Grid Project:   Springerville for NRAP-ARRA

Grid Team Member: Terry Miller
Grid Versions: Block Version 01
Data Source: Coblentz October 2011
Delivered To: George Zyvoloski
Date Completed: Oct 26 2011




Purpose of Project

Create a simple block mesh based on the Springerville GFM being developed for hydrologic and stress-permeability modeling as part of the National Risk Assessment Partnership NRAP-ARRA project. The model incorporates geologic stratigraphy and fault information.



PDF of GFM Quarterly Report by D. Coblentz


Block Mesh with Fault
Cells are colored by GFM stratigraphy, red nodes show fault location
Block mesh with fault




Project Summary

Dimension: 3D
Grid Element Type: hex and tet
Grid Software Used: LaGriT driven by Cshell script


material 8 = above topography
 ------------------------------ surface 7 = dem 
material 7 = congolmeration of upper units
 ------------------------------ surface 6 = san andreas
material 6 = san andreas
 ------------------------------ surface 5 = glorieta
material 5 = glorieta
 ------------------------------ surface 4 = upper supai
material 4 = upper supai
 ------------------------------- surface 3 = ft. apache
material 3 = ft. apache
 ------------------------------- surface 2 = low supai
material 2 = low supai
 -------------------------------- surface 1 = precambrian
material 1 = precambrian


The Fault zone  HAS  4723 POINTS

hex mesh
  Number of nodes:     44541
  Number of hex elements:  40000

define/XMIN/3.e4
define/YMIN/3.e4
define/ZMIN/2.5e2
define/XMAX/4.e4
define/YMAX/4.e4
define/ZMAX/2.5e3
define/NX/21
define/NY/21
define/NZ/101



Material           8 has     10782 nodes. #nodes/nnodes is   0.242069110000
Material           7 has      3771 nodes. #nodes/nnodes is   0.846635700000E-01
Material           6 has      1323 nodes. #nodes/nnodes is   0.297029700000E-01
Material           5 has       612 nodes. #nodes/nnodes is   0.137401500000E-01
Material           4 has      3009 nodes. #nodes/nnodes is   0.675557330000E-01
Material           3 has       636 nodes. #nodes/nnodes is   0.142789800000E-01
Material           2 has     20522 nodes. #nodes/nnodes is   0.460744020000
Material           1 has      3885 nodes. #nodes/nnodes is   0.872230100000E-01
Material           0 has         1 nodes. #nodes/nnodes is   0.224512240000E-04


tet mesh
  Number of nodes:     44540
  Number of tet elements:  239996

Material           8 has     10782 nodes. #nodes/nnodes is   0.242074530000    
Material           7 has      3771 nodes. #nodes/nnodes is   0.846654700000E-01 
Material           6 has      1323 nodes. #nodes/nnodes is   0.297036400000E-01 
Material           5 has       612 nodes. #nodes/nnodes is   0.137404600000E-01 
Material           4 has      3009 nodes. #nodes/nnodes is   0.675572520000E-01 
Material           3 has       636 nodes. #nodes/nnodes is   0.142793000000E-01 
Material           2 has     20522 nodes. #nodes/nnodes is   0.460754400000     
Material           1 has      3885 nodes. #nodes/nnodes is   0.872249700000E-01 

*** Construct and Compress Sparse Matrix:3D ***                                 
   *** Compress Area Coefficient Values ***                                     
AMatbld3d_stor: *****Zero Negative Coefficients ******                          
AMatbld3d_stor: Number of 'zero' (< 10e-8*max) coefs       0                    
AMatbld3d_stor: npoints =    44540  n connections =     302434                  
AMatbld3d_stor: Number of written coefs =        10                             
AMatbld3d_stor: Maximum num. connections to a node =         11                 
AMatbld3d_stor: Volume min =   7.0312500E+05                                    
AMatbld3d_stor: Volume max =   8.4318085E+06                                    
AMatbld3d_stor: Total Volume:   2.2500000E+11                                   
AMatbld3d_stor: abs(Aij/xij) min =   0.0000000E+00                              
AMatbld3d_stor: abs(Aij/xij) max =   1.1111111E+04                              
AMatbld3d_stor: (Aij/xij) max =   0.0000000E+00                                 
AMatbld3d_stor: (Aij/xij) min =  -1.1111111E+04                                 
Compression used for graph and coefficient values                               
The area coefficient values were written as scalar area/distance                
ascii STOR file written with name mesh_tet.stor        




Project Archive Location


FEHMN files for Block Mesh Version 01:

/scratch/fwo/tam/Springerville/gaz_vers01

-rw-r--r--  1 tam  _update_sharing   6082465 Oct 26 13:48 mesh_hex.fehmn
-rw-r--r--  1 tam  _update_sharing    490131 Oct 26 13:48 mesh_hex_material.zone
-rw-r--r--  1 tam  _update_sharing    134345 Oct 26 13:48 mesh_hex_multi_mat.zone

-rw-r--r--  1 tam  _update_sharing  13002230 Oct 26 13:48 mesh_tet.fehmn
-rw-r--r--  1 tam  _update_sharing   8432689 Oct 26 13:48 mesh_tet.stor
-rw-r--r--  1 tam  _update_sharing     84309 Oct 26 13:48 mesh_tet_interface.zone
-rw-r--r--  1 tam  _update_sharing    490104 Oct 26 13:48 mesh_tet_material.zone
-rw-r--r--  1 tam  _update_sharing    441344 Oct 26 13:48 mesh_tet_multi_mat.zone
-rw-r--r--  1 tam  _update_sharing    103197 Oct 26 13:48 mesh_tet_outside.zone
-rw-r--r--  1 tam  _update_sharing    595311 Oct 26 13:48 mesh_tet_outside_vor.area

Zone file for fault points fault.zone.pset

Mesh viewing files:
-rw-r--r--  1 tam  _update_sharing   4606292 Oct 26 13:48 mesh_hex.gmv
-rw-r--r--  1 tam  _update_sharing  15466725 Oct 26 13:48 mesh_hex.inp
-rw-r--r--  1 tam  _update_sharing  11244044 Oct 26 13:48 mesh_tet.gmv


 
 

Images and Project Details



Gallery of Mesh Images


 
 
 




Grid Team Contacts

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