CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Temperature/Energy in sonicFoam (http://www.cfd-online.com/Forums/openfoam-solving/78456-temperature-energy-sonicfoam.html)

Ruehri July 22, 2010 04:48

Temperature/Energy in sonicFoam
 
Dear Foamers,

I have seen several threads about getting wrong results with the solver sonicFoam and have unfortunately got a similar experience with this solver.

In my specific case, it seemed that the heat/energy flux was too high in the solution, that's why I searched the implementation for it in the code and I found 2 issues that bug me about this solver.

1. The internal energy and temperature are calculated in a predictor step, but not corrected after solving the pressure equation (at least I couldn"t find where). Wouldn't this result in a fairly bad calculated temperature field? It seems to me, that the predictor step should actually yield the corrected temperatures for the previous step.

2. The pressure equation in sonicFoam is derived from continuity, momentum and the equation of state. Wouldn't this mean that the corrector step actually doesn't solve for energy conservation? Even if the predictor obeys the energy conservation, couldn't there still could be huge issues with the corrector step?

It would be great if someone could have a look at it and confirm my assumptions, I think it would be at least fairly easy to build in a corrector for temperature in the solver.

Best,
Yu

KrisT July 24, 2010 03:19

Quote:

Originally Posted by Ruehri (Post 268466)
It would be great if someone could have a look at it and confirm my assumptions, I think it would be at least fairly easy to build in a corrector for temperature in the solver.

I too think it would be great if sonicFoam could be fixed - or removed from the distribution - it is kind of stupid to keep it while it gives the wrong results (even on the tutorial cases included with the distribution).


All times are GMT -4. The time now is 04:40.