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

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).