CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Drag Coefficient calculation for flow over a 2D Cylinder at High Reynodls Numbers

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

Like Tree3Likes
  • 2 Post By ghorrocks
  • 1 Post By evcelica

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 19, 2019, 09:51
Post Drag Coefficient calculation for flow over a 2D Cylinder at High Reynodls Numbers
  #1
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Good day everyone,

At the moment I'm carrying out simulations in CFX of flow over a 2D Cylinder in a rectangular domain immersed in a high reynolds flow to determine its drag coefficient.

The goal of the simulations is to obtain a Drag Coefficient vs Reynodls curve that resembles as closely as possible the graph attached in the first image, for Cd vs Re for flow over Cylinders.

I'm using an structured mesh created in ICEM as a 2D mesh and the extruded one element in the Z direction so it can be used in ICEM. This mesh is then used in CFX to run the simulations for flows in the range of 10^4 to 10^7 Reynodls number.

When I run the simulations between 10^4 and 10^6, the results I obtain for the drag force over the cylinder, using the following expression:

force_x()@CYLINDER

provide with acceptable values of the drag coefficient, almost equal in some cases to the graph I'm trying to reproduce. The formula I use to obtain the Drag Coefficient is as follows:

Cd = force_x@CYLINDER /(0.5*DENSITY*V^2*A)

The cylinder has 1 meter in diameter, and the fluid is air at 25ºC. The velocity is dependent on the Reynolds number at which I'm running the simulation.

The problem begins when I run the simulation for Reynodls numbers between 10^6 and 10^7. In this case, using the same expression to determine the drag force in the cylinder, and then the same expression to determine the drag coefficient, the values I get are around half and a third of what the actual values should be.

Mesh:
Structured mesh, with a high resolution around the cylinder wall. The first cell element thickness is determined for a y+ value of 1.
The domain size is 10 times the diameter length in the upstream direction, and 20 times in the downstream direction, to allow the flow to fully develop and avoid any backflows or undesired effects. The width is 7 diameters in each direction as well.

Turbulence model:
SST, with no transitional turbulence, as it can be seen in the second attached image.

Type of Simulation:
Steady State

Time Steps:
Local Time Factor: 30

Convergence Criteria:
1*10^-8
(The convergence criteria is reached)

All of this can be seen in the third picture attached.

It should be mentioned as well that this problem does not appear when meassuring the Lift Coefficient, as it remains close to 0 in all simulations, independently from the Reynodls number.

Can anybody provide me with any insight into why for the initial Reynodls values, until the Drag Coefficient drop in the graph, the Cd values are more than sattisfactory, but as the Reynodls increases in the range of 10^6 to 10^7, the results are completely off? Is there obvious mistake?

Has anybody attempted to do a similar simulation with succesfull results?

It's a very simple simulation and I cannot see where there migth be some errors in the setup.

Thank very much in advance for the replies.

Daniel Barreiro
Attached Images
File Type: png Cd vs Re for Cylinders.png (161.5 KB, 65 views)
File Type: jpg Mesh.JPG (170.8 KB, 40 views)
File Type: jpg Turbulence Model.JPG (44.9 KB, 36 views)
File Type: jpg Steady State Time Set Up.JPG (53.8 KB, 29 views)
File Type: jpg Drag Force Re 6x10^6.jpg (50.3 KB, 43 views)
DanielBarreiro is offline   Reply With Quote

Old   February 19, 2019, 10:56
Default
  #2
Senior Member
 
M
Join Date: Dec 2017
Posts: 640
Rep Power: 12
AtoHM is on a distinguished road
Resolution of the mesh seems to be high indeed, but still: you said the layer thickness is chosen so yplus is 1. What Reynolds number did you use to determine the neccessary cell height?

If you designed it for Re 10^4, you will get a higher yplus when simulating higher Re-number with the same mesh. This will in turn influence the boundary layer resolution and might lead to changes in the drag.


Check yplus values in post and see if its getting out of range for the higher Re-number simulations.
AtoHM is offline   Reply With Quote

Old   February 19, 2019, 11:30
Default
  #3
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Quote:
Originally Posted by AtoHM View Post
Resolution of the mesh seems to be high indeed, but still: you said the layer thickness is chosen so yplus is 1. What Reynolds number did you use to determine the neccessary cell height?

If you designed it for Re 10^4, you will get a higher yplus when simulating higher Re-number with the same mesh. This will in turn influence the boundary layer resolution and might lead to changes in the drag.


Check yplus values in post and see if its getting out of range for the higher Re-number simulations.
The y+ value of 1 has been calculated for a Reynolds number of 10^7, so when the simulations are run with a lower Reynodls number, the y+ value I get in CFX Post is always below 1.

Still, the results keep being innacurrate for Reynolds numbers above from 10^6.
DanielBarreiro is offline   Reply With Quote

Old   February 19, 2019, 11:57
Default
  #4
Senior Member
 
M
Join Date: Dec 2017
Posts: 640
Rep Power: 12
AtoHM is on a distinguished road
Alright, this was just the first thing coming to mind as it was missing from your besides that very clear explanation of your setup. Maybe someone else can point to a solution.
AtoHM is offline   Reply With Quote

Old   February 19, 2019, 12:11
Default
  #5
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Thank you very much for the reply anyways. Let's see is someone else is able to point out the problem.
DanielBarreiro is offline   Reply With Quote

Old   February 19, 2019, 17:48
Default
  #6
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,166
Rep Power: 23
evcelica is on a distinguished road
What are you comparing to? Red line and blue dots are what in your first graph?
Cylinders do have this dip around the Reynolds numbers like your red line shows, I think you are already aware of that though.

First off, I have found that drag force on a 2D cylinder is actually quite a difficult simulation @ any Reynolds number due to the difficulty in CFD predicting the flow separation point around the back of the cylinder. I always got a slightly lower drag coefficient than literature, and I don't know a way to correct for this. I feel CFD could better predict a more complex geometry than this simple cylinder, which turns out to be much more dificult than one would expect.

Second, are you using incompressible air? (You stated "air @ 25C" which is incompressible). If so compressibility effects could be playing a role.
Air ideal gas may be more appropriate for starters. One step further (probably not required, but a sensitivity analysis would determine this) may be using real gas with the total energy model instead of an isothermal model.
evcelica is offline   Reply With Quote

Old   February 19, 2019, 17:55
Default
  #7
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,665
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You are using local time scale factor. Quoting from the FAQ: https://www.cfd-online.com/Wiki/Ansy...gence_criteria

Use Local Timescale Factor. A factor of about 5 is a good guess to start with. If this is successful you should run the final few iterations to convergence with a physical timescale (not local timescale all the way to convergence).

So you need to do the final run to convergence with physical time scale, not local time scale factor.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   February 21, 2019, 10:22
Default
  #8
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
You are using local time scale factor. Quoting from the FAQ: https://www.cfd-online.com/Wiki/Ansy...gence_criteria

Use Local Timescale Factor. A factor of about 5 is a good guess to start with. If this is successful you should run the final few iterations to convergence with a physical timescale (not local timescale all the way to convergence).

So you need to do the final run to convergence with physical time scale, not local time scale factor.
I have run a few simulations with a Local Time Step of 5, as well as 10, since it's a very common value too, but I do not get consistent results, and the Cd value for Reynodls between 10^6 and 10^7 are still very far from what the actual values should be.

Any idea of what else migth be causing this problem?
DanielBarreiro is offline   Reply With Quote

Old   February 21, 2019, 10:26
Default
  #9
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Quote:
Originally Posted by evcelica View Post
What are you comparing to? Red line and blue dots are what in your first graph?
Cylinders do have this dip around the Reynolds numbers like your red line shows, I think you are already aware of that though.

First off, I have found that drag force on a 2D cylinder is actually quite a difficult simulation @ any Reynolds number due to the difficulty in CFD predicting the flow separation point around the back of the cylinder. I always got a slightly lower drag coefficient than literature, and I don't know a way to correct for this. I feel CFD could better predict a more complex geometry than this simple cylinder, which turns out to be much more dificult than one would expect.

Second, are you using incompressible air? (You stated "air @ 25C" which is incompressible). If so compressibility effects could be playing a role.
Air ideal gas may be more appropriate for starters. One step further (probably not required, but a sensitivity analysis would determine this) may be using real gas with the total energy model instead of an isothermal model.
Soryy I did not make that clear enough. I'm trying to reproduce the red curve in the graph as closely as possible to validate the mesh and the turbulence model I'm using for Reynodls numbers between 10^6 and 10^7.

I'm aware of the fact that the Cd it's always a bit less compared to the actual experimental value, but in my case, the Cd values I obtain with the simulations are about half, or a third of what they should be, which is not an acceptable difference.

As you well say, the energy mode I do not think is neccesary, as it a very simple simulation that shouldn't be affected by it.

I have gone once again over my set up parameters but I do not seem to find anything unusual. Any idea of what migth be causing this problem?
DanielBarreiro is offline   Reply With Quote

Old   February 21, 2019, 15:12
Default
  #10
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,166
Rep Power: 23
evcelica is on a distinguished road
So the blue dots are your results?
Perhaps the transient nature of the vortex shedding is not being resolved in your steady state model. That dip in the Cd lines up pretty well with a distinct vertex shedding mode, and that dip can also change with smooth vs rough walls.

Is your entrance and exit far enough away from the cylinder?

Are you are using ideal gas or constant properties?
evcelica is offline   Reply With Quote

Old   February 21, 2019, 19:29
Default
  #11
New Member
 
Daniel Barreiro Clemente
Join Date: Feb 2019
Location: Munich
Posts: 15
Rep Power: 7
DanielBarreiro is on a distinguished road
Quote:
Originally Posted by evcelica View Post
So the blue dots are your results?
Perhaps the transient nature of the vortex shedding is not being resolved in your steady state model. That dip in the Cd lines up pretty well with a distinct vertex shedding mode, and that dip can also change with smooth vs rough walls.

Is your entrance and exit far enough away from the cylinder?

Are you are using ideal gas or constant properties?
The blue dots are my results, and the tendency is a straigth line instead of increasing once again after the dip at around 8*10^5 Reynodls number.

The domain is big enough to allow for the flow to be fully developed. Around 30 times the cylinder diameter in the downstream direction and around 10 times in the upstream direction, as well as 10 times in each width direction.

I'm using constant density for the air at 25ºC.
DanielBarreiro is offline   Reply With Quote

Old   February 25, 2019, 11:09
Default
  #12
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,166
Rep Power: 23
evcelica is on a distinguished road
Well I would try ideal gas first.
Then run a transient model using time steps small enough to resolve the vortex shedding.
evcelica is offline   Reply With Quote

Old   February 25, 2019, 20:09
Default
  #13
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,665
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The dips in the CD curve are mainly due to boundary layer effects. The big drop at around 2-3E5 is due to the boundary layer becoming turbulent and staying attached longer. At Re less than this you are transitioning from Stokes flow to laminar boundary layer flow. At Re higher that 2-3E5 the separation point moves forward again so you get less pressure recovery and the CD goes up.

So if you want to model the bumps and wiggles in the CD curve the most important thing is to have a very accurate boundary layer model. If you are around the 2-3E5 point you will need a turbulence transition model to capture the turbulent boundary layer with everything else laminar. For Re above above this you will need a very good boundary layer mesh and turbulence model. For Re below this you will need to match your boundary layer mesh to the laminar boundary layer you get at that Re.

Vortex shedding is a smaller contributor to the CD curve variation. The big effect is boundary layer effects as I describe.
belier1988 and DanielBarreiro like this.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   February 26, 2019, 10:40
Default
  #14
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,166
Rep Power: 23
evcelica is on a distinguished road
Thanks for the info Glenn,

I am interested in this as well, and appreciate the knowledge you share.
DanielBarreiro likes this.
evcelica is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
2D cylinder drag coefficient miku11 OpenFOAM Running, Solving & CFD 0 June 28, 2016 08:36
doubt about drag coefficient (cylinder) CataV OpenFOAM Running, Solving & CFD 7 March 10, 2016 11:54
drag and lift calculation of a cylinder in wind tunnel fewgoodmen Main CFD Forum 4 February 24, 2014 16:05
Incorrect Drag and Drag Coefficient for flow over a cylinder ozzythewise Main CFD Forum 8 June 13, 2012 07:24
Circular cylinder drag at high Re Anton Lyaskin Main CFD Forum 0 January 14, 2003 05:48


All times are GMT -4. The time now is 22:23.