Grid Project:   Reverse Fault Benchmark  
                     
            
                    
           
             
               |   
 |   
 | 
             
               | 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. 
 | 
        
          |   
 |   
 | 
        
          |   
 |   
 | 
                                   
           
         
                     
 
 
                               
 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 |