Extending fields in a level set method by solving a biharmonic equation

Abstract

We present an approach for computing extensions of velocities or other fields in level set methods by solving a biharmonic equation. The approach differs from other commonly used approaches to velocity extension because it deals with the interface fully implicitly through the level set function. No explicit properties of the interface, such as its location or the velocity on the interface, are required in computing the extension. These features lead to a particularly simple implementation using either a sparse direct solver or a matrix-free conjugate gradient solver. Furthermore, we propose a fast Poisson preconditioner that can be used to accelerate the convergence of the latter. We demonstrate the biharmonic extension on a number of test problems that serve to illustrate its effectiveness at producing smooth and accurate extensions near interfaces. A further feature of the method is the natural way in which it deals with symmetry and periodicity, ensuring through its construction that the extension field also respects these symmetries.

Accepted by Journal of Computational Physics

Cited by in corpus (7)

• Numerical investigation of controlling interfacial instabilities in non-standard Hele-Shaw configurations
• Compression-driven viscous fingering in a radial Hele-Shaw cell
• A review of one-phase Hele-Shaw flows and a level-set method for non-standard configurations
• Viscous fingering patterns for Hele--Shaw flow in a doubly connected geometry driven by a pressure differential or rotation
• Moving boundary problems for quasi-steady conduction limited melting
• PDE-Based Multidimensional Extrapolation of Scalar Fields over Interfaces with Kinks and High Curvatures
• Interfacial dynamics and pinch-off singularities for axially symmetric Darcy flow