I need explanations about fixedFluxPressure
I had a look at the fixedFluxPressure boundary condition and I am not quite sure to well understand what this BC does.
I guess I have to use this BC if I want to evaluate the pressure according to a velocity input. Am I wrong ?
Furthermore, in the code source, one can read:
All enlightments are welcome!
buoyantPressure vs. fixedFluxPressure
I realize this is an old thread, but it is still very relevant.
I was having terrible problems trying to get buoyantBoussinesqSimpleFoam to converge. I changed the p_rgh BC's from buoyantPressure to fixedFluxPressure and the model started to converge.
Same question here. What is the difference between buoyantPressure and fixedFLuxPressure boundary conditions?
OpenFAOM 2.3.0 doesn't have buoyantPressure BC
OpenFAOM version 2.3.0 doesn't have buoyantPressure libraries. Am I right?
what is the difference between buoyantPressure and fixedFluxPressure?
why fixedFluxPressure converge better? is it the only difference between them?
OpenFOAM 2.2.2 has both boundary conditions. In newer Version, buoyantPressure boundary was dismissed. For clarification have a look into the code or into Doxygen.
From buoyantPressureFvPatchScalarField we learn that:
This boundary condition sets the pressure gradient appropriately for buoyant flow. If the variable name is one of pd, p_rgh or ph_rgh, we assume that the pressure variable is http://latex.codecogs.com/gif.latex?...%5Ccdot%20h%29 and the gradient set using:
Otherwise we assume that it is the static pressure, and the gradient calculated using:
From fixedFluxPressureFvPatchScalarField we see that:
This boundary condition adjusts the pressure gradient such that the flux on the boundary is that specified by the velocity boundary condition.
The predicted flux to be compensated by the pressure gradient is evaluated as http://latex.codecogs.com/gif.latex?...i_%7BH/A%7D%29, both of which are looked-up from the database, as is the pressure diffusivity used to calculate the gradient using:
Some additional remarks:
The fixedFluxPressure boundary is known do have a better convergence. This is in my opinion due to the more pressure related calculation of the gradient:
phi and phiHbyA (predicted flux field) are directly linked to pressure equation. On a wall, the flux difference tends to zero and so the boundary condition turns to zeroGradient.
The buoyanPressure boundary is more the physical method to describe the pressure gradient with Archimedes' principle.
I hope my explanation is understandable and correct.
Hi folks, I'm dealing with fixedFluxPressure BCs also and this thread was helpful. Checking the code things have changed at least in FOAM 2.3.0. The class definition is now more simple
Hope this helps for 2.3.0 users.
I would like to comment out a bit, why does OpenFOAM use the expression given by santiagomarquezd (gracias por el análisis del código fuente) for the pressure surface normal gradient in the boundaries.
First of all we need to write down the momentum equation discretized using the Rhie-Chow interpolation method. This reads:
Where the f subscript means interpolate to the face, is the variable p_rgh whose value is , and the vector u is the velocity. Now multiply the equation by the face surface vector where the vector n is the surface normal:
Thus the following equation for the pressure surface gradient can be obtained:
Hope this is useful for you foamers
|All times are GMT -4. The time now is 05:55.|