CFD Online Logo CFD Online URL
Home > Forums > Phoenics

How to realize fixed pressure in my codes?

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

LinkBack Thread Tools Display Modes
Old   March 13, 2003, 17:13
Default How to realize fixed pressure in my codes?
Posts: n/a

I am coding myself and using the method of Phoenics to deal with the bcs: source S=[cphi+cm(vm-p)][vphi-phi]

to fix outlet pressure:in phoenics it is suggested to make cm=1000*expected flow rate/T, vm=0.0 then the incell pressure can be set to fix pressure of 1e-3 order of the vm(actually zero here).

Should I put this mass source into P' equation then I will get incell pressure fixed?

Another question is: when outflow rate is known for example 1m3/s how can I introduce this "source" in the momentum equation by using S=[cphi+cm(vm-p)][vphi-phi]?

When outflow occurs the source S changed into S=cphi(vphi-phi) only and usally the cphi is zero so does that mean no source will be needed to add into the momentum equation for outlet cells?

Thanks a lot for your help and patience.

  Reply With Quote

Old   March 14, 2003, 06:26
Default Re: How to realize fixed pressure in my codes?
Posts: n/a
The expression dSm/dP should go into the ap (central coefficient) of the pressure-correction (p') equation, where Sm=cm*(vm-Pp)*T and T is is the source type, i.e. area, volume or cell(unity). Then a large value of Cm will act to suppress pressure corrections and so fix the pressure in the cell.

Yes, if there is outflow from the cell for the variable phi then nothing need be added to the phi equation. This follows from use of inflows=outflows (continuity) when deriving the finite-volume form of the equation for phi. This form ensures that -m*phi is taken from the cell when there is outflow.
  Reply With Quote

Old   March 16, 2003, 01:23
Default Re: How to realize fixed pressure in my codes?
Posts: n/a
Hi,cham SUPPORT:

Thanks so much!I have tried this for duct flow and it works pretty well.

But: When the pressure is "fixed", small changes resulting from dynamic- head variations can change inflow to outflow and vice versa, with consequent convergence difficulties.

How phoenics deal with this problem?Any magical tricks?(this problem is significant when fixed pressure is applied to large openings where inflow and outflow both exist.) During every iteration when incell pressure is smaller than the outside pressure then this cell should be taken as an inflow cell BC for this time iteration,right? To make the question clearer: For fixed pressure case how to deal with the situation where both inflow and outflow occurs in single opening in Phoenics?

Thanks again for your timely answer.

  Reply With Quote


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
how to realize the following data input? Steven CFX 4 March 24, 2008 16:53
Help: how to realize UDF on parallel cluster? Haoyin FLUENT 1 August 6, 2007 13:53
Has the pressure to be fixed at one node? Steve Song Main CFD Forum 5 February 9, 2003 14:14
fixed pressure alex Phoenics 6 October 6, 2001 11:29
fixed pressure john Main CFD Forum 4 August 6, 1999 00:34

All times are GMT -4. The time now is 08:45.