CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Problem with BC + ReactingFOAM + gravity (

Gloq November 27, 2012 06:29

Problem with BC + ReactingFOAM + gravity
Hi guys,

I searched for a problem looking like mine on the forum, but I could not find such a topic, so here is a new one !

My computation domain is a channel, in which air is flowing, while an other gas is injected in the middle.
To adress this problem, I use reactingFOAM.
For my first runs, I forgot to set the gravity vector to -9.81 u_z.
I obtained physical pressure and velocity fields.

But when I realized my mistake, I 'activated' gravity. Doing so, I got weird velocity and pressure fields for the corresponding runs.
Especially close to outlet BC where a recirculation zone appeared.
Big differences between theoretically very close runs.

(If useful : outlet BC are zeroGradient for U and fixedValue (1 atm.) for p)

How do you explain the differences ?
Which solution could you advise to me ?
I was considering the option to try to create a simplified version of reactingFOAM, as my case is non-reactive, but I'm not sure that this work would help me ...

Thanks in advance for your help ...


vainilreb November 27, 2012 06:48

Could you please post your files from the 0-directory and a screenshot of the strange flow behaviour you're talking about?

I'm dealing with a similar problem and perhaps I can help you. But I need a little more detailed information. ;)

Gloq November 27, 2012 08:49

I'll try to upload more precise elements. For the fields, I have roughly :
- a recirculation zone at the outlet
- a depression at the outlet. In the channel, pressure is not homogeneous. I mean, pressure is slightly higher in the lower part than in the upper part of the channel. For the case ran w/o gravity, pressure is homogeneous.

vainilreb November 27, 2012 09:08

The pressure gradient is typical of gravity-influenced cases.
It should be something like dp/dy = g*rho [Pa/m]. (Please correct me, if I'm misled.) You will find this pressure gradient in our atmosphere, too.

Concerning the recirculation zone I can't say anything until I've seen your BC files and a screenshot would really help my imagination, too. ;)

redhoax May 9, 2013 21:02


I am very interested in the solution to this matter. I am also doing a case on air flow and gas injection, using reactingFoam. Like Gloq I obtain this circulation zone at the outlet when gravity is set.

Did you discover a solution to the matter?

Could it be possible to use f.ex. fireFoam, or rewrite solver to solve for p_rgh?


All times are GMT -4. The time now is 14:33.