[cig-commits] r16316 - in short/3D/PyLith/trunk/examples/3d/hex8: . figures powerlaw
willic3 at geodynamics.org
willic3 at geodynamics.org
Mon Feb 22 20:01:00 PST 2010
Author: willic3
Date: 2010-02-22 20:00:59 -0800 (Mon, 22 Feb 2010)
New Revision: 16316
Added:
short/3D/PyLith/trunk/examples/3d/hex8/figures/
short/3D/PyLith/trunk/examples/3d/hex8/figures/dislocation.jpg
short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-gravity.jpg
short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-gravity_istress.jpg
short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-mesh.jpg
short/3D/PyLith/trunk/examples/3d/hex8/figures/savageprescott-t100.jpg
short/3D/PyLith/trunk/examples/3d/hex8/figures/shear.jpg
Removed:
short/3D/PyLith/trunk/examples/3d/hex8/dislocation.jpg
short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity.jpg
short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity_istress.jpg
short/3D/PyLith/trunk/examples/3d/hex8/hex8-mesh.jpg
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-t100.jpg
short/3D/PyLith/trunk/examples/3d/hex8/shear.jpg
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
Log:
Made some changes to make it easier to use power-law for Savage & Prescoot
example, and fixed some problems.
There will be some more changes after I finish debugging power-law.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -28,6 +28,16 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+
+# We specifically define the database info for the viscoelastic material
+# for each model, since the type of spatialdb can be different.
+
+[pylithapp.timedependent.materials.viscoelastic]
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+
+# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/dislocation.jpg
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/dislocation.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/dislocation.jpg)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-gravity.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity.jpg)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-gravity_istress.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity_istress.jpg)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/hex8-mesh.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/hex8-mesh.jpg)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/savageprescott-t100.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-t100.jpg)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/figures/shear.jpg (from rev 16301, short/3D/PyLith/trunk/examples/3d/hex8/shear.jpg)
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -30,6 +30,16 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+
+# We specifically define the database info for the viscoelastic material
+# for each model, since the type of spatialdb can be different.
+
+[pylithapp.timedependent.materials.viscoelastic]
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+
+# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
# Use the default FixedDOFDB, which has a uniform spatial distribution.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -30,6 +30,16 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+
+# We specifically define the database info for the viscoelastic material
+# for each model, since the type of spatialdb can be different.
+
+[pylithapp.timedependent.materials.viscoelastic]
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+
+# ----------------------------------------------------------------------
# initial stresses
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials.elastic]
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity.jpg
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/hex8-gravity_istress.jpg
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/hex8-mesh.jpg
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-23 04:00:59 UTC (rev 16316)
@@ -1,8 +1,8 @@
// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
//
// This spatial database defines a simple linear 1D temperature profile,
-// assuming a temperature of 200 degrees C (473 degrees K) at the surface,
-// increasing to 500 degrees C (773 degrees K) at 4 km depth.
+// assuming a temperature of 0 degrees C (273 degrees K) at the surface,
+// increasing to 400 degrees C (773 degrees K) at 4 km depth.
//
#SPATIAL.ascii 1
SimpleDB {
@@ -22,5 +22,5 @@
// (2) y coordinate (m)
// (3) z coordinate (m)
// (4) temperature (K)
-0.0 0.0 0.0 473.0
-0.0 0.0 -4.0 773.0
+0.0 0.0 0.0 273.0
+0.0 0.0 -4.0 673.0
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-23 04:00:59 UTC (rev 16316)
@@ -1,7 +1,7 @@
#SPATIAL.ascii 1
SimpleDB {
num-values = 3
- value-names = reference-stress reference-strain-rate powerlaw-exponent
+ value-names = reference-stress reference-strain-rate power-law-exponent
value-units = Pa 1/s none
num-locs = 9
data-dim = 1
@@ -11,12 +11,12 @@
space-dim = 3
}
}
- 0.000000e+00 0.000000e+00 0.000000e+00 3.595219e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -5.000000e+02 3.595219e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.000000e+03 3.595219e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.500000e+03 3.595219e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.000000e+03 3.595219e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.500000e+03 9.233654e+12 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.000000e+03 9.233654e+12 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.500000e+03 9.233654e+12 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -4.000000e+03 9.233654e+12 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 0.000000e+00 8.804741e+20 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -5.000000e+02 8.804741e+20 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.000000e+03 8.804741e+20 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.500000e+03 8.804741e+20 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.000000e+03 8.804741e+20 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.500000e+03 3.163259e+14 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.000000e+03 3.163259e+14 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.500000e+03 3.163259e+14 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -4.000000e+03 3.163259e+14 1.000000e-06 3.500000e+00
Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -15,7 +15,8 @@
faultcohesivekin = 1
fiatlagrange = 1
pylithapp = 1
-##problem = 1
+materials = 1
+powerlaw3d = 1
##implicit = 1
# ----------------------------------------------------------------------
@@ -58,10 +59,11 @@
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
+# We do not define the database info for the viscoelastic material, since
+# a different database type may be used for some examples.
[pylithapp.timedependent.materials.viscoelastic]
label = Viscoelastic material
id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-t100.jpg
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -23,12 +23,18 @@
# time stepping.
implicit.time_step = pylith.problems.TimeStepAdapt
+# Test option to turn off normalizing.
+normalizer.length_scale = 1.0*m
+# normalizer.pressure_scale = 1.0*Pa
+# normalizer.time_scale = 1.0*second
+# normalizer.density_scale = 1.0*kg/m**3
+
[pylithapp.timedependent.implicit.time_step]
# Define the total time for the simulation and the default time step size.
# Overwrite values in pylithapp.cfg
total_time = 700.0*year ; total time of simulation
max_dt = 50.0*year ; time step
-stability_factor = 2.0 ; use time step equal to stable value from materials
+stability_factor = 1.0 ; use time step equal to stable value from materials
[pylithapp.timedependent.implicit]
# Set the output to an array of 2 output managers.
@@ -38,23 +44,45 @@
# Set subdomain component to OutputSolnSubset (subset of domain).
output.subdomain = pylith.meshio.OutputSolnSubset
-[pylithapp.timedependent.materials.viscoelastic]
-output.cell_info_fields = [density,mu,lambda,maxwell_time] ; Maxwell model
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+# We specifically define the viscoelastic material for each model, since
+# the type of spatialdb can be different.
+# Comment out the lines below to avoid using a linear Maxwell rheology.
+# *** BEGIN MAXWELL SECTION ***
+# [pylithapp.timedependent.materials.viscoelastic]
+# db_properties.iohandler.filename = mat_maxwell.spatialdb
+# output.cell_info_fields = [density,mu,lambda,maxwell_time]
+# *** END MAXWELL SECTION ***
+
# ----------------------------------------------------------------------
# Uncomment the lines below to switch to the power-law rheology and
# nonlinear solver.
# ----------------------------------------------------------------------
+# *** BEGIN POWER-LAW SECTION ***
+[pylithapp.timedependent]
+implicit.solver = pylith.problems.SolverNonlinear
+implicit.time_step.max_dt = 10.0*year
+materials.viscoelastic = pylith.materials.PowerLaw3D
-# [pylithapp.timedependent]
-# implicit.solver = pylith.problems.SolverNonlinear
-# implicit.time_step.max_dt = 10.0*year
-# materials.viscoelastic = pylith.materials.PowerLaw3D
+[pylithapp.timedependent.materials.viscoelastic]
+db_properties = spatialdata.spatialdb.CompositeDB
+db_properties.db_A = spatialdata.spatialdb.SimpleDB
+db_properties.db_B = spatialdata.spatialdb.SimpleDB
+output.cell_info_fields = [density,mu,lambda,reference_strain_rate,reference_stress,power_law_exponent]
-# [pylithapp.timedependent.materials.viscoelastic]
-# db_properties.iohandler.filename = powerlaw_properties.spatialdb
-# output.cell_info_fields = [density,mu,lambda,reference_strain_rate,reference_stress,power_law_exponent]
+[pylithapp.timedependent.materials.viscoelastic.db_properties]
+values_A = [density,vs,vp]
+db_A.label = Elastic properties
+db_A.iohandler.filename = mat_elastic.spatialdb
+values_B = [reference-stress,reference-strain-rate,power-law-exponent]
+db_B.label = Power-law properties
+db_B.iohandler.filename = powerlaw_properties.spatialdb
+# *** END POWER-LAW SECTION ***
+
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/shear.jpg
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2010-02-23 02:41:02 UTC (rev 16315)
+++ short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2010-02-23 04:00:59 UTC (rev 16316)
@@ -26,6 +26,16 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+
+# We specifically define the database info for the viscoelastic material
+# for each model, since the type of spatialdb can be different.
+
+[pylithapp.timedependent.materials.viscoelastic]
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+
+# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
# We change spatial database for initial value from FixedDOFDB (which has
More information about the CIG-COMMITS
mailing list