Grid Project:   Generic Repository 3D octree

Version History

Octree 25m joined with conforming Tet slope	Octree 50m joined with conforming Tet slope
Detail of Zones: Repository, Fault, and Top	Detail of Zones: Repository, Fault, and Top
Gallery Images 25m Mesh	Gallery Images 50m Mesh

Project Description

This is a summary of work that started with the Generic Repository in 2D XY. The boundaries and criteria remain the same except to accomodate spacing restrictions.

This version uses the design criteria same as original, materials stay the same.

• Background mesh is 1000 x 1000 meters in the x,y. The top is Z upward with max elevation at 600 meters. The top is truncated by a slope that bottoms at 550 meters.
• repository 200m x 200m x 20m with finest resolution
• fault, deformation zone, is 20m wide along full mesh with finest resolution
• Background spacing is 25m and 50m (coarse).
• Octree Refinement is 3 levels to 3.125m and 6.35m (coarse).
• The top slope between X=300m and X=700m is maintained by building a seperate top piece that is connected with Delaunay triangulation and joined to the octree hex mesh.
• FHEM zone definitions are. 1 (Top left), 2 (Top right), 3 (Repository), 4 (Fault), 5 (Background)
• Exodus II file definitions are: 1 (Background), 2 (Fault), 3 (Repository), SS1 (Top left), SS2 (Top Right)

Project Details

Method Overview:
Create the hex mesh for the 1000x1000x600 box with the start spacing of 25m or 50m. Cut out the area where the slope piece will go then refine the fault and repository zones. Use the edges cut from the background grid to form the left and bottom sides of the slope piece. Triangulate these edges into a Delaunay triangulation, then stack into 3D. Join the 3D slope and the octree mesh and connect Delaunay into the final tet mesh.
The triangulated slope piece must be Delaunay (recon 1).
Start refinement of fault wider than the fault, otherwise the fault will stair step outside the smallest refine area.
Repository box is centered at Z=300, select height based on resolution.

-------

The 25m octree mesh:
number of nodes    =     634148        
number of elements =    3686223    

FEHM Material Zones
Material Top left     1 has      1722 nodes. #nodes/nnodes is   0.271545443684E-02
Material Top right    2 has      1296 nodes. #nodes/nnodes is   0.204368703999E-02
Material Repository   3 has     29575 nodes. #nodes/nnodes is   0.466373786330E-01
Material Fault        4 has    478130 nodes. #nodes/nnodes is   0.753972232342
Material Background   5 has    123425 nodes. #nodes/nnodes is   0.194631218910

Exodus Materials and Facesets:
Material 1 background num_el_in_blk1 = 1081173 ;
Material 2 repository num_el_in_blk2 = 2458156 ;
Material 3 fault      num_el_in_blk3 = 146894 ;
Faceset  1 Top left   num_side_ss1 = 3360 ;
Faceset  2 Top right  num_side_ss2 = 8240 ;

-------
The 50m octree mesh:
number of nodes    =  160583
number of elements =  921085

FEHM Material Zones
Material Top left     1 has     1417 nodes. #nodes/nnodes is   0.882409699261E-02
Material Top Right    2 has     1173 nodes. #nodes/nnodes is   0.730463396758E-02
Material Repository   3 has     5445 nodes. #nodes/nnodes is   0.339077003300E-01
Material Fault        4 has    46207 nodes. #nodes/nnodes is   0.287745267153
Material Background   5 has   106341 nodes. #nodes/nnodes is   0.662218272686

Exodus Materials and Facesets:
Material 1 background num_el_in_blk1 = 531503 ;
Material 2 repository num_el_in_blk2 = 365103 ;
Material 3 fault      num_el_in_blk3 = 24479 ;
Faceset  1 Top left   num_side_ss1 = 2780 ;
Faceset  2 Top right  num_side_ss2 = 3540 ;

Project Archive Locations

FEHMN files for 3D Octree 25m

/n/www_meshing_local/meshing/proj/UFD_GenericRepo_3D_octree

FEHM Files:
-rw-r--r-- 1 tam sft 153471776 Jun 19 14:30 tet_final_zones.gmv
-rw-r--r-- 1 tam sft 369245288 Jun 19 14:30 tet_final_zones.inp
-rw-r--r-- 1 tam sft 196159705 Jun 19 14:31 tet.fehmn
-rw-r--r-- 1 tam sft 134918966 Jun 19 14:33 tet.stor
-rw-r--r-- 1 tam sft   6975805 Jun 19 14:31 tet_material.zone
-rw-r--r-- 1 tam sft  10760950 Jun 19 14:32 tet_multi_mat.zone
-rw-r--r-- 1 tam sft    205198 Jun 19 14:32 tet_outside.zone
-rw-r--r-- 1 tam sft   1183723 Jun 19 14:32 tet_outside_vor.area

Exodus Files:
-rw-r--r-- 1 tam sft  47143 Jun 19 15:47 exo_1_top_left.faceset
-rw-r--r-- 1 tam sft 123703 Jun 19 15:47 exo_2_top_right.faceset
-rw-r--r-- 1 tam sft 103784288 Jun 19 15:47 tet_exodus.exo
-rw-r--r-- 1 tam sft 156752665 Jun 19 15:48 tet_exodus_ascii.exo

FEHMN files for 3D Octree 50m (coarse)

  /scratch/sft/tam/grid_gen/generic_repository/3Dgrid_octree_coarse

FEHM Files:
-rw-r--r-- 1 tam sftus 49165932 Jun 19 16:47 tet.fehmn
-rw-r--r-- 1 tam sftus 34312186 Jun 19 16:47 tet.stor
-rw-r--r-- 1 tam sftus  1766590 Jun 19 16:47 tet_material.zone
-rw-r--r-- 1 tam sftus  4705462 Jun 19 16:47 tet_multi_mat.zone
-rw-r--r-- 1 tam sftus    94923 Jun 19 16:47 tet_outside.zone
-rw-r--r-- 1 tam sftus   547136 Jun 19 16:47 tet_outside_vor.area
-rw-r--r-- 1 tam sftus 38467720 Jun 19 16:47 tet_final_zones.gmv
-rw-r--r-- 1 tam sftus 91938827 Jun 19 16:47 tet_final_zones.inp

Exodus Files:
-rw-r--r-- 1 tam sftus    36243 Jun 19 16:49 exo_1_top_left.faceset
-rw-r--r-- 1 tam sftus    49663 Jun 19 16:49 exo_2_top_right.faceset
-rw-r--r-- 1 tam sftus 26013176 Jun 19 16:49 tet_exodus.exo
-rw-r--r-- 1 tam sftus 36906997 Jun 19 16:49 tet_exodus_ascii.exo

Shaoping's work /scratch/sft/spchu/reference_case/
Dylan's work /home/dharp/projects/IntToolDev/2013/mesh/pflotran
Terry's work on debugging /scratch/sft/tam/grid_gen/generic_repository 

Grid Team Contacts