Grid Project:   Barrow Ice Wedge Polygon 40 V2 (fine scale)

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 2 is a modification of the Version 1 mesh which is a 3D mesh polygon based on single polygon #40 generated using LiDAR data set (area C). This mesh has changes to V1 based on the Mesh Input Parameters below. It has increased vertical resolution and a flat bottom. The convention of layer order was also changed so that numbering is from top to bottom. Two Exodus II mesh files are written, one multi-material with ice, and one single material.

Full Stacked Mesh with Ice Wedge	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 V2

General Description

Same as V01 using Polygon #40 with the following modifications:

• The live organic layer (uppermost layer) should be 2 cm thick with 2 cells.
• The dead organic layer should be 48 cm thick with 24 cells. (That would be materials 3 and 4 combined, so you could make each 24 cm think).
• The inorganic layer (layer 3 from the top) should have 2 cm thick cells near its top but then coarser resolution below - about 10 cm at the bottom.
• Materials 1 and 2 are both mineral soils - I need higher resolution at the top of the mineral soil.
• The domain bottom should be flat at Z = -1.

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 (exists in material 4 and 5)

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

Output Mesh Description

Output is a stacked prism mesh written as two Exodus files; a single material mesh and a multi-material Mesh with 6 materials (5 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

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     


SUMMARY VOLUME/AREA of ELEMENT COLORS (BLOCKS) 
Color      Num. Elements   Volume          Fractional Volume                    
    1       12544          0.7936425E+01     0.003490890                        
    2       75264          0.9523729E+02     0.041890767                        
    3       75264          0.9523710E+02     0.041890684                        
    4       95524          0.1464469E+03     0.064415680                        
    5      293616          0.1745817E+04     0.767909504                        
   30       49900          0.1827924E+03     0.080402475                        

SUMMARY ELEVATIONS of LAYER TOPS 
LAYER     MIN              MAX             DIFFERENCE 
  1     4.504300117E+00  4.992260933E+00 4.879608154E-01
  2     4.484300117E+00  4.972260933E+00 4.879608154E-01
  3     4.244300117E+00  4.732260933E+00 4.879608154E-01
  4     4.004300117E+00  4.492260933E+00 4.879608154E-01
  5     3.504300117E+00  3.992260933E+00 4.879608154E-01
  (Bottom is flat at -1.0)

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

* Materials and Cell Layers *
ID Name Cell Height Num Cells Mat Thick 1 moss (top layer) 0.01 2 0.20 2 upper organic 0.02 12 0.24 3 lower organic 0.02 12 0.24 4 upper mineral 0.025 20 0.5 5 lower mineral 0.09 50 2.2 - 4.8 30 ice 0.09 variable 0.00 - 2.3

* Boundary Side Sets *
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_v2 

tar file sent:
-rw-r--r-- 1 lucia sftus 164147118 Jun  3 11:33 mesh_3D_poly40_v02.tgz

Created from clean subdirectory with mesh files:
drwxr-xr-x 3 lucia sftus     12288 Jun  3 11:27 mesh_final/

LaGriT faceset files from mesh outside surface (elem number and face):
-rw-r--r--  1 lucia sftus     81639 Jun  3 11:12 fs_1_bot.faceset
-rw-r--r--  1 lucia sftus     69095 Jun  3 11:12 fs_2_top.faceset
-rw-r--r--  1 lucia sftus     86119 Jun  3 11:12 fs_3_rig.faceset
-rw-r--r--  1 lucia sftus     43111 Jun  3 11:12 fs_4_bac.faceset
-rw-r--r--  1 lucia sftus     21607 Jun  3 11:12 fs_5_lef.faceset
-rw-r--r--  1 lucia sftus     43111 Jun  3 11:12 fs_6_fro.faceset

Mesh and faceset files for 6 Materials:
-rw-r--r--  1 lucia sftus  26953248 Jun  3 11:12 mesh_fsets.exo
-rw-r--r--  1 lucia sftus  45822656 Jun  3 11:12 mesh_fsets.gmv
-rw-r--r--  1 lucia sftus  89405021 Jun  3 11:12 mesh_fsets.inp

Mesh and faceset files for 1 Material:
-rw-r--r--  1 lucia sftus  26951604 Jun  3 11:12 mesh_fsets_1mat.exo
-rw-r--r--  1 lucia sftus  42678672 Jun  3 11:12 mesh_fsets_1mat.gmv
-rw-r--r--  1 lucia sftus  89093748 Jun  3 11:12 mesh_fsets_1mat.inp

LaGriT output files contain information reported during the run.
-rw-r--r-- 1 lucia sftus  20338 Jun  3 11:12 interpolate_icewedge.out
-rw-r--r-- 1 lucia sftus 110772 Jun  3 11:12 make_poly.out
-rw-r--r-- 1 lucia sftus 261870 Jun  3 11:12 stack_full_poly.out
-rw-r--r-- 1 lucia sftus  77436 Jun  3 11:12 write_exo.out

Mesh images:
drwxr-x--x  2 lucia sftus      4096 Jun  3 11:18 images/
-rw-r--r--  1 lucia sftus     20338 Jun  3 11:12 interpolate_icewedge.out

Input Data files:

 /scratch/sft2/lucia/grid_gen/arctic/arctic_fine_v2/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_v2

Links to Images and Project Details

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