[cig-commits] r7744 -
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Jul 25 11:08:38 PDT 2007
Author: willic3
Date: 2007-07-25 11:08:37 -0700 (Wed, 25 Jul 2007)
New Revision: 7744
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_maxwell.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/timedep.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_lagrit.cfg
Log:
Added timedep.cfg and mat_maxwell.spatialdb to allow solution of
time-dependent problem (10 years) using Maxwell viscoelastic model, and
made necessary changes to other files to make things work correctly.
Also altered material databases to make elastic properties consistent
with benchmark definition (shear modulus = 30 GPa, Poisson's ratio =
0.25).
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README 2007-07-25 18:08:37 UTC (rev 7744)
@@ -6,6 +6,8 @@
pylithapp.cfg - General parameters for problem
+timedep.cfg - Parameters for time-dependent problem
+
strikeslip_cubit.cfg - Parameters common to CUBIT generated meshes
strikeslip_hex8_1000m.cfg - Parameters specific to hex8 1000m resolution sim
strikeslip_hex8_0500m.cfg - Parameters specific to hex8 500m resolution sim
@@ -44,3 +46,42 @@
tet4 @ 250m resolution
pylith strikeslip_lagrit.cfg strikeslip_tet4_0250m.cfg
+
+
+======================================================================
+All of the simulations above are for a single elastic time step.
+To run time-dependent simulations, you can simply append 'timedep.cfg'
+to the end of each example listed above.
+
+TO RUN TIME-DEPENDENT SIMULATIONS
+======================================================================
+
+hex8 @ 1000m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_1000m.cfg timedep.cfg
+
+hex8 @ 500m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_0500m.cfg timedep.cfg
+
+hex8 @ 250m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_0250m.cfg timedep.cfg
+
+tet4 @ 1000m resolution
+
+ pylith strikeslip_lagrit.cfg strikeslip_tet4_1000m.cfg timedep.cfg
+
+tet4 @ 500m resolution
+
+ pylith strikeslip_lagrit.cfg strikeslip_tet4_0500m.cfg timedep.cfg
+
+tet4 @ 250m resolution
+
+ pylith strikeslip_lagrit.cfg strikeslip_tet4_0250m.cfg timedep.cfg
+
+======================================================================
+NOTE: Any of the above problems can be run on more than one processor
+by appending -nodes=n, where n is the number of desired processes, to
+the end of any of the above example commands.
+======================================================================
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb 2007-07-25 18:08:37 UTC (rev 7744)
@@ -11,4 +11,4 @@
space-dim = 3
}
}
-0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_maxwell.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_maxwell.spatialdb 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_maxwell.spatialdb 2007-07-25 18:08:37 UTC (rev 7744)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = density vs vp viscosity
+ value-units = kg/m^3 m/s m/s Pa*s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 1.0e18
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg 2007-07-25 18:08:37 UTC (rev 7744)
@@ -42,8 +42,8 @@
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials]
# TEMPORARY
-viscoelastic_neg = pylith.materials.ElasticIsotropic3D
-viscoelastic_pos = pylith.materials.ElasticIsotropic3D
+viscoelastic_neg = pylith.materials.MaxwellIsotropic3D
+viscoelastic_pos = pylith.materials.MaxwellIsotropic3D
# ----------------------------------------------------------------------
# boundary conditions
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg 2007-07-25 18:08:37 UTC (rev 7744)
@@ -27,9 +27,9 @@
quadrature.cell.dimension = 3
[pylithapp.timedependent.materials.viscoelastic_pos]
-label = Elastic material +x
+label = Viscoelastic material +x
id = 3
-db.iohandler.filename = mat_elastic.spatialdb
+db.iohandler.filename = mat_maxwell.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
@@ -37,7 +37,7 @@
[pylithapp.timedependent.materials.viscoelastic_neg]
label = Viscoelastic material -x
id = 4
-db.iohandler.filename = mat_elastic.spatialdb
+db.iohandler.filename = mat_maxwell.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_lagrit.cfg 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_lagrit.cfg 2007-07-25 18:08:37 UTC (rev 7744)
@@ -27,9 +27,9 @@
quadrature.cell.shape = tetrahedron
[pylithapp.timedependent.materials.viscoelastic_pos]
-label = Elastic material +x
+label = Viscoelastic material +x
id = 3
-db.iohandler.filename = mat_elastic.spatialdb
+db.iohandler.filename = mat_maxwell.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
@@ -37,7 +37,7 @@
[pylithapp.timedependent.materials.viscoelastic_neg]
label = Viscoelastic material -x
id = 4
-db.iohandler.filename = mat_elastic.spatialdb
+db.iohandler.filename = mat_maxwell.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/timedep.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/timedep.cfg 2007-07-24 21:34:10 UTC (rev 7743)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/timedep.cfg 2007-07-25 18:08:37 UTC (rev 7744)
@@ -0,0 +1,16 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 10.0*year
+default_dt = 0.1*year
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.solutioniovtk]
+output_freq = time_step
+time_step = 1.0*year
More information about the cig-commits
mailing list