[cig-commits] r12584 - in long/3D/Gale/trunk: . input/benchmarks
walter at geodynamics.org
walter at geodynamics.org
Thu Aug 7 15:56:36 PDT 2008
Author: walter
Date: 2008-08-07 15:56:36 -0700 (Thu, 07 Aug 2008)
New Revision: 12584
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/benchmarks/isostasy.xml
Log:
r2298 at earth: boo | 2008-08-07 15:54:59 -0700
Make isostasy.xml use StressBC
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2297
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2298
Modified: long/3D/Gale/trunk/input/benchmarks/isostasy.xml
===================================================================
--- long/3D/Gale/trunk/input/benchmarks/isostasy.xml 2008-08-07 22:56:32 UTC (rev 12583)
+++ long/3D/Gale/trunk/input/benchmarks/isostasy.xml 2008-08-07 22:56:36 UTC (rev 12584)
@@ -122,6 +122,27 @@
<param name="DofLayout">pressureDofLayout</param>
<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>
@@ -130,7 +151,7 @@
<param name="Swarm">picIntegrationPoints</param>
<param name="Mesh">mesh-linear</param>
<param name="IC">stressICs</param>
- </struct>
+ </struct>
<struct name="ViscosityField">
<param name="Type">ViscosityField</param>
<param name="Context">context</param>
@@ -176,14 +197,16 @@
<param name="Type">ElementCellLayout</param>
<param name="Mesh">mesh-linear</param>
</struct>
- <struct name="discreteVoronoi">
- <param name="Type">CellularAutomataVoronoi</param>
- <param name="resolution">10</param>
- <param name="diagonalNeighbours">true</param>
- </struct>
<struct name="weights">
- <param name="Type">DiscreteVoronoiWeights</param>
- <param name="DiscreteVoronoi">discreteVoronoi</param>
+ <param name="Type">PCDVC</param>
+ <param name="resolutionX">10</param>
+ <param name="resolutionY">10</param>
+ <param name="resolutionZ">10</param>
+ <param name="lowerT">0.6</param>
+ <param name="upperT">25</param>
+ <param name="maxDeletions">3</param>
+ <param name="maxSplits">3</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
</struct>
<struct name="localLayout">
<param name="Type">MappedParticleLayout</param>
@@ -224,15 +247,9 @@
<param name="PeriodicBCsManager">periodicBCsManager</param>
<param name="allowFallbackToFirstOrder">True</param>
</struct>
- <struct name="removalRoutine">
- <param name="Type">DiscreteVoronoiRemove</param>
- <param name="idealParticleCount">particlesPerCell</param>
- <param name="DiscreteVoronoi">discreteVoronoi</param>
- </struct>
<struct name="splittingRoutine">
<param name="Type">ReseedSplitting</param>
<param name="idealParticleCount">particlesPerCell</param>
- <param name="DiscreteVoronoi">discreteVoronoi</param>
</struct>
<struct name="solutionVelocity">
<param name="Type">SolutionVector</param>
@@ -292,8 +309,8 @@
<param name="Type">Stokes_SLE_UzawaSolver</param>
<!-- <param name="velocitySolver">matrixSolver</param> -->
<param name="Preconditioner">preconditioner</param>
- <param name="tolerance">1.0e-5</param>
- <param name="maxIterations">5000</param>
+ <param name="tolerance">linearTolerance</param>
+ <param name="maxIterations">maxUzawaIterations</param>
</struct>
<struct name="stokesEqn">
<param name="Type">Stokes_SLE</param>
@@ -316,36 +333,17 @@
<param name="Type">BuoyancyForceTerm</param>
<param name="ForceVector">mom_force</param>
<param name="Swarm">picIntegrationPoints</param>
- <param name="gravity">gravity</param>
+ <param name="gravity">1.0</param>
</struct>
<struct name="backgroundShape">
<param name="Type">Everywhere</param>
</struct>
- <struct name="pdmsShape">
- <param name="Type">Box</param>
- <param name="startX">0.5</param>
- <param name="endX">1.5</param>
- <param name="startY">0.0</param>
- <param name="endY">0.05</param>
- <param name="startZ">minZ</param>
- <param name="endZ">maxZ</param>
- </struct>
- <struct name="crustShape">
- <param name="Type">Intersection</param>
- <list name="shapes">
- <param>backgroundShape</param>
- <param>!pdmsShape</param>
- </list>
- </struct>
- <struct name="crustViscosity">
- <param name="Type">MaterialViscosity</param>
- <param name="eta0">10.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="finalSofteningStrain">0.1</param>
<param name="initialDamageFraction">0.0</param>
<param name="initialDamageWavenumber">0.5</param>
<param name="initialDamageFactor">0.5</param>
@@ -358,11 +356,12 @@
<param name="MaterialPointsSwarm">materialSwarm</param>
<param name="Context">context</param>
<param name="StrainWeakening">strainWeakening</param>
- <param name="cohesion">10.0</param>
- <param name="cohesionAfterSoftening">0.0001</param>
+ <param name="cohesion">100000.0</param>
+ <param name="cohesionAfterSoftening">100000</param>
<param name="frictionCoefficient">0.0</param>
<param name="frictionCoefficientAfterSoftening">0.0</param>
</struct>
+
<struct name="pdmsViscosity">
<param name="Type">MaterialViscosity</param>
<param name="eta0">1.0</param>
@@ -375,26 +374,17 @@
<param name="Type">StoreStress</param>
<param name="MaterialPointsSwarm">materialSwarm</param>
</struct>
- <struct name="crust">
+ <struct name="pdms">
<param name="Type">RheologyMaterial</param>
- <param name="Shape">crustShape</param>
+ <param name="Shape">backgroundShape</param>
<param name="density">1.0</param>
<list name="Rheology">
- <param>crustViscosity</param>
+ <param>pdmsViscosity</param>
<param>yielding</param>
<param>storeViscosity</param>
<param>storeStress</param>
</list>
</struct>
- <struct name="pdms">
- <param name="Type">RheologyMaterial</param>
- <param name="Shape">pdmsShape</param>
- <param name="density">1.0</param>
- <list name="Rheology">
- <param>pdmsViscosity</param>
- <param>storeViscosity</param>
- </list>
- </struct>
<struct name="escapedRoutine">
<param name="Type">EscapedRoutine</param>
<param name="idealParticleCount">0</param>
@@ -425,24 +415,25 @@
<param name="Type"> DistortedOpGenerator </param>
<param name="fineVariable"> VelocityField </param>
</struct>
- <struct name="conditionFunctions">
- <param name="Type">StgFEM_StandardConditionFunctions</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>
- <param name="y_value">325</param>
+ <param name="y_value">1</param>
</struct>
+
+ <struct name="conditionFunctions">
+ <param name="Type">StgFEM_StandardConditionFunctions</param>
+ </struct>
</struct>
<list name="plugins">
<param>Underworld_EulerDeform</param>
- <param>Underworld_DumpSwarm</param>
<param>Underworld_VTKOutput</param>
</list>
- <param name="maxTimeSteps">2</param>
+ <param name="maxTimeSteps">10</param>
<param name="dumpEvery">1</param>
<param name="outputPath">./output.isostasy</param>
<param name="dim">2</param>
@@ -450,10 +441,10 @@
<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">1.0f</param>
+ <param name="maxY">1.0f</param>
<param name="maxZ">1.0f</param>
- <param name="elementResI">64</param>
+ <param name="elementResI">16</param>
<param name="elementResJ">16</param>
<param name="elementResK">8</param>
<param name="allowUnbalancing">True</param>
@@ -493,6 +484,7 @@
</struct>
</list>
</struct>
+
<struct name="velocityBCs">
<param name="type">CompositeVC</param>
<list name="vcList">
@@ -525,12 +517,12 @@
<struct>
<param name="name">vx</param>
<param name="type">double</param>
- <param name="value">0.0</param>
+ <param name="value">0</param>
</struct>
<struct>
<param name="name">vy</param>
<param name="type">double</param>
- <param name="value">0.0</param>
+ <param name="value">0</param>
</struct>
</list>
</struct>
@@ -546,35 +538,21 @@
<struct>
<param name="name">vy</param>
<param name="type">double</param>
- <param name="value">0.0</param>
- </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</param>
</struct>
</list>
</struct>
</list>
</struct>
- <param name="StepFunctionLowerOffset">0.8</param>
- <param name="StepFunctionUpperOffset">1.2</param>
- <param name="StepFunctionValue">1.0</param>
- <param name="StepFunctionDim">0</param>
- <param name="StepFunctionLessThan">False</param>
- <param name="checkpointEvery">1</param>
- <param name="mgLevels">7</param>
- <param name="gravity">1000</param>
+<!-- <param name="checkpointEvery">1</param> -->
+ <param name="mgLevels">2</param>
+ <param name="dtFactor">1.0</param>
+ <param name="pressureEquilibriumPressure">1.0</param>
+ <param name="pressureEquilibriumDensity">1.0</param>
+ <param name="pressureEquilibriumGravity">1.0</param>
+ <param name="linearTolerance">1.0e-7</param>
+ <param name="maxUzawaIterations">5000</param>
+ <param name="nonLinearTolerance">1.0e-5</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