Abstract:A second-order accurate and robust numerical scheme is developed for the Kapila model to simulate compressible multiphase flows. The scheme is formulated within the temporal-spatial coupling framework with the generalized Riemann problem (GRP) solver applied as the cornerstone. The use of the GRP solver enhances the capability of the resulting scheme to handle the stiffness of the Kapila model in two ways. Firstly, in addition to Riemann solutions, the time derivatives of flow variables at cell interfaces are obtained by the GRP solver. The coupled values, i.e. Riemann solutions and time derivatives, lead to a straightforward approximation to the velocity divergence at the next time level, enabling a semi-implicit time discretization to the volume fraction equation. Secondly, the use of time derivatives enables numerical fluxes to comprehensively account for the effect of the source term, which includes interactions between phases. The robustness of the resulting numerical scheme is therefore further improved. Several challenging numerical experiments are conducted to demonstrate the performance of the proposed finite volume scheme. In particular, a test case with a nonlinear smooth solution is designed to verify the numerical accuracy.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to develop an efficient and robust numerical scheme for simulating compressible multiphase flows. Specifically, it addresses the following two main issues:
1. **Stiffness problems in the Kapila model**:
- **Stiff source terms in the volume fraction equation**: The source terms in the volume fraction equation in the Kapila model are very stiff, which is due to the stiff limit of the relaxation process in the original seven - equation model. This stiffness makes it very difficult to approximate the source terms and maintain the boundedness of the volume fraction.
- **Non - monotonic behavior of Wood's sound speed**: The non - monotonic characteristics of Wood's sound speed related to the volume fraction can cause Mach number oscillations at the interface, resulting in instability in the numerical diffusion region.
2. **Improving the robustness and accuracy of the numerical scheme**:
- **Application of the time - space coupling framework**: By introducing the generalized Riemann problem (GRP) solver, this method enhances the ability to handle the stiffness of the Kapila model. In particular, the GRP solver not only provides the Riemann solution but also obtains the time derivatives of the flow variables at the interface, which helps to more accurately capture the inter - phase interactions and achieve semi - implicit time discretization.
- **Improving the numerical flux**: Using the time derivatives, the numerical flux can more comprehensively consider the influence of the source terms, thereby improving the robustness and accuracy of the numerical scheme.
### Specific contributions of the paper
- **High - resolution and robust numerical scheme**: A finite - volume scheme based on the time - space coupling framework has been developed, which is especially suitable for the Kapila model and can maintain high computational accuracy while dealing with stiffness problems.
- **Semi - implicit time discretization**: Through the time derivatives obtained by the GRP solver, semi - implicit time discretization of the volume fraction equation has been achieved, effectively alleviating the stiffness problem of the volume fraction.
- **Accurate analysis of multi - dimensional effects**: The GRP solver can better analyze multi - dimensional effects, especially in the presence of strong non - linear waves, effectively solving the evolution problem of basic physical information.
### Summary of mathematical formulas
- **Basic equations of the Kapila model**:
\[
\begin{aligned}
&\frac{\partial (\zeta_1 \rho)}{\partial t}+\nabla\cdot(\zeta_1 \rho\mathbf{u}) = 0,\\
&\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{u}) = 0,\\
&\frac{\partial (\rho\mathbf{u})}{\partial t}+\nabla\cdot(\rho\mathbf{u}\otimes\mathbf{u})+\nabla p = 0,\\
&\frac{\partial (\rho E)}{\partial t}+\nabla\cdot[(\rho E + p)\mathbf{u}] = 0,\\
&\frac{\partial \alpha_1}{\partial t}+\mathbf{u}\cdot\nabla\alpha_1=-\alpha_1(1 - \alpha_1)\left(\frac{\rho_1 c_1^2-\rho_2 c_2^2}{\alpha_1\rho_2 c_2^2+(1-\alpha_1)\rho_1 c_1^2}\right)\nabla\cdot\mathbf{u}.
\end{aligned}
\]
- **Sound speed of the mixture**:
\[
\frac{1}{\rho c^2}=\frac{\alpha_1}{\rho_1 c_1^2}+\frac{\alpha_2}{\rho_2 c_2^2}.
\]