|
[Sponsors] |
June 8, 2011, 23:25 |
Projection Methods
|
#1 |
New Member
Felipe Hernandez
Join Date: Apr 2011
Posts: 13
Rep Power: 15 |
I've got a incompressible Navier-Stokes code that works reasonably well (compared to the taylor-green vortex) that uses a pressure-correction time integration method. I want to implement a projection method now, both as a learning experience and to compare the accuracies.
I'm just trying to get a gist of how it works and how it compares to the pressure correction method. From what I understand, I get an intermediate velocity (probably from the old pressure) and then I project the velocity onto divergence-free vector field (by projecting the vector onto the null-space of another matrix?!), then I solve for pressure, and then I make it divergence-free again. Is this correct? If not, could you describe projection methods to me? Are there any good (and free?!) articles you could point me to? I have a strict $0 budget for learning CFD. EDIT: I also want the method to be accurate at the no-slip boundary conditions so that I can calculate the drag/lift forces. I was reading that projection methods may have a problem with this. Should I look at other options? Last edited by Petrov; June 8, 2011 at 23:53. |
|
June 10, 2011, 04:32 |
pjm
|
#2 |
Super Moderator
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19 |
Hi,
you got the projection almost right, except that the intermediate velocity is calculated with the pressure not included. This is the major difference to SIMPLE type methods: There you include the pressure gradient in the Navier-Stokes equation, so you only need to solve the poisson equation for the pressure correction. The procedure for the projection method is (with explicit time treatment for the velocities): - calculate intermediate velocity u* without the pressure gradient - calculate the divergence of the velocity field, this term will be the right hand side of the poisson equation for the pressure - discretize the poisson equation for the pressure (trivial for constant density) - solve the poisson equation iteratively - correct the velocities (explicitly) with the pressure gradient - go to next timestep Here is a link to a paper, where the projection method is used and explained with some detail: http://wissrech.ins.uni-bonn.de/rese.../croce/lsm.pdf Here is book, where the discretization of the projection method is described with even more detail: http://books.google.com/books?id=ehu...page&q&f=false Cheers Hans |
|
June 10, 2011, 10:50 |
|
#3 |
New Member
Felipe Hernandez
Join Date: Apr 2011
Posts: 13
Rep Power: 15 |
Wow, those are great resources!
Thanks for your help. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 08:36 |
projection method | sudhakar | Main CFD Forum | 4 | December 8, 2007 02:25 |
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM | kenn | Main CFD Forum | 2 | July 18, 2004 19:28 |
free surface flow and projection methods | sylvain | Main CFD Forum | 2 | April 26, 2002 10:45 |
Projection methods | ZY Sun | Main CFD Forum | 9 | September 15, 2000 19:00 |