Temperature/Energy in sonicFoam
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.
|All times are GMT -4. The time now is 19:40.|