Dear All! I need some help. I try to calculate shock tube problem when the high pressure chamber is filled with gas-particle dispersion. I use Mac-Cormack method with implicit account of interaction terms.

When I've used isothermal approach everything worked fine. When I've added energy equation I got negative temperatures and pressures after the shock wave at the beginning. Later on proper shock wave forms. Additionally I need very small CFL numbers ~ 0.05 whereas for isothermal case 0.2-0.3 was good. I am trying to simulate high pressure drops – 20-40 MPa.

Any suggestions?

Oleg

What kind of artificial viscosity are you using?

I use the term eps*(V_i+1 - 2*V_i + V_i-1) in momentum equation and play with eps. V is gas velocity.

Oleg

The problem is solved now. I've added TVD to Mac-Cormack method and temperature behaves properly.

Now there is one more issue. I have a closed end of the high pressure chamber. The material of the wall is metal with high conductivity. The behaviour of the pressure at the wall (which we measure) strongly depends on temperature BC. I've tried constant temperature and adiabatic wall. The real one should be in between. Could you help me with a proper BC for the temperature?

Best regards,

Oleg