[Zato_Ichi]

Greetings !

I tried to simulate a 1D shock wave generation in closed domain. Pressure field initially was set up with discontinuous conditions, from funkySetFields dictionary:

field p;
expression "3";
condition "(pos().x)>=(1)";
keepPatches 1;

So pressure equals 1 Pa on left and 3 Pa on right.

Boundary conditions for T and p are zeroGradient on both sizes, and U equal (0 0 0) on both sizes, so domain is closed by walls on right and left.
For calculating i used sonicFoam and rhoCentralFoam. And result of simulation was quite strange. Shock wave appears and propagates normally, but after reflection by the wall velocity U turns into zero, while pressure and temperature experience jumps, as it should be. As i know, in theory velocity must change sign and value (but not become a zero).

Maybe I'm putting wrong boundary conditions. I tried slip condition on velocity, but it didn't help.