[cig-commits] r18588 - in short/3D/PyLith/trunk/tests/2d: . refinefaulttip
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 10 16:11:40 PDT 2011
Author: brad
Date: 2011-06-10 16:11:39 -0700 (Fri, 10 Jun 2011)
New Revision: 18588
Added:
short/3D/PyLith/trunk/tests/2d/refinefaulttip/
short/3D/PyLith/trunk/tests/2d/refinefaulttip/pylithapp.cfg
short/3D/PyLith/trunk/tests/2d/refinefaulttip/quad4.mesh
Log:
Added refinefaulttip for testing refinement around fault tip.
Added: short/3D/PyLith/trunk/tests/2d/refinefaulttip/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/refinefaulttip/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/refinefaulttip/pylithapp.cfg 2011-06-10 23:11:39 UTC (rev 18588)
@@ -0,0 +1,147 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+meshiocubit = 1
+implicitelasticity = 1
+quadrature2d = 1
+fiatlagrange = 1
+faultcohesivedyn = 1
+
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+implicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOAscii
+debug = 1
+distributor.write_partition = True
+distributor.data_writer.filename = partition.vtk
+
+# SWITCH FOR REFINEMENT
+refiner = pylith.topology.RefineUniform
+
+
+[pylithapp.mesh_generator.reader]
+filename = quad4.mesh
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+formulation = pylith.problems.Implicit
+
+normalizer = spatialdata.units.NondimElasticQuasistatic
+normalizer.length_scale = 1.0*m
+normalizer.relaxation_time = 1.0*s
+
+bc = [xpos,xneg]
+interfaces = [fault]
+materials = [elastic]
+
+[pylithapp.timedependent.implicit.time_step]
+total_time = 0.0*s
+dt = 0.1*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 0
+db_properties = spatialdata.spatialdb.UniformDB
+db_properties.label = Elastic properties
+db_properties.values = [density, vp, vs]
+db_properties.data = [2500.0*km/m**3, 5.1962*km/s, 3.0*km/s]
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc]
+
+# Dirichlet BC on -x
+[pylithapp.timedependent.bc.xneg]
+bc_dof = [0, 1]
+label = xneg
+
+
+# Dirichlet BC on +x
+[pylithapp.timedependent.bc.xpos]
+bc_dof = [0, 1]
+label = xpos
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+id = 100
+label = fault
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.label = Final slip
+slip.values = [left-lateral-slip,fault-opening]
+slip.data = [0.1*m, 0.0*m]
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip initiation time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*s]
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+# Preconditioner settings.
+pc_type = asm
+sub_pc_factor_shift_type = nonzero
+
+# KSP
+ksp_rtol = 1.0e-12
+ksp_atol = 1.0e-15
+ksp_max_it = 500
+ksp_gmres_restart = 100
+
+ksp_monitor = true
+#ksp_view = true
+ksp_converged_reason = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+vertex_data_fields = [displacement]
+writer.filename = refinefaulttip.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_data_fields=[slip]
+writer.filename = refinefaulttip-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = refinefaulttip-statevars.vtk
Added: short/3D/PyLith/trunk/tests/2d/refinefaulttip/quad4.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/refinefaulttip/quad4.mesh (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/refinefaulttip/quad4.mesh 2011-06-10 23:11:39 UTC (rev 18588)
@@ -0,0 +1,101 @@
+// Original mesh
+//
+// 2 ----- 5 ----- 8 -----11 -----14
+// | | | | |
+// | 4 | 5 | 6 | 7 |
+// | | | | |
+// 1 ----- 4 ----- 7 -----10 -----13
+// | | | | |
+// | 0 | 1 | 2 | 3 |
+// | | | | |
+// 0 ----- 3 ----- 6 ----- 9 -----12
+//
+mesh = {
+ dimension = 2
+ use-index-zero = true
+ vertices = {
+ dimension = 2
+ count = 15
+ coordinates = {
+ 0 -2.0 -1.0
+ 1 -2.0 0.0
+ 2 -2.0 1.0
+ 3 -1.0 -1.0
+ 4 -1.0 0.0
+ 5 -1.0 1.0
+ 6 0.0 -1.0
+ 7 0.0 0.0
+ 8 0.0 1.0
+ 9 +1.0 -1.0
+ 10 +1.0 0.0
+ 11 +1.0 1.0
+ 12 +2.0 -1.0
+ 13 +2.0 0.0
+ 14 +2.0 1.0
+ }
+ }
+ cells = {
+ count = 8
+ num-corners = 4
+ simplices = {
+ 0 0 3 4 1
+ 1 3 6 7 4
+ 2 6 9 10 7
+ 3 9 12 13 10
+ 4 1 4 5 2
+ 5 4 7 8 5
+ 6 7 10 11 8
+ 7 10 13 14 11
+ }
+ material-ids = {
+ 0 0
+ 1 0
+ 2 0
+ 3 0
+ 4 0
+ 5 0
+ 6 0
+ 7 0
+ }
+ }
+ group = {
+ name = fault
+ type = vertices
+ count = 2
+ indices = {
+ 7 8
+ }
+ }
+ group = {
+ name = ypos
+ type = vertices
+ count = 5
+ indices = {
+ 2 5 8 11 14
+ }
+ }
+ group = {
+ name = yneg
+ type = vertices
+ count = 5
+ indices = {
+ 0 3 6 9 12
+ }
+ }
+ group = {
+ name = xpos
+ type = vertices
+ count = 3
+ indices = {
+ 12 13 14
+ }
+ }
+ group = {
+ name = xneg
+ type = vertices
+ count = 3
+ indices = {
+ 0 1 2
+ }
+ }
+}
More information about the CIG-COMMITS
mailing list