CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   CFX (https://www.cfd-online.com/Forums/cfx/)
-   -   Compressible amosphere using opening (https://www.cfd-online.com/Forums/cfx/160921-compressible-amosphere-using-opening.html)

Dino9111 October 15, 2015 13:52

Compressible amosphere using opening
 
3 Attachment(s)
Hi together,

I would like to simulate just a area (box with dx=50m x dy=50m x dz=25m) with compressible air using opening boundary at the sides and a wall at the bottom. Intention is to create a realistic atmosphere with a negative temperature gradient with increasing height without any flow (or almost).

Initial and boundary condition for temperature: T(y)=288K-(0.976 K/100m)*y
(using adiabatic temperature gradient of atmosphere g_ad)
ref. pressure: 1 atm at (0,0,0)
If I set "relative pressure" at "opening pressure and drn." at 0 Pa, I always get a flow downwards at the center of the box (when compressibility and byouyancy is enabled - I need both effects in later simulations).
I'm not sure but as far as I understood the doc, I should define hydrostatic pressure p(y) at my opening boundaries.
--> Initial and boundary condition for pressure: p(y)-pref with
p(y)=p0*(1-(g_ad*y/T0)^(M_air*g/(R*g_ad)) (barometric formula; p0,T0 at y=0)

But then, in the center there is a flow upwards.

Has anyone an idea how to set boundary conditions etc. correctly?
Thanks a lot in advance!
Attachment 42722

Attachment 42723

Attachment 42724

evcelica October 16, 2015 06:22

Translational periodic conditions might be suitable here.

Dino9111 October 20, 2015 05:17

Hi evcelica,

first of all thanks - it works qute well. But unfortunately then I can't define a wind profile at the sides (for later siimulations).
Isn't there any possibility to define relative pressure at an opening correctly (so that there is no flow initiated by pressure differences)?
And why doesn't relative pressure = 0 not work? No pressure difference-->no flow?!

Dino9111 October 21, 2015 07:15

OK i finally solved it...
for those who are interested:
You have to define for the realtive pressure at the opening boundary: p(y) - pref + density*g*y because p(y) is the absolute pressure, but the pressure you have to define at the opening boundary is interpreted to be without the hydrostatic part.


All times are GMT -4. The time now is 10:03.