CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

drichlet boundary condition for pressure

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 26, 2012, 21:11
Default drichlet boundary condition for pressure
  #1
Senior Member
 
Mehdi
Join Date: Jan 2011
Location: Iran
Posts: 128
Rep Power: 5
mb.pejvak is on a distinguished road
I plan to apply constant pressure boundary condition in outlet in collocated arrangement. unfortunately, I could not find any book or paper that explain how to apply it. does any one know how I do it or introduce a reference to me.

thanks in advance
__________________
Best Regards;
Mehdi
E-mail: mb.pejvak[at]Gmail[dot]com
mb.pejvak is offline   Reply With Quote

Old   May 26, 2012, 23:41
Default
  #2
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 386
Rep Power: 5
cdegroot is on a distinguished road
At volumes adjacent to the outlet, just take the integration point value of the pressure on the outlet faces to be equal to the desired pressure. Seems too simple. Is there some specific issue with this BC you are struggling with?
cdegroot is offline   Reply With Quote

Old   May 27, 2012, 13:11
Default
  #3
Senior Member
 
Mehdi
Join Date: Jan 2011
Location: Iran
Posts: 128
Rep Power: 5
mb.pejvak is on a distinguished road
thank you for your help but in simpler method, when momentum equation is descritized, value of pressure in cells appear and not in faces. so I guess I should assume a guest cells and value of these cells set constant pressure that I wanted to apply in outlet. but problem is I should have mass flux and velocity of these cells, because they are required for momentum equation. and I don't know how to deal with.
I hope I can explain the problem good.
__________________
Best Regards;
Mehdi
E-mail: mb.pejvak[at]Gmail[dot]com
mb.pejvak is offline   Reply With Quote

Old   May 28, 2012, 09:14
Default
  #4
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 386
Rep Power: 5
cdegroot is on a distinguished road
Yes ghost cells is one way to do it. Blazek's book "Computational Fluid Dynamics" has a detailed description of this approach. So does Versteeg and Malalasekera.

In my code I simply place a "volume" with its centroid at the integration point and calculate any quantities adjacent to the boundaries implicitly using this value. Then I apply a condition to the boundary volume and absorb it into the interior. It is similar to the ghost cell approach, but I place the centroid of ghost cell right at the integration point and the cell doesn't have any volume. It is more like a ghost point I guess. For example at the volume adjacent to the boundary you will have an equation like

AP*UP = AW*UW + AE*UE +B

in one dimension where 'E' happens to be a boundary volume. Then at the boundary volume, for a Dirichlet condition, you can write

AP*UP = B.

where AP=1 and B is your specified value. Then the second equation can be absorbed into the first. Hopefully this made sense and applies to your situation.
cdegroot is offline   Reply With Quote

Old   May 31, 2012, 20:15
Default
  #5
Senior Member
 
Mehdi
Join Date: Jan 2011
Location: Iran
Posts: 128
Rep Power: 5
mb.pejvak is on a distinguished road
Dear cdegroot
thank you for your answer. in fact, I think I could not find it out how to apply it. in collocated arrangement, for applying simpler method we have P and p' (correction part). the equation in two step with applying mometum interpolation (Rhie and Choe) in continuity equation:
ap*p=Aw*Pw+AE*Pe+B ; b=Uf*dy-vf*dx+ ...
ap*p'=Aw*P'w+AE*P'e+B ; b=Uf*dy-vf*dx+ ...
when we want to apply constant pressure in boundary, P'=0 and P=cont. so we will (assume E is guest cell):
ap*p=Aw*Pw+AE*C+B ; b=Uf*dy-vf*dx+ ... ; C=value of pressure in boundary
ap*p'=Aw*P'w+B ; b=Uf*dy-vf*dx+ ...
but we should calculate Ae in guest cell P equation that need geometric and flow character of the cell. in fact I dont know how to deal with this part.
__________________
Best Regards;
Mehdi
E-mail: mb.pejvak[at]Gmail[dot]com

Last edited by mb.pejvak; May 31, 2012 at 20:57.
mb.pejvak is offline   Reply With Quote

Old   June 1, 2012, 00:35
Default
  #6
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 386
Rep Power: 5
cdegroot is on a distinguished road
Quote:
Originally Posted by mb.pejvak View Post
Dear cdegroot
but we should calculate Ae in guest cell P equation that need geometric and flow character of the cell. in fact I dont know how to deal with this part.
If P is a ghost cell and it has a Dirichlet condition, it does not depend on neighbouring values, so Ae=Aw=0 in the ghost cell, true?

Last edited by cdegroot; June 1, 2012 at 00:37. Reason: Mistake
cdegroot is offline   Reply With Quote

Old   June 1, 2012, 00:41
Default
  #7
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 386
Rep Power: 5
cdegroot is on a distinguished road
With regards to my last comment, it is possible I misunderstood you. Are you talking about the equation in the ghost cell or the equation at a volume adjacent to the boundary and you are having trouble figuring out AE in that equation?
cdegroot is offline   Reply With Quote

Old   June 1, 2012, 03:26
Default
  #8
Member
 
Join Date: May 2012
Posts: 60
Rep Power: 4
Simbelmynė is on a distinguished road
I might have a similar problem and I have attached a working code in this thread:

Boundary condition, incompressible colocated code

Perhaps we can help each other?
Simbelmynė is offline   Reply With Quote

Old   June 1, 2012, 04:43
Default
  #9
Senior Member
 
Mehdi
Join Date: Jan 2011
Location: Iran
Posts: 128
Rep Power: 5
mb.pejvak is on a distinguished road
dear cdegroot
I think I can convey my purpose well. in fact my problem is applying Dirichlet pressure boundary condition and for velocities I applied Dirichlet boundary by set face of boundary cells as but in pressure equation we need to add a guest cell to apply constant pressure condition in outlet. so in continuity equation in simpler method, we should have AE in collocated arrangement. so I don't know how to find it?
If i made mistake about that I'd really appreciate it if you correct me.
__________________
Best Regards;
Mehdi
E-mail: mb.pejvak[at]Gmail[dot]com
mb.pejvak is offline   Reply With Quote

Old   June 1, 2012, 04:49
Default
  #10
Senior Member
 
Mehdi
Join Date: Jan 2011
Location: Iran
Posts: 128
Rep Power: 5
mb.pejvak is on a distinguished road
dear Simbelmynė
I'll be happy if I can help you. but in fact, I have not worked FDM. and I don't know I can help you. I think Hoffman book is based on FDM. Have you ever taken look at it?
__________________
Best Regards;
Mehdi
E-mail: mb.pejvak[at]Gmail[dot]com
mb.pejvak is offline   Reply With Quote

Old   June 1, 2012, 05:02
Default
  #11
Member
 
Join Date: May 2012
Posts: 60
Rep Power: 4
Simbelmynė is on a distinguished road
Dear mb.pejvak,

If you are using a structured mesh then I believe the two approaches are similar enough. I have a ton of different books where my favorites are the ones by C.A.J. Fletcher, and by Ferziger and Peric. I have managed to write working codes but I have some doubts about my boundary conditions (similar to what you describe, although I have solved it by using extrapolation from interior).

Ferziger and Peric suggests using extrapolation of pressure at boundaries if a colocated setup is used where the control volumes extend to the boundary.
Simbelmynė is offline   Reply With Quote

Old   June 1, 2012, 15:41
Default
  #12
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 386
Rep Power: 5
cdegroot is on a distinguished road
Quote:
Originally Posted by mb.pejvak View Post
dear cdegroot
I think I can convey my purpose well. in fact my problem is applying Dirichlet pressure boundary condition and for velocities I applied Dirichlet boundary by set face of boundary cells as but in pressure equation we need to add a guest cell to apply constant pressure condition in outlet. so in continuity equation in simpler method, we should have AE in collocated arrangement. so I don't know how to find it?
If i made mistake about that I'd really appreciate it if you correct me.
I think will have to bow out from this thread at this point. While I have used SIMPLE before I now use a fully coupled (not segregated) approach with unstructured grids. I don't recall any specific difficulties with pressure bcs using SIMPLE. Unfortunately I don't think I can offer any further assistance at this point. I hope you are able to figure out your issue.
cdegroot is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Can anyone give me some hint on how to make traction free boundary condition? poplar OpenFOAM 2 May 16, 2013 03:02
Boundary condition & reverse flows zhenglun.wei FLUENT 4 October 25, 2011 17:38
Setting outlet Pressure boundary condition using CAFFA code Mukund Pondkule Main CFD Forum 0 March 16, 2011 04:23
Domain Imbalance HMR CFX 3 March 6, 2011 20:10
How to set boundary condition in Fluent for the fo Peiyong FLUENT 1 November 10, 2006 11:44


All times are GMT -4. The time now is 09:38.