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	Full Stacked Mesh with Ice Wedge Cut to show internal materials
Mesh and Ice Wedge with moss and organic layers removed	Ice Wedge

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:

1 = bottom
2 = top
3 = East right
4 = North back
5 = West left
6 = South front
7 = Corner side
8 = Corner top

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

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

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

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