[cig-commits] [commit] williams/add-examples-grav2d: Small improvements to 2-D gravity example. (c5a9dbe)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Jun 20 12:19:51 PDT 2014
Repository : https://github.com/geodynamics/pylith
On branch : williams/add-examples-grav2d
Link : https://github.com/geodynamics/pylith/compare/99730884c86e4e5b96c472ce5a4e5094682941ee...c5a9dbe5776e87a83efa0040169d00fe3571af03
>---------------------------------------------------------------
commit c5a9dbe5776e87a83efa0040169d00fe3571af03
Author: Brad Aagaard <baagaard at usgs.gov>
Date: Fri Jun 20 12:19:44 2014 -0700
Small improvements to 2-D gravity example.
Fixed location of inner block. Adjusted solver tolerances (KSP smaller
than SNES).
>---------------------------------------------------------------
c5a9dbe5776e87a83efa0040169d00fe3571af03
examples/2d/gravity/README | 11 ++++++-----
examples/2d/gravity/geometry.jou | 7 ++++---
examples/2d/gravity/{mesh_tri3.jou => mesh.jou} | 2 +-
examples/2d/gravity/mesh_tri3.exo | Bin 866892 -> 0 bytes
examples/2d/gravity/pylithapp.cfg | 12 +++++++-----
5 files changed, 18 insertions(+), 14 deletions(-)
diff --git a/examples/2d/gravity/README b/examples/2d/gravity/README
index 333bdf9..5f09daa 100644
--- a/examples/2d/gravity/README
+++ b/examples/2d/gravity/README
@@ -9,7 +9,7 @@ different simulations, which should generally be performed in order:
1. pylith grav_stress.cfg
This will run a single elastic time step to compute the stresses due to
'turning on' gravity. These stresses will be used as initial stresses
- in one of the other simulations. To conver the output PyLith stresses
+ in one of the other simulations. To convert the stresses PyLith outputs
to a spatial database of initial stresses, run the stress2spatialdb.py
script:
./stress2spatialdb.py
@@ -36,7 +36,8 @@ different simulations, which should generally be performed in order:
This simulation uses a finite strain formulation, and should thus
produce less vertical deformation that the infinitesimal strain version.
-The results of all the different simulations may be viewed using ParaView,
-and a ParaView state file is contained in the output directory. By loading
-this state file into ParaView and changing the filename appropriately, the
-various simulation results may be viewed.
+The results of all the different simulations may be viewed using
+ParaView. We provide a ParaView state file (output/deform_view.pvsm),
+which you can load in ParaView to view the solution over the
+domain. The visualization is intended for steps 3 and 4 in which there
+is deformation.
diff --git a/examples/2d/gravity/geometry.jou b/examples/2d/gravity/geometry.jou
index 993ab75..efc9038 100644
--- a/examples/2d/gravity/geometry.jou
+++ b/examples/2d/gravity/geometry.jou
@@ -38,10 +38,10 @@ undo off
# Create outer box (50 km x 40 km x 10 km)..
# ----------------------------------------------------------------------
#{outerBlockX=40.0*km}
-#{outerBlockY=50.0*km}
+#{outerBlockY=30.0*km}
#{outerBlockZ=10.0*km}
brick x {outerBlockX} y {outerBlockY} z {outerBlockZ}
-volume 1 move x 0.0 y {-25.0*km} z 0.0
+volume 1 move x 0.0 y {-0.5*outerBlockY} z 0.0
# ----------------------------------------------------------------------
# Create inner box (2.5 km x 5 km x 10 km)..
@@ -50,7 +50,7 @@ volume 1 move x 0.0 y {-25.0*km} z 0.0
#{innerBlockY=5.0*km}
#{innerBlockZ=10.0*km}
brick x {innerBlockX} y {innerBlockY} z {innerBlockZ}
-volume 2 move x {-40.0*km} y {-2.5*km} z 0.0
+volume 2 move x 0.0 y {-0.5*innerBlockY} z 0.0
# ----------------------------------------------------------------------
# Webcut outer block with inner block and then form a midsurface
@@ -65,3 +65,4 @@ create midsurface volume 4 surface 20 22
delete volume 3 4
imprint all
merge all
+
diff --git a/examples/2d/gravity/mesh_tri3.jou b/examples/2d/gravity/mesh.jou
similarity index 98%
rename from examples/2d/gravity/mesh_tri3.jou
rename to examples/2d/gravity/mesh.jou
index 5451801..8a0b49e 100644
--- a/examples/2d/gravity/mesh_tri3.jou
+++ b/examples/2d/gravity/mesh.jou
@@ -72,7 +72,7 @@ nodeset 14 name "face_ypos"
# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
-export mesh "mesh_tri3.exo" dimension 2 overwrite
+export mesh "mesh.exo" dimension 2 overwrite
# End of file
diff --git a/examples/2d/gravity/mesh_tri3.exo b/examples/2d/gravity/mesh_tri3.exo
deleted file mode 100644
index 0899088..0000000
Binary files a/examples/2d/gravity/mesh_tri3.exo and /dev/null differ
diff --git a/examples/2d/gravity/pylithapp.cfg b/examples/2d/gravity/pylithapp.cfg
index 892e636..0d9b2b4 100644
--- a/examples/2d/gravity/pylithapp.cfg
+++ b/examples/2d/gravity/pylithapp.cfg
@@ -28,7 +28,7 @@ materials = 1
reader = pylith.meshio.MeshIOCubit
[pylithapp.mesh_generator.reader]
-filename = mesh_tri3.exo
+filename = mesh.exo
coordsys.space_dim = 2
# ----------------------------------------------------------------------
@@ -174,23 +174,25 @@ sub_pc_factor_shift_type = nonzero
# ksp_type = preonly
# Convergence parameters.
-ksp_rtol = 1.0e-8
+ksp_rtol = 1.0e-10
ksp_atol = 1.0e-12
ksp_max_it = 300
ksp_gmres_restart = 100
# Linear solver monitoring options.
ksp_monitor = true
-ksp_view = true
+#ksp_view = true
ksp_converged_reason = true
+ksp_error_if_not_converged = true
# Nonlinear solver monitoring options.
snes_rtol = 1.0e-10
-snes_atol = 1.0e-12
+snes_atol = 1.0e-11
snes_max_it = 100
snes_monitor = true
-snes_view = true
+#snes_view = true
snes_converged_reason = true
+snes_error_if_not_converged = true
# PETSc summary -- useful for performance information.
log_summary = true
More information about the CIG-COMMITS
mailing list