On Riemann solvers and kinetic relations for isothermal two-phase flows with surface tension (Q721334)

The compressible two-phase flow is studied in the paper. The open bounded domain \(\Omega\subset\mathbb{R}^m\) is filled with liquid and gas. So \(\Omega\) is portioned into the time-depended two open sets \(\Omega_l(t)\) (liquid) and \(\Omega_g(t)\) (gas). The surface \(\Gamma(t)\) is the common boundary of \(\Omega_l(t)\) and \(\Omega_g(t)\). The motion is governed by the equations for free boundary problem \[ \begin{cases}\frac{\partial \rho}{\partial t}+\operatorname{div}(\rho v)=0,\quad x\in \Omega_l(t)\cup\Omega_g(t),\;t\in(0,T), \\ \frac{\partial }{\partial t}(\rho v)+\operatorname{div}(\rho v\otimes v+p(\rho)\mathbf{I})=0,\quad x\in \Omega_l(t)\cup\Omega_g(t),\;t\in(0,T), \\ \left[\rho(v\cdot n-V_n)\right]=0\quad \text{on}\quad\Gamma(t),\\ \left[\rho(v\cdot n-V_n)v\cdot n+p(\rho)\right]=(m-1)\sigma K\quad\text{on}\;\Gamma(t),\\ \left[v\cdot\tau_i\right]=0\quad\text{on}\;\Gamma(t),\;i=1,2,\dots,m-1. \end{cases} \] Here \(\rho(x,t)\) is the unknown density, \(v(x,t)=(v_1,\dots,v_m)\) is the unknown velocity field, \(p(\rho)\) is a given function, \(\mathbf{I}\) is the unit \(m\times m\) matrix, \(n\) is the unit normal on \(\Gamma(t)\), \(V_n\) denotes the speed of \(\Gamma(t)\) in the normal direction \(n\), \(\sigma\) is the surface tension coefficient, \(K\) is the mean curvature of \(\Gamma(t)\), \(\tau_i\), \(i=1,\dots,m-1\) are unit vectors tangential to \(\Gamma(t)\). The symbol \([a]\) denotes the jump of the function \(a\) on the surface \(\Gamma(t)\). In addition the density fields \(\rho_l\) and \(\rho_g\) in liquid and gas regions accordingly and the velocity fields \(v_l\), \(v_g\) satisfy to the specific kinetic relation \[ F(\rho_l,v_l,\rho_g,v_g)=0. \] The authors investigate the free boundary problem analytically and numerically for physically relevant kinetic relations using the theory of monotone decreasing functions. It should be noted that the kinetic relation \(F\) has a very complex form. Nevertheless, according to the authors, the proposed method of numerical analysis of two-phase flow is very effective.