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 ... Gloq |
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. ;) |
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. |
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. ;) |
Hello.
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? Lasse |
All times are GMT -4. The time now is 21:59. |