RMSM Slice Trv2a 20m x 5m orthogonal mesh Source EarthVision HFM: From Tim Vogt (Navarro-HGL) TDR UGTA-4-1542 https://ugta.nv.doe.gov/sites/Document%20Center/Documents/RMSM_Prime_for%20TDR_20160624.zip Local: /scratch/nts/ev_models/RM_HFM_Prime_Nov_2016/RMSM_Prime/rmsmv2 Source Trv 2A fence line points (UTM NAD27): From Tim Vogt (Navarro-HGL) 575809.00 4125783.000 575165.00 4121503.000 572948.00 4119385.000 571721.00 4118975.000 570423.00 4117679.000 569075.93 4116988.994 568283.06 4115924.949 563113.05 4113274.550 Main Grid directory /scratch/nts/tam/rainier_mesa/RMSM/slice_2A_20x5m /FEHM_tet are FEHM files (tet) /EXO_hex are AMANZI ExodusII files (hex) Straight version 3D Slice along x-axis (3 cells 20m each): xic 0.000000000E+00 1.917600000E+04 1.917600000E+04 yic 0.000000000E+00 6.000000000E+01 6.000000000E+01 zic 0.000000000E+00 2.265000000E+03 2.265000000E+03 Material Zones 4-39 Fault Zones 51-54 Block HFM ZONE 1 MGCU 4 2 UCCU 9 3 LCA3 11 4 LCCU1 12 5 ATCU 14 6 LTCU1 15 7 RVA 17 8 OSBCU 18 9 TUBA 19 10 LTCU 20 11 LVTA1 21 12 BRA 22 13 LVTA2 24 14 UTCU1 25 15 SWA 26 16 PVTA 35 17 TCA 36 18 TMLVTA 38 19 TMWTA 39 20 Big Burn Valley Fault 51 21 RM1 Fault 52 22 NW T-Tunnel Fault 53 23 Mid T-Tunnel Fault 54 -------------------------------------------------------------- Hex Mesh for Amanzi:-------------------------------------------- /scratch/nts/tam/rainier_mesa/RMSM/slice_2A_20x5m/EXO_hex nodes = 1534188 hexs = 1145760 min volume = 3.3092105E+02 max volume = 6.6790123E+02 total volume = 0.2286049E+10 Material list for hex cells: Block Zone Hex Cells Volume Fractional Volume 1 4 43128 0.8589610E+08 0.037574041 2 9 260943 0.5208737E+09 0.227848881 3 11 178599 0.3558855E+09 0.155677135 4 12 82485 0.1645901E+09 0.071997630 5 14 21312 0.4248981E+08 0.018586569 6 15 10047 0.2005170E+08 0.008771335 7 17 95253 0.1901038E+09 0.083158242 8 18 141570 0.2826412E+09 0.123637437 9 19 6309 0.1258270E+08 0.005504124 10 20 81072 0.1616810E+09 0.070725101 11 21 2985 0.5936479E+07 0.002596829 12 22 70983 0.1417676E+09 0.062014256 13 24 5118 0.1021594E+08 0.004468818 14 25 27405 0.5485641E+08 0.023996166 15 26 5403 0.1078458E+08 0.004717561 16 35 2505 0.4995521E+07 0.002185221 17 36 561 0.1119589E+07 0.000489748 18 38 71526 0.1426391E+09 0.062395451 19 39 33345 0.6655354E+08 0.029112910 20 51 1164 0.2324000E+07 0.001016601 21 52 1347 0.2685215E+07 0.001174610 22 53 1371 0.2729345E+07 0.001193913 23 54 1329 0.2645918E+07 0.001157420 Outside side sets: sideset 1 top: 4716 sideset 2 bottom: 2883 sideset 3 left_w: 1047 sideset 4 front_s: 381920 sideset 5 right_e: 1110 sideset 6 back_n: 381920 Exodus II Mesh File: -rw-r--r-- 1 tamiller sft 88849460 Jan 3 2017 mesh_fsets.exo -------------------------------------------------------------- Tet Mesh for FEHM:-------------------------------------------- /scratch/nts/tam/rainier_mesa/RMSM/slice_2A_20x5m/FEHM_tet nodes = 1534188 tets = 6807154 min volume = 3.3092105E+02 max volume = 6.6790123E+02 total volume = 2.2860489E+09 Material list for tet nodes: 4 has 58040 nodes. #nodes/nnodes is 0.378310866654E-01 9 has 349776 nodes. #nodes/nnodes is 0.227987706661 11 has 236932 nodes. #nodes/nnodes is 0.154434785247 12 has 109772 nodes. #nodes/nnodes is 0.715505555272E-01 14 has 28432 nodes. #nodes/nnodes is 0.185322798789E-01 15 has 13336 nodes. #nodes/nnodes is 0.869254581630E-02 17 has 126712 nodes. #nodes/nnodes is 0.825922265649E-01 18 has 188204 nodes. #nodes/nnodes is 0.122673362494 19 has 8408 nodes. #nodes/nnodes is 0.548042356968E-02 20 has 107208 nodes. #nodes/nnodes is 0.698793083429E-01 21 has 3972 nodes. #nodes/nnodes is 0.258899177425E-02 22 has 94780 nodes. #nodes/nnodes is 0.617786087096E-01 24 has 6820 nodes. #nodes/nnodes is 0.444534840062E-02 25 has 36592 nodes. #nodes/nnodes is 0.238510537893E-01 26 has 7264 nodes. #nodes/nnodes is 0.473475223407E-02 35 has 3428 nodes. #nodes/nnodes is 0.223440676928E-02 36 has 792 nodes. #nodes/nnodes is 0.516233965755E-03 38 has 95104 nodes. #nodes/nnodes is 0.619897954166E-01 39 has 44680 nodes. #nodes/nnodes is 0.291228983551E-01 51 has 3112 nodes. #nodes/nnodes is 0.202843453735E-02 52 has 3600 nodes. #nodes/nnodes is 0.234651821665E-02 53 has 3668 nodes. #nodes/nnodes is 0.239084125496E-02 54 has 3556 nodes. #nodes/nnodes is 0.231783860363E-02 Outside face nodes (default with double on top stair-step): Face top 1 has 5580 nodes. Face bottom 2 has 3848 nodes. Face left_w 3 has 3532 nodes. Face right_e 5 has 3528 nodes. Face front_s 4 has 383547 nodes. Face back_n 6 has 383547 nodes. Outside face nodes (single on top stair-step): Face top 1 has 3848 nodes. Face bottom 2 has 3848 nodes. Face left_w 3 has 1816 nodes. Face right_e 5 has 1816 nodes. Face front_s 4 has 383547 nodes. Face back_n 6 has 383547 nodes. FEHM Files: -rw-r--r-- 1 tamiller sft 388020721 Dec 16 11:48 tet.fehmn -rw-r--r-- 1 tamiller sft 281106612 Dec 16 11:52 tet.stor -rw-r--r-- 1 tamiller sft 16877763 Dec 16 11:48 tet_material.zone -rw-r--r-- 1 tamiller sft 8619659 Dec 16 11:52 tet_outside.zone -rw-r--r-- 1 tamiller sft 49758111 Dec 16 11:52 tet_outside_vor.area -rw-r--r-- 1 tamiller sft 8562899 Dec 16 11:54 tet_sngl_outside.zone -rw-r--r-- 1 tamiller sft 49430451 Dec 16 11:54 tet_sngl_outside_vor.area -rw-r--r-- 1 tamiller sft 5947218 Dec 16 11:49 tet_multi_mat.zone -rw-r--r-- 1 tamiller sft 26803446 Dec 16 11:48 tet_interface.zone Method:------------------------------------------------------ Create quad surfaces along each section of tranverse 2A Fence section lengths (set NX = 2 to get section lengths) Divide equally into cells near 20 meter width Divide vertically into cells with 5m height grep max outx3dgen | grep length | grep = min max edge length = 0.1000E+03 max max edge length = 0.5810E+04 min max edge length = 0.1000E+03 max max edge length = 0.1335E+04 min max edge length = 0.1000E+03 max max edge length = 0.1509E+04 min max edge length = 0.1000E+03 max max edge length = 0.1834E+04 min max edge length = 0.1000E+03 max max edge length = 0.1294E+04 min max edge length = 0.1000E+03 max max edge length = 0.3066E+04 min max edge length = 0.1000E+03 max max edge length = 0.4328E+04 Join sections into fence line, this is line A A' Interpolate EV HFM RMSM PRIME zone materials on to fence. Intersect fault surfaces on quad fence. Copy materials on to quad surface along x-axis. Extrude quads into hex mesh and merge into slice 3 cells, 4 nodes wide. Truncate by removing air material. Write FEHM files. Write Exodus files. Done. ====================================================================== NOTES ====================================================================== From Tim Vogt traverse_2a_from_sections_001.trv # Field: 3 traverse non-numeric # Coordinate_System_Id: CRD189B1 # Coordinate_System_Name: NAD27 / UTM Zone 11 North # Units: meters 575809. 4125783. "2a" 575165. 4121503. "2a" 574747. 4120981. "2a" "Diamond Fortune" 574008. 4120384. "2a" 572948. 4119385. "2a" "Mission Ghost" 572477. 4119201. "2a" "Mighty Oak" 571721. 4118975. "2a" 570423. 4117679. "2a" "Hudson Seal" 569075.93 4116988.994 "2a" 568283.06 4115924.949 "2a" 567405.9 4115418.816 "2a" 563113.05 4113274.55 "2a" Begin forwarded message: From: "Vogt, Tim (Navarro-HGL)" Subject: RE: [EXTERNAL] Request for RM HFM Prime and Cross Section coordinates Date: November 29, 2016 at 1:46:55 PM MST To: "Birdsell, Kay" Cc: "Kwicklis, Edward Michael" , "Miller, Terry Ann" , "Kelkar, Sharad (Navarro-HGL)" , "Knop, Matthew (Navarro)" Kay, I’ve attached the traverse file for the cross-section. The traverse is longer than the section but that will be obvious. I believe the section you sent is cut from the ‘Prime’ version and of course anything you create from now on will be from the current version. If you have any questions or issues let me know. Tim From: Birdsell, Kay [mailto:khb@lanl.gov] Sent: Tuesday, November 29, 2016 12:42 PM To: Vogt, Tim (Navarro-HGL) Cc: Kwicklis, Edward Michael ; Miller, Terry Ann ; Kelkar, Sharad (Navarro-HGL) ; Knop, Matthew (Navarro) Subject: RE: [EXTERNAL] Request for RM HFM Prime and Cross Section coordinates Thanks, Tim! I believe it is ok if the cross section we use is based on HFM’ even if the figure is not. Thanks for looking for the coordinates for us. -- Kay H. Birdsell EES-16, MS T003 Los Alamos National Laboratory Los Alamos, NM 87545 office: (505) 665-0260 cell: (505) 690-5451 From: Vogt, Tim (Navarro-HGL) [mailto:Tim.Vogt@nv.doe.gov] Sent: Tuesday, November 29, 2016 12:36 PM To: Birdsell, Kay Cc: Kwicklis, Edward Michael ; Miller, Terry Ann ; Kelkar, Sharad (Navarro-HGL) ; Knop, Matthew (Navarro) Subject: RE: [EXTERNAL] Request for RM HFM Prime and Cross Section coordinates Kay, All of the EarthVision files for RMSM Prime are in the TDR as : UGTA-4-1542. https://ugta.nv.doe.gov/sites/Document%20Center/Documents/RMSM_Prime_for%20TDR_20160624.zip There is currently no additional documentation describing the ‘Prime’ version except what is in the documentation package (above). Although I believe there is some documentation being prepared describing the rational for the changes incorporated in the “Prime’ version. The cross section was done quite a bit earlier. I know I have the coordinates but I’m going to have to look for the them. Also we need to verify that the cross section you included was prepared using the most recent version. You can get started on the download and I’ll find the section coordinates and send them to you. Tim Tim Vogt, CPG UGTA Modeler HGL – Navarro Supporting DOE/NNSA Cell: 702-994-2489 Office: 702-295-2057 From: Birdsell, Kay [mailto:khb@lanl.gov] Sent: Monday, November 28, 2016 12:15 PM To: Vogt, Tim (Navarro-HGL) ; Knop, Matthew (Navarro) Cc: Kwicklis, Edward Michael ; Miller, Terry Ann ; Kelkar, Sharad (Navarro-HGL) Subject: [EXTERNAL] Request for RM HFM Prime and Cross Section coordinates Hi Tim and Matt, LANL will be doing some 2D calculations using the cross section A-A’ shown in the attached figure. This figure comes from the Draft chapter for the “Rainier Mesa/Shoshone Mountain CAU Flow and Transport Model, NNSS, Nye County, Nevada” dated Nov 6, 2014 (v7). Authors: Fenelon, Halford, Carle. To do the work, we need the RM HFM Prime Earthvision files and the coordinates that define the fence diagram A-A’. Will you please help us get these files and coordinates? Terry Miller will be making the grid. Tim – I believe you’re familiar with what Terry needs. Thanks very much, Kay -- Kay H. Birdsell EES-16, MS T003 Los Alamos National Laboratory Los Alamos, NM 87545 office: (505) 665-0260 cell: (505) 690-5451