
[Sponsors] 
September 19, 2014, 03:46 
drag coefficient calculation

#1 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
Hello everyone,
I need some assistance on understanding why in the drag coefficient i am getting the correct frequency according to the Re number but the amplitude is different from studies in papers. icoFoam is the solver and for different meshes i get different amplitudes but same frequency. the amplitude should oscillate around 1.44 and it oscillates around 3. Does anybody know about amplitude dependency on either numerical or physical factors that could help me understand why. Last edited by Konstantine; September 19, 2014 at 04:02. Reason: more info 

September 19, 2014, 05:54 

#2  
New Member
mohsen cheraghi
Join Date: Jun 2010
Location: Switzerland
Posts: 28
Rep Power: 16 
Quote:


September 19, 2014, 10:56 

#3 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
thank you my friend mohsen for your response.
So if i understand correctly you suggest that the calculation adds up to more than it should do, meaning the magnitude as it is calculated by the total drag force + the negative and positive oscillations so it adds up to more than just the magnitude. You could have a point there but how is it calculated? i am searching now from the forceCoeffs.C routine to understand better what it actually calculates and what formula it uses for the coefficient. 

September 19, 2014, 14:26 

#4 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
Dear mohsen,
i found how it calculates the coefficient. it is the usual formula of force/(dynamic pressure* area). then i found that it adds pressure, viscous and porous forces to calculate the total force. so it is something about the last three that is wrong, i am thinking that maybe porous is not zero so it gives different result than usually or pressure or viscous friction is overpredicted. i am not sure yet i will continue searching. Last edited by Konstantine; September 19, 2014 at 14:45. Reason: more info 

September 22, 2014, 12:45 
Higher Average CD

#5 
New Member
AI
Join Date: Jun 2014
Posts: 17
Rep Power: 12 
Hello,
I'm having what I assume is a similar issue, I'm getting higher average Cd than what I should be getting, however, not as much as yours. I started looking at sectional Cd in the vertical direction and i'm finding two distinct groups of sectional Cd curves, one group is higher than the other. I described the issue in the following thread http://www.cfdonline.com/Forums/ope...direction.html Please take a look and let me know what do you think. Thank you. 

September 23, 2014, 07:57 

#6 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
hello my friend sherif35
thank you for your concern on the topic. total drag= pressure drag + friction drag so what you see there may be the pressure drag and the friction drag , those will oscillate around a nonzero value in contrast to the lift which oscillates around zero. 

September 23, 2014, 09:29 

#7  
New Member
AI
Join Date: Jun 2014
Posts: 17
Rep Power: 12 
Quote:
The sectional Cd curves I show are for the total drag on the section. I'm still trying to figure out why certain sections will oscillate about certain value while the rest oscillate around another!! 

September 23, 2014, 10:03 

#8 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
i guess you talk about sections of the cylinder, sorry i did not understand the sectional forces well, it is reasonable that because of different pressures (due to vortex shedding) around the cylinder you will have different forces on the cylinder wall points, therefore it would oscillate around different values for each section but the frequency should be the same as vortex shedding is the phenomenon responsible for these pressures. if you have the same forces let's say on all sections then either you would have symmetric vortex shedding which can be the case ( i think only numerically achieved) some times but is not stable and will brake after a couple of shedding periods or symmetric laminar flow like on Re=3050 where there is a wake bubble stable that just circulates in two symmetric regions the water behind the cylinder and does not shed any vortex.
In general i believe you are interested in the total drag and not on sectional unless you want to compare with experimental on particular points of the cylinder. 

September 23, 2014, 10:31 

#9  
New Member
AI
Join Date: Jun 2014
Posts: 17
Rep Power: 12 
Quote:
You are correct, I'm interested in the total drag, but what I don't really understand is why there are two distinct group of sectional Cd curves. If they are all different I would have assumed that this is 3D effects. But in my case, it looks like certain sections produce certain results. If you look at the attached picture for Re=10,000 you will see section 3&6 Cd curves are exactly the same (until the shedding starts) and the rest of the sections are higher. I think if I understand what causes this issue for the sectional Cd, it will help me with the average drag for the whole cylinder. So I'm wondering if you have seen similar sectional Cd issue with your simulation? 

September 23, 2014, 13:55 

#10 
New Member
Join Date: Sep 2014
Posts: 10
Rep Power: 11 
everything is about pressure, i have seen different pressures around the cylinder wall yes ( therefore the forces and the coefficients are different), this is reasonable. on Re=200 you may have 3d effects and maybe transition to turbulence which makes the phenomena more complicated. depending on how many parts you have devided the cylinder there may be parts which have the same pressure curve e.g. the front parts of the cylinder before the seperation. after the separation, the pressure will oscillate due to vortex shedding or be more chaotic due to turbulent eddies which may appear because of transition to turbulence. you calculate the total force and then the coefficient so the average coefficient will be in the middle around which it oscillates. now if you want to relate the sectional coefficients on particular parts of the cylinder with the total average, think that calculate the total force you add the forces on each node where it is calculated so eventually the sectional average Cds would add up to the total average Cd. think if you really need that have you found any similar research, it has been many times when i was stuck with ideas and i was searching without a reason. you may be correct ok i am not questioning your idea but is it worth it as a result what do you want to prove.


Tags 
coefficient, cylinder flow, drag, forces and force coeff. 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Drag Coefficient for Ellipse Form (2D)  Ketut Utama  Main CFD Forum  8  December 11, 2014 11:03 
Question about heat transfer coefficient setting for CFX  Anna Tian  CFX  1  June 16, 2013 06:28 
problem with saving drag coefficient  colopolo  FLUENT  5  April 12, 2013 10:59 
Post Processing Drag Coefficient  squanto773  OpenFOAM PostProcessing  1  March 7, 2012 09:43 
Plot drag coefficient against Y(distance)  edwin  FLUENT  0  February 4, 2008 09:41 