|
[Sponsors] |
rhoReactingBuoyantFoam - Problems with gravity and Atmospheric BC |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 13, 2018, 19:06 |
rhoReactingBuoyantFoam - Problems with gravity and Atmospheric BC
|
#1 |
New Member
Henrik Hovrud
Join Date: Mar 2018
Location: Norway
Posts: 2
Rep Power: 0 |
Hello fellow foamers!
I am having trouble choosing the correct boundary-conditions for my gas dispersion case, where H2 gas and air is to flow freely out of the domain, into the atmosphere. Air should also be allowed to flow freely into the domain. With gravity set to a non-zero value in "constant/g" I experience a wall-like effect as the gas cloud nears the boundary, as if the H2 does not escape the domain. The flow velocity is non-zero, albeit small in y-direction on the boundary, so some gas/air clearly escapes. This problem is not present when gravity is set to "( 0 0 0 )", which I find odd... The effect seems limited to the y-normal boundary (the top) which is also normal to the direction of "g". I have tried different models for the 5 sides, including:
About the case:
The pictures show a jet of pure H2, an obstacle, H2 threshold between (0.01, 1), ground as a solid colour and the outline of the domain. "g" set to ( 0 -9.81 0 ) "g" set to ( 0 0 0 ) Vertical jet, to confirm gravity & buoyancy is working with g set to ( 0 9.81 0 ) Video, with better resolution and gravity enabled: https://www.dropbox.com/s/ug3j1qi24zjbxhc/sim.ogv?dl=0 If anyone has some insight into what I should do differently, it would be greatly appreciated Sincerely ~Hovrud Last edited by Hme; March 14, 2018 at 07:36. |
|
March 12, 2019, 17:16 |
|
#2 |
New Member
Join Date: Mar 2019
Posts: 1
Rep Power: 0 |
Hello Henrik,
I am having the same issue when applying these boundary conditions in buoyantBoussinesqPimpleFoam and twoLiquidMixingFoam. I am guessing the reason may be the zeroGradient boundary condition (including other boundary conditions that switch from zeroGradient and fixedValue depending on the flow direction at the boundary, e.g. freestream) that I used for temperature (in buoyantBoussinesqPimpleFoam)/ alpha (the volume fraction in twoLiquidMixingFoam). What boundary conditions did you use for your H2 and air? Do you have any progress in finding the reason solving this problem? Kind regards SL |
|
March 16, 2019, 15:12 |
|
#3 |
Member
Kirk Jarvis
Join Date: Mar 2009
Posts: 31
Rep Power: 17 |
Try type symmetry for the BC of the sides you should get the results you are expecting. You will have to set the sides to symmetry in the blockMesh also.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Boundary conditions when gravity is enabled | Cobra | STAR-CCM+ | 12 | January 13, 2020 07:38 |
Atmospheric BC in rhoReactingBuoyantFoam - Weird issues with gravity | Hme | OpenFOAM | 0 | March 12, 2018 09:06 |
buoyantBoussineqPimpleFoam Gravity Current | Bish | OpenFOAM | 0 | February 1, 2018 09:36 |
p_rgh outlet boundary condition for atmospheric simulation | Raphael_Santos | OpenFOAM Running, Solving & CFD | 0 | September 9, 2015 08:52 |
Gravity component and inlet velocity | Vidya | FLUENT | 8 | July 31, 2006 08:28 |