The need to follow fronts moving with curvature-dependent speed arises in the modeling of a wide class of physical phenomena, such as crystal growth, flame propagation and secondary oil recovery. In this paper, we show how to design numerical algorithms to follow a closed, non-intersecting hypersurface propagating along its normal vector field with curvature-dependent speed. The essential idea is an Eulerian formulation of the equations of motion into a Hamilton- Jacobi equation with parabolic right-hand side. This is in contrast to marker particle methods, which are rely on Lagrangian discretizations of a moving parameterized front, and suffer from instabilities, excessively small time step requirements, and difficulty in handling topological changes in the propagating front. In our new Eulerian setting, the numerical algorithms for conservation laws of hyperbolic systems may be used to solve for the propagating front In this form, the entropy-satisfying algorithms naturally handle singularities in the propagating front, as well as complicated topological changes such as merging and breaking. We demonstrate the versatility of these new algorithms by computing the solutions of a wide variety of surface motion problems in two and three dimensions showing sharpening, breaking and merging.
[1]
S. Osher,et al.
Uniformly high order accurate essentially non-oscillatory schemes, 111
,
1987
.
[2]
M. Grayson.
The heat equation shrinks embedded plane curves to round points
,
1987
.
[3]
P. Lions,et al.
Two approximations of solutions of Hamilton-Jacobi equations
,
1984
.
[4]
P. Concus,et al.
Variational Methods for Free Surface Interfaces
,
1987
.
[5]
J. Sethian.
AN ANALYSIS OF FLAME PROPAGATION
,
1982
.
[6]
J. Sethian,et al.
Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations
,
1988
.
[7]
S. Osher.
Riemann Solvers, the Entropy Condition, and Difference
,
1984
.
[8]
J. Sethian.
Curvature and the evolution of fronts
,
1985
.