[AnasCFD]

Dear respected researchers,

I am simulating a two phase flow problem in a rotating domain.
In this case there are liquid water, air, and a well mixed water-air regions. For this case I used compressibleInterDyMFoam.

Unfortunately I have faced very high velocities at the interfaces (in the order of 300 to 400 m/s). I am thinking of two reasons to explain this result:
1- Numerically. Due to the well known "parasitic currents" problem.
2- Physically. It is known that the speed of sound for a highly mixed two phase flow (like water and air) decreases dramatically (up to to 20 m/s) and hence this affects the compressibility of the "mixture" at this region.

The questions are now how to reduce the first effect (other than playing with the time step and grid size) and how to include the effect of the speed of sound on the model.