
[Sponsors] 
June 19, 2011, 00:10 
Turbulence model for flow over a cylinder at Re=10000

#1 
New Member
Mayank Ojha
Join Date: May 2011
Posts: 22
Rep Power: 7 
Hello,
Am trying to find the Lift coeff, Drag coeff and the Strouhal number for a cylinder in a cross flow at Re=10,000 Things I tried: 1. Created a mesh and ran cases for tstep=0.001, 0.0001 & 0.00001 using KW and SST turbulence model. The y+ for the mesh was around 3.54. So now I created another mesh. 2. For the new mesh, the y+ value over the cylinder surface was found to be 0.62. I ran with SST and KW at tstep=0.001 which is 2% of vortex shedding time. The time step is low enough to capture the vortex shedding. Results: 1. My St number is 0.22 while the expected value from DNS and experimentsis 0.2~0.21 2. My Cd (average) was found to be 1.3 and expected value is 1.1 3. Lift coeff was found to be ~1 while expected value is 0.5 (which is where the real problem is). The value is almost double. I would like to know if there is anything else I should be doing and where am I actually missing. Thank you very much. Ojha 

June 19, 2011, 06:51 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
Are you sure the lift coefficient is based on the correct reference area? Mixing up radius and diameter explains a factor of two.
Now that you have got a model which is pretty close (in everything except lift), if you wish to get more accurate you should do a proper convergence study to guide you into what options you have left. Consider Richardson extrapolation, grid convergence indexes and similar techniques to really squeeze the last drip of accuracy out of it. "Computational Fluid Dynamics" by Roache is the seminal textbook in this area, but a summary of some key concepts is here: http://www.cfdonline.com/Wiki/Ansys..._inaccurate.3F  and the reference to JFE is highly recommended reading. 

June 19, 2011, 12:01 

#3 
New Member
Mayank Ojha
Join Date: May 2011
Posts: 22
Rep Power: 7 
Thanks Glen. Appreciate your feedback.
I had another question. I am using 5% turbulence ( i.e medium intensity) at the inlet. Should this be causing the problem ??? I realize that this problem is slightly tricky because prior to separation, the flow is laminar while after separation its completely turbulent. How about the BSL model and the nonlinear RNG ke model ? Literature say that they have been proved to be a good estimate for near wall flows and flows with rotation (which is exactly my case). But literature also say that SST should be the best (which it does not, so far). Right now am running an LES for this case. Lets see how it goes. BTW what is Richardson extrapolation and grid convergence indexes ??? Looks like the library is out of Roache's txt book. 

June 19, 2011, 18:35 

#4  
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
Quote:
Quote:
Quote:
Quote:


June 20, 2011, 13:01 

#5 
New Member
Mayank Ojha
Join Date: May 2011
Posts: 22
Rep Power: 7 

June 20, 2011, 18:47 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
LES is a whole field of CFD in itself, your library should have textbooks on the topic. Otherwise start with general CFD books (eg Anderson) as they often introduce the concepts of LES.


June 21, 2011, 13:17 

#7 
Member
Paulo Alexandre Costa Rocha
Join Date: Mar 2009
Posts: 71
Rep Power: 9 
AFAIK a Reynolds number of 10000 is not fully turbulent for external flows. Pure RANS models are not proper. Maybe LES or a transition model can give better results.
My two cents. Best regards, Paulo Rocha 

February 3, 2012, 21:49 

#8 
New Member
cfd
Join Date: Oct 2011
Posts: 17
Rep Power: 6 
Hi,
I am modelling the turbulent flow over the cylinder by fluent(kepsilon model) I have a basic question about defining the boundary conditions for the domain. The domain is rectangular. for specifying velocity inlet for domain, I have to define the turbulent intensity as well as hydraulic diameter, I want to know I should use the diameter of the cylinder as the Hydraulic diameter or the cross length of the inlet face?! Thank you 

February 4, 2013, 04:44 

#10 
New Member
Join Date: Jan 2013
Posts: 3
Rep Power: 5 
Hi everyone,
I'm running a 3D cylinder simulation in FLUENT with a Re of 100,000. I've used a Transition SST model with a timestep size of 1e06 and 10,000 time steps. I'm trying to obtain the drag coefficient from this simulation The results I've obtained shows an oscillatory motion, with the Cd value hitting a peak of 1.21 midway before falling back to 0.726 at the end of the 10,000 time steps. Can anyone tell me why has this occurred? What do I need to do for the Cd values to be stablised? Thanks! Regards, Hee 

February 4, 2013, 05:10 

#11 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
Your question is about Fluent, not CFX and this is the CFX forum. But the answer is the same: http://www.cfdonline.com/Wiki/Ansys..._inaccurate.3F


May 18, 2015, 20:11 

#12  
Senior Member
Emre G
Join Date: May 2011
Location: Turkey
Posts: 120
Rep Power: 7 
Quote:


May 18, 2015, 21:03 

#13 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
My post #11 is a pretty clear response to this.


May 19, 2015, 02:09 

#14 
Senior Member
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 312
Rep Power: 7 
Is zonal RANS a thing? I.e. define two subdomains? A laminar one and a turbulent RANS one?


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Modified ke turbulence model UDF  Travis  Fluent UDF and Scheme Programming  6  October 7, 2015 13:54 
question about turbulence model selection and sensitivity  karananand  Main CFD Forum  1  February 26, 2010 05:41 
Centrifugal Pump and Turbulence Model  Michiel  CFX  12  January 25, 2010 04:20 
SSG Reynolds Turbulence Model  Georges  CFX  1  February 28, 2007 17:15 
Turbulence model  Herry  Phoenics  1  May 29, 2003 13:19 