Some numerical aspects of the level set method

    Mathias Moog Affiliation
    ; Rainer Keck Affiliation
    ; Aivars Zemitis Affiliation


Many practical applications imply the solution of free boundary value problems. If the free boundary is complex and can change its topology, it will be hard to solve such problems numerically. In recent years a new method has been developed, which can handle boundaries with complex geometries. This new method is called the level set method. However, the level set method also has some drawbacks, which are mainly concerning conservation of mass or numerical instabilities of the boundaries. Our aim is to analyze some aspects of the level set method on the basis of two‐phase flow in a Hele‐Shaw cell. We investigate instabilities of two‐phase flow between two parallel plates. A solution of the linearized problem is obtained analytically in order to check whether the numerical schemes compute reasonable results. The developed numerical scheme is based on finite difference approximations and the level set method. The equations of two‐phase Hele‐Shaw flow are written in a modified formulation using the one‐dimensional Dirac delta‐function. Since the level set function is not smooth enough after re‐initialization, special attention during the computation of curvature is needed. We propose a method that can solve the problems for two‐phase Hele‐Shaw flow with changing topology. The numerical solution shows good agreement with the analytical solution of the linearized problem. We describe the method below and analyze the results.

First Published Online: 14 Oct 2010

Keyword : -

How to Cite
Moog, M., Keck, R., & Zemitis, A. (1998). Some numerical aspects of the level set method. Mathematical Modelling and Analysis, 3(1), 140-151.
Published in Issue
Dec 15, 1998
Abstract Views
PDF Downloads
Creative Commons License

This work is licensed under a Creative Commons Attribution 4.0 International License.