[cig-commits] r4009 - long/3D/Gale/trunk
walter at geodynamics.org
walter at geodynamics.org
Wed Jul 12 11:30:42 PDT 2006
Author: walter
Date: 2006-07-12 11:30:41 -0700 (Wed, 12 Jul 2006)
New Revision: 4009
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/simple-extension.xml
Log:
r445 at earth: boo | 2006-07-12 11:29:28 -0700
Use updated input file from Luke
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:427
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:445
Modified: long/3D/Gale/trunk/simple-extension.xml
===================================================================
--- long/3D/Gale/trunk/simple-extension.xml 2006-07-11 20:38:41 UTC (rev 4008)
+++ long/3D/Gale/trunk/simple-extension.xml 2006-07-12 18:30:41 UTC (rev 4009)
@@ -1,13 +1,5 @@
<?xml version="1.0"?>
<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-
- <list name="plugins">
- <param>StG_FEM_StandardConditionFunctions</param>
- <param>StG_FEM_FrequentOutput</param>
- <param>StG_FEM_CPUTime</param>
- <param>Underworld_EulerDeform</param>
- </list>
-
<struct name="components">
<struct name="geometry">
<param name="Type">MeshGeometry</param>
@@ -34,7 +26,6 @@
<param name="sizeI">elementResI</param>
<param name="sizeJ">elementResJ</param>
<param name="sizeK">elementResK</param>
- <param name="dynamic">True</param>
</struct>
<struct name="bodyNodeLayout">
<param name="Type">BodyNL</param>
@@ -45,6 +36,7 @@
<param name="Type">HexaMD</param>
<param name="ElementLayout">elementLayout</param>
<param name="NodeLayout">bodyNodeLayout</param>
+ <param name="numPartitionedDims">numPartitionedDims</param>
</struct>
<struct name="meshLayout-constant">
<param name="Type">MeshLayout</param>
@@ -64,7 +56,6 @@
<param name="sizeJ">elementResJ</param>
<param name="sizeK">elementResK</param>
<param name="sizeShift">1</param>
- <param name="dynamic">True</param>
</struct>
<struct name="cornerNodeLayout">
<param name="Type">CornerNL</param>
@@ -75,6 +66,7 @@
<param name="Type">HexaMD</param>
<param name="ElementLayout">elementLayout</param>
<param name="NodeLayout">cornerNodeLayout</param>
+ <param name="numPartitionedDims">numPartitionedDims</param>
</struct>
<struct name="meshLayout-linear">
<param name="Type">MeshLayout</param>
@@ -169,11 +161,10 @@
<struct name="PressureField">
<param name="Type">FeVariable</param>
<param name="FEMesh">mesh-constant</param>
+ <param name="GeometryMesh">mesh-linear</param>
<param name="DofLayout">pressureDofLayout</param>
<param name="LinkedDofInfo">pressureLinkedDofs</param>
</struct>
-
- <!-- Cell layout. -->
<struct name="cellLayout">
<param name="Type">SingleCellLayout</param>
</struct>
@@ -204,7 +195,7 @@
</struct>
<struct name="timeIntegrator">
<param name="Type">TimeIntegrator</param>
- <param name="order">timeIntegratorOrder</param>
+ <param name="order">1</param>
<param name="simultaneous">t</param>
<param name="Context">context</param>
</struct>
@@ -212,30 +203,25 @@
<param name="Type">ElementCellLayout</param>
<param name="Mesh">mesh-linear</param>
</struct>
- <struct name="materialPointsParticleLayout">
- <param name="Type">RandomParticleLayout</param>
- <param name="seed">13</param>
- <param name="cellParticleCount">particlesPerCell</param>
- </struct>
<struct name="discreteVoronoi">
- <param name="Type">CellularAutomataVoronoi</param>
- <param name="resolution">10</param>
- <param name="diagonalNeighbours">true</param>
+ <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">DiscreteVoronoiWeights</param>
+ <param name="DiscreteVoronoi">discreteVoronoi</param>
</struct>
+ <struct name="localLayout">
+ <param name="Type">MappedParticleLayout</param>
+ </struct>
<struct name="picIntegrationPoints">
<param name="Type">IntegrationPointsSwarm</param>
<param name="CellLayout">elementCellLayout</param>
- <param name="ParticleLayout">materialPointsParticleLayout</param>
- <param name="WeightsCalculator">weights</param>
+ <param name="ParticleLayout">localLayout</param>
<param name="FiniteElement_Mesh">mesh-linear</param>
+ <param name="WeightsCalculator">weights</param>
<param name="TimeIntegrator">timeIntegrator</param>
- <!--param name="SplittingRoutine">splittingRoutine</param-->
- <!--param name="RemovalRoutine">removalRoutine</param-->
- <param name="EscapedRoutine">escapedRoutine</param>
<param name="IntegrationPointMapper">mapper</param>
</struct>
<struct name="mapper">
@@ -256,15 +242,24 @@
<param name="RemovalRoutine">removalRoutine</param>
<param name="EscapedRoutine">escapedRoutine</param>
</struct>
- <struct name="materialPointsAdvector">
+ <struct name="materialSwarmAdvector">
<param name="Type">SwarmAdvector</param>
<param name="Swarm">materialSwarm</param>
<param name="TimeIntegrator">timeIntegrator</param>
<param name="VelocityField">VelocityField</param>
<param name="PeriodicBCsManager">periodicBCsManager</param>
+ <param name="allowFallbackToFirstOrder">True</param>
</struct>
-
- <!-- More variables. -->
+ <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>
<param name="FeVariable">VelocityField</param>
@@ -277,19 +272,18 @@
<param name="Type">ForceVector</param>
<param name="FeVariable">VelocityField</param>
<param name="ExtraInfo">context</param>
- <param name="Swarm">gaussSwarm</param>
</struct>
<struct name="cont_force">
<param name="Type">ForceVector</param>
<param name="FeVariable">PressureField</param>
<param name="ExtraInfo">context</param>
- <param name="Swarm">gaussSwarm</param>
</struct>
<struct name="k_matrix">
<param name="Type">StiffnessMatrix</param>
<param name="RowVariable">VelocityField</param>
<param name="ColumnVariable">VelocityField</param>
<param name="RHS">mom_force</param>
+ <param name="allowZeroElementContributions">False</param>
</struct>
<struct name="constitutiveMatrix">
<param name="Type">ConstitutiveMatrixCartesian</param>
@@ -301,6 +295,7 @@
<param name="RowVariable">VelocityField</param>
<param name="ColumnVariable">PressureField</param>
<param name="RHS">cont_force</param>
+ <param name="allowZeroElementContributions">False</param>
</struct>
<struct name="gradientStiffnessMatrixTerm">
<param name="Type">GradientStiffnessMatrixTerm</param>
@@ -312,6 +307,7 @@
<param name="RowVariable">PressureField</param>
<param name="ColumnVariable">PressureField</param>
<param name="RHS">cont_force</param>
+ <param name="allowZeroElementContributions">True</param>
</struct>
<struct name="preconditionerTerm">
<param name="Type">UzawaPreconditionerTerm</param>
@@ -321,7 +317,6 @@
<struct name="uzawa">
<param name="Type">Stokes_SLE_UzawaSolver</param>
<param name="Preconditioner">preconditioner</param>
- <!--param name="tolerance">5.0e-4</param-->
<param name="tolerance">1.0e-5</param>
<param name="maxIterations">5000</param>
</struct>
@@ -332,6 +327,7 @@
<param name="StressTensorMatrix">k_matrix</param>
<param name="GradientMatrix">g_matrix</param>
<param name="DivergenceMatrix"></param>
+ <param name="CompressibilityMatrix">c_matrix</param>
<param name="VelocityVector">solutionVelocity</param>
<param name="PressureVector">solutionPressure</param>
<param name="ForceVector">mom_force</param>
@@ -340,70 +336,38 @@
<param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
<param name="nonLinearTolerance">nonLinearTolerance</param>
</struct>
- <struct name="c_matrix">
- <param name="Type">StiffnessMatrix</param>
- <param name="RowVariable">PressureField</param>
- <param name="ColumnVariable">PressureField</param>
- <param name="RHS">cont_force</param>
- </struct>
-
- <!-- Population control -->
- <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="escapedRoutine">
<param name="Type">EscapedRoutine</param>
<param name="idealParticleCount">0</param>
</struct>
-
- <!-- This is a box to cover the whole simulation -->
- <struct name="box">
+ <struct name="crustShape">
<param name="Type">Box</param>
<param name="startX">minX</param>
<param name="endX">maxX</param>
<param name="startY">minY</param>
- <param name="endY">maxY</param>
+ <param name="endY">0.35</param>
<param name="startZ">minZ</param>
<param name="endZ">maxZ</param>
</struct>
-
- <!-- This box just covers the the PDMS material -->
- <struct name="pdmsBox">
+ <struct name="weakZoneShape">
<param name="Type">Box</param>
- <param name="startX">pdmsMinX</param>
- <param name="endX">pdmsMaxX</param>
- <param name="startY">pdmsNinY</param>
- <param name="endY">pdmsMaxY</param>
- <param name="startZ">pdmsMinZ</param>
- <param name="endZ">pdmsMaxZ</param>
+ <param name="startX">0.5</param>
+ <param name="endX">1.5</param>
+ <param name="startY">minY</param>
+ <param name="endY">0.05</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
</struct>
-
- <!-- We subtract the PDMS from the box to make the sand's shape -->
- <struct name="sandShape">
+ <struct name="crustShape2">
<param name="Type">Intersection</param>
<list name="shapes">
- <param>box</param>
- <param>!pdmsBox</param>
+ <param>crustShape</param>
+ <param>!weakZoneShape</param>
</list>
</struct>
-
- <!-- Rheologies -->
- <struct name="StrainWeakening">
- <param name="Type">StrainWeakening</param>
- <param name="TimeIntegrator">timeIntegrator</param>
- <param name="MaterialPointsSwarm">materialSwarm</param>
- <param name="initialStrain">0.5</param>
- <param name="finalStrain">1.0</param>
- <param name="initialDamageFraction">0.0</param>
- <param name="initialDamageWavenumber">0.5</param>
- <param name="initialDamageFactor">0.5</param>
+ <struct name="crustViscosity">
+ <param name="Type">MaterialViscosity</param>
+ <param name="eta0">10.000</param>
</struct>
<struct name="director">
<param name="Type">Director</param>
@@ -415,50 +379,59 @@
<param name="initialDirectionZ">0.0</param>
<param name="dontUpdate">True</param>
</struct>
- <struct name="sandYielding">
+ <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="Director">director</param>
<param name="Context">context</param>
- <param name="StrainWeakening">StrainWeakening</param>
+ <param name="StrainWeakening">strainWeakening</param>
<param name="cohesion">10.0</param>
- <param name="cohesionAfterSoftening">10.0</param>
- <param name="frictionCoefficient">1.0/tan((36/180)*Pi)</param>
- <param name="frictionCoefficientAfterSoftening">1.0/tan((31/180)*Pi)</param>
+ <param name="cohesionAfterSoftening">0.0001</param>
+ <param name="frictionCoefficient">0.0</param>
+ <param name="frictionCoefficientAfterSoftening">0.0</param>
<param name="minimumYieldStress">0.00001</param>
<param name="ignoreOldOrientation">false</param>
</struct>
- <struct name="sandViscosity">
+ <struct name="mantleViscosity">
<param name="Type">MaterialViscosity</param>
- <param name="eta0">1.0</param>
+ <param name="eta0">1.000</param>
</struct>
- <struct name="sand">
+ <struct name="storeViscosity">
+ <param name="Type">StoreVisc</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+ <struct name="crust">
<param name="Type">RheologyMaterial</param>
- <param name="Shape">sandShape</param>
+ <param name="Shape">crustShape2</param>
<list name="Rheology">
- <param>sandViscosity</param>
- <!--param>sandYielding</param-->
+ <param>crustViscosity</param>
+ <param>yielding</param>
</list>
</struct>
- <struct name="pdmsViscosity">
- <param name="Type">MaterialViscosity</param>
- <param name="eta0">50000.0</param>
- </struct>
- <struct name="pdms">
+ <struct name="mantle">
<param name="Type">RheologyMaterial</param>
- <param name="Shape">pdmsBox</param>
- <param name="Rheology">pdmsViscosity</param>
+ <param name="Shape">weakZoneShape</param>
+ <param name="Rheology">mantleViscosity</param>
</struct>
-
<struct name="velocityRemesher">
<param name="Type">StripRemesher</param>
<param name="mesh">mesh-linear</param>
<param name="meshType">regular</param>
<list name="dims">
<param>true</param>
- <param>true</param>
+ <param>true</param>
</list>
</struct>
<struct name="pressureRemesher">
@@ -469,95 +442,82 @@
<param name="cellMesh">mesh-linear</param>
</struct>
</struct>
-
- <!-- Euler deformation parameters. -->
+ <list name="plugins">
+ <param>StG_FEM_FrequentOutput</param>
+ <param>Underworld_Vrms</param>
+ <param>StG_FEM_CPUTime</param>
+ <param>Underworld_EulerDeform</param>
+ <param>StG_FEM_StandardConditionFunctions</param>
+ <!--param>lucPlugin</param-->
+ </list>
<struct name="EulerDeform">
<list name="systems">
<struct>
+ <param name="mesh">mesh-constant</param>
<param name="remesher">pressureRemesher</param>
- <param name="velocityField">VelocityField</param>
- <param name="wrap">False</param>
- <list name="fields">
- <struct>
- <param name="field">PressureField</param>
- <param name="variable">pressure</param>
- </struct>
- </list>
+ <param name="velocityField">VelocityField</param>
+ <param name="wrap">False</param>
+ <list name="fields">
+ <struct>
+ <param name="field">PressureField</param>
+ <param name="variable">pressure</param>
+ </struct>
+ </list>
</struct>
<struct>
+ <param name="mesh">mesh-linear</param>
<param name="remesher">velocityRemesher</param>
- <param name="velocityField">VelocityField</param>
- <param name="wrap">True</param>
- <list name="fields">
- <struct>
- <param name="field">VelocityField</param>
- <param name="variable">velocity</param>
- </struct>
- </list>
+ <param name="velocityField">VelocityField</param>
+ <param name="wrap">True</param>
+ <list name="fields">
+ <struct>
+ <param name="field">VelocityField</param>
+ <param name="variable">velocity</param>
+ </struct>
+ </list>
</struct>
</list>
</struct>
-
- <!-- Overall simulation parameters. -->
- <param name="maxTimeSteps">200000</param>
+ <param name="maxTimeSteps">500</param>
<param name="outputPath">./output</param>
<param name="dumpEvery">1</param>
- <param name="saveEvery">1</param>
- <param name="nonLinearMaxIterations">500</param>
+ <param name="nonLinearMaxIterations">20</param>
<param name="journal.info.Context-verbose">True</param>
<param name="journal-level.info.Stg_ComponentFactory">2</param>
<param name="journal.debug">t</param>
<param name="journal-enable-branch.debug.StG_FEM">t</param>
<param name="journal-level-branch.debug.StG_FEM">2</param>
<param name="dim">2</param>
- <param name="elementResI">64</param>
- <param name="elementResJ">32</param>
- <param name="elementResK">16</param>
+ <param name="elementResI">32</param>
+ <param name="elementResJ">16</param>
+ <param name="elementResK">10</param>
<param name="shadowDepth">1</param>
<param name="minX">0.0</param>
<param name="minY">0.0</param>
<param name="minZ">0.0</param>
- <param name="maxX">4.0</param>
- <param name="maxY">2.0</param>
- <param name="maxZ">3.5</param>
- <param name="pdmsMinX">1.0</param>
- <param name="pdmsMinY">minY</param>
- <param name="pdmsMinZ">0</param>
- <param name="pdmsMaxX">3.0</param>
- <param name="pdmsMaxY">0.5</param>
- <param name="pdmsMaxZ">0.5</param>
+ <param name="maxX">2.0</param>
+ <param name="maxY">0.35</param>
+ <param name="maxZ">2.0</param>
<param name="buildNodeNeighbourTbl">True</param>
<param name="buildElementNeighbourTbl">True</param>
<param name="allowUnbalancing">True</param>
- <param name="particlesPerCell">40</param>
-
- <param name="StepFunctionOffset">2.0</param>
- <param name="StepFunctionValue">2.5</param>
-
-
- <!-- Boundary Conditions -->
+ <param name="particlesPerCell">30</param>
+ <param name="cellParticleCount">30</param>
+ <param name="seed">13</param>
<struct name="velocityBCs">
<param name="type">CompositeVC</param>
<list name="vcList">
- <!-- This is a split, moving, frictional boundary. -->
<struct>
<param name="type">WallVC</param>
<param name="wall">bottom</param>
<list name="variables">
- <struct>
- <param name="name">vx</param>
- <param name="type">func</param>
- <param name="value">StepFunction</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>
-
- <!-- Static, frictional left boundary -->
<struct>
<param name="type">WallVC</param>
<param name="wall">left</param>
@@ -565,39 +525,44 @@
<struct>
<param name="name">vx</param>
<param name="type">double</param>
- <param name="value">0.0</param>
+ <param name="value">-0.5</param>
</struct>
</list>
</struct>
-
- <!-- Moving right boundary -->
<struct>
<param name="type">WallVC</param>
- <param name="wall">right</param>
- <list name="variables">
- <struct>
- <param name="name">vx</param>
- <param name="type">double</param>
- <param name="value">2.5</param>
- </struct>
+ <param name="wall">right</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vx</param>
+ <param name="type">double</param>
+ <param name="value">0.5</param>
+ </struct>
</list>
</struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">front</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vz</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">back</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vz</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
</list>
</struct>
-
- <param name="gaussParticlesX">2</param>
- <param name="gaussParticlesY">2</param>
- <param name="gaussParticlesZ">2</param>
-
- <!-- Visualisation stuff -->
-<!-- <include>ParticleDensityVP.xml</include> -->
- <!--include>PressureVP.xml</include-->
-<!-- <include>./window.xml</include> -->
-<!-- <struct name="components" mergeType="merge"> -->
-<!-- <struct name="window" mergeType="merge"> -->
-<!-- <param name="Viewport">ParticleDensityVP</param> -->
-
-<!-- </struct> -->
-<!-- </struct> -->
-
+ <param name="interactive">t</param>
</StGermainData>
More information about the cig-commits
mailing list