Level-set method

If the curve ${\displaystyle \Gamma }$ moves in the normal direction with a speed ${\displaystyle v}$, then the level-set function ${\displaystyle \varphi }$ satisfies the level-set equation

${\displaystyle {\frac {\partial \varphi }{\partial t}}=v|\nabla \varphi |.}$

Here, ${\displaystyle |\cdot |}$ is the Euclidean norm (denoted customarily by single bars in PDEs), and ${\displaystyle t}$ is time. This is a partial differential equation, in particular a Hamilton–Jacobi equation, and can be solved numerically, for example, by using finite differences on a Cartesian grid.[2][3]

The numerical solution of the level-set equation, however, requires sophisticated techniques. Simple finite-difference methods fail quickly. Upwinding methods, such as the Godunov method, fare better; however, the level-set method does not guarantee the conservation of the volume and the shape of the level set in an advection field that does conserve the shape and size, for example, uniform or rotational velocity field. Instead, the shape of the level set may get severely distorted, and the level set may vanish over several time steps. For this reason, high-order finite-difference schemes are generally required, such as high-order essentially non-oscillatory (ENO) schemes, and even then the feasibility of long-time simulations is questionable. Further sophisticated methods to deal with this difficulty have been developed, e.g., combinations of the level-set method with tracing marker particles advected by the velocity field.[4]

Consider a unit circle in ${\textstyle \mathbb {R} ^{2}}$, shrinking in on itself at a constant rate, i.e. each point on the boundary of the circle moves along its inwards pointing normal at some fixed speed. The circle will shrink and eventually collapse down to a point. If an initial distance field is constructed (i.e. a function whose value is the signed euclidean distance to the boundary, positive interior, negative exterior) on the initial circle, the normalised gradient of this field will be the circle normal.

If the field has a constant value subtracted from it in time, the zero level (which was the initial boundary) of the new fields will also be circular and will similarly collapse to a point. This is due to this being effectively the temporal integration of the Eikonal equation with a fixed front velocity.

In combustion, this method is used to describe the instantaneous flame surface, known as the G equation.

The level-set method was developed in 1979 by Alain Dervieux,[5] and subsequently popularized by Stanley Osher and James Sethian. It has become popular in many disciplines, such as image processing, computer graphics, computational geometry, optimization, computational fluid dynamics, and computational biology.

A number of level-set data structures have been developed to facilitate the use of the level-set method in computer applications.

• Computational fluid dynamics
• Combustion
• Trajectory planning
• Optimization
• Image processing
• Computational biophysics
• Discrete complex dynamics: visulalisation of parameter plane and dynamic plane

To run a Math Model in the interface of two different fluids we need to soften the interactions between the fluids. Therefore we need to apply a specific function: Compact Level Set Method.

As a “spin off”, the CompactLSM is a complement of the LSM, that helps solving LSM equations. It can be used in numerical simulation of flow, for example, if we are working with discretization of the interface water-air, compacts at sixth order, ensures the accurate and fast calculation of the interface equations (Monteiro 2018).

The LSM uses a distance function to locate different fluids. A distance function is that whose value represents the smallest distance from the point where it is being analyzed to the interface. This distance function is identified by isolines (2D) or isosurfaces (3D), showing that  the negative values refer to one of the fluids, positive values refer to the other and the zero value corresponds to the position of the interface.

But, how is the Heaviside function inserted in the Compact Level Set Method?

Since the specific mass and viscosity are discontinuous at the interface, both excess diffusion problem (interface widening) and numerical oscillations are expected if there is no adequate treatment of the fluid near the interface. To minimize these problems, the Level Set method uses a smooth, cell-related Heaviside function that explicitly defines the interface position (${\displaystyle \varphi =0}$).

The transition in the interface is kept smooth, but with a thickness of the order of magnitude of the cell size, to avoid the introduction of disturbances with a length scale equal to that of the mesh, since the interface infers an abrupt jump property from one cell to the next (Unverdi and Tryggvason, 1992). To reconstruct the material properties of the flow, such as specific mass and viscosity, another marker function, ${\displaystyle I(\varphi )}$, of the Heaviside type is used:

${\displaystyle I(\varphi )={\begin{cases}0,&{\text{if }}\varphi <-\delta \Delta \\[8pt]{\dfrac {1}{2}}\left[1+{\dfrac {\varphi }{\delta \Delta }}+{\dfrac {1}{\pi }}\sin \left({\dfrac {\pi \varphi }{\delta \Delta }}\right)\right],&{\text{if }}|\varphi |\leq \delta \Delta \\[8pt]1,&{\text{if }}\varphi >\delta \Delta \end{cases}}}$

(1)

where ${\displaystyle \delta }$ is an empirical coefficient, usually equal to 1; 5 and ${\displaystyle \Delta }$ is the characteristic discretization of the problem, which varies according to the phenomenon to be simulated. The value of ${\displaystyle \delta }$ represents an interface with a thickness of three cells, and thus ${\displaystyle \delta \Delta }$ represents half the thickness of the interface. Note that in this method, the interface has a virtual thickness, as it is represented by a smooth function. Physical properties, such as specific mass and kinematic viscosity, are calculated as:

${\displaystyle \rho =(1-I)\rho _{1}+I\rho _{2}\qquad e\qquad v=(1-I)v_{1}+Iv_{2}}$

(2)

where ${\displaystyle \rho _{1}}$, ${\displaystyle \rho _{2}}$, ${\displaystyle v_{1}}$ and ${\displaystyle v_{2}}$ are the specific mass and kinematic viscosity of fluids 1 and 2. Equation 2 can be applied analogously to the other properties of the fluids.

• See Ronald Fedkiw's academic web page for many stunning pictures and animations showing how the level-set method can be used to model real-life phenomena, like fire, water, cloth, fracturing materials, etc.
• Multivac is a C++ library for front tracking in 2D with level-set methods.
• James Sethian's web page on level-set method.
• Stanley Osher's homepage.
• The Level Set Method. MIT 16.920J / 2.097J / 6.339J. Numerical Methods for Partial Differential Equations by Per-Olof Persson. March 8, 2005