A variational level set methodology without reinitialization for the prediction of equilibrium interfaces over arbitrary solid surfaces

Abstract

A B S T R A C TA robust numerical methodology to predict equilibrium interfaces over arbitrary solid surfaces is developed. The kernel of the proposed method is the distance regularized level set equations (DRLSE) with techniques to incorporate the no-penetration and mass-conservation constraints. In this framework, we avoid reinitialization typically used in traditional level set methods. This allows for a more efficient algorithm since only one advection equation is solved, and avoids numerical error associated with the re-distancing step. A novel surface tension distribution, based on harmonic mean, is prescribed such that the zero level set has the correct the liquid-solid surface tension value. This leads to a more accurate triple contact point location. The method uses second-order central difference schemes which facilitates easy parallel implementation, and is validated by comparing to traditional level set methods for canonical problems. The application of the method, in the context of Gibbs free energy minimization, to obtain liquid-air interfaces is validated against existing analytical solutions. The capability of our current methodology to predict equilibrium shapes over both structured and realistic rough surfaces is demonstrated.

Cooperative drag reduction in turbulent flows using polymer additives and superhydrophobic walls

“…In contrast, the average slip length on the superhydrophobic random texture is appreciably diminished (from b ≈ 25 to 18 μm) in the butanol solution. This may either be a consequence of adverse Marangoni stresses, arising from surfactant gradients that form adjacent to flow-obstructing texture asperities, or may also be due to the lowering and further impalement of the air-water interface into the random texture, caused by the reduction in the surface tension of the liquid phase [42].…”

A variational volume-of-fluid approach for front propagation

“…The simulation is stopped after 6000 iterations (t f inal = 300.0) and snapshots of the interface evolution are shown in Fig 13 for t = 0.0, 37.5, 75.0, 112.5, 150.0, 187.5, 225.0, 262.5, and 300.0. Problem parameters describing spiral geometry were chosen to match [13] for comparison. The spiral unwinds under curvature-driven motion, and its size decreases in the process until it eventually acquires a bean shape that later collapses rapidly.…”

“…This is more commonly known as the distance regularized level set evolution (DRLSE) where the regularization term penalizes the deviation of the level-set from an SDF, therefore, bypassing the need for an explicit reinitialization step. More recently, Alamé et al [13] devised a robust numerical methodology to predict equilibrium interfaces by utilizing a DRLSE approach. Their variational LS formulation comes in the context of Gibbs free energy minimization where they demonstrate comparisons between analytical solutions of liquid-air interface equilibria and the proposed variational LS method.…”

Inverse Dirichlet weighting enables reliable training of physics informed neural networks

“…using the finite-difference method with necessary reinitialization to preserve the sign distance property of the level-set function ϕ. The spiral parameters are chosen to be np = 400, D = 2.5, a = 3, m = 2 [57]. Time integration is performed using a second-order Runge-Kutta scheme with a step size of dt = 0.001.…”