|
[Sponsors] |
May 26, 2011, 07:36 |
Implicit elliptic pressure field equation.
|
#1 |
New Member
Tom
Join Date: Mar 2011
Posts: 12
Rep Power: 15 |
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. |
|
|
|
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 12:21 |
How to add correcting pressure equation in SIMPLE of MRFSimpleFOAM? | renyun0511 | OpenFOAM Programming & Development | 0 | November 4, 2010 01: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 |