CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLUENT (https://www.cfd-online.com/Forums/fluent/)
-   -   buoyancy driven flow (https://www.cfd-online.com/Forums/fluent/35606-buoyancy-driven-flow.html)

buoyancy driven flow

Hi, Does anybody know how to treat density when modeling flow with high difference of altitude (order of 800 m)? I have pressure inlet and pressure outlet BC's.The fluid is air. The hydrostatic pressure must be removed from the static pressure at the inlet and otulet and the changes of pressure with altitude have to be taken into account after getting the pressure field. Is it the same case with density? For incompressible ideal gas the density will depend on temperature only (the operating pressure is constatnt). Should I let fluent to calculate the density field in such a way and then just calculate the field again with pressure containing the hydrostatic term to see the impact of height? Is it ok to do it like that? Best regards

 Evan Rosenbaum January 31, 2005 14:12

Re: buoyancy driven flow

Instead of using incompressible ideal gas, try ideal gas. I think that should account for pressure, although I've never actually tried such a problem.

Re: buoyancy driven flow

thanks for respond, Ideal gas will be even worse... If I used ideal gas, fluent will calculate the density according to current pressure (density=(operating pressure+local pressure)/(R/M*T). This will take into account the local pressure which is lower due to the assumption of removed hydrostatic pressure and thus the density will be also much lower. Then the density will have impact on sollution of the Boussinesq term in momentum equation. I am just not sure wether ideal gas or incmopressible ideal gas will better fit to the problem. It is sad that I cannot find much about it in the UG. regards adam

 giorgio February 10, 2005 13:26

Re: buoyancy driven flow

hi,

I performed a similar simulation. If you set the operating pressure to zero you obtain that local pressure=absolute pressure, the density should be calculated in a rigth way. But after that I had a lot of troubles with the turbulence values, TKE especially.

regards,

giorgio

Re: buoyancy driven flow

Hi Giorgio, Yes, you're right but did you have large difference in altitute? If you do, you need to redefine the pressure due to the hydrostatic head. The redefined static pressure is p's=ps-rho0*g*h. It doesn't matter if it is absolute or gauge - it is lower and the calculated density is also lower. Probably it is necessary to define my own formula for density calculation by UDF. Regards

 giorgio February 11, 2005 05:55

Re: buoyancy driven flow

in my simulations the height of the domain was 1 km. The formula you've writted is valid only if rho is constant. In the hydrostatic approximation you have to integrate the hydrostatic law from the ground to the top of your domain. Inside the integral form you have rho=rho(pressure(z), temperature(z)) so calculate this integral and you get a vertical pressure profile p=p(z) with z the altitude. Build a pressure profile file to import in fluent. I applied this profile at the outlet (pressure outlet) while at the inlet I had set the inlet as velocity inlet.

It should work

Giorgio