[cig-commits] r16564 - short/3D/PyLith/trunk/doc/userguide/governingeqns
brad at geodynamics.org
brad at geodynamics.org
Mon Apr 19 17:46:53 PDT 2010
Author: brad
Date: 2010-04-19 17:46:53 -0700 (Mon, 19 Apr 2010)
New Revision: 16564
Modified:
short/3D/PyLith/trunk/doc/userguide/governingeqns/governingeqns.lyx
Log:
A few more edits.
Modified: short/3D/PyLith/trunk/doc/userguide/governingeqns/governingeqns.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/governingeqns/governingeqns.lyx 2010-04-19 23:16:23 UTC (rev 16563)
+++ short/3D/PyLith/trunk/doc/userguide/governingeqns/governingeqns.lyx 2010-04-20 00:46:53 UTC (rev 16564)
@@ -994,6 +994,16 @@
i.e., the same trial function is used for each component, and we employ quadratur
e to integrate the quantities over the domain of each finite-element cell.
We formulate the residual using
+\begin_inset Note Greyedout
+status open
+
+\begin_layout Plain Layout
+Need to refine the finite-element implementation stuff.
+\end_layout
+
+\end_inset
+
+
\begin_inset Formula \[
R=-\sum_{\text{vol cells}}\sum_{\text{quad pt}}\frac{1}{2}\sigma_{ij}(x_{q},t+\Delta t)(\phi_{,j}(x_{q})+\phi_{,i}(x_{q}))\: w_{q}|J_{m}(x_{q})|+\sum_{\text{vol cells}}\sum_{\text{quad pt}}f_{i}(x_{q},t+\Delta t)\phi(x_{q})\, w_{q}|J_{m}(x_{q})|+\sum_{\text{tract cells}}\sum_{\text{quad pt}}T_{i}(x_{q},t+\Delta t)\phi_{i}(x_{q})\, w_{q}|J_{m}(x_{q})|,\]
@@ -1104,6 +1114,16 @@
.
Applying the same finite-element formulation to the Jacobian that we did
to the residual, the expression for the Jacobian sparse matrix is
+\begin_inset Note Greyedout
+status open
+
+\begin_layout Plain Layout
+Need to refint the finite-element stuff.
+\end_layout
+
+\end_inset
+
+
\begin_inset Formula \[
J=\sum_{\text{vol cells }}\sum_{\text{quad pt}}\frac{1}{4}C_{ijkl}(\phi_{,l}(x_{q})+\phi_{,k}(x_{q}))(\phi_{,j}(x_{q})+\phi_{,i}(x_{q}))\ w_{q}|J_{m}(x_{q})|\]
More information about the CIG-COMMITS
mailing list