|
[Sponsors] |
August 3, 2012, 08:18 |
Problem with drag coefficient
|
#1 |
New Member
Mindaugas
Join Date: Jun 2012
Posts: 25
Rep Power: 14 |
Hi everyone, I'm modeling Turbulent flow around cylinder, the main aim is to dicover Cd dependence on Reynolds numeber.
The cylinder is in the infinite field, material air, Diameter of cylinder 1m, Re=10000. Model RNG,k-epsilon. I'm changing only tubulent kinetic energy and tubulent Dissipation Rate. And the resut is Cd=0.58, in case it should be approximatelly 1.1-1.2. Maybe someone had done similar work, and have some advices. Regards, Mindaugas |
|
August 3, 2012, 08:31 |
|
#2 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49 |
I encountered the same problem at Re=10^5. The best result I got was at about 70% of the experimental value.
Maybe you could try a k-omega based turbulence model or even RSM. Don't worry though, I wouldn't expect a RANS model to yield perfect results in the detached flow over a blunt body. |
|
August 3, 2012, 08:36 |
|
#3 |
New Member
Mindaugas
Join Date: Jun 2012
Posts: 25
Rep Power: 14 |
I thing the main problems are:
1) Bad mesh, at monday I'll try to change it. 2) Bad calculation of Turbulent kinetic energy, turbulent dissipation rate. Also I have tried to change turbulent intensity and hydraulic diameter, but as I understand hydraulic diameter is essiantial only for pipes. But project is 2D cylinder. Laminar flow was quit good 1-3% inaccurate. I thing for this exeriment k-epsilon model is good. |
|
August 3, 2012, 08:48 |
|
#4 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49 |
I already did a mesh dependency study for this case.
Changing the turbulent kinetic energy and its dissipation rate at the inlet is not the right way to fix the results. You can assume that the measurement was carried out with very low turbulent intensities upstream of the cylinder. So increasing k at the inlet in your simulation might yield a higher C_d, but then the simulation no longer corresponds to the experiment. |
|
August 3, 2012, 08:54 |
|
#5 |
New Member
Mindaugas
Join Date: Jun 2012
Posts: 25
Rep Power: 14 |
But you certainly must change parametres on the boundary conditions, because turbulence parametres dependent on Reynols number.
quate Note by default, the FLUENT GUI enters k=1 m²/s² and ε =1m²/s³. These values MUST be changed, they are unlikely to be correct for your simulation. |
|
August 3, 2012, 09:09 |
|
#6 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49 |
But which value do you want to use? Unless you have a model of the wind tunnel that was used for the experiment, it seems impossible to me to "guess" correct values.
Aren't wind tunnels usually equipped to produce a homogneous flow field with very little turbulent disturbances? |
|
August 3, 2012, 18:38 |
|
#7 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49 |
BTW: The turbulence parameters at the inlet do NOT depend on the Reynolds number of the flow around the cylinder.
Just imagine that in the experiment, instead of the cylinder with 1m diameter, a smaller with d=0.1m is used. The Reynlods number decreases by a factor of 10, but the flow in the wind tunnel is still the same with the same turbulent quantities. |
|
August 6, 2012, 02:45 |
|
#8 |
New Member
Mindaugas
Join Date: Jun 2012
Posts: 25
Rep Power: 14 |
All parametres can be defined:
http://hpce.iitm.ac.in/website/Manua...ug/node217.htm I don't know what is a problem? Maybe others have modeled turbulence flow and results were approximately equals to experimental ones? |
|
August 7, 2012, 03:50 |
|
#9 |
New Member
Mindaugas
Join Date: Jun 2012
Posts: 25
Rep Power: 14 |
Results become more accurate. Re=10000 , Cd=0.9. I have changed turbulent model Spalart-Allmaras, turbulent intensity 0.5 and turbulent length scale 0.07, but still the result should be 1.15-1.2.
Maybe some advices, how to improve results. |
|
August 7, 2012, 04:34 |
|
#10 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49 |
A wall-resolving LES should yield better results here.
Still there is the problem of the unknown Reynolds stresses and length scales at the inlet, which becomes even more crucial in a LES. In my opinion, 0.9 is quite a good result as an approximation of the drag coefficient. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Drag Coefficient Convergence Problem | John | FLUENT | 18 | June 24, 2023 10:22 |
problem with saving drag coefficient | colopolo | FLUENT | 5 | April 12, 2013 11:59 |
lift and drag ceofficient problem | icemaniac178 | CFX | 6 | August 17, 2011 19:40 |
Problem in plotting drag coefficient | Sarinagara | FLUENT | 2 | September 28, 2010 07:54 |
drag coefficient | pukiwawa | FLUENT | 10 | August 29, 2008 16:31 |