[cig-commits] r7731 - in long/3D/Gale/trunk: . input

walter at geodynamics.org walter at geodynamics.org
Fri Jul 20 23:08:06 PDT 2007


Author: walter
Date: 2007-07-20 23:08:05 -0700 (Fri, 20 Jul 2007)
New Revision: 7731

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/extension.xml
   long/3D/Gale/trunk/input/shortening.xml
Log:
 r1886 at earth:  boo | 2007-07-20 10:47:02 -0700
 Make the shortening and extension examples use numbers closer to the real sandbox experiment.  The viscosity is still wrong, because setting it to the right number causes the code to crash.



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1885
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1886

Modified: long/3D/Gale/trunk/input/extension.xml
===================================================================
--- long/3D/Gale/trunk/input/extension.xml	2007-07-21 06:07:51 UTC (rev 7730)
+++ long/3D/Gale/trunk/input/extension.xml	2007-07-21 06:08:05 UTC (rev 7731)
@@ -308,10 +308,10 @@
     </struct>
     <struct name="pdmsShape">
       <param name="Type">Box</param>
-      <param name="startX">0.5</param>
-      <param name="endX">1.5</param>
+      <param name="startX">0.05</param>
+      <param name="endX">0.15</param>
       <param name="startY">0.0</param>
-      <param name="endY">0.05</param>
+      <param name="endY">0.005</param>
       <param name="startZ">minZ</param>
       <param name="endZ">maxZ</param>
     </struct>
@@ -324,13 +324,13 @@
     </struct>
     <struct name="crustViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">10.0</param>
+      <param name="eta0">1.0</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="softeningStrain">0.01</param>
       <param name="initialDamageFraction">0.0</param>
       <param name="initialDamageWavenumber">0.5</param>
       <param name="initialDamageFactor">0.5</param>
@@ -345,13 +345,13 @@
       <param name="StrainWeakening">strainWeakening</param>
       <param name="cohesion">10.0</param>
       <param name="cohesionAfterSoftening">0.1</param>
-      <param name="frictionCoefficient">0.0</param>
-      <param name="frictionCoefficientAfterSoftening">0.0</param>
+      <param name="frictionCoefficient">0.404026225835</param>
+      <param name="frictionCoefficientAfterSoftening">0.363970234266</param>
       <param name="minimumYieldStress">0.00001</param>
     </struct>
     <struct name="pdmsViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0</param>
+      <param name="eta0">50000.0</param>
     </struct>
     <struct name="storeViscosity">
       <param name="Type">StoreVisc</param>
@@ -364,7 +364,7 @@
     <struct name="crust">
       <param name="Type">RheologyMaterial</param>
       <param name="Shape">crustShape</param>
-      <param name="density">1.0</param>
+      <param name="density">1560.0</param>
       <list name="Rheology">
         <param>crustViscosity</param>
         <param>yielding</param>
@@ -375,7 +375,7 @@
     <struct name="pdms">
       <param name="Type">RheologyMaterial</param>
       <param name="Shape">pdmsShape</param>
-      <param name="density">1.0</param>
+      <param name="density">965.0</param>
       <list name="Rheology">
         <param>pdmsViscosity</param>
         <param>storeViscosity</param>
@@ -430,8 +430,8 @@
   <param name="minX">0.0f</param>
   <param name="minY">0.0f</param>
   <param name="minZ">0.0f</param>
-  <param name="maxX">2.0f</param>
-  <param name="maxY">0.35f</param>
+  <param name="maxX">0.20f</param>
+  <param name="maxY">0.035f</param>
   <param name="maxZ">1.0f</param>
   <param name="elementResI">64</param>
   <param name="elementResJ">8</param>
@@ -511,11 +511,6 @@
             <param name="type">double</param>
             <param name="value">0.0</param>
           </struct>
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
         </list>
       </struct>
       <struct>
@@ -525,13 +520,8 @@
           <struct>
             <param name="name">vx</param>
             <param name="type">double</param>
-            <param name="value">1.0</param>
+            <param name="value">0.01</param>
           </struct>
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
         </list>
       </struct>
 
@@ -550,7 +540,7 @@
       <struct>
         <param name="type">FrictionVC</param>
         <param name="wall">bottom</param>
-        <param name="StaticFriction">0.1</param>
+        <param name="StaticFriction">0.33</param>
         <param name="includeUpperX">False</param>
         <param name="StressField">StressField</param>
         <param name="PressureField">PressureField</param>
@@ -581,13 +571,14 @@
     </list>
   </struct>
 
-  <param name="StepFunctionLowerOffset">0.99</param>
-  <param name="StepFunctionUpperOffset">1.01</param>
-  <param name="StepFunctionValue">1.0</param>
+  <param name="StepFunctionLowerOffset">0.099</param>
+  <param name="StepFunctionUpperOffset">0.101</param>
+  <param name="StepFunctionValue">0.010</param>
   <param name="StepFunctionDim">0</param>
   <param name="StepFunctionLessThan">False</param>
   <param name="checkpointEvery">1</param>
   <param name="mgLevels">2</param>
+  <param name="dtFactor">1.0</param>
 
 <!--   <param name="journal.info">True</param> -->
 <!--   <param name="journal.debug">True</param> -->

Modified: long/3D/Gale/trunk/input/shortening.xml
===================================================================
--- long/3D/Gale/trunk/input/shortening.xml	2007-07-21 06:07:51 UTC (rev 7730)
+++ long/3D/Gale/trunk/input/shortening.xml	2007-07-21 06:08:05 UTC (rev 7731)
@@ -325,8 +325,8 @@
       <param name="Type">StrainWeakening</param>
       <param name="TimeIntegrator">timeIntegrator</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="initialSofteningStrain">0.5</param>
-      <param name="finalSofteningStrain">1.0</param>
+      <param name="initialSofteningStrain">0.0</param>
+      <param name="finalSofteningStrain">0.01</param>
       <param name="initialDamageFraction">0.0</param>
       <param name="initialDamageWavenumber">0.5</param>
       <param name="initialDamageFactor">0.5</param>
@@ -340,21 +340,21 @@
       <param name="Context">context</param>
       <param name="StrainWeakening">sandStrainWeakening</param>
       <param name="cohesion">10.0</param>
-      <param name="cohesionAfterSoftening">1.0</param>
+      <param name="cohesionAfterSoftening">0.1</param>
       <param name="frictionCoefficient">0.726542528005</param>
       <param name="frictionCoefficientAfterSoftening">0.600860619027</param>
       <param name="minimumYieldStress">0.00001</param>
     </struct>
     <struct name="sandViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0e2</param>
+      <param name="eta0">1.0</param>
     </struct>
     <struct name="beadsStrainWeakening">
       <param name="Type">StrainWeakening</param>
       <param name="TimeIntegrator">timeIntegrator</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="initialSofteningStrain">0.5</param>
-      <param name="finalSofteningStrain">1.0</param>
+      <param name="initialSofteningStrain">0.0</param>
+      <param name="finalSofteningStrain">0.01</param>
       <param name="initialDamageFraction">0.0</param>
       <param name="initialDamageWavenumber">0.5</param>
       <param name="initialDamageFactor">0.5</param>
@@ -368,14 +368,14 @@
       <param name="Context">context</param>
       <param name="StrainWeakening">beadsStrainWeakening</param>
       <param name="cohesion">10.0</param>
-      <param name="cohesionAfterSoftening">1.0</param>
+      <param name="cohesionAfterSoftening">0.1</param>
       <param name="frictionCoefficient">0.404026225835</param>
       <param name="frictionCoefficientAfterSoftening">0.363970234266</param>
       <param name="minimumYieldStress">0.00001</param>
     </struct>
     <struct name="beadsViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0e2</param>
+      <param name="eta0">1.0</param>
     </struct>
     <struct name="storeViscosity">
       <param name="Type">StoreVisc</param>
@@ -461,7 +461,7 @@
   <param name="maxY">0.035f</param>
   <param name="maxZ">1.0f</param>
   <param name="elementResI">64</param>
-  <param name="elementResJ">16</param>
+  <param name="elementResJ">8</param>
   <param name="elementResK">4</param>
   <param name="allowUnbalancing">True</param>
   <param name="buildNodeNeighbourTbl">True</param>
@@ -592,4 +592,5 @@
   </struct>
 
   <param name="checkpointEvery">1</param>
+  <param name="dtFactor">1.0</param>
 </StGermainData>



More information about the cig-commits mailing list