CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Axisymmetry problem and XFOIL

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

LinkBack Thread Tools Display Modes
Old   February 26, 2008, 16:02
Default Axisymmetry problem and XFOIL
Posts: n/a
Dear All I would like to get the boundary layer evolution around a rotationally symmetrical body (kind of fuselage). The problem is axisymmetric, so that I can plot the shape of the fuselage in a 2d plan. However, does XFOIL is able to give an answer for this problem? The equation used by the code are still valid for axisymmetric shapes? Thanks
  Reply With Quote

Old   August 27, 2011, 02:13
Default Xfoil hack for drag of axisymmetric bodies of revolution
New Member
Keith Pickering
Join Date: Aug 2011
Posts: 1
Rep Power: 0
KAP is on a distinguished road
Here's how I did this. First, you need a version of Xfoil compiled from source. That's because you're going to change the code and recompile it. I'm also assuming you've got a Linux system, where make commands work. In Windows, you're on your own for compiling. And finally, of course, I'm assuming that the airfoil you loaded is symmetrical.

Assuming you've got that, navigate to the Xfoi/src directory.

1. In xfoil.f, find the line (about line 1178 in my version, 6.97) that computes CD. It reads like this:

CD = 2.0*THWAKE * (UEWAKE/QINF)**(0.5*(5.0+SHWAKE))

After that line (and before the ELSE) insert the following lines:

C------ compute axisymmetric drag coeff CDX for body of revolution in a
C similar manner, using frontal area FRS as reference area.
C First term is now area of wake as defined by theta.
FRS = PI*(THICKB*0.5)**2
CDX = 2.0*PI*THWAKE**2 * (UEWAKE/QINF)**(0.5*(5.0+SHWAKE)) / FRS

Theta is the boundary layer momentum thickness and THWAKE is theta in the downstream wake. For axisymmetric flow, instead of multiplying by 2 (which you do in 2-d to get total thickness of the wake) we use the area of circle to get the AREA of the wake, and convert to Cd by dividing by the cross sectional area FRS. The other terms of the equation are for the wake momentum and don't change. Easy, no? [Reference: J.L. Hess & R.M. James, "On the problem of shaping an axisymmetric body to obtain low drag at large Reynolds numbers," McDonnel Douglas report MDC J6791, (1975).]

2. In XFOIL.INC define the two new variables you just used. You can really do this
anywhere as long as you define them as real. To be consistent, I did this in the same
line where CD is defined (just added CDX and FRS to the line), like this:


3. In xoper.f, output your result. Find the line where results are written to the screen (not the plot) which is about line 2704 in my version. It looks line this:


Below this, add the line:

WRITE(*,2025) CDX

Then add the format line you just referenced. I put it just below line label 2020, though it really can go anywhere in subroutine VISCAL

2025 FORMAT ('Axisymmetric drag CDx =',F9.7 )

Finally, save all files and compile. In directory /Xfoil/bin issue a "make xfoil" to compile it all up.

KAP is offline   Reply With Quote

Old   August 27, 2011, 03:16
Senior Member
truffaldino's Avatar
Join Date: Jan 2011
Posts: 248
Blog Entries: 5
Rep Power: 10
truffaldino is on a distinguished road
Originally Posted by KAP View Post
The other terms of the equation are for the wake momentum and don't change.
Could you explain why the momentum thikness does not change?
Indeed, for an inviscid solution surface velocities on 2d airfoil and axisymmetric body with the same profile will be quite different, so bl equations too.

Is this explained in the reference given above?

truffaldino is offline   Reply With Quote


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
Pressure coeff for aerofoils problem Luk_Fiz CFX 6 June 18, 2010 07:26
[FloWorks] Request advice for an airfoil calculation problem Bogey Jammer Main CFD Forum 0 September 29, 2009 17:06
3D to 2D problem parvez2683 FLUENT 0 July 30, 2009 01:10

All times are GMT -4. The time now is 05:52.