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

NS on polar coordinate grid

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 31, 2005, 23:32
Default NS on polar coordinate grid
  #1
Brian
Guest
 
Posts: n/a
Hi,

I am trying to solve the incompressible NS equations on a polar coordinate grid. I'm using a modified "stable fluid" method and am getting stable but strange results. After the pressure correction step the fluid flow is biased towards the angular direction and flows that should go horizontal are being curved. There is no divergence, however, and without the pressure correction the flow is as expected, but obviously not divergence free. Therefore I believe the pressure correction step is not maintaining conservation of momentum correctly. I have looked into coordinate transformations to set up the pressure poisson equation but the results are the same as the method I'm using now (divergence of the pressure gradient).

Here's a rundown of the method:

U0 = U^(n-1)

U1 = U0 + dt*F // gravity, body forces

U2 = transport(U1) // convection term

U3 = U2 - (dt/rho)*del(p) // continuity

U^n = U3

where pressure is found by solving this poisson equation:

(dt/rho)del^2(p) = divergence(U2)

and the del^2(p) matrix is set up as the divergence of the pressure gradient

For the CFD experienced, does this sound like a resonable approach? I don't think I have everything quite right.

Thanks in advance, Brian
  Reply With Quote

Old   June 2, 2005, 05:35
Default Re: NS on polar coordinate grid
  #2
Alexander Starostin
Guest
 
Posts: n/a
I have not noticed a small detail maybe because it is undoubtely right in your scheme:

For the polar coordinate system the additional terms come to NS equation due to the curvature of the coordinates.
  Reply With Quote

Old   June 3, 2005, 20:39
Default Re: NS on polar coordinate grid
  #3
Brian
Guest
 
Posts: n/a
You're right, there are additional terms. However, the pressure correction step guarantees a divergence free field and not momentum. I'm wondering if I need to either find a momentum correction term or somehow tie this in when solving for pressure so the pressure gradient respects momentum?
  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
MapFields to New Grid For Extreme Grid Deformations due to Body Motion albcem OpenFOAM 0 May 5, 2009 14:17
Is polar coordinate available in Gambit? David FLUENT 1 July 20, 2007 16:08
GRID TO GRID INTERPOLATION in FLUENT calogero FLUENT 3 June 4, 2003 08:32
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 05:59
Troubles modelling flow through a grid Hans Klaufus CFX 1 June 28, 2000 16:43


All times are GMT -4. The time now is 10:34.