# Axisymmetry problem and XFOIL

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

 February 26, 2008, 16:02 Axisymmetry problem and XFOIL #1 JC Guest   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

 August 27, 2011, 02:13 Xfoil hack for drag of axisymmetric bodies of revolution #2 New Member   Keith Pickering Join Date: Aug 2011 Posts: 1 Rep Power: 0 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 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: COMMON/CR09/ ADEG,ALFA,AWAKE,MVISC,AVISC, & XCMREF,YCMREF, & CL,CM,CD,CDP,CDF,CL_ALF,CL_MSQ,CDX,FRS, & PSIO,CIRC,COSA,SINA,QINF, & GAMMA,GAMM1, & MINF1,MINF,MINF_CL,TKLAM,TKL_MSQ,CPSTAR,QSTAR, & CPMN,CPMNI,CPMNV,XCPMNI,XCPMNV 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: WRITE(*,2020) ALFA/DTOR, CL, CM, CD, CDF, CDP 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

August 27, 2011, 03:16
#3
Senior Member

Join Date: Jan 2011
Posts: 247
Blog Entries: 5
Rep Power: 9
Quote:
 Originally Posted by KAP The other terms of the equation are for the wake momentum and don't change. KAP
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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Luk_Fiz CFX 6 June 18, 2010 07:26 Bogey Jammer Main CFD Forum 0 September 29, 2009 17:06 parvez2683 FLUENT 0 July 30, 2009 01:10

All times are GMT -4. The time now is 13:21.