|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Join Date: Dec 2009
Posts: 33
Rep Power: 16 ![]() |
Hey guys,
I've been reading through the Versteeg CFD book and I'm having some trouble understanding the boundary conditions. My question is, if you don't know the value of the pressure or velocity at the boundary then do you simply put zeroGradient or is it more complex than that? |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
fixedGradient bassically means that a specific value shall be assigned, for instance, no-slip condition which means U=0 at wall (Def. of fluid). Moreover, at the inlet usually a velcotiy is prescirbed.
zeroGradient means that the gradient of respective quantity is zero, meaning that the actual value is constant. Ususally, pressure is taken to be zeroGradient at a inlet, because the pressure is assumed to be constant there. Hopfully that will help |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Join Date: Dec 2009
Posts: 33
Rep Power: 16 ![]() |
idrama thanks. But I have another question.
What boundary conditions would you put for the outlet? Because we wouldn't have any idea of the velocity and pressure values. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 ![]() |
U=0 would be fixedValue, not fixedGradient! (Dirichlet vs Neumann).
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 ![]() |
Since you are not interested in the exact pressure value for an incompressible flow in a channel, but only pressure differences, you'd set the pressure to be 0 at the outlet and zeroGradient on the inlet.
|
|
![]() |
![]() |
![]() |
![]() |
#6 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
Firstly, you are right; i wrote under pressure (ain't got no time).
Actually, you can set every constant value you want at the outlet. Fact is, that in the NS equations only the alteration of the pressure is presented, e.g. dp/dx, which means that only for pressure differences are being solved. In postprocessing is it just easier to handle the pressure distribution about zero. Check it out: Make a easy pipe case with 100000 Pa and 0 Pa at the outlet. Note, my statements are only for incompressible flows. I've never dealt compressible or something else. Good luck |
|
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21 ![]() |
I have provided a fairly self-explanatory derivation[1] for why the zero gradient B/C is applied for pressure at walls. The same logic can be extended to inlets as well.
References: [1] http://www.cfd-online.com/Forums/ope...ssure-b-c.html |
|
![]() |
![]() |
![]() |
![]() |
#8 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
Hello,
I think there is something wrong! You wrote: u_A = 0 => du_A/dn = 0 [1] Should I have understood everything correctly, then is [1] wrong. Just consider the steady-state solution (you can find this form in every fluid mechanic textbook): u(y)=-1/(2*mu) dp/dx (h-y)y which has roots at 0 and h. The derivative goes as follows u'(y)=-1/(2*mu) dp/dx (h-2*y) which is for u'(0)=-1/(2*mu) dp/dx h != 0. All quatntites are non-zero, even dp/dx due to the pressure drop. Furthremore, du/dn = 0 implies you would have no wall shear stresses. Note, the presents of root of a function f at a postion t, i.e. f(t)=0, does not implie that the slope is zero. |
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21 ![]() |
I think du/dx = 0 and dv/dx =0 are vaild as both u and v at the bottom wall do not change in the X-direction. Also, in the derivation, I have NOT assumed du/dy=0, which means that the wall shear stress is free to exist.
|
|
![]() |
![]() |
![]() |
![]() |
#10 |
Senior Member
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21 ![]() |
Ah, I think I know what happened. You have confused my "x" for "n". Just for clarification, nowhere in the derivation have I ever assumed/stated that du/dn = 0 at the wall. There are no derivatives with respect to 'n' in my derivation. That's why I clearly marked out the co-ordinate system as X and Y.
Last edited by msrinath80; August 7, 2010 at 15:10. |
|
![]() |
![]() |
![]() |
![]() |
#11 | ||
Member
Join Date: Jul 2012
Posts: 31
Rep Power: 14 ![]() |
hey guys,
i think i have to revive this thread, because i dont quite understand the following statements: Firstly, Quote:
Quote:
Please correct me if I'm wrong, but zeroGradient for p at inlet doesn't seem to make sense to me. |
|||
![]() |
![]() |
![]() |
![]() |
#12 | |
New Member
XU XIAOJUN
Join Date: Feb 2018
Location: Tokyo
Posts: 2
Rep Power: 0 ![]() |
Quote:
Thanks for your explanation, It's very helpful. But I still have a question about the gradient. Dose zeroGradient means the gradient of respective quantity to the time is zero? or space? Since you said the actual value is a constant. I suppose it to be the time. Am I right? Jo
__________________
&追梦人& |
||
![]() |
![]() |
![]() |
![]() |
#13 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 ![]() ![]() |
Quote:
No. In space. In the direction normal to the patch. "Constant in time" would be fixedValue
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
![]() |
![]() |
![]() |
![]() |
#14 |
New Member
XU XIAOJUN
Join Date: Feb 2018
Location: Tokyo
Posts: 2
Rep Power: 0 ![]() |
Dear gschaider
Thanks for your quick reply. I got it. Thank you! Jo
__________________
&追梦人& |
|
![]() |
![]() |
![]() |
![]() |
#15 |
Senior Member
qutadah
Join Date: Jun 2021
Location: USA
Posts: 101
Rep Power: 5 ![]() |
and what about zero with time?
|
|
![]() |
![]() |
![]() |
![]() |
#16 |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 ![]() ![]() |
You mean a boundary condition whose value doesn't change over time? That is called fixedValue
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|
![]() |
![]() |
![]() |
![]() |
#17 |
New Member
Guo Zongchao
Join Date: Nov 2021
Posts: 8
Rep Power: 5 ![]() |
Hello, everyone.
I meet a trouble while setting pressure gradient. My inlet pressure gradient is time-varying, and there is a formula. FixedValue and fixedGradient are all useless for me. How could I write this inlet condition? I would appreciate it if someone could give me some ideas. |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
interFoam | andrea.pasquali | OpenFOAM | 8 | February 23, 2021 13:06 |
fixedGradient BC type with the channelOodles solver | Bedotto | OpenFOAM Running, Solving & CFD | 5 | April 15, 2019 10:49 |
Pressure instability with rhoSimpleFoam | daniel_mills | OpenFOAM Running, Solving & CFD | 44 | February 17, 2011 18:08 |
ZeroGradient pressure Outlet BC prob (adjustPhi.C) | Stylianos | OpenFOAM | 3 | March 23, 2010 13:28 |
zeroGradient pressure and moving mesh | kassiotis | OpenFOAM Running, Solving & CFD | 0 | July 22, 2009 10:28 |