Grid Project:   RM SUBCAU Rhino 2010
      with Faults, N-Tunnel, and Cavities

Grid Team Member: Terry Miller
Grid Versions: RM SUBCAU Rhino 2010
Data Source: Rhino from Coblentz
Delivered To: Ed Kwicklis
Date Completed: May 5 2010




Gallery of FEHM Model, Faults, Tunnel, Cavities

Gallery of Mesh Details

Gallery of Material Zones




/scratch/nts/tam/rainier_mesa/HSU_basecase_SubCAU/fehm_octree_rhino_v01
sundog% cat README
This mesh is full octree with refinement down to 12 meters
along the faults and along layers ATCU, RVA, and BRA
The tunnel elements are refined to 6 meters

This mesh is truncated at ground surface.
   The point count is      872032
   The element count is   5028526

Work directories are under  
/scratch/nts/tam/rainier_mesa/HSU_basecase_SubCAU/grid_gen/
Final tet mesh is 
   /grid_gen/rhino_octree_flts_tunnel/tet_zones/tet_rhino.inp

pending - top areas by voronoi method
pending - Cavity Zone files - cavity radius 1x radius (zones equal LOCID in table)  
possible - list of nodes at octree refine interfaces for debugging infil problems

Regular FEHM files:
  tet.fehmn (real coordinates) and tet_zero.fehmn (coorindates at 0,0,0)
  tet_zero.stor
  tet_multi_mat.zone
  tet_outside.zone

Do not use these areas:  tet_outside_median.area (these are median areas, not voronoi)
use replacement file:    tet_outside_voronoi.area

Data file for top nodes of mesh. Includes attribute information including
voronoi volume, material(hsu_id), fault(flt_id), tunnel(nt_bot and nt_top),
the z vector voronoi area (zn_vor_area), and node number (id_node)
==> top_node_attributes.dat <==
# lineno, vor_volume, hsu_id,  layer_id,  flt_id,  nt_bot,  nt_top,  zn_vor_area, id_node 
      1  0.688860000000E+04   20   6    6  0    0  0.288000000000E+03     3534
      2  0.507840000000E+05   16   8    8  0    0  0.110400000000E+04     3565
      3  0.155520000000E+05   20   6    6  0    0  0.288000000000E+03     4118


Material Zone files (zones 4-43, see stats below)
  tet_rhino_material.zone  - HSU ID based on RM numbers 4 through 43
  tet_layer_material.zone  - Rhino layer ID numbered from top 1 to bottom 11

Fault Zone files (zones 100 - 117, see stats below)
  flt101_material.zone  flt105_material.zone  flt109_material.zone  flt113_material.zone
  flt102_material.zone  flt106_material.zone  flt110_material.zone  flt114_material.zone
  flt103_material.zone  flt107_material.zone  flt111_material.zone  flt116_material.zone
  flt104_material.zone  flt108_material.zone  flt112_material.zone  flt117_material.zone

Tunnel Zone files (top is 51, bottom is 52)
tunnel_bot_material.zone  
tunnel_top_material.zone

Cavity Zone files - cavity radius 1x radius (zones equal LOCID in table)  
Initial data is from NTS_nuclear_tests_RM-SM.xls
22 zones with locations and radius computed with corrections in cavity_vers1.dat
Huron (796) and Diamond Ace (795) are simultaneous
Mineral Quarry(802) and Randsburg (803) are simultaneous but seperate drifts
README in /scratch/nts/tam/rainier_mesa/HSU_basecase_SubCAU/grid_gen/N-tunnel_cavities

Cavity zones are written with radius and radius x 2
cav784_MMIST_material.zone	cav784_MMISTx2_material.zone
cav785_HUDSS_material.zone	cav785_HUDSSx2_material.zone
cav786_NMIST_material.zone	cav786_NMISTx2_material.zone
cav787_DMIST_material.zone	cav787_DMISTx2_material.zone
cav788_HUSKY_material.zone	cav788_HUSKYx2_material.zone
cav789_MINGB_material.zone	cav789_MINGBx2_material.zone
cav790_HYBLA_material.zone	cav790_HYBLAx2_material.zone
cav791_MEPIC_material.zone	cav791_MEPICx2_material.zone
cav792_DIABH_material.zone	cav792_DIABHx2_material.zone
cav793_MINIR_material.zone	cav793_MINIRx2_material.zone
cav794_MJADE_material.zone	cav794_MJADEx2_material.zone
cav795_DIACE_material.zone	cav795_DIACEx2_material.zone
cav796_HURLA_material.zone	cav796_HURLAx2_material.zone
cav797_MRAIN_material.zone	cav797_MRAINx2_material.zone
cav798_TOMME_material.zone	cav798_TOMMEx2_material.zone
cav799_BEECH_material.zone	cav799_BEECHx2_material.zone
cav800_MYARD_material.zone	cav800_MYARDx2_material.zone
cav801_MNOTE_material.zone	cav801_MNOTEx2_material.zone
cav802_MQUAR_material.zone	cav802_MQUARx2_material.zone
cav803_RANDS_material.zone	cav803_RANDSx2_material.zone
cav804_MECHO_material.zone	cav804_MECHOx2_material.zone
cav805_HUNTS_material.zone	cav805_HUNTSx2_material.zone

original input
# id   x             y               top     depth    radius	# event name
784	570299.5	4117383.4	2226.9	377.0	36.0	MIDIMIST
785	570423.8	4117679.1	2194.6	344.4	36.0	HUDSON SEAL
787	570531.7	4117294.5	2255.8	402.0	35.0	DIANA MIST
786	570213.7	4117998.6	2226.0	376.1	36.0	MISTY NORTH
788	570712.7	4117196.1	2264.7	413.3	35.0	HUSKY ACE
789	570342.8	4118293.5	2238.6	389.0	35.0	MING BLADE
790	570645.0	4117278.6	2253.6	404.0	35.0	HYBLA FAIR
791	569882.5	4118151.2	2250.7	398.0	35.0	MIGHTY EPIC
792	570032.3	4118123.6	2239.1	388.0	35.0	DIABLO HAWK
793	570508.5	4118400.4	2239.4	389.5	35.0	MINERS IRON
795	570384.6	4118484.1	-99.9	407.3	35.0	DIAMOND ACE
796	570384.6	4118484.1	-99.9	408.5	35.0	HURON LANDING
794	570541.6	4118119.5	1855.0	379.2	35.0	MINI JADE
798	570165.9	4118229.4	-99.9	404.8	35.0	TOMME/MIDNIGHT
797	570358.9	4117242.2	2238.7	388.6	35.0	MISTY RAIN
799	570092.0	4118220.6	-99.9	404.5	35.0	DIAMOND BEECH
800	570533.5	4118168.3	2230.4	375.0	36.0	MILL YARD
801	570226.3	4118284.8	2250.3	398.7	35.0	MIDDLE NOTE
804	570159.3	4117041.9	2259.4	400.2	68.0	MISTY ECHO
802	569724.6	4117907.0	2243.1	389.4	35.0	MINERAL QUARRY
803	569724.6	4117907.0	2243.1	389.4	35.0	RANDSBURG
805	570104.0	4117916.2	1853.9	385.3	35.0	HUNTERS TROPHY



Mesh attributes include:
 vorvol     2.160000000E+02  8.847360000E+05 8.845200000E+05    872032  
 hsu_id                   4               43              39    872032  
 layer_id                 1               12              11    872032  
 elayer_id                1               11              10   5028526  
 itetlev                  0                4               4   5028526  


Files appended with "zero" are at the Rhino model zero coordinates
which have the tunnel at zero. The other files are at real coordinates
which are in UTM meters.

define XTRANS 570500.0
define YTRANS 4117750.0
define ZTRANS 1850.0

ATTRIBUTE NAME              MIN               MAX         DIFFERENCE    LENGTH  
 xic               -1.400000000E+03  1.480000000E+03 2.880000000E+03    872032  
 yic               -1.048000000E+03  9.680000000E+02 2.016000000E+03    872032  
 zic               -5.500000000E+02  4.580000000E+02 1.008000000E+03    872032  
 
trans 1,0,0/0. 0. 0./XTRANS YTRANS ZTRANS                                       

ATTRIBUTE NAME              MIN               MAX         DIFFERENCE    LENGTH  
 xic                5.691000000E+05  5.719800000E+05 2.880000000E+03    872032  
 yic                4.116702000E+06  4.118718000E+06 2.016000000E+03    872032  
 zic                1.300000000E+03  2.308000000E+03 1.008000000E+03    872032  


==== STATS ================================================================


*********dump_material_lists********                                            
Minimum material ID value =      4                                              
Maximum material ID value =     43                                              
Total possible materials  =     40                                              
Layer   11  MGCU    HSU   4 has    33293 nodes. #nodes/nnodes is   0.381786448204E-01 
Layer   10  LCA3    HSU  11 has   135501 nodes. #nodes/nnodes is   0.155385352831     
Layer    9  LCCU1   HSU  12 has    67036 nodes. #nodes/nnodes is   0.768733257495E-01 
Layer    8  ATCU    HSU  16 has    63785 nodes. #nodes/nnodes is   0.731452515504E-01 
Layer    7  RVA     HSU  19 has    46141 nodes. #nodes/nnodes is   0.529120490991E-01 
Layer    6  OSBCU   HSU  20 has   148550 nodes. #nodes/nnodes is   0.170349253238     
Layer    5  LTCU    HSU  21 has   187593 nodes. #nodes/nnodes is   0.215121692782     
Layer    4  BRA     HSU  25 has    29825 nodes. #nodes/nnodes is   0.342017265421E-01 
Layer    3  TMLVTA  HSU  39 has   136234 nodes. #nodes/nnodes is   0.156225918315     
Layer    2  TMWTA   HSU  40 has    22459 nodes. #nodes/nnodes is   0.257547888151E-01 
Layer    1  AA      HSU  43 has     1615 nodes. #nodes/nnodes is   0.185199625702E-02 

*********Fault list********                                            
Fault    101 has     50310 nodes. #nodes/nnodes is   0.576928369601E-01 
Fault    102 has     19515 nodes. #nodes/nnodes is   0.223787659168E-01 
Fault    103 has     34695 nodes. #nodes/nnodes is   0.397863839859E-01 
Fault    104 has     50499 nodes. #nodes/nnodes is   0.579095721258E-01 
Fault    105 has     36032 nodes. #nodes/nnodes is   0.413195846024E-01 
Fault    106 has     41313 nodes. #nodes/nnodes is   0.473755550255E-01 
Fault    107 has     13714 nodes. #nodes/nnodes is   0.157264871014E-01 
Fault    108 has     15305 nodes. #nodes/nnodes is   0.175509614326E-01 
Fault    109 has     34569 nodes. #nodes/nnodes is   0.396418938755E-01 
Fault    110 has     12549 nodes. #nodes/nnodes is   0.143905269531E-01 
Fault    111 has     16089 nodes. #nodes/nnodes is   0.184500110088E-01 
Fault    112 has     23939 nodes. #nodes/nnodes is   0.274519742395E-01 
Fault    113 has     12357 nodes. #nodes/nnodes is   0.141703515467E-01 
Fault    114 has     16346 nodes. #nodes/nnodes is   0.187447249642E-01 
Fault    115 (removed) coincident and equal to 102 
Fault    116 has     39654 nodes. #nodes/nnodes is   0.454731019045E-01 
Fault    117 has     14226 nodes. #nodes/nnodes is   0.163136215185E-01 

*********Tunnel list********                                            
Tunnel top    51 has    9233 nodes. #nodes/nnodes is   0.105879142050E-01 
Tunnel bottom 52 has    9231 nodes. #nodes/nnodes is   0.105856207112E-01 

*********Cavity list********                                            
Cavity  784 has       276 nodes. #nodes/nnodes is   0.316502146710E-03 
Cavity  785 has       207 nodes. #nodes/nnodes is   0.237376610033E-03 
Cavity  787 has       221 nodes. #nodes/nnodes is   0.253431066750E-03 
Cavity  786 has       427 nodes. #nodes/nnodes is   0.489660929874E-03 
Cavity  788 has       198 nodes. #nodes/nnodes is   0.227055887857E-03 
Cavity  789 has       402 nodes. #nodes/nnodes is   0.460992257165E-03 
Cavity  790 has       212 nodes. #nodes/nnodes is   0.243110344575E-03 
Cavity  791 has       418 nodes. #nodes/nnodes is   0.479340207699E-03 
Cavity  792 has       460 nodes. #nodes/nnodes is   0.527503577850E-03 
Cavity  793 has       200 nodes. #nodes/nnodes is   0.229349381674E-03 
Cavity  795 has       333 nodes. #nodes/nnodes is   0.381866720487E-03 
Cavity  796 has       333 nodes. #nodes/nnodes is   0.381866720487E-03 
Cavity  794 has       327 nodes. #nodes/nnodes is   0.374986239037E-03 
Cavity  798 has       299 nodes. #nodes/nnodes is   0.342877325603E-03 
Cavity  797 has       210 nodes. #nodes/nnodes is   0.240816850758E-03 
Cavity  799 has       338 nodes. #nodes/nnodes is   0.387600455029E-03 
Cavity  800 has       368 nodes. #nodes/nnodes is   0.422002862280E-03 
Cavity  801 has       439 nodes. #nodes/nnodes is   0.503421892775E-03 
Cavity  804 has       813 nodes. #nodes/nnodes is   0.932305236505E-03 
Cavity  802 has       364 nodes. #nodes/nnodes is   0.417415874647E-03 
Cavity  803 has       364 nodes. #nodes/nnodes is   0.417415874647E-03 
Cavity  805 has       375 nodes. #nodes/nnodes is   0.430030090639E-03 


*********dump_multi_material_list********                                       
Layer   11  MGCU    HSU   4 has     19282 multi-material connections.                    
Layer   10  LCA3    HSU  11 has     60737 multi-material connections.                    
Layer    9  LCCU1   HSU  12 has     81671 multi-material connections.                    
Layer    8  ATCU    HSU  16 has    139555 multi-material connections.                    
Layer    7  RVA     HSU  19 has     81974 multi-material connections.                    
Layer    6  OSBCU   HSU  20 has    172218 multi-material connections.                    
Layer    5  LTCU    HSU  21 has    171528 multi-material connections.                    
Layer    4  BRA     HSU  25 has     64385 multi-material connections.                    
Layer    3  TMLVTA  HSU  39 has    155429 multi-material connections.                    
Layer    2  TMWTA   HSU  40 has     32899 multi-material connections.                    
Layer    1  AA      HSU  43 has      4246 multi-material connections.                    

*********dump_outside_list********                                              
Face top           1 has     15941 nodes.                                       
top       Area            Area_x          Area_y          Area_z                
          0.6732125E+07   0.1231884E+07   0.1238736E+07   0.5806080E+07         
Face bottom        2 has     12088 nodes.                                       
bottom    Area            Area_x          Area_y          Area_z                
          0.5931985E+07   0.1464240E+06   0.1494240E+06   0.5806080E+07         
Face left_w        3 has      3563 nodes.                                       
left_w    Area            Area_x          Area_y          Area_z                
          0.2650313E+07   0.2208984E+07   0.2557440E+06   0.5838540E+06         
Face right_e       5 has      7909 nodes.                                       
right_e   Area            Area_x          Area_y          Area_z                
          0.3296991E+07   0.2208936E+07   0.5788560E+06   0.1504758E+07         
Face back_n        6 has      6408 nodes.                                       
back_n    Area            Area_x          Area_y          Area_z                
          0.3678084E+07   0.3383040E+06   0.3016566E+07   0.9553740E+06         
Face front_s       4 has      9307 nodes.                                       
front_s   Area            Area_x          Area_y          Area_z                
          0.3748264E+07   0.4316460E+06   0.3015774E+07   0.1026630E+07         


*** 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 =   872032  n connections =    7214560                  
AMatbld3d_stor: Number of written coefs =       689                             
AMatbld3d_stor: Maximum num. connections to a node =         55                 
AMatbld3d_stor: Volume min =   2.7000000E+01                                    
AMatbld3d_stor: Volume max =   8.8473600E+05                                    
AMatbld3d_stor: Total Volume:   4.9753757E+09                                   
AMatbld3d_stor: abs(Aij/xij) min =   0.0000000E+00                              
AMatbld3d_stor: abs(Aij/xij) max =   1.0800000E+02                              
AMatbld3d_stor: (Aij/xij) max =   0.0000000E+00                                 
AMatbld3d_stor: (Aij/xij) min =  -1.0800000E+02                                 
Compression used for graph and coefficient values                               
The area coefficient values were written as scalar area/distance                
ascii STOR file written with name tet_zero.stor                                 
*** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL ***                             



Gallery of Full Model, Faults, Tunnel, Cavities

Gallery of Final Materials

 
 
 




Grid Team Contacts

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