Copy 312/ RM A54F28 o NACA RESEARCH MEMORANDUM ON THE RANGE OF APPLICABILITY OF THE TRANSONIC AREA RULE By John R. Spreiter Ames Aeronautical Laboratory Moffett Field, Calif. DOCUMENTS DEPARTMENT lOW"*" • — ^ 1 20 MARSTON SCIENCE UBRARY ^g^cB NO. 30?6 RO. BOX 11701) GAINESVILLE, FL 32611-7011 USA UNIVERSITY OF FLORIDA jteoRITY: J--'- CROWL^ dat: CLASSIFIED DOCUMENT This material contains information affecting the National Defense of the United States within the meani n g of the espionage laws, Title 18, U.S.C., Sees. 793 and 794, the transmission or revelation of which in any manner to an unauthorized person is prohibited by law. NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS WASHINGTON August 23, 1954 t^- — — — imm^ i» CONFIDENTIAL^ NACA RM A5UF28 CONFIDENTIAL NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS RESEARCH MEMORANDUM ON THE RANGE OF APPLICABILITY OF THE TRANSONIC AREA RULE By John R. Spreiter SUMMARY Some insight into the range of applicability of the transonic area rule has been gained by comparison with the appropriate similarity rule of transonic flow theory and with available experimental data for a large family of rectangular wings having NACA 63AXXX profiles. In spite of the small, number of geometric variables available for such a family, the range is sufficient that cases both compatible and incompatible with the area rule are included. INTRODUCTION A great deal of effort is presently being expended in correlating the zero-lift drag rise of wing-body combinations on the basis of their streamwise distribution of cross-section area. This work is based on the discovery and generalization announced by Whitcomb in reference 1 that "near the speed of sound, the zero-lift drag rise of thin low-aspect-ratio wing-body combinations is primarily dependent on the axial distribution of cross-sectional area normal to the air stream." It is further conjec- tured in reference 1 that this concept, known as the transonic area rule, is valid for wings with moderate twist and camber. Since an accurate pre- diction of drag is of vital importance to the designer, and since the use of such a simple rule is appealing, it is a matter of great and immediate concern to investigate the applicability of the transonic area rule to the widest possible variety of shapes of aerodynamic interest. The experimental data contained in reference 1 and many subsequent papers have shown that this simple rule is often remarkably successful for a wide variety of shapes ranging in complexity from simple bodies of revo- lution to models of complete airplanes. Furthermore, important reductions in the transonic drag of wing-body combinations have been realized by indenting the body so that the axial distribution of cross-section area corresponds to that of smooth bodies of revolution having low drag at CONFIDENTIAL CONFIDENTIAL NACA RM A5*4-F28 supersonic speeds. On the other hand, it is important not to overlook the fact that there are a number of test results on equivalent bodies for which the correlation of drag rise by the area rule is unsatisfactory. Inasmuch as the models tested are generally of complex geometry, and only the original model and an equivalent body of revolution are tested, it is difficult to ascertain whether these discrepancies are attributable to viscous phenomena or to the fact that the drag rise may depend on other geometric parameters than the axial distribution of cross-section area. It is the purpose of this note to examine in further detail the applicability of the area rule. Despite the fact that most of the emphasis in the tests relative to the area rule has been on wing -body combinations, there exists such a scarcity of experimental data of a sufficiently sys- tematic type that the present discussion will be confined to rectangular wings without bodies. Transonic drag data are available from bump tests in the Ames l6-foot high-speed wind tunnel for a large family of rectangu- lar wings having NACA 63AXXX sections, aspect ratios varying from -0.5 to 6.0, thickness ratios from 2 to 10 percent, and both symmetrical and cam- bered profiles. These results are reported in references 2, 3, and k and have been studied by McDevitt (refs. 3 and 5) who showed in a convincing manner that the experimental data can be correlated successfully by means of the transonic similarity rules. These same results will be used herein to evaluate one phase of the transonic area rule. The relationship between the two rules is naturally of interest and will also be explored. Since the transonic area rule is considered to apply equally to all low-aspect-ratio wing-body combinations, detailed examination of such a limited class of aerodynamic shapes as a family of rectangular wings is not without value inasmuch as limitations revealed in special applications must appear as a limitation in the general case. The restricted range of the investigation is compensated somewhat by the fact that the geometric simplicity increases the chances of understanding the underlying causes. Although the results can only be said to apply with surety to the specific cases investigated, the method of approach is not restricted and may be applied similarly to other cases as more data become available and as understanding increases. In this way, the present discussion may be con- sidered more as suggestive than definitive. PRINCIPAL SYMBOLS A aspect ratio A [(7+1) M 2 t] 1/3 A b wing span CONFIDENTIAL NACA RM A5UF28 CONFIDENTIAL C D^ °w <l s p D o AC Dn A C L c 2- c D D o D, ^w qc 2 3 1 ^o qSp qs p [Mo 2 ( 7 + I)]" L/3 T 2/3 C L ideal lift coefficient wing chord drag drag at zero lift wave drag at zero lift /D G Vqc 2 'Mo /D D N W 2 , Mref function of indicated variables dimensionless function describing thickness distribution g 2 dimensionless function describing camber distribution H h c h maximum camber of wing h t h L M Mach number M re f reference Mach number less than the critical q dynamic pressure, -£■ U Q 2 CONFIDENTIAL h CONFIDENTIAL NACA RM A5^F28 R Reynolds number ^ref reference Reynolds number corresponding to M re f S area S c cross-section area Sm maximum cross-section area Sp plan-form area s dimensionless area distribution function t maximum thickness U Q free-stream velocity x,y Cartesian coordinates in plane of wing where x extends in the direction of the free-stream velocity Y dimensionless function describing plan form ordinates of wing profiles Zy ordinates of upper surface of wing profiles Zj ordinates of lower surface of wing profiles a angle of attack a *o a t/c ratio of specific heats, for air 7 = 1.1+ Mq 2 [(7 + 1) M^tF' 3 P density t c CONFIDENTIAL NACA RM A54F28 CONFIDENTIAL DIMENSIONAL CONSIDERATIONS The forces on a body moving through an infinite fluid are dependent on a number of parameters equal to that necessary to describe the problem. Thus, it is well known from dimensional considerations that the drag D of a body can be expressed as the product of a pressure, say the dynamic P n o pressure q = -~ U , a characteristic area S, and some function of Mach number M , Reynolds number R, the geometry of the body, and gas proper- ties such as the ratio of the specific heats 7, etc. This can be written symbolically as follows: D = q_Sf(M , R, geometry, gas properties) (l) (Note that S does not necessarily refer to the plan-form area, but, rather, to any combination of geometric lengths having the dimensions of area.) The geometry of a body can be described by a number of dimension- less parameters, the number of which depends on the complexity of the shape. If these considerations are applied to the present family of rec- tangular wings having NACA 63AXXX profiles, the description is particularly simple since the wing plan form is determined by specifying the aspect ratio A, the profile by the thickness ratio t = t/c and the camber ratio H = h/c, and the inclination of the wing by the angle of attack a. In this way, equation (l) can be rewritten in the following more explicit form: -g = f(M , R, A, t, H, a, gas properties) (2) If, as in nearly all problems of aerodynamic practice, all measurements are made in the same fluid, air, the gas properties can be represented by merely a set of constants and therefore disappear as parameters in equa- tion (2) leaving only 2r = f(M , R, A, r, H, a) (3) qS If attention is confined to the drag at zero lift D , as is the case with the transonic area rule, further simplification results because a is then a dependent parameter and may be eliminated even though it is not a constant for variously cambered wings. Then 2g = f(M , R, A, x, H) (1+) qS The proof of the foregoing step is as follows: The relationship for lift analogous to equation (3) is C L =jg = f(Mo, R, A, t, H, a) (5) CONFIDENTIAL CONFIDENTIAL NACA RM A5^F28 where f is, of course, a different function of the indicated variables. This can be solved for a a = f(Mo, R, A, t, H, C L ) (6) and substituted in its place in equation (3) to produce ^ = f(M Q , R, A, T, H, C L ) (7) In this form, it is clear that the restriction to zero lift fixes the value of the last parameter, thereby leaving D /qS dependent only on the remaining parameters as in equation (k) . Although dimensional considerations such as the foregoing display the parameters or factors that influence the drag, they do not provide any information on the nature of the functional relations involved. Thus, it may be that some of the parameters are more important than others, or that a parameter is of importance over a certain range of values but of negli- gible importance in another range, or that a parameter is important or not depending on the value of another parameter, etc. Examples are that D /qS depends on Mq at transonic and supersonic speeds but not at sub- sonic speeds, D /qS depends on A for small A but not for large A (if S refers to the plan-form area), etc. Since there exists a great body of literature pertinent to the nature of these dependencies which is probably quite familiar to the readers, no further amplification appears necessary here. TRANSONIC AREA RULE Theoretical Considerations The customary verbal statement of the transonic area rule is given in the INTRODUCTION. It should be noted that the term "drag rise" does not refer to the increase in drag force D , but to the increase in the ratio D /q between a subcritical reference Mach number M re f and a transonic Mach number Mq. If Sc(x) represents the streamwise or axial distribution of cross-section area, the transonic area rule states that the variations A(D /q) with Mach number M Q are the same for all low- aspect-ratio wing-body combinations having the same Sc(x). Note that the rule in this form relates a very wide class of wings, bodies, etc., but that all related configurations must have the same longitudinal length c. This restriction can be removed easily by recasting the state- ment in dimensionless form using c as the characteristic length. Thus, we can say that the variation of A(D/qc 2 ) with Mach number Mq is the same for all low-aspect-ratio wing-body combinations having the same CONFIDENTIAL NACA RM A5UF28 CONFIDENTIAL S c (x/c)/c 2 . From a slightly different point of view, the area rule states that, near the speed of sound, A(D/qc 2 ) is a function of Mq and S c (x/c)/c 2 , or r 'D \ ^D Q \ (V \ [ s cU/c) (8) The axial distribution of cross-section area S c (x/c) can be written in the form S c (x/c) = S m s(x/c) (9) where S^ is the maximum cross-section area and s(x/c) is an area- distribution function. Application of the transonic area rule to the drag results of the present family of rectangular wings is facilitated partic- ularly by the fact that s(x/c) is the same for all wings. As a result, the area distribution S c (x/c) of any wing in this family is specified by stating the value of the maximum cross-section area S m , so that A(D /qc 2 ) = f(M , S m /c 2 ) (lO) Since the chord c and thickness t of each wing are constant across the span b, the maximum cross-section area is equal to the product of the span and thickness Sm = bt (11) and Sm b t , , q§ = __ = AT (12) c c c whence equation (10) can be rewritten for this family of airfoils as The remarkable simplicity of these statements is emphasized by com- parison with the functional relation revealed by dimensional consideration alone. Thus, if the unspecified area S is replaced by c 2 , equation (U) yields A (jj*) = A (j^>) = f ( M o, R, A, t, H) - f(M ref , Rref, A, t, H) (lk) where R r ef refers to the value of Reynolds number associated with the subcritical reference Mach number M re f, and the symbol f refers in each case to the appropriate function of the indicated variables. In the cus- tomary discussion of drag-rise data, M re f and R re f are constants and no CONFIDENTIAL 8 CONFIDENTIAL NACA RM A5I1F28 longer appear as parameters in equation (lU). Further simplification occurs in most cases because the wind-tunnel or flight test technique determines a specific relation between the Mach and Reynolds numbers. As a result, either Mo or R can be removed as parameters since the value of either is determined by that of the other. Since the present problem is more closely connected with effects of compressibility than of vis- cosity, it is appropriate to retain M as the significant parameter. In this way, equation (lh) reduces to A (qcV = f(M °' A ' T ' H) (15) Comparison of equations (13) and (15) highlights the fact that the area rule affirms, in dimensionless terms, that the drag-rise parameter A(D /qc 2 ) for the present family of wings depends on Mach number and the product of aspect ratio and thickness ratio At (or the maximum cross- section area parameter S m /c 2 ) but is independent of camber ratio H, and aspect ratio A or thickness ratio t taken separately. Comparison with Experiment Application to wings having identical area distributions . - The applicability of the transonic area rule to the present family of rectangu- lar wings can be examined in several ways. Perhaps the most obvious way is to actually compare the variation of A(D D /qc 2 ) with Mq for two or more wings having identical area distributions. For the present family of wings, this means comparing wings having the same S™/c 2 or At. The tran- sonic area rule predicts that the variation of A(D /qc 2 ) with M should be the same for all such wings. An example of such a comparison is shown in sketches (a) and (b) . The experimental data are from references 2 and 3- Both wings have ./Or .6 .7 .8 Sketch (b) CONFIDENTIAL NACA RM A5UF28 CONFIDENTIAL symmetrical sections (H = 0) and At = 0.l6 t>ut one has an aspect ratio of 2 and thickness ratio of 0.08; whereas, the other has an aspect ratio of k and thickness ratio of 0.0U. The first sketch shows the total drag and the second the drag rise determined by subtracting the value of D /qc 2 at Mq = 0.6. Although the two curves in sketch (b) are not identical as predicted by the area rule, they are closely related. Innumerable reasons could be advanced to explain the differences between the two curves; per- haps there are viscous effects which may significantly affect the drag rise, perhaps the measurements are not sufficiently accurate or the flow field sufficiently uniform, or perhaps the aspect ratio or thickness ratio is too large, etc. In any case, this particular comparison would probably be scored in favor of the transonic area rule . Another comparison, this time among three wings of aspect ratio 2, thickness ratio 0.06, but different amounts of camber is shown in sketches (c) and (d) ficient C z . designation. The amount of camber is specified by the ideal lift coef- in accordance with the NACA scheme for airfoil section Again, the basic data are presented the drag rise in the first sketch and in the second. In this case, however, the variation of A(D /qc 2 ) with Mq is definitely not the same for the three wings and the use of the area rule could lead to serious error. At a Mach number of unity, where the area rule is sup- posed to be most accurate, the drag rise of the wing with greatest camber is nearly twice that of the uncambered wing. Wings having similar area dis- tribution .- Although a certain num- ber of additional comparisons of the type described in the preceding sec- tion can be made using the data of references 2 through 5> the number is definitely limited because the test program was not designed to preserve a single value for At for all wings. It is furthermore not practical to carry out extensive programs of such a type because it necessitates the testing of very thin wings of high aspect ratio and thick wings of low aspect ratio. As mentioned in the derivation of equation (l3), however, all members of the present family of wings have similar distributions of CONFIDENTIAL Sketch (c) .08 .06 .04 .02 A =2 , r = .06 °$ o-gfe^=^ t r ^\ a-oo 2 —$°*ko .8 .9 1.0 I.I Sketch (d) 1.2 10 CONFIDENTIAL NACA RM A51+F28 cross-section area (a single area-distribution function s(x/c)), and the transonic area rule can be extended to include such cases by introducing Sm/c 2 or At as a second parameter. In exchange for being able to cor- relate the drag rise of bodies having not only identical, but also similar area distributions, we incur the complications of a dependence on two parameters rather than only one. Thus, whereas the curves representing the variation of A(Do/qc 2 ) with M Q for all wings having identical Sjj/c 2 or At should coincide to form a single line, those for a family of wings having similar area distributions should form a family of lines. Simplicity can be regained, however, by restricting attention to a single Mach number and ascertaining the variation of A(D /qc 2 ) with S m /c 2 or At. Since the customary statement of the area rule restricts attention to near-sonic speeds, the most appropriate single Mach number to select for such a comparison is unity. Sketch (e) shows the variation of A(Do/qc 2 ) with Sm/c 2 for all the uncambered wings of references 2 and 3« It can be seen that all these results fall near to a single curved line for wings of all aspect ratios up to 3 hut that those for wings of aspect ratios k and 6 depart from this line in a systematic manner. Sm/c sv .25 The same results are replotted in sketch (f) versus the square of rather than the first power. It can be seen that the curved line .44 .40 .36 .32 28 $.20 < .16 .12 .08 .04 P h to-l -0 / D A 1 / / / ^-A't / 1/ / p f -A*4 I # 1 Sketch (e) ./ .2 (S m /C)' . (At)' Sketch (f) CONFIDENTIAL NACA RM A5^F28 CONFIDENTIAL 11 of sketch (e) for wings having aspect ratios less than 3 is now a straight line, indicating that the sonic drag rise for these wings is proportional to the square of the maximum cross-section area. Such a dependence is consistent with the formulas of linearized compressible flow theory for the wave drag of nonlifting slender wings and bodies at supersonic speeds. These plots can be interpreted as showing that the sonic drag rise of the uncambered members of the present family of wings depends on the cross-section area in accordance with the area rule, provided the aspect ratio is about 3 or less. On the other hand, the results for wings of higher aspect ratio can only be interpreted as indicating that some parame- ter other than the cross-section area must be involved. Inasmuch as the geometry of wings of the present family is described completely by the two parameters, aspect ratio and thickness ratio, it is clear that these parameters must assume importance in some form other than their simple product At for wings of larger aspect ratio. One could seek the new relation empirically, but the transonic similarity rules provide a theo- retical basis for proceeding. Some properties of transonic similarity rules are reviewed in the following section, although the reader is refer- red to the original references for further details. TRANSONIC SIMILARITY RULES Statement of Rule Transonic similarity rules are derived from the nonlinear equation of inviscid flow theory and are known for thin wings (e.g., ref . 6) and slender bodies of revolution (ref. 7)> but not for wing-body combinations with pointed noses. In contrast to the transonic area rule which relates the zero-lift drag rise of families of bodies having identical or similar axial distributions of cross-section area, the transonic similarity rules relate the aerodynamic properties of much more highly restricted families of bodies. Even for wings alone, the restrictions imposed on the members of a single family are much more severe than for the area rule, since all members of a single family must have affinely related plan forms, affinely related thickness distributions, and affinely related camber distributions. To be more explicit: if the plan form is given by ^2 = Y(x / c) < 16 ) CONFIDENTIAL 12 CONFIDENTIAL NACA RM A5^F28 y-frW as indicated in sketch (g), it is required that Y(x/c) be a single function for all wings of a given family. Furthermore, if the ordinates of the wing surface are given by 2 = * Li c c ~2 * 1 h + t 1 I a x t/c c (IT) Sketch (g) where the plus sign is associated with the ordinates of the upper surface Zu and the minus sign with those of the lower surface Z z , it is required that g 1 (x/c, y/b) and g 2 (x/c, y/b) be single functions for all wings of a given family. The first of these restrictions requires that related plan forms be obtainable one from another by a differential lengthening of lateral and longitudinal dimensions. Thus rectangular plan forms consti- tute one family, triangular plan forms with straight trailing edges another, etc. Examples of related plan forms are illustrated in sketch (h) . The relationship expressed in equation (l7) requires that the thickness distribution and camber □ -> Sketch (h) /\ variations must be the same for all wings of a particular family. The magnitudes of the maximum thickness A ratio t/c and camber thickness ratio h/c, as well as the angle of attack a, may be different for vari- ous members of a single family. From the present point of view, one of the most significant proper- ties of the family of wings described above is that all members have the same dimensionless area-distribution function s(x/c). This can be shown as follows: +Y(x/c) +Y S(x/c) = J (Z u - Z z )dy = I J -Y(x/c) J -Y tfc (§,*)* = bt +Y/b r -Y/b f'f = btg(x/c) = SmsU/c) (18) CONFIDENTIAL NACA RM A5^F28 CONFIDENTIAL 13 where g represents the integral of the preceding expression and is pro- portional to s(x/c), since Sm is proportional to bt for a family of wings. Because of this fact, a family of wings suitable for correlation by the transonic similarity rules is also part of a family of bodies suit- able for correlation by the area rule. For wings of such a family, the similarity rules of inviscid, small- disturbance, transonic flow theory provide that the wave drag LV be given by T 5/3 = ^P [Mo % + l)]l/3 f (t > K > £ ' S) U9) where T t/c Mo 2 & [( 7+ l)M 2 T ]2/3 A [(7+l)MoS-] 1/& A h h/t a a/r and S„ is the plan-form area of the wing. Equation (19) is not in the proper form for comparison with the area rule since the latter is con- cerned with conditions at zero lift and, hence, imposes an indirect requirement on a. We can proceed toward the desired form, however, by introducing the transonic similarity rule for the lift coefficient C L = L/qSp ° L = [M^l)]^ f( *o' l > *> °> (20) and defining the reduced lift coefficient [M g (y + l)] 1/3 T 2 / 3 C L - L "° " + i" C L = f (* , A, h, a) (21) Now since C^ is a function of the same four variables that appear in equation (l9)> cc can be replaced with Cl i- n the latter equation, whence °w = *Sp [^g^Jl/a f (t , A, h, C L ) (22) C0NFIDENTIA1 Ik CONFIDENTIAL NACA RM A5UF28 The condition of zero lift eliminates the last parameter, leaving only D .5/3 Ov = qS p i^Li)]!'* f{ to> l > 5) qS -5/3 r Mq 2 P [Mo 2 (7+Dl 1/3 [[(7+-)Mo 2 t] 2/ [(7 + l)Mo 2 T ] 1/3 A, h/t! (23) Application to Rectangular Wings and Comparison with Transonic Area Rule The restrictions introduced in the derivation of the transonic simi- larity rule for zero-lift wave drag are such that they permit the direct application of equation (23) to the present family of rectangular wings having NACA 63AXXX profiles. It is natural to compare this functional relation with the corresponding relations of equations (lO) and (l3) given by the transonic area rule. A G?) = f(M °' Sm/c2) = f(M °' at) (2k) At first glance, the two sets of relationships appear to bear only slight resemblance. It can be seen upon closer examination, however, that some of the apparent differences are superficial and of little or no signifi- cance. For instance, equation (23) is concerned with wave drag Dq^. whereas equation (2k) is concerned with drag rise. It is evident, how- ever, that the two rules are actually concerned with the same quantity, since the drag rise can be considered to be an approximation for the wave drag under the assumption that the friction drag coefficient Is independent of Mach number. Another point of apparent lack of resemblance is that equation (23) does not show an explicit relationship between wave drag and maximum cross-section area, since the latter is not a function of SpT 5 ^ 3 nor of | , A, or h taken separately. This does not mean that the transonic similarity rule is incompatible with the area rule because there are sev- eral possibilities for making such a dependence visible. Two permissible procedures are to multiply the right side of equation (23) by either i Q or A. The first procedure is of no help for the discussion of conditions at M = 1, since the quantity M 2 - 1 appears in two places and an inde- terminate form ensues. The second is perfectly acceptable, however, and produces the following relationships: CONFIDENTIAL NACA RM A5^F28 CONFIDENTIAL 15 Dow = qSpAr 2 f(| Q , A, h) qC \^y) f( ^0' K > 5) ( 2 5) Since the maximum cross-section area Sjq of the present family of rec- tangular wings is equal to the product bt, equation (25) can be rewritten as follows: D ov = ^ c2 (^i) 2 f U , A, h) c 2 / 1[(7+i)Mo 2 t] 2/3 ' [(7+1)Mo 2 t] 1/3 A, h/tl (26) Equation (2*0 can likewise be rewritten by multiplying f by the square of Sm/c 2 , A &) " (¥J f(M °' "^ ' (?J »(*. *T) (27) This appears to be the closest that the two rules can be brought together without introducing additional simplifications. Both rules are now concerned with essentially the same quantity, D^/qc 2 and A(D /qc 2 ) . Each rule states that this quantity is proportional to the square of Sm/ c2 times some unknown function of certain specified parameters. The two rules disagree completely, however, as to the nature of the parameters. The transonic similarity rule specifies three parameters i Q , A, and h, whereas the area rule specifies only two, Mq and At. Since neither set can be transformed into the other, it is apparent that the only way in which both rules can be universally correct is for the function f to be actually a constant and not to depend on the value of any of the five parameters. It is obvious that such is not the case, however, since it requires, for instance, that the wave drag be independent of Mach number. There is a way out of this apparent impasse, however, if the range of validity of one or both of the rules is restricted sufficiently that f is independent of the remaining variables. Comparison of the two rules shows that the drag depends in both cases on the Mach number of the stream and the geometry of the wing and, in the case of the transonic similarity rule, on the ratio of the specific heats 7. Since nearly all problems of aerodynamic interest are concerned solely with air, however, y is a con- stant and need not be retained as a parameter. The resulting simplified equation can be written in full in either of the following forms: CONFIDENTIAL 16 CONFIDENTIAL NACA RM A5^F28 r »/r 2 £ss . ra»Y f r,»i^. (vm^ *, vt 2 V 2 2 T \2/3 ; = (AT) 2 f (Mq 2 T) Mo 2 - 1 ,„ , ,1/3 / („ 2 T )2/3^ (Mo T) A, h/t (28) With regard to further simplifications, it should be noted that the derivation of the transonic similarity rule requires that the thickness and camber be small with respect to the chord, but does not restrict the Mach number or aspect ratio. The transonic area rule, on the other hand, restricts itself to Mach numbers near the speed of sound and to wings having low aspect ratios. The functional relation given in equation (27) representing the area rule simplifies if attention is fixed on any given Mach number, since Mq appears as an isolated parameter, thus Do \ /Sm On the other hand, the only Mach number at which the corresponding rela- tion given by the transonic similarity rule simplifies, irrespective of the thickness ratio, is Mq = 1. The resulting expression is 22g = (J|> f(A T 1/3 , h/t) (30) This expression simplifies further for wings having symmetrical sections, since h/t is zero. For such cases, the transonic similarity rule for wave drag at Mq = 1 reduces to ASmY f ( AT i/3) ( 31) qc 2 WJ The foregoing analysis has developed certain equations relating to the drag rise or wave drag of a family of wings having relatively simple geometry. Despite the great restrictions imposed by the selection of such a family of bodies, the analysis has disclosed a number of significant points which complement those discussed in connection with the experimental data presented in sketches (a) through (f). First of all, it has been shown that the only way in which the wave drag at Mq = 1 can depend on the axial distribution of cross-section area (defined in the present fam- ily of wings by the value of the maximum cross-section area S m ) and still CONFIDENTIAL NACA RM A5^F28 CONFIDENTIAL 17 be compatible with the transonic similarity rule for drag is for Dov/qc 2 to be proportional to the square of S m /c 2 . This is precisely the depend- ence disclosed by the experimental data for low-aspect-ratio wings having symmetrical sections and shown in sketch (f). The transonic similarity rule states that the wave drag may depend on the camber, whereas the area rule states that it does not. The experimental data in sketch (d) show that camber has a significant effect on drag. It is evident, both from a priori considerations and from the experi- mental results shown in sketch (f), that some change must occur in the relation between wave drag and maximum cross-section area as the aspect ratio becomes very large. The only possibility permitted by the transonic area rule is that A(D /qc 2 ) varies with S m /c 2 in some other manner than as the square. It can be seen from sketch (e), however, that the data for wings of the present family cannot be correlated on this basis if the aspect ratio is greater than about 3- The transonic similarity rule for the zero-lift wave drag of uncambered wings at Mq = 1 provides a differ- ent dependence by stating that Dc^/qc 2 is equal to the square of Sm/c 2 times some function of At 1/3 . This statement is compatible with the area rule if Dov/qc 2 is independent of At 1/3 for small values of the latter. If wave-drag data for wings of larger aspect ratio can be cor- related successfully by considering the value of this quantity, a theoret- ical basis for a limit to the range of applicability of the transonic area rule has been found. Comparison with Experiment The foregoing discussion has served to focus attention on the fact that the parameter At /3 , already familiar from prior papers on transonic flow (e.g., ref. 3)> may be of importance in defining the limit of appli- cability of the transonic area rule as applied to a family of af finely related wings. The functional relation of equation (31) suggests that if D^c^qSm 2 is plotted as a function of At 1/3 , all the data for the uncambered wings should fall on a single curve. Moreover, the values of D(v. c 2 /qS m 2 should be independent of At 1 /3 over whatever range the area rule applies. Although this method of plotting is conceptually simple, it imposes severe requirements on the accuracy of the experimental deter- mination of wave drag because any errors are magnified as a result of dividing by the square of the maximum cross -section area. For results such as the present in which the wave drag is not actually measured, but is inferred from measurements of total drag by subtracting an estimated friction drag, greatest difficulties are experienced if the wings are thin and of low aspect ratio so that the wave drag is only a small fraction of the total drag. Consequently, data points for which the wave drag is less than half the total drag are omitted in the plot of A(D c 2 /qS m 2 ) versus CONFIDENTIAL 18 CONFIDENTIAL NACA RM A5UF28 3.0 25 20 15 1.0 AT 1 'a shown in 8ketch (i). Even with this precaution, the data evidences considerable scatter for thin, low- aspect-ratio wings. (The symbols refer to the same wings as in sketch (e).) The principal points of interest in this plot are threefold. First, except for the smallest values of At 3 where the scatter is too large to provide any positive conclusions, the points determine essentially a single line, indicating that the sonic wave -drag characteristics of the uncam- bered wings of the present family can be correlated successfully by the tran- sonic similarity rule (this has been shown previously by McDevitt, ref. 3); second, A(D c 2 /qSm 2 ) is, at best, independent of At 1/3 for values of the latter up to about unity, indi- cating that the drag in this range varied in accordance with the area rule as well; third, A(D c 2 /qS m 2 ) varies appreciably with At 1/3 at larger values of the latter, indicating that the range of applicability of the area rule was exceeded. 1 1 1 /-A*l <b y, a-K* — A=4 * /'? L A=. p M„ = l r A=6 h--0 1.2 1.6 At* 20 24 28 Sketch (i) Before leaving the discussion of the correlation of drag data at Mq = 1 by use of the transonic similarity rule, it is worthwhile to call attention to another method of plotting which illustrates the same factors although in a slightly different manner. This procedure, which is described more completely in references 3> 5> and 6, is based on equa- tion (23) rather than equation (26) and consists, for the zero-lift drag of a family of uncambered wings in an air stream with Mq = 1, of plotting the variation of Cj)~ / t5/3 w ith At 1/3 where Cj)~ re P resen "t s ^W/qSp. Sketch (j) shows the data of sketch (i) replotted in this manner. (Once again, the lack of wave-drag data requires the substitution of drag-rise information and the introduction of the drag-rise coefficient ^D n > defined as equal to A(D /qS p ).) The curve formed by the data points on this type of plot is perhaps somewhat simpler than that of sketch (i) since it is asymptotic to straight 3.5 3.0 25 2.0 1.5 1.0 .5 ^ h*r*- -4 -A = 6 / / 1/ ! - 4=2 Mo'l /i--0 / SI --/ / lines at both large and small At 1/3 At small At 1/3 1.2 1.6 Ari 2.0 24 2.8 the points define a Sketch (j) straight line passing through the origin. Such a line is in accordance CONFIDENTIAL NACA RM A51+F28 CONFIDENTIAL 19 with the area rule. At values of At 1 / 3 larger than about unity, how- ever, the line determined by the data points departs from this initial trend and turns toward the horizontal. This trend is contrary to the area rule but consistent with the fact that the results for wings of high aspect ratio must tend toward those for wings of infinite aspect ratio. It cannot be emphasized too much that the critical value of unity for At 1/3 is determined solely on the basis of data for a very special family of rectangular wings having symmetrical profiles. Other families of wings would be represented by different curves on such a plot. It is interesting to consider for a moment the nature of this limit and to compare it with the verbal restriction of the area rule to low- aspect-ratio wings. Although the one requires that AtI^ 3 t> e small, and the other that A be small, these two statements are in better agreement, insofar as engineering applications are concerned, than might appear at first glance. This agreement results from the fact that other consider- tions, such as designing for structural strength or for the avoidance of excessive drag, tend to preserve a rather narrow range for the values of t likely to be met in practice. The effectiveness of t as a parameter is further diminished by the fact that only its cube root is involved. Consequently, the restriction to small At*' 3 represents a limitation primarily on the aspect ratio and only secondarily on the thickness ratio. A=I, T = 04, Ar*=342 A =2, r'./O, At*= 928 It is also interesting to exam- ine the shapes of equivalent bodies of revolution having the same axial dis- tribution of cross-section area as wings lying on either side of the limit. Accordingly, sketch (k) has been prepared showing the shapes of bodies of revolution equivalent to one of the wings having At 1 / 3 much less than unity, to the two wings having At 1 /3 nearest unity, and to the wing having the largest At 1/3 of any tested. It can be seen that the equivalent bodies are blunt and stubby rather than pointed and slender. Thus, although it has been shown that the drag-rise characteristics of the var- ious members of the present family of uncambered rectangular wings are related to one another in the manner predicted by the transonic area rule, provided At 1 3 is less than about unity, it would appear that the drag- rise characteristics of the equivalent bodies of revolution might be con- siderably different. It is shown in sketch (d) that camber has an effect on the zero-lift drag rise. The transonic similarity rule suggests that the ratio of maxi- mum camber to maximum thickness h/t might be an appropriate parameter to A=3,r= 04, Ar*= 1.026 A=6,r- 10, A T *=2784 Sketch (k) CONFIDENTIAL 20 CONFIDENTIAL NACA RM A5^F28 use in addition to At 1/3 to correlate the sonic zero-lift wave-drag characteristics of a family of cambered wings. Accordingly, sketch (z) has been prepared showing the influence of h/t for wings of various At 1 '3. For reasons of simplicity, only three values for h/t, namely, 0, 0.222, and O.Mj-U, are included on this plot. The important effects of camber are readily evident from this graph. It is appar- ent that the area rule is not applica- ble to wings having different camber- thickness ratios, even if the values of At 3 are sufficiently small to permit successful correlation of the wave drag of uncambered wings. On the other hand, it is permissible in the theory and seems to be indicated by the experimen- tal data that the area rule is applica- ble to families of cambered wings pro- viding h/t is maintained constant. This result is recognized in sketch (Z) by the fact that A(C D /t 5 /3) for the wings of constant h/t is approximately proportional to At 1/3 . 4.0 3.5 30 25 \T 20 4. < 1.5 10 .5 h/t -- 444 (63 A 406} c / V 1 ^Q / □ / // I V 1 "\_ h/t =22 2 {63 A 206 J / ) -h/t-- c / ) / M --/ 12 1.6 Ar* 2.0 24 23 Sketch ( Z) SUMMARY OF RESULTS The range of applicability of the transonic area rule has been inves- tigated by comparison with the appropriate similarity rule of transonic flow theory and with available experimental data for a large family of rectangular wings having NACA 63AXXX profiles. These wings are of af finely related geometry and are hence immediately amenable to analysis by the transonic similarity rules. On the other hand, the axial distributions of cross-section area are not identical, in most cases, but merely similar. (The ratio of the local cross section to the maximum cross section is a given function.) It is shown, however, how the transonic area rule can also be used to correlate the sonic drag-rise data for such a family of wings. It is found that the sonic zero-lift drag-rise data for the present family of wings can be successfully correlated on the basis of the area rule, provided the wing profiles are symmetrical and the product of the aspect ratio and the cube root of the thickness ratio is less than about unity. Within this range, the sonic drag rise varied as the square of the maximum cross-section area, all wings having equal chords. It is demon- strated that this is the only dependence of drag on maximum cross-section area for a family of wings like the present that is compatible with both the area rule and the transonic similarity rule. CONFIDENTIAL NACA RM A5UF28 CONFIDENTIAL 21 It was found that the addition of camber greatly increased the sonic drag rise and that the application of the transonic area rule to a family of wings, some of which are cambered and others not, could lead to serious error. On the other hand, it is indicated by the transonic similarity rules and the experimental data that the area rule is applicable to fami- lies of cambered wings, provided the camber distribution, as well as the area distribution, are similar and that the ratio of the maximum ordinates of the camber and thickness distribution is maintained constant. Ames Aeronautical Laboratory National Advisory Committee for Aeronautics Moffett Field, Calif., June 28, 195 1 *- 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. Nelson, Warren H., and McDevitt, John B.: The Transonic Characteris- tics of 17 Rectangular, Symmetrical Wing Models of Varying Aspect Ratio and Thickness. NACA RM A51A12, 1951. 3. McDevitt, John B.: A Correlation by Means of the Transonic Similarity Rules of the Experimentally Determined Characteristics of 22 Rec- tangular Wings of Symmetrical Profile. NACA RM A51L17b, 1952. h. Nelson, Warren H., and Krumm, Walter J.: The Transonic Characteris- tics of 38 Cambered Rectangular Wings of Varying Aspect Ratio and Thickness as Determined by the Transonic-Bump Technique. NACA RM A52D11, 1952. 5. McDevitt, John B.: A Correlation by Means of Transonic Similarity Rules of the Experimentally Determined Characteristics of 18 Cam- bered Wings of Rectangular Plan Form. NACA RM A53G31, 1953. 6. Spreiter, John R.: On the Application of Transonic Similarity Rules to Wings of Finite Span. NACA Rep. 1153, 1953- 7. Oswatitsch, K., and Berndt, S. B.: Aerodynamic Similarity at Axi sym- metric Transonic Flow Around Slender Bodies. Kungl. Tekniska Hogskolan, Stockholm. Institutionen for Flygteknik. Tech. Note 15, 1950. CONFIDENTIAL NACA-Langley - 8-23-54 - 325 •J J J e -*r O <-> « >i - -j 01 3 01 3 a feiS H X s « a < a. u OT z co i— i £3 < < CM ^ -1 J e Z o u o u <->K ~ - >, .-^ _ -co t-. m u - 5 0D O 00 01 ' O 3 01 3 Q, jh c •* c IO o •-a <: S « ■jl < QJ ! ) < w z «• a W SS H 3 u. O >- 01 H <] ►J oi r; oi y S -J ■v co O < -a <C g < K 05-3 W < s x < rt z z z o H z w 2 a Si b. a> z ago OS m < O — CO 153-2 o •S ° — . CJ rt n a -a 9- c 01 Q, 01 OD - u X «g 01 i to . < W U -I < P z K" < • K CM < U -3" ~ m Z os w — Z » <: a K 3 H < c -Q rt yi 3 « < ^ -ay a - z m *j OD ™ 3 3 s S > 5 <u oj en s-2 a CO ■" s-, 2 3 -3 o u a ■" en is 5-3 u B z o u J J.J 3 t! ■ * b •d y- 01 3 a £ ir> S -2 o <• a s u >. ^ . 10 ■- « O"* 1 o u m a CO -J- IS OJ M ^5 oj y " j S Oh %< o u. 3 > z < -c < g < K K «S ^ c X 5-Sc- < « z z z o H z w oi 9 i" „z MgO • T OS in c <c Is . < w u a z K CM < U -a-' O — CO a ■- <" £ U — ' o >- -— 01 cs rt ^ b£ o c O) 01 3 a < - Z 31 u n] OO •z: co 5< o .S a ^ at £ oo ■« an S !5 a jz o — t - "?3 Z as 'to o Z " aj X OT t- CM »-l (x, 3 -^ Z o u >< O O rt SU Uk - CO t- CO 01 > 00 O 00 01 o c ai c a — ^ -3 l™ « ^s < s 01 K ~ < 01 ( ) < OT Z j a H Z z o u < CONFIDENTIAL UNIVERSITY OF FLORIDA UNIVERSITY OF FLORIDA DOCUMENTS DEPARTMENT 1 20 MARSTON SCIENCE LIBRARY i. BOX 117011 SVILLE, FL 32611-7011 USA CONFIDENTIAL