[MisterAceituna]

Hi everyone!

I'm currently working on a case of an engine with the piston moving according to crank rod angle.

The simulation works fine with rhoPimpleFoam and similar solvers, though I'm encountering this issue:

when using both k-epsilon and k-omega, if I put a fixedValue at inlet, the simulations run nicely, while if I put the specification of turbulentMixingLengthDissipationRateInlet (or turbulentMixingLengthFrequencyInlet) , which is what I really want to use, the solvers crashes in few steps, generating the error of negative temperature.

What do you think can be the main problem in setting this BC?

Tell me if uploading the BC field or some dicts may help.

Thank you very much!