
[Sponsors] 
September 12, 2019, 18:48 
Under Predicting Drag on Cyclist (Fluent)

#1 
New Member
Jacob Flanagan
Join Date: Sep 2019
Posts: 5
Rep Power: 2 
Hi,
I am having some trouble with a CFD simulation I am running on fluent. I am trying to simulate the drag force on a cyclist travelling at 11.5m/s. I would expect a drag force of around 1822N (definitely no less that 17N) however after running multiple simulations with different turbulence models I am only getting drag forces around 1113N, which is far to low. Of the drag I am getting 12N of force from pressure drag and 1N from viscous drag, which seems a little low. I have examined the flow field and everything seems to be looking right in terms of flow structure. I was wondering if anyone would be able to shed some light on this issue as I am not really sure what to change in my model. My setup parameters are as follows: Domain size: 40m x 15.5m x 10m Mesh sizing: Polyhexcore mesh 10 prism layers with 1.2 growth rate and 0.05mm first layer Global min=0.5mm max=1000mm far boi min=0.5mm max=20mm max wheel boi min=0.5mm max=10mm surface mesh of rider min=0.5mm max=10mm max skewness of volume mesh 0.91 number of cells 7.8million Boundary conditions: velocity inlet 11.5m/s pressure outlet 0Pa moving ground boundary 11.5m/s (to match inlet velocity) Symmetry walls No slip wall condition on rider and Frame Rotating boundary condition on wheels Turbulence model: I am using the SSTKw model primarily to solve and using second order SIMPLE scheme. If anyone has any idea on why I am getting such low values for drag I would be extremely grateful. Thanks, Jacob NearWallYPlus.PNG VelocityMagClose.PNG VelocityMagFar.PNG YPlus Far.PNG Streamline 4.PNG 

September 12, 2019, 19:05 

#2  
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 2,092
Rep Power: 32 
Quote:
https://www.cfdonline.com/Wiki/Best...stage_analysis The drag values you are comparing to: where do they come from? Wind tunnel testing? The CFD model seems very simplified, missing lots of small parts a regular bicycle would have. And humans tend to have a different, more complex shape.
__________________
Please do not send me CFDrelated questions via PM 

September 12, 2019, 20:21 

#3 
New Member
Jacob Flanagan
Join Date: Sep 2019
Posts: 5
Rep Power: 2 
Thanks for the advice flotus1,
For the rotating wheel I also simulated a solution with no wheel rotation and just a no slip wall condition and the result was extremely similar. My drag estimations are based off research papers such as https://www.sciencedirect.com/scienc...67610518305762 which have a CdA of around 0.25 for a cyclist similar to mine. There are also other similar studies which have made the simplifications that I have made which get more accurate results. Would those minor simplifications cause that much decrease in drag? Thanks. 

September 13, 2019, 02:26 

#4  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 807
Rep Power: 22 
Quote:
It could you never know with turbulence. Here in this video where the tyres are rotated you could see the flow pattern at front tyre is much different than the back trye. Rotation seems to matter https://youtu.be/jAfJnUZAMco 

September 13, 2019, 03:56 

#5 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 4,894
Rep Power: 52 
At the first glance, I would advice that computing "accurately" the drag requires a very very fine grid close to the walls. Have you estimated your y+ distribution around the walls?


September 13, 2019, 04:37 

#6  
New Member
Jacob Flanagan
Join Date: Sep 2019
Posts: 5
Rep Power: 2 
Quote:
I did some rough y+ calculations but I wasn't 100% sure as the Reynolds numbers was hard to get right (not sure what the characteristic length should be) and then I also wasn't 100% sure what equation to use for the skin friction coefficient (Cf) used to determine the wall shear stress and therefore y+. Any advice on the y+, Re and Cf calculations would also be greatly appreciated. Thanks. 

September 13, 2019, 04:49 

#7  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 4,894
Rep Power: 52 
Quote:
Viscous drag prediction requires to use a wallresolution, that is at least 34 nodes having y+<=1. 

September 13, 2019, 04:56 

#8  
New Member
Jacob Flanagan
Join Date: Sep 2019
Posts: 5
Rep Power: 2 
Quote:
Thanks. 

September 13, 2019, 05:04 

#9  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 4,894
Rep Power: 52 
Quote:
You have to compute explicitly the values of y+, that requires to evaluate the u_tau velocity. An estimation of the bulk based Reynolds number is of O(10^6) therefore you need a well refined grid, a rough estimation would suggest you need to put the nodes at a distance lesser than 10^3 m. 

September 13, 2019, 05:09 

#10  
New Member
Jacob Flanagan
Join Date: Sep 2019
Posts: 5
Rep Power: 2 
Quote:
I am currently have my cell thickness at 5 x 10^(5)m (50 micrometers) which would suggest that the boundary layer would be resolved in that case? 

September 13, 2019, 05:50 

#11  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 4,894
Rep Power: 52 
Quote:
That depends not only on the node closest to the wall but on the description of the BL, that is on the next nodes normal to the wall. Again, you need to explicitly compute the y+ distribution around the body to ensure that is fully described by at least 34 nodes within y+. To do that you need to evaluate u_tau from your computation. 

September 13, 2019, 09:58 

#12 
Senior Member
Joern Beilke
Join Date: Mar 2009
Location: Dresden
Posts: 305
Rep Power: 14 
I would do a transient simulation at first.


Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculating drag coefficient from Fluent  Richard  FLUENT  12  April 11, 2018 15:35 
Problem in using parallel process in fluent 14  Tleja  FLUENT  3  September 13, 2013 10:54 
Drag force calculation in FLUENT  pychear  FLUENT  0  August 20, 2013 10:35 
direction vectors for lift and drag in fluent  nauman55  FLUENT  6  September 10, 2012 10:34 
How does FLUENT calculate lift and drag?  xTamx420  FLUENT  0  May 30, 2011 13:35 