[cig-commits] commit: Add deformable lower boundary to cookbook
Mercurial
hg at geodynamics.org
Thu Jul 15 12:08:04 PDT 2010
changeset: 294:c86df357143d
user: Walter Landry <wlandry at caltech.edu>
date: Tue Jul 13 19:38:17 2010 -0700
files: doc/gale.lyx doc/images/Paraview_bottom.png input/cookbook/viscous_bottom.xml
description:
Add deformable lower boundary to cookbook
diff -r 3750621515f0 -r c86df357143d doc/gale.lyx
--- a/doc/gale.lyx Thu Jun 17 11:00:52 2010 -0700
+++ b/doc/gale.lyx Tue Jul 13 19:38:17 2010 -0700
@@ -1,4 +1,4 @@
-#LyX 1.6.5 created this file. For more info see http://www.lyx.org/
+#LyX 1.6.6 created this file. For more info see http://www.lyx.org/
\lyxformat 345
\begin_document
\begin_header
@@ -775,16 +775,11 @@ reference "eq:body force"
\begin_layout Standard
\begin_inset Formula \begin{equation}
-\frac{DT}{Dt}\equiv\frac{\partial T}{\partial t}+v\cdot\nabla T=\kappa\nabla^{2}T+Q,\label{eq:energy}\end{equation}
+\frac{\partial T}{\partial t}+v\cdot\nabla T=\kappa\nabla^{2}T+Q,\label{eq:energy}\end{equation}
\end_inset
where
-\begin_inset Formula $\frac{D}{Dt}$
-\end_inset
-
- is the material time derivative (taken at a point moving with respect to
- the fluid),
\begin_inset Formula $T$
\end_inset
@@ -1821,23 +1816,11 @@ Gale works on a variety of computational
\end_layout
\begin_layout Itemize
-Mac OS X 10.4.6 (G4, G5, and Intel)
+Mac OS X 10.4.6 (Intel)
\end_layout
\begin_layout Itemize
-Windows 2000 and XP SP2
-\end_layout
-
-\begin_layout Itemize
-RedHat Fedora Core 5 (x86)
-\end_layout
-
-\begin_layout Itemize
-OpenSuse 10.0 (x86)
-\end_layout
-
-\begin_layout Itemize
-Gentoo (x86)
+Windows Vista
\end_layout
\begin_layout Itemize
@@ -1845,8 +1828,8 @@ Debian stable (x86 and AMD64), testing (
\end_layout
\begin_layout Standard
-Gale has also been tested on clusters running Redhat 7.2 (x86) and RedHat
- Enterprise Linux 3 (EM64T).
+Gale has also been tested on clusters running RedHat Enterprise Linux 3
+ (EM64T).
\end_layout
\begin_layout Subsection
@@ -1874,7 +1857,8 @@ You must also have python 2.2.1 or great
You must also have python 2.2.1 or greater installed.
If you do not already have MPI, then in many cases PETSc can install a
version for you.
- Installing PETSc also requires a Blas/Lapack implementation.
+ Installing PETSc also requires a Blas/Lapack implementation, which, again,
+ PETSc can install for you.
\end_layout
\begin_layout Subsection
@@ -2345,7 +2329,6 @@ In parallel, the analogous approach woul
\end_inset
when configuring.
- The PETSc installation on TACC Lonestar has been installed with Mumps.
\end_layout
\begin_layout Standard
@@ -2393,7 +2376,7 @@ You can also change the default values o
You can also change the default values of the input file without modifying
that file by appending arguments.
For example, to change only the number of time steps from the default value
- 10 to value 20, use the following command
+ of 10 to 20, use the following command
\end_layout
\begin_layout LyX-Code
@@ -2416,63 +2399,6 @@ 64
--elementResJ=64 --elementResK=64 --particlesPerCell=60 --dumpEvery=10
\end_layout
-\begin_layout Subsection
-Parallel Issues
-\end_layout
-
-\begin_layout Standard
-When writing
-\family typewriter
-.dat
-\family default
- files in parallel (see Sections
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Checkpoint Files"
-
-\end_inset
-
- and
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:.dat-Binary-Files"
-
-\end_inset
-
-), all of the processors have to cooperate.
- This is because the
-\family typewriter
-.dat
-\family default
- files are single files that have information from every processor.
- Gale uses file locking to ensure this cooperation.
- However, on some parallel file systems, this can be problematic.
- The only known solution at this point is to turn off checkpointing by setting
- the
-\family typewriter
-checkpointEvery
-\family default
- parameter to 0.
-
-\end_layout
-
-\begin_layout Standard
-The VTK files do not suffer from this problem.
- However, the VTK files may contain a large amount of data.
- Writing all of that information to disk may slow down your computation
- significantly.
- The easiest way to tell if this is the case is to remove the line
-\end_layout
-
-\begin_layout LyX-Code
-<param>Underworld_VTKOutput</param>
-\end_layout
-
-\begin_layout Standard
-If this makes a big difference, then you can output less frequently or output
- fewer particles.
-\end_layout
-
\begin_layout Section
Output and Visualization
\end_layout
@@ -2598,7 +2524,7 @@ name "sub:Visualization-with-ParaView"
\begin_layout Standard
\color none
-These instructions are for Paraview version 3.6.0.
+These instructions are for Paraview version 3.6.1.
To visualize step 10 of
\family typewriter
input/cookbook/yielding.xml
@@ -3483,26 +3409,11 @@ particlesPerCell
\end_layout
\begin_layout Standard
-In addition, you need to vary the scaling factor for time steps (
+In addition, you may need to vary the scaling factor for time steps (
\family typewriter
dtFactor
\family default
).
- For problems with inflow boundaries (see Sections
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Viscous-In-Outflow"
-
-\end_inset
-
- and
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sub:Flux-Boundary-Conditions"
-
-\end_inset
-
-), you must vary the size of the boundary box.
Depending on the boundary conditions, you may need to vary the size of
the box (
\family typewriter
@@ -4781,7 +4692,7 @@ reference "sec:Viscous-Material-in"
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
+ Then, when material piles up, gravity forces 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.
@@ -5172,6 +5083,87 @@ Strain rate invariant and velocities of
\end_layout
\begin_layout Section
+Viscous Material with Deformable Bottom Boundary
+\end_layout
+
+\begin_layout Standard
+The previous example can be modified so that, instead of having material
+ flow through the bottom boundary, the boundary itself deforms.
+ You can do this by changing the one line
+\end_layout
+
+\begin_layout LyX-Code
+<param name="staticBottom">True</param>
+\end_layout
+
+\begin_layout Standard
+to
+\end_layout
+
+\begin_layout LyX-Code
+<param name="wrapBottom">True</param>
+\end_layout
+
+\begin_layout Standard
+A worked example is in
+\family typewriter
+input/cookbook/viscous_bottom.xml
+\family default
+.
+ Figure
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "fig:Strain-deformable-bottom"
+
+\end_inset
+
+ shows the strain rate invariant and velocity.
+\end_layout
+
+\begin_layout Standard
+\begin_inset Float figure
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\align center
+\begin_inset Graphics
+ filename images/Paraview_bottom.png
+ lyxscale 75
+ width 100col%
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Plain Layout
+\begin_inset Caption
+
+\begin_layout Plain Layout
+\begin_inset CommandInset label
+LatexCommand label
+name "fig:Strain-deformable-bottom"
+
+\end_inset
+
+Strain rate invariant and velocities of viscous material with a deformable
+ bottom boundary
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Section
\begin_inset CommandInset label
LatexCommand label
name "sec:Viscous-Material-with"
@@ -7150,7 +7142,7 @@ reference "fig:Temperature-and-velocity"
\end_inset
- is mostly because the viscosity is higher everywhere.
+ are mostly because the viscosity is higher everywhere.
\end_layout
\begin_layout LyX-Code
@@ -7586,9 +7578,9 @@ Software Components of Gale
\begin_layout Standard
Gale makes use of several physics libraries, including StGermain, StgFEM,
PiCellerator, and UnderWorld.
- These are proven, capable, open-source finite element method libraries
- written by the Victorian Partnership for Advanced Computing (VPAC) and
- Louis Moresi's group at Monash University (see Fig
+ These are open-source finite element method libraries written by the Victorian
+ Partnership for Advanced Computing (VPAC) and Louis Moresi's group at Monash
+ University (see Fig
\begin_inset CommandInset ref
LatexCommand ref
reference "fig:Mapping-between-MicroFEM"
diff -r 3750621515f0 -r c86df357143d doc/images/Paraview_bottom.png
Binary file doc/images/Paraview_bottom.png has changed
diff -r 3750621515f0 -r c86df357143d input/cookbook/viscous_bottom.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/input/cookbook/viscous_bottom.xml Tue Jul 13 19:38:17 2010 -0700
@@ -0,0 +1,516 @@
+<?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">timeIntegratorOrder</param>
+ <param name="simultaneous">f</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="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="Swarm">picIntegrationPoints</param>
+ <param name="gravity">gravity</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">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="FieldVariablesToCheckpoint">
+ <param>VelocityField</param>
+ <param>PressureField</param>
+ </list>
+ <param name="timeIntegratorOrder">1</param>
+ <param name="maxTimeSteps">10</param>
+ <param name="outputEvery">1</param>
+ <param name="dumpEvery">1</param>
+ <param name="outputPath">./output</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-linear</param>
+ <param name="remesher">velocityRemesher</param>
+ <param name="velocityField">VelocityField</param>
+ <param name="wrapTop">True</param>
+ <param name="wrapBottom">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">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>
+ <param name="gravity">1.0</param>
+</StGermainData>
More information about the CIG-COMMITS
mailing list