[cig-commits] commit: Add dike example
Mercurial
hg at geodynamics.org
Sun Jul 25 09:09:31 PDT 2010
changeset: 303:fefcce820596
tag: tip
user: Walter Landry <wlandry at caltech.edu>
date: Sun Jul 25 09:06:59 2010 -0700
files: input/examples/dike.xml
description:
Add dike example
diff -r e50f25eaf548 -r fefcce820596 input/examples/dike.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/input/examples/dike.xml Sun Jul 25 09:06:59 2010 -0700
@@ -0,0 +1,978 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+ <import>
+ <toolbox>Underworld</toolbox>
+ </import>
+ <plugins>
+ <struct>
+ <param name="Type">Underworld_EulerDeform</param>
+ <param name="Context">context</param>
+ </struct>
+ <struct>
+ <param name="Type">Underworld_VTKOutput</param>
+ <param name="Context">context</param>
+ </struct>
+ <struct>
+ <param name="Type">StgFEM_StandardConditionFunctions</param>
+ <param name="Context">context</param>
+ </struct>
+ </plugins>
+ <struct name="components">
+ <struct name="context">
+ <param name="Type">UnderworldContext</param>
+ </struct>
+ <struct name="mesh-linear">
+ <param name="Type">FeMesh</param>
+ <param name="elementType">linear</param>
+ </struct>
+ <struct name="linearMesh-generator">
+ <param name="Type">CartesianGenerator</param>
+ <param name="mesh">mesh-linear</param>
+ <param name="dim">dim</param>
+ <param name="shadowDepth">shadowDepth</param>
+ <param name="regular">False</param>
+ <list name="size">
+ <param>elementResI</param>
+ <param>elementResJ</param>
+ <param>elementResK</param>
+ </list>
+ <list name="minCoord">
+ <param>minX</param>
+ <param>minY</param>
+ <param>minZ</param>
+ </list>
+ <list name="maxCoord">
+ <param>maxX</param>
+ <param>maxAirY</param>
+ <param>maxZ</param>
+ </list>
+ </struct>
+ <struct name="velocity">
+ <param name="Type">MeshVariable</param>
+ <param name="mesh">mesh-linear</param>
+ <param name="Rank">Vector</param>
+ <param name="DataType">Double</param>
+ <param name="VectorComponentCount">dim</param>
+ <list name="names">
+ <param>vx</param>
+ <param>vy</param>
+ <param>vz</param>
+ </list>
+ </struct>
+ <struct name="velocityBCs">
+ <param name="Type">CompositeVC</param>
+ <param name="Data">mesh-linear</param>
+ </struct>
+ <struct name="velocityICs">
+ <param name="Type">CompositeVC</param>
+ <param name="Data">mesh-linear</param>
+ </struct>
+ <struct name="velocityDofLayout">
+ <param name="Type">DofLayout</param>
+ <param name="mesh">mesh-linear</param>
+ <param name="BaseVariableCount">dim</param>
+ <list name="BaseVariables">
+ <param>vx</param>
+ <param>vy</param>
+ <param>vz</param>
+ </list>
+ </struct>
+ <struct name="VelocityField">
+ <param name="Type">FeVariable</param>
+ <param name="FEMesh">mesh-linear</param>
+ <param name="DofLayout">velocityDofLayout</param>
+ <param name="BC">velocityBCs</param>
+ <param name="IC">velocityICs</param>
+ <param name="LinkedDofInfo">velocityLinkedDofs</param>
+ </struct>
+ <struct name="VelocityMagnitudeField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">Magnitude</param>
+ <param name="FeVariable">VelocityField</param>
+ </struct>
+ <struct name="VelocityGradientsField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">Gradient</param>
+ <param name="FeVariable">VelocityField</param>
+ </struct>
+ <struct name="VelocityGradientsInvariantField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">TensorInvariant</param>
+ <param name="FeVariable">VelocityGradientsField</param>
+ </struct>
+ <struct name="StrainRateField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">TensorSymmetricPart</param>
+ <param name="FeVariable">VelocityGradientsField</param>
+ </struct>
+ <struct name="VorticityField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">TensorAntisymmetricPart</param>
+ <param name="FeVariable">VelocityGradientsField</param>
+ </struct>
+ <struct name="StrainRateInvariantField">
+ <param name="Type">OperatorFeVariable</param>
+ <param name="Operator">SymmetricTensor_Invariant</param>
+ <param name="FeVariable">StrainRateField</param>
+ </struct>
+ <struct name="pressure">
+ <param name="Type">MeshVariable</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-linear</param>
+ <list name="BaseVariables">
+ <param>pressure</param>
+ </list>
+ </struct>
+ <struct name="PressureField">
+ <param name="Type">FeVariable</param>
+ <param name="FEMesh">mesh-linear</param>
+ <param name="DofLayout">pressureDofLayout</param>
+ <param name="LinkedDofInfo">pressureLinkedDofs</param>
+ </struct>
+ <struct name="StressField">
+ <param name="Type">StressField</param>
+ <param name="StrainRateField">StrainRateField</param>
+ <param name="Context">context</param>
+ <param name="ConstitutiveMatrix">constitutiveMatrix</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="Mesh">mesh-linear</param>
+ <param name="IC">stressICs</param>
+ </struct>
+ <struct name="ViscosityField">
+ <param name="Type">ViscosityField</param>
+ <param name="Context">context</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="Mesh">mesh-linear</param>
+ <param name="ConstitutiveMatrix">constitutiveMatrix</param>
+ </struct>
+ <struct name="cellLayout">
+ <param name="Type">SingleCellLayout</param>
+ </struct>
+ <struct name="particleLayout">
+ <param name="Type">GaussParticleLayout</param>
+ </struct>
+ <struct name="gaussSwarm">
+ <param name="Type">IntegrationPointsSwarm</param>
+ <param name="CellLayout">cellLayout</param>
+ <param name="ParticleLayout">particleLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="IntegrationPointMapper">gaussMapper</param>
+ </struct>
+ <struct name="gaussMapper">
+ <param name="Type">GaussMapper</param>
+ <param name="IntegrationPointsSwarm">gaussSwarm</param>
+ <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+ </struct>
+ <struct name="backgroundLayout">
+ <param name="Type">BackgroundParticleLayout</param>
+ </struct>
+ <struct name="gaussMSwarmMovementHandler">
+ <param name="Type">ParticleMovementHandler</param>
+ </struct>
+ <struct name="gaussMaterialSwarm">
+ <param name="Type">MaterialPointsSwarm</param>
+ <param name="CellLayout">cellLayout</param>
+ <param name="ParticleLayout">backgroundLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ <list name="ParticleCommHandlers">
+ <param>gaussMSwarmMovementHandler</param>
+ </list>
+ </struct>
+ <struct name="timeIntegrator">
+ <param name="Type">TimeIntegrator</param>
+ <param name="order">1</param>
+ <param name="simultaneous">t</param>
+ <param name="Context">context</param>
+ </struct>
+ <struct name="elementCellLayout">
+ <param name="Type">ElementCellLayout</param>
+ <param name="Mesh">mesh-linear</param>
+ </struct>
+ <struct name="weights">
+ <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>
+ <param name="Inflow">True</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">localLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ <param name="WeightsCalculator">weights</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="IntegrationPointMapper">mapper</param>
+ </struct>
+ <struct name="mapper">
+ <param name="Type">CoincidentMapper</param>
+ <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+ <struct name="materialSwarmParticleLayout">
+ <param name="Type">MeshParticleLayout</param>
+ <param name="cellParticleCount">particlesPerCell</param>
+ <param name="mesh">mesh-linear</param>
+ </struct>
+ <struct name="pMovementHandler">
+ <param name="Type">ParticleMovementHandler</param>
+ </struct>
+ <struct name="pShadowSync">
+ <param name="Type">ParticleShadowSync</param>
+ </struct>
+ <struct name="materialSwarm">
+ <param name="Type">MaterialPointsSwarm</param>
+ <param name="CellLayout">elementCellLayout</param>
+ <param name="ParticleLayout">materialSwarmParticleLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ <list name="ParticleCommHandlers">
+ <param>pMovementHandler</param>
+ <param>pShadowSync</param>
+ </list>
+ <param name="SplittingRoutine">splittingRoutine</param>
+ <param name="RemovalRoutine">removalRoutine</param>
+ <param name="EscapedRoutine">escapedRoutine</param>
+ </struct>
+ <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>
+ <struct name="solutionVelocity">
+ <param name="Type">SolutionVector</param>
+ <param name="FeVariable">VelocityField</param>
+ </struct>
+ <struct name="solutionPressure">
+ <param name="Type">SolutionVector</param>
+ <param name="FeVariable">PressureField</param>
+ </struct>
+ <struct name="mom_force">
+ <param name="Type">ForceVector</param>
+ <param name="FeVariable">VelocityField</param>
+ <param name="ExtraInfo">context</param>
+ </struct>
+ <struct name="cont_force">
+ <param name="Type">ForceVector</param>
+ <param name="FeVariable">PressureField</param>
+ <param name="ExtraInfo">context</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>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="StiffnessMatrix">k_matrix</param>
+ </struct>
+ <struct name="g_matrix">
+ <param name="Type">StiffnessMatrix</param>
+ <param name="RowVariable">VelocityField</param>
+ <param name="ColumnVariable">PressureField</param>
+ <param name="RHS">mom_force</param>
+ <param name="transposeRHS">cont_force</param>
+ <param name="allowZeroElementContributions">False</param>
+ </struct>
+ <struct name="gradientStiffnessMatrixTerm">
+ <param name="Type">GradientStiffnessMatrixTerm</param>
+ <param name="Swarm">gaussSwarm</param>
+ <param name="StiffnessMatrix">g_matrix</param>
+ </struct>
+ <struct name="preconditioner">
+ <param name="Type">StiffnessMatrix</param>
+ <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>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="StiffnessMatrix">preconditioner</param>
+ </struct>
+ <struct name="uzawa">
+ <param name="Type">Stokes_SLE_UzawaSolver</param>
+ <param name="velocitySolver">matrixSolver</param>
+ <param name="Preconditioner">preconditioner</param>
+ <param name="tolerance">linearTolerance</param>
+ <param name="monitor">false</param>
+ <param name="maxIterations">5000</param>
+ <param name="minIterations">1</param>
+ </struct>
+ <struct name="stokesEqn">
+ <param name="Type">Stokes_SLE</param>
+ <param name="SLE_Solver">uzawa</param>
+ <param name="Context">context</param>
+ <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>
+ <param name="ContinuityForceVector">cont_force</param>
+ <param name="killNonConvergent">false</param>
+ <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+ <param name="nonLinearMinIterations">nonLinearMinIterations</param>
+ <param name="nonLinearTolerance">nonLinearTolerance</param>
+ <param name="makeConvergenceFile">false</param>
+ <param name="isNonLinear">true</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>
+ <param name="allowZeroElementContributions">True</param>
+ </struct>
+ <struct name="mixedStabiliser">
+ <param name="Type">MixedStabiliserTerm</param>
+ <param name="Swarm">gaussSwarm</param>
+ <param name="picSwarm">picIntegrationPoints</param>
+ <param name="storeVisc">storeViscosity</param>
+ <param name="StiffnessMatrix">c_matrix</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="HydrostaticTerm">hydrostaticTerm</param>
+ <param name="gravity">9.81</param>
+ </struct>
+ <struct name="background">
+ <param name="Type">Everywhere</param>
+ </struct>
+ <struct name="escapedRoutine">
+ <param name="Type">EscapedRoutine</param>
+ <param name="idealParticleCount">0</param>
+ </struct>
+ <struct name="velocityRemesher">
+ <param name="Type">RegularRemesherCmpt</param>
+ <param name="mesh">mesh-linear</param>
+ <list name="remeshDims">
+ <param>0</param>
+ <param>1</param>
+ <param>2</param>
+ </list>
+ </struct>
+ <struct name="storeViscosity">
+ <param name="Type">StoreVisc</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+ <struct name="storeStress">
+ <param name="Type">StoreStress</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ </struct>
+
+<!--============== END "BASIC COMPONENTS" =================-->
+ <struct name="hydrostaticTerm">
+ <param name="Type">HydrostaticTerm</param>
+ <param name="upperDensity">2300</param>
+ <param name="upperAlpha">3.0e-5</param>
+ <param name="lowerDensity">2300</param>
+ <param name="lowerAlpha">3.0e-5</param>
+ <param name="height">maxY</param>
+<!-- <param name="materialBoundary">maxY</param>-->
+ <param name="T_0">TemperatureProfileTop</param>
+ <param name="linearCoefficient">TemperatureProfileLinearCoefficient</param>
+ <param name="exponentialCoefficient1">TemperatureProfileExponentialCoefficient1</param>
+ <param name="exponentialCoefficient2">TemperatureProfileExponentialCoefficient2</param>
+ <param name="gravity">9.81</param>
+ </struct>
+
+ <struct name="backgroundShape">
+ <param name="Type">Everywhere</param>
+ </struct>
+
+ <struct name="boxShape">
+ <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="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ <!--<param name="beta">0.785398163398</param>-->
+ </struct>
+
+<!-- ============== Crust ================ -->
+<!-- <struct name="crustShape"> -->
+<!-- <param name="Type">Box</param> -->
+<!-- <param name="startX">minX</param> -->
+<!-- <param name="endX">maxX</param> -->
+<!-- <param name="startY">10000</param> -->
+<!-- <param name="endY">maxY</param> -->
+<!-- <param name="startZ">minZ</param> -->
+<!-- <param name="endZ">maxZ</param> -->
+<!-- </struct> -->
+<!-- <struct name="crustViscosity"> -->
+<!-- <param name="Type">FrankKamenetskii</param> -->
+<!-- <param name="TemperatureField">TemperatureField</param> -->
+<!-- ref. visc is 1e22 so this times eta0 gives -->
+<!-- viscosity=1e25, at T=273 & viscosity=1e20 at T=1473 -->
+<!-- <param name="eta0">2.31739465e4</param> -->
+<!-- <param name="theta">0.011512925465</param> -->
+<!-- </struct> -->
+<!-- <struct name="crust"> -->
+<!-- <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.0e16</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> -->
+
+<!-- ============== Strain Weakening ================ -->
+ <struct name="strainWeakening">
+ <param name="Type">StrainWeakening</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="initialSofteningStrain">0.0</param>
+ <param name="finalSofteningStrain">0.1</param>
+ <param name="initialDamageFraction">0.9</param>
+ <param name="initialDamageWavenumber">0.0</param>
+ <param name="initialDamageFactor">0.9</param>
+ <param name="strainLimitedShape">dikeShape</param> -->
+ <param name="initialStrainShape">initStrainBox</param> -->
+ <param name="healingRate">6.0</param>
+ </struct>
+ <struct name="initStrainBox">
+ <param name="Type">Box</param>
+ <param name="startX">29500</param>
+ <param name="endX">45000.0</param>
+ <param name="startY">0.0</param>
+ <param name="endY">maxY</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ </struct>
+ <struct name="yielding">
+ <param name="Type">DruckerPrager</param>
+ <param name="PressureField">PressureField</param>
+ <param name="VelocityGradientsField">VelocityGradientsField</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="Context">context</param>
+ <param name="StrainRateField">StrainRateField</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>
+ <param name="HydrostaticTerm">hydrostaticTerm</param>
+ </struct>
+
+<!-- ============== Mantle ================ -->
+ <struct name="mantleShape">
+ <param name="Type">Intersection</param>
+ <list name="shapes">
+ <param>boxShape</param>
+<!-- <param>!crustShape</param> -->
+ </list>
+ </struct>
+
+ <struct name="airShape">
+ <param name="Type">Intersection</param>
+ <list name="shapes">
+ <param>backgroundShape</param>
+ <param>!mantleShape</param>
+ </list>
+ </struct>
+
+ <struct name="airViscosity">
+ <param name="Type">MaterialViscosity</param>
+ <param name="eta0">1e-4</param>
+ </struct>
+
+ <!-- Define the compressibility function -->
+ <struct name="compresibility">
+ <param name="Type">Compressible</param>
+ <param name="GeometryMesh">mesh-linear</param>
+ <param name="StiffnessMatrix">c_matrix</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="oneOnLambda">10.0</param> <!-- this is the 'compressibility factor, 1/lambda, which must be set between 0 and 10. The higher, the more compressible -->
+ </struct>
+
+ <struct name="air">
+ <param name="Type">RheologyMaterial</param>
+ <param name="Shape">airShape</param>
+ <param name="density">0</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">0.0</param>
+ <param name="lambda">0.0</param>
+ </struct>
+ </list>
+ <list name="Rheology">
+ <param>airViscosity</param>
+ <param>storeViscosity</param>
+ <param>storeStress</param>
+ </list>
+ <param name="Compressible">compresibility</param>-->
+ </struct>
+
+ <struct name="nonnewViscosity">
+ <param name="Type">NonNewtonian</param>
+ <param name="StrainRateInvariantField">StrainRateInvariantField</param>
+ <param name="TemperatureField">TemperatureField</param>
+ <param name="n">3.0</param>
+ <param name="T_0">1473</param>
+ <param name="T_melt">1480</param>
+ <param name="A">1.0e-12</param>
+ <param name="refStrainRate">1.0e22</param>
+ <param name="minViscosity">1.0e-2</param>
+ <param name="maxViscosity">1.0e3</param>
+ </struct>
+
+ <struct name="mantleViscous">
+ <param name="Type">RheologyMaterial</param>
+ <param name="Shape">mantleShape</param>
+ <param name="density">2300</param> <!--3300-1030 kg/m^3-->
+ <param name="alpha">3.0e-5</param>
+ <param name="diffusivity">1.0e16</param> <!--1.0e-6*1.0e22-->
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">8.48484848485e-13</param>
+ <param name="lambda">0.0</param>
+ </struct>
+ </list>
+ <list name="Rheology">
+ <param>nonnewViscosity</param>
+ <param>yielding</param>
+ <param>storeViscosity</param>
+ <param>storeStress</param>
+ </list>
+ </struct>
+
+<!-- ============== Diking Zone ================ -->
+ <struct name="dikeShape">
+ <param name="Type">Box</param>
+ <param name="startX">29000</param>
+ <param name="endX">31000</param>
+ <param name="startY">0</param>
+ <!-- <param name="startY">8000</param> -->
+ <param name="endY">16000</param>
+ <param name="startZ">minZ</param>
+ <param name="endZ">maxZ</param>
+ </struct>
+ <struct name="nonDikeShape">
+ <param name="Type">Intersection</param>
+ <list name="shapes">
+ <param>boxShape</param>
+ <param>!dikeShape</param>
+ </list>
+ </struct>
+
+
+ <struct name="divergenceForce">
+ <param name="Type">DivergenceForce</param>
+ <param name="ForceVector">cont_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="DomainShape">dikeShape</param>
+ <param name="GeometryMesh">mesh-linear</param>
+ <param name="force_type">double</param>
+ <param name="force_value">-2.0e9</param>
+ </struct>
+
+<!-- ============== Stress BC on bottom ===================== -->
+ <struct name="stressBC">
+ <param name="Type">StressBC</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="wall">top</param>
+<!-- <param name="wall">bottom</param> -->
+ <param name="y_type">HydrostaticTerm</param>
+ <param name="y_value">hydrostaticTerm</param> <!-- 16000*3300*(1-3e-5*1473/2)*9.81 -->
+<!-- <param name="y_type">double</param>-->
+<!-- <param name="y_value">5.04e+8</param>--> <!-- 16000*3300*(1-3e-5*1473/2)*9.81=5.07e+8 -->
+ </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>
+
+<!-- ====the end of component section==== -->
+
+<!-- ============== Mesh ==================== -->
+ <struct name="EulerDeform">
+ <list name="systems">
+ <struct>
+ <param name="mesh">mesh-linear</param>
+ <param name="remesher">velocityRemesher</param>
+ <param name="velocityField">VelocityField</param>
+ <param name="wrapTop">False</param>
+ <param name="staticTop">True</param>-->
+ <param name="staticBottom">True</param> -->
+ <param name="staticFront">True</param>
+ <param name="staticBack">True</param>
+ <param name="staticLeft">True</param>
+<!-- <param name="staticLeftTop">True</param>-->
+<!-- <param name="floatLeftTop">True</param>-->
+<!-- <param name="xLeftCoord">minX</param>-->
+ <param name="staticLeftBottom">True</param> -->
+ <param name="staticRight">True</param>
+<!-- <param name="staticRightTop">True</param>-->
+<!-- <param name="floatRightTop">True</param>-->
+<!-- <param name="xRightCoord">maxX</param>-->
+ <param name="staticRightBottom">True</param> -->
+ <list name="fields">
+ <struct>
+ <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>
+ </struct>
+
+ <!-- ============== Velocity BCs ==================== -->
+ <struct name="velocityBCs">
+ <param name="type">CompositeVC</param>
+ <list name="vcList">
+ <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>
+ <struct>
+ <param name="type">WallVC</param>
+ <param name="wall">left</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vx</param>
+ <param name="type">double</param>
+ <param name="value">-4e12</param> <!-- -4e-10*1e25 -->
+ </struct>
+ <struct>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
+
+ <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">4e12</param>
+ </struct>
+ <struct>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
+ </list>
+ </struct>
+
+<!-- ============== Temperature BCs =============== -->
+ <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>
+ <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">func</param>
+ <param name="value">TemperatureProfile</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>
+ <struct>
+ <param name="type">MeshShapeVC</param>
+ <param name="Shape">nonDikeShape</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">MeshShapeVC</param>
+ <param name="Shape">dikeShape</param>
+ <list name="variables">
+ <struct>
+ <param name="name">temperature</param>
+ <param name="type">double</param>
+ <param name="value">1500</param>
+ </struct>
+ </list>
+ </struct>
+ </list>
+ </struct>
+
+<!-- ============== Temperature ICs ===================== -->
+ <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>
+
+<!-- ================= PARAMETERS ======================= -->
+ <list name="FieldVariablesToCheckpoint">
+ <param>VelocityField</param>
+ <param>PressureField</param>
+ </list>
+ <param name="maxTimeSteps">15000</param>
+ <param name="outputEvery">1</param>
+ <param name="dumpEvery">2</param>
+ <param name="outputPath">output_dike</param>
+ <param name="dim">2</param>
+ <param name="shadowDepth">1</param>
+ <param name="minX">0.0f</param>
+ <param name="minY">0.0f</param>
+ <param name="minZ">0.0f</param>
+ <param name="maxX">60000.0f</param>
+ <param name="maxY">16000.0f</param>
+ <param name="maxAirY">18000.0f</param>
+ <param name="maxZ">2500.0f</param>
+ <param name="elementResI">120</param>
+ <param name="elementResJ">36</param>
+ <param name="elementResK">5</param>
+ <param name="allowUnbalancing">True</param>
+ <param name="buildNodeNeighbourTbl">True</param>
+ <param name="buildElementNeighbourTbl">True</param>
+ <param name="particlesPerCell">30</param>
+ <param name="seed">13</param>
+ <param name="gaussParticlesX">2</param>
+ <param name="gaussParticlesY">2</param>
+ <param name="gaussParticlesZ">2</param>
+ <param name="defaultDiffisivity">1.0e16</param> <!--1.0e-6*1.0e25-->
+
+<!-- === PARAMETERS for temperature profile rotation ============ -->
+ <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.10909</param> <!--1200K/depth[1200K]=0.102176 0.10909-->
+<!-- Heating rate = 2.8e-6 W/m^3, h_r=10 km, thermal conductivity 3 W/m K -->
+ <param name="TemperatureProfileExponentialCoefficient1">0</param>
+ <!-- 1/h_r = 10^-4 -->
+ <param name="TemperatureProfileExponentialCoefficient2">0</param>
+
+
+
+ <param name="TemperatureProfileMax">1473.0</param>
+ <param name="TemperatureProfileH0">7000.0</param>
+ <param name="TemperatureProfiledH">27000.0</param>
+ <param name="TemperatureProfileHm">16000.0</param>
+ <param name="ExtensionCentreX">30000</param>
+
+ <!--<param name="SolidBodyRotationCentreX">15000</param>-->
+ <!--<param name="SolidBodyRotationCentreY">0</param>-->
+ <!--<param name="SolidBodyRotationCentreZ">0</param>-->
+ <!--<param name="TemperatureRotationOmega">0</param>--><!--counterclockwise rotation from x axis-->
+<!-- === PARAMETERS for temperature profile rotation end========== -->
+<!-- <param name="GaussianHeight">5.25e8</param> -->
+<!-- <param name="GaussianWidth">5000</param> -->
+<!-- <param name="GaussianCenter">30000</param> -->
+<!-- <param name="GaussianDim">0</param> -->
+
+ <!--<param name="dt">0</param>-->
+ <param name="dtFactor">0.5</param>
+ <param name="linearTolerance">1.0e-5</param>
+ <param name="nonLinearTolerance">1.0e-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