Flow past cylinder (Re=10^5)
Hello !
I am simulating incompressible unsteady flows past a cylinder with StarCCM+. The Reynolds number is around 10^5, assuming that diameter (D), density and viscosity are 1 (to make things simpler for me, as long as nondimensional values are correct, it should not influence on the computation, I am wrong?) and upstream velocity U=10^5 m/s. The turbulence model is realizable kepsilon, with standard values of turbulence intensity and turbulence length scale (various models and parameters have been tried, results are roughly the same). The grid is a conventional 2D grid. The distance from the velocity inlet to the cylinder is 15*D, 25*D from the cylinder to the pressure outlet and 15*D from the cylinder to the sides, supposed to be symmetry boundaries. The mesh is structured with quadcells. The timestep size is 10^6 s, because the Strouhal number (St) will be around 0.2, so the vortex shedding frequency will be around f=St*U/D=0.2*10^5/1=2*10^4 Hz so the period about 1/(2*10^4)=5*10^5 s. I expect about 50 points per period of oscillation of the aerodynamic coefficients. Given this timestep size and the refinement of my mesh, the maximum Courant number remains below 50. The maximum wall Y+ is also maintained below 50. The number of inner iterations is 10. It seems to be enough. With the same geometry, laminar simulations give satisfactory results for lower Reynolds numbers. The problem is the value of the drag coefficient, known to be around 1.2 at Re=10^5 according to many documents. No matter the changes in any parameter I have tried, my mean value of the coefficient is 0.255, which is much too far from the result I expect. So my first questions are what could be the reasons of such a difference? Is it due to the grid? Should I enlarge the domain (I’ve given a try: it increases the coefficient, but not enough)? Which parts of the domain should I extend (distance from the velocity inlet to the cylinder, from the cylinder to the sides, etc)? What is the influence of the turbulence intensity and turbulence length scale that are specified in the inlet boundary? I have tried to decrease the timestep size to 10^7 s too, but the improvement is not that great: I get a mean drag coefficient of about 0.32. The vortex shedding frequency and the amplitude of the lift coefficient oscillation also vary. But the calculations take so much time, I can’t reduce the timestep once more, it would last days, I can’t afford it. My second set of questions is why does the drag coefficient vary with timestep size? How could I know my timestep is small enough? And what could I do to get more accurate results? What else should I try? Any help is welcome. If you need any further detail, let me know. Thank you, Laurent http://www.cfdonline.com/Forums/mem...ure25mesh.png >> Mesh << http://www.cfdonline.com/Forums/mem...oefficient.png >> Effect of a change in the timestep size on drag coefficient << http://www.cfdonline.com/Forums/mem...oefficient.png >> Effect of a change in the timestep size on lift coefficient << 
Hello,
I would actually like to ask you something about mesh generation. I am going to do similar things, but the cylinder would rotate somehow. Do you or anyone have any suggestions on how to generate the mesh (2D or 3D) for this type of problem in StarCCM+ and StarCD (I am using at the moment)? Thank you. B Ken 
http://www.grc.nasa.gov/WWW/wind/val...l/lam_cyl.html
check this site, though it is for a Laminar case Good Luck 
This is just a shot in the dark.
"Theory" gives 1.2 x 10^4. You get 0.255x10^4. The ratio is 4. Now the ratio of diameter^2/radius^2 = 4. (radius & diameter are for your cylinder). Would it be worth checking to see that the theory and your results are on the same basis? Hope this is some help  even if it just removes one potential problem. 
I'm not sure I understand what you mean.
First, the expected result of the drag coefficient is 1.2, not 1.2x10^4... Then, for the calculation of the drag coefficient, I use the diameter as reference length, radius shouldn't be used. Anyway, I don't think the problem comes from the calculation of the coefficient, because laminar simulations give good results. The basis of my results and the one of the expected results are the same as long as Reynolds numbers are the same. In my case, "drag coefficient [is] only dependent on the Reynolds number" (from "Boundary layer theory" by Hermann Schlichting, p. 9). Here is an example of the document where you can check the expected result: http://www.cfdonline.com/Forums/mem...hlichting2.png >> Evolution of Cd vs. Re, from "Boundary layer theory" by Hermann Schlichting, p. 19 (preview available on http://books.google.com/books?id=8YugVtom1y4C) << Good try, thank you. Any other idea ? Laurent 
The mesh needs to be remeshed. What is needed is a ring around the cylinder that has a dense distribution, these four walls have a lot to do with these wild oscillations in your results

Actually I have tried to add a boundary layer of small cells.
Perhaps not small enough... What should be the size of the first row? I thought it was enough refined, as my max y+ is between 1 and 100 (about 3050). If I put smaller cells, max y+ will be lower, how should I handle the wall treatment? Thanks for the tip, Laurent 
Laurent
only I can ask you to review the hypersonic cylinder case on the validation server given before (there is a shot for the mesh used), you need a high density ring around the cylinder, simply said, remove these walls from your domain. the y+ values seem to be OK Good Luck 
flow around cylinder
From my experience, you will not get the correct results using a 2D mesh. Flow around a cylinder is 3D for Re>200.
For Re< 3 x 10⁵, the boundary layer is laminar, so you need to use a lowRe model. Probably the kepsilon model that you are using is a highRe one, so you are getting the drag values for turbulent boundarylayer. The best option for your flow is to use the SST komega model, but without 3D effects you are going to get high Cd values. Your Cd probably will be around 1.7. regards, Fabio 
Hi,
Thank you very much, Fabio, for your advices. This study is in fact a test case, a kind of preliminary work to the study of a much more complex geometry (a 3D vertical axis wind turbine). The Reynolds number will be about 5x10^5, but I know I can't get results for the flow over a cylinder at this critical point, so I first study Re=10^5. In reality, I guess flow will be fully turbulent over a wind turbine (To me, lowRe models doesn't seem to be appropriate, is it ?), perhaps I should run tests at Re=10^6 or 10^7 ? My goal is to train and get necessary knowledge. But I feel a little worried, because it is so difficult to get accurate results on what I thought was a simple test case that I just can't guarantee usable results on a complex geometry. If I don't get good results for the flow past a cylinder, how could I get good results for another geometry ? I'm still gonna work on the cylinder following your help, but does anyone have any advice on how I should work to study a wind turbine ? Am I going the right way (I mean: study the flow over a cylinder (fixed and rotating), be sure to get accurate coefficients, and then study the flow over a wind turbine) ? By the way, can anyone give me reasons why drag coefficient is so influenced by the timestep size ? If you have an answer for any of those questions, any help is welcome. Thanks once again, Regards, Laurent 
All times are GMT 4. The time now is 11:05. 