[CIG-LONG] Thermally driven counterflow
Magali Billen
mibillen at ucdavis.edu
Thu Oct 14 13:04:20 PDT 2010
Walter - can you explain whether there was an actual bug in the code
that could
affect other people's modeling?
- Magali
On Oct 14, 2010, at 11:29 AM, Walter Landry wrote:
> Walter Landry <walter at geodynamics.org> wrote:
>> Karen Paczkowski <karen.paczkowski at yale.edu> wrote:
>>> Hi,
>>>
>>> Thank you for all the help. I am unfortunately still getting a very
>>> unexpected flow field. The only velocities I have set are:
>>>
>>> Vx = 0 and Vy = -Vslab on the left side
>>> Vy=0 across the top
>>>
>>> My input file is producing a flow field that flows in the opposite
>>> direction as Vslab (up and to the right) and has a magnitude 2-3
>>> orders of magnitude larger than Vslab. This was not occurring prior
>>> to adding the thermal field, so I am assuming it has to do with
>>> that,
>>> but I'm really not sure what is causing this exactly. I have
>>> attached
>>> a picture of the velocity field and a copy of the input file. It is
>>> written for Gale 1.5.0. I would really appreciate any advice.
>>
>> The problem is that there is no maximum temperature in the
>> HydrostaticTerm. I am on vacation until September, so I will not be
>> able to get to it until then.
>
> Sorry for the long delay, but I finally implemented this. You can get
> the code from the repository. I am attaching a new input file and a
> picture of the velocity and temperature. Please let me know if you
> have any more questions.
>
> Cheers,
> Walter Landry
> walter at geodynamics.org
> <?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>maxY</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>
> </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="mesh">mesh-linear</param>
> <param name="cellParticleCount">particlesPerCell</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="nonLinearTolerance">nonLinearTolerance</param>
> <param name="makeConvergenceFile">false</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="gravity">gravity</param>
> <param name="HydrostaticTerm">hydrostaticTerm</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="bottomBoundaryShape">
> <param name="Type">Box</param>
> <param name="startX">minX</param>
> <param name="endX">maxX</param>
> <param name="startY">minY</param>
> <param name="endY">10.00</param>
> <param name="startZ">minZ</param>
> <param name="endZ">maxZ</param>
> </struct>
>
> <!-- <struct name="fixedCornerShape"> -->
> <!-- <param name="Type">Box</param> -->
> <!-- <param name="startX">1999000</param> -->
> <!-- <param name="endX">maxX</param> -->
> <!-- <param name="startY">439000</param> -->
> <!-- <param name="endY">maxY</param> -->
> <!-- <param name="startZ">minZ</param> -->
> <!-- <param name="endZ">maxZ</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>
> </struct>
> <struct name="backgroundViscosity">
> <param name="Type">MaterialViscosity</param>
> <param name="eta0">1.0</param>
> </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>
> <struct name="viscous">
> <param name="Type">RheologyMaterial</param>
> <param name="Shape">boxShape</param>
> <param name="density">3000.01</param>
> <param name="alpha">3e-5</param>
> <param name="diffusivity">1e15</param>
> <list name="heatingElements">
> <struct>
> <param name="Q">0</param>
> <param name="lambda">0</param>
> </struct>
> </list>
> <list name="Rheology">
> <param>backgroundViscosity</param>
> <param>storeViscosity</param>
> <param>storeStress</param>
> </list>
> </struct>
>
> <struct name="hydrostaticTerm">
> <param name="Type">HydrostaticTerm</param>
> <param name="upperDensity">3000</param>
> <param name="upperAlpha">3e-5</param>
> <param name="height">maxY</param>
> <param name="T_0">TemperatureProfileTop</param>
> <param name="T_max">TemperatureProfileMax</param>
> <param name="T_max_depth">100000</param>
> <param
> name="linearCoefficient">TemperatureProfileLinearCoefficient</param>
> <param name="TemperatureProfileMax">TemperatureProfileMax</param>
> <param name="gravity">gravity</param>
> </struct>
>
> <struct name="stressBCTop">
> <param name="Type">StressBC</param>
> <param name="ForceVector">mom_force</param>
> <param name="Swarm">picIntegrationPoints</param>
> <param name="wall">top</param>
> <param name="y_type">HydrostaticTerm</param>
> <param name="y_value">hydrostaticTerm</param>
> </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="FieldVariablesToCheckpoint">
> <param>VelocityField</param>
> <param>PressureField</param>
> </list>
> <param name="maxTimeSteps">60</param>
> <param name="outputEvery">10</param>
> <param name="dumpEvery">10</param>
> <param name="outputPath">output_3</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">1320000.0f</param>
> <param name="maxY">440000.0f</param>
> <param name="maxZ">0.5f</param>
> <param name="elementResI">128</param>
> <param name="elementResJ">64</param>
> <param name="elementResK">4</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="particlesPerCell">30</param>
> <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="staticRight">True</param>
> <param name="staticLeft">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>
> <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">0.0</param>
> </struct>
> <struct>
> <param name="name">vy</param>
> <param name="type">double</param>
> <param name="value">-1.0e12</param>
> </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>
> <!-- <struct>-->
> <!-- <param name="type">WallVC</param> -->
> <!-- <param name="wall">right</param>-->
> <!-- <list name="variables">-->
> <!-- <struct> -->
> <!-- <param name="name">vy</param>-->
> <!-- <param name="type">double</param>-->
> <!-- <param name="value">0.0</param>-->
> <!-- </struct>-->
> <!-- </list>-->
> <!-- </struct>-->
> <!-- <struct> -->
> <!-- <param name="name">vx</param> -->
> <!-- <param name="type">double</param> -->
> <!-- <param name="value">0.0</param> -->
> <!-- </struct> -->
> <!-- </list> -->
> <!-- </struct> -->
> <!-- <struct> -->
> <!-- <param name="type">MeshShapeVC</param> -->
> <!-- <param name="Shape">fixedCornerShape</param> -->
> <!-- <list name="variables"> -->
> <!-- <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">double</param>
> <param name="value">0.0</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>
>
>
> <!-- ============== 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>
>
>
> <param name="TemperatureProfileTop">273.0</param>
> <param name="TemperatureProfileLinearCoefficient">0.0135</param>
> <param name="TemperatureProfileExponentialCoefficient1">0</param>
> <param name="TemperatureProfileExponentialCoefficient2">0</param>
> <param name="TemperatureProfileMax">1623.0</param>
> <param name="TemperatureProfileH0">100000.0</param>
>
>
> <!-- <param name="checkpointEvery">1</param>-->
> <param name="dtFactor">0.1</param>
> <param name="checkpointEvery">10</param>
> <param name="gravity">10.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>
>
> </StGermainData>
> <vertical_slab.png>_______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
-----------------------------
Associate Professor, U.C. Davis
Department of Geology/KeckCAVES
2129 Earth & Physical Sciences Bldg
Davis, CA 95616
-----------------
mibillen at ucdavis.edu
(530) 754-5696
--------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20101014/7c5e0a41/attachment-0001.htm
More information about the CIG-LONG
mailing list