
[Sponsors] 
how to calculate drag force coefficient of cylinder in oscillating flow 

LinkBack  Thread Tools  Display Modes 
August 24, 2014, 09:59 
how to calculate drag force coefficient of cylinder in oscillating flow

#1 
New Member
congvu
Join Date: Apr 2012
Posts: 12
Rep Power: 7 
Hello,
I try to calculate the drag force coefficient on cylinder in case with oscillating flow at inlet boundary. i.e U=Uosin(2pi*f*t) How to set velocity in Reference values? any ideas ? Thank you 

September 17, 2014, 07:32 

#2 
Senior Member
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 6 
Hi,
I'm also looking for the same question. Could you find the answer? Thanks and regards. Maimouna 

September 23, 2014, 01:06 

#3 
New Member
congvu
Join Date: Apr 2012
Posts: 12
Rep Power: 7 
Sorry, I did not figure out yet. How about you?


September 23, 2014, 04:27 

#4 
Senior Member
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 6 
Sorry, not yet too. Which type of oscillating are you doing: forced or natural?


September 23, 2014, 04:41 

#5 
New Member
congvu
Join Date: Apr 2012
Posts: 12
Rep Power: 7 

September 23, 2014, 04:52 

#6 
Senior Member
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 6 
Unfortunately, I don't know what you exactly mean, if you are interested you case google ''how to accurately flow around cylinder'' there are some attached folders might help.
All the best. 

September 23, 2014, 05:11 
Help

#7 
New Member
yes
Join Date: Sep 2014
Posts: 2
Rep Power: 0 
Hello everyone
If anyone can help me to simulate the problem attached using Ansys, and if you can explain how I can treat with the roughness parameter, please I need it in my project, and that will be so helpful. I have been working on it too may time but it gave me a different results than the analytic solution especially with the friction part. Please if you can contact me with my Email that will be so nice of you. Haider_skran@yahoo.com 2.jpg 

September 23, 2014, 05:59 

#8  
Senior Member
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 996
Rep Power: 17 
Quote:
I think you need to take care on 2 things in this problem: 1 you have a developed flow: this means the fluid is fully developed, so in fluent you need to write an udf to have a fully developed flow at the inlet (there is an example in udf manual, velocity inlet + k + epsilon (<turbulent dissipation)); otherwise, you can make the pipe longer, monitor the velocity profile in the axial direction, identify where the fluid fully develops, and start your calculations (pressure drop, etc.) from that point. Outlet boundary should be set to outflow, because the fluid is fully developed. 2 Roughness parameter: the Moody's chart was obtained by coating the internal surfaces of tested pipes with a monolayer of sand, and the pipe wall roughness was defined as the average diameter of the sand grains (epsilon); the sentence "the roughness of stainless steel is 0,002 mm" is quite ambiguous: which roughness? I think this refers to that of the Moody's chart, so you can input directly this value into the fluent "roughness heigh" box (take care to units). Roughness constant is a parameter of the uniformity of sand grains: a value of 0,5 means all the grains are of equal diameter, so I think you should set this value to 0,5. Several algorithms were developed to link various roughness to the sand grain roughness, I usually use the followings: epsilon=Ra*5,863 epsilon=Rrms*3,100 epsilon=Rzd*0,978 epsilon=sand grain roughness Ra=average roughness (this is the average value of the heights, + or , measured from a mean line of a profile) Rrms=root mean squared roughness Rzd=peak to valley roughness Daniele 

September 23, 2014, 07:45 

#9  
Senior Member
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 996
Rep Power: 17 
Quote:
You should investigate about how to calculate the drag force; Maybe you could use:  for pressure forces: F_P(f,t1) * A[0] where F_P is the face pressure and A[0] is the projected area in x direction;  for viscous forces: F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR)[0] this is not documented in udf manual; Then you should add: Domain *d = Get_Domain(1); to get domain pointer, which should be 1, as you have only 1 phase; Thread *t1 = Lookup_Thread(d,xxx); where xxx is the id of the wall of the cylinder (you read it in the fluent gui) Once you have the total force (pressure+viscous), you can easily compute the drag coefficient with time varying velocity. Daniele 

September 25, 2014, 04:53 

#10 
New Member
yes
Join Date: Sep 2014
Posts: 2
Rep Power: 0 
First of all I would like to thank you Daniele
I have done what you said , but I did not get the results that I looking for , If you can simulate the problem that will be so helpful 

September 25, 2014, 04:59 

#11 
Senior Member
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 996
Rep Power: 17 
Rayder,
you have 2 posts on cfdonline, I will not simulate your problem and nobody on this forum will do it for you. Describe what you have done, you didn't describe anything about your setup, and I will help you with suggestions. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Drag coefficient too high at flow around a cyclinder  Gunni  OpenFOAM Running, Solving & CFD  16  October 19, 2015 19:21 
how to calculate darg force on cylinder in laminar flowCOMSOL  newmax  Main CFD Forum  0  January 24, 2014 07:28 
Drag coefficient of flow past cylinder vs time  pedroxramos  FLUENT  0  January 14, 2013 13:39 
drag force in flow over a cylinder  student  Main CFD Forum  1  December 13, 2008 17:59 
Lift force coefficient for 2phase flow  Summer  FLUENT  0  April 29, 2008 17:58 