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

fractional step method in cylindrical coordinates, help please?

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
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: 14
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

 


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
SIMPLE algorithm in 3D cylindrical coordinates zouchu Main CFD Forum 1 January 20, 2014 17: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 05:45
Fractional Step, PRESTO! & LES Paolo Lampitella FLUENT 2 December 1, 2007 01:46
Multigrid method for cylindrical coordinates chinthakindi Main CFD Forum 0 October 30, 2003 04:00


All times are GMT -4. The time now is 09:08.