Top Banner
5-11-2005 Lamps of Aladdin 1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah
25

5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

Dec 22, 2015

Download

Documents

Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 1

Moving Mesh Adaptation Techniques

Todd Phillips

Gary Miller

Mark Olah

Page 2: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 2

Introduction

• The Meshing Problem– Discretize a Spatial Domain– Minimize Size (Number of Triangles)– Maximize Quality (Shape of Triangles)

• Mesh Adaptation– Introduce More Elements– Base Decisions on the Function being Interpolated

• How do we make such a decision?

Page 3: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 3

The Basic Setup

• A Function Exists

Page 4: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 4

The Basic Setup

• A Mesh Also Exists

Page 5: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 5

The Basic Setup

• We Approximate the Function with the Mesh

Page 6: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 6

The Basic Setup

• We can adapt the Mesh to Better Approximate the Function

Page 7: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 7

Background: Local Feature Size

• What is Local Feature Size?– Roughly, the Size of Triangles We Should Be Using– The Smallest Distance to Two Geometric Domain Features

• Local Feature Size Should be k-Lipschitz– Doesn’t change too fast– For all p,q lfs(p) <= lfs(q) + k*dist(p,q)– This is like having derivative bounded by k.

Page 8: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 8

Optimal Meshing and Adaptive Meshing

• Optimal Meshing– Consider the LFS of the Input Domain– Output a Mesh

• Triangle Size is within a constant of LFS• Triangle Aspect Ratio is bounded

– Note the Impossabilty of this if LFS is not Lipschitz• Adaptive Meshing

– LFS Accounts for Other ‘Features’ Provided by Some Oracle

– Introduce Geometric Features to Accommodate– Perturb the Original Mesh to obtain a new Optimal Mesh

Page 9: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 9

Example of Modifying Local Feature Size

• A Sample Geometric Domain

Page 10: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 10

• LFS of the Sample Domain

Example of Modifying Local Feature Size

Page 11: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 11

• Modify the Geometric Domain to Capture a New Feature

Example of Modifying Local Feature Size

Page 12: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 12

• LFS of the Modified Domain

Example of Modifying Local Feature Size

Page 13: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 13

How can we Define an Oracle?

• Function-Angle Based Refinement– Some Function is Approximated by the Mesh– This Embeds the Mesh as a surface in one-higher

dimension– Observe the angle between faces of this surface– Introduce Features Where This Angle is Small

• Insert Circumcenters, Split Edges, etc.

Page 14: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 14

An Example in One Dimension: Function

Page 15: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 15

An Example in One Dimension: Function and Uniform Mesh

Page 16: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 16

An Example in One Dimension: Function, Mesh, and Error

Page 17: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 17

An Example in One Dimension: Function and Adaptive Mesh

Page 18: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 18

An Example in One Dimension: Adaptive Mesh LFS

Page 19: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 19

Will this work?

• Original Feature Size must be Small Enough

• Actual Function must be Differentiable

Page 20: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 20

Practicalities

• Real Mesh Adaptation Requires Coarsening As Well– Sometimes Functional Features Move or Disappear– Thrashing is a dirty word

• Real Solutions aren’t always differentiable– Give up after a Minimum Feature Size

Page 21: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 21

Why is this so great for Moving Meshes?

• Unlike the One Dimensional Movie– Fixed Mesh ( Grid points didn’t move left/right )– Constantly ‘Chasing’ the Feature

• Mesh moves with the Velocity Function– Hence, Mesh moves with the Velocity features– After features are captured, very little Mesh Adaptation

Page 22: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 22

How does this compare to other methods?

• Sensitivity Analysis and a posteriori Error Estimates– Double the Mesh Size Everywhere, compare solutions– Require multiple meshes of different sizing be maintained– Require multiple solves at every iteration– Don’t take advantage of Moving Mesh

• (Persistance of Features)

Page 23: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 23

Some 2-D Examples

Page 24: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 24

Some 2-D Examples

Page 25: 5-11-2005Lamps of Aladdin1 Moving Mesh Adaptation Techniques Todd Phillips Gary Miller Mark Olah.

5-11-2005 Lamps of Aladdin 25

Conclusion

• Summary– Use an oracle based on Function-Angles– Introduce new geometric features– Use existing geometrically driven meshing algorithms

• Future– Tighter Coupling of Oracle with Meshing Algorithms– Adaptive Time Refinement