[cig-commits] r16249 - in short/3D/PyLith/trunk/examples/3d/hex8: . powerlaw
willic3 at geodynamics.org
willic3 at geodynamics.org
Tue Feb 9 16:21:06 PST 2010
Author: willic3
Date: 2010-02-09 16:21:05 -0800 (Tue, 09 Feb 2010)
New Revision: 16249
Added:
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_params.spatialdb
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
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
Log:
Put in directory to compute power-law properties from parameters and
temperature.
Still need to finish fixing .cfg file to get elastic and power-law
properties from different spatial databases.
Added: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_gendb.cfg 2010-02-10 00:21:05 UTC (rev 16249)
@@ -0,0 +1,28 @@
+# -*- Python -*-
+[powerlaw_gendb]
+
+reference_value = strain_rate
+
+[powerlaw_gendb.db_exponent]
+label = Power-law exponent
+iohandler.filename = powerlaw_params.spatialdb
+
+[powerlaw_gendb.db_activation_energy]
+label = Activation energy
+iohandler.filename = powerlaw_params.spatialdb
+
+[powerlaw_gendb.db_powerlaw_coefficient]
+label = Experimentally derived power-law coefficient
+iohandler.filename = powerlaw_params.spatialdb
+
+[powerlaw_gendb.db_temperature]
+label = Temperature
+iohandler.filename = temperature.spatialdb
+
+[powerlaw_gendb.geometry]
+data_dim = 1
+reader = spatialdata.utils.PointsStream
+reader.filename = powerlaw_points.txt
+
+[powerlaw_gendb.iohandler]
+filename = powerlaw_properties.spatialdb
Added: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_params.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_params.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_params.spatialdb 2010-02-10 00:21:05 UTC (rev 16249)
@@ -0,0 +1,38 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database defines the power-law parameters for two materials:
+// Wet granite and dry olivine. The properties are from Strehlau and
+// Meissner (1987). Wet granite is specified from zero to 2 km depth,
+// while dry olivine is specified from 2 to 4 km depth.
+// Note that the flow constant is expressed as log10(flow-constant), where
+// the units are GPa^(-n)/s. The activation energy (Q) is expressed in
+// kJ/mol. To accommodate the strange units of the flow constant, the
+// given units of are 'None', and we specify a new parameter
+// (flow-constant-multiplier), with a value of 9, to indicate that the
+// underlying units are GPa. A value of 6 would indicate MPa, etc.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = log-flow-constant activation-energy power-law-exponent flow-constant-scale // names of the material property values
+ value-units = None kJ/mol None None // units
+ num-locs = 4 // number of locations
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1000.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) flow-constant (None)
+// (5) activation-energy (kJ/mol)
+// (6) power-law-exponent (None)
+// (7) activation-energy-multiplier (None)
+0.0 0.0 0.0 2.0 137.0 1.5 9.0
+0.0 0.0 -1.9 2.0 137.0 1.5 9.0
+0.0 0.0 -2.1 15.5 535.0 3.5 9.0
+0.0 0.0 -4.0 15.5 535.0 3.5 9.0
Added: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/powerlaw_points.txt 2010-02-10 00:21:05 UTC (rev 16249)
@@ -0,0 +1,10 @@
+# Locations at which we want to output power-law properties.
+0.0 0.0 0.0
+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 -2500.0
+0.0 0.0 -3000.0
+0.0 0.0 -3500.0
+0.0 0.0 -4000.0
Added: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-10 00:21:05 UTC (rev 16249)
@@ -0,0 +1,26 @@
+// -*- 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.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1 // number of material property values
+ value-names = temperature // names of the material property values
+ value-units = K // units
+ num-locs = 2 // number of locations
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1000.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (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
Added: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-10 00:21:05 UTC (rev 16249)
@@ -0,0 +1,22 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = reference-stress reference-strain-rate powerlaw-exponent
+ value-units = Pa 1/s none
+ num-locs = 9
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1
+ 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
Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-10 00:10:20 UTC (rev 16248)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-10 00:21:05 UTC (rev 16249)
@@ -46,14 +46,14 @@
# nonlinear solver.
# ----------------------------------------------------------------------
-#[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.iohandler.filename = mat_powerlaw.spatialdb
-#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]
# ----------------------------------------------------------------------
# boundary conditions
More information about the CIG-COMMITS
mailing list