CFD Online URL
[Sponsors]
Home > Forums > Main CFD Forum

Pressure correction in cylindrical coordinates

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 12, 2012, 02:47
Default Pressure correction in cylindrical coordinates
  #1
Member
 
Peter
Join Date: Oct 2011
Posts: 52
Rep Power: 5
new_at_this is on a distinguished road
I am trying to write a code for the equivalent of a lid driven cavity in cylindrical coordinates. As a reference, I have adapted a 2D cartesian lid driven cavity code written by Seibold which can be found at: http://www-math.mit.edu/cse/codes/mi...navierstokes.m

The non linear terms are calculated using an explicit upwind method and the viscosity terms are implicit. I think I am having problems with the pressure corrections step in cylindrical coordinates.

Right now, I am finding P by solving the following equation where ** denotes the velocity field after a partial step for the nonlinear and viscosity terms.:

\frac{1}{r}\frac{\partial P}{\partial r} + \frac{\partial^2 P}{\partial r^2}+\frac{\partial^2 P}{\partial z^2} = \frac{u_r^{**}}{r}+\frac{\partial u_r^{**}}{\partial r} +  \frac{\partial u_z^{**}}{\partial z}

Could some one give me some tips? My u_z velocity profiles look like this and I expect the profile at 0.5L to be similar to a poiseuille profile.

http://s1118.photobucket.com/albums/...nt=axiprof.jpg
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
Divergence detected in AMG solver: pressure correction thisisit FLUENT 7 April 1, 2014 07:54
Pressure Outlet Guage pressure Mohsin FLUENT 34 January 12, 2012 06:45
Help! Divergence detected in AMG solver: pressure correction user0314 FLUENT 3 September 2, 2011 17:12
The correction on pressure equation of SIMPLE algorithm in MRFSimpleFOAM solver renyun0511 OpenFOAM Running, Solving & CFD 0 November 10, 2010 02:47
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13


All times are GMT -4. The time now is 05:00.