[hannes]

Hello all,

I've tried for some time now to set up a steady compressible flow simulation with rhoSimplecFoam (tried rhoSimpleFoam at first, but the SIMPLEC variant seems to be more stable).

I have chosen the converging channel validation case for the compressible SIMPLE algorithm from the Ferziger&Peric Book pp. 323 (You can find the case attached to this message).
Prescribing the flowrate seemed to be quite shaky, so I prescribed the pressure ratio, which resulted in a better convergence. The inlet mach number is then a bit higher than in the case presented in the book.

According to the literature, there should be some shocks originating at the throat, but they do not appear although the case converges (see attached picture).
The picture does not change qualitatively if I use higher order schemes instead of upwind for convection-like terms.

Since there are pretty nice looking results presented in the Ferziger&Peric Book for a compressible SIMPLE solver, I expected to see something similar with the OpenFOAM solver. Do I expect too much?
Are there any suggestions on the scheme selection for this type of flows?
Any hints are welcome.

Thanks in advance, Hannes