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

walter at geodynamics.org walter at geodynamics.org
Wed Mar 4 09:09:49 PST 2009


Author: walter
Date: 2009-03-04 09:09:49 -0800 (Wed, 04 Mar 2009)
New Revision: 14225

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/benchmarks/geomod_2008_brittle.xml
Log:
 r2527 at dante:  boo | 2009-03-04 08:57:40 -0800
 Make geomod brittle use a transform'd boundary layer



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

Modified: long/3D/Gale/trunk/input/benchmarks/geomod_2008_brittle.xml
===================================================================
--- long/3D/Gale/trunk/input/benchmarks/geomod_2008_brittle.xml	2009-03-04 17:09:44 UTC (rev 14224)
+++ long/3D/Gale/trunk/input/benchmarks/geomod_2008_brittle.xml	2009-03-04 17:09:49 UTC (rev 14225)
@@ -114,26 +114,6 @@
       <param name="LinkedDofInfo">pressureLinkedDofs</param>
     </struct>
 
-    <struct name="nodalPressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Scalar</param>
-      <param name="DataType">Double</param>
-    </struct>
-    <struct name="nodalPressureDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <list name="BaseVariables">
-        <param>nodalPressure</param>
-      </list>
-    </struct>
-    <struct name="NodalPressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">nodalPressureDofLayout</param>
-      <param name="LinkedDofInfo">nodalPressureLinkedDofs</param>
-    </struct>
-
     <struct name="StressField">
       <param name="Type">StressField</param>
       <param name="StrainRateField">StrainRateField</param>
@@ -368,22 +348,32 @@
     </struct>
     <struct name="quartzYielding">
       <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">NodalPressureField</param>
+      <param name="PressureField">PressureField</param>
       <param name="VelocityGradientsField">VelocityGradientsField</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
       <param name="Context">context</param>
       <param name="StrainWeakening">quartzStrainWeakening</param>
       <param name="cohesion">30.0</param>
       <param name="cohesionAfterSoftening">30</param>
-      <!-- Angle of internal friction = 35 degrees -->
-      <param name="frictionCoefficient">0.70020753821</param>
-      <param name="frictionCoefficientAfterSoftening">0.70020753821</param>
+      <param name="frictionCoefficient">1.05737126344</param>
+      <param name="frictionCoefficientAfterSoftening">0.75168293455</param>
+<!--       <param name="frictionCoefficient">0.726542528005</param> -->
+<!--       <param name="frictionCoefficientAfterSoftening">0.600860619027</param> -->
+      <param name="boundaryCohesion">30</param>
+      <param name="boundaryCohesionAfterSoftening">30</param>
+      <param name="boundaryFrictionCoefficient">0.299314553072</param>
+      <param name="boundaryFrictionCoefficientAfterSoftening">0.257458783778</param>
+<!--       <param name="boundaryFrictionCoefficient">0.286745385759</param> -->
+<!--       <param name="boundaryFrictionCoefficientAfterSoftening">0.249328002843</param> -->
+      <param name="boundaryBottom">True</param>
+      <param name="boundaryRight">True</param>
     </struct>
-
     <struct name="quartzViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0e12</param>
+      <param name="eta0">1.0e2</param>
+<!--       <param name="eta0">1.0e12</param> -->
     </struct>
+
     <struct name="corundumStrainWeakening">
       <param name="Type">StrainWeakening</param>
       <param name="TimeIntegrator">timeIntegrator</param>
@@ -397,20 +387,32 @@
     </struct>
     <struct name="corundumYielding">
       <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">NodalPressureField</param>
+      <param name="PressureField">PressureField</param>
       <param name="VelocityGradientsField">VelocityGradientsField</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
       <param name="Context">context</param>
       <param name="StrainWeakening">corundumStrainWeakening</param>
       <param name="cohesion">30.0</param>
       <param name="cohesionAfterSoftening">30</param>
-      <param name="frictionCoefficient">0.726542528005</param>
-      <param name="frictionCoefficientAfterSoftening">0.726542528005</param>
+      <param name="frictionCoefficient">1.05737126344</param>
+      <param name="frictionCoefficientAfterSoftening">0.75168293455</param>
+<!--       <param name="frictionCoefficient">0.726542528005</param> -->
+<!--       <param name="frictionCoefficientAfterSoftening">0.600860619027</param> -->
+      <param name="boundaryCohesion">30</param>
+      <param name="boundaryCohesionAfterSoftening">30</param>
+      <param name="boundaryFrictionCoefficient">0.497230605386</param>
+      <param name="boundaryFrictionCoefficientAfterSoftening">0.468805151222</param>
+<!--       <param name="boundaryFrictionCoefficient">0.445228685309</param> -->
+<!--       <param name="boundaryFrictionCoefficientAfterSoftening">0.42447481621</param> -->
+      <param name="boundaryBottom">True</param>
+      <param name="boundaryRight">True</param>
     </struct>
     <struct name="corundumViscosity">
       <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0e12</param>
+      <param name="eta0">1.0e2</param>
+<!--       <param name="eta0">1.0e12</param> -->
     </struct>
+
     <struct name="storeViscosity">
       <param name="Type">StoreVisc</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
@@ -441,6 +443,7 @@
         <param>storeStress</param>
       </list>
     </struct>
+
     <struct name="escapedRoutine">
       <param name="Type">EscapedRoutine</param>
       <param name="idealParticleCount">0</param>
@@ -461,7 +464,6 @@
 
   </struct>
   <list name="plugins">
-<!--     <param>Gale_SurfaceProcess</param> -->
     <param>Underworld_EulerDeform</param>
     <param>Underworld_VTKOutput</param>
   </list>
@@ -496,6 +498,7 @@
         <param name="remesher">velocityRemesher</param>
         <param name="velocityField">VelocityField</param>
         <param name="wrapTop">True</param>
+        <param name="floatRightTop">True</param>
         <list name="fields">
           <struct>
             <param name="field">VelocityField</param>
@@ -509,9 +512,6 @@
       </struct>
     </list>
   </struct>
-  <struct name="SurfaceProcess">
-    <param name="mesh">mesh-linear</param>
-  </struct>
   <struct name="velocityBCs">
     <param name="type">CompositeVC</param>
     <list name="vcList">
@@ -545,7 +545,8 @@
             <param name="name">vx</param>
             <param name="type">double</param>
             <!-- 2.5 cm/hr -->
-            <param name="value">6.94444444444e-6</param>
+            <param name="value">6.94444444444e4</param>
+<!--             <param name="value">6.94444444444e-6</param> -->
           </struct>
           <struct>
             <param name="name">vy</param>
@@ -556,12 +557,15 @@
       </struct>
       <struct>
         <param name="type">WallVC</param>
-        <param name="wall">right</param>
+        <param name="wall">bottom</param>
         <list name="variables">
           <struct>
             <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
+            <param name="type">func</param>
+            <param name="value">StepFunction</param>
+<!--             <param name="type">double</param> -->
+<!--             <param name="value">6.94444444444e4</param> -->
+<!--             <param name="value">6.94444444444e-6</param> -->
           </struct>
           <struct>
             <param name="name">vy</param>
@@ -572,12 +576,12 @@
       </struct>
       <struct>
         <param name="type">WallVC</param>
-        <param name="wall">bottom</param>
+        <param name="wall">right</param>
         <list name="variables">
           <struct>
             <param name="name">vx</param>
-            <param name="type">func</param>
-            <param name="value">StaticFriction</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
           </struct>
           <struct>
             <param name="name">vy</param>
@@ -589,50 +593,17 @@
     </list>
   </struct>
 
-  <struct name="stressICs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">AllNodesVC</param>
-        <list name="variables">
-          <struct>
-            <param name="name">stress</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct> 
-    </list>
-  </struct>
+  <param name="linearTolerance">1.0e-7</param>
+  <param name="nonLinearTolerance">1.0e-3</param>
 
-<!--   <param name="checkpointEvery">1</param> -->
-<!--   <param name="dtFactor">1.0</param> -->
-<!--   <param name="dt">3.6</param> -->
+  <param name="StepFunctionLowerValue">6.94444444444e4</param>
+<!--   <param name="StepFunctionLowerValue">6.94444444444e-6</param> -->
+  <param name="StepFunctionUpperValue">0</param>
+  <param name="StepFunctionLowerOffset">0.23</param>
+  <param name="StepFunctionUpperOffset">0.23</param>
 
-  <param name="maxUzawaIterations">20000</param>
-  <param name="linearTolerance">1.0e-5</param>
-
-  <param name="StaticFrictionCoefficient">0.267949192431</param>
-  <param name="StaticFrictionNormalDim">1</param>
-  <param name="StaticFrictionTangentDim">0</param>
-
-  <param name="StepFunctionProduct1Start">-100</param>
-  <param name="StepFunctionProduct1End">0.23</param>
-  <param name="StepFunctionProduct1Value">6.94444444444e-6</param>
-  <param name="StepFunctionProduct1Dim">0</param>
-  <param name="StaticFrictionLowerV">6.94444444444e-6</param>
-  <param name="StaticFrictionUpperV">0</param>
-
-
-
-  <param name="StaticFriction2Coefficient">0.267949192431</param>
-  <param name="StaticFriction2NormalDim">0</param>
-  <param name="StaticFriction2TangentDim">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> -->
-  <param name="diffusionCoefficient">1.0e-8</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