Grid Project:   Barrow Ice Wedge Polygon 40 V3.0.1 (fine scale)

Grid Team Member: Lucia Short
Grid Version: Version 03.0.1:
single polygon #40 with corners
Data Source: LiDAR from Area C
Delivered To: Scott Painter
Date Completed: July 19, 2013
V3.0.1 Mesh (6 layers & corners): 311,273 Nodes and 602,112 Cells
V3 Mesh (6 layers) : 311,273 Nodes and 602,112 Cells
V2 Mesh (5 layers) : 311,273 Nodes and 602,112 Cells
V1 Mesh (prototype):   51,344 Nodes and  94,080 Cells
   

Purpose of Project

These meshes were developed in support of the Arctic Next Generation Ecosystem Experiment (NGEE-Arctic) Project and LANL Laboratory Directed project LDRD20120068DR. These two projects share a need for high quality computational meshes for use in fine-scale modeling of thermal hydrology in ice-wedge polygon landscapes.

Version 3.0.1 is a modification of the Version 3.0 mesh which adds corner faces to the exodus file.

This version has a correction to file write_exo.lgi and now has two additional layers for corner and top faces.

The following are additional faces that were added to the V3.0.0:

  -----------------
  layer 7 -> corner side faces
  layer 8 -> corner top faces
  -----------------

Full Stacked Mesh with Ice Wedge
6 material layers plus ice material

LaGriT Finite Element Mesh

Full Stacked Mesh with Ice Wedge
Cut to show internal materials

LaGriT Finite Element Mesh


Mesh and Ice Wedge
with moss and organic layers removed

LaGriT Finite Element Mesh


Ice Wedge

LaGriT Finite Element Mesh

Mesh Input Parameters V3.0.1

General Description

Same as V03 using Polygon #40.

Material Ice Wedge

We want the top of the ice wedge to be at 50 cm depth (corresponds to the top of the inorganic layer).

Top of wedge starts at top of mineral layer (replaces cells only in materials 5 and 6)

Top wedge width (distance from polygon) is aprox 1.2

Depth from top of ice to tip is aprox 2.8

Outside Boundaries (Exodus side sets)

The side sets defined are top, bottom, and sides (corresponding to the 4 sectors within 45 degrees of each direction). The exodus file contains following surface facesets:

Faces 7 and 8 are new to V3.0.1



Output Mesh Description

Output is a stacked prism mesh written as two Exodus files; a single material mesh and a multi-material Mesh with 7 materials (6 layers plus the Ice Wedge). Files are written in formats; Exodus II (.exo), AVS (.inp), and GMV (.gmv)
Final Mesh files: mesh_fsets.exo, mesh_fsets.gmv, mesh_fsets.inp
Final Mesh with single material: mesh_fsets_1mat.exo, mesh_fsets_1mat.gmv, mesh_fsets_1mat.inp

Mesh Element Type:             prism
Mesh Software Used: LaGriT, ExodusII

Total nodes:       311273       
Total elements:    602112     
Total Volume:      0.2273467E+04      

from file write_exo.out.txt:
Mesh Extents (Meters):
        MIN              MAX             DIFFERENCE    
 x    4.634709930E+01  7.183999634E+01 2.549289703E+01     
 y    1.062200012E+02  1.326065216E+02 2.638652039E+01     
 z   -1.000000000E+00  4.992260933E+00 5.992260933E+00     

from file write_exo.out.txt:
SUMMARY VOLUME/AREA METERS of ELEMENT COLORS (BLOCKS) 
Color    Num. Elements   Volume          Fractional Volume                    
    1       12544      0.7936425E+01     0.003490890                        
    2       56448      0.7142782E+02     0.031418013                        
    3       56448      0.7142801E+02     0.031418096                        
    4       37632      0.4761855E+02     0.020945342                        
    5       95524      0.1464469E+03     0.064415680                        
    6      293616      0.1745817E+04     0.767909504                        
    7/30   49900       0.1827924E+03     0.080402475        

from file stack_full_poly.out.txt:
SUMMARY ELEVATIONS at LAYER TOPS (METERS)
        MIN          MAX       DIFF MAX below (thickness)   
    1 4.504300117  4.992260933   .02
    2 4.484300117  4.972260933   .18
    3 4.304300117  4.792260933   .18
    4 4.124300117  4.612260933   .12
    5 4.004300117  4.492260933   .50
    6 3.504300117  3.992260933  4.99 
      (Bottom is flat at -1.0)

Mesh Layer Definitions (meters)

ID Layer Name Layer
Thick
Layer
Vertical Cells
Layer
Volume
Min Top
Elevation
Max Top
Elevation
1 moss(top) 0.02 2 x .01 7.936425 4.50 4.99
2 upper organic 0.18 9 x .02 79.36425 4.48 4.97
3 upper mineral 0.18 9 x .02 71.42801 4.30 4.79
4 lower organic 0.12 6 x .02 47.61855 4.12 4.61
5 lower mineral 0.50 20 x .025 146.4469 4.00 4.49
6 deep mineral 2.2 -
4.8
50 x .09 1745.817 3.50 3.99
7/30 ice 0.00 -
2.3
variable 182.7924 4.49 4.00
LaGriT Finite Element Mesh

LaGriT Finite Element Mesh

Boundary Side Sets

Corners:
purple side (ss_c_side) has 10 faces              black  top  (ss_c_top)  has 13 faces
Top and Bottom:
neon-pink  top    (num_side_ss2) has 6272 faces
neon-green bottom (num_side_ss1) has 6272 faces
 
Sides:
red    east right  (num_side_ss3) has 6144 faces
yellow north back  (num_side_ss4) has 3072 faces
green  west left   (num_side_ss5) has 1536 faces
blue   south front (num_side_ss6) has 3072 faces
LaGriT Finite Element 
Mesh LaGriT Finite Element 
Element Mesh


Project Archive Location

Full Set of files:

/scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v3_1_poly40/ 

tar file was not sent, but was created (*tgz):
-rw-r--r--   1 lucia  sftus  170910794 Aug 26 16:48 mesh_fine_poly40_v03_0_1.tgz

Created from clean subdirectory with mesh files:
drwxr-xr-x   3 lucia  sftus       8192 Aug 26 16:46 arctic_v3_0_1

LaGriT faceset files with elem number and face (*.faceset):
-rw-r--r--  1 lucia  sftus  81639 Jul 19 10:03 fs_1_bot.faceset
-rw-r--r-- 1 lucia sftus 69095 Jul 19 10:03 fs_2_top.faceset
-rw-r--r-- 1 lucia sftus 86119 Jul 19 10:03 fs_3_rig.faceset
-rw-r--r-- 1 lucia sftus 43111 Jul 19 10:03 fs_4_bac.faceset
-rw-r--r-- 1 lucia sftus 21607 Jul 19 10:03 fs_5_lef.faceset
-rw-r--r-- 1 lucia sftus 43111 Jul 19 10:03 fs_6_fro.faceset
-rw-r--r-- 1 lucia sftus 223 Jul 19 10:03 fs_7_Csides.faceset
-rw-r--r-- 1 lucia sftus 246 Jul 19 10:03 fs_8_Ctop.faceset Mesh and faceset files for 6 Materials (mesh_fsets.*):
-rw-r--r-- 1 lucia sftus 26954056 Jul 19 10:04 mesh_fsets.exo
-rw-r--r-- 1 lucia sftus 46374600 Jul 19 10:04 mesh_fsets.gmv
-rw-r--r-- 1 lucia sftus 89405021 Jul 19 10:04 mesh_fsets.inp Mesh and faceset files for 1 Material (mesh_fsets_1mat.*):
-rw-r--r-- 1 lucia sftus 26952084 Jul 19 10:04 mesh_fsets_1mat.exo
-rw-r--r-- 1 lucia sftus 42678672 Jul 19 10:04 mesh_fsets_1mat.gmv
-rw-r--r-- 1 lucia sftus 89093748 Jul 19 10:04 mesh_fsets_1mat.inp LaGriT output files contain information reported during the run (*.out.txt):
-rw-r--r-- 1 lucia sftus 4963 Jul 19 10:02 get_boundary_pts.out.txt
-rw-r--r-- 1 lucia sftus 4961 Jul 19 10:02 get_poly_surf.out.txt
-rw-r--r-- 1 lucia sftus 20338 Jul 19 10:03 interpolate_icewedge.out.txt
-rw-r--r-- 1 lucia sftus 110772 Jul 19 10:02 make_poly10040.out.txt
-rw-r--r-- 1 lucia sftus 5943 Jul 19 10:02 merge_ice.out.txt
-rw-r--r-- 1 lucia sftus 269338 Jul 19 10:03 stack_full_poly.out.txt
-rw-r--r-- 1 lucia sftus 105639 Jul 19 10:04 write_exo.out.txt Mesh images: /scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v3_1_poly40/images_v3_1_p40

Input Data files:

 /scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v3_1_poly40/input/

 -rw-r--r--  2 lucia sftus 30928912 Apr 18 13:59 lidar_site3.inp
 -rw-r--r--  1 lucia sftus      123 Apr 18 13:59 poly10001.avs
 -rw-r--r--  1 lucia sftus      178 Apr 18 13:59 poly10002.avs 
						...
 -rw-r--r--  1 lucia sftus      155 Apr 18 13:59 poly10071.avs
 -rw-r--r--  1 lucia sftus      190 Apr 18 13:59 poly10072.avs

Main Work Directory for Polygon Meshing Method:

/scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v3_1_poly40/

 

Links to Images and Project Details


Gallery of All Project Images

Method Description Overview (from V01)

Input and Output Files Running LaGriT to create the Mesh

These are LaGriT input (.lgi) and output (outx3dgen) files used in creating this mesh. The output files provide good information about the data and mesh, including min and max values of mesh attributes.
The complete set of scripts and input files are in /scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v3_1_poly40. The output files for various mesh versions are in /scratch/sft2/lucia/grid_gen/arctic/.

Step 1: Create Ice Object from Polygon (make_poly.lgi)     Input       Output

Step 2: Create and stack layers into full mesh (stack_full_poly.lgi)     Input       Output

Step 3: Interpolate Ice Material on to Full Stacked Mesh (interpolate_icewedge.lgi)     Input       Output

Step 4: Create facesets and write final Exodus files (write_exo.lgi)     Input       Output

 




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