rhoPimpleFoam with fixedFluxPressure + flowRateInletVelocity --> p diverges
 

Hello!

I'm trying to specify mass flow at compressible subsonic inlets, and I have understood that this can be done using fixedFluxPressure and flowRateInletVelocity. This seems to work fine for rhoCentralFoam, but not for rhoPimpleFoam; unless I assign a very low mass flow (flowRate), p diverges (to negative and positive value) usually on the first time step. I'm using OF 2.1.1.

Any ideas how to overcome the p divergence? For fixedFluxPressure I need to store the matrix rAU = 1/UEqn.A(). I guess I should calculate rAU for fixedFluxPressure after the U is completely solved for that time step, i.e. after line 85 in pEqn.H, right?

By the way, in 2.1.x there is a new feature in fixedFluxPressure: pressure diffusivity. I guess this is meant to "smear" the pressure gradient, or is it? Could this be of help for my case? Also, instead of dividing the phi with density (in compressible solvers), the solver seems to use the field phi*UEqn.H()/UEqn.A() for adjusting the pressure gradient.

Tatu