[cig-commits] r16344 - in short/3D/PyLith/trunk/examples/3d/hex8: . powerlaw
willic3 at geodynamics.org
willic3 at geodynamics.org
Thu Feb 25 16:45:50 PST 2010
Author: willic3
Date: 2010-02-25 16:45:50 -0800 (Thu, 25 Feb 2010)
New Revision: 16344
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt
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/savageprescott.cfg
Log:
Cleaned up savageprescott.cfg. Power-law is now the default.
Altered temperatures to give more viscous strain.
Changed points definitions to have a clean break at 2 km depth.
Changed temperature interpolation to linear.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg 2010-02-25 23:07:29 UTC (rev 16343)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg 2010-02-26 00:45:50 UTC (rev 16344)
@@ -17,6 +17,7 @@
[powerlaw_gendb.db_temperature]
label = Temperature
+query_type = linear
iohandler.filename = temperature.spatialdb
[powerlaw_gendb.geometry]
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt 2010-02-25 23:07:29 UTC (rev 16343)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt 2010-02-26 00:45:50 UTC (rev 16344)
@@ -3,7 +3,8 @@
0.0 0.0 -500.0
0.0 0.0 -1000.0
0.0 0.0 -1500.0
-0.0 0.0 -2000.0
+0.0 0.0 -1900.0
+0.0 0.0 -2100.0
0.0 0.0 -2500.0
0.0 0.0 -3000.0
0.0 0.0 -3500.0
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-25 23:07:29 UTC (rev 16343)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-26 00:45:50 UTC (rev 16344)
@@ -2,7 +2,7 @@
//
// This spatial database defines a simple linear 1D temperature profile,
// assuming a temperature of 100 degrees C (373 degrees K) at the surface,
-// increasing to 650 degrees C (923 degrees K) at 4 km depth.
+// increasing to 1000 degrees C (1273 degrees K) at 4 km depth.
//
#SPATIAL.ascii 1
SimpleDB {
@@ -23,4 +23,4 @@
// (3) z coordinate (m)
// (4) temperature (K)
0.0 0.0 0.0 373.0
-0.0 0.0 -4.0 923.0
+0.0 0.0 -4.0 1273.0
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-25 23:07:29 UTC (rev 16343)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-26 00:45:50 UTC (rev 16344)
@@ -3,7 +3,7 @@
num-values = 3
value-names = reference-stress reference-strain-rate power-law-exponent
value-units = Pa 1/s none
- num-locs = 9
+ num-locs = 10
data-dim = 1
space-dim = 3
cs-data = cartesian {
@@ -12,11 +12,12 @@
}
}
0.000000e+00 0.000000e+00 0.000000e+00 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -5.000000e+02 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.000000e+03 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.500000e+03 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.000000e+03 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.500000e+03 1.935443e+11 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.000000e+03 1.935443e+11 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.500000e+03 1.935443e+11 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -4.000000e+03 1.935443e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -5.000000e+02 1.977176e+13 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.000000e+03 2.801905e+11 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.500000e+03 1.528518e+10 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.900000e+03 2.687566e+09 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.100000e+03 1.201397e+12 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -2.500000e+03 1.483187e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.000000e+03 1.798919e+10 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.500000e+03 3.284422e+09 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -4.000000e+03 8.099233e+08 1.000000e-06 3.500000e+00
Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-25 23:07:29 UTC (rev 16343)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-26 00:45:50 UTC (rev 16344)
@@ -23,17 +23,11 @@
# 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
+max_dt = 10.0*year ; time step
stability_factor = 1.0 ; use time step equal to stable value from materials
[pylithapp.timedependent.implicit]
@@ -48,9 +42,9 @@
# 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.
+# Users can switch between linear Maxwell rheology and power-law
+# rheology by commenting and uncommenting the lines in the appropriate
+# section.
# *** BEGIN MAXWELL SECTION ***
# [pylithapp.timedependent.materials.viscoelastic]
# db_properties.iohandler.filename = mat_maxwell.spatialdb
@@ -58,10 +52,6 @@
# output.cell_data_fields = [total_strain,stress,viscous_strain]
# *** 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
@@ -171,7 +161,7 @@
output_freq = time_step
time_step = 5.0*year
writer.filename = savageprescott.vtk
-writer.time_format = %04.0f
+writer.time_format = %06.1f
writer.time_constant = 1.0*year
# Give basename for VTK domain output of solution over ground surface.
@@ -179,14 +169,14 @@
label = 17 ; nodeset for subdomain
skip = 0
writer.filename = savageprescott-groundsurf.vtk
-writer.time_format = %04.0f
+writer.time_format = %06.1f
writer.time_constant = 1.0*year
# Give basename for vtk fault rupture output.
[pylithapp.timedependent.interfaces.fault.output]
skip = 0
writer.filename = savageprescott-fault.vtk
-writer.time_format = %04.0f
+writer.time_format = %06.1f
writer.time_constant = 1.0*year
# Give basename for VTK output of state variables.
@@ -195,7 +185,7 @@
time_step = 5.0*year
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = savageprescott-elastic.vtk
-writer.time_format = %04.0f
+writer.time_format = %06.1f
writer.time_constant = 1.0*year
[pylithapp.timedependent.materials.viscoelastic.output]
@@ -203,5 +193,5 @@
time_step = 5.0*year
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = savageprescott-viscoelastic.vtk
-writer.time_format = %04.0f
+writer.time_format = %06.1f
writer.time_constant = 1.0*year
More information about the CIG-COMMITS
mailing list