CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Flow Past a Cylinder with Non-Zero Average Lift Coefficient

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

Like Tree1Likes
  • 1 Post By innoaloe

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 21, 2020, 22:21
Default Flow Past a Cylinder with Non-Zero Average Lift Coefficient
  #1
New Member
 
Innocentio Loe
Join Date: May 2020
Posts: 4
Rep Power: 2
innoaloe is on a distinguished road
Hello,

I've been working on simple flow past a cylinder simulation with Re=100.
Boundary condition is flat velocity profile on the inlet, 0 on the outlet, and slip top and bottom wall.
The domain is 20d in width and 10d in height, with the cylinder centered at the origin, 5d to both left and bottom wall

I'm trying to validate the result, compared to literature data.
Expected lift and drag coefficient (CL and CD) around the cylinder wall are CL = +/-0.33 and CD = 1.34 +/-0.009

The simulation is working fine. However, what I notice is that even though my CL oscillates at the correct amplitude, it is not oscillating around 0 average but instead has small negative average value instead.
Also my CD is oscillating around 1.42 instead, so I guess there is the same DC drift happening as well here.

For the CL, it is easy to remove this additional drift, but for CD I cannot really do so, since there is no baseline of judgement.

My question is whether this is commonly encountered? If so, how to validate such a case then? I'm using COMSOL Multiphysics at the moment. Might try to compare with OpenFOAM later on.
innoaloe is offline   Reply With Quote

Old   May 22, 2020, 04:07
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,424
Rep Power: 58
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Have you checked that the grid is perfectly symmetric?
FMDenaro is offline   Reply With Quote

Old   May 22, 2020, 07:41
Default
  #3
New Member
 
Innocentio Loe
Join Date: May 2020
Posts: 4
Rep Power: 2
innoaloe is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Have you checked that the grid is perfectly symmetric?
Hi, thank you for your input. Certainly I haven't thought about the mesh symmetry. Also, currently I am using an unstructured mesh which maybe will make it harder to create symmetrical mesh. I will try using a structured mesh.
innoaloe is offline   Reply With Quote

Old   May 24, 2020, 06:20
Default
  #4
Senior Member
 
Gerry Kan's Avatar
 
Gerry Kan
Join Date: May 2016
Posts: 224
Rep Power: 7
Gerry Kan is on a distinguished road
Dear innoaloe:

One thing you should also check is whether the CL oscillates about some value as the solver iterates. It could be that the solution is oscillating, as it happens physically. Look at your solution field to verify.

Gerry

PS I am sorry I did not read your entire message. I would second with FMDenaro and use a perfectly symmetric mesh. I don't know if you can do this natively with COMSOL but I would be surprised if you couldn't. Otherwise you need a very fine mesh to attenuate any asymmetric effects.
Gerry Kan is offline   Reply With Quote

Old   May 24, 2020, 08:19
Default
  #5
New Member
 
Innocentio Loe
Join Date: May 2020
Posts: 4
Rep Power: 2
innoaloe is on a distinguished road
Quote:
Originally Posted by Gerry Kan View Post
Dear innoaloe:

One thing you should also check is whether the CL oscillates about some value as the solver iterates. It could be that the solution is oscillating, as it happens physically. Look at your solution field to verify.

Gerry

PS I am sorry I did not read your entire message. I would second with FMDenaro and use a perfectly symmetric mesh. I don't know if you can do this natively with COMSOL but I would be surprised if you couldn't. Otherwise you need a very fine mesh to attenuate any asymmetric effects.
Dear Gerry,
thank you for your input as well. I did as FMDenaro's suggested to create symmetrical mesh, also ensuring the mesh around the cylinder to be more refined. And yes, it is possible to do so in COMSOL.

I managed to get an acceptable result by doing so. There is still a very tiny average value on the oscillation, but the value is negligible if compared with the CL oscillation itself, so in my opinion such small numerical error is to be expected.

Regards
Gerry Kan likes this.
innoaloe is offline   Reply With Quote

Reply

Tags
comsol, drag coefficient, lift and drag, lift coefficient

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
Simulating Flow past Circular Cylinder parthiv1991 OpenFOAM Running, Solving & CFD 4 June 3, 2018 06:19
Drag force coefficient too high for a flow past a cylinder using komega sst Scabbard OpenFOAM Running, Solving & CFD 37 March 21, 2016 17:16
Drag Coefficient for flow past a Cylinder o_mars_2010 Main CFD Forum 0 April 18, 2013 07:17
How to calculate the drag coefficient for flow past cylinder o_mars_2010 Tecplot 0 April 18, 2013 02:26
Drag coefficient of flow past cylinder vs time pedroxramos FLUENT 0 January 14, 2013 13:39


All times are GMT -4. The time now is 03:37.