
[Sponsors] 
October 12, 2010, 09:48 
SIMPLE pressure correction on unstructured mesh

#1 
New Member
abcdef123
Join Date: Dec 2009
Posts: 17
Rep Power: 10 
Hi, could anyone provide some details or references on the derivation of this equation using the finite volume method for discretization?
Thanks in advance. 

October 14, 2010, 09:25 

#2 
New Member
abcdef123
Join Date: Dec 2009
Posts: 17
Rep Power: 10 
So I can get as far as this,
The continuity equation is defined as and using Gauss' Theorem we can say that Approximating this integral for each control volume surface we get and the SIMPLE approximations for the velocity correction (collocated grid) is where and p is the pressure. Normally (on a structured mesh) we would be able to approximate the pressure gradient term to get an equation with coefficients on each neighboring point. The problem is how is this done for an unstructured mesh? I suspect the thing I am missing is a pressure gradient approximation that would allow be to construct the pressure correction equation. Any ideas? Thanks again in advance. Last edited by abcdef123; October 14, 2010 at 09:40. Reason: fixed latex equations 

October 14, 2010, 09:37 

#3 
New Member
abcdef123
Join Date: Dec 2009
Posts: 17
Rep Power: 10 
I am also having trouble getting the latex equations to show up in my post.. sorry


October 14, 2010, 14:00 

#4 
Member
ganesh
Join Date: Mar 2009
Posts: 40
Rep Power: 10 
Dear abcdef,
Well, I maynot be able to directly help you out in writing down the pressure correction equation for SIMPLE, but I can point out how it is done on unstructured meshes. Please note that I am not writing down equations here either, so kindly excuse me. The basic steps are as follows. 1. The momentum equation gives you an auxiliary velocity, call it u* 2. Thsi auxiliary velocity doesnot satisfy continuity and this leads us to a pressure correction equation. 3. The PCE is of the form div u* = div(grad(dp))=Lap(dp) The idea then is to integrate this equation over a control volume and apply the divergence theorem on both sides. The LHS is merely \sum(u*_f s_f) where u*_f is the auxiliary velocity at the face (the velocities are at the face centroids in a staggered approach; if a collocated approach is used these are interpolated from the centroidal velocities with a suitable correction to avoid pressurevelocity decoupling). This term acts as the source because the auxiliary velocity is known by solving the momentum equation. The RHS leads to the normal pressure gradient, which can be approximated in different ways. The way that I use, and would recommend (although there exist other approaches) is to use an orthogonal correction to the normal pressure gradient; which consists of two parts: one which is obtained using a finite difference of pressures in cells sharing the face and the other which depends on the gradients of pressures in these cells. For simplicity, in a collocated approach, as a first try, you can neglect the second term. The result is an implicit sytem of equations for the pressures that can be solved using any suitable linear solver. For more details, you could refer to the following paper (and references therein) or also google for related sources. A numerical method for largeeddy simulation in complex geometries, Mahesh et.al., JCP 2005. Hope this helps. Regards, Ganesh 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
2D Mesh Generation Tutorial for GMSH  aeroslacker  Open Source Meshers: Gmsh, Netgen, CGNS, ...  12  January 19, 2012 04:52 
Neumann pressure BC and velocity field  Antech  Main CFD Forum  0  April 25, 2006 02:15 
BCs for Pressure Correction Equation (SIMPLE)  Bharath Somayaji  Main CFD Forum  1  March 1, 2006 07:12 
what the result is negatif pressure at inlet  chong chee nan  FLUENT  0  December 29, 2001 06:13 
unstructured vs. structured grids  Frank Muldoon  Main CFD Forum  1  January 5, 1999 11:09 