[David1]

Hello,

Does anyone know why for the transonic case in rhoSimpleFoam the flux is calculated as:

phi == pEqn.flux(); ?

If transonic is turned off, then the flux is calculated in the 'normal' way, i.e. first calculting the flux without the contribution of the pressure gradient:

phi = fvc::interpolate(rho)*(fvc::interpolate(U) & mesh.Sf() )

And after the calulation of the pressure, the old flux value gets corrected using the new pressure:

phi -= pEqn.flux();

So I do not understand why the flux in the transonic case is not calculated as:

phi == fvc::interpolate(rho)*(fvc::interpolate(U) & mesh.Sf() ) - pEqn.flux()

as well.

I hope somebody can help me with this.
Thank you in advance!