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

Pressure correction

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

Like Tree2Likes
  • 1 Post By FMDenaro
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 23, 2015, 16:40
Default Pressure correction
  #1
New Member
 
Join Date: Dec 2015
Posts: 2
Rep Power: 0
Harr is on a distinguished road
Hi everybody,

I am not a cfd expert, I am writing a code that can solve the flow of a lid driven cavity. I am doing this in in-compressible steady-state, using the simple algorithm. When I inverse and multiply my pressure correction matrix with my source terms, I get a near singularity. My equations are:

P'=pressure correction.
Ap*P'p=Ae*P'e+Aw*P'w+An*P'n+As*P's+Su
Su=(density*u_w-density*u_e)*area+(density*v_s-density*v_n)*area

I do not know why I keep getting near singularity.
Any help would be appreciated.
Harr is offline   Reply With Quote

Old   December 23, 2015, 16:46
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,726
Rep Power: 60
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The pressure equation Div Grad p = q with proper BC.s leads to an algebric system that has indeed a singular matrix.
However, when the BC.s satisfies the compatibility relation, you get a solution apart from a constant value
Harr likes this.
FMDenaro is offline   Reply With Quote

Old   December 23, 2015, 17:00
Default
  #3
New Member
 
Join Date: Dec 2015
Posts: 2
Rep Power: 0
Harr is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
However, when the BC.s satisfies the compatibility relation, you get a solution apart from a constant value
Firstly thank you. Sorry I do not know to much about CFD yet, what do you mean by "when the BC.s satisfies the compatibility relation"?
Harr is offline   Reply With Quote

Old   December 23, 2015, 17:15
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,726
Rep Power: 60
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Harr View Post
Firstly thank you. Sorry I do not know to much about CFD yet, what do you mean by "when the BC.s satisfies the compatibility relation"?

Actually, this is not a CFD issue but a general mathematical aspect of linear algebra.
Consider the equation for pressure and integrate both side over the computational domain

Int [V] Div Grad p dV = Int [V] q dV

use Gauss

Int [S] n.Grad p dS = Int[V] q dV


The compatibility relation must be fulfilled by the boundary condition you prescribe on the frontier of the domain

n.Grad p = f

so that you have to fulfill the condition

Int [S] f dS = Int[V] q dV

The above relation must be satisfied in discrete sense after discretization.
Only this way, despite the singular matrix, you get one solution (apart a constant).

For further details you can see in literature, for example

https://books.google.it/books?id=lXB...matrix&f=false
FMDenaro is offline   Reply With Quote

Old   December 25, 2015, 18:49
Default
  #5
Member
 
Kaya Onur Dag
Join Date: Apr 2013
Posts: 94
Rep Power: 10
kaya is on a distinguished road
Quote:
Originally Posted by Harr View Post
Hi everybody,

I am not a cfd expert, I am writing a code that can solve the flow of a lid driven cavity. I am doing this in in-compressible steady-state, using the simple algorithm. When I inverse and multiply my pressure correction matrix with my source terms, I get a near singularity. My equations are:

P'=pressure correction.
Ap*P'p=Ae*P'e+Aw*P'w+An*P'n+As*P's+Su
Su=(density*u_w-density*u_e)*area+(density*v_s-density*v_n)*area

I do not know why I keep getting near singularity.
Any help would be appreciated.

I am very bad at these finite volume notations to me they all are gibberish so do excuse me

is this poisson equation, is [Ap] the laplace or grad.div. operator?
Quote:
Ap*P'p=Ae*P'e+Aw*P'w+An*P'n+As*P's+Su
Su=(density*u_w-density*u_e)*area+(density*v_s-density*v_n)*area
you may need to set a mean for your correction pressure such as setting 1 single value manually by modifying your matrix and your source (right hand side) otherwise to my knowledge your matrix system will be under determined.
kaya is offline   Reply With Quote

Old   December 26, 2015, 04:13
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,726
Rep Power: 60
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by kaya View Post
I am very bad at these finite volume notations to me they all are gibberish so do excuse me

is this poisson equation, is [Ap] the laplace or grad.div. operator?


you may need to set a mean for your correction pressure such as setting 1 single value manually by modifying your matrix and your source (right hand side) otherwise to my knowledge your matrix system will be under determined.

In a FV formulation, A*p should be the discretization of (1/|V|) Int [S] n.Grad p dS.
However, you could find also several researchers considering Div Grad p. For the pressure equation it is not suitable to discretize the laplacian operator.

The correction to the matrix entry and to the source term is consequent to the proper setting of the BC.s in such a way they satisfy the compatibility constraint
kaya likes this.
FMDenaro is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Periodic flow using Cyclic - comparison with Fluent nusivares OpenFOAM Running, Solving & CFD 30 December 12, 2017 05:35
Fluidized Bed: Error: Divergence detected in AMG solver: pressure correction Error Ob Mole89 FLUENT 5 April 12, 2014 09:32
Timestep and Pressure Correction Relationship in SIMPLE rks171 Main CFD Forum 23 May 4, 2012 01:04
changing the coefficients of pressure correction Noel Phoenics 1 April 7, 2009 08:54
Does star cd takes reference pressure? monica Siemens 1 April 19, 2007 11:26


All times are GMT -4. The time now is 23:27.