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

Projection Methods

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By Petrov
  • 2 Post By valgrinda

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 8, 2011, 23:25
Default Projection Methods
  #1
New Member
 
Felipe Hernandez
Join Date: Apr 2011
Posts: 13
Rep Power: 15
Petrov is on a distinguished road
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?
rajibroy likes this.

Last edited by Petrov; June 8, 2011 at 23:53.
Petrov is offline   Reply With Quote

Old   June 10, 2011, 04:32
Default pjm
  #2
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
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
rajibroy and mostanad like this.
valgrinda is offline   Reply With Quote

Old   June 10, 2011, 10:50
Default
  #3
New Member
 
Felipe Hernandez
Join Date: Apr 2011
Posts: 13
Rep Power: 15
Petrov is on a distinguished road
Wow, those are great resources!

Thanks for your help.
Petrov is offline   Reply With Quote

Reply


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
[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


All times are GMT -4. The time now is 04:07.