[cig-commits] r8034 - in long/3D/Gale/trunk: . input/examples
walter at geodynamics.org
walter at geodynamics.org
Wed Sep 26 13:12:18 PDT 2007
Author: walter
Date: 2007-09-26 13:12:18 -0700 (Wed, 26 Sep 2007)
New Revision: 8034
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/examples/tibet.xml
Log:
r1941 at earth: boo | 2007-09-26 13:07:07 -0700
Make tibet use realistic viscosity and cohesion. Unfortunately, it required turning off normal stress boundaries and making the velocity of the left side 100 times slower
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1940
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1941
Modified: long/3D/Gale/trunk/input/examples/tibet.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet.xml 2007-09-26 20:12:15 UTC (rev 8033)
+++ long/3D/Gale/trunk/input/examples/tibet.xml 2007-09-26 20:12:18 UTC (rev 8034)
@@ -292,7 +292,7 @@
<param name="ContinuityForceVector">cont_force</param>
<param name="killNonConvergent">false</param>
<param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
- <param name="nonLinearTolerance">nonLinearTolerance</param>
+ <param name="nonLinearTolerance">5e-3</param>
<param name="makeConvergenceFile">false</param>
</struct>
<struct name="buoyancyForceTerm">
@@ -337,15 +337,16 @@
<param name="endZ">maxZ</param>
</struct>
- <struct name="stressBC">
- <param name="Type">StressBC</param>
- <param name="ForceVector">mom_force</param>
- <param name="Swarm">picIntegrationPoints</param>
- <param name="wall">bottom</param>
- <param name="y_type">double</param>
+<!-- <struct name="stressBC"> -->
+<!-- <param name="Type">StressBC</param> -->
+<!-- <param name="ForceVector">mom_force</param> -->
+<!-- <param name="Swarm">picIntegrationPoints</param> -->
+<!-- <param name="wall">bottom</param> -->
+<!-- <param name="y_type">double</param> -->
<!-- This number is just adding up rho*g*h at the edge -->
- <param name="y_value">3077200000</param>
- </struct>
+<!-- <param name="y_value">3077200000</param> -->
+<!-- </struct> -->
+
<struct name="surfaceAdaptor">
<param name="Type">SurfaceAdaptor</param>
<param name="mesh">mesh-linear</param>
@@ -402,10 +403,35 @@
<struct name="crustViscosity">
<param name="Type">FrankKamenetskii</param>
<param name="TemperatureField">TemperatureField</param>
- <param name="eta0">2.37419189702e22</param>
- <param name="theta">3.16724221869e-3</param>
+ <!-- at T=273, viscosity=1e25, at T=1273, viscosity=1e20 -->
+ <param name="eta0">2.31739465e26</param>
+ <param name="theta">0.011512925465</param>
</struct>
+ <struct name="strainWeakening">
+ <param name="Type">StrainWeakening</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="softeningStrain">0.1</param>
+ <param name="initialDamageFraction">0.0</param>
+ <param name="initialDamageWavenumber">0.5</param>
+ <param name="initialDamageFactor">0.5</param>
+ <param name="healingRate">0.0</param>
+ </struct>
+ <struct name="yielding">
+ <param name="Type">MohrCoulomb</param>
+ <param name="PressureField">PressureField</param>
+ <param name="VelocityGradientsField">VelocityGradientsField</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="Context">context</param>
+ <param name="StrainWeakening">strainWeakening</param>
+ <param name="cohesion">44e6</param>
+ <param name="cohesionAfterSoftening">4e6</param>
+ <param name="frictionCoefficient">0.57735026919</param>
+ <param name="frictionCoefficientAfterSoftening">0.57735026919</param>
+ <param name="minimumYieldStress">0.0001</param>
+ </struct>
+
<struct name="crustViscous">
<param name="Type">RheologyMaterial</param>
<param name="Shape">crustShape</param>
@@ -420,6 +446,7 @@
</list>
<list name="Rheology">
<param>crustViscosity</param>
+ <param>yielding</param>
<param>storeViscosity</param>
<param>storeStress</param>
</list>
@@ -437,8 +464,8 @@
<struct name="mantleViscosity">
<param name="Type">FrankKamenetskii</param>
<param name="TemperatureField">TemperatureField</param>
- <param name="eta0">2.37419189702e22</param>
- <param name="theta">3.16724221869e-3</param>
+ <param name="eta0">2.31739465e26</param>
+ <param name="theta">0.011512925465</param>
</struct>
<struct name="mantleViscous">
@@ -557,8 +584,8 @@
<param name="maxX">1000000.0f</param>
<param name="maxY">100000f</param>
<param name="maxZ">1000000f</param>
- <param name="elementResI">32</param>
- <param name="elementResJ">16</param>
+ <param name="elementResI">64</param>
+ <param name="elementResJ">8</param>
<param name="elementResK">4</param>
<param name="allowUnbalancing">True</param>
<param name="buildNodeNeighbourTbl">True</param>
@@ -646,7 +673,8 @@
<param name="name">vx</param>
<param name="type">double</param>
<!-- 1 cm/yr -->
- <param name="value">3.17097919838e-10</param>
+ <param name="value">3.17097919838e-12</param>
+<!-- <param name="value">3.17097919838e-10</param> -->
</struct>
<struct>
<param name="name">vy</param>
@@ -655,6 +683,19 @@
</struct>
</list>
</struct>
+
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">bottom</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>
@@ -762,4 +803,10 @@
<param name="TemperatureProfileExponentialCoefficient2">1.0e-4</param>
<param name="checkpointEvery">1</param>
+ <param name="dtFactor">0.1</param>
+
+<!-- <param name="journal.info">True</param> -->
+<!-- <param name="journal.debug">True</param> -->
+<!-- <param name="journal-level.info">2</param> -->
+<!-- <param name="journal-level.debug">2</param> -->
</StGermainData>
More information about the cig-commits
mailing list