Grid Project:   Reverse Fault Benchmark

Grid Team Member: Carl Gable
Grid Versions: v01
Data Source: Benchmark Definitions:
http://www.geodynamics.org:8080/cig/workinggroups/eq/workarea/benchmarks
http://www.geodynamics.org:8080/cig/workinggroups/eq/workarea/benchmarks/benchmark-rs-nog/description-rs-nog
Delivered To: SCEC Community
Date Completed: July 2006, update Sept 2006 (add GeoFEST mesh files)

Reverse Fault Benchmark, Variable Resolution, Slip Distribution
Reverse Fault Benchmark, 1000m resolution, slip distribution
Close up view of variable resolution mesh. The finest cells have edge length of ~125m. The fault tip has increased resolution. The slip distribution is mapped onto fault nodes as a floating point attribute.
Cut away view of uniform resolution, 1000m mesh.  The slip distribution is mapped onto fault nodes as a floating point attribute.
no comment
no comment
no comment
no comment



Purpose of Project

Build grids of various resolution (1000m, 500m, 250m, variable) for reverse fault benchmark calculation.

Project Summary

 
Dimension: 3D


Grid Element Type: tetrahedral


Grid Software Used: LaGriT


Uniform Grids



Uniform node spacing of 3 different resolutions. Like a finite difference mesh node distribution however the points are connected into tetrahedral elements.
    tet_1000m_res   tet_0500m_res  tet_0250m_res
Number of Nodes: 15625 117649 912673
Number of Elements: 82944 663552 4962216
Refine Faults, No Smoothing



This starts with a uniform mesh of 1000 m resolution and then refines just the elements adjacent to the fault to 125 m (res_03).
    tet_1000m_res   tet_0500m_res  tet_0250m_res
Number of Nodes: 22973 52793 160100
Number of Elements: 121918 287347 888750
Refine Faults, With Smoothing



This starts with a uniform mesh of 1000 m resolution and then refines just the elements adjacent to the fault to 125 m (res_03). In addition, smoothing is applied to interior nodes to produce a more isotropic mesh.
    tet_var_res_01   tet_var_res_02  tet_var_res_03
Number of Nodes: 22973 59093 187746
Number of Elements: 121918 321645 1040841

Project Archive Files:

Output files for calculations:

UPDATE: Sept 2006, Add GeoFEST mesh files for all cases.

Individual GeoFEST mesh files:
tet_1000m_res.gft.gz
tet_0500m_res.gft.gz
tet_0250m_res.gft.gz
tet_var_res_01_smooth.gft.gz
tet_var_res_02_smooth.gft.gz
tet_var_res_03_smooth.gft.gz
tet_var_res_01_nosmooth.gft.gz
tet_var_res_02_nosmooth.gft.gz
tet_var_res_03_nosmooth.gft.gz

All GeoFEST mesh files in a single tar file:
tar_geofest_all_cases_gzipped.tar

Uniform Grids

tet_0250m_res.gmv
tet_0250m_res.inp
tet_0250m_res.table
tet_0250m_res_outside.zone
tar_0250all.tar.gz

tet_0500m_res.gmv
tet_0500m_res.inp
tet_0500m_res.table
tet_0500m_res_outside.zone
tar_0500all.tar.gz

tet_1000m_res.gmv
tet_1000m_res.inp
tet_1000m_res.table
tet_1000m_res_outside.zone
tar_1000all.tar.gz

---------------------------------------------

Refine Elements Bounding Fault, NO smoothing of node positions.
tet_var_res_01.gmv
tet_var_res_01.inp
tet_var_res_01.table
tet_var_res_01_outside.zone
tar_bm_reverse_var_res_no_smooth_01_all.tar.gz

tet_var_res_02.gmv
tet_var_res_02.inp
tet_var_res_02.table
tet_var_res_02_outside.zone
tar_bm_reverse_var_res_no_smooth_02_all.tar.gz

tet_var_res_03.gmv
tet_var_res_03.inp
tet_var_res_03.table
tet_var_res_03_outside.zone
tar_bm_reverse_var_res_no_smooth_03_all.tar.gz

---------------------------------------------

Refine Elements Bounding Fault, with smoothing of node positions for an isotropic mesh.
tet_var_res_01.gmv
tet_var_res_01.inp
tet_var_res_01.table
tet_var_res_01_outside.zone
tar_bm_reverse_var_res_smoothed_01_all.tar.gz

tet_var_res_02.gmv
tet_var_res_02.inp
tet_var_res_02.table
tet_var_res_02_outside.zone
tar_bm_reverse_var_res_smoothed_02_all.tar.gz

tet_var_res_03.gmv
tet_var_res_03.inp
tet_var_res_03.table
tet_var_res_03_outside.zone
tar_bm_reverse_var_res_smoothed_03_all.tar.gz

---------------------------------------------


123315569 tar_0250all.tar.gz
15897375 tar_0500all.tar.gz
1920735 tar_1000all.tar.gz

2979898 tar_bm_reverse_var_res_no_smooth_01_all.tar.gz
7441716 tar_bm_reverse_var_res_no_smooth_02_all.tar.gz
24413184 tar_bm_reverse_var_res_no_smooth_03_all.tar.gz

2923133 tar_bm_reverse_var_res_smoothed_01_all.tar.gz
8322019 tar_bm_reverse_var_res_smoothed_02_all.tar.gz
28808343 tar_bm_reverse_var_res_smoothed_03_all.tar.gz

What is in the files?
*.gmv - this is a standard GMV file. Mostly for graphics if you use GMV to view grids.
*.inp - this is an AVS format file of the node coordinate list and the tet connectivity.
*.table - this is a table of node attributes for each mesh.

The first column is the node number.
The second column is a node attribute, 'if_fault'. It is 0 for nodes not on the fault
plane, it is > 0 for nodes on the fault plane.
The third column is a node attribute, 'slip', that is 1.0 and tapers to 0.0 according to the specifications of the benchmark description.
The fourth column is a node attribute, 'top', = 1 if a top node, = -5 if not a top node
The fifth column is a node attribute, 'bottom', = 2 if a bottom node, = -5 if not a bottom node
The sixth column is a node attribute, 'left_w', = 3 if a left side node, = -5 if not a left side node
The seventh column is a node attribute, 'right_e', = 5 if a right side node, = -5 if not a right side node
The eight column is a node attribute, 'back_n', = 6 if a back side node, = -5 if not a back side node
The ninth column is a node attribute, 'front_s', = 4 if a front side node, = -5 if not a front side node

*.zone - this is the same information as provide in the *.table file for the top, bottom, left, ... except
the list is not all nodes. There is a list of just top nodes (by node number), a list of just bottom nodes, etc.

NOTE: For *.table and *.zone there are some nodes that appear in more than one list. For example a corner node
may appear in the top,left and front list.

LaGriT Control Files to Build Mesh:

Uniform Grids:
tet_uniform_grid_4_res.lgi
identify_fault_nodes.mlgi

Refine Faults, No Smoothing:
tet_refine_flt_nosmooth.lgi
geometry.mlgi
distance_field.lgi
identify_fault_nodes.mlgi

Refine Faults, With Smoothing:
tet_refine_flt_smooth.lgi
geometry.mlgi
distance_field.lgi
identify_fault_nodes.mlgi


Work Directory:

/scratch/ymp/gable/grid_gen/crustal_deformation/bm_reverse/bm_reverse_uniform/ver_01
/scratch/ymp/gable/grid_gen/crustal_deformation/bm_reverse/bm_reverse_unstructured/refine_fault_only
/scratch/ymp/gable/grid_gen/crustal_deformation/bm_reverse/bm_reverse_unstructured/with_smoothing
 

Images and Project Details


Grid Statistics

Gallery of Project Images

 




Grid Team Contacts

Carl Gable gable -at- lanl -dot- gov 505-665-3533
Terry Miller tamiller -at- lanl -dot- gov 505-667-8009