[Piotr95]

Hello,
I have a problem with the multi-phase transonic simulation. It diverges. I try to make a simulation of DeLaval nozzle, where shock waves appear.

Firstly I did a simulation of single-phase, it works perfectly fine, as you can see on the picture below.

Settings:
Inlet
Total Pressure = 98880 Pa
Total Temperature 297.15 K
Outlet
Static Pressure = 55000 Pa

Operating pressure = 0Pa
Ideal-gas
k-w SST, y+ ~1,



Then I tried to move it to Euler-Euler framework. The secondary phase is water droplets. The dispersed turbulence model is on. The volume fraction of water is 1e10-6 and the droplet diameter is 1e10-8. Although it start to compute easily, it diverges after some time.

I have already tried transient and steady approach, all discretization schemes, Simple and Couple methods, Implicit and Explicit formulation for volume fraction. None of them work. What is happening:
-The solution is getting shifted away from the nozzle center line
-The pressure between shock wave and outlet starts to fluctuate

-Mass imbalance starts to appear on the shock wave

Either The Fluent is not capable of solving the interaction of droplets with shock waves or I am doing something wrong. I would be honestly grateful for any suggestions.