Geometrical level set reinitialization using closest points method and kink detection for thin filaments, topology changes and two-phase flows

Félix Henri (1) , Mathieu Coquerelle (1) , Pierre Lubin (1)

1 I2M - Institut de Mécanique et d'Ingénierie

We introduce a robust and high order strategy to perform the reinitialization in a level set framework. The reinitialization by closest-points (RCP) method is based on geometric considerations. It relies on a gradient descent to find the closest points at the interface in order to solve the eikonal equation and thus reinitializing the level set field. Furthermore, a new algorithm, also based on a similar geometric approach, is introduced to detect precisely all the ill-defined points of the level set. These points, also referred to as kinks, can mislead the gradient descent and more widely impact the accuracy of level set methods. This algorithm coupled with the precise computation of the closest points of the interface, permits the novel method to be robust and accurate even when performing the reinitialization every time step after solving the advection equation. Furthermore, they both require very few given parameters with the advantage of being based on a geometrical approach and independent of the application. The proposed method was tested on various benchmarks, and demonstrated equivalent or even better results compared to solving the Hamilton-Jacobi equation.

Mots clés: geometrical approach; level set; reinitialization; closest points; kink detection; medial axis

Domaines: Modélisation et simulation; Mécanique des fluides [physics.class-ph]

Félix Henri, Mathieu Coquerelle, Pierre Lubin. Geometrical level set reinitialization using closest points method and kink detection for thin filaments, topology changes and two-phase flows. Journal of Computational Physics, 2022. ⟨hal-03159919v2⟩