,HM<^HX.^n i^^pjpEiNTIAL- ^' c^f°^ ^^«^^:^'^^'/^^?<^:>i;J;s'^SJi0 NACA :f RESEARCH MEMORANDUM SOME, EXAMPLES OF THE APPLICATIONS OF THE TRANSO^C AND SUPERSONIC AREA RULES TO THE # ^ # PREDICTION OF WAVE DRAG ^ ^ By Robert L. Nelson and Clement J. Welsh ^ ^ ^^U) - (x - xof (Al) CONFIDENTIAL l6 CONFIDENTIAL NACA RM L56DII Then, A = \ P^7p2r2(x) _ (x"- xo)2 dx (A2) Differentiating the expression for A gives dA 2 r^u (^ - ^0) dXQ n2 P ^^^ ^pVCx) - (x - Xq)' ^x (A3) where x, and x^ are the roots of x = xq - pR(x) and x = xq + pR(x), respectively. CONFIDENTIAL MCA RM L56DII CONFIDENTIAL 17 APPENDIX B THE NET FORCE ACTING ON THE OBLIQUE AREA OF A BODY OF REVOLUTION AT ZERO ANGLE OF ATTACK For a body of revolution at zero angle of attack, the net force is independent of roll angle. The derivation will be made for a roll angle of 90°. The net force in the 9 direction (the z-direction for a roll angle of 90°) (see sketch (b)) can be written as f Cp sin a dS /■,..2 . .„2 Since dS = Ifdy + dz and sin a dy \/dy2 + dz" The pressure coefficient at zero lift is a function of x only. The equation for y in terms of x is given by equation (l) of appen- dix A as follows: = Wp^'(x) - U -xo)' CONFIDENTIAL 18 CONFIDENTIAL MCA RM L56DII and p2R(x)g - (x - xo) dy - — -—=z=^=:^:^i=^=r dx Then, the net force can be written as K. I = (-2y^0 cp dy^ 4- (2 pO Cp dy /upper surface \ -^ /lower siirface •Xo S[p'^^-)i - (^ - -0)] ^ 2 rxo S[p'^(^)i - ^^ - ^)] dx + ^ / — 5= =- dx f 2 rxu ^p[p'^(^)i - (x - xo)| — = — / dx ^ '^^2 yp2H2(x) - (x - xo)2 Integrating the expression by parts gives XI dC^ If — ii is essentially constant between the limits of integrations, dx Then, from equation for the frontal projection of the oblique area in appendix A (eq. (A2)) -F dCr, t = _BA —2- q ^ dx CONFIDENTIAL MCA RM L56DII CONFIDENTIAL 19 APPENDIX C METHOD FOR DETERMINING WING-AREA DISTRIBUTION AND AREA-DISTRIBUTION-CURVE SLOPE This method assumes that the wing is thin and that the oblique Mach plane can be replaced by a plane perpendicular to the wing chord plane. The method also assumes the wing has straight leading and trailing edges and constant thickness ratio. The method is developed first for pointed-tip wings. Then, correc- tions are made for curved-wing — body junctirres and finite wing tips. In addition, the right- and left-hand wing panels are considered separately. Pointed-tip wings , lowing sketch: Consider the right-wing panel shown in the fol- max X* = Xq + y(p cos 9 - tan A) (c) CONFIDENTIAL 20 CONFIDENTIAL NA.CA RM L56DII The frontal projection of the area of one wing panel cut by the Mach plane is given by A = 2 Tz dy and — can be written as c <= ' ^ -max " 2 c ^^^ Then, for — constant. A = ^o^ol j^^^)^^'^ The value of t\ and v are related by the intersection line of the Mach plane and the wing chord plane for the right-wing panel given by the equation x' = Xp^ + y(p cos 9 - tan A) and for the left-wing panel by the equation x' = Xq + y(p cos e + tan A) Then, V = v„ + -^(p cos e ± tanA)T] With ^ = 1 - Tl ^0 ' V = — i — Vq + — (P cos e ± tanA)T] 1 - iL ^0 and Let Ti = vq - V (tan A ± p cos 9)-^ - v ^0 m = -T^ tan A ^0 CONFIDENTIAL NACA RM L56DII CONFIDENTIAL 21 In terms of tapered wing geometry, m is given by Then, Let m = A LJlA tan A 4 1 - A = ^0 - ^ "^ " Ji ± £_co^\ _ ^ \ tan A / K = m 1 + L^2L_i \ tan A ; for the left-wing panel and K = m 1 - P-^HlJ \ tan A for the right-wing panel. Then, ' K - V dT] K - Vq '^ ~~ (K - V) 2 c , K - V ^0 ' - K - V ^ a. = . OLlM dTi = — ^ ^ dv ^0 (K - v)5 The equation for the area can then be written as ^ ^ower (K - "y CONFIDENTIAL 22 CONFIDENTIAL NACA RM L56DII The slope of the ELrea-distribution curve is obtained by differentiating the expression for A dA _ t >(k - Vq) V, upper f(v) , _ ^ ' dv + ^o^.. (K - v)5 lower (k - vq)2 dv^^^3^ (^ - ^of '^Vper -.. (k-v N5 a ^^°"^^^"/K-v n5 dv, ^^PP^^' \ lower j ^ upper) or ^= ^0^0 |h(k,vo) Curves of G(K^vo) and H(K,vo) have been made up for a 65A series airfoil and are given in figures I6 and I7 for values of K from to 2.k and vq from to 1. In evalxiating / ^ — - — dv, f (v) was ^ (K - v)^ assumed to vary linearly between airfoil ordinate stations. Figure I8 gives a plot of f(v) for this ass\jmption. For K and Vq greater than 1, G(K,vo) and H(K, vq) are given by the expressions G(K,vo) = G(K,1) il_l_!2)! (K - 1)2 H(K,vo) =H(K,1) k:L-ml (K - 1) For K and Vq less than 0, G(K, vq) and H(K,vo) are given by the expressions g(k,vo) = G(K,0) (^ -J^f h(k,Vo) = H(k,0) ^^^.^ K CONFIDENTIAL NACA RM L56DII CONFIDENTIAL 25 Correction for curved-ving — body juncture .- The following sketch shows a pointed-tip wing mounted on a curved body. Line of intersection of Mach plane and wing chord plane The areas and slopes will be referred to the actual wing geometry (cp, s, and a). The areas and slopes, however, will be for the exposed pointed-wing tip. In sketch (d) consider one point of intersection of the wing panel with the body. The area of the wing cut by the Mach plane through this point is determined only by the product of SqC^ of the exposed wing through the point and the value of Vq for the exposed wing at the point of intersection. As the point of intersection changes, s 0^ '0^ and ^0 change and account for the intersection line. Expressed in terms of the actual wing-body characteristics. ^o'^o is given by ^0^0 I^[-(-^>iI The quantity — is related to v^ by the expression s u ^0 = (1 - A)m| 1 - (1 - A)| The area of the exposed wing panel cut by the Mach plane can be written as CONFIDENTIAL 2k CONFIDENTIAL NACA RM L56DII c^s 4 A=^[l- (1 -A)|] G(K,vo) The area is calculated for given value of Vq. The center-line value of V is given by Vr = V5[l - (1 - A)|] 1- K(l - A)| The slope is obtained by differentiating the expression for A and is given by ^ crs i[l - (1 - A)|j|[l - (1 - A)|]h(k,vo) - 2(1 - A)^ g(k,vo) dVv If — §• = 0, dVQ 1 - A 1 - (1 - A)^+ (k - vo)(l -A) if dVr dA dv. = rTf[i- (1 - ^)flH(K,vo) Correction for finite vlng tip .- In order to correct the pointed- tip wing panel and slopes for the finite wing tip, the areas and slopes outboard of the wing tip are subtracted. Intersection line of Mach plane and wing From sketch (e) (e) s„c^ = c^s A'^ ^ 1 _ A CONFIDENTIAL NACA RM L56DII CONFIDENTIAL 25 Then the areas and slopes are given by Atip = -^ _ ^ g(k,vo) dAtip c^s|a The center-line value of v is given by v^ = VqA + K(l - A) CONFIDENTIAL 26 CONFIDENTIAL NACA RM L56DII REFERENCES 1. Whitcomb, Richard T.: A Study of the Zero-Lift Drag-Rise Character- istics of Wing-Body Combinations Near the Speed of Sound. NACA RM L52H08, 1952. 2. Hall, James Rudyard: Comparison of Free-Flight Measurements of the Zero-Lift Drag Rise of Six Airplane Configurations and Their Equivalent Bodies of Revolution at Transonic Speeds. NACA RM L53J21a, 195^. 3. Jones, Robert T.: Theory of Wing-Body Drag at Supersonic Speeds. NACA RM A53Hl8a, 1953. k. Holdaway, George H.: Comparison of Theoretical and Experimental Zero- Lift Drag-Rise Characteristics of Wing-Body-Tail Combinations Near the Speed of Sound. NACA RM A55H17, 1953. 5. Alksne, Alberta: A Comparison of Two Methods for Computing the Wave Drag of Wing -Body Combinations. NACA RM A55A06a, 1955. 6. Lomax, Harvard: The Wave Drag of Arbitrary Configurations in Line- arized Flow As Determined by Areas and Forces in Oblique Planes. NACA RM A55A18, 1955. 7. Sears, William R.: On Projectiles of Minimum Wave Drag. Quarterly Appl. Math., vol. IV, no. k, Jan. 19^7, pp. 36I-366. 8. Staff of the Computing Section, Center of Analysis (Under Direction of Zdenek Kopal): Tables of Supersonic Flow Around Cones. Tech. Rep. No. 1, M.I.T., 19^7. 9. Fraenkel, L. E.: The Theoretical Wave Drag of Some Bodies of Revolu- tion. Rep. No. Aero. 2^4-20, British R.A.E., May 1951. 10. Van Dyke, Milton D.: Application of Hypersonic Small-Disturbance Theory. Jour. Aero. Sci., vol. 21, no. 3, Mar. 195^4-, pp. 179-186. 11. Van Dyke, Milton D.: Practical Calculation of Second-Order Super- sonic Flow Past Nonlifting Bodies of Revolution. NACA TN 2^hk, 1952. 12. Hart, Roger G., and Katz, Ellis R. : Flight Investigations at High- Subsonic, Transonic, and Supersonic Speeds To Determine Zero-Lift Drag of Fin-Stabilized Bodies of Revolution Having Fineness Ratios of 12.5^ 8.91j and 6.04 and Varying Positions of Maximum Diameter. NACA RM L9I3O, 19^+9. CONFIDENTIAL MCA RM L56DII CONFIDENTIAL 27 13. Wallskog, Harvey A., and Hart, Roger G. : Investigation of the Drag of Blunt -Nosed Bodies of Revolution in Free Flight at Mach Numbers From 0.6 to 2.3. NACA RM L53D1W, I953. ih. Chapman, Dean R., and Kester, Robert H.: Turbulent Boundary-Layer and Skin-Friction Measurements in Axial Flow Along Cylinders at Mach Numbers Between O.5 and 3.6. NACA TN 3097, 1954. 15. Stevens, Joseph E., and Purser, Paul E.: Flight Measurements of the Transonic Drag of Models of Several Isolated External Stores and Nacelles. NACA RM L5I+LO7, 1955. 16. Morrow, John D., and Nelson, Robert L. : Large-Scale Flight Measure- ments of Zero-Lift Drag of 10 Wing-Body Configurations at Mach Numbers From 0.8 to 1.6. NACA RM L52Dl8a, I953. 17. Wallskog, Harvey A., and Morrow, John D. : Large-Scale Flight Measure- ments of Zero-Lift Drag and Low-Lift Longitudinal Characteristics of a Diamond-Wing — Body Combination at Mach Numbers From O.725 to 1.5^. NACA RM L53C17, 1953- 18. Welsh, Clement J., Wallskog, Harvey A., and Sandahl, Carl A.: Effects of Some Leading-Edge Modifications, Section and Plan-Form Variations, and Vertical Position on Low-Lift Wing Drag at Transonic and Super- sonic Speeds. NACA RM L54K01, 1955. 19. Sandahl, Carl A., and Stoney, William E.: Effect of Some Section Modifications and Protuberances on the Zero-Lift Drag of Delta Wings at Transonic and Supersonic Speeds. NACA RM L53L2i^-a, 195^. 20. Holdaway, George H., and Mersman, William A.: Application of Tchebichef Form of Harmonic Analysis to the Calculation of Zero-Lift Wave Drag of Wing-Body-Tail Combinations. NACA RM A55J28, 1956. 21. Beane, Beverly: The Characteristics of Supersonic Wings Having Biconvex Sections. Jour. Aero. Sci., vol. I8, no. 1, Jan. 1951? pp. 7-20. CONFIDENTIAL 28 COKFIDEOTIAL MCA RM L56DII (a) e = 90*= (b) e = 0°. Flg-ure 1.- The areas and pressures which influence the drag of configura- tions at supersonic speeds. COKFIDENTIAL NACA RM L56DII CONFIDENTIAL 29 (a) Low-flneness-ratlo body; — = 6.04. d /.2 (b) High-fineness-ratio body; — = 10. d Figvire 2.- The effect of Mach number on the area-distribution-curve slope of bodies of revolution. CONFIDENTIAL 50 CONnDEOTIAL NACA RM L56DII <9 M 5" /. 15- 4.68 D 1.5 /. 07- 3.94 10 1.12-3.01 A /2.5 /. 16-2.47 t^ /5 7.22-2. 12 G^ 17.5 728-7.88 Q 20 /. 33-7. 70 .5V .4 - .3 - 7ir?e of perfe>ct agreemerff — v^ /q .2 .3 .4 Cff 7'£xact f/?eory) Figi:ire 3.- Comparison of the drag of cones calculated with the area rule with the drag from exact theory (ref . 8) . CONFIDENTIAL CONFIDENTIAL 51 £xactt/?eory/'fi = 5''J Area- rule t^eort/ 5/e/7der-boofi/ tJ76>ori/ .2 .4 . .6 .d /.O Figure \.- The drag of cones in collapsed form. Ci)UlrJ Slender-body t/?eory Area-ru/e f/7eory Area- rule ^J7eory "^ pressure term O Secorfd -order theory J I L 4 , 6 .3 J.O Figure 5-- The drag of parabolic noses in collapsed form. CONFIDENTIAL 32 CONFIDENTIAL NACA RM L56DII I 'I !fe » 5^ V N QJ ^ -p A -P t Td (D +3 . 03 >> H U Ti 0) H ^ tU +3 >i CO TJ (U •H £1 • Ti 1 CO ^ • ja Q) Ti fU tIJ •H 0) § CQ XI •^ cd -ri MD ^ fi cd S 6 ft -P •\ QH £3 -4- (U • 6 bO -H 03 Jh II fl ro X -3 0) 0) s: -p ^ * s -p cd ^•^ S3 (U g^ •H fH u oj oj n 0) e M 03 1 1 ■<6 H) ^ M •H p^ CONFIDENTIAL NACA RM L56DII CONFIDENTIAL 55 l.0r- / / \ y f /^ Area-ru/e /heory y^ 1 U/t -err' — -D M .0/ " 7 III III /.O /./ ^^ /.2 /.3 /.4 /.S /.(, (f) Model C-5 (ref. l8); A = 3; A = 0.2; Ac/J+ = ^5°; NACA 65AOO3 air- foil section; — = O.O606. Sw Exper/ment Area-ru/e iheory £quiya/er>t body (g) Model 6 (ref. 16); A = i^-; A = 0.6; A^/ij. = ^5°; NACA 65AOO6 airfoil Sb section: — = 0.0606. Sw Figure 11.- Continued. CONFIDENTIAL 1^0 CONFIDENTIAL NACA RM L56DII <^o Area-ru/e t/Jeory £qu/^a/enf body (h) Model 2 (ref. 16); A = 5.0^+; A = 0.59^; Ajc/if = 0°; - = 0.0^5; .04 ^i> .02 — = 0.0606. Sw Figure 11.- Concluded. Exppr/znent /.S /J> Figure 12.- Comparison of equivalent-body drag and configuration drag for an airplane configuration. CONFIDENTIAL NACA RM L56DII CONFIDENTIAL kl /6 14 /2 /O (t/£)2far,f- 8 ./ ^rea - ru/e theori^ ■3 .4 /3 ^ane .& .8 Figure 13 •- Comparison of the drag of delta wings calculated with two versions of the linearized theory. CCNFIDENTIAL k2 CONFIDEFriAL MCA RM L56DII .a/6 .0/2 Co .008 .004 - Figure 1^4-.- Effect of wing-panel separation on wing drag. 6o° delta wing; NACA 65AOO3 airfoil section. .0/6 .0/2 ^D .008 .004- Figure 15-- Comparison of the sum of component drags with the configxira- tion drag. 6o° delta wing; NACA 65AOO3 airfoil section. CONFIDENTIAL NACA RM L56DII CONFIDENTIAL h-i .S2r Figure 16.- Area distribution parameter G(K, vq) for 65A series airfoil. COI^i'IDENTIAL kk CONFIDEETIAL MCA RM L56DII (a) vq from to 0.^4-5. Figure 17-- Area-distribution- slope parameter H(K, vq) for 65A series airfoil. CONFIDENTIAL NACA RM L56DII CONFIDENTIAL h^ (b) Vq from O.5 to 1.0. Figure 17-- Concluded. CONFIDENTIAL 46 CONFIDENTIAL NACA EM L56DII O Airfo// ord/haie sfcrf/^/?s^ /^v) Figure I8.- Approximation of 65A series airfoil for the calculation of G(K,Vo) and H(K,Vo). CONFIDENTIAL NACA - Lajigley Field, Va. CONFIDENTIAL UNIVERSITY OF FLORIDA 3 1262 08106 568 1 UN!\/ERSITY OF FLORIDA DOCUMENTS DEPARTMENT 1 20 MARSTON SCIENCE UBRARY P.O. BOX 117011 GAINESVILLE, PL 32511-7011 USA CONFIDENTIAL