[cig-commits] commit: Small edits.

Mercurial hg at geodynamics.org
Thu Feb 7 09:12:14 PST 2013


changeset:   168:443f4b2a5a50
tag:         tip
user:        Brad Aagaard <baagaard at usgs.gov>
date:        Thu Feb 07 09:12:09 2013 -0800
files:       faultRup.tex
description:
Small edits.


diff -r f66e7a77e6a8 -r 443f4b2a5a50 faultRup.tex
--- a/faultRup.tex	Thu Feb 07 09:00:04 2013 -0800
+++ b/faultRup.tex	Thu Feb 07 09:12:09 2013 -0800
@@ -351,7 +351,7 @@ this effective plastic strain and the el
 this effective plastic strain and the elastic strain. The fault
 tractions are associated with the total strain, not the effective
 plastic strain. This illustrates a key difference between this
-approach and the domain decomposition approach, in which the Lagrange
+approach and the domain decomposition approach in which the Lagrange
 multipliers and the constraint equation directly relate the fault slip
 to the fault tractions (Lagrange multipliers). One implication of this
 difference is that when using double couple point forces, the body
@@ -413,7 +413,7 @@ We want to solve these equations for the
 We want to solve these equations for the coefficients $\bm{u}_n$
 and $\bm{l}_p$ subject to $\bm{u} = \bm{u}_0 \text{ on
 }S_u$. When we prescribe the slip, we specify $\bm{d}$ on $S_f$,
-and when we use a fault constitutive model, we specify how the
+and when we use a fault constitutive model we specify how the
 Lagrange multipliers $\bm{l}$ depend on the fault slip, slip rate,
 and state variables.
 
@@ -445,9 +445,9 @@ time-dependence only enters through the 
 time-dependence only enters through the constitutive models and the
 loading conditions. As a result, the quasi-static simulations are a
 series of static problems with potentially time-varying physical
-properties and boundary conditions. The stability of the solution is
-limited to resolving these temporal variations. Considering the
-deformation at time $t+\Delta t$,
+properties and boundary conditions. The temporal accuracy of the
+solution is limited to resolving these temporal
+variations. Considering the deformation at time $t+\Delta t$,
 \begin{linenomath*}\begin{gather}
   \label{eqn:quasi-static:residual:elasticity}
   \begin{split}
@@ -479,7 +479,7 @@ with parallel processing
 with parallel processing
 \citep{PETSC:efficient,PETSc:manual}. In solving the
 system, we compute the residual (i.e., $\bm{r} = \bm{b} -
-\bm{A} \cdot \bm{u}$ and the Jacobian of the system
+\bm{A} \cdot \bm{u}$) and the Jacobian of the system
 ($\bm{A}$). In our case the solution is $\bm{u} =
 \left( \begin{smallmatrix} \bm{u}_n \\
     \bm{l}_n \end{smallmatrix} \right)$, and the residual is
@@ -515,7 +515,7 @@ This matches the tangent stiffness matri
 This matches the tangent stiffness matrix in conventional solid
 mechanics finite-element formulations. In computing the residual, we
 use the expression given in equation~(\ref{eqn:residual:elasticity})
-with an implementation for infinitesimal strain and an implementation
+with one implementation for infinitesimal strain and another implementation
 for small strain and rigid body motion. Following a similar procedure,
 we find the portion of the Jacobian associated with the constraints,
 equation~(\ref{eqn:quasi-static:residual:fault}), is
@@ -556,14 +556,13 @@ associated with the Lagrange multipliers
 
 The matrix $\bm{L}$ defined in
 equation~(\ref{eqn:jacobian:constraint}) is spectrally equivalent to
-the identity, because it just involves integration of products of the
-basis functions. This makes the traditional LBB stability criterion
+the identity, because it involves integration of products of the basis
+functions. This makes the traditional LBB stability criterion
 \citep{Brenner:Scott:2008} trivial to satisfy by choosing the space of
 Lagrange multipliers to be exactly the space of displacements,
-restricted to the fault. In simple terms to specify the problem we
-need to know the distance between any pair of vertices spanning the
-fault, which can be expressed as a relative displacement (i.e., fault
-slip).
+restricted to the fault. This means we simly need to know the distance
+between any pair of vertices spanning the fault, which can be
+expressed as a relative displacement, i.e., fault slip.
 
 % ------------------------------------------------------------------
 \subsection{Dynamic Simulations}



More information about the CIG-COMMITS mailing list