[cig-commits] r22249 - short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction

[brad at geodynamics.org]

Author: brad
Date: 2013-06-13 09:53:51 -0700 (Thu, 13 Jun 2013)
New Revision: 22249

Modified:
   short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc
Log:
Remove effective slip rate stuff from rate-state friction model (not necessary in current formulation).

Modified: short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-13 15:07:31 UTC (rev 22248)
+++ short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-13 16:53:51 UTC (rev 22249)
@@ -419,22 +419,18 @@
   // thetaTpdt = thetaT * exp(-slipRate/L * dt)
   //             + dt - 0.5*(sliprate/L)*dt**2 + 1.0/6.0*(slipRate/L)*dt**3;
 
-  // Since regulatized friction -> 0 as slipRate -> 0, limit slip
-  // rate to some minimum value
-  const PylithScalar slipRateEff = std::max(1.0e-2*_linearSlipRate, slipRate);
-
   const PylithScalar dt = _dt;
   const PylithScalar thetaTVertex = stateVars[s_state];
   const PylithScalar L = properties[p_L];
-  const PylithScalar vDtL = slipRateEff * dt / L;
+  const PylithScalar vDtL = slipRate * dt / L;
   const PylithScalar expTerm = exp(-vDtL);
 
   PylithScalar thetaTpdtVertex = 0.0;
   if (vDtL > 1.0e-20) {
-    thetaTpdtVertex = thetaTVertex * expTerm + L / slipRateEff * (1 - expTerm);
+    thetaTpdtVertex = thetaTVertex * expTerm + L / slipRate * (1 - expTerm);
     PetscLogFlops(7);
   } else {
-    thetaTpdtVertex = thetaTVertex * expTerm + dt - 0.5 * slipRateEff/L * dt*dt;
+    thetaTpdtVertex = thetaTVertex * expTerm + dt - 0.5 * slipRate/L * dt*dt;
     PetscLogFlops(9);
   } // if/else