A note on hyperbolic relaxation of the Navier-Stokes-Cahn-Hilliard system for incompressible two-phase flows

¹ Institute of Aerodynamics and Gas Dynamics, University of Stuttgart, Wankelstraße 3, 70563 Stuttgart, Germany
² Institute of Applied Analysis and Numerical Simulation, University of Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany

^* e-mail: jens.keim@iag.uni-stuttgart.de
^** e-mail: Hasel-Cicek.Konan@mathematik.uni-stuttgart.de
^*** e-mail: christian.rohde@mathematik.uni-stuttgart.de

Abstract

We consider the two-phase dynamics of two incompressible and immiscible fluids. As a mathematical model we rely on the Navier-Stokes-Cahn-Hilliard system that belongs to the class of diffuse-interface models. Solutions of the Navier-Stokes-Cahn-Hilliard system exhibit strong non-local effects due to the velocity divergence constraint and the fourth-order Cahn-Hilliard operator. We suggest a new first-order approximative system for the inviscid sub-system. It relies on the artificialcompressibility ansatz for the Navier-Stokes equations, a friction-type approximation for the Cahn- Hilliard equation and a relaxation of a third-order capillarity term. We show under reasonable assumptions that the first-order operator within the approximative system is hyperbolic; precisely we prove for the spatially one-dimensional case that it is equipped with an entropy-entropy flux pair with convex (mathematical) entropy. For specific states we present a numerical characteristic analysis. Thanks to the hyperbolicity of the system, we can employ all standard numerical methods from the field of hyperbolic conservation laws. We conclude the paper with preliminary numerical results in one spatial dimension.