[cig-commits] r11976 - in long/3D/Gale/trunk: . input/benchmarks
walter at geodynamics.org
walter at geodynamics.org
Thu May 15 12:07:11 PDT 2008
Author: walter
Date: 2008-05-15 12:07:11 -0700 (Thu, 15 May 2008)
New Revision: 11976
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/benchmarks/shortening.xml
Log:
r2156 at earth: boo | 2008-05-15 12:00:19 -0700
Mostly working shortening model. There still may be a problem with large slopes at later times causing large velocities and inhibiting convergence.
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2155
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2156
Modified: long/3D/Gale/trunk/input/benchmarks/shortening.xml
===================================================================
--- long/3D/Gale/trunk/input/benchmarks/shortening.xml 2008-05-15 19:07:08 UTC (rev 11975)
+++ long/3D/Gale/trunk/input/benchmarks/shortening.xml 2008-05-15 19:07:11 UTC (rev 11976)
@@ -283,7 +283,7 @@
<struct name="uzawa">
<param name="Type">Stokes_SLE_UzawaSolver</param>
<param name="Preconditioner">preconditioner</param>
- <param name="tolerance">1.0e-5</param>
+ <param name="tolerance">linearTolerance</param>
<param name="maxIterations">maxUzawaIterations</param>
</struct>
<struct name="stokesEqn">
@@ -395,7 +395,7 @@
</struct>
<struct name="beadsViscosity">
<param name="Type">MaterialViscosity</param>
- <param name="eta0">100000.0</param>
+ <param name="eta0">1.0e12</param>
</struct>
<struct name="cornerViscosity">
<param name="Type">MaterialViscosity</param>
@@ -431,17 +431,6 @@
<param>storeStress</param>
</list>
</struct>
-<!-- <struct name="corner"> -->
-<!-- <param name="Type">RheologyMaterial</param> -->
-<!-- <param name="Shape">cornerShape</param> -->
-<!-- <param name="density">1560</param> -->
-<!-- <list name="Rheology"> -->
-<!-- <param>sandViscosity</param> -->
-<!-- <param>sandYielding</param> -->
-<!-- <param>storeViscosity</param> -->
-<!-- <param>storeStress</param> -->
-<!-- </list> -->
-<!-- </struct> -->
<struct name="surfaceAdaptor">
<param name="Type">SurfaceAdaptor</param>
<param name="mesh">mesh-linear</param>
@@ -477,8 +466,8 @@
<struct name="conditionFunctions">
<param name="Type">StgFEM_StandardConditionFunctions</param>
</struct>
- <struct name="kineticFriction">
- <param name="Type">StressBC</param>
+ <struct name="kineticFrictionBottom">
+ <param name="Type">KineticFriction</param>
<param name="ForceVector">mom_force</param>
<param name="Swarm">picIntegrationPoints</param>
<param name="wall">bottom</param>
@@ -490,6 +479,21 @@
<param name="vx_type">double</param>
<param name="vx_value">0.0</param>
</struct>
+
+ <struct name="kineticFrictionRight">
+ <param name="Type">KineticFriction</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="wall">right</param>
+ <param name="StressField">StressField</param>
+ <param name="PressureField">PressureField</param>
+ <param name="VelocityField">VelocityField</param>
+ <param name="friction_type">double</param>
+ <param name="friction_value">0.34432761329</param>
+ <param name="vy_type">double</param>
+ <param name="vy_value">0.0</param>
+ </struct>
+
</struct>
<list name="plugins">
<param>Underworld_EulerDeform</param>
@@ -605,21 +609,42 @@
</struct>
</list>
</struct>
-<!-- <struct> -->
-<!-- <param name="type">StaticFrictionVC</param> -->
-<!-- <param name="wall">bottom</param> -->
-<!-- <param name="StaticFriction">0.34432761329</param> -->
-<!-- <param name="includeUpperX">False</param> -->
-<!-- <param name="StressField">StressField</param> -->
-<!-- <param name="PressureField">PressureField</param> -->
-<!-- <list name="variables"> -->
-<!-- <struct> -->
-<!-- <param name="name">vx</param> -->
-<!-- <param name="type">double</param> -->
-<!-- <param name="value">0.0</param> -->
-<!-- </struct> -->
-<!-- </list> -->
-<!-- </struct> -->
+ <struct>
+ <param name="type">StaticFrictionVC</param>
+ <param name="wall">bottom</param>
+ <param name="StaticFriction">0.34432761329</param>
+ <param name="includeUpperX">False</param>
+ <param name="StickOnFirstStep">False</param>
+ <param name="StressField">StressField</param>
+ <param name="PressureField">PressureField</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vx</param>
+<!-- <param name="type">func</param> -->
+<!-- <param name="value">StepFunction</param> -->
+ <param name="type">double</param>
+ <param name="value">0.0</param>
+ </struct>
+ </list>
+ </struct>
+
+ <struct>
+ <param name="type">StaticFrictionVC</param>
+ <param name="wall">right</param>
+ <param name="StaticFriction">0.34432761329</param>
+ <param name="includeLowerY">False</param>
+ <param name="StickOnFirstStep">False</param>
+ <param name="StressField">StressField</param>
+ <param name="PressureField">PressureField</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0.0</param>
+ </struct>
+ </list>
+ </struct>
+
</list>
</struct>
@@ -639,17 +664,17 @@
</list>
</struct>
- <param name="StepFunctionLowerOffset">0.0</param>
- <param name="StepFunctionUpperOffset">0.035</param>
- <param name="StepFunctionValue">-1.0</param>
- <param name="StepFunctionDim">0</param>
- <param name="StepFunctionLessThan">False</param>
-
- <param name="checkpointEvery">1</param>
+<!-- <param name="checkpointEvery">1</param> -->
<param name="dtFactor">1.0</param>
<param name="maxUzawaIterations">20000</param>
+ <param name="linearTolerance">1.0e-5</param>
+ <param name="StepFunctionLowerOffset">0.38</param>
+ <param name="StepFunctionUpperOffset">0.40</param>
+ <param name="StepFunctionValue">-6.94444444444e-6</param>
+ <param name="StepFunctionDim">0</param>
+ <param name="StepFunctionLessThan">False</param>
<param name="journal.info">True</param>
<param name="journal.debug">True</param>
More information about the cig-commits
mailing list