[Taendyr]

Hi,



I'm trying to simulate a mix of gas and dust against a disc using twoPhaseEulerFoam. After some iterations, the solver stops without any error message, but we can see some negative temperature for the air:
Courant Number mean: 0.122993 max: 0.401974
Max Ur Courant Number = 0.0456224
Time = 0.195

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 1.22398e-08, Final residual = 2.15733e-15, No Iterations 1
alpha.particles volume fraction = 0.199907 Min(alpha.particles) = 0.00159117 Max(alpha.particles) = 0.992272
Constructing momentum equations
smoothSolver: Solving for e.particles, Initial residual = 0.0643472, Final residual = 2.2115e-12, No Iterations 3
smoothSolver: Solving for e.air, Initial residual = 0.252945, Final residual = 2.55661e-13, No Iterations 5
min T.particles 119.995
min T.air -18672.9
GAMG: Solving for p_rgh, Initial residual = 0.999183, Final residual = 3.63007e-18, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.40899e-18, Final residual = 1.40899e-18, No Iterations 0
PIMPLE: iteration 2
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 1.07889e-08, Final residual = 2.02332e-15, No Iterations 1
alpha.particles volume fraction = 0.199984 Min(alpha.particles) = 0.00158145 Max(alpha.particles) = 0.996012
Constructing momentum equations.


I tried to change the pressure tolerance, to change deltaT,etc.. I always have this behaviour.


How can I fix it ? Imposing a limit to the temperature doens't work.



I joined the files case.