
[Sponsors] 
komega SST simulation of turbulent flow around a circular cylinder 

LinkBack  Thread Tools  Display Modes 
September 23, 2010, 05:25 
komega SST simulation of turbulent flow around a circular cylinder

#1 
New Member
Dan M
Join Date: Sep 2010
Location: Munich, Germany
Posts: 20
Rep Power: 7 
Hi everyone,
I have been analyzing the performance of the kepsilon and komega SST turbulence models using a test case of a 2D circular cylinder. The numerical model was set up for Re = 3900 in order to compare the results with experimental data. As comparison parameter the drag coefficient Cd was calculated which according to experiments should be at around 0.99 for Re = 3900. The results significantly deviate from this value, especially for the komega SST model. With OpenFoam, an average value of about 1.6 was calculated using komega SST. The following boundary conditions were used for the model: At the wall k = 0 (and kqRWallFunction) omega = 10*6*nu/(0.075*y^2) (and omegaWallFunction) nut = 0 (and nutWallFunction) At the inlet k = 3/2*(U*I)^2 omega = Cmu^(0.25)*k^(0.5)/L nut = calculated uniform 0 internalField at timestep 0 k = 3/2*(U*I)^2 omega = Cmu^(0.25)*k^(0.5)/L nut = 0 The calculation resulted in an average yPlus value of around 2 with the maximum just over 4. Both mesh refinement and a test case without the use of wall functions did not help to significantly improve the results. Further, a second test case, which made use of the kepsilon model, seemed to improve the results significantly using the exact same boundary conditions as before which leads to the following values for epsilon: At the wall epsilon = 0 (and epsilonWallFunction) At the inlet and internalField at timestep 0 epsilon=Cmu^(0.75)*k^(1.5)/L Here, Cd averaged about 0.9. I find this rather surprising, as the literature suggests that the komega SST model should perform better for this particular case (or at least just as good). Has anyone experienced similar difficulties or can tell me if there is a mistake in the boundary conditions? I would like to use komega SST for future calculations and need to make sure it is set up correctly. Any suggestions for a possible source of error is much appreciated. Thanks a lot, Dan 

September 27, 2010, 03:19 

#2 
New Member
Dan M
Join Date: Sep 2010
Location: Munich, Germany
Posts: 20
Rep Power: 7 
I haven’t received an answer to my post yet and am still facing the same problems with my model. I ran some further tests and by now am fairly certain that my solution is neither mesh dependent nor does the use of wall function significantly change the results. Maybe I should post my question a bit more general. I am interested what boundary conditions are typically used at the inlet and walls when making use of the komega SST model. Is there an alternative approach to the boundary conditions I have listed in my fist post? I would very much appreciate any input on this.
Thanks a lot, Dan 

October 17, 2010, 17:18 
reply

#3 
Member
joegi
Join Date: Nov 2009
Location: genoa
Posts: 67
Rep Power: 8 
Hi,
I had a similar problem but with a different geometry. For the yplus value you reported you should not use wall functions, if you want to use WFs you should have 30<yplus<300 (if i well recall). If you still want to use that mesh, you should use a near wall treatment (do not use wall functions in your BCs) and an appropriate lowre turbulence model. In my personal experience, I found that is extremely tricky to get good results. Hope this helps, joegi 

October 18, 2010, 09:58 

#4 
New Member
Michel
Join Date: Jun 2010
Posts: 7
Rep Power: 7 
Hi,
i have the same problem with the komegasst model. I have a Renumber of 5000 and komegast isnt working, with and without bcs. I tested the kepsilon model with the same case and it is working, but the problem is, that kepsilon can not compute seperation that good. Perhaps there is a bug in the komegasstmodel? 

October 18, 2010, 10:25 

#5 
New Member
Dan M
Join Date: Sep 2010
Location: Munich, Germany
Posts: 20
Rep Power: 7 
Thank you very much for the answers. I have spent some more time with this particular problem and suspect that one of the problems is that the Reynolds number is in the transitional regime between fully laminar and fully turbulent boundary flow. This means that while the flow is turbulent in the regions behind the cylinder, the boundary around the cylinder itself is still laminar (the wall has a damping effect on the development of turbulence). Only for high Reynolds flow the near wall regions start becoming turbulent as well. For komega SST calculations with Reynolds numbers in the transitional regime this means that a turbulence model is used near the wall to predict a behavior that in practice is actually laminar, which might be one source of error. I have started to use the komega SST model for high Reynolds flows, but have found that it is very difficult to get satisfactory results here as well.


October 18, 2010, 10:43 

#6 
Member
joegi
Join Date: Nov 2009
Location: genoa
Posts: 67
Rep Power: 8 
Hi Dan,
There are some discussions about transitional turbulence models, but so far I think they are not implemented in openfoam. Let's see what the gurus have to say. In any case I think kwsst will give you way much better results than kepsilon (if you manage to get it working for your case), just pay attention to your y+. Also try to change to a lowre model (near wall treatment), LaunderSharmaKE LamBremhorstKE LienCubicKELowRe LienLeschzinerLowRe qZeta SpalartAllmaras I like SpalartAllmaras, but is often used for high speed external aerodynamics. Another thing, I would rather prefer to use a small value for epsilon (108) instead of 0, there is a discussion of this somewhere in the forum. Let us know if you get your simulations working. Have a nice day, Joel 

October 18, 2010, 11:03 

#7 
Senior Member
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 8 
There are blended wall functions in OpenFOAM. omegaWallFunction (as recommended by Menter) and the spalding wall functions for nut. Both should give good results for any y+ value. However, I also believe that the kw SST & the kw can be run as lowre turbulence models without any additional damping functions for the near wall region (these are normally required for the ke models). Wilcox discusses this in his book, Turbulence Modelling for CFD. However, the y+ isn't the only important element. You should make sure you have enough layers off the wall to properly resolve your boundary layer, otherwise you will be wasting your time. I wouldn't recommend anything less than 10 (for a cylinder, you could easily get away with 2030).


October 18, 2010, 11:14 

#8 
New Member
Dan M
Join Date: Sep 2010
Location: Munich, Germany
Posts: 20
Rep Power: 7 
I find interesting what was written about setting epsilon equal to 108. In my experience I have found that I seem to get better results when setting k equal to zeroGradient at the wall. Of course the turbulent kinetic energy should theoretically be equal to zero right at the wall (just like epsilon), but I suspect that the boundary conditions for k, epsilon and omega are in fact applied to the first cell center in OpenFOAM, rather than to the wall itself (unlike for the velocity for example). According to Pope, k is equal to zero right at the wall, but constant (and not equal to zero) in the loglaw region adjacent to the wall. If in fact they are applied to the first cell center in OpenFOAM it would mean that the boundary condition k=0 for the walls, as I had previously suggested, would therefore not be correct. Does anybody know how the boundary conditions for k, omega and epsilon are treated in OpenFOAM? By the way: Mesh refinement (and therefore reduction of yPlus) did not significantly improve my results.
Last edited by DanM; October 18, 2010 at 12:13. 

October 19, 2010, 10:01 
omegaWallFunction blended wall function

#9 
Member
joegi
Join Date: Nov 2009
Location: genoa
Posts: 67
Rep Power: 8 
Hi,
Regarding to the omega blended wall function (omegaWallFunction) for lowre, how should I setup my case?. I guess I have to impose the omegaWallFunction in the viscous walls. But what about k and nut?. Should I use kqrwallfuntion and nutwallfunction on the viscous walls as well? or should I assign to k and nut a fixedValue (1e10) or a zeroGradient?. Finally, I wonder if this model will suffer of some kind of stability problems when using simplefoam?. From my previous experience, I had terrible convergence problems using ke with simplefoam, so I ran all my simulations using pisofoam. Regards jg 

October 19, 2010, 10:56 

#10 
Senior Member
Olivier
Join Date: Jun 2009
Location: France, grenoble
Posts: 252
Rep Power: 9 

July 6, 2012, 22:30 

#11  
Senior Member
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 12 
Quote:
Interesting! My understanding is this: fluctuation reaches to their peak around y+~15~25, so once you put the first cell center in the viscous layer, the k value is small if its y+<5, so setting wall condition of k as fixedvalue to a very small value still make sense.
__________________
~ Daniel WEI  Boeing Research & Technology  China Beijing, China 

December 8, 2012, 06:10 
flow over a circular cylinder.. @ Re=1000

#12 
New Member
Rishitosh Ranjan
Join Date: Dec 2012
Location: Trichirapalli, Tamilnadu, India
Posts: 18
Rep Power: 5 
hiii..
i m very new to in CFD area... i have worked on flow over a single cylinder, two clinder in tandem and tow cylinder in side by side @ diff Re with diff ... i got good results.... but when i started working for Re=1000 .. for single cylinder... i didnt get gud result.... i have use used mesh wid y~min <1 to Y+max<5.. wid Kwsst model... Cd wich i m getting is 1.4 and as per journals and different thesis, i went through, Cd should be 1~ or < 1. i hv tried wid 0.5% turbulence @ inlet and hydraulic dia equal to dia of cylinder.... and @ outlet 0, 0.5 and 0.7 back flow turbulence and same hyd dia... but still... 

December 8, 2012, 09:53 

#13 
Super Moderator

y+ is not only the important factor. Placing first cell in viscous sublayer will give you the Y+ around 1, but it does not mean you have properly resolved the boundary layer. You should ensure 1540 layers in boundary layer with atleast 35 cells in viscous sublayer.
Use second order scheme. If flow is transitional then use the transition model and use the appropriate mesh resolution. In this case your stream wise meshing is also important. Resolving the normal and streamwise in the immediate vicinity will give you the correct values of Cd. Hope it help you. 

December 8, 2012, 13:16 
fow past circular cylinder...

#14 
New Member
Rishitosh Ranjan
Join Date: Dec 2012
Location: Trichirapalli, Tamilnadu, India
Posts: 18
Rep Power: 5 
hello Far sir,
thnx 4 showing interest.. sir i m using 2nd order scheme.. and i hv also tried wid 45 layers of boundary... wid 1st row height 0.007D and 1.1 increment.. but still i dint get result.. dis time y+min was <.4.... Rishi 

December 11, 2012, 10:38 
flow past a circular cylinder.. @Re=1000

#16 
New Member
Rishitosh Ranjan
Join Date: Dec 2012
Location: Trichirapalli, Tamilnadu, India
Posts: 18
Rep Power: 5 
Hello.. Far sir..
Here is detail of d mesh which i m using... @Re1000 * No of nodal grid point on cylinder=100 * Domain size: 16DX35D, D=1cm * No of nodes=12298 * No of cells=19800 * with boundary layer: * 1st row: .004D * Growth factor: 1.05 * No of rows: 45 *Y+min=.104 *y+max=.6115 and while solving (wid 2nd order) BC which i m giving: *velocity inlet: u=1m/sec with turbulent intensity (I)= 0.6% and turbulent length scale(l)= 0.07D *pressure outlet: gauge pr=0 with back flow turbulent intensity = 0.7% and turbulent length scale(l)= 0.007 *upper and lower boundary as symmetry. Rishi 

December 12, 2012, 14:39 

#17 
Member
Aathavan
Join Date: Nov 2012
Posts: 70
Rep Power: 5 
Hi Far and Others,
I am using kepsilon turbulence model, I ran the case many time with different mesh. unfortunately I dint touch system and solution folder, but non of my case shows good results. Now I have changed the turbulence model in to kwsst, under system: Grad  fourth; div  p  Gauss linear, k  upwind omega  upwind and rest of them are Gauss linear, Under Solution: I am using GAMG Solver and GaussSeidel  smoother for p, remaining all smoothSolver. Mesh details: first layer is 3 mm growth factor is 1.05 and 15 layers. the total mesh size is 3.1 millions I have noticed moment, there is no change in that. actual torque should be 70 Nm but the simulation shows 57Nm. as well the pressure raise is much less tan the experimental result. my geometry is an axial flow fan passage my questions are: is there any way to check the origin of my geometry ( I am suspecting this), the torque from simulation is it right? the out show : force ( x y z) viscous (x y z) moment (x y z) viscous (x y z) here is the moment is torque of the fan? please help me to get clear Thanks a lot, Aadhavan 

Tags 
boundary conditions, cylinder, epsilon, omega, sst 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Flow past 2 smooth circular cylinder  slip  FLUENT  0  July 8, 2010 18:45 
Flow over a cylinder  Anna  Main CFD Forum  9  March 24, 2006 15:32 
vortex flow over circular cylinder  Sheila  FLUENT  7  March 21, 2006 15:28 
flow around a circular cylinder  Hanson He  Main CFD Forum  0  November 17, 1998 11:39 
Simulation of the Flow past a circular cylinder using STARCD  M. S. GUEROUACHE  Main CFD Forum  0  October 1, 1998 10:51 