[cig-commits] r14539 - in long/3D/Gale/trunk: . input/examples
walter at geodynamics.org
walter at geodynamics.org
Mon Mar 30 02:43:59 PDT 2009
Author: walter
Date: 2009-03-30 02:43:59 -0700 (Mon, 30 Mar 2009)
New Revision: 14539
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/examples/rifting.xml
Log:
r2633 at dante: boo | 2009-03-30 02:39:29 -0700
Add the final version of rifting
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2632
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2633
Modified: long/3D/Gale/trunk/input/examples/rifting.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/rifting.xml 2009-03-30 09:43:44 UTC (rev 14538)
+++ long/3D/Gale/trunk/input/examples/rifting.xml 2009-03-30 09:43:59 UTC (rev 14539)
@@ -299,7 +299,7 @@
<param name="ContinuityForceVector">cont_force</param>
<param name="killNonConvergent">false</param>
<param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
- <param name="nonLinearTolerance">1.0e-2</param>
+ <param name="nonLinearTolerance">nlinearTolerance</param>
<param name="makeConvergenceFile">false</param>
</struct>
<struct name="c_matrix">
@@ -316,13 +316,6 @@
<param name="storeVisc">storeViscosity</param>
<param name="StiffnessMatrix">c_matrix</param>
</struct>
- <struct name="buoyancyForceTerm">
- <param name="Type">BuoyancyForceTerm</param>
- <param name="ForceVector">mom_force</param>
- <param name="TemperatureField">TemperatureField</param>
- <param name="Swarm">picIntegrationPoints</param>
- <param name="gravity">9.81</param>
- </struct>
<struct name="background">
<param name="Type">Everywhere</param>
</struct>
@@ -340,7 +333,37 @@
<param>true</param>
</list>
</struct>
+ <struct name="storeViscosity">
+ <param name="Type">StoreVisc</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+ <struct name="storeStress">
+ <param name="Type">StoreStress</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+<!--============== END "BASIC COMPONENTS" =================-->
+ <struct name="hydrostaticTerm">
+ <param name="Type">HydrostaticTerm</param>
+ <param name="upperDensity">2800</param>
+ <param name="upperAlpha">3.0e-5</param>
+ <param name="lowerDensity">3300</param>
+ <param name="lowerAlpha">3.0e-5</param>
+ <param name="height">maxY</param>
+ <param name="materialBoundary">10000</param>
+ <param name="T_0">TemperatureProfileTop</param>
+ <param name="linearCoefficient">TemperatureProfileLinearCoefficient</param>
+ <param name="exponentialCoefficient1">TemperatureProfileExponentialCoefficient1</param>
+ <param name="exponentialCoefficient2">TemperatureProfileExponentialCoefficient2</param>
+ <param name="gravity">gravity</param>
+ </struct>
+ <struct name="buoyancyForceTerm">
+ <param name="Type">BuoyancyForceTerm</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="TemperatureField">TemperatureField</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="gravity">9.81</param>
+ </struct>
<struct name="boxShape">
<param name="Type">Box</param>
<param name="startX">minX</param>
@@ -351,54 +374,65 @@
<param name="endZ">maxZ</param>
</struct>
- <struct name="dikeShape">
- <param name="Type">Box</param>
- <param name="startX">17500</param>
- <param name="endX">22500</param>
- <param name="startY">15000</param>
- <param name="endY">maxY</param>
- <param name="startZ">minZ</param>
- <param name="endZ">maxZ</param>
- </struct>
+<!-- ============== Crust ================ -->
+ <struct name="crustShape">
+ <param name="Type">Box</param>
+ <param name="startX">minX</param>
+ <param name="endX">maxX</param>
+ <param name="startY">10000</param>
+ <param name="endY">maxY</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ </struct>
+ <struct name="crustViscosity">
+ <param name="Type">FrankKamenetskii</param>
+ <param name="TemperatureField">TemperatureField</param>
+<!-- ref. visc is 1e22 so this times eta0 gives -->
+<!-- viscosity=1e25, at T=1273 & viscosity=1e20 at T=1273 -->
+ <param name="eta0">2.31739465e4</param>
+ <param name="theta">0.011512925465</param>
+ </struct>
+ <struct name="crust">
+ <param name="Type">RheologyMaterial</param>
+ <param name="Shape">crustShape</param>
+ <param name="density">2800</param>
+ <param name="alpha">3.0e-5</param>
+ <param name="diffusivity">1.0e-6</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">1.0e-12</param>
+ <param name="lambda">0.0</param>
+ </struct>
+ </list>
+ <list name="Rheology">
+ <param>crustViscosity</param>
+ <param>yielding</param>
+ <param>storeViscosity</param>
+ <param>storeStress</param>
+ </list>
+ </struct>
- <struct name="crustShape">
- <param name="Type">Box</param>
- <param name="startX">minX</param>
- <param name="endX">maxX</param>
- <param name="startY">15000</param>
- <param name="endY">maxY</param>
- <param name="startZ">minZ</param>
- <param name="endZ">maxZ</param>
- </struct>
-
- <struct name="storeViscosity">
- <param name="Type">StoreVisc</param>
- <param name="MaterialPointsSwarm">materialSwarm</param>
- </struct>
- <struct name="storeStress">
- <param name="Type">StoreStress</param>
- <param name="MaterialPointsSwarm">materialSwarm</param>
- </struct>
-
- <struct name="crustViscosity">
- <param name="Type">FrankKamenetskii</param>
- <param name="TemperatureField">TemperatureField</param>
-<!-- at T=273, viscosity=1e25, at T=1273, viscosity=1e20 -->
- <param name="eta0">2.31739465</param>
-<!-- <param name="eta0">2.31739465e26</param> -->
- <param name="theta">0.011512925465</param>
- </struct>
-
+<!-- ============== Strain Weakening ================ -->
<struct name="strainWeakening">
<param name="Type">StrainWeakening</param>
<param name="TimeIntegrator">timeIntegrator</param>
<param name="MaterialPointsSwarm">materialSwarm</param>
<param name="finalSofteningStrain">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>
+ <param name="initialDamageFraction">0.5</param>
+ <param name="initialDamageWavenumber">0.0</param>
+ <param name="initialDamageFactor">0.9</param>
+ <param name="initialStrainShape">initStrainBox</param>
+ <param name="healingRate">6.0</param>
</struct>
+ <struct name="initStrainBox">
+ <param name="Type">Box</param>
+ <param name="startX">29500.0</param>
+ <param name="endX">45000.0</param>
+ <param name="startY">0000.0</param>
+ <param name="endY">maxY</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ </struct>
<struct name="yielding">
<param name="Type">DruckerPrager</param>
<param name="PressureField">PressureField</param>
@@ -412,55 +446,19 @@
<param name="frictionCoefficientAfterSoftening">0.57735026919</param>
</struct>
- <struct name="crust">
- <param name="Type">RheologyMaterial</param>
- <param name="Shape">crustShape</param>
- <param name="density">2800</param>
- <param name="alpha">3.0e-5</param>
- <param name="diffusivity">1.0e-6</param>
- <list name="heatingElements">
- <struct>
- <param name="Q">1.0e-12</param>
- <param name="lambda">0.0</param>
- </struct>
- </list>
- <list name="Rheology">
- <param>crustViscosity</param>
- <param>yielding</param>
- <param>storeViscosity</param>
- <param>storeStress</param>
- </list>
- </struct>
-
-
- <struct name="divergenceForce">
- <param name="Type">DivergenceForce</param>
- <param name="ForceVector">cont_force</param>
- <param name="Swarm">picIntegrationPoints</param>
- <param name="DomainShape">dikeShape</param>
- <param name="GeometryMesh">mesh-linear</param>
- <param name="force_type">func</param>
- <param name="force_value">StepFunction</param>
-<!-- <param name="force_type">double</param> -->
-<!-- <param name="force_value">-1.4e-13</param> -->
- </struct>
-
-
-
- <!-- Mantle -->
+<!-- ============== Mantle ================ -->
<struct name="mantleShape">
<param name="Type">Intersection</param>
<list name="shapes">
<param>boxShape</param>
- <param>!crustShape</param>
+ <param>!crustShape</param>
</list>
</struct>
<struct name="mantleViscosity">
<param name="Type">FrankKamenetskii</param>
<param name="TemperatureField">TemperatureField</param>
- <param name="eta0">2.31739465</param>
-<!-- <param name="eta0">2.31739465e26</param> -->
+ <param name="eta0">2.31739465e3</param> <!--gives max vis=1e3 -->
<param name="theta">0.011512925465</param>
</struct>
@@ -478,22 +476,43 @@
</list>
<list name="Rheology">
<param>mantleViscosity</param>
+<!-- <param>yielding</param> -->
<param>storeViscosity</param>
<param>storeStress</param>
</list>
</struct>
-
+<!-- ============== Diking Zone ================ -->
+ <struct name="dikeShape">
+ <param name="Type">Box</param>
+ <param name="startX">29500</param>
+ <param name="endX">30500</param>
+ <param name="startY">9000</param>
+ <param name="endY">17000</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ </struct>
+ <struct name="divergenceForce">
+ <param name="Type">DivergenceForce</param>
+ <param name="ForceVector">cont_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="DomainShape">dikeShape</param>
+ <param name="GeometryMesh">mesh-linear</param>
+ <param name="force_type">double</param>
+ <param name="force_value">-0.8e9</param> <!--0.62*0.25.*(2*4e-10)/1000*1e22 -->
+ </struct>
+
+<!-- ============== Stress BC on bottom ===================== -->
<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>
- <param name="y_value">622935000</param>
+ <param name="y_value">5.0852e+8</param> <!-- 16000*3300*(1-3e-5*1473/2)*9.81 -->
</struct>
- <!-- Temperature components -->
+<!-- ============== Temperature Components ===================== -->
<struct name="temperature">
<param name="Type">MeshVariable</param>
<param name="Rank">Scalar</param>
@@ -521,7 +540,8 @@
<param name="DofLayout">temperatureDofLayout</param>
<param name="BC">temperatureBCs</param>
<param name="IC">temperatureICs</param>
- <param name="LinkedDofInfo">temperatureLinkedDofs</param>
+ <param name="LinkedDofInfo">temperature
+LinkedDofs</param>
</struct>
<struct name="TemperatureGradientsField">
<param name="Type">OperatorFeVariable</param>
@@ -530,7 +550,6 @@
</struct>
<!-- Energy Equation -->
-
<!-- <struct name="residual"> -->
<!-- <param name="Type">ForceVector</param> -->
<!-- <param name="FeVariable">TemperatureField</param> -->
@@ -572,34 +591,8 @@
<!-- </struct> -->
</struct>
- <list name="plugins">
- <param>Underworld_EulerDeform</param>
- <param>Underworld_VTKOutput</param>
- </list>
- <param name="maxTimeSteps">10</param>
- <param name="outputEvery">1</param>
- <param name="dumpEvery">1</param>
- <param name="outputPath">./output.rifting</param>
- <param name="dim">2</param>
- <param name="shadowDepth">1</param>
- <param name="minX">0.0f</param>
- <param name="minY">0.0f</param>
- <param name="minZ">0.0f</param>
- <param name="maxX">40000.0f</param>
- <param name="maxY">20000f</param>
- <param name="maxZ">1000000f</param>
- <param name="elementResI">40</param>
- <param name="elementResJ">20</param>
- <param name="elementResK">4</param>
- <param name="allowUnbalancing">True</param>
- <param name="buildNodeNeighbourTbl">True</param>
- <param name="buildElementNeighbourTbl">True</param>
- <param name="particlesPerCell">30</param>
- <param name="seed">13</param>
- <param name="gaussParticlesX">2</param>
- <param name="gaussParticlesY">2</param>
- <param name="gaussParticlesZ">2</param>
- <param name="particlesPerCell">30</param>
+
+<!-- ============== Mesh ==================== -->
<struct name="EulerDeform">
<list name="systems">
<struct>
@@ -608,19 +601,27 @@
<param name="velocityField">VelocityField</param>
<param name="wrapTop">True</param>
<param name="staticBottom">True</param>
- <list name="fields">
+ <param name="staticFront">True</param>
+ <param name="staticBack">True</param>
+ <param name="staticLeft">True</param>
+<!-- <param name="staticLeftTop">True</param> -->
+ <param name="staticRight">True</param>
+<!-- <param name="staticRightTop">True</param> -->
+ <list name="fields"> -->
+ <struct>
+ <param name="field">VelocityField</param>
+ <param name="variable">velocity</param>
+ </struct>
<struct>
- <param name="field">VelocityField</param>
- <param name="variable">velocity</param>
- </struct>
- <struct>
<param name="field">PressureField</param>
<param name="variable">pressure</param>
</struct>
- </list>
+ </list>
</struct>
</list>
</struct>
+
+ <!-- ============== Velocity BCs ==================== -->
<struct name="velocityBCs">
<param name="type">CompositeVC</param>
<list name="vcList">
@@ -653,8 +654,7 @@
<struct>
<param name="name">vx</param>
<param name="type">double</param>
- <param name="value">-4e16</param>
-<!-- <param name="value">-4e-10</param> -->
+ <param name="value">-4e12</param> <!-- -4e-10*1e22 -->
</struct>
<struct>
<param name="name">vy</param>
@@ -670,9 +670,7 @@
<struct>
<param name="name">vx</param>
<param name="type">double</param>
- <!-- 1 cm/yr -->
- <param name="value">4e16</param>
-<!-- <param name="value">4e-10</param> -->
+ <param name="value">4e12</param>
</struct>
<struct>
<param name="name">vy</param>
@@ -687,16 +685,22 @@
<param name="Shape">dikeShape</param>
<list name="variables">
<struct>
- <param name="name">vx</param>
- <param name="type">func</param>
- <param name="value">Velocity_Extension</param>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0.0</param>
</struct>
+ <struct>
+ <param name="name">vz</param>
+ <param name="type">double</param>
+ <param name="value">0.0</param>
+ </struct>
</list>
</struct>
</list>
</struct>
+<!-- ============== Temperature BCs =============== -->
<struct name="temperatureBCs">
<param name="type">CompositeVC</param>
<list name="vcList">
@@ -717,8 +721,8 @@
<list name="variables">
<struct>
<param name="name">temperature</param>
- <param name="type">double</param>
- <param name="value">1373</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
</struct>
</list>
</struct>
@@ -769,6 +773,7 @@
</list>
</struct>
+<!-- ============== Temperature ICs ===================== -->
<struct name="temperatureICs">
<param name="type">CompositeVC</param>
<list name="vcList">
@@ -785,31 +790,52 @@
</list>
</struct>
- <!-- These numbers have to be kept in sync with the top and bottom BC's. The
- parser can not read the variables and use them there, so they have
- to be specified in two places. -->
-
+<!-- ================= PARAMETERS ======================= -->
+ <list name="plugins">
+ <param>Underworld_EulerDeform</param>
+ <param>Underworld_VTKOutput</param>
+ </list>
+ <param name="maxTimeSteps">100</param>
+ <param name="dumpEvery">1</param>
+ <param name="outputPath">./output.rifting</param>
+ <param name="dim">3</param>
+ <param name="shadowDepth">1</param>
+ <param name="minX">0.0f</param>
+ <param name="minY">0.0f</param>
+ <param name="minZ">0.0f</param>
+ <param name="maxX">60000.0f</param>
+ <param name="maxY">16000.0f</param>
+ <param name="maxZ">6000.0f</param>
+ <param name="elementResI">240</param>
+ <param name="elementResJ">64</param>
+ <param name="elementResK">12</param>
+ <param name="allowUnbalancing">True</param>
+ <param name="buildNodeNeighbourTbl">True</param>
+ <param name="buildElementNeighbourTbl">True</param>
+ <param name="particlesPerCell">30</param>
+ <param name="seed">13</param>
+ <param name="gaussParticlesX">2</param>
+ <param name="gaussParticlesY">2</param>
+ <param name="gaussParticlesZ">2</param>
+ <param name="particlesPerCell">30</param>
<param name="TemperatureProfileTop">273</param>
- <!-- Assume thermal conductivity of 3 W/m K and heat flux at 0.029 W/m^2 -->
- <param name="TemperatureProfileLinearCoefficient">0.055</param>
-<!-- <param name="checkpointEvery">1</param> -->
- <param name="dtFactor">1.0</param>
+ <param name="TemperatureProfileLinearCoefficient">0.075</param>
- <param name="ExtensionFactor">1.4e13</param>
-<!-- <param name="ExtensionFactor">1.4e-13</param> -->
- <param name="ExtensionCentreX">20000</param>
+ <param name="TemperatureProfileMax">1473.0</param>
+ <param name="TemperatureProfileH0">7000.0</param>
+ <param name="TemperatureProfiledH">27000.0</param>
+ <param name="TemperatureProfileHm">16000.0</param>
+ <param name="ExtensionCentreX">30000.0</param>
- <param name="linearTolerance">1e-5</param>
- <param name="StepFunctionLowerValue">-1.4e-13</param>
- <param name="StepFunctionUpperValue">100000000</param>
- <param name="StepFunctionLowerOffset">-100000001</param>
- <param name="StepFunctionUpperOffset">-100000000</param>
+ <param name="dtFactor">0.5</param>
+ <param name="linearTolerance">1.0e-5</param>
+ <param name="nlinearTolerance">1.0e-2</param>
-
<param name="journal.info">True</param>
+ <param name="journal-level.info">2</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