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

fractional step method in cylindrical coordinates, help please?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 9, 2012, 13:31
Default fractional step method in cylindrical coordinates, help please?
  #1
Member
 
Peter
Join Date: Oct 2011
Posts: 52
Rep Power: 5
new_at_this is on a distinguished road
I am trying to solve developing internal flow in a circular duct using a staggered grid and cylindrical navier stokes equations.

At the moment I am using a fractional step method with primitive variables but my code is currently wrong and I think it has to do with the pressure correction step and I think I might be using the wrong equations. Please let me know if you see anything wrong. My process is listed below.

fractional step for advection terms
U^* = U^n-\Delta t(\frac{\partial UV}{\partial r}+\frac{\partial U^2}{\partial z})
V^* = V^n-\Delta t(\frac{\partial V^2}{\partial r}+\frac{\partial UV}{\partial z})

fractional step for diffusion terms
U^{**} = U^*+\frac{\Delta t}{Re}(\frac{1}{r}\frac{\partial U^{**}}{\partial r}+\frac{\partial^2 U^{**}}{\partial r^2} + \frac{\partial^2 U^{**}}{\partial z^2})
V^{**} = V^*+\frac{\Delta t}{Re}(\frac{1}{r}\frac{\partial V^{**}}{\partial r}+\frac{\partial^2 V^{**}}{\partial r^2} + \frac{\partial^2 V^{**}}{\partial z^2}-\frac{V^{**}}{r^2})

pressure correction
U^{n+1} = U^{**} - \Delta t(\frac{\partial P}{\partial z})
V^{n+1} = V^{**} - \Delta t(\frac{\partial P}{\partial r})

where Pressure is calculated by
\frac{1}{r}\frac{\partial P}{\partial r}+\frac{\partial^2 P}{\partial r^2}+\frac{\partial^2 P}{\partial z^2} = \Delta t(\frac{U^{**}}{r}+\frac{\partial U^{**}}{\partial r}+\frac{\partial V^{**}}{\partial z})
new_at_this is offline   Reply With Quote

Old   May 10, 2012, 16:00
Default
  #2
Member
 
Peter
Join Date: Oct 2011
Posts: 52
Rep Power: 5
new_at_this is on a distinguished road
so I'm pretty sure i've narrowed it down to be an error with the poisson pressure solver for cylindrical coordinates.

Is there anyone out there who is familiar marker and cell cfd in cylindrical? or with the mit18086_navierstokes code written by benjamin seibold? I modified his code to cylindrical coordinates.

Also is there anyone out there who might consider looking at my code? I have been staring at it for so long that I might be missing something because I am so used to it.

If someone knows of a simple test case that I can run to debug my poisson solver, please let me know.
new_at_this is offline   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
SIMPLE algorithm in 3D cylindrical coordinates zouchu Main CFD Forum 1 January 20, 2014 18:02
Pressure constant B.C. in fractional step method amin144 Main CFD Forum 10 April 10, 2012 17:58
Time Accuracy of the NITA - FRACTIONAL STEP method Paolo Lampitella FLUENT 1 March 19, 2008 06:45
Fractional Step, PRESTO! & LES Paolo Lampitella FLUENT 2 December 1, 2007 02:46
Multigrid method for cylindrical coordinates chinthakindi Main CFD Forum 0 October 30, 2003 05:00


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