[SimulatorMaster]

Hey all,

An air-blast atomizer is being analyzed by utilizing DPM along with wall-film models activated on 3 different walls one of which a jet is impinging.

The activate impingement/splashing model is Stanton-Rutland while the chosen separation model is Foucart with specified critical angle and weber constants.

Moreover, activated DPM models are
-- Saffman Lift Force
-- Two Way Turbulence Coupling
--Breakup
--Stochastic Collision with Coalescence

Breakup Model: KHRT

DPM iterations are kept as low as possible: 5
Unsteady Tracking with 10^-6s time step and 10 number of steps

Turbulence Model: Kw-SST with curvature correction, and Production Limiter

The flow is steady. Unsteady is tried as well but did not work.

While running the simulation with these settings, there seems to be no convergence problems. However, whenever "Particle Stripping" model is activated, almost all of the residuals shoot up, try to converge for a while, then give errors.

At some point, the temperature limit is reached (5000K but decreased down to 400K) and turbulent viscosity limit is reached as well. At this point, the continuity and omega residuals go up to the moon.

A reason could be sudden increase in particels in mesh cells. So, the critical shear value is increased but still no use. Volume Displacement model is tried as well but this did not work. The only thing working was deactivating two-way turbulence coupling model which somewhat increased omega convergence yet the problems still persist. I believe my mesh is decent as well.

What could be the reason?