IDEAS Jaramillo Watershed

Version History

Project Overview

DR project on creating a "disturbance science" demonstration on the Upper Jaramillo River using ATS. A method is developed to create a triangulation conforming to the watershed boundary and river banks. The DEM elevations are interpolated on to this triangulation to create the top of the mesh. Layers are translated below then filled with prism elements.
Mesh Directory: /scratch/tundra/tam/IDEAS/Jaramillo

General Mesh Images

Source DEM and Polygons

DEM Polygon Images and Description

DEM and Polygon Directory: /scratch/tundra/tam/IDEAS/Jaramillo/data

From: "Middleton, Richard" 
 The boundary of the Upper Jaramillo watershed
 The DEM
 The Jaramillo stream channel as a single line (Channel_single_part.txt)
 The stream channel as several individual lines (Channel_multipart.txt)
 A single polygon that is a 1m buffer around the Jaramillo River

DEM V1.2 Triangulation

This triangulation starts with the same polygons and triangulation as V1
Smooth, massage, and massage2 are used to optimize distribution.
A call to recon 1 creates a delaunay mesh.

DEM Triangulation V1.2 Images and Description  

Sep 14 tri_surf_dem_v1.2/tri_full_recon.inp
elements =  128395    nodes  =   91301 

                  MIN            MAX             DIFFERENCE 
edgemax  7.482054217E-01  2.974543903E+01   2.899723361E+01
area     1.389550553E-01  2.550517326E+02   2.549127775E+02
aratio   9.721532198E-02  9.999992999E-01   9.027839779E-01 

DEM V1 Triangulation

DEM Triangulation V1 Images and Description  

Aug 25 tri_surf_dem_v1/tri_flip.inp
elements =  259420    nodes  =  130169

                  MIN            MAX             DIFFERENCE
edgemax  7.483000780E-01  5.506619899E+01   5.431789891E+01
area     1.393946118E-01  9.571184741E+02   9.569790795E+02
aratio   9.272161037E-02  9.999995584E-01   9.072779480E-01

Exodus Mesh V3 Sep 2015

V3 Mesh Images and Description  
-rw-r--r-- 1 tamiller sft  62031716 Sep 14 17:26 EXO_full_v3/prism_fs.exo

Same as V2  with Richard's smooth 2 meter DEM interpolated to top
The DEM is a very smooth version of the original 1m DEM.
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.

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


Title: LAGRIT TO EXODUS                                                         
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

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   13627 faces
  6 top not river     107950 faces

Exodus Mesh V2.1 and V2.2 Jan 2016

V2.1 Mesh Images and Description  

V2.2 is V2.0 with modified layer spacing:
/scratch/tundra/tam/IDEAS/Jaramillo/EXO_full_v2.2

Title: LaGriT to Exodus 
number of nodes: 581373 
number of elements: 1027160 
number of element blocks: 3 
number of side sets: 6 

Materials (Total 8 cells vertical):
2 1 cell  .15 thick
3 3 cells .75 thick
4 4 cells 2.0 thick

V2.1 is V2.0 with modified layer spacing:
/scratch/tundra/tam/IDEAS/Jaramillo/EXO_full_v2.1

Title: LaGriT to Exodus 
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

Exodus Mesh V2 Sep 2015

V2 Mesh Images and Description  
-rw-r--r-- 1 tamiller sft  62031716 Sep 15 10:04 EXO_full_v2/prism_fs.exo

Mesh V2 with new tri surface and original 1 meter DEM
This mesh is translated to zero coordinate: XTRANS 360600   YTRANS 3973000.

The river is represented by a single cell width of 1 meter (material 1)
The river 1 meter wide and has a 1 meter land buffer (material 2)
The River node elevations are interpolated from a center line DEM
Land is material 3

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

Title: LAGRIT TO EXODUS                                                         
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

Exodus Mesh V1b Aug 2015

V1b Mesh Images and Description  

This fixes negative volumes by flipping the triangluated surface cells up.
Also modified top material 1 to have 2 cells instead of 3.

/scratch/tundra/tam/IDEAS/Jaramillo/EXO_full_v1b
-rw-r--r-- 1 tamiller sft 112248404 Aug 28 07:56 prism_fs.exo


Title: LAGRIT TO EXODUS
number of nodes:             1301690
number of elements:          2334780
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

Exodus side sets:
 1 bottom            259420 faces
 2 top full             259420 faces
 3 all sides               8244 faces
 4 top river               6815 faces
 5 top river buffer   13634 faces
 6 top not river     238981 faces

Exodus Mesh V1 Aug 2015

V1 Mesh Images and Description  

Mesh is in /scratch/tundra/tam/IDEAS/Jaramillo/EXO_full_v1
-rw-r--r-- 1 tamiller sft      123681228 Aug 27 08:35 prism_fs.exo

Unstructured stacked prism Mesh
River 1 meter wide
River buffer 1 meter wide
The River node elevations are interpolated from a center line DEM
Triangulation from .5 to 60 meters away from river

This is a first iteration triangulation using nicely formed polygons.
Improvements on the triangulation will be based on the model results.
Issues to look at:
 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 can be improved; remove fan effects along parts of the river,
break triangles inside branch intersects, improve smooth transition from small
1 meter cells near river to 60 meter at boundary, improve triangulation along
boundary to avoid ponding or other mesh efects.
See if triangulation of river can be improved to ensure single cells across width.


Title: LAGRIT TO EXODUS
number of nodes:             1431859
number of elements:          2594200
number of element blocks:          2
number of side sets:               6

Materials:
1 3 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            259420 faces
 2 top full          259420 faces
 3 all sides           9160 faces
 4 top river           6815 faces
 5 top river buffer   13634 faces
 6 top not river     238981 faces

Grid Team Contacts