Grid Project: SALT SDDI Heaters

Grid Team Member:	Carl Gable, Dylan Harp, Terry Miller
Version:	3D Prism Mesh
Delivered To:	Phil Stauffer
Date:	March 2013 to present

Version History

Version	Scope	Finished	Total
V2	Half with 5 Heaters	Aug 1 2016	258,800 Nodes	1,476,328 Tets
V1	Quarter with 2.5 Heaters	Jun 1 2013	117,670 Nodes	667,192 Tets
V0	Quarter first version	March 2013	24,496 Nodes	133,476 Tets

Purpose of Project

Build a 3D mesh representing heater canisters in a salt repository for FEHM simulations. Because of symmetry of the room and panels, a quarter of the model is represented in V1 of the mesh. This is extended along the X-axis for V2 half model representation.

Quarter Mesh Geometry with 2.5 Heaters

Half Mesh Geometry with 5 Heaters

Mesh V2 Half Mesh with 5 Heaters

Use the same LaGriT command files that built V1 Quarter Mesh, but update geometry to extend the X-axis into half model representation of heaters in salt repository.

Images of V1 and V2 Mesh

README_v2.txt for V2 Half Mesh


Number Nodes:    258800
Number Elements: 1476328

Mesh Material Zones for FEHM (Heaters left to right):

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

V2 FEHM Mesh Files:

Mesh Directory: /scratch/sft/tam/SALT/SDDI_meshV2_half

-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 V1 Quarter Mesh with 2.5 Heaters

Refined tet Mesh based on modifications from an early gridder mesh to a refined tet mesh using LaGriT.

Images of V2 that include V1 images


Number Nodes:    258800
Number Elements: 1476328

Mesh Material Zones for FEHM (Heaters left to right):

Material Heater 1  1 has       225 nodes. #nodes/nnodes is   0.191212713253E-02 
Material Heater 2  2 has       396 nodes. #nodes/nnodes is   0.336534366943E-02 
Material Heater 3  3 has       396 nodes. #nodes/nnodes is   0.336534366943E-02 
Material Salt      4 has     32770 nodes. #nodes/nnodes is   0.278490692377     
Material Air       5 has     41271 nodes. #nodes/nnodes is   0.350735098124     
Material Room      6 has     36554 nodes. #nodes/nnodes is   0.310648411512     
Material Outside   7 has      6058 nodes. #nodes/nnodes is   0.514829605818E-01

V1 Quarter Mesh with 2.5 Heaters

Mesh Directory: /scratch/er/dharp/projects/salt/SDDI/mesh3d13

-rw-r--r-- 1 dharp sft 35709496 Jun 11  2013 grid.fehmn
-rw-r--r-- 1 dharp sft  1294613 Jun 11  2013 grid_material.zone
-rw-r--r-- 1 dharp sft  2556434 Jun 11  2013 grid_multi_mat.zone
-rw-r--r-- 1 dharp sft    78632 Jun 11  2013 grid_outside.zone
-rw-r--r-- 1 dharp sft   453092 Jun 11  2013 grid_outside_vor.area
-rw-r--r-- 1 dharp sft 25521138 Jun 11  2013 grid.stor
-rw-r--r-- 1 dharp sft 27130148 Jun 11  2013 grid_tet01.gmv

Mesh V0 Quarter Mesh with 2.5 Heaters

Development of a refined tet mesh based on modifications from an early gridder mesh to a refined tet mesh using LaGriT.


Number Nodes:     24496
Number Elements: 133476

Mesh Material Zones for FEHM:

Material  1 has       540 nodes. #nodes/nnodes is   0.220444146544E-01 
Material  2 has       940 nodes. #nodes/nnodes is   0.383736118674E-01 
Material  3 has       940 nodes. #nodes/nnodes is   0.383736118674E-01 
Material  4 has      7357 nodes. #nodes/nnodes is   0.300334751606     
Material  5 has      2374 nodes. #nodes/nnodes is   0.969137847424E-01 
Material  6 has     12345 nodes. #nodes/nnodes is   0.503959834576

V0 Quarter Mesh

Mesh Directory: /home/gable/people/harp/mesh3d02/ 

-rw-r--r-- 1 gable sft   7211786 Mar 27  2013 grid.fehmn
-rw-r--r-- 1 gable sft    269666 Mar 27  2013 grid_material.zone
-rw-r--r-- 1 gable sft    515348 Mar 27  2013 grid_multi_mat.zone
-rw-r--r-- 1 gable sft     40880 Mar 27  2013 grid_outside.zone
-rw-r--r-- 1 gable sft    235160 Mar 27  2013 grid_outside_vor.area
-rw-r--r-- 1 gable sft   5142791 Mar 27  2013 grid.stor

Data Sources


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 187, no. 3 (2014): 294-307.

Mesh Parameters for LaGriT command files:

   Y extent of EDZ into wall 3.9624
   Z extent of EDZ below floor -1.542
   Z extent of EDZ above ceiling  4.572

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

Meshing Method:


Create the Level 0 base hex grid

  lagrit < grid_L0.lgi

edit geometry.mlgi or geometry_half.mlgi
test result with lagrit < test_geometry.lgi

  lagrit < refine_main.lgi

Aux LaGriT infiles:
  geometry_half.mlgi, geometry.mlgi
  build_eltsets.mlgi
  refine_ESET.mlgi

Grid Team Contacts

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