[fportela]

Hello all,

I have been trying to get a solution for a supersonic turbulent channel flow.

I managed to get sonicFoam working for a Mach number (based on wall temperature) of about 1.5 - t1.png. To achieve this I simply add a momentum source term to the UEqn and eEqn.

When I try to get to M ~ 3 - t2.png - the simulation crashes during laminar-turbulent transition with the following message:

[15] --> FOAM FATAL ERROR:
[15] Maximum number of iterations exceeded
[15]
[15]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const

I am not sure what can be going on, so if you need more information to help me sort this out, just let me know!

Cheers,
Felipe