[cig-commits] r13423 - in long/3D/Gale/trunk: . input/examples
walter at geodynamics.org
walter at geodynamics.org
Sat Nov 29 22:50:41 PST 2008
Author: walter
Date: 2008-11-29 22:50:41 -0800 (Sat, 29 Nov 2008)
New Revision: 13423
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/examples/tibet.xml
Log:
r2374 at dante: boo | 2008-11-20 14:31:54 -0800
unrealistic hacks to get tibet working with linear pressure elements
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2372
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2374
Modified: long/3D/Gale/trunk/input/examples/tibet.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet.xml 2008-11-30 00:29:41 UTC (rev 13422)
+++ long/3D/Gale/trunk/input/examples/tibet.xml 2008-11-30 06:50:41 UTC (rev 13423)
@@ -4,15 +4,6 @@
<struct name="conditionFunctions">
<param name="Type">StgFEM_StandardConditionFunctions</param>
</struct>
- <struct name="mesh-constant">
- <param name="Type">FeMesh</param>
- <param name="elementType">constant</param>
- </struct>
- <struct name="constantMesh-generator">
- <param name="Type">C0Generator</param>
- <param name="mesh">mesh-constant</param>
- <param name="elementMesh">mesh-linear</param>
- </struct>
<struct name="mesh-linear">
<param name="Type">FeMesh</param>
<param name="elementType">linear</param>
@@ -108,42 +99,23 @@
</struct>
<struct name="pressure">
<param name="Type">MeshVariable</param>
- <param name="mesh">mesh-constant</param>
+ <param name="mesh">mesh-linear</param>
<param name="Rank">Scalar</param>
<param name="DataType">Double</param>
</struct>
<struct name="pressureDofLayout">
<param name="Type">DofLayout</param>
- <param name="mesh">mesh-constant</param>
+ <param name="mesh">mesh-linear</param>
<list name="BaseVariables">
<param>pressure</param>
</list>
</struct>
<struct name="PressureField">
<param name="Type">FeVariable</param>
- <param name="FEMesh">mesh-constant</param>
+ <param name="FEMesh">mesh-linear</param>
<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>
@@ -330,13 +302,13 @@
<param name="nonLinearTolerance">1e-2</param>
<param name="makeConvergenceFile">false</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="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>
@@ -354,14 +326,21 @@
<param>true</param>
</list>
</struct>
- <struct name="pressureRemesher">
- <param name="Type">CellRemesher</param>
- <param name="mesh">mesh-constant</param>
- <param name="meshType">regular</param>
- <param name="dim">3</param>
- <param name="cellMesh">mesh-linear</param>
- </struct>
+ <element type="struct" name="c_matrix">
+ <element type="param" name="Type">StiffnessMatrix</element>
+ <element type="param" name="RowVariable">PressureField</element>
+ <element type="param" name="ColumnVariable">PressureField</element>
+ <element type="param" name="RHS">cont_force</element>
+ <element type="param" name="allowZeroElementContributions">True</element>
+ </element>
+ <element type="struct" name="mixedStabiliser">
+ <element type="param" name="Type">MixedStabiliserTerm</element>
+ <element type="param" name="Swarm">gaussSwarm</element>
+ <element type="param" name="materialSwarm">materialSwarm</element>
+ <element type="param" name="StiffnessMatrix">c_matrix</element>
+ </element>
+
<struct name="bottomBoundaryShape">
<param name="Type">Box</param>
<param name="startX">minX</param>
@@ -382,21 +361,21 @@
<!-- <param name="y_value">3080340000</param> -->
<!-- </struct> -->
- <struct name="surfaceAdaptor">
- <param name="Type">SurfaceAdaptor</param>
- <param name="mesh">mesh-linear</param>
- <param name="sourceGenerator">linearMesh-generator</param>
- <param name="surfaceType">plateau</param>
- <param name="x1">300000</param>
- <param name="x2">350000</param>
- <param name="x3">650000</param>
- <param name="x4">700000</param>
- <param name="z1">200000</param>
- <param name="z2">250000</param>
- <param name="z3">750000</param>
- <param name="z4">800000</param>
- <param name="height">3030.3030303</param>
- </struct>
+<!-- <struct name="surfaceAdaptor"> -->
+<!-- <param name="Type">SurfaceAdaptor</param> -->
+<!-- <param name="mesh">mesh-linear</param> -->
+<!-- <param name="sourceGenerator">linearMesh-generator</param> -->
+<!-- <param name="surfaceType">plateau</param> -->
+<!-- <param name="x1">300000</param> -->
+<!-- <param name="x2">350000</param> -->
+<!-- <param name="x3">650000</param> -->
+<!-- <param name="x4">700000</param> -->
+<!-- <param name="z1">200000</param> -->
+<!-- <param name="z2">250000</param> -->
+<!-- <param name="z3">750000</param> -->
+<!-- <param name="z4">800000</param> -->
+<!-- <param name="height">3030.3030303</param> -->
+<!-- </struct> -->
<struct name="boxShape">
<param name="Type">Box</param>
@@ -421,89 +400,91 @@
</struct>
<!-- Crust -->
- <struct name="crustShape">
- <param name="Type">PolygonShape</param>
- <list name="verticies">
- <asciidata>
- <columnDefinition name = "x" type="double"/>
- <columnDefinition name = "y" type="double"/>
- 0 68000
- 300000 68000
- 350000 51030.3030303
- 650000 51030.3030303
- 700000 68000
- 1000000 68000
- 1000000 100000
- 0 100000
- </asciidata>
- </list>
- </struct>
+<!-- <struct name="crustShape"> -->
+<!-- <param name="Type">PolygonShape</param> -->
+<!-- <list name="verticies"> -->
+<!-- <asciidata> -->
+<!-- <columnDefinition name = "x" type="double"/> -->
+<!-- <columnDefinition name = "y" type="double"/> -->
+<!-- 0 68000 -->
+<!-- 300000 68000 -->
+<!-- 350000 51030.3030303 -->
+<!-- 650000 51030.3030303 -->
+<!-- 700000 68000 -->
+<!-- 1000000 68000 -->
+<!-- 1000000 100000 -->
+<!-- 0 100000 -->
+<!-- </asciidata> -->
+<!-- </list> -->
+<!-- </struct> -->
- <struct name="crustViscosity">
- <param name="Type">FrankKamenetskii</param>
- <param name="TemperatureField">TemperatureField</param>
+<!-- <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.31739465e26</param>
- <param name="theta">0.011512925465</param>
- </struct>
+<!-- <param name="eta0">2.31739465e26</param> -->
+<!-- <param name="theta">0.011512925465</param> -->
+<!-- </struct> -->
- <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>
- </struct>
- <struct name="yielding">
- <param name="Type">MohrCoulomb</param>
- <param name="PressureField">NodalPressureField</param>
- <param name="VelocityGradientsField">VelocityGradientsField</param>
- <param name="MaterialPointsSwarm">materialSwarm</param>
- <param name="Context">context</param>
- <param name="StrainWeakening">strainWeakening</param>
- <param name="cohesion">44e6</param>
- <param name="cohesionAfterSoftening">4e6</param>
- <param name="frictionCoefficient">0.57735026919</param>
- <param name="frictionCoefficientAfterSoftening">0.57735026919</param>
- </struct>
+<!-- <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> -->
+<!-- </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="Context">context</param> -->
+<!-- <param name="StrainWeakening">strainWeakening</param> -->
+<!-- <param name="cohesion">44e6</param> -->
+<!-- <param name="cohesionAfterSoftening">4e6</param> -->
+<!-- <param name="frictionCoefficient">0.57735026919</param> -->
+<!-- <param name="frictionCoefficientAfterSoftening">0.57735026919</param> -->
+<!-- </struct> -->
- <struct name="crustViscous">
- <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="crustViscous"> -->
+<!-- <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> -->
<!-- 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.31739465e26</param>
- <param name="theta">0.011512925465</param>
+ <param name="Type">MaterialViscosity</param>
+ <param name="eta0">1.0</param>
+<!-- <param name="Type">FrankKamenetskii</param> -->
+<!-- <param name="TemperatureField">TemperatureField</param> -->
+<!-- <param name="eta0">2.31739465e26</param> -->
+<!-- <param name="theta">0.011512925465</param> -->
</struct>
<struct name="mantleViscous">
@@ -637,17 +618,6 @@
<struct name="EulerDeform">
<list name="systems">
<struct>
- <param name="mesh">mesh-constant</param>
- <param name="remesher">pressureRemesher</param>
- <param name="velocityField">VelocityField</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>
@@ -658,6 +628,10 @@
<param name="field">VelocityField</param>
<param name="variable">velocity</param>
</struct>
+ <struct>
+ <param name="field">PressureField</param>
+ <param name="variable">pressure</param>
+ </struct>
</list>
</struct>
</list>
@@ -832,8 +806,8 @@
<!-- <param name="checkpointEvery">1</param> -->
<param name="dtFactor">1.0</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="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