Grid Project:  ASCEM F-Basin

Grid Team Member: Terry Miller
Version: 3D Prism Mesh
Data Source: Haruko Wainwright and Jeff Thibault and Greg Flach
Delivered To: David Moulton and Haruko Wainwright
Date: May 22 2014 to present

Version History

Version Scope Finished Total
V3.1 V3 New Feature Data
from Haruko
Jul 23 2015 982,998 Nodes 1,849,039 Cells
V2.2 V2 with new Surfaces
from Haruko
May 19 2015 199,494 Nodes 371,518 Cells
V2.2b V2.2 converted to tet
for Rao
May 19 2015 199,494 Nodes 1,114,554 Tets
V2.1 3 layers, no pinchout
SRS Elevations
Nov 12 2014 199,494 Nodes 371,518 Cells
V1.0 3 layers, pinchout, 3 barriers May 22 2014 228,713 Nodes 428,930 Cells
V1.60 V1 with no pinchout (6 meters) Aug 25 2014 232743 Nodes 437080 Cells
V1.61 V1.60 with barrier cells removed Aug 26 2014 227153 Nodes 414600 Cells
V1.30 V1 with no pinchout (3 meters) Aug 25 2014 232743 Nodes 437080 Cells
V1.20 V1 with no pinchout (2 meters) Aug 25 2014 232743 Nodes 437080 Cells

Mesh Overview

This project is for ASCEM F-Basin Modeling. This mesh starts a triangluation based on the F-Basin shape and features; F-Basin and Barriers.The triangulated surface is used to create layer surfaces by assigning elevations equal to or relative to the DEM elevations. The surface layers are stacked and connected such that trianglular prism elements are formed between layers. Material elements (blocks) are identified to represent layers and features such as the basin and barriers. Exodus side sets are identified as boundary faces along bottom, sides, and top.

Future work may use a different boundary polygon and consider strategies for upper aquifer refinement so that elements at outlet are not so thin compared to element at outlet where layer is much thicker.

Mesh Files are available on local servers or you can download the compressed ".tar.gz" archive.
See mesh details for each version below.

Image pages for V3:

Image pages for V2 and V1:


V3 Prism Meshes

Source Data (SRS feet):
/scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/data/data_polys_v3
  barrier 1=A, 2=B, 3=C, 4=D, 5=E
  basins 1=West, 2=Middle, 3=fbasin box
  New Data, see 20150704domain_fig.tiff
  Archive:  domain_spec20150704.zip replaced with corrections

  Surfaces are from V2
  V2.2 SRS elevation data updates by Haruko:
  /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/data_surfs/Mesh_interfaces20150415
  Provided by Haruko Wainwright hmwainwright@lbl.gov
  Mesh_interfaces20150415.zip

  Boundary polygon from version 1 and 2 mesh work.





The following definitions V3 meshes with added feature data:

Element type: triangular prism formed from stacked triangulated surfaces

Materials (Exodus Blocks):

16  E barrier top layer one cell thick
15  D barrier top layer one cell thick
14  C barrier top layer one cell thick
13  B barrier top layer one cell thick
12  A barrier top layer one cell thick

11  E barrier 5 (east side)
10  D barrier 4 (intersect)
9   C barrier 3
8   B barrier 2
7   A barrier 1 (west)

6  box1 one cell thick
5  box2 one cell thick
4  box3 (old basin) one cell thick
3 Upper Layer - Upper Aquifer (to ground surface)
2 Middle Layer - Tan Clay Layer
1 Bottom Layer - Lower Aquifer


Boundary Faces (Exodus Side Sets):
 1 bottom 
 2 east 
 3 north inlet 
 4 west side 
 5 river outlet south face 
 6 top 
 7 river at top layer  
 8 top basin  3 (old fbasin) 
 9 top basin  2 (middle basin) 
10 top basin  1 (east basin) 
11 top barrier A  
12 top barrier B  
13 top barrier C  
14 top barrier D  
15 top barrier E  

Mesh information can be found in mesh directory in README.txt, RUNS, *.out
Mesh Files are written in formats; Exodus II (.exo), AVS (.inp), and GMV (.gmv)


V3.1 Prism Mesh with New box and barrier data features

number of nodes = 982998
number of elements = 1849039

Horizontal edge lengths from 12 m to 2 m near features.
Box features include Basin plus 2 more to the East
Barrier features are A, B, C, D, E (west to east)
Barriers are 2 elements wide, edge lengths of .61 m (2 ft)
No pinchout, upper aquifer has 1.5m thickness at outlet.
Upper Aquifer has 6 layers of proportional refinement.

V3.1 Exodus Mesh Overview

V3.1 All Mesh and Data images


MESH FILES for V3.1 Prism Mesh
  Jul 23 2015
   81M Jul 23 14:33 prism_m_ss.exo (Exodus with facesets)
  163M Jul 23 14:33 prism_m_ss.gmv (GMV binary)
   39K Jul 23 08:13 mesh_info.txt
    4K Jul 23 08:15 README.txt
    2M images directory

   Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v03.1/exodus 

   Download Archive:  fb_prismv3.1.tar.gz  (116M)  


V2 Prism Meshes


The following definitions apply to all the V1 and V2 meshes:

Element type: triangular prism formed from stacked triangulated surfaces

Materials (Blocks) for all grids:
1 Bottom Layer - Lower Aquifer
2 Middle Layer - Tan Clay Layer
3 Upper Layer - Upper Aquifer (to ground surface)
4 F-Basin (one element deep at top of Upper layer)
5 Barrier 1 (2 elements wide, depth is full Upper layer)
6 Barrier 2 (2 elements wide, depth is full Upper layer)l
7 Barrier 3 (2 elements wide, depth is full Upper layer)

Boundary Side Sets (Facesets) for all grids:
faceset ss1 = bottom
faceset ss2 = east side
faceset ss3 = north inlet
faceset ss4 = west side
faceset ss5 = south side river outlet
faceset ss6 = top (minus step and features)
faceset ss7 = top step feature
faceset ss8 = top of F-Basin feature
faceset ss9 = top of Barrier 1 feature
faceset ss10 = top of Barrier 2 feature
faceset ss11 = top of Barrier 3 feature

Mesh information can be found in mesh directory in README.txt, RUNS, *.out
Mesh Files are written in formats; Exodus II (.exo), AVS (.inp), and GMV (.gmv)


Source Data:

Barrier data:
  /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/data
  Provided by Haruko Wainwright hmwainwright@lbl.gov

V2.2 SRS elevation data updates by Haruko:
  /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/data_surfs/Mesh_interfaces20150415
  Provided by Haruko Wainwright hmwainwright@lbl.gov
  Mesh_interfaces20150415.zip
  V2.2 New_Interfaces_20150415.pdf 

V2 SRS elevation data:
  /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/data_surfs/haruko_email_data
  Provided by Haruko Wainwright hmwainwright@lbl.gov
  Uses data from ASCEM_transmittal.zip to convert to SRS coordinates.

2010 Elevation data (not in SRS coordinates):
  SRNL-L6200-2010-00025_plusData.zip
  /scratch/nts/gable/grid_gen/ascem/f_area_demo/data_download/Flow/ASCEM_transmittal
       /ASCEM_transmittal/Hydrostrat/cubit_quads
       /ASCEM_transmittal/reformat_polygon
       /ASCEM_transmittal/Topo/

Source elevation data:
Report SRNL-L6200-2010-00025.pdf
Greg Flach
Savannah River National Laboratory
773-42A, Savannah River Site, Aiken, SC 29808
803-725-5195
gregory.flach@srnl.doe.gov



V2 Prism Mesh with Basin and 3 Barriers Using SRS Elevations

V2.2 Input Elevations are from Haruko. Sent by email May 5 2015 in file Mesh_interfaces20150415.zip. Images showing comparisons of V2.2 and previous V2.1 surfaces are included in the Mesh and Surface QA Page.

These surface elevations do not pinchout the upper aquifer at outlet side of mesh.

V2.1 Input Elevations are from Haruko. They come from Greg's 2010 data set but are subset and converted to SRS coordinates. This conversion results in ground surface elevations that do not pinchout the upper aquifer at outlet side of mesh.






number of nodes = 199494
number of elements = 371518

Horizontal edge lengths from 50 m to 2 m near features.
F-Basin edge lengths near 7 m
Barriers are 2 elements wide, edge lengths of .61 m (2 ft)
No pinchout, upper aquifer has 1.5m thickness at outlet.
Upper Aquifer has 6 layers of proportional refinement.


MESH FILES for V2.2 Prism Mesh with Elevations Haruko Mesh_interfaces20150415
   May 12 2015
   README.txt mesh_info.out.txt
   prism_m_v2.20.exo prism_m_v2.20.gmv
   images/ 03_Gordon_elev.png 04_LA_elev.png 05_TanClay_elev.png 06_Top_elev.png prism_5x.png

   Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v02/3D_prism_v02.2/exodus

   Download Archive:  fb_mesh_prismv2.20.tar.gz  (24M)  

MESH FILES for V2.2b Prism Mesh Converted to Tet Mesh (for Rao)

This converts the V2.2 Prism mesh to a tet mesh. There are two versions, mat3 has 3 material layers as defined by the SRS Elevation surfaces. The mesh mat7 adds the 4 materials defining the F Basin box, and each of the 3 barriers.
number of nodes = 199494
number of tets = 1114554
Images of the tet mesh materials: V2.2 Tet Mesh Images

  May 19 2015
  prism_to_tet.out.txt
  tet_mat3.inp tet_mat3_ascii.gmv tet_mat3.gmv
  tet_mat7.inp tet_mat7_ascii.gmv tet_mat7.gmv
  images/ mesh_3mat_5x.png mesh_3mat_clr1_5x.png mesh_3mat_clr2_5x.png 
          mesh_3mat_clr3_5x.png mesh_3mat.png
          mesh_7mat_5x.png mesh_7mat_clr3_5x.png mesh_7mat_clrs124567_5x.png 
          mesh_7mat_clrs12567_5x.png mesh_7mat_clrs4567_5x.png mesh_7mat.png

   Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_tet_v02.2/output_files

   Download Archive:  fbasin_tet.tar.gz  (109M)  

MESH FILES for V2.1 Prism Mesh
  Nov 20 2014
  17M Nov 12 14:33 prism_m_ss.exo (Exodus with facesets)
  34M Nov 12 14:33 prism_m_ss.gmv (GMV binary)
  28K Nov 20 08:13 mesh_info.out.txt
  10K Nov 20 08:15 README.txt

   Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v02/mesh_top_ref6

   Download Archive (4 files):  fb_mesh_prismv2.1.tar.gz  (23M)  


V1 Prism Meshes


V1 Prism Mesh with Basin and 3 Barriers with Pinchout top layer

V3 Prism




number of elements = 428930
number of nodes = 228713

Horizontal edge lengths from 50 m to 2 m near features.
F-Basin edge lengths near 7 m
Barriers are 2 elements wide, edge lengths of .61 m (2 ft)


MESH FILES for V1.1 Prism Mesh
    May 28 2014
    19M May 22 14:49 prism_rmedge_m_ss.exo (Exodus with facesets)
    35M May 28 16:42 prism_stack_rmedge_m.gmv (GMV binary)
    87M May 21 10:52 prism_stack_rmedge_m.inp (AVS ascii format)
   6.6K May 28 12:46 README.txt

   Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v01/mesh_full

   Download Archive (4 files):  fb_mesh_prismv1.1.tar.gz  (38M)   



V1 Prism Mesh with Basin and 3 Barriers with No Pinchout in top layer

V3 Prism




number of elements = 437080
number of nodes = 232743

number of elements (barriers removed) = 414600
number of nodes (barriers removed) = 227153

Top surface (DEM) does not pinchout lower layer.
Top surface raised where layer is less than MIN_THICK

Horizontal edge lengths from 50 m to 2 m near features.
F-Basin edge lengths near 7 m
Barriers are 2 elements wide, edge lengths of .61 m (2 ft)

All Images for V1 Mesh with No Pinchout

Images V1.6 Mesh 6 meter thick

Images V1.3 Mesh 3 meter thick

Images V1.2 Mesh 2 meter thick



MESH FILES for V1.60 Prism Mesh plus 6 meters
 Aug 25 2014
 20M Aug 25 09:35 prism_m_ss.exo
 39M Aug 25 09:35 prism_m_ss.gmv
 25K Aug 25 07:12 mesh_info.out.txt
 11K Aug 25 11:21 README.txt

 Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v01_p6m/mesh_full

 Download Archive (4 files):  fb_mesh_prismv1.60.tar.gz  (27M)  

MESH FILES for V1.61 Prism Mesh plus 6 meters (Barrier cells removed)

Aug 26 2014
 19M Aug 26 10:58 prism_m_ss_rmB.exo
 37M Aug 26 10:58 prism_m_ss_rmB.gmv
 21K Aug 26 10:37 mesh_info.out.txt
 11K Aug 26 10:55 README.txt
172K Aug 26 10:57 prism_v01_p6m_rmB_mat1-4_2x.png
241K Aug 26 10:30 prism_v01_p6m_rmB_ss6_8_10.png
273K Aug 26 10:30 prism_v01_p6m_rmB_ss8_10.png
544K Aug 26 10:30 prism_v01_p6m_rmB_ss9.png
288K Aug 26 10:30 prism_v01_p6m_rmB_ss_sides_3x.png
2.3M Aug 26 10:36 prism_v01_p6m_rmB_ssBars.png

Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v01_p6m/mesh_rmBarriers

 Download Archive (10 files):  fb_mesh_prismv1.61.tar.gz  (29M)  

MESH FILES for V1.30 Prism Mesh plus 3 meters

Aug 25 2014
 20M Aug 25 11:30 prism_m_ss.exo
 39M Aug 25 11:30 prism_m_ss.gmv
 25K Aug 25 11:29 mesh_info.out.txt
 11K Aug 25 11:29 README.txt
168K Aug 25 12:49 prism_v01_p3m_detail_2x.png
327K Aug 25 12:49 prism_v01_p3m_E_2x.png
302K Aug 25 12:49 prism_v01_p3m_matall_2x.png
 32K Aug 25 12:49 prism_v01_p3m_ss_sides_3x.png
309K Aug 25 12:49 prism_v01_p3m_SW_2x.png

 Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v03/3D_prism_v01_p3m/mesh_full

 Download Archive (9 files):  fb_mesh_prismv1.30.tar.gz  (28M)  


MESH FILES for V1.20 Prism Mesh plus 2 meters
Aug 25 2014
 20M Aug 25 12:04 prism_m_ss.exo
 39M Aug 25 12:04 prism_m_ss.gmv
 25K Aug 25 12:06 mesh_info.out.txt
 11K Aug 25 12:09 README.txt
129K Aug 25 12:46 prism_v01_p2m_detail_2x.png
281K Aug 25 12:46 prism_v01_p2m_E_2x.png
316K Aug 25 12:46 prism_v01_p2m_matall_2x.png
287K Aug 25 12:46 prism_v01_p2m_SW_2x.png

 Local dir: /scratch/nts/tam/ascem/f_area_demo/pri_unstruc_v02/3D_prism_v01_p2m/mesh_full

 Download Archive (8 files):  fb_mesh_prismv1.20.tar.gz  (28M)  





Grid Team Contacts

Carl Gable <gable -at- lanl -dot- gov> 505-665-3533
Terry Miller <tamiller -at- lanl -dot- gov> 505-667-8009
Lucia Short <lucia -at- lanl -dot- gov> 505-665-7957