[cig-commits] r7473 - in long/3D/Gale/trunk: . documentation
documentation/images input/cookbook
walter at geodynamics.org
walter at geodynamics.org
Mon Jun 25 03:06:15 PDT 2007
Author: walter
Date: 2007-06-25 03:06:14 -0700 (Mon, 25 Jun 2007)
New Revision: 7473
Added:
long/3D/Gale/trunk/documentation/images/Paraview_normal_stress.png
long/3D/Gale/trunk/input/cookbook/viscous_normal_stress.xml
Removed:
long/3D/Gale/trunk/documentation/images/Paraview_viscous_open_boundary.png
long/3D/Gale/trunk/input/cookbook/viscous_open_boundary.xml
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/documentation/gale.lyx
long/3D/Gale/trunk/documentation/images/Paraview_thermal.png
long/3D/Gale/trunk/input/cookbook/thermal.xml
Log:
r1846 at earth: boo | 2007-06-25 03:01:34 -0700
Add in radiogenic to the manual and thermal example. Change open boundary to normal stress.
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1842
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1846
Modified: long/3D/Gale/trunk/documentation/gale.lyx
===================================================================
--- long/3D/Gale/trunk/documentation/gale.lyx 2007-06-25 09:18:26 UTC (rev 7472)
+++ long/3D/Gale/trunk/documentation/gale.lyx 2007-06-25 10:06:14 UTC (rev 7473)
@@ -1,4 +1,4 @@
-#LyX 1.4.3 created this file. For more info see http://www.lyx.org/
+#LyX 1.4.4 created this file. For more info see http://www.lyx.org/
\lyxformat 245
\begin_document
\begin_header
@@ -3345,12 +3345,8 @@
\end_layout
-\begin_layout Standard
-
-\end_layout
-
\begin_layout Section
-Viscous Material in Extension with Open Boundaries
+Viscous Material in Extension with Normal Stress Boundaries
\end_layout
\begin_layout Standard
@@ -3359,11 +3355,14 @@
\end_inset
- to use an open boundary on the bottom, instead of specifying the velocity.
- A stress boundary condition simulates the effect of material below the
- material pushing up, supporting the material in the box.
- This will allow material to flow in or out of the simulation as material
- is thinned out or piled up.
+ to use a stress boundary normal to the bottom surface, instead of specifying
+ the velocity.
+ A normal stress boundary condition simulates the effect of material below
+ the material pushing up, supporting the material in the box.
+ Then, when material piles up, material will overcome the stress boundary
+ and flow out of the simulation.
+ Conversely, if material is thinned out, the stress boundary will push new
+ material into the simulation.
This kind of boundary is often more relevant for geological simulations.
\end_layout
@@ -3686,7 +3685,7 @@
\begin_layout Standard
A worked example is at
\family typewriter
-input/cookbook/viscous_open_boundary.xml
+input/cookbook/viscous_normal_stress.xml
\family default
.
A visualization of the strain rate invariant and velocity (see Section
@@ -3696,7 +3695,7 @@
\end_inset
) is shown in Figure
-\begin_inset LatexCommand \ref{fig:Strain-viscous-open-boundary}
+\begin_inset LatexCommand \ref{fig:Strain-viscous-normal-stress}
\end_inset
@@ -3723,12 +3722,12 @@
\end_layout
\begin_layout Caption
-\begin_inset LatexCommand \label{fig:Strain-viscous-open-boundary}
+\begin_inset LatexCommand \label{fig:Strain-viscous-normal-stress}
\end_inset
Strain rate invariant and velocities of viscous material in extension with
- an open boundary
+ a normal stress boundary
\end_layout
\end_inset
@@ -4595,8 +4594,7 @@
\end_inset
- and heats it on the bottom.
-
+, heats it on the bottom, and adds in radiogenic heating throughout.
\end_layout
\begin_layout Enumerate
@@ -4975,8 +4973,8 @@
\end_deeper
\begin_layout Enumerate
-Specify the background material's coefficient of thermal expansivity and
- diffusivity by adding after
+Specify the background material's thermal expansivity, thermal diffusivity,
+ radiogenic heating rate, and radiogenic decay time scale by adding after
\end_layout
\begin_deeper
@@ -5008,6 +5006,30 @@
<param name="diffusivity">1.0</param>
\end_layout
+\begin_layout LyX-Code
+ <list name="heatingElements">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Q">1.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="lambda">1.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
\begin_layout Standard
For the sphere, after the lines
\end_layout
@@ -5040,8 +5062,32 @@
<param name="diffusivity">10.0</param>
\end_layout
+\begin_layout LyX-Code
+ <list name="heatingElements">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Q">1000.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="lambda">10.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
\begin_layout Standard
-This makes the sphere more expansive and conductive.
+This makes the sphere more expansive, conductive, and radioactive.
\end_layout
\end_deeper
@@ -5114,7 +5160,7 @@
\end_layout
\begin_layout LyX-Code
- <param name="value">0.0001</param>
+ <param name="value">0.01</param>
\end_layout
\end_deeper
@@ -7789,6 +7835,30 @@
</struct>
\end_layout
+\begin_layout LyX-Code
+<struct name="internalHeatingTerm">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Type">RadiogenicHeatingTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="ForceVector">residual</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Swarm">picIntegrationPoints</param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
\begin_layout Standard
You should never need to modify these components.
\end_layout
@@ -7832,7 +7902,12 @@
\end_inset
-)
+) and radiogenic heating (see Section
+\begin_inset LatexCommand \ref{sec:Materials}
+
+\end_inset
+
+).
\end_layout
\begin_layout Standard
@@ -9722,16 +9797,103 @@
\end_layout
\begin_layout Standard
-You can also specify a density, a coefficient of thermal expansivity (
+You can specify a density, a coefficient of thermal expansivity (
\begin_inset Formula $\alpha$
\end_inset
-), and a thermal diffusivity for each material.
+), and a thermal diffusivity, for each material.
\end_layout
\begin_layout Standard
+The density and expansivity are used by the BuoyancyForceTerm component
+ (see Section
+\begin_inset LatexCommand \ref{sub:BouyancyForceTerm}
+
+\end_inset
+
+) to create buoyancy forces.
+ The diffusivity is used by the temperature solver (see Section
+\begin_inset LatexCommand \ref{sec:Temperature-components}
+
+\end_inset
+
+).
+\end_layout
+
+\begin_layout Standard
+You can also specify multiple radiogenic heating rates (Q) and radiogenic
+ timescales (
+\begin_inset Formula $\lambda$
+\end_inset
+
+).
+ This simulates the action of multiple radioactive materials with different
+ half-lives.
+ To enable this, you must provide a list of Q's and
+\begin_inset Formula $\lambda$
+\end_inset
+
+'s.
+ For example, to specify two different radioactive species, add something
+ like
+\end_layout
+
+\begin_layout LyX-Code
+<list name="heatingElements">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Q">1.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="lambda">1.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Q">2.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="lambda">2.0</param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+</list>
+\end_layout
+
+\begin_layout Standard
+At a given time, each radioactive element will generate
+\end_layout
+
+\begin_layout Standard
+\begin_inset Formula \[
+Qe^{-\lambda t}\]
+
+\end_inset
+
+ units of energy.
+\end_layout
+
+\begin_layout Standard
\begin_inset Tabular
-<lyxtabular version="3" rows="4" columns="2">
+<lyxtabular version="3" rows="6" columns="2">
<features>
<column alignment="center" valignment="top" leftline="true" width="0">
<column alignment="center" valignment="top" leftline="true" rightline="true" width="0">
@@ -9785,7 +9947,7 @@
\end_inset
</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
\begin_inset Text
\begin_layout Standard
@@ -9805,7 +9967,7 @@
\end_inset
</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
\begin_inset Text
\begin_layout Standard
@@ -9815,25 +9977,51 @@
\end_inset
</cell>
</row>
-</lyxtabular>
+<row bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+\begin_layout Standard
+Q
+\end_layout
+
\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+\begin_layout Standard
+0
+\end_layout
-\newline
-The density and expansivity are used by the BuoyancyForceTerm component
- (see Section
-\begin_inset LatexCommand \ref{sub:BouyancyForceTerm}
+\end_inset
+</cell>
+</row>
+<row bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+\begin_layout Standard
+lambda
+\end_layout
+
\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
-) to create buoyancy forces.
- The diffusivity is used by the temperature solver (see Section
-\begin_inset LatexCommand \ref{sec:Temperature-components}
+\begin_layout Standard
+0
+\end_layout
\end_inset
+</cell>
+</row>
+</lyxtabular>
-).
+\end_inset
+
+
\end_layout
\begin_layout Subsection
@@ -11935,6 +12123,13 @@
\begin_layout Standard
Gales supports two types of Buoyancy forces.
+ The first one,
+\family typewriter
+BuoyancyForceTerm
+\family default
+, is more general, allowing you to specify buoyancy properties for each
+ material.
+
\end_layout
\begin_layout Subsection
Added: long/3D/Gale/trunk/documentation/images/Paraview_normal_stress.png
===================================================================
(Binary files differ)
Property changes on: long/3D/Gale/trunk/documentation/images/Paraview_normal_stress.png
___________________________________________________________________
Name: svn:mime-type
+ image/x-png
Modified: long/3D/Gale/trunk/documentation/images/Paraview_thermal.png
===================================================================
(Binary files differ)
Deleted: long/3D/Gale/trunk/documentation/images/Paraview_viscous_open_boundary.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/input/cookbook/thermal.xml
===================================================================
--- long/3D/Gale/trunk/input/cookbook/thermal.xml 2007-06-25 09:18:26 UTC (rev 7472)
+++ long/3D/Gale/trunk/input/cookbook/thermal.xml 2007-06-25 10:06:14 UTC (rev 7473)
@@ -370,6 +370,12 @@
<param name="density">1.0</param>
<param name="alpha">1.0</param>
<param name="diffusivity">1.0</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">1.0</param>
+ <param name="lambda">1.0</param>
+ </struct>
+ </list>
<list name="Rheology">
<param>backgroundViscosity</param>
<param>storeViscosity</param>
@@ -388,6 +394,12 @@
<param name="density">1.0</param>
<param name="alpha">10.0</param>
<param name="diffusivity">10.0</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">1000.0</param>
+ <param name="lambda">10.0</param>
+ </struct>
+ </list>
<list name="Rheology">
<param>sphereViscosity</param>
<param>storeViscosity</param>
@@ -468,6 +480,11 @@
<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>
@@ -570,7 +587,7 @@
<struct>
<param name="name">vx</param>
<param name="type">double</param>
- <param name="value">0.0001</param>
+ <param name="value">0.01</param>
</struct>
</list>
</struct>
Added: long/3D/Gale/trunk/input/cookbook/viscous_normal_stress.xml
===================================================================
--- long/3D/Gale/trunk/input/cookbook/viscous_normal_stress.xml 2007-06-25 09:18:26 UTC (rev 7472)
+++ long/3D/Gale/trunk/input/cookbook/viscous_normal_stress.xml 2007-06-25 10:06:14 UTC (rev 7473)
@@ -0,0 +1,498 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+ <struct name="components">
+ <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>
+ </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>
+ <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-constant</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>
+ <list name="BaseVariables">
+ <param>pressure</param>
+ </list>
+ </struct>
+ <struct name="PressureField">
+ <param name="Type">FeVariable</param>
+ <param name="FEMesh">mesh-constant</param>
+ <param name="DofLayout">pressureDofLayout</param>
+ <param name="LinkedDofInfo">pressureLinkedDofs</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="gaussMaterialSwarm">
+ <param name="Type">MaterialPointsSwarm</param>
+ <param name="CellLayout">cellLayout</param>
+ <param name="ParticleLayout">backgroundLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ </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="materialSwarm">
+ <param name="Type">MaterialPointsSwarm</param>
+ <param name="CellLayout">elementCellLayout</param>
+ <param name="ParticleLayout">materialSwarmParticleLayout</param>
+ <param name="FeMesh">mesh-linear</param>
+ <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="splittingRoutine">
+ <param name="Type">ReseedSplitting</param>
+ <param name="idealParticleCount">particlesPerCell</param>
+ <param name="originalParticleShape">bottomBoundaryShape</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">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="Preconditioner">preconditioner</param>
+ <param name="tolerance">1.0e-5</param>
+ <param name="maxIterations">5000</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="buoyancyForceTerm">
+ <param name="Type">BuoyancyForceTerm</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="gravity">1.0</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">StripRemesher</param>
+ <param name="mesh">mesh-linear</param>
+ <param name="meshType">regular</param>
+ <list name="dim">
+ <param>true</param>
+ <param>true</param>
+ <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>
+
+ <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">0.05</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">1.0</param>
+ <list name="Rheology">
+ <param>backgroundViscosity</param>
+ <param>storeViscosity</param>
+ <param>storeStress</param>
+ </list>
+ </struct>
+
+ <struct name="stressBC">
+ <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">double</param>
+ <param name="y_value">0.35</param>
+ </struct>
+
+ </struct>
+ <list name="plugins">
+ <param>Underworld_EulerDeform</param>
+ <param>Underworld_DumpSwarm</param>
+ <param>Underworld_VTKOutput</param>
+ </list>
+ <param name="maxTimeSteps">10</param>
+ <param name="outputEvery">1</param>
+ <param name="dumpEvery">1</param>
+ <param name="outputPath">./output.template</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">2.0f</param>
+ <param name="maxY">0.35f</param>
+ <param name="maxZ">0.3f</param>
+ <param name="elementResI">32</param>
+ <param name="elementResJ">16</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-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>
+ <param name="wrapTop">True</param>
+ <param name="staticBottom">True</param>
+ <list name="fields">
+ <struct>
+ <param name="field">VelocityField</param>
+ <param name="variable">velocity</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">0.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">1.0</param>
+ </struct>
+ <struct>
+ <param name="name">vy</param>
+ <param name="type">double</param>
+ <param name="value">0.0</param>
+ </struct>
+ </list>
+ </struct>
+ </list>
+ </struct>
+ <param name="checkpointEvery">1</param>
+</StGermainData>
Deleted: long/3D/Gale/trunk/input/cookbook/viscous_open_boundary.xml
===================================================================
--- long/3D/Gale/trunk/input/cookbook/viscous_open_boundary.xml 2007-06-25 09:18:26 UTC (rev 7472)
+++ long/3D/Gale/trunk/input/cookbook/viscous_open_boundary.xml 2007-06-25 10:06:14 UTC (rev 7473)
@@ -1,498 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
- <struct name="components">
- <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>
- </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>
- <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-constant</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>
- <list name="BaseVariables">
- <param>pressure</param>
- </list>
- </struct>
- <struct name="PressureField">
- <param name="Type">FeVariable</param>
- <param name="FEMesh">mesh-constant</param>
- <param name="DofLayout">pressureDofLayout</param>
- <param name="LinkedDofInfo">pressureLinkedDofs</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="gaussMaterialSwarm">
- <param name="Type">MaterialPointsSwarm</param>
- <param name="CellLayout">cellLayout</param>
- <param name="ParticleLayout">backgroundLayout</param>
- <param name="FeMesh">mesh-linear</param>
- </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="materialSwarm">
- <param name="Type">MaterialPointsSwarm</param>
- <param name="CellLayout">elementCellLayout</param>
- <param name="ParticleLayout">materialSwarmParticleLayout</param>
- <param name="FeMesh">mesh-linear</param>
- <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="splittingRoutine">
- <param name="Type">ReseedSplitting</param>
- <param name="idealParticleCount">particlesPerCell</param>
- <param name="originalParticleShape">bottomBoundaryShape</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">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="Preconditioner">preconditioner</param>
- <param name="tolerance">1.0e-5</param>
- <param name="maxIterations">5000</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="buoyancyForceTerm">
- <param name="Type">BuoyancyForceTerm</param>
- <param name="ForceVector">mom_force</param>
- <param name="Swarm">picIntegrationPoints</param>
- <param name="gravity">1.0</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">StripRemesher</param>
- <param name="mesh">mesh-linear</param>
- <param name="meshType">regular</param>
- <list name="dim">
- <param>true</param>
- <param>true</param>
- <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>
-
- <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">0.05</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">1.0</param>
- <list name="Rheology">
- <param>backgroundViscosity</param>
- <param>storeViscosity</param>
- <param>storeStress</param>
- </list>
- </struct>
-
- <struct name="stressBC">
- <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">double</param>
- <param name="y_value">0.35</param>
- </struct>
-
- </struct>
- <list name="plugins">
- <param>Underworld_EulerDeform</param>
- <param>Underworld_DumpSwarm</param>
- <param>Underworld_VTKOutput</param>
- </list>
- <param name="maxTimeSteps">10</param>
- <param name="outputEvery">1</param>
- <param name="dumpEvery">1</param>
- <param name="outputPath">./output.template</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">2.0f</param>
- <param name="maxY">0.35f</param>
- <param name="maxZ">0.3f</param>
- <param name="elementResI">32</param>
- <param name="elementResJ">16</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-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>
- <param name="wrapTop">True</param>
- <param name="staticBottom">True</param>
- <list name="fields">
- <struct>
- <param name="field">VelocityField</param>
- <param name="variable">velocity</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">0.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">1.0</param>
- </struct>
- <struct>
- <param name="name">vy</param>
- <param name="type">double</param>
- <param name="value">0.0</param>
- </struct>
- </list>
- </struct>
- </list>
- </struct>
- <param name="checkpointEvery">1</param>
-</StGermainData>
More information about the cig-commits
mailing list