CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Temperature field error in compressible LES (

aljo August 21, 2013 03:42

Temperature field error in compressible LES
3 Attachment(s)
Dear Foamers,

I encountered a serious problem in a high speed gas jet simulation (LES) with sprayFoam.

The jet velocity is 177 m/s with random perturbations added. The temperature at the inlet is set to 600K (!!). The velocity field (attachment) and pressure field are looking good. For the RANS simulation I get reasonable temperature profiles. For the LES (smagorinsky) the temperature in the field remains constant with a non-smooth transition.
Chemistry and so on are switched off.
I tried various fvSchemes (upwind, linear, ..). Also the transsonic option.
My guess is that it has something to do with the thermodynamic settings (same as for aachenBomb-Tutorial).

The only post I found that is dealing with a similar problem is

Can anyone give me a hint of how to solve this problem

Best regards


cosimobianchini August 22, 2013 04:12

Can you explain better what is plotted in the first two pictures (please provide scales)?
Can you provide a picture of the mesh for picture 2?


For the LES (smagorinsky) the temperature in the field remains constant with a non-smooth transition.
Can you explain better what you expect and why?

aljo August 22, 2013 10:39


Some more information about the case:
The jet diameter is 2.2 mm and the cell size at the outlet is about 0.05 mm. An o-block mesh has been used. For the contours the non-interpolated values have been plotted. So it is visible that the mesh is relatively fine at the inlet (velocity field).
The temperature field however is not so so well resolved. I would have expected a temperature field similar to the gas velocity field. The strangest thing is, that the temperature field is static for all time steps.

Hope the case is clear now


cosimobianchini August 23, 2013 03:40

From what you say I guess that first picture is velocity and second is temperature on the same cutting plane of the same mesh. In this case I agree with you temperature field looks wired.

If temperature field is the same at every iteration it might be that your fvSolution file is set so that energy is not actually solved.
Can you plot your fvSolution dictionary and a log of few iterations?
Also what solver are you using?
Try to use constant transport properties first (mu and Prandtl) instead of Sutherland as I'm not sure in which cases the formula implemented to estimate alpha (thermal diffusion coefficient) provides reasonable results.

aljo December 10, 2013 04:30

Problem solved
Dear Foamers,

I was able to solve this problem a while ago. It seems to be important for this specific chase to initialise the flow field with a streamwise velocity component u_field>5 m/s (u_jet = 240 m/s). If u_field is in the order of 1 the jet cannot develop well.
Still I think the core is somewhat short in comparison with fluent 13.0.0 .


All times are GMT -4. The time now is 18:22.