A high order continuous Galerkin spectrally stabilized level-set approach for incompressible two-phase flows

Read full text
Authors: Nadish Saini, Ananias G. Tomboulides, Dillon R. Shaver, Paul Fischer, Elia Merzari, Aleksandr V. Obabko, Donna Post Guillen
Cited by: 0
References: 53

Abstract

This paper introduces a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation.The two-phase framework introduced herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al [1], is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The methodology is implemented in the high-order CG spectral element method based code Nek5000.

Fields of Study & Subjects

Lattice Boltzmann Simulation Studies Computational Fluid Dynamics and Aerodynamics Fluid Dynamics and Turbulent Flows

Publication details

Journal / source
Journal of Computational Physics
Publisher
Elsevier BV
Volume
561
Pages
114961–114961
Publication date
24 April 2026
Language
EN