Abstract:We extend the unstructured LEvel set / froNT tracking (LENT) method for handling two-phase flows with strongly different densities (high-density ratios) by providing the theoretical basis for the numerical consistency between the mass and momentum conservation in the collocated Finite Volume discretization of the single-field two-phase Navier-Stokes equations. Our analysis provides the theoretical basis for the mass conservation equation introduced by Ghods and Herrmann [3] and used in [4, 5, 6, 7, 8]. We use a mass flux that is consistent with mass conservation in the implicit Finite Volume discretization of the two-phase momentum convection term, and solve the single-field Navier-Stokes equations with our SAAMPLE segregated solution algorithm [2]. The proposed $\rho$LENT method recovers exact numerical stability for the two-phase momentum advection of a spherical droplet with density ratios ranging in $[1, 10^4]$. Numerical stability is demonstrated for in terms of the relative $L_\infty$ velocity error norm, for density-ratios in the range of $[1, 10^4]$, dynamic viscosity-ratios in the range of $[1, 10^4]$ and very strong surface tension forces, for challenging mercury/air and water/air fluid pairings. In addition, the solver performs well in cases characterized by strong interaction between two phases, i.e., oscillating droplets and rising bubbles. The proposed $\rho$LENT method is applicable to any other two-phase flow simulation method that discretizes the single-field two-phase Navier-Stokes Equations using the collocated unstructured Finite Volume Method but does not solve an advection equation for the phase indicator using a flux-based approach, by adding the proposed geometrical approximation of the mass flux and the auxiliary mass conservation equation to the solution algorithm.
What problem does this paper attempt to address?
This paper aims to solve the serious challenges encountered in the numerical simulation of two - phase fluids under the condition of high - density ratio (that is, the density difference between the two fluids is extremely large). Specifically, the paper focuses on how to maintain the consistency of mass conservation and momentum conservation numerically, especially when using the unstructured finite volume method (Finite Volume Method, FVM) for two - phase fluid simulation. High - density ratios can lead to a deterioration in the condition number of the pressure Poisson equation and the accumulation of numerical errors due to inconsistent mass and momentum transport methods. These problems are particularly prominent when dealing with gas - liquid two - phase flows in natural and industrial processes such as fuel injection, water tank sloshing, mold filling, and water - flooding - oil.
To address these challenges, the paper proposes an extended unstructured LENT (Level Set / Front Tracking) method, called the ρLENT method. This method solves the problem by providing the theoretical basis for numerical consistency between mass and momentum conservation in the collocated finite - volume discretization of the single - field two - phase Navier - Stokes equations on unstructured grids. Specific improvement measures include:
1. **Consistency of mass flux**: Use a mass flux consistent with the mass conservation equation to ensure mass conservation of the two - phase momentum convection term during the implicit finite - volume discretization process.
2. **Solution of the momentum equation**: Use the SAAMPLE piecewise solution algorithm to solve the single - field Navier - Stokes equations.
3. **Numerical stability**: The proposed ρLENT method can achieve accurate numerical stability in the case where the density ratio ranges from 1 to 10 ^ 4, and is suitable for cases where the dynamic viscosity ratio is also in the range of 1 to 10 ^ 4 and there is a strong surface tension effect. In addition, this method also performs well in handling complex situations of strong interaction between the two phases, such as oscillating droplets and rising bubbles.
In conclusion, through proposing a new numerical method, the paper effectively solves the numerical stability and consistency problems in the two - phase fluid simulation under the condition of high - density ratio, providing important technical support for research in related fields.