Grid Project:   3D 1m Generic Repository (Extruded 2D)

Project Description

This is a summary of work that started with the Generic Repository 2D in XY. The 3D mesh is rotated so Z points up. Exodus files and facesets are written. Note the rotation may experience precision issues. A version that avoids precision issues - switches axis and connects the points by Delaunay did not work because connect left a triangle shaped hole in middle of mesh.

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

• For the repository: 1x1m resolution; repository size is 200m x 10m
• For the deformation zone, use 2x2m resolution
• For the rest of the area: 10m x 10m
• The grid resolution is higher near the repository and then gradually coarsens to the surrounding area of 10 m x 10m grid. The max Y axis is the top of the grid and is flat at 600, has a conforming slope from X=300 to X=700, then is flat at 550.

Project Details

Date: May 13, 2013 3:04:33 PM MDT
Subject: 3D Repo 1m with Exodus mesh and facesets

Read the tet_exodus.exo file with gmv, 
then select surfaces under Display menu to toggle the facesets on or off.
The first image are the materials and the top facesets.
The second image is just the materials, no facesets.

/scratch/sft/tam/grid_gen/generic_repository/
First version starting with Shaoping grids: 3Dgrid1m_v1_exodus
Final version starting with Quan grids: 3Dgrid1m_quan_rot

Quan's original 3D 1m XY mesh meshxy_quan.gmv ->  
 /home/quanb/src/grid_gen/generic_repository2D/v1_extrude_1m/mesh.gmv 

number of dimension: 3 
number of nodes:     36112 
number of elements: 107268 

This mesh (rotated zone) counts:
Material           1 has       132 nodes. #nodes/nnodes is   0.365529465489E-02 
Material           2 has        66 nodes. #nodes/nnodes is   0.182764732745E-02 
Material           3 has      4848 nodes. #nodes/nnodes is   0.134249001741     
Material           4 has      7890 nodes. #nodes/nnodes is   0.218486934900     
Material           5 has     23176 nodes. #nodes/nnodes is   0.641781151295  

Both rotated and un-rotated meshes have 4 neg ccoefs.
This mesh:

*** Construct and Compress Sparse Matrix:3D ***                                 
   *** Compress Area Coefficient Values ***                                     
AMatbld3d_stor: *****Negative Coefficients ******                               
AMatbld3d_stor: Total Number of Negative Coefficients         4                 
AMatbld3d_stor: Number of Significant Negative Coefs         4                  
AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs         0                 
Negative coef at row     1509  Column     2054 with value  -3.5372844E-06       
Negative coef at row     1512  Column     2055 with value  -8.3614028E-06       
Negative coef at row    19565  Column    20110 with value  -3.5372844E-06       
Negative coef at row    19568  Column    20111 with value  -8.3614028E-06       
AMatbld3d_stor: npoints =    36112  ncoefs =     272812                         
AMatbld3d_stor: Number of unique coefs =     48096                              
AMatbld3d_stor: Maximum num. connections to a node =         15                 
AMatbld3d_stor: Volume min =   3.1879425E-01                                    
AMatbld3d_stor: Volume max =   5.7019814E+01                                    
AMatbld3d_stor: Total Volume:   5.7500000E+05                                   
AMatbld3d_stor: abs(Aij/xij) min =   0.0000000E+00                              
AMatbld3d_stor: abs(Aij/xij) max =   1.1403963E+02                              
AMatbld3d_stor: (Aij/xij) max =   8.3614028E-06                                 
AMatbld3d_stor: (Aij/xij) min =  -1.1403963E+02                                 
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 tet.stor                            

Images

Exodus Materials and Facesets	Exodus Materials
tet	tet_zones_pts
tet_repo_pts	tet_zones
4 Nodes with Neg Ccoefs	tet_hole (version not used)

Project Archive Locations

FEHMN files for 3D XZ (rotated from XY):

Final Mesh Used with FEHM WINDOWS version.
The Linux version quits after about 20 time-steps due to numerical confusion.

/scratch/sft/tam/grid_gen/generic_repository/3Dgrid1m_quan_rot
uses meshxy_quan.gmv 
  -> /home/quanb/src/grid_gen/generic_repository2D/v1_extrude_1m/mesh.gmv

Final Files:

-rw-r--r-- 1 tam sft       2572 May 29 14:27 rotate_fehm.lgi
-rw-r--r-- 1 tam sft    7497267 May 29 14:28 meshxz_geom.inp
-rw-r--r-- 1 tam sft    6961994 May 29 14:28 tet.fehmn
-rw-r--r-- 1 tam sft     397409 May 29 14:28 tet_material.zone
-rw-r--r-- 1 tam sft     229846 May 29 14:28 tet_multi_mat.zone
-rw-r--r-- 1 tam sft     405365 May 29 14:28 tet_outside.zone
-rw-r--r-- 1 tam sft    2339231 May 29 14:28 tet_outside_vor.area
-rw-r--r-- 1 tam sft    8371597 May 29 14:28 tet.stor
-rw-r--r-- 1 tam sft   13605488 May 29 14:28 tet_final.inp
-rw-r--r-- 1 tam sft    5450872 May 29 14:28 tet_final.gmv
-rw-r--r-- 1 tam sft    3442564 May 29 14:28 tet_final.exo

Files with outside attributes for Dylan and Satish:
-rw-r--r-- 1 tam sft   20609199 May 29 16:33 tet_out_attrib.inp
-rw-r--r-- 1 tam sft    9495328 May 29 16:33 tet_out_attrib.gmv

Original FEHMN files for 2D XY :

/home/quanb/src/grid_gen/generic_repository2D

-rw-r----- 1 quanb altern 2461978 Jan 25 10:57 mesh.fehmn
-rw-r----- 1 quanb altern  198793 Jan 25 10:57 mesh_material.zone
-rw-r----- 1 quanb altern   76822 Jan 25 10:57 mesh_multi_mat.zone
-rw-r----- 1 quanb altern   23176 Jan 25 10:57 mesh_outside_vor.area
-rw-r----- 1 quanb altern    4116 Jan 25 10:57 mesh_outside.zone
-rw-r----- 1 quanb altern 6019974 Jan 25 10:57 mesh.stor


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