[cig-commits] commit: Some fixes to wars model
Mercurial
hg at geodynamics.org
Mon Jul 26 13:07:20 PDT 2010
changeset: 307:76d3ba1044f8
user: Walter Landry <wlandry at caltech.edu>
date: Mon Jul 26 13:02:15 2010 -0700
files: input/examples/wars.xml
description:
Some fixes to wars model
diff -r 25aa4928d042 -r 76d3ba1044f8 input/examples/wars.xml
--- a/input/examples/wars.xml Mon Jul 26 12:59:20 2010 -0700
+++ b/input/examples/wars.xml Mon Jul 26 13:02:15 2010 -0700
@@ -488,7 +488,7 @@
<param name="A">13.3220431236</param> <!-- 6.18354466027e-8*(1e20)^(1/2.4) -->
<!-- <param name="refStrainRate">1.5e-16</param> -->
<param name="refStrainRate">1.5e4</param>
- <!-- <param name="minViscosity">1.0e-2</param> -->
+ <param name="minViscosity">1.0e3</param>
<!-- <param name="maxViscosity">1.0e20</param> -->
</struct>
@@ -548,11 +548,11 @@
<param name="TemperatureField">TemperatureField</param>
<param name="n">4.5</param>
<param name="T_0">59927.7978339</param>
- <param name="A">3.78647900946e-6</param>
+ <!-- <param name="A">3.78647900946e-6</param> -->
<param name="A">0.105361027689</param>
<!-- <param name="refStrainRate">1.5e-16</param> -->
<param name="refStrainRate">1.5e4</param>
- <!-- <param name="minViscosity">1.0e-2</param> -->
+ <param name="minViscosity">1.0e3</param>
<!-- <param name="maxViscosity">1.0e20</param> -->
</struct>
@@ -636,11 +636,12 @@
<param name="frictionCoefficient">0.427954752593</param>
<param name="frictionCoefficientAfterSoftening">0.427954752593</param>
<param name="HydrostaticTerm">hydrostaticTerm</param>
+ <param name="minimumViscosity">minimumViscosity</param>
</struct>
-<!-- ============== Stress BC on bottom ===================== -->
- <struct name="stressBC">
+<!-- ============== Stress BC on top and bottom ===================== -->
+ <struct name="topStressBC">
<param name="Type">StressBC</param>
<param name="ForceVector">mom_force</param>
<param name="Swarm">picIntegrationPoints</param>
@@ -648,6 +649,17 @@
<param name="y_type">HydrostaticTerm</param>
<param name="y_value">hydrostaticTerm</param>
</struct>
+
+ <!-- <struct name="bottomStressBC"> -->
+ <!-- <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">HydrostaticTerm</param> -->
+ <!-- <param name="y_value">hydrostaticTerm</param> -->
+ <!-- <param name="bottomDensity">3139.0821</param> -->
+ <!-- </struct> -->
+ <!-- 3139.0821 = 3300 * (1- 3.1e-5 * 1573) -->
<!-- ============== Temperature Components ===================== -->
<struct name="temperature">
@@ -766,11 +778,11 @@
<param name="type">double</param>
<param name="value">0</param>
</struct>
- <struct>
- <param name="name">vy</param>
- <param name="type">double</param>
- <param name="value">0</param>
- </struct>
+ <!-- <struct> -->
+ <!-- <param name="name">vy</param> -->
+ <!-- <param name="type">double</param> -->
+ <!-- <param name="value">0</param> -->
+ <!-- </struct> -->
</list>
</struct>
@@ -807,6 +819,18 @@
<!-- </struct> -->
<!-- </list> -->
<!-- </struct> -->
+
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">top</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
</list>
@@ -863,7 +887,7 @@
<param>VelocityField</param>
<param>PressureField</param>
</list>
- <param name="maxTimeSteps">150</param>
+ <param name="maxTimeSteps">10000</param>
<param name="outputEvery">1</param>
<param name="dumpEvery">2</param>
<param name="outputPath">output_wars</param>
@@ -903,9 +927,10 @@
<param name="dtFactor">0.5</param>
<param name="linearTolerance">1.0e-5</param>
<param name="nonLinearTolerance">1.0e-2</param>
+ <param name="minimumViscosity">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