Optimization-based Level-Set Re-initialization: A Robust Interface Preserving Approach in Multiphase Problems

A. Hashemi,M.R. Hashemi,P. Ryzhakov,R. Rossi

2023-09-21

Abstract:In spite of its overall efficiency and robustness for capturing the interface in multiphase fluid dynamics simulations, the well-known shortcoming of the level-set method is associated with the lack of a systematic approach for preserving the regularity of the distance function. This is mainly due to the stretching (or compressing) effect of the strain rate especially in the vicinity of the liquid-gas interface. Level-set re-initialization is an effective treatment for this issue. However, the traditional approach based on the hyperbolic Hamilton-Jacobi equation is a computationally expensive procedure. Crucially, due to the hyperbolic nature of the formulation, the accuracy of the results hinges significantly on the specialized handling of blind spots near the liquid-gas interface intersecting the substrate. The present work proposes a two-step elliptic level-set re-initialization approach that strictly preserves the location of zero level-set via incorporation of an element splitting process. The primary initialization step helps remove any non-smoothness in the to-be regularized level-set function dramatically improving the efficiency of the secondary optimization step. Geometric representation of the boundary conditions is utilized in the initialization step, while the optimization step minimizes the reliance of the results on the treatment of the blind spots. The performance of the proposed method is examined for free and sessile three-dimensional droplets.

Fluid Dynamics,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the instability issue in interface capturing of the Level - Set Method in multiphase fluid dynamics simulations. Specifically, the distance property of the level - set function is easily affected by the strain rate near the fluid interface and becomes irregular, which will lead to a decline in the accuracy of calculation results. Although the traditional level - set re - initialization method can improve this problem, it has a high computational cost and is not effective in dealing with certain complex situations (such as blind spots near the contact line between a droplet and a solid surface).

To solve the above problems, this paper proposes an optimized elliptic - type level - set re - initialization method. This method is achieved through a two - step strategy:

1. **Elliptic - type Initial Reconstruction**: First, use the elliptic equation to preliminarily reconstruct the level - set function to remove non - smoothness and ensure the efficiency of subsequent optimization steps.
2. **Optimization Step**: Further optimize the level - set function by minimizing an objective function, and introduce a penalty term to keep the position of the zero - level set unchanged, thereby accurately maintaining the interface configuration.

In addition, this method also includes geometric treatment of boundary conditions. Especially when dealing with blind spots near the contact line, the normal gradient is approximately calculated by assuming a constant curvature. These improvements enable the new method to effectively handle complex situations that are difficult for traditional methods to deal with while maintaining the zero - level - set configuration.

The effectiveness of this method has been verified through a series of static and dynamic droplet simulations. The results show that this method performs excellently in maintaining the unit - norm of the level - set function gradient, accurately maintaining the interface and the contact line, and has high computational efficiency and robustness.