an introduction to mesh generation algorithms part 1 · pdf fileoverview • the mesh...

Post on 06-Mar-2018

238 Views

Category:

Documents

1 Downloads

Preview:

Click to see full reader

TRANSCRIPT

An Introduction to Mesh Generation An Introduction to Mesh Generation AlgorithmsAlgorithms

Part 1Part 1

Mesh Generation Tutorial 1Instituto Tecnológico de Veracruz 21-25 April 2008

Part 1Part 1

Fathi El-YafiProject and Software Development Manager

Engineering Simulation

Overview

• The Mesh Generation Process• Meshing Algorithms�Tri / Tet Methods

�Quad / Hex Methods

Mesh Generation Tutorial 2Instituto Tecnológico de Veracruz 21-25 April 2008

�Quad / Hex Methods

�Hybrid Methods

�Surface Meshing

• Algorithm Characteristics

Mesh Generation

Geometry / Mesh

Mesh Generation Tutorial 3Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Engine

Mesh Generation Process

Set up sizing function for

Mesh surface

Smooth/Cleanup surface mesh

Verify Quality

For each surface

Mesh Vertices

Mesh Curves

Verify/correct for sizing criteria on

curves

Apply Manual Sizing, Match

Intervals

Mesh Generation Tutorial 4Instituto Tecnológico de Veracruz 21-25 April 2008

function for surface

surface surface mesh Quality

Set up sizing function for

volume

Mesh volume

Verify Quality

Smooth/Cleanup volume mesh

For each volume

The Mesh Generation Process

Meshing Algorithms

Mesh Generation Tutorial 5Instituto Tecnológico de Veracruz 21-25 April 2008

Tri/Tet Methods

� Octree� Advancing Front� Delaunay

Meshing Algorithms

Mesh Generation Tutorial 6Instituto Tecnológico de Veracruz 21-25 April 2008

http://www.ansys.com

Octree / Quadtree

Mesh Generation Tutorial 7Instituto Tecnológico de Veracruz 21-25 April 2008

� Define initial bounding box (root of Quadtree)� Recursively break into 4 leaves per root to resolve geometry� Find intersections of leaves with geometry boundary� Mesh each leaf using corners, side nodes and intersections with geometry� Delete Outside� (Yerry and Shephard, 84), (Shepherd and Georges, 91)

Octree / Quadtree

Mesh Generation Tutorial 8Instituto Tecnológico de Veracruz 21-25 April 2008

QMG, Cornell University

Octree / Quadtree

Mesh Generation Tutorial 9Instituto Tecnológico de Veracruz 21-25 April 2008

QMG, Cornell University

Advancing Front

C

Mesh Generation Tutorial 10Instituto Tecnológico de Veracruz 21-25 April 2008

A B

C

� Begin with boundary mesh - define as initial front� For each edge (face) on front, locate ideal node C based on front AB

C D

Advancing Front

Mesh Generation Tutorial 11Instituto Tecnológico de Veracruz 21-25 April 2008

A B

Cr

� Determine if any other nodes on current front are within search radius r of ideal location C (Choose D instead of C)

D

D

Advancing Front

Mesh Generation Tutorial 12Instituto Tecnológico de Veracruz 21-25 April 2008

� Book-Keeping: New front edges added and deleted from front as triangles are formed� Continue until no front edges remain on front

D

Advancing Front

Mesh Generation Tutorial 13Instituto Tecnológico de Veracruz 21-25 April 2008

� Book-Keeping: New front edges added and deleted from front as triangles are formed� Continue until no front edges remain on front

Advancing Front

Mesh Generation Tutorial 14Instituto Tecnológico de Veracruz 21-25 April 2008

� Book-Keeping: New front edges added and deleted from front as triangles are formed� Continue until no front edges remain on front

Advancing Front

Mesh Generation Tutorial 15Instituto Tecnológico de Veracruz 21-25 April 2008

� Book-Keeping: New front edges added and deleted from front as triangles are formed� Continue until no front edges remain on front

Cr

Advancing Front

Mesh Generation Tutorial 16Instituto Tecnológico de Veracruz 21-25 April 2008

A

B

C

�Where multiple choices are available, use best quality (closest shape to equilateral)�Reject any that would intersect existing front�Reject any inverted triangles (|AB X AC| > 0)�(Lohner,88;96)(Lo,91)

Advancing Front

Mesh Generation Tutorial 17Instituto Tecnológico de Veracruz 21-25 April 2008

Ansys, Inc.www.ansys.com

Delaunay

Mesh Generation Tutorial 18Instituto Tecnológico de Veracruz 21-25 April 2008

TriangleJonathon Shewchukhttp://www-2.cs.cmu.edu/~quake/triangle.html

Tetmesh-GHS3DINRIA, Francehttp://www.simulog.fr/tetmesh/

Circumcircle

Delaunay

Mesh Generation Tutorial 19Instituto Tecnológico de Veracruz 21-25 April 2008

Empty Circle (Sphere) Property: No other vertex is contained within the circumcircle (circumsphere) of any triangle (tetrahedron)

Delaunay

Mesh Generation Tutorial 20Instituto Tecnológico de Veracruz 21-25 April 2008

Delaunay Triangulation: Obeys empty-circle (sphere) property

Delaunay

Mesh Generation Tutorial 21Instituto Tecnológico de Veracruz 21-25 April 2008

Non-Delaunay Triangulation

X

Given a Delaunay Triangulation of n nodes, How do I insert node n+1 ?

Delaunay

Mesh Generation Tutorial 22Instituto Tecnológico de Veracruz 21-25 April 2008

Lawson Algorithm�Locate triangle containing X�Subdivide triangle�Recursively check adjoining triangles to ensure empty-circle property. Swap diagonal if needed�(Lawson,77)

X

Given a Delaunay Triangulation of n nodes, How do I insert node n+1 ?

Delaunay

Mesh Generation Tutorial 23Instituto Tecnológico de Veracruz 21-25 April 2008

Lawson Algorithm�Locate triangle containing X�Subdivide triangle�Recursively check adjoining triangles to ensure empty-circle property. Swap diagonal if needed�(Lawson,77)

Bowyer-Watson Algorithm�Locate triangle that contains the point�Search for all triangles whose circumcircle contain the point (d<r)

Xd

Delaunay

Mesh Generation Tutorial 24Instituto Tecnológico de Veracruz 21-25 April 2008

�Delete the triangles (creating a void in the mesh)�Form new triangles from the new point and the void boundary�(Watson,81)

r cd

Given a Delaunay Triangulation of n nodes, How do I insert node n+1 ?

X

Bowyer-Watson Algorithm�Locate triangle that contains the point�Search for all triangles whose circumcircle contain the point (d<r)

Delaunay

Mesh Generation Tutorial 25Instituto Tecnológico de Veracruz 21-25 April 2008

Given a Delaunay Triangulation of n nodes, How do I insert node n+1 ?

�Delete the triangles (creating a void in the mesh)�Form new triangles from the new point and the void boundary�(Watson,81)

Delaunay

Mesh Generation Tutorial 26Instituto Tecnológico de Veracruz 21-25 April 2008

�Begin with Bounding Triangles (or Tetrahedron)

Delaunay

Mesh Generation Tutorial 27Instituto Tecnológico de Veracruz 21-25 April 2008

�Insert boundary nodes using Delaunay method (Lawson or Bowyer-Watson)

Delaunay

Mesh Generation Tutorial 28Instituto Tecnológico de Veracruz 21-25 April 2008

�Insert boundary nodes using Delaunay method (Lawson or Bowyer-Watson)

Delaunay

Mesh Generation Tutorial 29Instituto Tecnológico de Veracruz 21-25 April 2008

�Insert boundary nodes using Delaunay method (Lawson or Bowyer-Watson)

Delaunay

Mesh Generation Tutorial 30Instituto Tecnológico de Veracruz 21-25 April 2008

�Insert boundary nodes using Delaunay method (Lawson or Bowyer-Watson)

Delaunay

Mesh Generation Tutorial 31Instituto Tecnológico de Veracruz 21-25 April 2008

�Insert boundary nodes using Delaunay method (Lawson or Bowyer-Watson)

Delaunay

Mesh Generation Tutorial 32Instituto Tecnológico de Veracruz 21-25 April 2008

�Recover boundary�Delete outside triangles�Insert internal nodes

h

Delaunay: Node Insertion

Mesh Generation Tutorial 33Instituto Tecnológico de Veracruz 21-25 April 2008

Grid Based� Nodes introduced based on a regular lattice� Lattice could be rectangular, triangular, quadtree, etc…� Outside nodes ignored

h

Delaunay: Node Insertion

Mesh Generation Tutorial 34Instituto Tecnológico de Veracruz 21-25 April 2008

Grid Based� Nodes introduced based on a regular lattice� Lattice could be rectangular, triangular, quadtree, etc…� Outside nodes ignored

Delaunay: Node Insertion

Mesh Generation Tutorial 35Instituto Tecnológico de Veracruz 21-25 April 2008

Centroïd�Nodes introduced at triangle centroïds�Continues until edge length, hl ≈

Delaunay: Node Insertion

Mesh Generation Tutorial 36Instituto Tecnológico de Veracruz 21-25 April 2008

l

Centroïd�Nodes introduced at triangle centroïds�Continues until edge length, hl ≈

Delaunay: Node Insertion

Mesh Generation Tutorial 37Instituto Tecnológico de Veracruz 21-25 April 2008

α

Circumcenter (“Guaranteed Quality”)� Nodes introduced at triangle circumcenters� Order of insertion based on minimum angle of any triangle� Continues until minimum angle > predefined minimum )30( o≈α

(Chew, Ruppert, Shewchuk)

Delaunay: Node Insertion

Mesh Generation Tutorial 38Instituto Tecnológico de Veracruz 21-25 April 2008

Circumcenter (“Guaranteed Quality”)� Nodes introduced at triangle circumcenters� Order of insertion based on minimum angle of any triangle� Continues until minimum angle > predefined minimum )30( o≈α

(Chew, Ruppert, Shewchuk)

A B

C

Delaunay: Node Insertion

Mesh Generation Tutorial 39Instituto Tecnológico de Veracruz 21-25 April 2008

C

Advancing Front� “Front” structure maintained throughout� Nodes introduced at ideal location from current front edge

(Marcum,95)

Delaunay: Node Insertion

Mesh Generation Tutorial 40Instituto Tecnológico de Veracruz 21-25 April 2008

Advancing Front� “Front” structure maintained throughout� Nodes introduced at ideal location from current front edge

(Marcum,95)

Delaunay: Node Insertion

Mesh Generation Tutorial 41Instituto Tecnológico de Veracruz 21-25 April 2008

Voronoï-Segment� Nodes introduced at midpoint of segment connecting the circumcircle centers of two adjacent triangles

(Rebay,93)

Delaunay: Node Insertion

Mesh Generation Tutorial 42Instituto Tecnológico de Veracruz 21-25 April 2008

Voronoï-Segment� Nodes introduced at midpoint of segment connecting the circumcircle centers of two adjacent triangles

(Rebay,93)

h

h

Delaunay: Node Insertion

Mesh Generation Tutorial 43Instituto Tecnológico de Veracruz 21-25 April 2008

h

h

Edges� Nodes introduced at along existing edges at l=h� Check to ensure nodes on nearby edges are not too close

(George,91)

Delaunay: Node Insertion

Mesh Generation Tutorial 44Instituto Tecnológico de Veracruz 21-25 April 2008

Edges� Nodes introduced at along existing edges at l=h� Check to ensure nodes on nearby edges are not too close

(George,91)

Delaunay: Boundary Constrained

Mesh Generation Tutorial 45Instituto Tecnológico de Veracruz 21-25 April 2008

Boundary Intersection� Nodes and edges introduced where Delaunay edges intersect boundary

Delaunay: Boundary Constrained

Mesh Generation Tutorial 46Instituto Tecnológico de Veracruz 21-25 April 2008

Boundary Intersection� Nodes and edges introduced where Delaunay edges intersect boundary

Delaunay: Boundary Constrained

Mesh Generation Tutorial 47Instituto Tecnológico de Veracruz 21-25 April 2008

(George,91)(Owen,99)

Local Swapping�Edges swapped between adjacent pairs of triangles until boundary is maintained

Delaunay: Boundary Constrained

Mesh Generation Tutorial 48Instituto Tecnológico de Veracruz 21-25 April 2008

(George,91)(Owen,99)

Local Swapping�Edges swapped between adjacent pairs of triangles until boundary is maintained

Delaunay: Boundary Constrained

Mesh Generation Tutorial 49Instituto Tecnológico de Veracruz 21-25 April 2008

(George,91)(Owen,99)

Local Swapping�Edges swapped between adjacent pairs of triangles until boundary is maintained

Delaunay: Boundary Constrained

Mesh Generation Tutorial 50Instituto Tecnológico de Veracruz 21-25 April 2008

(George,91)(Owen,99)

Local Swapping�Edges swapped between adjacent pairs of triangles until boundary is maintained

Delaunay: Boundary Constrained

Mesh Generation Tutorial 51Instituto Tecnológico de Veracruz 21-25 April 2008

(George,91)(Owen,99)

Local Swapping�Edges swapped between adjacent pairs of triangles until boundary is maintained

D C

VS

Delaunay: Boundary Constrained

Mesh Generation Tutorial 52Instituto Tecnológico de Veracruz 21-25 April 2008

Local Swapping ExampleRecover edge CD at vector Vs

VS

D C

E2

E3

E4E5E8

Delaunay: Boundary Constrained

Mesh Generation Tutorial 53Instituto Tecnológico de Veracruz 21-25 April 2008

E1

E2

E6E7

E8

Local Swapping ExampleMake a list (queue) of all edges Ei, that intersect Vs

D CE1

E2

E3

E4E5E8

Delaunay: Boundary Constrained

Mesh Generation Tutorial 54Instituto Tecnológico de Veracruz 21-25 April 2008

E2

E6E7

E8

Local Swapping ExampleSwap the diagonal of adjacent triangle pairs for each edge in the list

D CE3

E4E5E8

Delaunay: Boundary Constrained

Mesh Generation Tutorial 55Instituto Tecnológico de Veracruz 21-25 April 2008

E2

E5

E6E7

E8

Local Swapping ExampleCheck that resulting swaps do not cause overlapping triangles. If they do, then place edge at the back of the queue and try again later

D C

E3

E4E5E8

Delaunay: Boundary Constrained

Mesh Generation Tutorial 56Instituto Tecnológico de Veracruz 21-25 April 2008

E5

E6E7

E8

Local Swapping ExampleCheck that resulting swaps do not cause overlapping triangles. I they do, then place edge at the back of the queue and try again later

D C

E6

Delaunay: Boundary Constrained

Mesh Generation Tutorial 57Instituto Tecnológico de Veracruz 21-25 April 2008

E6

Local Swapping Example�Final swap will recover the desired edge.�Resulting triangle quality may be poor if multiple swaps were necessary�Does not maintain Delaunay criterion!

D E

BEdge Recovery• Force edges into triangulation by performing 2-3 swap transformation

DE = edge to be recovered

Delaunay: Boundary Constrained

Mesh Generation Tutorial 58Instituto Tecnológico de Veracruz 21-25 April 2008

A

C

E

ABC = non-conforming face

3D Local Swapping�Requires both boundary edge recovery and boundary face recovery

(George,91;99)(Owen,00)

B

D E2-3 Swap

Edge Recovery• Force edges into triangulation by performing 2-3 swap transformation

DE = edge to be recovered

Delaunay: Boundary Constrained

Mesh Generation Tutorial 59Instituto Tecnológico de Veracruz 21-25 April 2008

A

C

E

ABCEACBD

2-3 Swap

ABC = non-conforming face

3D Local Swapping�Requires both boundary edge recovery and boundary face recovery

(George,91;99)(Owen,00)

B

D E

Edge Recovery• Force edges into triangulation by performing 2-3 swap transformation

DE = edge recovered

2-3 Swap

Delaunay: Boundary Constrained

Mesh Generation Tutorial 60Instituto Tecnológico de Veracruz 21-25 April 2008

A

C

E

BAEDCBEDACED

ABCEACBD

2-3 Swap

3D Local Swapping�Requires both boundary edge recovery and boundary face recovery

(George,91;99)(Owen,00)

D E

BDE = edge recovered

Edge Recovery• Force edges into triangulation by performing 2-3 swap transformation

2-3 Swap

Delaunay: Boundary Constrained

Mesh Generation Tutorial 61Instituto Tecnológico de Veracruz 21-25 April 2008

A

C

E

3D Local Swapping�Requires both boundary edge recovery and boundary face recovery

(George,91;99)(Owen,00)

BAEDCBEDACED

ABCEACBD

2-3 Swap

A BS

Edge AB to be recovered

Delaunay: Boundary Constrained

Mesh Generation Tutorial 62Instituto Tecnológico de Veracruz 21-25 April 2008

A BS

3D Edge Recovery�Form queue of faces through which edge AB will pass�Perform 2-3 swap transformations on all faces in the list�If overlapping tets result, place back on queue and try again later�If still cannot recover edge, then insert “Steiner” point

Exploded view of tets intersected by AB

a

b

c

d e

a

bc

d

e

a

b

c

de

n2

a

b

n3n4

n5

n1

a

b

n1

n2

n3

n4

n5

2-3 Swap 2-2 Swap Face Split Edge Split Edge Suppress

abce, acbdaceb, adcb abce, acbd

Delaunay

Mesh Generation Tutorial 63Instituto Tecnológico de Veracruz 21-25 April 2008

a

b

c

d e

a

bc

d

e

a

b

c

def

n1

a

b

n2

n3n4

n5

n1

c

a

b

n3,i

12 3

n3,j

n3,k

n2,i

n1,i

n1,j

n1,kn2,j

n2,k

abde, bcde, cade

aceb, adcb

adeb, edcb

abce, acbd

abfe, bcfe, cafebafd, cbfd, acfd

abnini+1 {i=1…N}

abnini+1, cbnini+1 {i=1…N}

abnini+1 {i=1…N}

nm,knm.jnm.ia, nm,knm.jnm.ib {m=1…M}

N = no. adj. tets at edge abM = no. unique trias in polygon P={n1,n2,n3}

Quad/HexMethods

Structured Unstructured

Mesh: Structured / Unstructured

Mesh Generation Tutorial 64Instituto Tecnológico de Veracruz 21-25 April 2008

Structured Unstructured�No specific requirements for geometry�quads/hexes placed to conform to geometry.�No connectivity requirement (although optimization of connectivity is beneficial)

� Requires geometry to conform to specific characteristics� Regular patterns of quads/hexes formed based on characteristics of geometry� Internal nodes always attached to same number of elements

63

Structured: Mapped Meshing

Mesh Generation Tutorial 65Instituto Tecnológico de Veracruz 21-25 April 2008

6

6

3� 4 topological sides� Opposite sides must have similar intervals

Geometry Requirements

Algorithm

�Trans-finite Interpolation (TFI)� Maps a regular lattice of quads onto polygon

(Thompson,88;99), (Cook,82)

Structured: 3D Mapped Meshing

Mesh Generation Tutorial 66Instituto Tecnológico de Veracruz 21-25 April 2008

� 6 topological surfaces� opposite surfaces must have similar mapped meshes

Geometry Requirements

Block-Structured

Structured: Mapped Meshing

Mesh Generation Tutorial 67Instituto Tecnológico de Veracruz 21-25 April 2008

http://www.gridpro.com/gridgallery/tmachinery.html http://www.pointwise.com/case/747.htm

i

j

+3i∑ = 0iInterval

∑ = 0jInterval

Structured: Sub-Mapping

Mesh Generation Tutorial 68Instituto Tecnológico de Veracruz 21-25 April 2008

+3i

+2j

q -1j

� Blocky-type surfaces (principally 90 degree angles)

Geometry Requirements

(White,95)

Structured: Sub-Mapping

Mesh Generation Tutorial 69Instituto Tecnológico de Veracruz 21-25 April 2008

� Automatically decomposes surface into mapable regions based on assigned intervals

(White,95)

Structured: Sweeping

Mesh Generation Tutorial 70Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

target

linking surfaces

Structured: Sweeping

Mesh Generation Tutorial 71Instituto Tecnológico de Veracruz 21-25 April 2008

source

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 72Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 73Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 74Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 75Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 76Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 77Instituto Tecnológico de Veracruz 21-25 April 2008

Geometry Requirements� Source and target surfaces topologically similar� Linking surfaces mapable or submapable

Structured: Sweeping

Mesh Generation Tutorial 78Instituto Tecnológico de Veracruz 21-25 April 2008

Cubit, Sandia National Labs

Gambit, Fluent Inc.

Sweep Direction

1-to-1 sweepable

Sweeping

Mesh Generation Tutorial 79Instituto Tecnológico de Veracruz 21-25 April 2008

Sweep Direction

n-to-1 sweepable

Sweep Direction

Sweeping

Mesh Generation Tutorial 80Instituto Tecnológico de Veracruz 21-25 April 2008

Sweep Direction

n-to-m sweepableMulti-Sweep

Sweep Direction

Sweeping

Mesh Generation Tutorial 81Instituto Tecnológico de Veracruz 21-25 April 2008

Sweep Direction

Sweep Direction

Sweeping

Mesh Generation Tutorial 82Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 83Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 84Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 85Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 86Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 87Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 88Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 89Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Sweeping

Mesh Generation Tutorial 90Instituto Tecnológico de Veracruz 21-25 April 2008

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

The fundamental strategy of multi-sweep is to convert an n-to-m sweepable volume into a

number of n-to-1 sweepable volumes.

Sweep Direction

Decomposition / Sweep Overview

Mesh Generation Tutorial 91Instituto Tecnológico de Veracruz 21-25 April 2008

Sweep Direction

Decomposition / Sweep Overview

Mesh Generation Tutorial 92Instituto Tecnológico de Veracruz 21-25 April 2008

Sweeping

Mesh Generation Tutorial 93Instituto Tecnológico de Veracruz 21-25 April 2008

(White, 2004)CCSweep

Medial Axis

Mesh Generation Tutorial 94Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 95Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 96Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 97Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 98Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 99Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 100Instituto Tecnológico de Veracruz 21-25 April 2008

(Price, 95;97)(Tam,91)

� Medial Object - Roll a Maximal circle or sphere through the model. The center traces the medial object � Medial Object used as a tool to automatically decompose model into simpler mapable or sweepable parts

Medial Axis

Mesh Generation Tutorial 101Instituto Tecnológico de Veracruz 21-25 April 2008

Medial Axis + Midpoint Subdivision (Price, 95) (Sheffer, 98)

Embedded Voronoï Graph

Meshing Algorithms

Mesh Generation Tutorial 102Instituto Tecnológico de Veracruz 21-25 April 2008

Unstructured: Indirect Quad

Mesh Generation Tutorial 103Instituto Tecnológico de Veracruz 21-25 April 2008

Triangle splitting� Each triangle split into 3 quads� Typically results in poor angles

Unstructured: Indirect Hex

Mesh Generation Tutorial 104Instituto Tecnológico de Veracruz 21-25 April 2008

(Taniguchi, 96)

Tetrahedron splitting� Each tetrahedron split into 4 hexahedra� Typically results in poor angles

Unstructured: Indirect Hex

Mesh Generation Tutorial 105Instituto Tecnológico de Veracruz 21-25 April 2008

(Taniguchi, 96)

Tetrahedron splitting� Each tetrahedron split into 4 hexahedra� Typically results in poor angles

Unstructured: Indirect Hex

Mesh Generation Tutorial 106Instituto Tecnológico de Veracruz 21-25 April 2008

(Taniguchi, 96)

Tetrahedron splitting� Each tetrahedron split into 4 hexahedra� Typically results in poor angles

Unstructured: Indirect Hex

Mesh Generation Tutorial 107Instituto Tecnológico de Veracruz 21-25 April 2008

Example of geometry meshed by tetrahedron splitting

Unstructured: Indirect Quad

Mesh Generation Tutorial 108Instituto Tecnológico de Veracruz 21-25 April 2008

Triangle Merge� Two adjacent triangles combined into a single quad � Test for best local choice for combination� Triangles can remain if attention is not paid to order of combination

Unstructured: Indirect Quad

Mesh Generation Tutorial 109Instituto Tecnológico de Veracruz 21-25 April 2008

Triangle Merge� Two adjacent triangles combined into a single quad � Test for best local choice for combination� Triangles can remain if attention is not paid to order of combination

Unstructured: Indirect Quad

Mesh Generation Tutorial 110Instituto Tecnológico de Veracruz 21-25 April 2008

Triangle Merge� Two adjacent triangles combined into a single quad � Test for best local choice for combination� Triangles can remain if attention is not paid to order of combination

Directed Triangle Merge� Merging begins at a boundary� Advances from one set of triangles to the next� Attempts to maintain even number of

Unstructured: Indirect Quad

Mesh Generation Tutorial 111Instituto Tecnológico de Veracruz 21-25 April 2008

� Attempts to maintain even number of intervals on any loop� Can produce all-quad mesh� Can also incorporate triangle splitting(Lee and Lo, 94)

Triangle Merge w/ local transformations (“Q-Morph”)� Uses an advancing front approach� Local swaps applied to improve resulting quad� Any number of triangles

Unstructured: Indirect Quad

Mesh Generation Tutorial 112Instituto Tecnológico de Veracruz 21-25 April 2008

A B

NB

e e

NA

e e C

D� Any number of triangles merged to create a quad� Attempts to maintain even number of intervals on any loop� Produces all-quad mesh from even intervals(Owen, 99)

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 113Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 114Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 115Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 116Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 117Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph

Unstructured: Indirect Quad

Mesh Generation Tutorial 118Instituto Tecnológico de Veracruz 21-25 April 2008

Unstructured: Indirect Quad

Mesh Generation Tutorial 119Instituto Tecnológico de Veracruz 21-25 April 2008

Q-Morph Lee,Lo Method

AB

CD

AB

CD

E

AB

CD

EF

Indirect

Mesh Generation Tutorial 120Instituto Tecnológico de Veracruz 21-25 April 2008

B B B

AB

CD

EF

G

AB

CD

EF

GH

AB

CD

EF

GH

Tetrahedral Merge w/ local transformations (“H-Morph ”)

Unstructured-Hex

Mesh Generation Tutorial 121Instituto Tecnológico de Veracruz 21-25 April 2008

H-Morph

“Hex-Dominant Meshing”(Owen and Saigal, 00)

Unstructured-Hex

Mesh Generation Tutorial 122Instituto Tecnológico de Veracruz 21-25 April 2008

H-Morph

“Hex-Dominant Meshing”(Owen and Saigal, 00)

Unstructured-Hex

Mesh Generation Tutorial 123Instituto Tecnológico de Veracruz 21-25 April 2008

Meshing Algorithms

Mesh Generation Tutorial 124Instituto Tecnológico de Veracruz 21-25 April 2008

Unstructured-Hex

Mesh Generation Tutorial 125Instituto Tecnológico de Veracruz 21-25 April 2008

Grid-Based

� Generate regular grid of quads/hexes on the interior of model� Fit elements to the boundary by projecting interior faces towards the surfaces� Lower quality elements near boundary� Non-boundary conforming

Unstructured-Hex

Mesh Generation Tutorial 126Instituto Tecnológico de Veracruz 21-25 April 2008

Grid-Based

� Generate regular grid of quads/hexes on the interior of model� Fit elements to the boundary by projecting interior faces towards the surfaces� Lower quality elements near boundary� Non-boundary conforming

Unstructured-Hex

Mesh Generation Tutorial 127Instituto Tecnológico de Veracruz 21-25 April 2008

Grid-Based

� Generate regular grid of quads/hexes on the interior of model� Fit elements to the boundary by projecting interior faces towards the surfaces� Lower quality elements near boundary� Non-boundary conforming

Unstructured-Hex

Mesh Generation Tutorial 128Instituto Tecnológico de Veracruz 21-25 April 2008

Grid-Based

� Generate regular grid of quads/hexes on the interior of model� Fit elements to the boundary by projecting interior faces towards the surfaces� Lower quality elements near boundary� Non-boundary conforming

(Schneiders,96)

Unstructured-Hex

Mesh Generation Tutorial 129Instituto Tecnológico de Veracruz 21-25 April 2008

Grid-Based

Unstructured-Hex

Mesh Generation Tutorial 130Instituto Tecnológico de Veracruz 21-25 April 2008

http://www.numeca.be/hexpress_home.html

Grid-Based

Gambit, Fluent, Inc.

Direct Quad

Mesh Generation Tutorial 131Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Unstructured-Quad

Mesh Generation Tutorial 132Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Unstructured-Quad

Mesh Generation Tutorial 133Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Form new row and check for overlap

Unstructured-Quad

Mesh Generation Tutorial 134Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Insert “Wedge”

Unstructured-Quad

Mesh Generation Tutorial 135Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Seams

Unstructured-Quad

Mesh Generation Tutorial 136Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Close Loops and smooth

Unstructured-Quad

Mesh Generation Tutorial 137Instituto Tecnológico de Veracruz 21-25 April 2008

Paving� Advancing Front: Begins with front at boundary� Forms rows of elements based on front angles� Must have even number of intervals for all-quad mesh

(Blacker,92)(Cass,96)

Unstructured-Hex

Mesh Generation Tutorial 138Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 139Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 140Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 141Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 142Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 143Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Unstructured-Hex

Mesh Generation Tutorial 144Instituto Tecnológico de Veracruz 21-25 April 2008

Plastering

(Blacker, 93)� 3D extension of “paving”� Row-by row or element-by-element

Exterior Hex mesh Remaining Void

Unstructured-Hex

Mesh Generation Tutorial 145Instituto Tecnológico de Veracruz 21-25 April 2008

Ford Crankshaft

Plastering + Tet Meshing

“Hex-Dominant Meshing”

Direct

Mesh Generation Tutorial 146Instituto Tecnológico de Veracruz 21-25 April 2008

Whisker Weaving

� First constructs dual of the quad/hex mesh� Inserts quad/hex at the intersections of the dual chords

Direct

Mesh Generation Tutorial 147Instituto Tecnológico de Veracruz 21-25 April 2008

Whisker Weaving�Spatial Twist Continuum - Dual of a 3D hex mesh (Murdoch, 96)�Hexes formed at intersection of twist planes

Twist Plane

Direct

Mesh Generation Tutorial 148Instituto Tecnológico de Veracruz 21-25 April 2008

Whisker Weaving�Spatial Twist Continuum - Dual of a 3D hex mesh (Murdoch, 96)�Hexes formed at intersection of twist planes

Twist Plane

Direct

Mesh Generation Tutorial 149Instituto Tecnológico de Veracruz 21-25 April 2008

Whisker Weaving�Spatial Twist Continuum - Dual of a 3D hex mesh (Murdoch, 96)�Hexes formed at intersection of twist planes

Twist Planes

Direct

Mesh Generation Tutorial 150Instituto Tecnológico de Veracruz 21-25 April 2008

Whisker Weaving�Spatial Twist Continuum - Dual of a 3D hex mesh (Murdoch, 96)�Hexes formed at intersection of twist planes

Whisker Weaving

� Define the topology of the twist planes using whisker diagrams� Each whisker diagram represents a closed loop of the surface dual

Direct

Mesh Generation Tutorial 151Instituto Tecnológico de Veracruz 21-25 April 2008

the surface dual� Each boundary vertex on the diagram represents a quad face on the surface� Objective is to resolve internal connectivity by “weaving” the chords following a set of basic rules

(Tautges,95;96)Whisker diagrams used to resolve hex mesh above

Hex Meshing Research

Unconstrained Paving

Remove constraint of pre-meshed

Mesh Generation Tutorial 152Instituto Tecnológico de Veracruz 21-25 April 2008

Remove constraint that we must define number of quad when row is advanced.

This constrains only 1 DOF.

Remove constraint of pre-meshed boundary.

Unconstrained PavingEach Row Advancement Constrains Only 1 DOF

Hex Meshing Research

Mesh Generation Tutorial 153Instituto Tecnológico de Veracruz 21-25 April 2008

Quads are only completely defined when 2 unconstrained rows cross

Quad Elements

Unconstrained PavingEach Row Advancement Constrains Only 1 DOF

Hex Meshing Research

Mesh Generation Tutorial 154Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 155Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering(DOF = 2)

Hex Meshing Research

Mesh Generation Tutorial 156Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering(DOF = 1)

Hex Meshing Research

Mesh Generation Tutorial 157Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 158Instituto Tecnológico de Veracruz 21-25 April 2008

(DOF = 0)

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 159Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 160Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 161Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 162Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 163Instituto Tecnológico de Veracruz 21-25 April 2008

Unconstrained Plastering

Hex Meshing Research

Mesh Generation Tutorial 164Instituto Tecnológico de Veracruz 21-25 April 2008

CFD Meshing

Hybrid Methods

Mesh Generation Tutorial 165Instituto Tecnológico de Veracruz 21-25 April 2008

Image courtesy of acelab, University of Texas, Austin, http://acelab.ae.utexas.edu

Image courtesy of Roy P. Koomullil, Engineering Research Center, Mississippi State University, http://www.erc.msstate.edu/~roy/

Advancing Layers Method

Hybrid Methods

Mesh Generation Tutorial 166Instituto Tecnológico de Veracruz 21-25 April 2008

Advancing Layers Method

Hybrid Methods

Mesh Generation Tutorial 167Instituto Tecnológico de Veracruz 21-25 April 2008

Discretized Boundary

(Pirzadeh, 1994)Advancing Layers Method

Hybrid Methods

Mesh Generation Tutorial 168Instituto Tecnológico de Veracruz 21-25 April 2008

Define Normals at boundary nodes

distribution function

Advancing Layers Method

Hybrid Methods

Mesh Generation Tutorial 169Instituto Tecnológico de Veracruz 21-25 April 2008

Generate nodes along normals according to distribution function Form layer

Distance from wall

Ele

men

t si

ze

Advancing Layers Method

distribution function

Hybrid Methods

Mesh Generation Tutorial 170Instituto Tecnológico de Veracruz 21-25 April 2008

Distance from wall

Ele

men

t si

ze

Generate nodes along normals according to distribution function Form layer

Advancing Layers Method

distribution function

Hybrid Methods

Mesh Generation Tutorial 171Instituto Tecnológico de Veracruz 21-25 April 2008

Distance from wall

Ele

men

t si

ze

Generate nodes along normals according to distribution function Form layer

Advancing Layers Method

Hybrid Methods

Mesh Generation Tutorial 172Instituto Tecnológico de Veracruz 21-25 April 2008

Define new boundary for triangle mesher

Hybrid Methods

Mesh Generation Tutorial 173Instituto Tecnológico de Veracruz 21-25 April 2008

Mesh with triangles

Hybrid Methods

Mesh Generation Tutorial 174Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Hybrid Methods

Mesh Generation Tutorial 175Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Hybrid Methods

Mesh Generation Tutorial 176Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Hybrid Methods

Mesh Generation Tutorial 177Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Blend Regions

Hybrid Methods

Mesh Generation Tutorial 178Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Blend Regions

Hybrid Methods

Mesh Generation Tutorial 179Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Blend Regions

Hybrid Methods

Mesh Generation Tutorial 180Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Smoothed Normals

Hybrid Methods

Mesh Generation Tutorial 181Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Smoothed Normals

Hybrid Methods

Mesh Generation Tutorial 182Instituto Tecnológico de Veracruz 21-25 April 2008

Convex Corner Concave Corner

Smoothed Normals

Define Normals

Hybrid Methods

Mesh Generation Tutorial 183Instituto Tecnológico de Veracruz 21-25 April 2008

Multiple Normals

αDefine Normals every a degrees

Hybrid Methods

Mesh Generation Tutorial 184Instituto Tecnológico de Veracruz 21-25 April 2008

Multiple Normals

Hybrid Methods

Mesh Generation Tutorial 185Instituto Tecnológico de Veracruz 21-25 April 2008

Intersecting Boundary Layers

Hybrid Methods

Mesh Generation Tutorial 186Instituto Tecnológico de Veracruz 21-25 April 2008

Intersecting Boundary Layers

Hybrid Methods

Mesh Generation Tutorial 187Instituto Tecnológico de Veracruz 21-25 April 2008

Intersecting Boundary Layers

Delete overlapping elements

Hybrid Methods

Mesh Generation Tutorial 188Instituto Tecnológico de Veracruz 21-25 April 2008

Intersecting Boundary Layers

Hybrid Methods

Mesh Generation Tutorial 189Instituto Tecnológico de Veracruz 21-25 April 2008

Image courtesy of SCOREC, Rensselaer Polytechnic Institute, http://www.scorec.rpi.edu/ (Garimella, Shephard, 2000)

Non-Conforming Diagonal Edge

Hex-Tet Interface

Mesh Generation Tutorial 190Instituto Tecnológico de Veracruz 21-25 April 2008

Conforming quad-triangle Conforming hex-tet?

Non-Conforming Node

Solutions� Free Edge (Non-conforming)

Hex-Tet Interface

Mesh Generation Tutorial 191Instituto Tecnológico de Veracruz 21-25 April 2008

� Free Edge (Non-conforming)� Multi-point Constraint� Pyramid

Hex-Tet Interface

Mesh Generation Tutorial 192Instituto Tecnológico de Veracruz 21-25 April 2008

Heat sink meshed with hexes, tets and pyramids

Hex-Tet Interface

Mesh Generation Tutorial 193Instituto Tecnológico de Veracruz 21-25 April 2008

Pyramid Elements for maintaining compatibility between hex and tet elements (Owen,00)

A,B N1

N2

N3N4

N5 A,B N1

N2

N3N4

N5

(a) (b)

N4 N3

N2

Hex-Tet Interface

Mesh Generation Tutorial 194Instituto Tecnológico de Veracruz 21-25 April 2008

Tetrahedral transformations to form Pyramids�Use 2-3 swaps to obtain 2 tets at diagonal�combine 2 tets to form pyramid

A,B N1

N3N4

N5 A,B N1

N3

N5

(a) (b)

(c)A

B

N5 N1

(d)

A,B N1

N2

N3N4

N5 A,B N1

N2

N3N4

N5

(a) (b)

N4 N3

N2

Hex-Tet Interface

Mesh Generation Tutorial 195Instituto Tecnológico de Veracruz 21-25 April 2008

A,B N1

N3N4

N5 A,B N1

N3

N5

(a) (b)

(c)A

B

N5 N1

(d)

Tetrahedral transformations to form Pyramids�Use 2-3 swaps to obtain 2 tets at diagonal�combine 2 tets to form pyramid

A,B N1

N2

N3N4

N5 A,B N1

N2

N3N4

N5

(a) (b)

N4 N3

Hex-Tet Interface

Mesh Generation Tutorial 196Instituto Tecnológico de Veracruz 21-25 April 2008

A,B N1

N3N4

N5 A,B N1

N3

N5

(a) (b)

(c)A

B

N5 N1

(d)

Tetrahedral transformations to form Pyramids�Use 2-3 swaps to obtain 2 tets at diagonal�combine 2 tets to form pyramid

A,B N1

N2

N3N4

N5 A,B N1

N2

N3N4

N5

(a) (b)

N3

Hex-Tet Interface

Mesh Generation Tutorial 197Instituto Tecnológico de Veracruz 21-25 April 2008

(d)A,B N1

N3N4

N5 A,B N1

N3

N5

(a) (b)

(c)A

B

N5 N1

Tetrahedral transformations to form Pyramids�Use 2-3 swaps to obtain 2 tets at diagonal�combine 2 tets to form pyramid

A,B N1

N2

N3N4

N5 A,B N1

N2

N3N4

N5

(a) (b)

N3

Hex-Tet Interface

Mesh Generation Tutorial 198Instituto Tecnológico de Veracruz 21-25 April 2008

(d)A,B N1

N3N4

N5 A,B N1

N3

N5

(a) (b)

(c)A

B

N5 N1

Tetrahedral transformations to form Pyramids�Use 2-3 swaps to obtain 2 tets at diagonal�combine 2 tets to form pyramid

N1N

A

N2N3

Pyramid Open Method

N1A

C

N2N3

N

Hex-Tet Interface

Mesh Generation Tutorial 199Instituto Tecnológico de Veracruz 21-25 April 2008

N1N4

B

Non-Conforming Condition:Tets at quad diagonal A-B N

N

N14

B

A

C

N23

� Insert C at midpoint AB:� Split all tets at edge AB

B

N1N4

� Move C to average N1,N2…Nn� Create New Pyramid A,Nn,B,N1,C

Surface Meshing

Direct 3D Meshing Parametric Space Meshing

Mesh Generation Tutorial 200Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

� Elements formed in 3D using actual x-y-z representation of surface

� Elements formed in 2D using parametric representation of surface� Node locations later mapped to 3D

A

B

Surface Meshing

Mesh Generation Tutorial 201Instituto Tecnológico de Veracruz 21-25 April 2008

3D Surface Advancing Front� Form triangle from front edge AB

A

B

C

NC Tangent plane

Surface Meshing

Mesh Generation Tutorial 202Instituto Tecnológico de Veracruz 21-25 April 2008

3D Surface Advancing Front� Define tangent plane at front by averaging normals at A and B

A

B

C

NC

D

Surface Meshing

Mesh Generation Tutorial 203Instituto Tecnológico de Veracruz 21-25 April 2008

3D Surface Advancing Front� define D to create ideal triangle on tangent plane

A

B

C

NC

Surface Meshing

Mesh Generation Tutorial 204Instituto Tecnológico de Veracruz 21-25 April 2008

D

3D Surface Advancing Front� project D to surface (find closest point on surface)

3D Surface Advancing Front� Must determine overlapping or intersecting triangles in 3D. (Floating point robustness issues)� Extensive use of geometry evaluators (for

Surface Meshing

Mesh Generation Tutorial 205Instituto Tecnológico de Veracruz 21-25 April 2008

� Extensive use of geometry evaluators (for normals and projections)� Typically slower than parametric implementations� Generally higher quality elements� Avoids problems with poor parametric representations (typical in many CAD environments)( Lo,96;97); (Cass,96)

Parametric Space Mesh Generation

� Parameterization of the NURBS provided by the CAD model can be used to reduce the mesh generation to 2D

Surface Meshing

Mesh Generation Tutorial 206Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

u

v

� Isotropic: Target element shapes are equilateral triangles� Equilateral elements in parametric space may be distorted when mapped to 3D space.� If parametric space resembles 3D space v

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 207Instituto Tecnológico de Veracruz 21-25 April 2008

� If parametric space resembles 3D space without too much distortion from u-vspace to x-y-zspace, then isotropic methods can be used. u

v

u

v

� Parametric space can be “customized” or warped so that isotropic methods can be used.� Works well for many cases.� In general, isotropic mesh generation does not work well for parametric meshing

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 208Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

u

v

Warped parametric space

v

� Anisotropic: Triangles are stretched based on a specified vector field�Triangles appear stretched in 2d (parametric space), but are near equilateral in 3D

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 209Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

u

v

space), but are near equilateral in 3D

� Stretching is based on field of surface derivatives

=∆z

v

y

v

x

v

δδ

δδ

δδ

,,v

=∆z

u

y

u

x

u

δδ

δδ

δδ

,,u

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 210Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

x

y

z

z

u

y

u

x

u

δδ

δδ

δδ

,,

z

v

y

v

x

v

δδ

δδ

δδ

,,

uu ∆⋅∆=E vu ∆⋅∆=F vv ∆⋅∆=G

=GF

FE)(XM

� Metric, M can be defined at every location on surface. Metric at location X is:

� Distances in parametric space can now be measured as a function of direction and location on the surface. Distance from point X to Q is defined as:

XQXQXQl T )()( XM≈

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 211Instituto Tecnológico de Veracruz 21-25 April 2008

u

v

x

y

z

X

Q

)(XQl

u

v

X Q)(XQl

M(X)

� Use essentially the same isotropic methods for 2D mesh generation, except distances and angles are now measured with respect to the local metric tensor M(X).� Can use Delaunay (George, 99) or Advancing Front Methods (Tristano,98)

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 212Instituto Tecnológico de Veracruz 21-25 April 2008

� Is generally faster than 3D methods� Is generally more robust (No 3D intersection calculations)� Poor parameterization can cause

Parametric Space Mesh Generation

Surface Meshing

Mesh Generation Tutorial 213Instituto Tecnológico de Veracruz 21-25 April 2008

� Poor parameterization can cause problems� Not possible if no parameterization is provided� Can generate your own parametric space (Flatten 3D surface into 2D) (Marcum, 99) (Sheffer,00)

1. Conforming Mesh� Elements conform to a prescribed surface mesh

2. Boundary Sensitive� Rows/layers of elements roughly conform to the contours of the boundary

3. Orientation Insensitive

Algorithm Characteristics

Mesh Generation Tutorial 214Instituto Tecnológico de Veracruz 21-25 April 2008

3. Orientation Insensitive� Rotating/Scaling geometry will not change the resulting mesh

4. Regular Node Valence� Inherent in the algorithm is the ability to maintain (nearly) the same number of elements adjacent each node)

5. Arbitrary Geometry� The algorithm does not rely on a specific class/shape of geometry

6. Commercial Viability (Robustness/Speed)� The algorithm has been used in a commercial setting

Ma

ppe

dS

ub-m

ap

Tri

Spl

itT

ri M

erg

eQ

-Mor

phG

rid-B

ase

dM

edi

al A

xis

Pa

ving

Ma

ppe

dS

ub-m

ap

Sw

ee

ping

Tet S

plit

Tet M

erg

eH

-Mor

phG

rid-B

ase

dM

edi

al S

urf.

Pla

ste

ring

Whi

ske

r W

.

Tris Hexes

Oct

ree

De

laun

ay

Adv

. Fro

nt

QuadsTets

Conforming 1

Algorithm Characteristics

Mesh Generation Tutorial 215Instituto Tecnológico de Veracruz 21-25 April 2008

Conforming Mesh

Boundary Sensitive

Orientation Insensitive

Commercially Viability

Regular Node Valence

Arbitrary Geometry

1

2

3

4

5

6

Blacker, Ted D. and R. J. Myers,. “Seams and Wedges in Plastering: A 3D Hexahedral Mesh Generation Algorithm,” Engineering With Computers, 2, 83-93 (1993)

Blacker, Ted D., “The Cooper Tool”, Proceedings, 5th International Meshing Roundtable, 13-29 (1996)

Blacker, Ted D., and Michael B. Stephenson. “Paving: A New Approach to Automated Quadrilateral Mesh Generation”, International Journal for Numerical Methods in Engineering, 32, 811-847 (1991)

Canann, Scott A., Joseph R. Tristano and Matthew L. Staten, “An approach to combined Laplacian and optimization-based smoothing for triangular, quadrilateral and quad-dominant meshes,” Proc. 7th International Meshing Roundtable, pp.479-494 (1998)

Canann, S. A., S. N. Muthukrishnan, and R. Phillips, “Topological refinement procedures for triangular finite element meshes,” Engineering with Computers, 12 (3/4), 243-255 (1996)

Cass, Roger J, Steven E. Benzley, Ray J. Meyers and Ted D. Blacker. “Generalized 3-D Paving: An Automated Quadrilateral Surface Mesh Generation Algorithm”, International Journal for Numerical Methods in Engineering, 39, 1475-1489 (1996)

References

Mesh Generation Tutorial 216Instituto Tecnológico de Veracruz 21-25 April 2008

Algorithm”, International Journal for Numerical Methods in Engineering, 39, 1475-1489 (1996)

Chew, Paul L., "Guaranteed-Quality Triangular Meshes", TR 89-983, Department of Computer Science, Cornell University, Ithaca, NY, (1989)

Cook, W.A., and W.R. Oakes, “Mapping Methods for Generating Three-Dimensional Meshes”, Computers in Mechanical Engineering, August, 67-72 (1982)

CUBIT Mesh Generation Toolkit, URL:http://endo.sandia.gov/cubit (2001)

Edelsbrunner, H. and N. Shah, “Incremental topological flipping works for regular triangular,” Proc. 8th Symposium on Computational Geometry, pp. 43-52 (1992)

Field, D. “Laplacian smoothing and Delaunay triangulations,” Communications in Numerical Methods in Engineering, 4, 709-712, (1988)

Freitag, Lori A. and Patrick M. Knupp, “Tetrahedral element shape optimization via the Jacobian determinant and the condition number,” Proc. 8th International Meshing Roundtable, pp. 247-258 (1999)

Freitag, Lori, “On Combining Laplacian and optimization-based smoothing techniques,” AMD Trends in Unstructured Mesh Genertation, ASME, 220, 37-43, July 1997

Garmilla, Rao V. and Mark S. Shephard, “Boundary layer mesh generation for viscous flow simulations,” International Journal for Numerical Methods in Engineering, 49, 193-218 (2000)

George, P. L., F. Hecht and E. Saltel, "Automatic Mesh Generator with Specified Boundary", Computer Methods in Applied Mechanics and Engineering, 92, 269-288 (1991)

George, Paul-Louis and Houman Borouchaki, “Delaunay Triangulation and Meshing”, Hermes, Paris, 1998

Joe, B. Geompack90 software URL: http://tor-pw1.attcanada.ca/~bjoe/index.htm (1999)

Joe, B., "GEOMPACK - A Software Package for the Generation of Meshes Using Geometric Algorithms", Advances in Engineering Software, 56(13), 325-331 (1991)

Joe, B., “Construction of improved quality triangulations using local transformations,” SIAM Journal on Scientific Computing, 16, 1292-1307 (1995)

Lau, T.S. and S.H. Lo, “Finite Element Mesh Generation Over Analytical Surfaces”, Computers and Structures, 59(2), 301-309 (1996)

Lau, T.S., S. H. Lo and C. K. Lee, “Generation of Quadrilateral Mesh over Analytical Curved Surfaces,” Finite Elements in Analysis and Design, 27, 251-272 (1997)

References

Mesh Generation Tutorial 217Instituto Tecnológico de Veracruz 21-25 April 2008

Khawaja, Aly and Yannis Kallinderis, “Hybrid grid generation for turbomachinery and aerospace applications,” International Journal for Numerical Methods in Engineering, 49, 145-166 (2000)

Knupp, P. “Algebraic Mesh Quality Metrics for Unstructured Initial Meshes,” SAND 2001-1448J, Sandia National Laboratories, submitted 2001

Lawson, C. L., "Software for C1 Surface Interpolation", Mathematical Software III, 161-194 (1977)

Lee, C.K, and S.H. Lo, “A New Scheme for the Generation of a Graded Quadrilateral Mesh,” Computers and Structures, 52, 847-857 (1994)

Liu, Anwei, and Barry Joe, “Relationship between tetrahedron shape measures,” BIT, 34, 268-287 (1994)

Liu, Shang-Sheng and Rajit Gadh, "Basic LOgical Bulk Shapes (BLOBS) for Finite Element Hexahedral Mesh Generation", 5th International Meshing Roundtable, 291-306 (1996)

Lo, S. H., “Volume Discretization into Tetrahedra - II. 3D Triangulation by Advancing Front Approach”, Computers and Structures, 39(5), 501-511 (1991)

Lo, S. H., “Volume Discretization into Tetrahedra-I. Verification and Orientation of Boundary Surfaces”, Computers and Structures, 39(5), 493-500 (1991)

Lohner, R. K. Morgan and O.C. Zienkiewicz, “Adaptive grid refinement for thr compressible Euler equations, “in I. Babuska, O. C. Zienkiewicz, J. Gago, and E.R. Oliviera (Eds.) Accuracy estimates and adaptive Refinements in Finite Element Computations, Wiley, p. 281-297 (1986)

Lohner, R., “Progress in Grid Generation via the Advancing Front Technique”, Engineering with Computers, 12, 186-210 (1996)

Lohner, Rainald, Paresh Parikh and Clyde Gumbert, "Interactive Generation of Unstructured Grid for Three Dimensional Problems", Numerical Grid Generation in Computational Fluid Mechanics `88, Pineridge Press, 687-697 (1988)

Lohner, Rainald and Juan Cebral, “Generation of non-isotropic unstructured grids via directional enrichment,” International Journal for Numerical Methods in Engineering, 49, 219-232 (2000)

Marcum, David L. and Nigel P. Weatherill, “Unstructured Grid Generation Using Iterative Point Insertion and Local Reconnection”, AIAA Journal, 33(9),1619-1625 (1995)

Marcum, David L. and J. Adam Gaither “Unstructured Surface Grid Generation using Global Mapping and Physical Space Approximation”, Proceedings 8th International Meshing Roundtable, 397-406 (1999)

Mitchell, Scott A., "High Fidelity Interval Assignment", Proceedings, 6th International Meshing Roundtable, 33-44 (1997)

References

Mesh Generation Tutorial 218Instituto Tecnológico de Veracruz 21-25 April 2008

Murdoch, Peter, and Steven E. Benzley, “The Spatial Twist Continuum”, Proceedings, 4th International Meshing Roundtable, 243-251 (1995)

Owen, Steve J. and Sunil Saigal, "H-Morph: An Indirect Approach to Advancing Front Hex Meshing". International Journal for Numerical Methods in Engineering, Vol 49, No 1-2, (2000) pp. 289-312

Owen, Steven J. "Constrained Triangulation: Application to Hex-Dominant Mesh Generation", Proceedings, 8th International Meshing Roundtable, 31-41 (1999)

Owen, Steven J. and Sunil Saigal, "Formation of Pyramid Elements for Hex to Tet Transitions", Computer Methods in Applied Mechanics in Engineering, Accepted for publication (approx. November 2000)

Owen, Steven J. and Sunil Saigal, "Surface Mesh Sizing Control", International Journal for Numerical Methods in Engineering, Vol. 47, No. 1-3, (2000) pp.497-511.

Owen, Steven J., M. L. Staten, S. A. Canann, and S. Saigal "Q-Morph: An Indirect Approach to Advancing Front Quad Meshing", International Journal for Numerical Methods in Engineering, Vol 44, No. 9, (1999) pp. 1317-1340

Pirzadeh S. “Viscous unstructured grids by the advancing-layers method”, In proceedings 32nd Aerospace Sciences Meeting and Exhibit, AIAA-94-0417, Reno, NV 1994.

Price, M.A. and C.G. Armstrong, “Hexahedral Mesh Generation by Medial Surface Subdivision: Part I”, International Journal for Numerical Methods in Engineering. 38(19), 3335-3359 (1995)

Price, M.A. and C.G. Armstrong, “Hexahedral Mesh Generation by Medial Surface Subdivision: Part II,” International Journal for Numerical Methods in Engineering, 40, 111-136 (1997)

Rebay, S., “Efficient Unstructured Mesh Generation by Means of Delaunay Triangulation and Bowyer-Watson Algorithm”, Journal Of Computational Physics, 106,125-138 (1993)

Ruppert, Jim, “A New and Simple Algorithm for Quality 2-Dimensional Mesh Generation”, Technical Report UCB/CSD 92/694, University of California at Berkely, Berkely California (1992)

Salem, Ahmed Z. I., Scott A. Canann and Sunil Saigal, “Mid-node admissible spaces for quadratic triangular arbitrarily curved 2D finite elements,” International Journal for Numerical Methods in Engineering, 50, 253-272 (2001)

Schneiders, Robert, “A Grid-Based Algorithm for the Generation of Hexahedral Element Meshes”, Engineering With Computers, 12, 168-177 (1996)

Shephard, Mark S. and Marcel K. Georges, “Three-Dimensional Mesh Generation by Finite Octree Technique”, International Journal for Numerical Methods in Engineering, 32, 709-749 (1991)

Sheffer, Alla and E. de Sturler, “Surface Parameterization for Meshing by Triangulation Flattening”, Proceedings, 9th International Meshing Roundtable (2000)

References

Mesh Generation Tutorial 219Instituto Tecnológico de Veracruz 21-25 April 2008

Sheffer, Alla and E. de Sturler, “Surface Parameterization for Meshing by Triangulation Flattening”, Proceedings, 9th International Meshing Roundtable (2000)

Sheffer, Alla, Michal Etzion, Ari Rappoport, Michal Bercovier, “Proceedings 7th International Meshing Roundtable,” pp.347-364 (1998)

Shewchuk, Jonathan Richard, “Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator”, URL: http://www.cs.cmu.edu/~quake/triangle.html (1996)

Shimada K., “Physically-based Mesh Generation: Automated triangulation of surfaces and volumes via bubble packing,” Ph.D. Thesis, MEDept., MIT, Cambridge , MA (1993)

Shimada, K., “Anisotropic triangular meshing of parametric surfaces via close packing ofellipsoidal bubbles,” Proc. 6th International Meshing Roundtable, pp. 375-390(1997)

Staten, Matthew L. and Scott A. Canann, “Post Refinement Element Shape Improvement for Quadrilateral Meshes,” AMD-220, Trends in Unstructured Mesh Generation, ASME 1997 pp.9-16

Staten, M. L., S. A. Canann and S. J. Owen, “BMSweep: Locating interior nodes during sweeping”, Engineering With Computers, vol 5, No. 3. (1999) pp.212-218

Staten, Mathew. L., Steven J. Owen, Ted D. Blacker, “Unconstrained Plastering”, Proceedings 14th International meshing Roundtable (2005)

Talbert, J.A., and A.R. Parkinson, “Development of an Automatic, Two Dimensional Finite Element Mesh Generator using Quadrilateral Elements and Bezier Curve Boundary Definitions”, International Journal for Numerical Methods in Engineering, 29, 1551-1567 (1991)

Tam, T. K. H. and C. G. Armstrong. “2D Finite Element Mesh Generation by Medial Axis Subdivision”, Advances in Engineering Software, 13, 313-324 (1991)

Taniguchi, Takeo, Tomoaki Goda, Harald Kasper and Werner Zielke, "Hexahedral Mesh Generation of Complex Composite Domain", 5th International Conference on Grid Generation in Computational Field Simmulations, Mississippi State University. 699-707 (1996)

Tautges, Timothy J., Shang-sheng Liu, Yong Lu, Jason Kraftcheck and Rajit Gadh, "Feature Recognition Applications in Mesh Generation", Trends in Unstructured Mesh Generation, ASME, AMD-220, 117-121 (1997)

Tautges, Timothy J., Ted Blacker and Scott Mitchell, “The Whisker-Weaving Algorithm: A Connectivity Based Method for Constructing All-Hexahedral Finite Element Meshes,” International Journal for Numerical Methods in Engineering, 39, 3327-3349 (1996)

Tautges, Timothy J., “The Common Geometry Module (CGM): A Generic Extensible Geometry Interface,” Proceedings, 9th International Meshing Roundtable, pp. 337-347 (2000)

Thompson, Joe F, Bharat K. Soni, Nigel P. Weatherill Eds. Handbook of Grid Generation, CRC Press, 1999

References

Mesh Generation Tutorial 220Instituto Tecnológico de Veracruz 21-25 April 2008

Thompson, Joe F, Bharat K. Soni, Nigel P. Weatherill Eds. Handbook of Grid Generation, CRC Press, 1999

Watson, David F., “Computing the Delaunay Tesselation with Application to Voronoi Polytopes”, The Computer Journal, 24(2) 167-172 (1981)

Weatherill, N. P. and O. Hassan, “Efficient Three-dimensional Delaunay Triangulation with Automatic Point Creation and Imposed Boundary Constraints”, International Journal for Numerical Methods in Engineering, 37, 2005-2039 (1994)

White, David R. and Paul Kinney, “Redesign of the Paving Algorithm: Robustness Enhancements through Element by Element Meshing,” Proceedings, 6th International Meshing Roundtable, 323-335 (1997)

White, David R.,. “Automated Hexahedral Mesh Generation by Virtual Decomposition”, Proceedings, 4th International Meshing Roundtable, pp.165-176 (1995)

Yerry, Mark A. and Mark S, Shephard, “Three-Dimensional Mesh Generation by Modified Octree Technique”, International Journal for Numerical Methods in Engineering, 20, 1965-1990 (1984)

top related