[cig-commits] r19419 - in short/3D/PyLith/branches/v1.6-stable/tests/3d: . refine
brad at geodynamics.org
brad at geodynamics.org
Mon Jan 23 15:40:44 PST 2012
Author: brad
Date: 2012-01-23 15:40:43 -0800 (Mon, 23 Jan 2012)
New Revision: 19419
Added:
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/createbc.jou
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/geometry.jou
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.cfg
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.exo.gz
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.jou
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/output/
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/pylithapp.cfg
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.cfg
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.exo.gz
short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.jou
Log:
Added global refinement test.
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/createbc.jou
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/createbc.jou (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/createbc.jou 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,77 @@
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 3
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 6
+nodeset 11 group face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 4
+nodeset 12 group face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 12
+group "face_ypos" add node in surface 14
+nodeset 13 group face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 10
+group "face_yneg" add node in surface 17
+nodeset 14 group face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 11
+group "face_zneg" add node in surface 15
+nodeset 15 group face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 9
+group "face_zpos" add node in surface 16
+nodeset 17 group face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+nodeset 20 group fault
+nodeset 20 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault (interior)
+# ----------------------------------------------------------------------
+group "fault_interior" add node in fault_surface
+group "fault_interior" remove node in face_ypos
+group "fault_interior" remove node in face_yneg
+group "fault_interior" remove node in face_zneg
+nodeset 21 group fault_interior
+nodeset 21 name "fault_interior"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face w/o fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 group face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/geometry.jou
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/geometry.jou (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/geometry.jou 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,33 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 6 km x 6 km x 3 km
+# -3 km <= x <= 3 km
+# -3 km <= y <= 3 km
+# -3 km <= z <= 0 km
+reset
+brick x 6000 y 6000 z 3000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -1500
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+volume 1 3 name "elastic"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+delete body 2
+
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.cfg 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,41 @@
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.implicit.time_step]
+total_time = 0.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8-statevars.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/hex8-fault.vtk
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.jou
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.jou (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/hex8.jou 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,26 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all scheme map
+volume all size 500
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8.exo" dimension 3 overwrite
+
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/pylithapp.cfg 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,146 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+fiatlagrange = 1
+fiatsimplex = 1
+pylithapp = 1
+materials = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOCubit
+reorder_mesh = True
+
+refiner = pylith.topology.RefineUniform
+
+distributor.partitioner = parmetis
+distributor.write_partition = True
+distributor.data_writer.filename = output/tet4_partition.vtk
+
+[pylithapp.mesh_generator.reader]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s ; total time of simulation
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties = spatialdata.spatialdb.UniformDB
+db_properties.label = Elastic material properties
+db_properties.values = [density, Vs, Vp]
+db_properties.data = [2500.0*kg/m**3, 3000.0*m/s, 5291.502622129181*m/s]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_pos,x_neg,z_neg]
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0, 1]
+label = face_xpos
+db_initial.label = Dirichlet BC on +x
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0, 1]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg_nofault
+db_initial.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+interfaces.fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 10
+label = fault
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.label = Final slip
+slip.values = [left-lateral-slip, reverse-slip, fault-opening]
+slip.data = [1.0*m, 0.0*m, 0.0*m]
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip time
+slip_time.values = [slip-time]
+slip_time.data = [0*s]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation]
+output = [domain]
+
+[pylithapp.problem.formulation.output.domain]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_data_fields = [slip,traction_change]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+
+# Preconditioner settings.
+pc_type = asm
+sub_pc_factor_shift_type = nonzero
+
+# Convergence parameters.
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 100
+ksp_gmres_restart = 50
+
+# Linear solver monitoring options.
+ksp_monitor = true
+#ksp_view = true
+ksp_converged_reason = true
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.cfg 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,41 @@
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.implicit.time_step]
+total_time = 0.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4-statevars.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tet4-fault.vtk
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.jou
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.jou (rev 0)
+++ short/3D/PyLith/branches/v1.6-stable/tests/3d/refine/tet4.jou 2012-01-23 23:40:43 UTC (rev 19419)
@@ -0,0 +1,26 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all scheme tetmesh
+volume all size 500
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4.exo" dimension 3 overwrite
+
More information about the CIG-COMMITS
mailing list