[cig-commits] r8032 - in long/3D/Gale/trunk: . input/examples
walter at geodynamics.org
walter at geodynamics.org
Wed Sep 26 13:12:12 PDT 2007
Author: walter
Date: 2007-09-26 13:12:12 -0700 (Wed, 26 Sep 2007)
New Revision: 8032
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/examples/tibet.xml
Log:
r1939 at earth: boo | 2007-09-26 02:33:43 -0700
Add temperature to tibet
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1938
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1939
Modified: long/3D/Gale/trunk/input/examples/tibet.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet.xml 2007-09-26 20:12:09 UTC (rev 8031)
+++ long/3D/Gale/trunk/input/examples/tibet.xml 2007-09-26 20:12:12 UTC (rev 8032)
@@ -298,6 +298,7 @@
<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.8</param>
</struct>
@@ -354,7 +355,7 @@
<param name="x2">350000</param>
<param name="x3">650000</param>
<param name="x4">700000</param>
- <param name="height">5000</param>
+ <param name="height">3030.3030303</param>
</struct>
<struct name="boxShape">
@@ -388,8 +389,8 @@
<columnDefinition name = "y" type="double"/>
0 68000
300000 68000
- 350000 53000
- 650000 53000
+ 350000 51030.3030303
+ 650000 51030.3030303
700000 68000
1000000 68000
1000000 100000
@@ -401,6 +402,14 @@
<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>backgroundViscosity</param>
<param>storeViscosity</param>
@@ -420,6 +429,14 @@
<param name="Type">RheologyMaterial</param>
<param name="Shape">mantleShape</param>
<param name="density">3300</param>
+ <param name="alpha">3.0e-5</param>
+ <param name="diffusivity">1.0e-6</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">8.48484848485e-13</param>
+ <param name="lambda">0.0</param>
+ </struct>
+ </list>
<list name="Rheology">
<param>backgroundViscosity</param>
<param>storeViscosity</param>
@@ -427,6 +444,85 @@
</list>
</struct>
+
+ <!-- Temperature components -->
+ <struct name="temperature">
+ <param name="Type">MeshVariable</param>
+ <param name="Rank">Scalar</param>
+ <param name="DataType">Double</param>
+ <param name="mesh">mesh-linear</param>
+ </struct>
+ <struct name="temperatureBCs">
+ <param name="Type">CompositeVC</param>
+ <param name="Data">mesh-linear</param>
+ </struct>
+ <struct name="temperatureICs">
+ <param name="Type">CompositeVC</param>
+ <param name="Data">mesh-linear</param>
+ </struct>
+ <struct name="temperatureDofLayout">
+ <param name="Type">DofLayout</param>
+ <param name="mesh">mesh-linear</param>
+ <list name="BaseVariables">
+ <param>temperature</param>
+ </list>
+ </struct>
+ <struct name="TemperatureField">
+ <param name="Type">FeVariable</param>
+ <param name="FEMesh">mesh-linear</param>
+ <param name="DofLayout">temperatureDofLayout</param>
+ <param name="BC">temperatureBCs</param>
+ <param name="IC">temperatureICs</param>
+ <param name="LinkedDofInfo">temperatureLinkedDofs</param>
+ </struct>
+ <struct name="TemperatureGradientsField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">Gradient</param>
+ <param name="FeVariable">TemperatureField</param>
+ </struct>
+
+ <!-- Energy Equation -->
+
+ <struct name="residual">
+ <param name="Type">ForceVector</param>
+ <param name="FeVariable">TemperatureField</param>
+ </struct>
+ <struct name="massMatrix">
+ <param name="Type">ForceVector</param>
+ <param name="FeVariable">TemperatureField</param>
+ </struct>
+ <struct name="predictorMulticorrector">
+ <param name="Type">AdvDiffMulticorrector</param>
+ </struct>
+ <struct name="EnergyEqn">
+ <param name="Type">AdvectionDiffusionSLE</param>
+ <param name="SLE_Solver">predictorMulticorrector</param>
+ <param name="Context">context</param>
+ <param name="PhiField">TemperatureField</param>
+ <param name="Residual">residual</param>
+ <param name="MassMatrix">massMatrix</param>
+ <param name="courantFactor">0.25</param>
+ </struct>
+ <struct name="lumpedMassMatrixForceTerm">
+ <param name="Type">LumpedMassMatrixForceTerm</param>
+ <param name="Swarm">gaussSwarm</param>
+ <param name="ForceVector">massMatrix</param>
+ </struct>
+ <struct name="defaultResidualForceTerm">
+ <param name="Type">AdvDiffResidualForceTerm</param>
+ <param name="Swarm">gaussSwarm</param>
+ <param name="ForceVector">residual</param>
+ <param name="ExtraInfo">EnergyEqn</param>
+ <param name="VelocityField">VelocityField</param>
+ <param name="defaultDiffusivity">defaultDiffusivity</param>
+ <param name="UpwindXiFunction">Exact</param>
+ </struct>
+ <struct name="internalHeatingTerm">
+ <param name="Type">RadiogenicHeatingTerm</param>
+ <param name="ForceVector">residual</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ </struct>
+
</struct>
<list name="plugins">
<param>Underworld_EulerDeform</param>
@@ -544,5 +640,109 @@
</struct>
</list>
</struct>
+
+
+ <struct name="temperatureBCs">
+ <param name="type">CompositeVC</param>
+ <list name="vcList">
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">top</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">double</param>
+ <!-- Remember to keep in sync with TemperatureProfileTop -->
+ <param name="value">273</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">bottom</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">double</param>
+ <!-- Remember to keep in sync with TemperatureProfile parameters -->
+ <param name="value">1332.998763</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">left</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">right</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">front</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
+ </struct>
+ </list>
+ </struct>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">back</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
+ </struct>
+ </list>
+ </struct>
+ </list>
+ </struct>
+
+ <struct name="temperatureICs">
+ <param name="type">CompositeVC</param>
+ <list name="vcList">
+ <struct>
+ <param name="type">AllNodesVC</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">func</param>
+ <param name="value">TemperatureProfile</param>
+ </struct>
+ </list>
+ </struct>
+ </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. -->
+
+ <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.0096667</param>
+ <!-- Heating rate = 2.8e-6 W/m^3, h_r=10 km, thermal conductivity 3 W/m K -->
+ <param name="TemperatureProfileExponentialCoefficient1">93.333</param>
+ <!-- 1/h_r = 10^-4 -->
+ <param name="TemperatureProfileExponentialCoefficient2">1.0e-4</param>
+
<param name="checkpointEvery">1</param>
</StGermainData>
More information about the cig-commits
mailing list