Prism Mesh V2.1 Jan 21 2016 /scratch/tundra/tam/IDEAS/Jaramillo/EXO_full_v2.1 Modified from V2.0 from mesh_unstruc_stack.lgi:# define / L1 / .25 mesh_unstruc_stack.lgi:# define / L2 / .75 mesh_unstruc_stack.lgi:# define / L3 / 2.0 to mesh_unstruc_stack.lgi:define / L1 / .10 mesh_unstruc_stack.lgi:define / L2 / .15 mesh_unstruc_stack.lgi:define / L3 / 0.75 mesh_unstruc_stack.lgi:define / L4 / 2.0 This mesh is translated to zero coordinate: define XTRANS 360600. define YTRANS 3973000. (re-run Sep 17 to remove duplicate nodes) Uses Richard's 1 meter DEM interpolated on to surface for top. Additional layering formed by translating this top surface. Unstructured stacked prism Mesh River 1 meter wide (same elev across river elements) River buffer 1 meter wide The River node elevations are interpolated from a center line DEM This is an improved iteration triangulation V1 bndry and river polygons. The surface triangluation is constrained to a max edge length of 30m. This triangulation is Delaunay and has better element ratios and distribution. Otherwise the 3D mesh is setup the same as V1b. Surface V1 tri count: 2,594,200 Surface V2 and V3 tri count: 1,155,555 Improvements on the triangulation will be based on the model results. Issues to look at: Capturing 1 meter feature for large region was difficult computational for Delaunay Is the 1 meter feature too small for modeling? Are the river faces oriented correctly for flow? Is the river positioned within the topology correctly for flow? What is the effect of river 1 cell wide, impact where there are more? The triangulation for this version was improved by: Use smarter refinement to reduce the total number of elements. Remove fan effects along parts of the river Break triangles inside branch intersects Improve smooth transition from small 1 meter cells near river to the boundary Improve triangulation along boundary to avoid mesh effects, allow smoothing V2.1--------------------------------------------------------------------- Title: LaGriT to Exodus number of dimension: 3 number of nodes: 645970 number of elements: 1155555 number of element blocks: 4 number of side sets: 6 Materials (Total 9 cells vertical): 1 1 cell .10 thick 2 1 cell .15 thick 3 3 cells .75 thick 4 4 cells 2.0 thick sideset tag: 1 nfaces: 128395 sideset tag: 2 nfaces: 128395 sideset tag: 3 nfaces: 7173 sideset tag: 4 nfaces: 6815 sideset tag: 5 nfaces: 13630 sideset tag: 6 nfaces: 107950 V2.0--------------------------------------------------------------------- Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 645970 number of elements: 1155555 number of element blocks: 2 number of side sets: 6 Materials (Total 9 cells vertical): 1 2 cells .25 thick 2 3 cells .75 thick 2 4 cells 2.0 thick ExodusII: Done writing to ExodusII file: prism_fs.exo Exodus side sets: 1 bottom 128395 faces 2 top full 128395 faces 3 all sides 7173 faces 4 top river 6815 faces 5 top river buffer 13630 faces 6 top not river 107950 faces Mesh stats: MIN MAX DIFFERENCE X -4.975466309E+02 2.124468018E+03 2.622014648E+03 Y -6.834048462E+02 1.823483887E+03 2.506888733E+03 Z 2.726891950E+03 3.323518000E+03 5.966260499E+02 Cell Volume 1.736938191E-02 1.275258663E+02 1.275084969E+02 Edge Min 1.250000000E-01 5.000000000E-01 3.750000000E-01 Edge Max 7.483734194E-01 3.225586790E+01 3.150749448E+01 Check for zero volumes: rmpoint/zero_volume/1.e-7 0 points removed and 0 elements removed. The point count is 645970. The element count is 2334780. -------------------------------------------- element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 297414 element norm edge length b/w .05 and .1 : 230508 element norm edge length b/w .1 and .2 : 235341 element norm edge length b/w .2 and .5 : 246609 element norm edge length b/w .5 and 1. : 145683 min max edge length = 0.7484E+00 max max edge length = 0.3226E+02 ----------------------------------------------------------- element norm min edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 513580 element norm edge length b/w .5 and 1. : 641975 min min edge length = 0.1250E+00 max min edge length = 0.5000E+00 ----------------------------------------------------------- Stacked Layers: surface name layer color type llcorner zic surf3.inp 1 2 -1 1) 2.726910E+03 refine 2 2 2 64598) 2.727410E+03 refine 3 2 2 129195) 2.727910E+03 refine 4 2 2 193792) 2.728410E+03 surf4.inp 5 2 0 258389) 2.728910E+03 refine 6 2 2 322986) 2.729160E+03 refine 7 2 2 387583) 2.729410E+03 surf5.inp 8 1 0 452180) 2.729660E+03 refine 9 1 2 516777) 2.729785E+03 tri_top.inp 10 1 -2 581374) 2.729910E+03 Elements per layer: 128395 stacked total: 1283950 Nodes per layer: 64597 stacked total: 645970 -----------------------------------------------------------