CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Matching cylinder drag

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 2, 2010, 02:16
Default Matching cylinder drag
  #1
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
Hello all,

I am beginning to try to match the drag coefficients seen here:
http://scienceworld.wolfram.com/phys...inderDrag.html

My case setup for 1e5 Reynolds number can be downloaded from here:
http://files.me.com/ben.racine/0g7e83

All you need to do still is unzip and execute "simpleFoam > log.simpleFoam"
because the grids are already made with blockMesh and snappyHexMesh.

After varying Reynolds I get the following results... where the columns are as defined in the Wolfram link above...

Re Drag / L Drag (N?) C_D
1 0.0 0.00007 33091
10 0.0 0.00687 32621
100 0.0 0.00617 293
1000 0.0 0.06170 29
10000 0.1 0.62 2.95
100000 1.5 11.91 0.57
233333 6.3 50.71 0.44
1000000 106.2 849.85 0.40
10000000 10060.3 80482.00 0.38

I don't find to this very strong agreement with the Wolfram reference.

Any assistance is greatly appreciated.

Regards,
Ben Racine
i3enhamin@gmail.com

Last edited by bjr; March 2, 2010 at 03:35.
bjr is offline   Reply With Quote

Old   March 2, 2010, 03:24
Default I've attached some images...
  #2
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
I seem to be getting good results when it is inertia/pressure dominated and rather terrible results when it's viscous dominated... probably due to my very poorly resolved boundary layer mesh?
Attached Images
File Type: jpg 1e2Re_slice_pressure_contours.jpg (65.9 KB, 32 views)
File Type: jpg 1e7Re_slice_pressure_contours.jpg (63.0 KB, 22 views)
File Type: jpg 1e2Re_slice_velocity_contours.jpg (59.6 KB, 23 views)
File Type: jpg 1e7Re_slice_velocity_contours.jpg (61.0 KB, 23 views)
File Type: jpg data.jpg (45.2 KB, 28 views)
bjr is offline   Reply With Quote

Old   March 2, 2010, 05:30
Default
  #3
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Ben

I have looked at your setup, and two things might cause the problems:

1. Looking at the discretization of the cylinder at the upper and lower wall they are completely different, and it appears that the mesh lines have a twist. Why are you making 3D computations when using simpleFoam, would that not be rather waist? The difference in mesh gives a non-uniformity in the direction of the cylinder.

2. As you pointed out yourself, the mesh resolution is rather coarse, and my best guess is that you will not be able to model the boundary layer and the boundary layer separation correctly using such a coarse mesh. Hence the prediction of drag coefficient will be difficult.

Best regards,

Niels
ngj is offline   Reply With Quote

Old   March 2, 2010, 05:41
Default
  #4
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
Hi Niels,

1. This can be corrected by changing the line 43 of snappyHexMesh from

" max (4.0 0.7 2.5)"
to
" max (4.0 0.7 8.0)"

2. I've noticed that I am getting reasonable pressure drag. It seems the viscous drag is the suspicious part. This all seems to hang together that it's boundary layer resolution. But, I saw on the external aerodynamics thread...

Automotive test case

... some talk of turbDragCoefficient... anybody know if that needs to be addressed in conjunction?

Regards,
Ben Racine
bjr is offline   Reply With Quote

Old   March 2, 2010, 15:45
Default
  #5
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
I think a relevant follow up question is... has anyone had much success setting up low Re grids with snappyHexMesh... i.e. with adequate near-body boundary layer resolution?
bjr is offline   Reply With Quote

Old   March 2, 2010, 16:28
Default
  #6
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
files.me.com/ben.racine/ny44nk contains an updated case setup... again just issue "simpleFoam > log.simpleFoam"
bjr is offline   Reply With Quote

Old   March 3, 2010, 03:39
Default
  #7
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Ben

It is fairly easy to make a mesh with good mesh resolution near the cylinder surface using blockMesh. I have made some tests using that in the past and it works quite well. My suggestion is as you have problems obtaining sufficient resolution using snappyHexMesh, then try to use blockMesh and see the difference.

Best regards,

Niels
ngj is offline   Reply With Quote

Old   March 3, 2010, 21:20
Default
  #8
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
Fair enough... for academic investigations... but, I'm hoping to find an algorithmic tet or hex-mesher that is also capable of creating a boundary layer. I'm half tempted to write the boundary layer part myself and then let one of these other tools treat that surface as the body itself... and combine the results.
bjr is offline   Reply With Quote

Old   March 16, 2010, 13:47
Default
  #9
bjr
Member
 
Ben Racine
Join Date: Mar 2009
Location: Seattle, WA, USA
Posts: 62
Rep Power: 8
bjr is on a distinguished road
Send a message via AIM to bjr Send a message via Skype™ to bjr
Turns out I didn't know about the Low Re version of the k-omega SST model. This was all failing at low Re and succeeding at high Re. There is a thread about the low Re version "out there".
bjr is offline   Reply With Quote

Reply

Tags
cylinder, drag coefficient, openfoam, simplefoam, validation

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Drag Coefficient Verification around Cylinder MFerris CFX 44 May 11, 2015 18:27
Coefficient DRAG and LIFT - Cylinder 3D Jwolf CFX 11 February 4, 2010 06:08
Drag on non-infinite cylinder between parallel walls deltaquattro Main CFD Forum 0 September 10, 2009 13:01
Solver and geometry choose for drag coefficient calculation around circular cylinder at large Re lin OpenFOAM Running, Solving & CFD 3 April 16, 2009 10:50
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


All times are GMT -4. The time now is 01:44.