CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   rhoPimpleFoam with fixedFluxPressure + flowRateInletVelocity --> p diverges (https://www.cfd-online.com/Forums/openfoam-solving/113661-rhopimplefoam-fixedfluxpressure-flowrateinletvelocity-p-diverges.html)

 tatu February 24, 2013 15:20

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

 styleworker March 21, 2013 16:10

In alphaContactAngleFvPatchScalarField.H is a note:

"
Note that if any of the first three options are used the boundary conditionon p_rgh must set to guarantee that the flux is corrected to be zero at the wall" (fixedFluxPressure)

Does it mean the flux is corrected to be zero across the wall, or along the wall?

 All times are GMT -4. The time now is 07:27.