V2 half mesh with 5 heaters in salt pile in room Same setup and refinement as V1 mesh extended on other side of zero along X axis Mesh directory /scratch/sft/tam/SALT/SDDI_meshV2_half Nodes: 258800 Tets : 1476328 Mesh material zones: Material Heater 1 1 has 481 nodes. #nodes/nnodes is 0.185857806355E-02 Material Heater 2 2 has 481 nodes. #nodes/nnodes is 0.185857806355E-02 Material Heater 3 3 has 481 nodes. #nodes/nnodes is 0.185857806355E-02 Material Heater 4 4 has 481 nodes. #nodes/nnodes is 0.185857806355E-02 Material Heater 5 5 has 481 nodes. #nodes/nnodes is 0.185857806355E-02 Material Salt 6 has 68526 nodes. #nodes/nnodes is 0.264783620834 Material Air 7 has 91441 nodes. #nodes/nnodes is 0.353326886892 Material Room 8 has 84820 nodes. #nodes/nnodes is 0.327743440866 Material Outside 9 has 11608 nodes. #nodes/nnodes is 0.448531694710E-01 FEHM Files: -rw-r--r-- 1 tamiller sft 78904307 Aug 1 2016 grid_half.fehmn -rw-r--r-- 1 tamiller sft 2847109 Aug 1 2016 grid_half_material.zone -rw-r--r-- 1 tamiller sft 5102606 Aug 1 2016 grid_half_multi_mat.zone -rw-r--r-- 1 tamiller sft 119068 Aug 1 2016 grid_half_outside.zone -rw-r--r-- 1 tamiller sft 686518 Aug 1 2016 grid_half_outside_vor.area -rw-r--r-- 1 tamiller sft 55709005 Aug 1 2016 grid_half.stor Mesh viewing files: -rw-r--r-- 1 tamiller sft 57482988 Aug 1 2016 grid_tet_half.gmv -rw-r--r-- 1 tamiller sft 123353207 Aug 1 2016 grid_tet_half.inp -rw-r--r-- 1 tamiller sft 5924260 Aug 1 2016 tet_surfmesh.gmv heater centers left to right define / X_H1 / -1.8288 define / X_H2 / -.9144 define / X_H3 / 0. define / X_H4 / 0.9144 define / X_H5 / 1.8288 mine-run salt pile covering heaters with slope X end of mine-run salt slope 4.6507272561576825 X upper top of slope 3.5171272561576825 X Past end of room 13. Y back wall of room 2.4384 Z top of crushed salt 1.8288 Y extent of EDZ into wall 3.9624 Z extent of EDZ below floor -1.542 Z extent of EDZ above ceiling 4.572 AMatbld3d_stor: *****Zero Negative Coefficients ****** AMatbld3d_stor: Number of 'zero' (< compress_eps) coefs 1 AMatbld3d_stor: npoints = 258800 ncoefs = 2086254 AMatbld3d_stor: Number of unique coefs = 108 AMatbld3d_stor: Maximum num. connections to a node = 36 AMatbld3d_stor: Volume min = 1.1061268E-04 AMatbld3d_stor: Volume max = 1.8122782E+00 AMatbld3d_stor: Total Volume: 1.6310504E+04 AMatbld3d_stor: abs(Aij/xij) min = 1.3360982E-15 AMatbld3d_stor: abs(Aij/xij) max = 1.3716000E+00 AMatbld3d_stor: (Aij/xij) max = -1.3360982E-15 AMatbld3d_stor: (Aij/xij) min = -1.3716000E+00 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name grid_half.stor Face top 1 has 336 nodes. top Sum Voronoi Area_x Area_y Area_z 0.2229673E+02 0.2972897E+02 0.4459346E+03 Face bottom 2 has 336 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2229673E+02 0.2972897E+02 0.4459346E+03 Face left_w 3 has 1245 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.6689019E+03 0.4100943E+02 0.2229673E+02 Face right_e 5 has 1245 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.6689019E+03 0.4100943E+02 0.2229673E+02 Face back_n 6 has 651 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.4459346E+02 0.8918692E+03 0.2972897E+02 Face front_s 4 has 6988 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.3742541E+02 0.8918692E+03 0.2972897E+02 RUNS: Expand the Level 0 base hex grid lagrit < grid_L0.lgi cp outx3dgen grid_L0.out.txt Create geometry for half-size mesh copy geometry.mlgi and edit into geometry_half.mlgi lagrit < test_geometry.lgi Refine and write FEHM files edit build_eltsets.mlgi to use new geom defs lagrit < refine_main.lgi cp outx3dgen refine_main.out.txt ====== V1 quarter mesh version ========================== Start with mesh3d from Dylan which uses mesh3d from Carl Gable. /scratch/er/dharp/projects/salt/SDDI/mesh3d13 uses grid_L0.inp -> /home/gable/people/harp/mesh3d02/grid_L0.inp refine_main.lgi geometry.mlgi refine_ESET.mlgi build_eltsets.mlgi ====== V0 quarter mesh version ========================== FROM 00README in /home/gable/people/harp/mesh3d02 Build a fine regular mesh and set material (itetclr) properties lagrit < test_geometry.lgi INPUT: grid_L0.inp Refine a mesh to conform to geometry boundaries and use the fine mesh to set itetclr values lagrit < refine_interp.lgi Necessary aux. input files infile / box_pset.mlgi infile / refine_pset.mlgi infile / build_inter.mlgi infile / geometry.mlgi OUTPUT: grid.stor grid_outside_vor.area grid_outside.zone grid_multi_mat.zone grid_material.zone grid.fehmn grid_tet.gmv ======= Source Data ============================================================== Early Gridder version mesh with model descriptions: Report describing early version models with Gridder mesh descriptions: Harp, Dylan R., Philip H. Stauffer, Phoolendra K. Mishra, Daniel G. Levitt, and Bruce A. Robinson. "Thermal Modeling of High-Level Nuclear Waste Disposal in a Salt Repository." Nuclear Technology 18 7, no. 3 (2014): 294-307. From: Phillip Stauffer April 10 2012 New mesh ideas. . . since the spacing is 1 ft along the drift, I think we should go to 0.3 m spacing along the drift (0.984 ft close enough eh!). Our initial runs will be with 1 ft between canisters (1 node) then a canister which will be 2 nodes wide and 1 or 2 nodes high depending on how we do the vertical resolution. We may be able to do 0.3 vertical so that we hit the 4ft and 6ft of salt cover pretty close. Also, the end caps of 50ft are pretty loose now, so we may want to pull those in to more like 20 ft on each end for backfill, which is conservative. P