Numerical Solution to the Riemann Problem for a Liquid-Gas Two-phase Isentropic Flow Model

TL;DR

This paper numerically solves the Riemann problem for a two-phase liquid-gas flow model using finite-difference schemes, highlighting the effectiveness of the TVD FORCE scheme in capturing complex wave interactions.

Contribution

It introduces a numerical approach to solve the Riemann problem for a specific two-phase flow model, comparing multiple schemes and demonstrating the advantages of the FORCE scheme.

Findings

The FORCE scheme is TVD and monotone, avoiding oscillations.

Numerical solutions show smooth profiles, shock waves, and rarefaction waves.

The scheme comparison highlights the performance differences among Lax--Friedrichs, Lax--Wendroff, and FORCE.

Abstract

A recently introduced two-phase flow model by Chun Shen is studied in this work. The model is derived to describe the dynamics of immersed water bubbles in liquid water as carrier. Several assumptions are made to obtain a reduced form of the mathematical model. The established model consists of nonlinear coupled PDEs in which the unknowns are the densities of the liquid and gas phases and the velocity of the liquid phase; these depend on space and time. For numerical purposes a one-dimensional space--time coordinate system $(x,t)$ is considered. Using the Python programming framework, several Riemann-type initial value problems for the two-phase flow model are solved numerically. A comparison of three finite-difference schemes is presented in order to examine their performance: the Lax--Friedrichs scheme, the Lax--Wendroff scheme, and the FORCE scheme. The FORCE scheme is total-variation diminishing (TVD) and monotone and does not create oscillations. As expected, the numerical solution of the Riemann problems consists of combinations of smooth profiles, shock waves, and rarefaction waves.