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

Implicit elliptic pressure field equation.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 26, 2011, 07:36
Default Implicit elliptic pressure field equation.
  #1
New Member
 
Tom
Join Date: Mar 2011
Posts: 12
Rep Power: 6
TheBoyce is on a distinguished road
I am trying to model stratified spin-up from rest. My pressure update equation is implicit and therefore the result is large set of linear equations to be solved. The structure of the coefficient matrix is tridiagonal with fringes which is a very sparse matrix. To solve this system of equations I am using the iterative biconjugate gradient method. My problem is that in order to get a solution an arbitrary point needs to be set to an arbitrary value, however, the point at which I define the pressure is acting as a sink. The value of the pressure I set is simply added to the entire pressure field and the pressure gradients are not changed. The location of where I set the pressure does make a difference to the pressure gradients and therefore alters the strength of the sink as it is the pressure gradients that are used in the velocity update equations. Any ideas as to why this sink is occurring and any ideas as to how I can remove it would be greatly appreciated.

Thanks in advance.
TheBoyce is offline   Reply With Quote

Old   June 3, 2011, 10:12
Default
  #2
New Member
 
Tom
Join Date: Mar 2011
Posts: 12
Rep Power: 6
TheBoyce is on a distinguished road
Any ideas? Please? I'm desperate.
TheBoyce is offline   Reply With Quote

Old   June 3, 2011, 13:08
Default pressure sink
  #3
Senior Member
 
CFDtoy
Join Date: Mar 2009
Location: United States
Posts: 145
Blog Entries: 2
Rep Power: 8
CFDtoy is on a distinguished road
If I were you, I will look closely at my pressure correction eqn and see how the af terms (mass flux term) is modeled. Be careful implementing your rhie-chow interpolation and since you have body force, make sure you are absorbing rho*g within p, pstatic should be treated as effectively pcell-rho*g...body force involved.

I think your treatment of the pressure for evaluating face vel will be crucial..take a look at that.

good luk.

CFDtoy

Quote:
Originally Posted by TheBoyce View Post
Any ideas? Please? I'm desperate.
__________________
CFDtoy
CFDtoy 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
Pressure work term in turbulent K.E. equation lost.identity Main CFD Forum 0 March 8, 2011 13:21
How to add correcting pressure equation in SIMPLE of MRFSimpleFOAM? renyun0511 OpenFOAM Programming & Development 0 November 4, 2010 02:38
FEM pressure poisson equation (Implicit) cfd-beginner Main CFD Forum 0 August 9, 2005 13:32
pressure poisson equation on staggerred grid chandra Main CFD Forum 9 May 5, 2005 06:01
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02


All times are GMT -4. The time now is 05:11.