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

Grid Team Member: Quan Bui, Carl Gable
Grid Versions: XY and rotate XZ V1_extrude_1m
Delivered To: Shaoping Chu, Scott Painter
Date Completed: 04/14/2013
Rotated versions (Shaoping and Terry): 05/29/2013

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.

Material Zones
5 - Red Background
4 - Yellow Deformation Zone
3 - Dark Green Repository
2 - Cyan Top Layer 2
1 - Blue Top Layer 1


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.

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 ->  

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                            


Exodus Materials and Facesets
LaGriT Finite Element Mesh
Exodus Materials
LaGriT Finite Element Mesh
LaGriT Finite Element Mesh
LaGriT Finite Element Mesh
LaGriT Finite Element Mesh
LaGriT Finite Element Mesh
4 Nodes with Neg Ccoefs
LaGriT Finite Element Mesh
tet_hole (version not used)
LaGriT Finite Element Mesh

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.

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
-rw-r--r-- 1 tam sft     229846 May 29 14:28
-rw-r--r-- 1 tam sft     405365 May 29 14:28
-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 :


-rw-r----- 1 quanb altern 2461978 Jan 25 10:57 mesh.fehmn
-rw-r----- 1 quanb altern  198793 Jan 25 10:57
-rw-r----- 1 quanb altern   76822 Jan 25 10:57
-rw-r----- 1 quanb altern   23176 Jan 25 10:57 mesh_outside_vor.area
-rw-r----- 1 quanb altern    4116 Jan 25 10:57
-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

Carl Gable gable -at- lanl -dot- gov 505-665-3533
Quan Bui quanb -at- lanl -dot- gov 505-667-7362