
[Sponsors] 
Fluid flow through wall U %3d 0 ignored when set pressure gradient for buoyancy 

LinkBack  Thread Tools  Display Modes 
February 29, 2008, 17:40 
Hi there,
I added buoyant f

#1 
Senior Member
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8 
Hi there,
I added buoyant force (1 + alpha*T)*g to an extended icoFoam solver, g = (0 9.81 0), then set BCs for p as pressure gradient such that eq grad(p) = g is satisfied. And what happens is  fluid goes through a boundary where BC U=(0 0 0) is set! Maybe you have an idea why? K. 

March 1, 2008, 01:27 
Karlis,
I'm guessing you'r

#2 
Senior Member
Michael Jaworski
Join Date: Mar 2009
Location: Champaign, IL, USA
Posts: 126
Rep Power: 8 
Karlis,
I'm guessing you're setting this BC at the bottom of the domain. Have you tried setting it to a fixedValue boundary with a value of rho*g*h where h is your fluid depth? That strikes me as a more direct approach. Regards, Mike J. 

March 1, 2008, 07:33 
the flow you are seeing is bec

#3 
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,769
Rep Power: 21 
the flow you are seeing is because of your PRESSURE boundary condition. A nonzero pressure gradient means there must be flow through the boundary because of the pressure Laplacian. You did not account for this in the pressure equation and hence you are getting nonzero flux.
Remember, boundary conditions on p and U are not independent  as I am sure you already know. Enjoy, Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk 

March 1, 2008, 12:34 
Harvoje,
I wouldn't ask and

#4 
Senior Member
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8 
Harvoje,
I wouldn't ask and make such useless experiments if I knew about BC nonindependence etc. Decided to implement gravitational force and temperature dependent density in icoFoam. Looking at the buoyantFoam solver I see some dynamic pressure introduced and only that is used for laplace( p ). Don't understand that stuff yet... > Where to read more about how to build buoyant solver for incompressible fluid? I read your thesis, Harvoje, but I would appreciate if someone point me in the right direction for buoyancy stuff.. > When I finish that solver, I would be happy to share it. Also I haven't found a BC patch type for convective + radiation heat transfer; BC is this: grad(T) =  surfaceNormalVector/k * ( h*(T  Treference) + sigma*epsilon(T^4  Treference^4) ) For sake of good usage constants h, Treference, epsilon must be read from 0/T file. I already begin a thread about this, but there is no reply so far. I wasn't able to access T field defined in solver from separate *.C file which defined class temperatureOpeningFvPatchScalarField : public fixedGradientFvPatchScalarField. Best regards, Kārlis 

March 1, 2008, 13:25 
Harvoje,
I found in an old

#5 
Senior Member
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8 
Harvoje,
I found in an old thread ( http://www.cfdonline.com/OpenFOAM_D...ges/1/815.html ), there was boussinesqBuoyantFoam by you, but link doesn't work... Maybe you could post it in again if you consider it useful? K. 

March 1, 2008, 23:53 
Karlis,
I do not believe t

#6 
Senior Member
Michael Jaworski
Join Date: Mar 2009
Location: Champaign, IL, USA
Posts: 126
Rep Power: 8 
Karlis,
I do not believe the pressure gradient term is correct as a boundary condition in the way you have used it. Namely, the gradient of the pressure is necessary *within* the fluid because of momentum balance. At a physical wall, however, pressure has a specific value. If you write out the equations for a momentum balance of the wall itself, the wallfluid force is the (pressure)*(area), and has nothing to do with the gradient. There already is a buoyant solver in OpenFOAM if I recall correctly. What does it specify for pressures? Regards, Mike J. 

March 2, 2008, 06:42 
Michael,
first about BC: ph

#7 
Senior Member
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8 
Michael,
first about BC: physically there is gradient. Is there any reason to say it cannot be implemented in OpenFOAM? Harvoje said: "You did not account for this in the pressure equation and hence you are getting nonzero flux." So, my problem (one of..) appears to be writing correct pressure eq! Second  I don't like specifying pressure because it is an integral quantity. There is static pressure and dynamic, but in case when density changes and flow is present, I'm not sure if it's still that simple. So I decided to use gradient as BC! About buoyantFoam  it is for compressible flow and tries to separate p = pd + rho*gh + pRef; But I currently cannot accept that. How to show that it solves the same NS equation with same conditions as for full pressure p? K. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
How can I get the pressure gradient along the direction normal to the local wall surface  qtian  OpenFOAM Running, Solving & CFD  3  April 4, 2014 05:05 
Pressure gradient for channel flow??  cfdIsMad  CFX  1  July 3, 2008 22:22 
Pressure boundary & buoyancy  Bouke  Main CFD Forum  2  October 22, 2002 11:37 
pressure gradient term in low speed flow  Atit Koonsrisuk  Main CFD Forum  2  January 10, 2002 11:52 
Pressure boundary on buoyancydriven flow  raymond  CDadapco  1  September 20, 2001 05:25 