There is a gas-water two-phase compressible multifluid Riemann problem, governed by the stiffened gas equation of state:

p(\rho,e)=(\gamma-1)\rho e - \gamma p_\infty

For gas, \gamma=1.4, p_\infty=0; For water, \gamma=4.4, p_infty=6.E+8

And the initial conditions are For gas, \rho=1.2 Kg/m^3, p=1.E+5 Pa, velocity=0; For water,\rho=1000 Kg/m^3, p=1.E+9 Pa, velocity=0.

I have used a simple mixture-type Eulerian numerical method proposed by Saurel/Abgrall(SIAM 1998) or Shyue(JCP 1999) to simulate this problem, but the result is very poor because the pressure curve is very different from the exact solution. It is very obvious in the pressure plot for using log curve.

I think this problem is very stiff and the density ratio and pressure ratio is very high(1000/1.2 and 10000/1).

Has someone done some work in this problem and to tell me how to solve this problem?

Ma, Dong-Jun

It seems to me it is a numerical problem. why not to have a look to this book:

Riemann solvers and numerical methods for fluid dynamics.

By: E.F. Toro

Springer.

Hi,

Are you using an explicit or implicit scheme for time integration. I used an implicit method to solve the 2-phase problem in K. M. Shyue, JCP 1998. But in this case the density ratio is not as high as yours. but the pressure ratio is high.

Do you know how to compute temperature on the water side.

Buvana

I'm working with isentropic two-fluid model and using exact Riemann solver, explicit scheme.

Hi,

I was also doing this problem last year, just in my work, for the water \gamma is taken as 7.0. I am thinking the major problem in your calculation may be the treatment of flux through air-water interface. Please check whether you keep the interface like the wave behaves itself.

Good luck!