A Second Order Fully-discrete Linear Energy Stable Scheme for a Binary Compressible Viscous Fluid Model
This work provides a provably stable and efficient numerical method for simulating compressible fluid mixtures, which is important for computational fluid dynamics and materials science.
The authors developed a linear, second-order fully discrete numerical scheme for a binary compressible viscous fluid model that preserves mass conservation and energy dissipation. The scheme was validated through mesh refinement tests and examples of spinodal decomposition and interface evolution, demonstrating convergence and applicability.
We present a linear, second order fully discrete numerical scheme on a staggered grid for a thermodynamically consistent hydrodynamic phase field model of binary compressible fluid flow mixtures derived from the generalized Onsager Principle. The hydrodynamic model not only possesses the variational structure, but also warrants the mass, linear momentum conservation as well as energy dissipation. We first reformulate the model in an equivalent form using the energy quadratization method and then discretize the reformulated model to obtain a semi-discrete partial differential equation system using the Crank-Nicolson method in time. The numerical scheme so derived preserves the mass conservation and energy dissipation law at the semi-discrete level. Then, we discretize the semi-discrete PDE system on a staggered grid in space to arrive at a fully discrete scheme using the 2nd order finite difference method, which respects a discrete energy dissipation law. We prove the unique solvability of the linear system resulting from the fully discrete scheme. Mesh refinements and two numerical examples on phase separation due to the spinodal decomposition in two polymeric fluids and interface evolution in the gas-liquid mixture are presented to show the convergence property and the usefulness of the new scheme in applications.