[yschati]

Hello
I am doing a multiphase flow simulation, with water and air (main phase) as the phases.
I have a duct. Air enters the duct (I have specified total temperature and total pressure as the inlet condition, setting air volume fraction = 1 and water volume fraction = 0). Inside the duct, I inject water through small holes (typically 2 mm in diameter). At the holes I give an inlet boundary condition of a lower total temperature but higher total pressure (water volume fraction = 1, air volume fraction = 0). Basically I want the water to cool the air. At the outlet, I have specified bulk mass flow rate. At the remaining surfaces I have the adiabatic wall (no slip for air, free slip for water droplets) conditions and the rotating periodic conditions. I am using k-epsilon turbulence model, buoyant flow, etc.
The problem is that the simulation is crashing. Either it gives a notice of very high Mach number or sometimes it just crashes without displaying any error. Microsoft gives an error that solver.hpmpi has crashed. I have tried changing the mesh, boundary conditions, turbulence models, but to no avail.
Can someone please advise what to do and what can be the possible reasons of failure?
Thanks.