Property-preserving numerical approximation of a Cahn-Hilliard-Navier-Stokes model with variable density and degenerate mobility

Daniel Acosta-Soba,Francisco Guillén-González,J. Rafael Rodríguez-Galván,Jin Wang

DOI: https://doi.org/10.1016/j.apnum.2024.11.005

2024-11-28

Abstract:In this paper, we present a new computational framework to approximate a Cahn-Hilliard-Navier-Stokes model with variable density and degenerate mobility that preserves the mass of the mixture, the pointwise bounds of the density and the decreasing energy. This numerical scheme is based on a finite element approximation for the Navier-Stokes fluid flow with discontinuous pressure and an upwind discontinuous Galerkin scheme for the Cahn-Hilliard part. Finally, several numerical experiments such as a convergence test and some well-known benchmark problems are conducted.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is how to preserve the key physical properties of the Cahn - Hilliard - Navier - Stokes (CHNS) model in numerical simulations, especially when the model involves variable density and degenerate mobility. Specifically, the authors hope to develop a numerical approximation method that can ensure the following three important properties:

1. **Mass conservation**: The total mass of the mixture remains unchanged throughout the simulation process.
2. **Point - wise bounds on density**: Ensure that the phase - field function ($\phi$) always stays within its physically allowed range (usually $[- 1,1]$) in the numerical solution, thereby guaranteeing that the mixture density $\rho(\phi)$ is within its physically reasonable range (i.e., $[\rho_1,\rho_2]$).
3. **Energy stability**: The total energy of the system does not increase with time, that is, it satisfies the discrete energy dissipation law.

To achieve these goals, the authors propose a new numerical framework based on the finite element method (for the fluid flow part in the Navier - Stokes equation) and the discontinuous Galerkin (DG) upwind scheme (for the Cahn - Hilliard equation part). Through this method, they can not only accurately simulate complex fluid mixing behaviors but also ensure the numerical fidelity of the above - mentioned physical properties.

In addition, the paper also verifies the effectiveness of this method through a series of numerical experiments, including convergence tests and some classic benchmark problems, such as the mixing of two bubbles and unstable phenomena under the action of gravity. These experimental results are consistent with the theoretical analysis, further demonstrating the reliability and superiority of the proposed numerical method.

In summary, this research aims to provide a numerical tool that is both accurate and stable for simulating multiphase fluid systems with complex physical properties, especially in cases involving fluids of different densities.