[Gloq]

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 ...

Gloq