Level Set Method in a Finite Element Setting
|
|
- Stuart Stephens
- 6 years ago
- Views:
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
The Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationWETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES
November 20, 2018 WETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES Bishal Bhattarai and Nikolai V. Priezjev Department of Mechanical and Materials Engineering
More informationTheoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization. M2DO Lab 1,2. 1 Cardiff University
Theoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization M2DO Lab 1,2 1 Cardiff University 2 University of California, San Diego November 2017 A brief description of theory
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationMath 690N - Final Report
Math 690N - Final Report Yuanhong Li May 05, 008 Accurate tracking of a discontinuous, thin and evolving turbulent flame front has been a challenging subject in modelling a premixed turbulent combustion.
More informationSimulation Details for 2D
Appendix B Simulation Details for 2D In this appendix we add some details two-dimensional simulation method. The details provided here describe the method used to obtain results reported in Chapters 3
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More informationCS 231. Fluid simulation
CS 231 Fluid simulation Why Simulate Fluids? Feature film special effects Computer games Medicine (e.g. blood flow in heart) Because it s fun Fluid Simulation Called Computational Fluid Dynamics (CFD)
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationShape Modeling. Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell. CS 523: Computer Graphics, Spring 2011
CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell 2/15/2011 1 Motivation Geometry processing: understand geometric characteristics,
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationModeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures
University of California, Davis Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures N. Sukumar UC Davis Rutgers University,
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationA Grid Based Particle Method for Evolution of Open Curves and Surfaces
A Grid Based Particle Method for Evolution of Open Curves and Surfaces Shingyu Leung a,, Hongkai Zhao b a Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Hong
More informationInverse and Implicit functions
CHAPTER 3 Inverse and Implicit functions. Inverse Functions and Coordinate Changes Let U R d be a domain. Theorem. (Inverse function theorem). If ϕ : U R d is differentiable at a and Dϕ a is invertible,
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationOverview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationGeometrical Modeling of the Heart
Geometrical Modeling of the Heart Olivier Rousseau University of Ottawa The Project Goal: Creation of a precise geometrical model of the heart Applications: Numerical calculations Dynamic of the blood
More informationMeasuring Lengths The First Fundamental Form
Differential Geometry Lia Vas Measuring Lengths The First Fundamental Form Patching up the Coordinate Patches. Recall that a proper coordinate patch of a surface is given by parametric equations x = (x(u,
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationWater. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things:
Notes Added a 2D cross-section viewer for assignment 6 Not great, but an alternative if the full 3d viewer isn t working for you Warning about the formulas in Fedkiw, Stam, and Jensen - maybe not right
More informationCharacteristic Aspects of SPH Solutions
Characteristic Aspects of SPH Solutions for Free Surface Problems: Source and Possible Treatment of High Frequency Numerical Oscillations of Local Loads. A. Colagrossi*, D. Le Touzé & G.Colicchio* *INSEAN
More informationApplications* X3D Moving Grid Methods for Semiconductor INTRODUCTION ANDREW KUPRAT *, DAVID CARTWRIGHT, J. TINKA GAMMEL,
VLSI DESIGN 1998, Vol. 8, Nos. (1-4), pp. 117-121 Reprints available directly from the publisher Photocopying permitted by license only (C) 1998 OPA (Overseas Publishers Association) N.V. Published by
More informationSnakes, Active Contours, and Segmentation Introduction and Classical Active Contours Active Contours Without Edges
Level Sets & Snakes Snakes, Active Contours, and Segmentation Introduction and Classical Active Contours Active Contours Without Edges Scale Space and PDE methods in image analysis and processing - Arjan
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More informationPARAMETRIC SHAPE AND TOPOLOGY OPTIMIZATION WITH RADIAL BASIS FUNCTIONS
PARAMETRIC SHAPE AND TOPOLOGY OPTIMIZATION WITH RADIAL BASIS FUNCTIONS Michael Yu Wang 1 and Shengyin Wang 1 Department of Automation and Computer-Aided Engineering The Chinese University of Hong Kong
More informationDr. Ulas Bagci
Lecture 9: Deformable Models and Segmentation CAP-Computer Vision Lecture 9-Deformable Models and Segmentation Dr. Ulas Bagci bagci@ucf.edu Lecture 9: Deformable Models and Segmentation Motivation A limitation
More informationFrom CAD surface models to quality meshes. Patrick LAUG. Projet GAMMA. INRIA Rocquencourt. Outline
From CAD surface models to quality meshes Patrick LAUG Projet GAMMA INRIA Rocquencourt Tetrahedron II Oct. 2007 1 Outline 1. Introduction B-Rep, patches 2. CAD repair ant topology recovery 3. Discretization
More informationLevel Set Models for Computer Graphics
Level Set Models for Computer Graphics David E. Breen Department of Computer Science Drexel University Ross T. Whitaker School of Computing University of Utah Ken Museth Department of Science and Technology
More informationThe Level Set Method THE LEVEL SET METHOD THE LEVEL SET METHOD 203
The Level Set Method Fluid flow with moving interfaces or boundaries occur in a number of different applications, such as fluid-structure interaction, multiphase flows, and flexible membranes moving in
More informationOverview. Applications of DEC: Fluid Mechanics and Meshing. Fluid Models (I) Part I. Computational Fluids with DEC. Fluid Models (II) Fluid Models (I)
Applications of DEC: Fluid Mechanics and Meshing Mathieu Desbrun Applied Geometry Lab Overview Putting DEC to good use Fluids, fluids, fluids geometric interpretation of classical models discrete geometric
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationLevel Sets Methods in Imaging Science
Level Sets Methods in Imaging Science Dr. Corina S. Drapaca csd12@psu.edu Pennsylvania State University University Park, PA 16802, USA Level Sets Methods in Imaging Science p.1/36 Textbooks S.Osher, R.Fedkiw,
More informationDroplet collisions using a Level Set method: comparisons between simulation and experiments
Computational Methods in Multiphase Flow III 63 Droplet collisions using a Level Set method: comparisons between simulation and experiments S. Tanguy, T. Ménard & A. Berlemont CNRS-UMR6614-CORIA, Rouen
More informationDiscrete representations of geometric objects: Features, data structures and adequacy for dynamic simulation. Part I : Solid geometry
Discrete representations of geometric objects: Features, data structures and adequacy for dynamic simulation. Surfaces Part I : Solid geometry hachar Fleishman Tel Aviv University David Levin Claudio T.
More informationTechnical Report TR
Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,
More informationA splitting method for numerical simulation of free surface flows of incompressible fluids with surface tension
A splitting method for numerical simulation of free surface flows of incompressible fluids with surface tension Kirill D. Nikitin Maxim A. Olshanskii Kirill M. Terekhov Yuri V. Vassilevski Abstract The
More informationSurfaces and Integral Curves
MODULE 1: MATHEMATICAL PRELIMINARIES 16 Lecture 3 Surfaces and Integral Curves In Lecture 3, we recall some geometrical concepts that are essential for understanding the nature of solutions of partial
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
More informationA hybrid level-set method in two and three dimensions for modeling detonation and combustion problems in complex geometries
and Applied Mechanics University of Illinois at Urbana-Champaign TAM Report No. 1040 UILU-ENG-2004-6001 ISSN 0073-5264 A hybrid level-set method in two and three dimensions for modeling detonation and
More informationAdaptive numerical methods
METRO MEtallurgical TRaining On-line Adaptive numerical methods Arkadiusz Nagórka CzUT Education and Culture Introduction Common steps of finite element computations consists of preprocessing - definition
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationData Visualization. Fall 2017
Data Visualization Fall 2017 Vector Fields Vector field v: D R n D is typically 2D planar surface or 2D surface embedded in 3D n = 2 fields tangent to 2D surface n = 3 volumetric fields When visualizing
More informationThe Transient Modeling of Bubble Pinch-Off Using an ALE Moving Mesh
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston The Transient Modeling of Bubble Pinch-Off Using an ALE Moving Mesh Christopher J. Forster, Marc K. Smith * School of Mechanical Engineering,
More informationMULTIPHASE LEVEL SET EVOLUTION WITH APPLICATIONS TO AUTOMATIC GENERATIONAL TRACKING OF CELL DIVISION OF ESCHERICHIA COLI. A Thesis.
MULTIPHASE LEVEL SET EVOLUTION WITH APPLICATIONS TO AUTOMATIC GENERATIONAL TRACKING OF CELL DIVISION OF ESCHERICHIA COLI A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment
More informationMotivation: Level Sets. Input Data Noisy. Easy Case Use Marching Cubes. Intensity Varies. Non-uniform Exposure. Roger Crawfis
Level Set Motivation: Roger Crawfi Slide collected from: Fan Ding, Charle Dyer, Donald Tanguay and Roger Crawfi 4/24/2003 R. Crawfi, Ohio State Univ. 109 Eay Cae Ue Marching Cube Input Data Noiy 4/24/2003
More informationPrinciples and Goals
Simulations of foam rheology Simon Cox Principles and Goals Use knowledge of foam structure to predict response when experiments are not easy to isolate causes of certain effects to save time in designing
More information2.11 Particle Systems
2.11 Particle Systems 320491: Advanced Graphics - Chapter 2 152 Particle Systems Lagrangian method not mesh-based set of particles to model time-dependent phenomena such as snow fire smoke 320491: Advanced
More informationHigh Resolution Sharp Computational Methods for Elliptic and Parabolic Problems in Complex Geometries
High Resolution Sharp Computational Methods for Elliptic and Parabolic Problems in Complex Geometries Frédéric Gibou Chohong Min Ron Fedkiw November 2, 2012 In honor of Stan Osher s 70 th birthday Abstract
More informationEvolution, Implementation, and Application of Level Set and Fast Marching Methods for Advancing Fronts
Evolution, Implementation, and Application of Level Set and Fast Marching Methods for Advancing Fronts J.A. Sethian Dept. of Mathematics University of California, Berkeley 94720 Feb. 20, 2000 Abstract
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More informationDevelopment of Hybrid Fluid Jet / Float Polishing Process
COMSOL Conference - Tokyo 2013 Development of Hybrid Fluid Jet / Float Polishing Process A. Beaucamp, Y. Namba Dept. of Mechanical Engineering, Chubu University, Japan Zeeko LTD, United Kingdom Research
More informationLaplacian Meshes. COS 526 Fall 2016 Slides from Olga Sorkine and Yaron Lipman
Laplacian Meshes COS 526 Fall 2016 Slides from Olga Sorkine and Yaron Lipman Outline Differential surface representation Ideas and applications Compact shape representation Mesh editing and manipulation
More informationLecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More informationThe Design and Implementation of a Modeling Package
The Design and Implementation of a Modeling Package p.1/20 The Design and Implementation of a Modeling Package Hossein ZivariPiran hzp@cs.toronto.edu Department of Computer Science University of Toronto
More informationVerification of Moving Mesh Discretizations
Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions
More informationAn Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE)
An Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE) Shingyu Leung Department of Mathematics, Hong Kong University of Science and Technology masyleung@ust.hk May, Shingyu Leung (HKUST)
More informationRobotics. Chapter 25-b. Chapter 25-b 1
Robotics Chapter 25-b Chapter 25-b 1 Particle Filtering Particle filtering uses a population of particles (each particle is a state estimate) to localize a robot s position. This is called Monte Carlo
More informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationPredicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity
Predicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity November 8th, 2006 Outline Motivation Motivation Motivation for Modelling Breast Deformation Mesh Generation
More informationAdvances in 3-D Forging Process Modeling
Advances in 3-D Forging Process Modeling W J Slagter, C J L Florie and A C J Venis MSC.Software Corporation, Groningenweg 6, 2803 PV Gouda, The Netherlands MSC/SuperForge is a new code developed for performing
More informationTopic 5-6: Parameterizing Surfaces and the Surface Elements ds and ds. Big Ideas. What We Are Doing Today... Notes. Notes. Notes
Topic 5-6: Parameterizing Surfaces and the Surface Elements ds and ds. Textbook: Section 16.6 Big Ideas A surface in R 3 is a 2-dimensional object in 3-space. Surfaces can be described using two variables.
More informationGeometric Modeling in Graphics
Geometric Modeling in Graphics Part 10: Surface reconstruction Martin Samuelčík www.sccg.sk/~samuelcik samuelcik@sccg.sk Curve, surface reconstruction Finding compact connected orientable 2-manifold surface
More informationFree Surface Lattice-Boltzmann fluid simulations. with and without level sets.
ree Surface Lattice-Boltzmann fluid simulations with and without level sets Nils Thürey, Ulrich Rüde University of Erlangen-Nuremberg System Simulation roup Cauerstr. 6, 91054 Erlangen, ermany Email: Nils.Thuerey@cs.fau.de
More informationVector Visualization. CSC 7443: Scientific Information Visualization
Vector Visualization Vector data A vector is an object with direction and length v = (v x,v y,v z ) A vector field is a field which associates a vector with each point in space The vector data is 3D representation
More informationParameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia
Parameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia 2008 1 Non-Convex Non Convex Boundary Convex boundary creates significant distortion Free boundary is better 2 Fixed
More informationCoping with the Ice Accumulation Problems on Power Transmission Lines
Coping with the Ice Accumulation Problems on Power Transmission Lines P.N. Shivakumar 1, J.F.Peters 2, R.Thulasiram 3, and S.H.Lui 1 1 Department of Mathematics 2 Department of Electrical & Computer Engineering
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationGeometric Acoustics in High-Speed Boundary Layers
Accepted for presentation at the 9th International Symposium on Shock Waves. Madison, WI. July -9,. Paper #8 Geometric Acoustics in High-Speed Boundary Layers N. J. Parziale, J. E. Shepherd, and H. G.
More information8 Special Models for Animation. Chapter 8. Special Models for Animation. Department of Computer Science and Engineering 8-1
Special Models for Animation 8-1 L-Systems 8-2 L-Systems Branching Structures Botany Display geometric substitution turtle graphics Animating plants, animals 8-3 Plant examples http://algorithmicbotany.org/papers/#abop
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationLos Alamos National Laboratory
Los Alamos National Laboratory is operated by the University of California for the United States Department of Energy under contract W-7405-ENG-36 TITLE: AUTHOR(S). X3D MOVING GRID METHODS FOR SEMICONDUCTOR
More informationAn Embedded Boundary Method with Adaptive Mesh Refinements
An Embedded Boundary Method with Adaptive Mesh Refinements Marcos Vanella and Elias Balaras 8 th World Congress on Computational Mechanics, WCCM8 5 th European Congress on Computational Methods in Applied
More informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
More informationMedical Image Segmentation using Level Sets
Medical Image Segmentation using Level Sets Technical Report #CS-8-1 Tenn Francis Chen Abstract Segmentation is a vital aspect of medical imaging. It aids in the visualization of medical data and diagnostics
More informationLevel Set Techniques for Tracking Interfaces; Fast Algorithms, Multiple Regions, Grid Generation, and Shape/Character Recognition
Level Set Techniques for Tracking Interfaces; Fast Algorithms, Multiple Regions, Grid Generation, and Shape/Character Recognition J.A. Sethian Abstract. We describe new applications of the level set approach
More informationSIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH
Fifth International Conference on CFD in the Process Industries CSIRO, Melbourne, Australia 3-5 December 006 SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH Manasa Ranjan Behera and Murali.
More informationMass-Spring Systems. Last Time?
Mass-Spring Systems Last Time? Implicit Surfaces & Marching Cubes/Tetras Collision Detection & Conservative Bounding Regions Spatial Acceleration Data Structures Octree, k-d tree, BSF tree 1 Today Particle
More informationThe Eikonal Equation
The Eikonal Equation Numerical efficiency versus computational compleity Shu-Ren Hysing III Institute of Applied Mathematics LSIII University of Dortmund Level set - methodology By embedding an interface
More informationIan Mitchell. Department of Computer Science The University of British Columbia
CPSC 542D: Level Set Methods Dynamic Implicit Surfaces and the Hamilton-Jacobi Equation or What Water Simulation, Robot Path Planning and Aircraft Collision Avoidance Have in Common Ian Mitchell Department
More informationLectures on Higher Gauge Theory III
Lectures on Higher Gauge Theory III John Baez joint with Toby Bartels, Alissa Crans, Alex Hoffnung, Aaron Lauda, Chris Rogers, Urs Schreiber, Danny Stevenson and Konrad Waldorf Courant Research Center
More informationMarker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 1, pp. 1-41 Commun. Comput. Phys. January 2009 REVIEW ARTICLE Marker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
More informationCOMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS
COMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS Tayfun Tezduyar tezduyar@rice.edu Team for Advanced Flow Simulation and Modeling (T*AFSM) Mechanical Engineering and Materials Science Rice University
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationImagery for 3D geometry design: application to fluid flows.
Imagery for 3D geometry design: application to fluid flows. C. Galusinski, C. Nguyen IMATH, Université du Sud Toulon Var, Supported by ANR Carpeinter May 14, 2010 Toolbox Ginzburg-Landau. Skeleton 3D extension
More informationMOTION OF A LINE SEGMENT WHOSE ENDPOINT PATHS HAVE EQUAL ARC LENGTH. Anton GFRERRER 1 1 University of Technology, Graz, Austria
MOTION OF A LINE SEGMENT WHOSE ENDPOINT PATHS HAVE EQUAL ARC LENGTH Anton GFRERRER 1 1 University of Technology, Graz, Austria Abstract. The following geometric problem originating from an engineering
More informationCurve Sampling and Geometric Conditional Simulation
Curve Sampling and Geometric Conditional Simulation Ayres Fan Joint work with John Fisher, William Wells, Jonathan Kane, and Alan Willsky S S G Massachusetts Institute of Technology September 19, 2007
More informationAs a consequence of the operation, there are new incidences between edges and triangles that did not exist in K; see Figure II.9.
II.4 Surface Simplification 37 II.4 Surface Simplification In applications it is often necessary to simplify the data or its representation. One reason is measurement noise, which we would like to eliminate,
More informationSIMULATION OF A DETONATION CHAMBER TEST CASE
SIMULATION OF A DETONATION CHAMBER TEST CASE Daniel Hilding Engineering Research Nordic AB Garnisonen I4, Byggnad 5 SE-582 10 Linköping www.erab.se daniel.hilding@erab.se Abstract The purpose of a detonation
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationModeling Human Embryos Using a Variational Level Set Approach
Modeling Human Embryos Using a Variational Level Set Approach Uffe Damgaard Pedersen June 1, 2004 Abstract At fertility clinics, living human embryos are evaluated by visual inspection of light microscopy
More information