Level Set Method in a Finite Element Setting

Transcription

1 Level Set Method in a Finite Element Setting John Shopple University of California, San Diego November 6, 2007

2 Outline 1 Level Set Method 2 Solute-Solvent Model 3 Reinitialization 4 Conclusion

3 Types of problems Geometric based motion Motion based on curvature Normal direction motion Simple constant motion Physical based motion Fluid dynamics Two-phase flow Combustion Liquid-gas interactions Biology Biomolecular surfaces Membranes Material science Solidification Epitaxial growth of thin films John Shopple Level Set Method in a Finite Element Setting 3/ 28

4 Front Tracking Methods 1 Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

5 Front Tracking Methods Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

6 Front Tracking Methods Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

7 Front Tracking Methods Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

8 Front Tracking Methods Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

9 Front Tracking Methods Place markers on interface 2 Calculate velocity for each marker 3 Propagate markers one time step 4 Check for topological changes of interfaces add/remove markers split/join interfaces if necessary 5 Go to step 2 John Shopple Level Set Method in a Finite Element Setting 4/ 28

10 Level Set Method The interface is represented as the zero level set of a function φ(x, y, t) = 0 Topological changes will occur naturally as time progresses. John Shopple Level Set Method in a Finite Element Setting 5/ 28

11 Derivation of the Level Set Evolution Equation Let ( x(t), y(t) ) parameterize a particle that stays on the moving interface. φ(x(0), y(0), t = 0) = 0 φ(x(1), y(1), t = 1) = 0 φ(x(2), y(2), t = 2) = 0 John Shopple Level Set Method in a Finite Element Setting 6/ 28

12 Derivation of the Level Set Evolution Equation 0 1 Let ( x(t), y(t) ) parameterize a particle that stays on the moving interface. φ(x(0), y(0), t = 0) = 0 φ(x(1), y(1), t = 1) = 0 φ(x(2), y(2), t = 2) = 0 John Shopple Level Set Method in a Finite Element Setting 6/ 28

13 Derivation of the Level Set Evolution Equation Let ( x(t), y(t) ) parameterize a particle that stays on the moving interface. φ(x(0), y(0), t = 0) = 0 φ(x(1), y(1), t = 1) = 0 φ(x(2), y(2), t = 2) = 0 John Shopple Level Set Method in a Finite Element Setting 6/ 28

14 Derivation of the Level Set Evolution Equation Let ( x(t), y(t) ) parameterize a particle that stays on the moving interface. φ(x(0), y(0), t = 0) = 0 φ(x(1), y(1), t = 1) = 0 φ(x(2), y(2), t = 2) = 0 dφ dt = φ t + φ x x (t) + φ y y (t) = φ t + φ (x (t), y (t) ) = φ t + φ V = 0 John Shopple Level Set Method in a Finite Element Setting 6/ 28

15 Level Set Evolution Equation For externally generated velocity fields: φ t + φ V = 0 For internally generated velocity fields: φ t + φ V n = 0 where V n = V n = V φ φ Interface = Γ(t) = {(x, y) φ(x, y, t) = 0} John Shopple Level Set Method in a Finite Element Setting 7/ 28

16 Basic Types of Interface Motion Vector field V Scalar field V n giving speed in level sets normal direction ) V n proportional to level set curvature κ = ( φ φ John Shopple Level Set Method in a Finite Element Setting 8/ 28

17 Local Method to Solve Level Set Equation Evolution equation φ t + V n φ = 0 1 Approximate φ at each vertex v i by the average of φ over a patch around v i 2 Use a Runga-Kutta type method using the values of V n at each vertex. v i John Shopple Level Set Method in a Finite Element Setting 9/ 28

18 Lennard-Jones Potential Lennard-Jones potential [ (σ ) 12 ( σ ) ] 6 U(r) = 4ε r r σ r U(r = σ) = 0 Minimum: U(r = σ σ) = ε r > σ attracts r < σ repels ε John Shopple Level Set Method in a Finite Element Setting 10/ 28

19 Solute-Solvent Model Molecules consist of atoms at locations x 1,...,x N Solvent region Ω w Solute region Ω m A solute-solvent interface Γ Ω w Ω m Γ John Shopple Level Set Method in a Finite Element Setting 11/ 28

20 Solute-Solvent Model Molecules consist of atoms at locations x 1,...,x N Solvent region Ω w Solute region Ω m A solute-solvent interface Γ Ω w Ω m Γ A free energy model G[Γ; x 1,..., x N ] = Γ γ ds + ρ w N i=1 Ω w U sw ( x x i ) dv John Shopple Level Set Method in a Finite Element Setting 11/ 28

21 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

22 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

23 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

24 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

25 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

26 Algorithm 1 Define a regular grid 2 Define atoms (with radii σ i ) 3 Construct initial surface 4 At each time step: 1 body-fit mesh to interface Γ 2 calculate L-J potential on Γ 3 calculate curvature κ 4 on Γ let V n = γ 0 κ+ρ 0 N i=1 U( x x i ) 5 extend V n to regular grid 6 evolve LS function on regular grid 5 Stop when free energy stops decreasing John Shopple Level Set Method in a Finite Element Setting 12/ 28

27 Body Fitted Meshes Refine mesh into interior and exterior regions. Add edges at zero level set Add edges to make conforming mesh John Shopple Level Set Method in a Finite Element Setting 13/ 28

28 Body Fitted Meshes Can calculate length of interface easily Explicitly separates interior and exterior May need extra refinement in some applications John Shopple Level Set Method in a Finite Element Setting 14/ 28

29 Reinitialization What, when, and why? John Shopple Level Set Method in a Finite Element Setting 15/ 28

30 Reinitialization What is reinitialization? Redefining the level set function to make it nicer without moving the location of the interface (too much). John Shopple Level Set Method in a Finite Element Setting 16/ 28

31 Reinitialization What is reinitialization? Redefining the level set function to make it nicer without moving the location of the interface (too much). Why reinitialize? Over time the level set function can become flat, causing errors in calculations of certain quantities. John Shopple Level Set Method in a Finite Element Setting 16/ 28

32 Reinitialization What is reinitialization? Redefining the level set function to make it nicer without moving the location of the interface (too much). Why reinitialize? Over time the level set function can become flat, causing errors in calculations of certain quantities. When do we reinitialize? Often enough that the level set function never becomes bad. Not so often to introduce more errors. John Shopple Level Set Method in a Finite Element Setting 16/ 28

33 Ideal Reinitialization Doesn t move the interface φ 1 (especially near the interface) Efficient Preserves key quantities near the interface (curvature, normal) Signed distance function? John Shopple Level Set Method in a Finite Element Setting 17/ 28

34 Ideal Reinitialization Doesn t move the interface φ 1 (especially near the interface) Efficient Preserves key quantities near the interface (curvature, normal) Signed distance function? John Shopple Level Set Method in a Finite Element Setting 17/ 28

35 Ideal Reinitialization Doesn t move the interface φ 1 (especially near the interface) Efficient Preserves key quantities near the interface (curvature, normal) Signed distance function? John Shopple Level Set Method in a Finite Element Setting 17/ 28

36 Ideal Reinitialization Doesn t move the interface φ 1 (especially near the interface) Efficient Preserves key quantities near the interface (curvature, normal) Signed distance function? John Shopple Level Set Method in a Finite Element Setting 17/ 28

37 Ideal Reinitialization Doesn t move the interface φ 1 (especially near the interface) Efficient Preserves key quantities near the interface (curvature, normal) Signed distance function? John Shopple Level Set Method in a Finite Element Setting 17/ 28

38 Reinitialization Examples A model starting problem John Shopple Level Set Method in a Finite Element Setting 18/ 28

39 Reinitialization, Poisson s Equation Method 1 Method 1 1 Body-fit mesh to interface 2 Solve φ = G in Ω φ = G + in Ω + φ = 0 on Γ φ n = 0 on Ω John Shopple Level Set Method in a Finite Element Setting 19/ 28

40 Reinitialization, Poisson s Equation Method φ = 1 in Ω φ = 1 in Ω φ = 0 on Γ φ n = 0 on Ω John Shopple Level Set Method in a Finite Element Setting 20/ 28

41 Reinitialization, Poisson s Equation Method 1 Black: original interface Red: new interface John Shopple Level Set Method in a Finite Element Setting 21/ 28

42 Reinitialization, Poisson s Equation Method 2 Method 2 1 Body-fit mesh to interface 2 Solve φ = G in Ω φ = 0 in Ω + φ = 0 on Γ φ = 1 on Ω John Shopple Level Set Method in a Finite Element Setting 22/ 28

43 Reinitialization, Poisson s Equation Method φ = 10 in Ω 0.1 φ = 0 in Ω φ = 0 on Γ φ = 1 on Ω John Shopple Level Set Method in a Finite Element Setting 23/ 28

44 Reinitialization, Poisson s Equation Method 2 John Shopple Level Set Method in a Finite Element Setting 24/ 28

45 Constrained Gradient Jump Across Interface Method Constraint Add additional constraint φ = 0 n e for all new body-fitted edges e, so unrefinement won t move interface. Constraint is enforced by basis test functions along interface Overdetermined system can only be solved in a least squares way Constraint can be weighted more John Shopple Level Set Method in a Finite Element Setting 25/ 28

46 Constrained Gradient Jump Across Interface Method Constraint Add additional constraint φ = 0 n e for all new body-fitted edges e, so unrefinement won t move interface. Constraint is enforced by basis test functions along interface Overdetermined system can only be solved in a least squares way Constraint can be weighted more John Shopple Level Set Method in a Finite Element Setting 25/ 28

47 Constrained Gradient Jump Across Interface Method Constraint Add additional constraint φ = 0 n e for all new body-fitted edges e, so unrefinement won t move interface. Constraint is enforced by basis test functions along interface Overdetermined system can only be solved in a least squares way Constraint can be weighted more John Shopple Level Set Method in a Finite Element Setting 25/ 28

48 Constrained Gradient Jump Across Interface Method Constraint Add additional constraint φ = 0 n e for all new body-fitted edges e, so unrefinement won t move interface. Constraint is enforced by basis test functions along interface Overdetermined system can only be solved in a least squares way Constraint can be weighted more John Shopple Level Set Method in a Finite Element Setting 25/ 28

49 Constrained Gradient Jump Comparison Black: original interface Red: new interface John Shopple Level Set Method in a Finite Element Setting 26/ 28

50 Solute-Solvent Model Examples 2 atom systems: 1.6σ, 2.0σ, 2.4σ, 2.8σ. sigma = 1.6 energy = 3.111e 005 sigma = 2.0 energy = 2.543e 005 sigma = 2.4 energy = 2.172e 005 sigma = 2.8 energy = 3.283e 005 John Shopple Level Set Method in a Finite Element Setting 27/ 28

51 Conclusion Future work: Implement 3 dimensional algorithm Speed up algorithm (various ways) Add more terms to free energy functional John Shopple Level Set Method in a Finite Element Setting 28/ 28