mesh generation for geological applications · mesh generation for geological applications carl w ....
TRANSCRIPT
![Page 1: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/1.jpg)
Mesh Generation for Geological Applications
Carl W. Gable
Earth & Environmental Sciences Division
http://meshing.lanl.gov
http://lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting
Golden CO June 23-27
![Page 2: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/2.jpg)
Slide 2
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Modeling and Simulation Workflow
Conceptual ModelSolid Model: Geometry & Topology
Mesh Generation
Mesh Optimization
Set Up: BC and IC
Physics Model
Results, Interpretation, Decision
![Page 3: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/3.jpg)
Slide 3
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Modeling and Simulation Workflow
Conceptual ModelSolid Model: Geometry & Topology
Mesh Generation
Mesh Optimization
Set Up: BC and IC
Physics Model
Results, Interpretation, Decision
![Page 4: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/4.jpg)
Slide 4
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Earth Science Related Mesh Generation
• Flow and Transport in Porous Media
• Tectonics, Stress and Strain, Faulting, Earthquake Studies
• Yucca Mountain Project (YMP)
• Heat & Mass Transport Studies of Nuclear Fuel In Salt Repository (GNEP)
• Nevada Test Site Underground Test Area (UGTA)
• Los Alamos Environmental Programs (LANL EP)
• Oil Shale and Water Resources
• Hard and Deeply Buried Targets (HDBT)
• WFO– Navy Coso Hydrothermal– Southern California Earthquake Center (SCEC)– Rhine Valley Fault Systems (NSF)
![Page 5: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/5.jpg)
Slide 5
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Conceptual Model
• Imagine a Spherical Cow
![Page 6: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/6.jpg)
Slide 6
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Conceptual Model
• Imagine a Spherical Cow
![Page 7: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/7.jpg)
Slide 7
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
• Geologic1. Geometry Poorly Known2. EarthVision, Stratamodel, gOcad,
GIS3. Very little FE integration4. Development driven by oil patch
(high end), environmental (low end).
5. No symmetry6. Often high aspect ratio
• Engineering1. Geometry Know w/ High Precision2. CAD Model (ACIS, Pro/E, Autocad)3. FE application well integrated4. Development driven by automotive,
aircraft, semiconductor, mechanical engineering, fluid dynamics
5. Often symmetry6. Generally unit aspect ratio
Solid Model: Geometry and Topology
![Page 8: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/8.jpg)
Slide 8
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Solid Model: Geometry and Topology
• Geologic1. Geometry Poorly Known2. EarthVision, Stratamodel, gOcad,
GIS3. Very little FE integration4. Development driven by oil patch
(high end), environmental (low end).
5. No symmetry6. Often high aspect ratio
• Engineering1. Geometry Know w/ High Precision2. CAD Model (ACIS, Pro/E, Autocad)3. FE application well integrated4. Development driven by automotive,
aircraft, semiconductor, mechanical engineering, fluid dynamics
5. Often symmetry6. Generally unit aspect ratio
![Page 9: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/9.jpg)
Slide 9
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Model: Geometry, No Topology
![Page 10: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/10.jpg)
Slide 10
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
San Andreas Fault
![Page 11: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/11.jpg)
Slide 11
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
What do we do if the geologic framework model is not in CAD format?
Option 11. Derive CAD solid model (ACIS, Pro/E, etc.) from
geologic framework model (EarthVision, gOcad, Stratamodel, …).
2. Use mesh generation tools that require CAD solid model input.
Option 25. Use mesh generation tools that do not require CAD
solid model input.
![Page 12: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/12.jpg)
Slide 12
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT: Los Alamos Grid Toolbox
meshing.lanl.gov lagrit.lanl.gov
![Page 13: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/13.jpg)
Slide 13
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT
• Mesh Generation
• Mesh Optimization
• Geometric Computing
• Element and Node Set Manipulation
• IO, Boundary Conditions, Initial Conditions
![Page 14: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/14.jpg)
Slide 14
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT: Mesh Generation
Delaunay point connection
– 2D triangulation
– 3D tetrahedralization
![Page 15: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/15.jpg)
Slide 15
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT: Mesh Generation
• 2D arbitrary (concave) polygon triangulation
![Page 16: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/16.jpg)
Slide 16
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT: Mesh Generation
• Block structured i,j,k connectivity
• 2.5D Stacking of z(x,y) surfaces
2x5
3x2
![Page 17: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/17.jpg)
Slide 17
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Orthogonal Finite Difference
• Logical Structured
• Block Structured
• Quadtree, Octree
• Unstructured (quad, tri, hex, tet, prism, pyramid, polyhedra)
• 2.5 D Unstructured
![Page 18: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/18.jpg)
Slide 18
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Orthogonal Finite Difference
![Page 19: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/19.jpg)
Slide 19
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Quadtree, Octree
![Page 20: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/20.jpg)
Slide 20
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Unstructured (quad, tri, hex, tet, prism, pyramid, polyhedra)
![Page 21: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/21.jpg)
Slide 21
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Same Geometry, Different Mesh Method
![Page 22: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/22.jpg)
Slide 22
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Unstructured 3D (hex, tet, hybrid)
![Page 23: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/23.jpg)
Slide 23
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Unstructured polyhedra (control volume methods)Voronoi Polygon
Median Polygon
![Page 24: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/24.jpg)
Slide 24
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
• Unstructured polyhedra3D Voronoi Polyhedra
![Page 25: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/25.jpg)
Slide 25
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Generation Depends on the Solver and Physics
![Page 26: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/26.jpg)
Slide 26
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Optimization
• Quality: good and bad is application dependent.
• Flexible tools to characterize and improve mesh quality– Volume, area, length
– Gradient of volume, area, length
– Aspect Ratio
– Angle, min/max angle, angle ratio, solid angle, dihedral angle
– Error estimate
– Solution value, solution gradient (d/dx, d/dx2 ), solution error
– User defined
![Page 27: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/27.jpg)
Slide 27
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Optimize: Refine, Derefine, Smooth
Original Block Triangles Filtered Block Triangles: Remove small area and high aspect ratio elements while maintaining geometry.
![Page 28: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/28.jpg)
Slide 28
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Optimize: Adapt mesh to user defined function
![Page 29: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/29.jpg)
Slide 29
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Boundary and Initial Conditions
• Mesh to Mesh Interpolation– Set material properties to nodes and/or elements (e.g. density)
– Set boundary conditions from another model (e.g. boundary flux)
• Mesh-Object Intersections– Point (injection), line (well bore), surface (fault), volume (tunnel)
• Point, Element, Face Set Operations (Intersection, Union, Complement)
![Page 30: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/30.jpg)
Slide 30
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Intersect fault surfaces with mesh to select elements to be refined with quadtree type mesh refinement.
Boundary and Initial Conditions
Mesh to Object Intersections
![Page 31: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/31.jpg)
Slide 31
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Landers/Hector Mine, Cut Away View
Aspect Ratio
DistanceFrom Faults
![Page 32: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/32.jpg)
Slide 32
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT Mesh Optimization
• Refine
– edge, face, element
– Rivara
• Refine
– field value, field gradient, aspect ratio, volume, arbitrary point
![Page 33: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/33.jpg)
Slide 33
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT Mesh Optimization
• Derefine
–edge length, volume, aspect
Original Elements Filter: Remove small area and high aspect ratio elements while maintaining geometry.
![Page 34: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/34.jpg)
Slide 34
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008LaGriT Mesh Optimization
• Smooth
– elliptic, laplace, random
• Massage - refine-derefine-smooth
![Page 35: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/35.jpg)
Slide 35
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008LaGriT Geometry
• Constructive Solid Geometry
surface/s_a/intrface/cylinder/x1,y1,z1/x2,y2,z2/radius surface/s_b/intrface/cylinder/x1,y1,z1/x2,y2,z2/radius surface/s_c/intrface/box/xmin,ymin,zmin/xmax,ymax,zmax/ region/r_ring/gt s_c and gt s_b and le s_a
![Page 36: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/36.jpg)
Slide 36
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT Element and Node Manipulation Commands
• Translate
• Scale
• Rotate: rotatept, rotateln
• Filter
• Perturb
• Remove
• Multi-Key Sort
• Reorder
• KDTree Search
![Page 37: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/37.jpg)
Slide 37
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT Grid Attributes
• Real and integer node and element attributes
• Element area, element volume
• aspect ratio, dihedral angle, solid angle, min/max edge length ratio
• Identify Sliver, Wedge, Needle, Cap
• Unit normal, area normal vector
• Synthetic normal to surface nodes
• Volume/Area integration of floating point node or cell attributes
• Dual mesh connectivity
• Voronoi and median volume and face area
![Page 38: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/38.jpg)
Slide 38
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
LaGriT Output Options
• PyLith
• GeoFEST
• FEHM
• AVS, GMV, Tecplot, (netCDF)
• STL, FLAG, X3D
![Page 39: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/39.jpg)
Slide 39
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
Mesh Manipulation
• Extract Lower D - 3D – 2D – 1D– Volume, Face, Line
• Extract 2D Surface (plane, isosurface, arbitrary triangulation) from 3D mesh
• Extract Line (well bore) from 3D mesh
![Page 40: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/40.jpg)
Slide 40
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008Mesh Manipulation
• Extract 2D Surface (plane) from 3D mesh.
![Page 41: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/41.jpg)
Slide 41
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
METIS Interface
• Supports METIS mesh partition algorithm calls
Partition:metis /partition/ metis_partmeshnodal / node / 32 metis /partition/ metis_partmeshdual / dual / 32
Reorder: metis / reorder / metis_edgend / dual metis / reorder / metis_nodend / node
![Page 42: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/42.jpg)
Slide 42
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
What is the LaGriT Interface?
• Command Line Driven
• Control File Driven
• Call from Fortran, C, C++
• Data structures can be accessed and manipulated by user code.
• Developer interface for extension and user modules
• Platforms: Linux, Mac, Sun, SGI
![Page 43: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/43.jpg)
Slide 43
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
What LaGriT Is Not
• No GUI interface
• No advancing front algorithm
• No interface for ACIS, Autocad, … CAD
• Not unstructured hex mesh (see Cubit)
![Page 44: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/44.jpg)
Slide 44
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008Finite Element Mesh With CFM Landers & Hector Mine Faults
Tet
Asp
ect R
atio
• Tetrahedral Finite Element MeshTetrahedral Finite Element Mesh
• 340km EW, 360km NS, 124km Z340km EW, 360km NS, 124km Z
• Element size ~50km far fieldElement size ~50km far field
• Element size ~500m near faultsElement size ~500m near faults
![Page 45: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/45.jpg)
Slide 45
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
No Silver Bullet
• With a flexible tool kit and the expertise necessary to know how and when to use different tools, we can adjust our approach depending upon the specific situation.
• A cookbook and the kitchen full of tools and ingredients does not make one a gourmet chef.
![Page 46: Mesh Generation for Geological Applications · Mesh Generation for Geological Applications Carl W . Gable Earth & Environmental Sciences Division ... • Point, Element, Face Set](https://reader030.vdocument.in/reader030/viewer/2022040920/5e98441f21658006b14c0dcb/html5/thumbnails/46.jpg)
Slide 46
Demo Problems: http://meshing.lanl.gov Manual & Documentation: lagrit.lanl.gov
Numerical Modeling of Crustal Deformation and Earthquake Faulting Golden 2008
•Demo problems with all input and LaGriT contol files can be found at the URL http://meshing.lanl.gov •Documentation and user manuals for LaGrit can be found at the URL http://lagrit.lanl.gov