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

pressure outlet BC for incompressible flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 10, 2005, 06:23
Default pressure outlet BC for incompressible flow
  #1
khaiching
Guest
 
Posts: n/a
Hi: I have implemented the SIMPLE method by Patankar (1980) to compute incompressible flow. To correct the pressure field, the so-called Pressure-Correction Equation has been employed.

I am in the midst of searching the correct implementation of BC at the boundaries (inlet, outlet, wall) for pressure correction terms (p'). At outlet faces, Constant pressure (p) has been used and thus the p' is set to 0 while extrapolating the 3 velocity components; for other boundaries, p and p' are extrapolated.. Unfortunately, the solution diverges !!

Is my treatment correctly implemented? Could anyone help to advise on the correct implementation of pressure outlet BCs (fixing the P) for SIMPLE algorithm??

I have tried to scale the mass flow rate at outlet boundary. It works; however, this method needs prior knowledge of the outlet mass flow rate; which I think it is impractical for unsteady flow..What you guys think??

-khai ching-
  Reply With Quote

Old   October 10, 2005, 06:41
Default Re: pressure outlet BC for incompressible flow
  #2
Quarkz
Guest
 
Posts: n/a
hi,

maybe u can try using under-relaxation mtds for your pressure & maybe velocities as well. this should prevent divergence. recommended values from fluent are 0.3 & 0.7 respectively. hope it helps
  Reply With Quote

Old   October 10, 2005, 06:49
Default Re: pressure outlet BC for incompressible flow
  #3
khaiching
Guest
 
Posts: n/a
Thanks for the reply..

I have implemented 2 methods to treat the outlet BCs

(1) Pressure outlet: fixing the pressure (p' = 0) and extrapolate the velocities.. I have tried to reduce the u/relaxation factor to about 0.05 for pressure and 0.10 for velocities; However, it fails...

(2) mass scaling at outlet: scale the velocity vectors at outlet as accordance to the expected mass flow rate.. This works..I have extrapolated p' and P at all boundaries..

I am wondering whether the Pressure Correction Equation can be implemented together with the pressure outlet BC.. From Davidson (1996), he has extrapolated p' and P at ALL boundaries..

-khai ching-
  Reply With Quote

Old   October 10, 2005, 16:00
Default Re: pressure outlet BC for incompressible flow
  #4
Adrin Gharakhani
Guest
 
Posts: n/a
>> I have tried to scale the mass flow rate at outlet boundary. It works; however, this method needs prior knowledge of the outlet mass flow rate; which I think it is impractical for unsteady flow

If I understand you correctly, you extrapolate your velocities to the boundaries and then scale them using conservation of mass. The extrapolation is probably ok (there really is no good exit boundary condition, without having to make some sort of simplification/assumption). My response is particularly related to the scaling issue. There is nothing wrong with scaling, and indeed you may have to do it anyway to make sure you conserve mass. Your concern about the need to know the mass flow rate in advance seems unjustified since you're dealing with incompressible flow. Therefore, unsteady or not, you already know the exist mass flow rate based on your inlet (and wall) boundary conditions.

Adrin Gharakhani
  Reply With Quote

Old   October 13, 2005, 23:48
Default Re: pressure outlet BC for incompressible flow
  #5
khaiching
Guest
 
Posts: n/a
Dear Adrin: Does that mean I must scale the mass flow rate at outlet even I impose the constant Pressure BC at outlet?

I am using the Pressure Equation right now instead of the pressure correction equations. I am still in the midst of verifying my implementation..

Please advise.

-khai ching-
  Reply With Quote

Old   October 14, 2005, 15:42
Default Re: pressure outlet BC for incompressible flow
  #6
Adrin Gharakhani
Guest
 
Posts: n/a
Just want to make myself clear what I mean by scaling: you have an input mass flow rate m_i, and you get a solution that has an exit velocity distribution u_e. Now if you multiply u_e by the local density and cross-sectional area, chances are you'll get an exit mass flow rate m_e, which is not exactly equal to m_i. (hopefully the two are not too far off; if so, then you have an error some place). Then you want to scale your u_e such that m_e=m_i. This "scaling" has nothing to do with your BC type; it is there to make sure that you don't have source accumulation/depletion (by violating continuity)

I don't know the details of your approach, but make sure you don't over/under-specify your boundary conditions (i.e., have more BCs imposed than is necessary as dictated by the pde)

Adrin Gharakhani
  Reply With Quote

Old   October 15, 2005, 03:58
Default Re: pressure outlet BC for incompressible flow
  #7
khaiching
Guest
 
Posts: n/a
thanks..

My code works right now with the pressure outlet BC implemented and continuity is satisfied using the Pressure Equation..(Thanks to Dr Hrvoje Jasak)

I still have some doubts about the implementation of Pressure Correction Equations(PCE). Anyone has used the PCE in concert with the constant pressure outlet BC? What will be the boundary conditions (p') at the boundary faces? So far, I have extrapolated p' at all boundaries while fixing p'=0 at outlet(since I fix P_out). This method dint work unless I extrapolated p' and P at outlet as well. By doing that, I have to scale the mass flow rate at outlet to ensure global mass conservation..

-khai ching-
  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
reverse flow in pressure outlet BC shahzeb irfan Main CFD Forum 2 August 10, 2011 00:40
Disturbed flow field at outlet boundary (Multiphase flow through pipe) Michiel CFX 17 April 21, 2010 11:14
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 22:31
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 20:41.