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 |