Stabilized formulation for phase-transforming flows with special emphasis on cavitation inception

Despite significant advancement in the field of cavitation, the inception of vapor bubbles from pure liquid is still not well understood. In this paper, we propose a numerical algorithm to solve the Navier-Stokes-Korteweg equations, which are phasefield equations that govern the dynamics of multiphase compressible flows. These equations can predict cavitation inception without phenomenological assumptions about mass transfer. We present a modification to the bulk free energy of the fluid that can help us simulate cavitation at much larger length scales than previously possible. We propose a numerical scheme based on a Taylor-Galerkin discretization, a residual-based discontinuity capturing operator and a spatial filtering scheme. Our numerical method can be used for high speed flows with large pressure gradients, which are distinctive features of flows where we observe cavitation inception. We present numerical examples that show the accuracy, stability and robustness of our scheme. Finally, we use the proposed numerical scheme to study benchmark problems in the field of cavitation inception such as flow past a wedge and flow past a bluff body. The proposed numerical scheme opens opportunities to understand different aspects of cavitation in general and inception in particular, in much more detail. & COPY; 2023 Elsevier B.V. All rights reserved.

Automated summary

In this paper, a numerical algorithm is proposed to solve the Navier-Stokes-Korteweg equations, which can predict cavitation inception without assumptions about mass transfer. The proposed numerical scheme uses a modified bulk free energy and a Taylor-Galerkin discretization. It is capable of simulating high-speed flows with large pressure gradients, and is demonstrated to be accurate, stable, and robust.