# cylindrical coordinates at thita=0

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

 April 30, 2006, 20:27 cylindrical coordinates at thita=0 #1 George Guest   Posts: n/a Hi everybody! I am writing a code for 3D unsteady,turbulent,compressible flow in cylindrical coordinates. The problem that I am facing for many months (!!!) is that at thita=0 and 360 degrees the results that I am having are not the same (or not almost the same), when the variables are small. Actually they are very different. The bigest problem is with the pressure correction, because the values that it takes are quite small (10E-3). Beacuse of this fact the velocity correction at 0 and 360 degrees is very different and also the d(P)/d(theta) is very different (for example: at 0 degrees: 0.001 and at 360 degrees: -0.0005). These values are when the code has converged. During the early iterations the difference is very large (5 and -12). The way that I am implemending this, is that I use ghost cells at 0 (I=1) and 360 (I=NI) degrees that take the values of I=NI-1 and I=2 respectively. These ghost cells are actually the boundary cells and have exactly tha same characteristics (volume, dimensions etc.) Do you find this way correct,or is there a better one? Has anybody faced the same problem with cylindrical or spherical coordinates? I have also tried periodic or cyclic TDMA but the results were exactly the same. Any thoughts? Thank you in advance!

 May 1, 2006, 04:16 Re: cylindrical coordinates at thita=0 #2 ramp Guest   Posts: n/a You might try to have one more ghost cells at theta=360 and go for I=0 (theta=0), NI+1 (theta=360+2*dtheta) and then implement the conditions as NcellI=NI+2 phi(k,j,NcellI)=0.5*(phi(k,j,2)+phi(k,j,NcellI-1)) I am using 3-D FVM in cartesian system as a 2-D solver and its working for me for 2-D flow over cylinder.

 May 1, 2006, 10:45 Re: cylindrical coordinates at thita=0 #3 George Guest   Posts: n/a Now I have one ghost cell at 0 degrees and one at 360 degrees. You suggest me to have two ghost cells at 0 degrees and two at 360 degrees? The solver (TDMA) for the cylindrical coordinates in comparison to the solver for cartesian coordinates will change or will be the same? The flow that I am simulating is inside a cylinder and not around a cylinder. Waiting for your answer! Thank you

 May 1, 2006, 15:04 Re: cylindrical coordinates at thita=0 #4 ramp Guest   Posts: n/a I solved for external flow with the two ghost cells.... you might try it too internal also... I guess there wont be any change in the solver but only implemenetation of the condition will change.... Best regards.... ramp

 May 2, 2006, 06:40 Re: cylindrical coordinates at thita=0 #5 George Guest   Posts: n/a I have sent you my solver to see if it is the same with yours. If not, can you send me the solver with the changes for the cyclic conditions and its implementation that you have? Thank you

 May 2, 2006, 18:24 Re: cylindrical coordinates at thita=0 #6 George Guest   Posts: n/a Sorry, I forgot. My mail is gkosmadakis@gmail.com

 May 4, 2006, 08:40 Re: cylindrical coordinates at thita=0 #7 yousef Guest   Posts: n/a hi george i work on the problem of laminar mixed convection in the entrance region of a semicircular duct. so my problem is three dimenstinal and i work with SIMPLER AND FVM i do not understand your equation and i will asked to my supervisor. can i ask you george in which is your university? because your problem is exating and i really want work with turbelent and compressible flows. really if tell me i will try to join your reaserch. yours yousef

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Randheer Yadav FLUENT 2 January 18, 2015 17:25 zouchu Main CFD Forum 1 January 20, 2014 18:02 falopsy Main CFD Forum 44 January 15, 2014 03:17 T.D. OpenFOAM Running, Solving & CFD 17 September 26, 2010 15:37 Franz Wingelhofer CFX 0 December 28, 1999 08:46

All times are GMT -4. The time now is 13:57.