UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Multi-element thin airfoil theory Watt, George Donald 1984

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata

Download

Media
831-UBC_1985_A1 W37_7.pdf [ 8.18MB ]
Metadata
JSON: 831-1.0096761.json
JSON-LD: 831-1.0096761-ld.json
RDF/XML (Pretty): 831-1.0096761-rdf.xml
RDF/JSON: 831-1.0096761-rdf.json
Turtle: 831-1.0096761-turtle.txt
N-Triples: 831-1.0096761-rdf-ntriples.txt
Original Record: 831-1.0096761-source.json
Full Text
831-1.0096761-fulltext.txt
Citation
831-1.0096761.ris

Full Text

Multi-Element Thin Airfoil Theory  by George Donald Watt B.A.Sc, The University of British Columbia, 1977  A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of D O C T O R O F PHILOSOPHY in  The Faculty of Graduate Studies Department of Mechanical Engineering  We accept this thesis as conforming to the required standard  The University of British Columbia December 1984 ©  George Donald Watt, 1984  In p r e s e n t i n g  t h i s t h e s i s i n p a r t i a l f u l f i l m e n t of  requirements f o r an advanced degree a t the  the  University  of B r i t i s h Columbia, I agree t h a t the L i b r a r y s h a l l make it  f r e e l y a v a i l a b l e f o r reference  and  study.  I further  agree t h a t p e r m i s s i o n f o r e x t e n s i v e copying of t h i s t h e s i s f o r s c h o l a r l y purposes may  be  department or by h i s or her  granted by the head o f  representatives.  my  It i s  understood t h a t copying or p u b l i c a t i o n o f t h i s t h e s i s f o r f i n a n c i a l gain  s h a l l not be  allowed without my  permission.  Department of  Mechanical  Engineering  The U n i v e r s i t y of B r i t i s h Columbia 1956 Main Mall Vancouver, Canada V6T  1Y3  Date  DE-6  (.3/81)  December 10,  1984  written  ii  Abstract A linearized, two-dimensional, potential flow analysis of multi-element airfoil configurations is attempted. Arguments are presented that suggest the accuracy of the linearized theory should be as good as or better than that of the well-known one-element thin airfoil theory — despite the large mean line curvature common to multi-element configurations. The results obtained in this thesis tend to support this expectation. Good success is achieved in developing a general two-element tandem airfoil linearized theory.  The effects of incidence, leading or trailing edge flap deflection, and  camber on the overall and localized lifts and moments are summarized in integrals which are the tandem airfoil versions of the historical Munk integrals. Analytical solutions for the forces on tandem N A C A airfoils are then obtained from these integrals. The expressions for the forces reduce, when the two airfoil elements come together, to the familiar one-element thin airfoil theory formulas for an airfoil with a simple flap. All the forces are calculated on a small hand-held computer and the results are compared with exact potential flow theory. The overall force results for the tandem airfoil incidence problem are particularly simple and, in fact, this simplicity allows the solution of the incidence problem for an arbitrary number of in-line airfoils. Although thickness has no effect on tandem thin airfoil theory forces, it does affect pressures. The general thickness distribution analysis for tandem airfoils, an exceedingly simple analysis involving only elementary functions, is presented. A similar type of analysis enables the design of tandem airfoil camber lines. Examples are given. Modelling the effects of overlap of a staggered two-element airfoil configuration is a much more complicated analysis than that of tandem airfoil elements. Nevertheless, substantial success is achieved. This includes the development of the equivalent Munk integrals that give the overall and localized forces on airfoil elements of arbitrary chord length, arbitrarily positioned relative to one another (providing the chord lines are approximately parallel to each other and to the flow at infinity — the linearization requirement). Although analytical solutions are obtained for the overall forces of the incidence and flap deflection problems, they are in terms of parameters which can only be obtained through a trial and error solution of the staggered parallel slit conformal mapping.  ui  Contents Abstract Tables  v  Figures  vi  Symbols  vii  Acknowledgement  x  Classical Thin Airfoil Theory  1 1.1 1.2 1.3 1.4  Introduction The potential flow - real flow correlation The linearizing procedure One-element thin airfoil theory functions and forces . . .  1 1 5 12 18  Tandem Symmetrical Airfoils  2 2.1 2.2 2.3 2.4  The mappings The flow functions The forces N in-line uncambered airfoils  24 24 29 32 37  General Tandem Thin Airfoil Theory  3 3.1 3.2 3.3 3.4 3.5 3.6  The thickness distribution analysis The camber line design problem Camber line analysis: general functions Camber line analysis: special functions Thickness distribution design and other problems . . . . Munk's integrals extended  39 41 47 52 60 64 65  Tandem Airfoil Forces an Easy Way  4 4.1 4.2 4.3 4.4  The overall forces The localized forces Segmented boundaries: the segment selectors Tandem N A C A airfoils, an example  66 66 72 76 84  Staggered Two-Element Thin Airfoil Theory  5 5.1 5.2 5.3 6  ii  The staggered parallel slit mapping The overall forces The localized forces  100 100 Ill 119  Conclusions and Recommendations  123  References  126  Contents  iv  APPENDICES A  Exact-Numerical Potential Flow Computations A.l A.2 A. 3  B  C  . 127 129 134  .  141  General 141 Elliptic function series representations in powers of q* . . 145 Elliptic function series representations in powers of q 147 16  The Selector Function Expansion Coefficients Cl C. 2 C.3  D  The flow field formulation The forces Computational accuracy Elliptic Functions  B. l B.2 B. 3  127  £{z) and F{z) A(u) and ${OJ)  and T(w)  Computer Program: The Tandem Airfoil Forces  153 153 154 . 156 159  V  Tables 1 2 3  Potential flow predictions of the forces on the tandem airfoils of Figure 8  50  Comparisons of staggered slit mapping parameter values at the tandem and biplane slit extremes  109  Comparisons of the exact and present numerical predictions of the lift on Williams' configuration A slotted flap arrangement  134  vi  Figures 1  Lift characteristics of a NACA 23012 airfoil with a 20% plain flap  2  N arbitrary airfoils  13  3  Joukowski transformation for the one-element problem  19  4  The tandem airfoil mappings  25  5  N in-line uncambered airfoils  37  6  Solution planes for the thickness and camber problems  41  7  Thin airfoil theory predictions of the pressure distributions over tandem 15% thick ellipses  44  Potential flow predictions of tandem airfoils designed by thin airfoil theory to give uniform chordwise load distributions at a = 0° .  51  The elementary mappings for the second camber function  55  8 9  .  .  7  Figures 10 to 17: The aerodynamic characteristics of tandem NACA 23012 airfoils. 10  Lift curve slope characteristics  89  11  Flap deflection effects on lift  90  12  Camber line lift characteristics  91  13  The aerodynamic centers  92  14  Flap deflection effects on moment  93  15  Camber line moment characteristics  94  16  C . vs. a for different 7/  97  17  C . - vs. a for different ij  98  18  The staggered parallel slit mappings  19  Doubletflowin the A-plane corresponding to zero incident uniform flow in the 2-plane  20  t  m  The control point Zj in the middle of the ;  101  t h  segment  102 131  Figures 21 to 25: The accuracy of the exact-numerical calculations. NACA 23012 airfoil elements in a slotted flap configuration at a = 0°. 21  Overall lift and drag  136  22  Localized lift  137  23  Localized drag  138  24  Overall leading edge moment  139  25  Localized leading edge moment  140  vii  Symbols a  p o l y n o m i a l b o u n d a r y c o n d i t i o n coefficients.  n  A  F o u r i e r series c o e f f i c i e n t s f o r t h i c k n e s s flow f u n c t i o n .  A{OJ)  airfoil selector function.  b  p o l y n o m i a l b o u n d a r y c o n d i t i o n coefficients.  n  n  B  F o u r i e r series c o e f f i c i e n t s f o r c a m b e r flow f u n c t i o n .  n  c  c h o r d l e n g t h o f o v e r a l l c o n f i g u r a t i o n : (c^ H  Cj  chord length of airfoil ; .  Ci = -r—z—  o v e r a l l lift coefficient f o r airfoil c o n f i g u r a t i o n .  Ci. = - j — ^ —  l o c a l i z e d l i f t c o e f f i c i e n t f o r a i r f o i l j.  Cm* ~  o v e r a l l m o m e n t coefficient ( a b o u t the o r i g i n ) .  i ^1°  h  c ). N  M Cfju'n — i — ^ 7 tt  C  i' 2o<7 , C -  mac  mj  , A C  C  p  C  t  l o c a l i z e d m o m e n t coefficient f o r airfoil airfoil;).  overall a n d localized a e r o d y n a m i c center p i t c h i n g m o m e n t coefficients; E q s . 4.72 a n d 4.73. t h e s h i f t i n a c o e f f i c i e n t ' s v a l u e a t a = rj = 0 ° d u e t o t h e p r e s e n c e o f a camber line.  m  = n—^p-  pressure coefficient.  Cj  closed contours i n the complex physical plane about the overall configuration  d  F o u r i e r series c o e f f i c i e n t s f o r t h e 2  n  e  n d  camber function.  p o l y n o m i a l b o u n d a r y c o n d i t i o n coefficients.  n  E,  a n d a i r f o i l j.  p o l y n o m i a l b o u n d a r y c o n d i t i o n coefficients.  n  D  ( a b o u t t h e l e a d i n g edge of  E'  c o m p l e t e elliptic integrals of the 2  n d  kind.  £  error i n selector f u n c t i o n expansion coefficients.  /  linearized t a n d e m airfoil physical plane geometry, F i g . 4.  F(z) 7{(JJ)  flow flap  function for solving the t a n d e m airfoil a - p r o b l e m . selector function.  viii  Symbols  constants, Eq. 2.37. h  the gap in a staggered two-element airfoil configuration, Fig. 18.  H  constant, Eq. 5.16.  U*)  selector functions for the localized airfoil force solutions, Eq. 4.68.  k, k'  modulus and complementary modulus of elliptic functions,  K,  complete elliptic integrals of the 1 kind, st  K'  length of segment  Fig. 20.  m  the boundary condition selector.  M  number of airfoil elements in a configuration (Appendix A only).  N  number of airfoil elements in a configuration (except in Appendix A). number of segments in the configuration's boundary and in the boundary of the k  th  airfoil (Appendix A only).  q, q'  theta function parameter, Appendix B.  R  radius to outer annulus circle representing the flap, Fig. 4. flap chord to total chord ratio: (1—f)/{s + l — /). linearized tandem airfoil physical plane geometry, Fig. 4.  3  airfoil segment selector. <• *j  thickness to chord ratios.  i-plane  "tandemized" staggered slits, Chapter 5, Fig. 18.  T(«)  flap segment selector,  w(z) = u  —  tu  w'{z) = u' — iv' Xae  X  >  X  j>  Jac  X  y(x),  jj  y'(x)  z = x + iy  complex velocity. nondimensionalized complex disturbance velocity, section 1.3. aerodynamic centers, Eqs. 4.72 and 4.73. leading and trailing edge x-coordinates of airfoil j. airfoil boundaries, section 1.3. linearized complex physical plane.  angle of attack of airfoil configuration.  P  6  orientation of doublet in the rectangle, Fig. 19. vorticity (Appendix A only). staggered slit mapping parameter, Eq. 5.21. "special* camber flow functions. thin airfoil theory solution planes.  Symbols  ix  rj  flap deflection angle.  6  angular variable in ^-planes.  i>i, t ? , t?3, #4  the four theta functions, Appendix B.  A = fi + iv  the rectangle solution plane, Figs. 4 and 18.  p  fluid density.  2  S , E 3  4  flow function constants.  r  theta function parameter,Appendix B.  <j>  velocity potential.  V>  stream function.  ip  value of stream function on & airfoil element, Appendix A. th  k  u = x + *7  modified rectangle solution plane, Fig. 4.  u' = x' + *V  modified rectangle solution plane, Fig. 4.  fl = <f> + i0  complex potential.  X  Acknowledgement I would like to thank Dr. G.V. Parkinson for the opportunity he has given me here at UBC and for the enthusiasm and support he has shown over the years. A major factor in the successful solution of a complicated problem is the problem's initial formulation: Dr. Parkinson's insightful suggestion to linearize the doubly connected flow about a slotted flap onto tandem slits provided this graduate student with a problem with extremely simple solutions just waiting to be discovered. The research was carried out under a grant from the Natural Sciences and Engineering Research Council of Canada. Additional support was recieved from the University Graduate Fellowship and H.R. MacMillan Family Fellowship funds. My thanks to my wife, Ying, who translated this entire manuscript into English. Also, my thanks to Mr. Feng Chen who donated a large part of his summer holidays towards the TgXing of this thesis. T^jX is a computerized typesetting process (sponsored by the American Mathematical Society) that is particularly suited to mathematically oriented dissertations such as this. My thanks to Dr. Larry Roberts for his work in getting TjjjX on-line in time for the processing of this thesis and for his help in achieving that end.  CHAPTER 1 Classical Thin Airfoil Theory  1.1  INTRODUCTION  In recent years there has been a great deal of interest in short take-off and landing aircraft. Their success is, to a large extent, attributable to the low-speed, high lift airfoil configurations that have been developed over the years. These multi-element configurations often employ multiple slotted flaps and a leading edge slotted slat to achieve large mean line curvature, and thereby high lift, while minimizing boundary layer growth and separation effects which lead to stall. Naturally, the development of these airfoils has been closely followed by the development of theories capable of predicting the aerodynamic forces associated with them. In keeping with what seems to be a general trend in engineering today, these theories are almost all numerical. They exploit to the fullest the power of big, but fast, large memory, mainframe computers. , There are several numerical theories available that can accurately predict the potential (i.e. inviscid)flowsabout multi-element airfoils, and in conjunction with boundary layer theory they can also give good predictions of the real flow lifts and moments. These potentialflowtheories are basically of two types. First, there are the more analytically oriented theories which use numerical methods to conformally map arbitrary airfoil shapes to circles, the flow about which is solved as accurately as required by superposing doublets and vortices on the flow in the manner of Williams (1971). The forerunner 1  Classical Thin Airfoil Theory /  1.1  2  of this approach was a method by Garrick (1936) who worked out a two-element configuration analogue to Theodoreen's (1931) analytical theory for arbitrary single-element airfoils. A modern example is seen in a paper by Halsey (1979). Secondly, there are the surface singularity panel methods, such as Hess and Smith (1967), which represent airfoil surfaces as a series of discrete line segments over which source strength and vorticity are appropriately distributed. A particularly efficient version of this approach which uses only vorticity, employs the stream function as fundamental variable, and uses a trailing point Kutta condition is given by Kennedy (1977). Of the two approaches to the numerical potential flow problem, the surface singularity methods, especially Kennedy's method, are the simplest and are probably as fast as the others. It is a version of Kennedy's method, presented in Appendix A, which is used in this thesis to calculate exact-numericalf potentialflowtheory predictions of flows. The success of the above numerical approaches depends on access to a large computer and, without expensive parameter studies, they give solutions to just one problem— trends are not easily established. The object of this thesis is to attempt an analytical description of these multiply connected flow fields, with the aim of obtaining a better qualitative understanding of the problem and simplified but realistic aerodynamic force predictions that can be calculated on today's inexpensive, readily available programmable hand calculators or small microcomputers. Such attempts have been made in the past, especially during the era of the slide rule dependent engineer who had no choice but to develop theories that were as analytical as possible. However, these earlier multi-element airfoil theories were not motivated by small gap slotted flap type configurations but rather by the large gap biplane configurations prevalent at the time. In order to get simple solutions, either approximations were made which are incompatible with a small gap analysis (Munk, 1922b; Millikan, 1930; Glauert, 1926) or some geometric considerations, such as camber and different airfoil element chord lengths, were ignored (Glauert, 1926). Those that were general (Garrick, 1936)  f By which is meant a numerical approximation to exact.  Classical Thin Airfoil Theory / 1.1  3  are complicated enough that they do not compare favourably with the modern numerical methods. The main reason for this complexity is that the flow field around multielement airfoils is multiply connected. A one-element airfoil flow is singly connected^ so that, at worst, it is described in terms of the singly periodic trigonometric functions. O n the other hand, a two-element airfoil flow, such as one about a slotted flap or leading edge slotted slat airfoil, is doubly connected and a general description of it would require use of doubly periodic elliptic functions which, numerically, are much more difficult to evaluate. As will be seen, however, elliptic functions have properties that allow them to be manipulated in much the same way that trigonometric functions are. In particular, the theory of the theta functions (Whittaker and Watson, 1927) used in conjunction with a programmable hand calculator allows rapid numerical evaluation of most elliptic functions. The desire to obtain a simple, useful, and general multi-element analytical theory is also motivated to a large extent by the success of the well-known single-element linearized airfoil theory, or thin airfoil theory. It originated in the early part of this century when Max M . Munk (1922a) first used a linearized boundary condition approach on the ideal flow.about a simple airfoil. He was able to describe the airfoil's overall lift and moment in terms of simple integrals. Subsequent contributions by Birnbaum (1923) and Glauert (1926) resulted in a much improved general formulation of the problem. Others applied the theory to special situations such as a simple attached flap, and calculated localized quantities such as flap lift or flap hinge moments. Much of this work is summarized by Abbott and von Doenhoff (1959). Stewart (1942) gives a particularly efficient formulation of the general problem using complex analysis and solving for a complex disturbance velocity. He shows how it is just as easy in thin airfoil theory to specify a pressure distribution and work out a camber line or thickness profile as it is to do the converse. His methods are similar to the ones used in this thesis. |  The definition of connectivity used is given in Woods (1961): a region is singly connected if it has one contour  separating it from its complement domain, doubly connected if it has two, etc.  Classical Thin Airfoil Theory  /  1.1  4  One-element thin airfoil theory lift and moment predictions are typically simple analytical expressions which agree quite well with experiment if boundary layer growth and separation effects are minimal. In fact, this agreement is usually better than that of exact potential flow theory which tends to overestimate loadings. Pressure distribution predictions, while generally less accurate and often containing singularities at points where the linearization process breaks down (such as leading and trailing edges), do show good overall qualitative agreement with experiment and are also given by simple expressions. The reason for the theory's simplicity is that application of the linearized boundary condition results in the geometric amplitude parameters for airfoil angle of attack, camber, and thickness appearing linearly in theflowsolutions. These three geometric effects, which completely determine the airfoil profile and its orientation in the flow, can then be shown to contribute independently to the aerodynamic characteristics of the airfoil. The effects are solved for in separate problems, the solutions of which are linearly superposed as required. One-element thin airfoil theory then is an efficient analytical tool that has been and continues to be a very useful aid to both the student and practising aerodynamicist. Its simplicity is one very good reason for using it as a starting point for developing an analytical multi-element thin airfoil theory. However, for this one would also want to anticipate that a multi-element extension of the theory would give predictions as realistic as those of the one-element theory. At first this seems unlikely since multi-element configurations generally operate with large mean line curvature, as in Figure 2 for example (page 13), so that the linearization requirement, that deviations of the boundary from some average straight line through the airfoil be small, is violated. Undoubtedly, neglecting the second order effects of large angles of attack or flap deflections would eventually be a limiting restriction on the use of a multi-element theory. However, as will be seen in the next section, these effects are not the limiting constraint on one-element thin airfoil theory and, in some cases, may not even be in a multi-element extension of it.  Classical Thin Airfoil Theory  /  1.2  5  This thesis then is devoted to the extension of classical thin airfoil theory into the multi-element regime. Only partial success is achieved but this includes a general tandem two-element thin airfoil theory, extensions of its overall lift and moment results to an arbitrary number of in-line airfoils, and some two-element staggered airfoil results. Also, a method is developed that greatly facilitates the calculation of multi-element thin airfoil theory lifts and moments. Based on Munk's integrals (section 4 of this chapter), in principle the method is applicable to any multi-element configuration. Comparisons are made between the tandem wing thin airfoil theory and exact-numerical potential flow theory. Variations of the aerodynamic forces with slot size and with flap chord size are examined. All of the aerodynamic forces for the tandem wing thin airfoil theory are calculated on a small, portable, inexpensive microcomputer. Some of these results are as simple as those from one-element thin airfoil theory.  1.2  T H E POTENTIAL FLOW - REAL FLOW CORRELATION  As an example of the relationship between one-element thin airfoil theory, exact potential flow theory, and the real flow consider Figure 1. Here, the nondimensionalized lift characteristics of a NACA 23012 airfoil with a 20% plain flap are shown for various angles of attack and flap deflections. The experimentally measured points in the figure are from Abbott and Greenberg (1939). The exact-numerical potentialflowcalculations were made using the theory described in Appendix A. Airfoil geometry was obtained from Abbott and von Doenhoff (1959). The predictions of one-element thin airfoil theory shown in the figure are described in section 1.4. There are several characteristics worth noting in Figure 1, as follows. a) The relationship between thin airfoil theory and exact potential theory. The mathematical relationship between the exact and linearized theories is more formally presented in the next section of this chapter. They are both two-dimensional,  Classical Thin Airfoil Theory /  1.2  6  incompressible, potential-flow theories. Essentially, the difference is that whereas the exact theory applies the tangential velocity boundary condition exactly on the airfoil boundary, the linearized theory applies it instead on some average straight line through the airfoil parallel to the flow at infinity. This is the same as assuming that deviations of the airfoil boundary from the average line are small, so that second order terms in a and r} for example, can be neglected. t  The errors that this assumption gives rise to are seen in Figure 1.  This  N A C A airfoil is 12% thick (maximum thickness relative to chord length) and it is this thickness which is primarily responsible for the difference between the exact and linearized theories at low values of a and T}. If the thickness were reduced to, say, 1% then the exact theory would collapse down virtually on top of the linearized theory (camber only being about 2% for this airfoil). For high values of a and rj, this collapse would not be as complete. Note that the angle of zero lift predictions, a^, are virtually identical for the two theories, so that differences between the theories are due to different lift curve slopes, dCi/da.  Abbott and von Doenhoff (1959,  section 3.5) show how the second order effect of thickness on the lift curve slope for a Joukowski airfoil is given by the factor ( l + A £ ) , where (t/c) is the thickness to chord ratio and A w 0.77. For typical airfoil section thicknesses of 12% to 18% then, one would expect that the exact predictions would be greater than the linearized ones by approximately 9% to 14%, which is consistent with Figure 1. Thus, there is a consistent relationship between one-element thin airfoil theory and exact potential flow theory which is, of course, completely independent of real flow effects such as boundary layer growth and separation.  Note too, from  Figure 1, that second order effects of o and rj are minimal for the ranges considered since the variations of the exact theory's predictions with both a and r\ are almost linear. For larger values of a and rj, these effects eventually result in the exact theory being overpredicted by the linearized theory.  FIGURE 1  Lift characteristics of a NACA 23012 airfoil with a 20% plain Bap.  Classical Thin Airfoil Theory /  1.2  8  b) The correlation between potential flow theory and experiment. Abbott and von Doenhbff (1959) give experimental lift and moment results for many different airfoil section profiles that the United States National Advisory Committee on Aeronautics (NACA) has studied over the years. These results show that, for thicknesses of 12% to 18% (and with noflapdeflections), experimental values of lift curve slopes (measured in the minimum drag regime of the profile — i.e. well below the stall) are all scattered about the thin airfoil theory prediction with most of them within 5% of that prediction and 95% of them within 10% of it. This accuracy, when combined with less accurate predictions of  , results in thin airfoil theory  overpredicting lifts at the point of maximum lift to drag (typically 2° < a < 6°) by 5% to 15% — although some poorly designed profiles can be overpredicted by as much as 20% to 30%. It is very rare indeed to see thin airfoil theory underpredict realflowresults (those profiles that have their lift curve slopes underpredicted invariably have the magnitudes of their zero lift angles overpredicted) so that one can conclude that one-element thin airfoil theory predictions are usually closer to realflowresults than exact potentialflowpredictions are. This is, of course, purely coincidental. It comes from the errors introduced by the linearization of the exact potential flow partially compensating for the necessary boundary layer corrections to this flow. It is these boundary layer corrections ,among other things, that airfoil profile designers try to minimize. By so doing they can delay the onset of stall, minimize drag, and maximize the lift curve slope, all of which also happen to improve the correlation between thin airfoil theory predictions and the real flow results. In practice, however, one must use caution when relying on this correlation since, as is obvious from Figure 1, it breaks down without warning as angle of attack increases. This breakdown begins gradually but is completed abruptly as the airfoil stalls (dC(/da becomes negative). The stall is caused byflowseparation from the  Classical Thin Airfoil Theory / 1.2  9  airfoil's upper surface well up from the trailing edge. An airfoil that develops high lift gets that lift primarily from a large negative pressure peak on its upper surface just behind the leading edge. Theflowover this surface, then, experiences large adverse pressure gradients which promote boundary layer growth and eventually flow separation. The separation begins at the trailing edge and, as the lift increases, so too does the adverse pressure gradient, the separation point moves upstream, the drag increases dramatically while the lift curve slope decreases, and eventually the airfoil stalls. Stall can be gradual, which would be preferable, or it can be sudden as in the case of the 23012 profile. These boundary layer related flow characteristics can be controlled to some extent by choosing a profile with a thickness distribution and camber line that minimizes upper surface adverse pressure gradients, at least over the primary operating range of a. In other words it is possible to use inviscid theory, in conjunction with a designer's knowledge of how inviscid flow characteristics affect boundary layer growth, to minimize degradation of the real flow. Here one-element thin airfoil theory can and, in fact, has played an important role. The NACA 6-series airfoils, an improved design that followed the 230 profiles (Abbott and von Doenhoff, 1959), are constructed from camber lines calculated by thin airfoil theory from designer specified uniform load type pressure distributions. These camber lines are simple analytical expressions, c) T h e effect o f flap deflection o n the potential flow - real flow correlation.  Figure 1 shows how this correlation breaks down as the plain flap is deflected through 10 and 20 degrees. For deflections greater than about 10 or 15 degrees the experimental results of Abbott and Greenberg (1939) also show that the drag increases dramatically, due to flow separation from the flap, as the lift slowly increases up to a maximum AC7 (C/ at a = 0°) of 1.5 at i] = 60°. Abbott and von Doenhoff (1959) conclude that one-element thin airfoil theory . . .permits cala  culation of the angle of zero lift, the pitching moment coefficient, and the chordwise  Classical Thin Airfoil Theory /  1.2  10  load distribution with reasonable accuracy . . . for most commonly used wing sections . . . for flap deflections not over 10 or 15 degrees." But even then, thin airfoil theory predictions of quantities such as flap lift and hinge moment are "relatively poor" since they are localized results from a part of the flow experiencing pronounced boundary layer separation. Other one-element configuration high lift devices, such as splitflapsand leading edge flaps, have also been used but, while they do increase the lift, none of them can match the high lift to drag ratios obtained by the multi-element configurations.  Thus, it appears that one-element thin airfoil theory correlates well with realflowresults unless angle of attack or mean line curvature are large enough that boundary layer effects appreciably affect the overall flow. The linearized theory's neglect of the second order geometric effects of these large deviations is not the reason for the breakdown of the potentialflow- real flow correlation. One would expect then that an airfoil configuration that dramatically reduced boundary layer effects on the overallflowwould also increase both the degree of this correlation and the range of deviations over which it would apply. One of the reasons the flow separates so easily from the attached flap of Figure 1 is that, due to the adverse pressure gradient over the back of the main part of the airfoil, boundary layer growth and momentum loss in that boundary layer are already appreciable when the flow reaches the flap. For even small flap deflections the further adverse pressure gradient over the flap (resulting from the pressure peak at the deflected flap leading edge) easily causes the flow to separate. However, if the flap were detached from the airfoil enough that the relatively cleanflowfrom below the airfoil came through the slot and swept the wake coming off the back of the airfoil away from the flap, so that boundary layer growth over the flap had to begin anew, then separation from the flap should not be as imminent. This is what appears to happen.  Classical Thin Airfoil Theory /  1.2  11  Foster, Irwin, and Williams (1970) made an extensive investigation of the flow around an airfoil with a slotted flap. They measured the lift and drag on the two airfoil elements and compared the lift with exact-numerical potentialflowpredictions. They also investigated the size and interactions of the boundary layers developed over the airfoil and flap. The nose of their flap was tucked under the trailing edge of their airfoil, as is usually done in order to drop the flap below the wake of the airfoil, and the effect of variation of overlap (the distance the flap leading edge is ahead of the airfoil trailing edge) on the aerodynamic forces was examined. The effect was minimal. Large effects were found, however, for variation in gap (the vertical distance between the lower surface of the airfoil and the upper surface of the flap). It was found that for very small gaps and relatively large flap deflections (TJ = 30° for their 31% flap), the wake from the airfoil was merged with the boundary layer over the flap, the drag was high, and the exact-numerical inviscid lift greatly overpredicted the experimental lift. Although the wing had not necessarily stalled, the flap was experiencing a great deal of separation. As the gap was increased, the airfoil wake and flap boundary layer began to separate from each other and eventually there was a region offlowbetween them that attained the free stream energy head of the flow. As this happened,flowseparation from the flap was reduced, drag was reduced, and lift increased. After the wake and boundary layer were effectively separated by the small region of inviscid flow, the experimental lift variation once again correlated well with the inviscidflowpredictions. The report concluded that, while both boundary layer and inviscid effects were important in obtaining the optimum flap position (which for this configuration was about a 2% gap), once in that position the airfoil wake and flap boundary layer would be separated and the overall real flow would be strongly dependent on the inviscid solution for the flow.  In conclusion, exact potentialflowtheory and, even more so, thin airfoil theory give realistic predictions of the real flow forces on a one-element airfoil if boundary layer  Classical Thin Airfoil Theory /  12  1.3  effects on the overall flow are small. Multi-element airfoil configurations are specifically designed to minimize boundary layer growth while operating with greater mean line curvature, and so it is not surprising that exact potential flow predictions of their forces also give good agreement with experiment, but for even higher lift configurations. It is reasonable to anticipate that the consistent relationship between one-element thin airfoil theory and exact potential flow theory would extend to multi-element configurations where one could then anticipate that a multi-element thin airfoil theory would give better predictions of the real flow forces than exact potential theory does, and that a close correlation with the realflowwould apply for even higher lift configurations than it does with one-element thin airfoil theory.  1.3  T H E LINEARIZING P R O C E D U R E  Consider the N airfoils in the steady, two-dimensional, incompressible, potential flow of Figure 2. The airfoils all have arbitrary positions, chord lengths, angles of attack, thickness, and camber. The boundary of airfoil j, say, is given by yy(x), its chord by Cy, and its angle of attack to the horizontalflowat infinity, u^,, by cry. The complex velocity, w[z) = u — iv, representing the exact potential flow about these airfoils can be written in terms of the flow at infinity and a nondimensional complex disturbance velocity, xv'(z) = u' — iv': w(z)  = u (l + 0O  V J \ Z ) ) .  The conditions which determine w\z) exactly and uniquely are: 1) tangent flow on all airfoil boundaries  2)  VJ'(Z)  v  dx  1 + u'  vanishes at infinity.  3) the Kutta condition applied to the airfoil trailing edges. The forces on the airfoils can be obtained through their surface pressures which are known via Bernoulli's equation and the pressure coefficient:  Classical Thin Airfoil Theory  /  1.3  13  z-PLANE  FIGURE 2  N arbitrary airfoils.  where p is the fluid density. In multi-element thin airfoil theory it is now assumed that the airfoils are thin and that they operate at small angles of attack. Thus, if y  Ja  is some average horizontal  line through a i r f o i l a s in Figure 2 , it is assumed that:  <  (i.i)  I.  Classical thin airfoil theory then makes the following corollary assumption:  u'l.kl  du'  d(y/cj)  du' d(x/cj)  dv'  <  IT  (1.2)  Although there are localized violations of this last assumption (at stagnation points and leading edges for example) that can lead to localized bad predictions (infinite pressures at stagnation points and leading edges), overall the assumption is good so that, as was previously stated, the overall force predictions and the overall qualitative nature of the  1.3  Classical Thin Airfoil Theory /  14  pressure distribution have close relationships with the exact theory. It is to be expected, however, that the linearized theory's predictions of more localized characteristics will be less accurate. A good example of a badly modelled localized characteristic is the Kutta condition, which is needed to determine the circulation around an airfoil in potential flow. Descriptions of it that are often used in exact potential flow are: a) finite velocities are required at trailing edges. b) finite pressures are required at trailing edges. and for exact-numerical potential flow a favourite is: c) there can be no load on the trailing edge. In exact potential flow these descriptions all amount to the same thing. In thin airfoil theory they can all be different and, in fact, none of them are correct Kutta conditions for all situations. This is a consequence of applying Equation 1.2 to the pressure coefficient and tangentflowboundary condition. Neglecting second order products in the disturbance velocities, as thin airfoil theory does, one gets: C = -2u'  (1.3)  p  (M. so that on the surface of airfoil ; : W  '(z) = u ' - i V = - 3 » - A . V  '  2  dx  (i.5) v  '  Thus, even though Equation 1.1 may be properly satisfied, in thin airfoil theory: A) if the airfoil is a thin curved plate with a weak infinite slope at the trailing edge, then a) and c) are violated while b) is true. Here, a logarithmic singularity at the trailing edge satisfies the boundary condition on v' but makes u' discontinuous across the trailing edge—hence the load. An example of this is the uniform load camber line used for some of the NACA 6-series airfoils.  Classical Thin Airfoil Theory / 1.3  15  B) if the airfoil is a long thin ellipse and so has an infinite slope at the trailing edge, then only a) is violated. Here, thin airfoil theory requires a simple pole at the trailing edge. The pole is entirely contained within the imaginary part of w'{z) and, in fact, u' is constant everywhere on the surface of the ellipse. C) if the airfoil is conventional and, due to thickness, has an acute nonzero trailing edge angle, then a) and b) are violated while c) is true.  In this case the boundary  condition requires that v' have a discrete jump at the trailing edge. This is provided by a logarithmic singularity which results in an infinite but continuous u \ D) if the airfoil's trailing edge angle is zero and if it has a non-infinite slope at the trailing edge, as for a flat plate or the well known Joukowski airfoils, then a), b), and c) are all true. There are no singularities at the trailing edge. Perhaps the only statement that can be made to describe a general linearized Kutta condition is that there can be no infinite loads at the trailing edge. Of course, this statement would not prevent violation of the Kutta condition in an exact potential flow.  Consider now the Taylor series expansion of u' and v' about y = yy for a given z:  ^ ' J ^ ^ J + ^ - ' i . ) ^ ) * -  ( 1 6 )  Thus, a further consequence of the small deviation assumption is that the boundary condition for airfoil  can be applied on the horizontal slit y = y  dx/j/dx =  an error of order |(yy airfoil's surface, C  p  2  y )/cj\ Ja  J a  to within  (using Eqs. 1.1 and 1.2). Similarly, the pressure on the  = —2u', is taken as equivalent to the pressure on this slit to within  the same small error. In this multi-element linearization each airfoil elementhas its own independent y . As in one-element thin airfoil theory this procedure breaks down at only Ja  a few localized points where u' or v' are singular. Also, in practice it is found that the  Classical Thin Airfoil Theory / 1.3  16  linearized theory will often give sensible answers even though the deviations are not very small (in the strict mathematical sense). Now, if y' j[x') and yjy(as') give the camber line and thickness distribution e  respectively for airfoil j (written in terms of a localized coordinate system, Fig. 2), then the airfoil's upper and lower surfaces are: y'u = y'c + y't I  I  (1.7)  I  VL =  y -vt c  so that the linearized boundary condition on the surface of airfoil j is: v ' - ^ - - a + ^ ± ^ i ~ dx ~ dx' dx' • V  +  ±  (  (18)  L  8  )  This shows that amplitude parameters for angle of attack, camber, and thickness that appear linearly in the boundary condition will also appear linearly in the flow solution, w'(z). In fact, if w' (z), w' -(z) and w' (z) represent functions that solve the three problems for aj  ej  t  tj  airfoil j independently of both each other and the other airfoils, then the linearity of Laplace's equation ensures that the effect of airfoil j on the overallflowis simply: Wji*) = <M)  + 'cj< > + h< ) W  Z  W  Z  ( ) L9  and that the overall flow solution will then be: N u/(z) =  ] > > ; ( * ) .  (I.IO)  The linearity of Equation 1.10 is transferred through to the airfoil forces as well. Neglecting second order products in the disturbance velocities, the Blasius equations for lift (L), drag {D) and moment (Af) become: t  D-iL D-iL  - j — —  n  ./ . , <b w{z) dz  = 2i  tt  }/> oo  M n  =2*< j> zw'{z) dz C  where M is the moment about the origin of the coordinate system and C is a curve Q  enclosing either the airfoil element of interest for localized forces or the configuration of interest for overall forces.  Classical Thin Airfoil Theory  /  1.3  IT  However, from exact potential flow theory it is known that the circulation (r) around some closed curve C and the net source strength (Q) within it is given by: -r + iQ = u j> w\z) dz 00  C so that in multi-element thin airfoil theory the drag on airfoil j is zero, assuming that it is a closed body. The lift and moment are:  CV-UA.)*  ( U I ,  The overall lift and moment expressions are the same as Equations 1.11 and 1.12 except that the nondimensionalizing length is taken instead as the sum of the individual chord lengths. It is interesting that in thin airfoil theory the lift and drag on an airfoil, regardless of the presence of another airfoil, are just proportional to its circulation and zero, respectively: D = 0 3  if  Qj =  0  so that the Kutta-Joukowski law applies on a localized basis. In an exact potential flow, however, this is only true for an isolated body or for the overall forces on an isolated group of bodies, as can be seen from the localized lift and drag expressions developed in Appendix A.  Perhaps the greatest advantage of thin airfoil theory is its ability to describe complicated general problems in terms of a linear superposition of independently solved simpler problems, each of which deals with only one aspect of the boundary condition. The individual linearly superposable force solutions are, for a given orientation of the airfoils of  Classical Thin Airfoil Theory / 1.4  18  Figure 2, simply constants multiplied by their respective geometric amplitude parameters. This permits a rapid investigation and easy understanding of the relative effects of these geometries on the forces. An understanding of the effects of the relative positions of airfoil elements in a configuration (interference effects) is obtained by simply looking at the variation of these force solution "constants" with position. For example, the variation of C /a with e  gap for a slotted flap airfoil would be of interest. Investigations of these types are presented in Chapter 4.  1.4  O N E - E L E M E N T T H I N AIRFOIL T H E O R Y FUNCTIONS A N D FORCES  What follows is a brief summary of the one-element thin airfoil theory problem. The functions used to solve the multi-element problems are just extensions of the simple functions presented here. Figure 3a shows the slit in the z-plane on which the boundary conditions for the airfoil of Figure 1 would be applied. Some problems are more easily solved in the f-plane where the flow about the slit has been conformally mapped into the flow about a unit circle using the well-known Joukowski transformation: 4z-2 On the surface of the slit f = e  t9  = c + ^.  (1.13)  and so: 2 z - l = cos0.  (1.14)  This transformation separates the upper and lower surfaces of the slit (they become, respectively, the upper and lower semi-circular arcs of the unit circle) with the result that simpler functions will solve a problem in the f-plane than would be possible in the z-plane.  Classical Thin Airfoil Theory  /  1.4  19  i  z-PLANE  f- PLANE  , dyu dx  dyL  f  1  dx  a) Linearized  physical  FIGURE 3  b) Solution  plane.  Joukowski  transformation  for the one-element  plane.  problem.  As previously mentioned, the general problem can be broken up and solved as a series of independent simpler problems: T h e a-Problem If a is the angle of attack of the airfoil, then:  dx  dx  = —a  and the solution is: 2ia  w' {z) = ta Q  7+T  (1.15)  which is correct and unique since: 1) on the slit surfaces y/z — 1 is imaginary so that v' = —a. 2) w' (z) -* 0 as z —* oo. a  3) the Kutta condition is satisfied. On the surface of the slit: w' (z) = a tan - + ia a  It  so that the pressure is singular at the leading edge—positive infinite on the bottom surface [9 = — ir) and negative infinite on the top surface (9 = TT).  /  Classical Thin Airfoil Theory  1.4  20  The lift and moment can be evaluated using the residue theorem and the expansion: ,, . x  iot ia  ^_ (/ 11 \\  ia t'a  as z —• co.  Letting the closed curve C expand out to infinity in Equations 1.11 and 1.12 gives:  C = 2ira t  ir  (1-16)  The aerodynamic centre of the airfoil, z , is that point about which the pitching moment ac  does not vary with a. That is, dC /da  = 0; and since:  mac  Cmac — m C  t  h  e  n  _  :  X  a  c  dC  ~ ~  + t ac C  0  m  x  Jda  dCt/da  _  7  (l-l ) 1 (  ~~A  L  1  8  )  The r;-Problem When w'(z) is written in terms of the f rather than the z variable it can be much simpler, as seen in Eq. 1.15. This is useful when writing down the solution to the more complicated flap deflection problem, where: v' =  0  for  0 <x  v' =  —T)  for  / < x < 1.  < f  The solution is: 2i9  c-e f i0  f  IT  j.+  l  '  (1.19)  i0  $-e- J.  which satisfies all the necessary conditions. The discrete jump in v' at the flap leading edge, z = / , is provided by logarithmic singularities in u' on the airfoil's top and bottom surfaces at this location. Note that when the flap is deflected the strongest resulting singularity in the flow is at the airfoil leading edge and not the flap leading edge, so that most of the lift from the flap deflection is borne by the front of the airfoil.  Classical Thin Airfoil Theory / 1.4  21  Now: as f —• oo  and:  ri °(i) •  as z —• co  +  so that:  C  t  = 2(9 f + sin9 f)t}  = 2 [cos- (2/ - 1) + 2^/(1 - /)] r,. 1  (1.20)  Similarly, using higher order terms from the above expansions: C  mo  = -^{9 + 2sm9 f  (1.21)  + lsm29 )  f  f  so that: C  mac  + cos9 ) = -2r//v7(l - f) .  = -\^9 {\  f  f  (1.22)  The Camber and Thickness Problems Both camber lines and thickness distributions can be represented in a general fashion through the use of Fourier series. The problems are orthogonal in the sense that where the camber boundary condition is even in 9: dyv  dy  dx  dx  c  dy  dy  dx  dx  L  '  c  resulting in pressures that are odd in 9, the thickness boundary condition is odd in 9: dy _ dy dx dx v  t  _ dy^ _ _dy ' dx dx  L  which results in pressures being even in 9, so that thickness has no effect on lift. If for camber:  then the boundary condition becomes:  ^ =f  +  tt—1  f>"C0Sfl*  U-24)  n=l  so that:  n  B = - f ^cosn9d9 ir J dx n  ;  n = 0,l,2,...  (1.25)  Classical Thin Airfoil Theory  / 1.4  22  This Laurent expansion accommodates virtually all conventional camber lines, including those requiring logarithmic singularities in w' (z). Equation 1.19 is just a special case of c  Equation 1.23. Lift and moment are: = -ff(i? +  AC  tfj  0  t  ^m  f l c  =  4(^1 + ^2).  (1.26) (1.27)  If the thickness flow function is: A,  0 0  *)  (1-28)  = E^ n=l  then the boundary condition is: ±-f-  = ^A smn&  with:  (1.29)  n  1  r»=l K  •"1/  , sinnfldtf dx  ;  n = 1,2,3,...  (1.30)  The net source strength contributed to the flow by a closed body is zero so that, from Eq. 1.28, Ai = 0. Equation 1.28 is not completely general since it will not account for strong singularities in dyjdx  which are usually present at the leading edge (at the trailing  edge as well for an ellipse). These singularities are added if necessary but they make the evaluation of the A 's somewhat more difficult. Since thickness has no effect on lift in the n  one-element theory, however, modelling it is usually unnecessary. Note that with the camber and thickness problems it is just as easy to specify a pressure distribution and work out what the corresponding boundary should be as it is to do the converse. This is particularly true when the specified u' function is simple, as it usually is, because this allows use of closed form analytic functions such as were used in Equation 1.19.  Classical  Thin  Airfoil  Theory  /  23  1.4  Finally, then, the total force is just a linear superposition of all the individual forces:  dC, d  =  and:  CL  -n(B  a  d  +  (1.31)  V  + B{) + lira  0  = ACL  mac  dC,  + 2(9  f  + sin  9 )rj f  dCmae  n  mac '  a _ '  >  (1.32)  dj  = ^B  1  +  j  9 ) - 2fvm-f) 2  ri.  The M u n k Integrals If only lift and moment are required then it turns out that one can bypass finding the general flow solution and solve directly for the forces in terms of the boundary condition, as first presented by Munk (1922a). Using Eq. 1.25 in Eqs. 1.26 and 1.27 and changing variables one gets: l 0  0  Note that these expressions give the general one-element thin airfoil theory force solutions, including the a and r) problem solutions, since to£(z) (Eq. 1.23) is a generalflowfunction capable of handling the a and 17 problem boundary conditions. It is by understanding and generalizing these simple expressions that thin airfoil theory is successfully extended to the multi-element problem.  CHAPTER 2 Tandem Symmetrical Airfoils  This multi-element airfoil linearization is initially restricted to the simple geometry of tandem airfoils, as in Figure 4a. In this chapter only the effects of angle of attack and flap deflection are considered. Camber and thickness are dealt with in later chapters. The linearized physical plane, Fig. 4b, has the airfoil and flap boundary conditions applied on the same average line which is taken as the real axis. The airfoil slit stretches between 0 and s and the flap slit between / and 1. The only conditions on a and / are: s and / are real  ;  0 < 5< / < 1  so that this could also be the problem of an airfoil with a leading edge slotted slat. This linearization will not model the effects of overlap. Note that for in-line airfoils, thickness will again have no effect on lift. This is because source strength distributed along the slits, with which thickness could be modelled, would have exactly the same effect on u' (the pressure) on the tops of the slits as it would on the bottoms. This would not be true for staggered slits.  2.1  T H E MAPPINGS  The conformal mapping from the z-plane to the A-plane, Fig. 4c, is a special case of a more general staggered slit mapping used by Garrick (1936) for developing biplane airfoil 24  Tandem Symmetrical Airfoils /  2.1  25  V  z - PLANE  k  5  a) Physical plane  /  b) Linearized physical plane  0) -PLANE  'i  (c)  (e)  (d) The rectangle solution planes.  f - PLANE  f ) The annulus solution plane. FIGURE 4  oi - PLANE  The tandem airfoil mappings.  Tandem Symmetrical Airfoils  /  2.1  26  potential flow solutions. Consider the small rectangle A ^ , iK', 0, K, K + iK' A ^ of the y  A-plane. A Schwarz-Christoffel transformation is used to map the interior of this rectangle to the upper half of the z-plane, to corresponding points —oo, 0, a, / , 1, +oo. The rectangle below the /i-axis then automatically transforms into the bottom half of the z-plane in a similar way. The left and right sides of the big rectangle form both surfaces of the airfoil and flap, respectively. The mapping is given by the elliptic integral: A  =  _y7(i^) f  2  dz dz  /J x/z(z-s)(z-f)(z-l) ^/ziz-sXz-fKz-l) »  ''  ( 2 - 1 )  It can be reduced to a standard form by the transformation:  f(z — s)  o  •  ° = z(f-s) < I f  so that:  . A= /  2.2  ^  a  °  =sn~ (o,fc) 1  (2.3)  or. z  =  7Z-  l_L— 2 Bn  2  •  3  4  (-)  {Xfk)  k is the modulus of the Jacobian elliptic function sn(A,A;):  *  2 =  7[[_1)  > 0<*<1.  (2.5)  K is a complete elliptic integral of the first kind:  f  ,  K=f  —  a  .  (2.6)  0<Jfc'<l  (2.7)  If k' is the complementary modulus: fc' = l - * 2  = i^j  2  ;  then K' is the same function of k' as K is of k. When z goes to infinity A = A ^ , and from E q . 2.4:  Sn2  A  °° = J~  (2-8)  Tandem Symmetrica] Airfoils  27  / 2.1  where snA is written for sn(A,A:). If A,^ = /x^, ± iK' then the periodicity of the elliptic function gives: » snA^ =  1  so that, using Eqs. 2.5 and 2.8:  snVoo = l - s  5  c n  = «  2  dnVoo = y  J  (2.9)  where, in general, the Jacobian elliptic functions, sn, cn, and dn, are related via: sn A + cn A = 1 2  ;  2  dn A + k sn A = 1. 2  2  (2.10)  2  Thus:  The w-plane, Fig. 4d, is a simple scaling of the A-plane: u=  x  + i  1  =  ^  (2.12)  ITT  Woo = Xco ± y iK' T = — .  (2-13) , \ (2.14)  Similarly, the w'-plane, Fig. 4e, is just: u' = x' + ii = OJT' a 4 = ± ! + »7*o  T  '  =  %  '  (2.15) (2.16)  (  2  -  1  7  )  By definition the Jacobian elliptic functions use the A variable as argument. However, theta functions, which will be used later and whose properties are summarized in Appendix B, are more naturally defined with the u and w' variables. The rectangle is unusual as a solution plane since it separates not only the upper and lower surfaces of the slits as desired, but also regions in the flow across which the  Tandem Symmetrical Airfoils / 2.1  28  flow characteristics must be continuous. For example, the real axis in the z-plane (z < 0 , z > 1) occurs twice in the rectangle, along u = ±iK'. Therefore, the single valuedness of the flow must be ensured by requiring of any flow function, F that: t  F{\) = F(\ ± 2iK').  (2.18)  This periodicity is consistent with the mapping, Eq. 2.4. Thus, theflowabout the slits in the z-plane is entirely contained within the infinite strip 0 < / i < K of the A-plane where it is periodic in 2iK'. Although the Jacobian elliptic functions are doubly periodic functions, so that the mapping is also periodic in the fi direction, there is no reason theflowshould be doubly periodic and, as will be seen, in general it is not. Any problem with maintaining theflow'ssingle valuedness is avoided in the f-plane, Fig. 4f, where: c = g'V*"  (2.19)  R = */' e  (2.20)  K K  and letting ^ = — d: d= t^t  K  '.  (2.21)  It is particularly convenient that: f = e~ '. 2,w  (2.22)  Here the entire flow is in the annulus 1 < r < R. The airfoil slit has now become the unit circle with its leading and trailing edges at f = — 1 and f = 1 respectively; the slot s < x < / corresponds to 1 < $ < R; ( = R and f = —R are, respectively, the leading and trailing edges of the flap; 1 < z < oo corresponds to —R < f < — d and —oo < z < 0 to — d < f < — 1. The upper half-plane semi-circular arcs are the upper surfaces of the airfoil and flap slits. The £-plane will be the starting point for constructing the necessary twoelement flow functions from the elementary functions of one-element thin airfoil theory.  Tandem Symmetrical  2.2  Airfoils  /  2.2  29  T H E FLOW FUNCTIONS  Tangent flow on the airfoil and flap surfaces requires that: v' = —a  on the airfoil slit  v = —a — rj  on the flap slit.  (2.23)  From the one-element theory it is known that the function i/(f + 1) has its imaginary part constant on the airfoil. It is not, however, constant on the circle r = R. This unwanted effect on the flap can be negated by using Milne-Thomson's circle theorem to add an image singularity at f = — R . The function: 2  (+ 1  $+ R  .  2  has its imaginary part constant (and, in fact, zero) on the flap but the new singularity has now created an unwanted effect on the airfoil — which is again negated using an image singularity, this time at $ = — l/R . This process continues indefinitely and leads to an 2  infinite series which converges everywhere in the annulus. The function:  '« = { ? T T M  +  §[s^TT-?+W]}  -  (2 24)  has its imaginary part constant and equal to 1 on the airfoil and constant and equal to. 0 on the flap. A similar function, F  2l  can be constructed to deal with the flap boundary con-  dition while having no effect on the airfoil boundary condition. Since w'(z) must be zero a  * f = $<x» these functions are not yetflowfunctions. However, this is easily remedied by  subtracting appropriate constants. Rather than working with infinite series, such as Eq. 2.24, it is more convenient to work with conventional elliptic functions. Rewriting Equation 2.24 using the u/ variable, it becomes: ;2n  Fi=i-  _ i  tana;'-4^ q sin 2or ^-f l + 2g' + q( 7 ' 4 n ^-f 1 + 2q' cos2w' cos 2w' +  1  (- ) 2  25  2n 2n  fizz 1 where:  •  q' = " ' = l/R. e  iT  * (2.26)  Tandem Symmetrical Airfoils / 2.2  30  This can now be summed in terms of the second theta function and its derivative: *i=« +  (2.27)  £V,g').  An invaluable aid to the user of elliptic functions is Whittaker and Watson (1927) where the four theta functions, their properties, and their relationships with the Jacobian elliptic functions are clearly presented. A brief summary is provided in Appendix B and this also serves to define the elliptic functions as they are used here. 0' One of the important properties of the -^(w',?') function is that, while its v  2  series representations, Eqs. 2.24 and 2.25, converge rapidly as the slot (/—s) becomes large (so that R —*• oo, q' —• 0), they lose their ability to converge as the slot becomes small (R,q' —* 1). This characteristic is common to the flow functions used in other analytical theories (Williams, 1973; Halsey, 1979), where it is overcome by using the computer to evaluate a large enough number of terms. However, the problem can be avoided entirely through the use of Jacobi's imaginary transformation, Eqs. B6, which allow F to be 1  rewritten:  2u  K'0\,  ir  K *i  ,  (2.28)  where: (2.29)  e* = e~* 'l . iT  K  K  As the slot goes to zero, K —* JT/2, K' —* oo, and q —* 0 so that the series representation of Fi in terms of the u variable and q parameter is ideal for small slot configurations. The four theta functions are very similar to the trigonometric functions sinw and cosy, which  and i?2(w) behave like as q —» 0. Indeed,  — #i(w + |) as  well. i?i(w) has a simple zero at the origin in the w-plane. t? (w) has this same zero moved 2  to a/ = rr/2.  $3(0/)  = i?4(w + f) has a simple zero at w = | +  and t? (u/) moves it to 4  w = jrr/2. Thus, the -^( ) functions are analogous to tanw and cotw. w  They have  simple poles in thew-plane at the four critical points of the tandem slit mapping. In the  Tandem Symmetrical Airfoils  / 2.2  31  z-plane they have square-root singularities at these locations. Also, as can be seen from Equations B l and B7, all four of these functions are purely imaginary on the surfaces of both the airfoil and flap (i.e. for u = »7 or w = | + »7), so that it's the second term of Equation B4 that enables the imaginary part of F , Eq. 2.28, to change from 1 on the x  airfoil to 0 on the flap. It is now possible to find theflowfunctions entirely in terms of theta functions in the w-plane. F has a singularity at the leading edge of the flap; its imaginary part is 2  equal to 0 on the airfoil and equal to 1 on the flap:  K' ;  2w  (2.30)  T h e a-Problem  v = -a a  on the surfaces of both slits.  As will be seen, a simpler solution is obtained by considering the overall angle of attack of the configuration as a single problem rather than by dealing separately with each airfoil element's boundary condition. The form of the solution will be: <{*) = C [F + F - [F + x  x  2  x  F) 2  ]  z=00  where C is some constant. But, using Eqs. B15: x  K' # 2 , F +F = i + i— #2 and using Eqs. 2.4 and 2.10: x  2  F + F = i-i x  ^  2  ^ 4 ,  ^  = 1  V  2K' 1-3  z{z - f)  .2K' snAdnA i  :—  JT  cn A  (2.31)  The square-bracketed expression is purely imaginary on the surface of either slit, so that: ii/ ( )=iajl0  2  \z-s)(z-\) z{z - f)  (2.32)  This satisfies all the required conditions and, as the slot goes to zero, it reduces in a simple manner to the corresponding one-element flow function, Eq. 1.15. In the last section of this chapter, this flow function is extended to include an arbitrary number of in-line airfoils.  Tandem Symmetrical Airfoils  /  2.3  32  T h e rj' P r o b l e m  on the airfoil  ;=o  v  on the flap. The flow function will have the form:  F2) + CZF2 - [C^Fi + F2) + C F 1  «;(*) = C (F + 2  3  X  If:  _2  2  X  (2.33)  = {Xoo) then, using Equations 87 and B3, 3 it can that T F^u^) ^3> ~jTbe shown ~K 3 =  l  8 0  that:  n  •(»-,)(«-I)"  <{*)  +  ^ ( 7  I ^  z{z - f)  (2.34)  T H E FORCES  2.3  Lift and moment are obtained with Blasius' equations, Eqs. 1.11 and 1.12. In order to evaluate the integrals using the residue theorem the following expansions are necessary. As z —• 00: 1  3 + 1  I  z{z - f)  l-(/-s)(3/+s-2) 8z  -f 2z  2  t  •(A)  (2.35)  and using Equations B17 and 2.1: , G -G  2u K'0' , M  2  2  l  , G ( +l+/)-G (3/+3-l) , 8z 2  3  1  2  (2.36)  where: „  l  t  2K'  S  rr  (2.37)  E' is a complete elliptic integral of the second kind (Appendix B). Calculating the overall forces is now straightforward. For the forces on the airfoil alone, an airfoil selector function, A{u), is used: 2u Alui) = 1 1:  K' — - r M w ) = 1 + t'odd(7) K V = 0 +1 odd(7)  on the airfoil (2.33) on the flap.  Tandem Symmetrical Airfoils /  2.3  33  The notation simply indicates that the imaginary part of A(u>) is odd in 7 on the slit surfaces. As z —• 00:  ^ , ~ I -  B  4  +  (^)GI-OJ W L  (~)  G,(3 +/+l) - G,(«+/+l) _ Q_J 3  +  J  0  (2.39)  where:  , , *  »  ^  + f  ^(Xoo) = E  3  + ^  .  (2.40)  Now, for any function F(z), Cauchy's Theorem gives: jF{z)A{u)dz C  = j F(z)dz + jF(z)[A(u)-l]dz Ct  j  +  Ci  F{z)A{u)dz  (2.41)  Cf  where C encloses the tandem airfoil configuration, C encloses just the airfoil, and Cy t  encloses just the flap. For the a-problem the airfoil lift is:  2  C, =  --Jw' (z)dz= -^.JF{z)dz  t  a  c.  c,  where: F(z)  z(z - f)  = 0 + 1 odd(7)  on the airfoil (2.42)  = 0 + »odd(7)  on the flap.  Since A(OJ) has a square-root pole only at z = s, each contour integral in Equation 2.41 can collapse onto the slit surfaces — i.e. there are no simple poles in the integrands at the slit leading or trailing edges (or anywhere else) that need to be accounted for. Therefore, since dz is also odd in 7 on the slits, the last two integrals of Equation 2.41 are zero due to the oddness of their integrands and:  *-*/™">* C  which is easily evaluated using Equations 2.35 and 2.39. C by replacing F(z) with zF[z) in the above procedure.  miQ  is handled in the same way  Tandem Symmetrical Airfoils  /  2.3  34  For the ^-problem this procedure will not work because the real part of V — -r-[u)  =  say, is not zero on both slit surfaces. By integrating by parts  first, however, a function is obtained which has this necessary characteristic. Since 7(w) is periodic in JTT, one gets:  j  -  7(u) dz =  C.  j(z - f)^- dz C  t  where / dz is taken as z - f to handle the singularity in 7(u) at z = / . Now: ( z - f ) ^ = -g^* " ^ " ^ " 7 ) = 0 + iodd( ) "dz 2 /z(z - s){z - f)(z - 1 )  on the airfoil  7  (2.44)  K l )  k  s  = 0 + iodd(7)  on the flap.  Identifying this function with F(z) from Equation 2.41 results in the last contour integral of that equation being zero. The second last integral is not zero since the square-root pole at z = s in A(ui) combines with that in Equation 2.44 to form a simple pole with residue: v , M ^d7 , f-s 2K' ^(z-s)(z-f)-A(u) = - — — '2K' (2.45) 2E'' A l  IT  so that for the r/-problem: (z-1)-  a  G (z - f) 2  ^(z - s)(z - f)(z - 1)  s [J  + 2iri(f-s) which is now solved by residues. C  2K'  A(UJ) dz  2K'  2E'  }  (2.46)  is evaluated the same way by taking / zdz as  mso  (z — f )/2 when integrating by parts. 2  2  Note that A(ui) is not the only airfoil selector that could be used in Equation 2.41. The function: 1  2UJ  K'  IT  -M  also has the necessary characteristics and has a square-root pole at z = 1. Analogous functions constructed from t ? 2 (  w  )  a  n  <  ^ $i( ) u  a r en o  * used, however, since they have square-  root poles at the flap and airfoil leading edges, respectively, which would combine with those in F(z) (Eq. 2.42) to produce simple poles in the £-plane.  Tandem Symmetrica/ Airfoils  a) b)  dC  d(h  G  S + 3  dr, -  da  d)  da dr,  0  ^  h)  - G\  «+W..  s  3  _  3 + 1 - / dC  2  l  i dr,  - G  1  ) ( 1 - E  3  ) -  7  ^ G  2  1-f  da  t  _ 3 + 1 - / dC  t  2  t. da ' 1-f da  1-f  da  3  t  dr,~T=J  dC  tt  dr,  4(/-3)(l-/) 2  da  (s+1-/)  +  (2.47)  2  dC, mo dr,  3  »0  m  G -'^-J-G  +  z  E —^ + — ( G oa a  j  dC i)  2  ( +l-f){l-V )  da m o g)  35  3  dC,  dC  2.3  = 2TT  t  da  c)  e)  /  _  da  2(3+l-/)  da  ^{[1  2  [(3+l + / ) G - ( 3 / + 3 - l ) C ? ] 2  1  - ( / - ) ( 3 / + 3 - 2 ) ] ( l - E ) + ( 3 / + 3 - 5 ) Gl a  3  + (3/- -l)G } S  dC„  2  dC„  dr,  j)  -^0\- G - f) {  2  Gi  da  5a  da dC, 1)  m  /o  _  dr,  /3+l-A 5G 2  V  1-/  ;  m o  dr,  /  3 \  2  1-/  5a  dc,  » 0 _j_ J v \ l - f j dr, • 1-f dr, m  Equations 2.47 summarize the forces for the a and r, problems. They are all in terms of the three constants G j , G , and E which are easily calculated with Equa2  tions 2.48. G  m / o  3  is the flap moment about its own leading edge.  Equations 2.48 are found in or derived from Whittaker and Watson (1927) and are all presented in Appendix B. As the slot in Figure 4b goes to zero, so too does q. For most airfoil flap configurations, q < 0.1 and the q terms in Equations 2.48 can be A  Tandem Symmetrical Airfoils /  k  b)  2e =  =  2  ft  it  (l + ^ C l + N/ib )  3  0  2K  f)  2  g = e~*£- = c + 2eo + 15*0 + 150e^ + O(l000cJ ) = (1 + 2q + 2g + 2g + 2? 4  d) e)  2  2  7  c)  36  Jfc -f-*' = l  a)  2  2.3  n \ n J  ir  22?'  9  16  7  + ••)  2  q  2  -|(a)+T k-(s)v,(«+v+v+^+-) :  (2.48)  g) e =  i V/ifc'  cos 2XCO =  1  *[l  -  ?4(2 -  - q (6 n  k)  -  20^2 +  32^)  7 6 £ 2 + 272cT 4 ~ 3 2 0 £ 6 ) +  0(g )] 16  0<2 <n Xoo  2K'  i) j)  4 ^ 2 ) + <? 8 (3 -  1-5  2 F"' Jr  2Xc /  VV7A  7+  7  +  g  1 2  sinl2  X o o +  0(  g  16  )]  neglected. For large slot configurations, evaluating the series up to and including the q  12  terms should give sufficient accuracy. If extremely large slots are of interest, then it  would be better to replace the q series with either q' series or the series of Appendix B which converge in powers of q . 16  Equations 2 . 4 7 and 2 . 4 8 have been programmed into a TI Programmable 5 9 hand calculator. It takes only half a minute to evaluate them to within an error O(l00(V ). 6  The results for dC /da and dC /da t  mo  are particularly simple, in keeping  Tandem Symmetrical Airfoils  / 2.4  37  z - PLANE  x„ x  X  2  5  FIGURE  22 *3  N in-line  X  33  NN  uncambered  airfoils.  with the simplicity of their flow function, Eq. 2.32, and they reduce to the one-element thin airfoil theory predictions, section 1.4, when f — s = 0 . The expression for the tandem airfoils' aerodynamic centre: dC Jda m  t  (2.49)  dC /dc t  is also simple while the pitching moment about the aerodynamic centre: dC ^mac  ~  mo  [  dC  dn  x  t  dr,  ac  s+l-f G  f  -( i-/)2{ 4^=7 5+  - Gi  .  2s s+l-f  - 1  (2.50)  is independent of E . 3  The rest of Equations 2.47 have also been shown to reduce to the one-element theory predictions for both f—s = 0 and f—s= 1. Comparisons between Equations 2.47 and exact-numerical potential flow theory are presented in Chapter 4 .  2.4  N IN-LINE U N C A M B E R E D AIRFOILS  The a-problem for the N in-line airfoils of Figure 5 is solved by: w' (z) = ta j ; a  (z - x ){z - z ) • • [z - x ) n  (z-x ){z-x )---(z-x ) l  2  i  NN  22  N  }  (2.51)  since, again, the square bracketed expression is purely imaginary on the slit surfaces, it goes to 1 as z —*• oo, and it is zero at all the trailing edges.  Tandem Symmetrica] Airfoils  / 2.4  38  The overall lift is: — Ct/ / =  r <coefficient of - in w'(z) as z —> oo z  —  7  which gives: C = 2rca.  (2.52)  t  The overall pitching moment (about the origin) coefficient is given by:  r  °  m o _  4  n  coefficient of -=• in w'(z) as z —• oo  i  (c ...+c„)» 1+  If the slots are defined by: s = x 0  x  s = x 1  s  2  -x  2  =  n  etc.  x  x  3 ~ 22  then: Cm  0  = ~Y  j + ( ., 1  Ci+  4  +  C j y  ) 2 [ 0 l + («0 + 8  C  Sl)c  2  + • • • + («0  + * • - + SiV-l)C/v]  J.  (2.53) If the origin in Figure 5 is at the aerodynamic centre, then C  mQ  = C  mac  = 0, a = 0  —  x  oc>  and:  * = (ci+'-'+c^) ac  1  3 i C + ( 3 i + 3 ) c + - - - + (3i+--- + 3 - ) c ^ 2  2  4  3  N  (ci+-+cNy  1  (2.54)  Surprisingly, these simple, general results for the overall lift and moment do not seem to have been presented before, although Glauert (1926) shows that the exact inviscid lift for tandem flat plates of equal chord is C = 2jrsina. t  Richardson (1981)  also arrived at the results Equations 2.52 and 2.53 at about the same time as the group here at U.B.C. He also found A C ; for a simple parabolic camber line — a result which is duplicated and generalized in the next chapter.  CHAPTER 3 General Tandem Thin Airfoil Theory  For tandem airfoil elements, the thickness distribution analysis and camber line design problems are solved with the same elementary functions and transformations used in oneelement thin airfoil theory. The first two sections of this chapter present this part of the theory. The camber line analysis and thickness distribution design problems, however, are not so easily solved. They require the use of elliptic integrals or functions for a satisfactory description of their general properties. This is dealt with in sections 3.3 through 3.5. In the final section of this chapter, the general camber functions of section 3.3 are used to obtain generalized expressions for overall lift and moment which are the tandem airfoil extensions of Munk's integrals. In the sections that follow, basically two classifications of functions are presented. The first consists of general functions which are analogous to the one-element camber and thickness functions of section 1.4 (Eqs. 1.23 and 1.28) in that they use simple Laurent expansions to describe theflowfields.They provide good qualitative descriptions and are useful for handling numerically specified airfoil profiles. Although they can be used for any type of airfoil profile, the necessity of evaluating a large number of coefficients makes their use undesirable when a boundary condition is specified as an analytic function or, perhaps, as a group of analytic functions each of which applies over some discrete segment of the airfoil. Here, the second classification special functions is more useful. 39  Generai Tandem Thin Airfoil Theory /  3.0  40  F o r e x a m p l e , t h e ry-problems of b o t h one-element a n d two-element t h i n airfoil theory are solved w i t h special functions that satisfy one particular type of b o u n d a r y c o n d i t i o n , w i t h the o n l y variable o r coefficient b e i n g the a m p l i t u d e of that p a r t i c u l a r t y p e of d e v i a t i o n . O t h e r t h a n t r y i n g t o s u m t h e series o f a g e n e r a l f u n c t i o n a f t e r h a v i n g e v a l u a t e d its coefficients, there is no general w a y of finding the special f u n c t i o n for a particular boundary condition.  Experience, simplifying the p r o b l e m to the analogous one-  e l e m e n t t h e o r y p r o b l e m , a n d , m o s t i m p o r t a n t l y , u n d e r s t a n d i n g t h e n a t u r e of t h e s i n g u l a r ities t h a t m u s t p r o d u c e t h e desired result are t h e best aids t o d o i n g it. T h e g e n e r a l r u l e t o follow w h e n s e a r c h i n g f o r a n a p p r o p r i a t e s i n g u l a r i t y is t h a t a d i s c o n t i n u i t y i n v' m u s t b e a s s o c i a t e d w i t h a s i n g u l a r i t y i n u ' , o r v i c e v e r s a . T h i s i s a c o n s e q u e n c e o f t h e a n a l y t i c i t y o f t h e flow f i e l d . U s u a l l y , t h e d i s c o n t i n u i t y i n v' c a n a l s o be associated w i t h a discontinuity i n the b o u n d a r y condition,  dy/dx, a n d t h e s i n g u l a r i t y  in u ' w i t h a singularity i n the airfoil element's surface pressure. T h e i m p o r t a n t  exception  to these a p p a r e n t l y s t r a i g h t f o r w a r d associations o c c u r s at l e a d i n g a n d t r a i l i n g edges. Iu t h e a - p r o b l e m f o r e x a m p l e , E q . 1.15, v' — —a i s p e r f e c t l y c o n t i n u o u s e v e r y w h e r e o n t h e b o u n d a r y , even a r o u n d the l e a d i n g edge. edge.  H o w e v e r , u ' h a s a s i n g u l a r i t y at t h e l e a d i n g  I n t h i s c a s e , t h e d i s c o n t i n u i t y i n v' o c c u r s b e t w e e n a p o i n t o n t h e a i r f o i l at t h e  l e a d i n g edge a n d a point i m m e d i a t e l y i n front of the l e a d i n g edge, while t h e associated singularity i n u ' is entirely o n the surface of t h e airfoil.  There is, of course, a physical  r e a s o n f o r t h e d i s c o n t i n u i t y i n v' a n d if t h i s c a n b e a n t i c i p a t e d t h e n t h e p r o b l e m r e s o l v e s itself. T h i n airfoil theory pressure singularities c a n also be thought  o f as b e i n g  a s s o c i a t e d w i t h e x t r e m e s i n t h e e x a c t p o t e n t i a l flow p r e s s u r e d i s t r i b u t i o n , s u c h as t h o s e w h i c h o c c u r at s u c t i o n p e a k s o r s t a g n a t i o n p o i n t s . O n e m u s t b e careful here t o o , however, since not a l l s t a g n a t i o n points p r o d u c e pressure singularities  that exist on the airfoil's  surface. A g o o d e x a m p l e o f t h i s i s s e e n i n t h e flow a r o u n d a n e l l i p s e , a s p r e s e n t e d i n t h e following section.  General Tandem Thin Airfoil Theory / 3.1  41  iy n  z - PLANE  •PLANE  FIGURE 6  3.1  Solution planes for the thickness and camber problems.  T H E THICKNESS DISTRIBUTION ANALYSIS  F i g u r e 6 shows the m a p p i n g s to the  and  solution planes where the general thickness  p r o b l e m functions are c o n s t r u c t e d . In the ^,-plane, a J o u k o w s k i t r a n s f o r m a t i o n opens the a i r f o i l s l i t u p i n t o t h e u n i t c i r c l e w h i l e l e a v i n g t h e flap a s t r a i g h t l i n e s l i t o n t h e p o s i t i v e real axis:  (3.1)  (3.2)  -(1 + cos 0 )  on the airfoil surface.  S  (3.3)  S i m i l a r l y , i n t h e f - p l a n e t h e flap b e c o m e s t h e u n i t c i r c l e w h i l e t h e a i r f o i l is left a s t r a i g h t r  line slit — this t i m e o n the negative real axis:  1 + /  ^r x-/  1-f  f  1\  = ,-I^±V(*-/)(*-l) =  1-/  (1 + c o s Of)  o n t h e flap s u r f a c e .  (3.4)  (3.5)  (3.6)  Genera/ Tandem Thin Airfoii Theory  / 3.1  42  The general thickness problem function is: w[(z) =  A  (  *0 4-  fo  1  n  C  ,  A  A  ,X^ ±  L  +  +  A  fn  (3.7)  1  On the airfoil, the second group of functions in this equation is purely real so that the airfoil boundary condition is simply: C O  0. ± ^ = ^ t a n ^ + £A, sinn0. dx 2 2 '  (3.8)  n  n=l  with: =  i / te _ ir J  (3.9)  sin n9 d9 t  [dx  t  exactly as in one-element thin airfoil theory when the leading edge singularity in dy /dx t  is accounted for. A  tQ  other A 's ln  must be determined from the nature of this singularity before the  can be calculated. If the airfoil elements are closed bodies, then the net source  strength they produce is zero and A  90  + A  tl  — 0 = Af -f A^, which is precisely what 0  Eq. 3.9 gives for n = 1. The pressure on the airfoil is:  (3.10)  C„ = "2 Similar simple expressions are obtained for the flap.  For conventional airfoil shapes, Equation 3.7 is completely general. w[(z) is just a linear superposition of the functions which satisfied each airfoil element's boundary condition in one-element thin airfoil theory.  Since these functions are purely real  everywhere on the real axis except on the slit whose boundary condition they are required to satisfy, any number of them can be linearly superposed to solve the general thickness problem for multi-element in-line airfoils. Note, too, that Equation 3.8 is odd in 9 since ±dy /dx t  C  pti  t  must be, and that  Eq. 3.10, is even in 9 confirming that thickness does not affect lift for in-line airfoils. t  General Tandem Tbin Airfoil Theory / 3.1  43  In the following example, both leading and trailing edge singularities are used to illustrate the magnitude of thickness interference effects for tandem bodies. Consider the tandem ellipses: y = t y/x(s — x)  for the "airfoil"  y = t s/(x- / ) ( l - x )  for the "flap"  tt  tj  9  f  (3.11)  where tj is the appropriate thickness to chord ratio. Now: s  dy , dx t  -t  ~2 *\/x(s - x) X  and: dx  V(*  "/HI"*)  so that theflowfunction is simply 3  w[(z) = t.  ~2 \Jz(z - s) Z  1-  1  z—  1+ / (3.12)  which goes to zero as z —» oo. The singularities at the trailing edges of the ellipses are entirely contained within their boundary conditions, so that the pressures are: _ C ,= - 2 * , - 2 t p  ;  1+ / 0< x<s  1+ V(l-x)(/-x)  C =-2t -2t, p  f  s X 2 1 \/x(x - s)  (3.13) f<  x< 1.  These equations are plotted in Figure 7 for a 33|% flap, for t„ = tj = 0.15, and for various slot sizes (an effective midpoint in slot sizes, between the one-element thin airfoil theory configurations at /—s = 0 and /—s = 1, is suggested by the rectangle geometry, Fig. 4c, when K = K'). In one-element thin airfoil theory, the pressure distribution over an ellipse is just C = — 2t, which is consistent with Equations 3.13 when f—s = l. It is p  the square bracketed expressions in these equations which describe the interference effects.  General Tandem Thin Airfoil Theory / 3.1  FIGURE 7  44  Thin airfoil theory predictions of the pressure distributions 15% thick ellipses, ( l - / ) / ( s + l - / ) = 1/3.  When /—s = 0, the strong discontinuity in dyjdx  over tandem  at z = / results in the relatively strong  singularities in Eqs. 3.13 at this location. Note that despite the fact that these ellipses have stagnation points at their leading and trailing edges, the pressures on the surfaces of the separated ellipses are not singular here. The singularities do in fact exist in u' but they are immediately ahead of or behind the ellipses, as can be seen from Equation 3.12 when, for example, z —> 1 along the real axis from z > 1. A more useful example is the thickness problem for in-line airfoils with the same thickness distributions as those of the NACA four and five digit airfoils (Abbott and von Doenhoff, 1959). That is: y (x) t  = t (Ay/x + Bx + Cx + Dx + Ex*) 2  3  ,  0< x< 1  (3.14)  for an airfoil of chord length 1. If the conventional NACA distribution is adjusted so  General Tandem Thin Airfoil Theory / 3.1  45  that y (l) = 0 (instead of 0.01*), then A = 1.47792, B = -0.63000, C = -1.77370, t  D = 1.44700, and E = -0.52122. Now, scaling to the airfoil coordinate system: dyi  k . « .\A [ ^Is  dx  + fl + 1 0  f  , ( £ ) /' r „\ 2( £ )/' r] \ 3  r. + u  +  (  ,  1 5 )  and, from Equation 3.7, since: as x —» 0 on the upper surface then: = At,  K  and the evaluation of the remaining A, 's can be accomplished with Equation 3.9. This is n  not, however, the most efficient method. Since Equation 3.14 is composed of simple algebraic functions, it is easiest to deal directly with the particular flow functions associated with each of those algebraic functions. For example, the function: 1 TT  ,, Z—  z ln 3  1  z  o Z 1 + z2 + - + 2 3  has its imaginary part equal to +z on the upper and lower surfaces of a slit stretching 3  between 0 and 1, it goes to zero as z —* oo, and, as a one-element thin airfoil theory thickness function, it cannot violate the Kutta condition. The flow function for the boundary of Equation 3.14 is: w\(z)  =  -±  -4=  In ^-4  2  + (B + 2Cz + 3Dz  '  V 2y/z  y/z+1 +  2C +  In  z  + AEz )  —  Z  (3.16)  3D  The first function in this equation is interesting. It was obtained in the  y/z-plane  where the  necessary logarithmic singularities are compatible with the requirement that theflowfunction be single-valued in theflowfield.  General Tandem Thin AirfoiJ Theory  /  3.1  46  The thickness flow function for tandem NACA airfoils is obtained by scaling Equation 3.16 to the airfoil and flap geometries, and then linearly superposing the two resulting functions. The airfoil pressure distribution becomes:  s/s  —  yfi  B+  +  p*  2C --rZD^)\AE(?) ]\» X  Z  S X  + 2C + ZD  7r - 2tan  +  1  If-x  (3.17)  1-/ In  +  1- x ~x  + 2C + 3D  and theflappressure is:  ir |  2V  *  y/x + ^/s  [  W  s  \sJ  J  x  —  s  *1 , n \f \ + — + + 2C + 3D - + -\+AE\[-) s 2] [\3/ 2s 3 x  IA / ! - /  +  .  +  2  x  2  x  1  (3.18)  y ^ T - y ^ T  ^  + 2C + 3D  +  3  1-/  . ( ^ )  +  2 j  +  4  2  i  i  +  4  , ( ^ i )  K l - / y  3  ] l n l ^  2(1-/)  +  3  These pressure distributions are graphed in an example in the next section.  General Tandem Thin Airfoil Theory / 3.2  3.2  47  T H E C A M B E R LINE DESIGN P R O B L E M  Here, pressure distributions which are odd in 0, and 0j are specified and the camber lines needed to create them are then worked out. The general function to work with is, again, a simple linear combination of one-element flow functions: (3.19) The first grouping in this equation is purely imaginary everywhere on the real axis except on the airfoil slit and thus does not affect the flap's pressure distribution. Similarly, the second group of functions has no effect on the airfoil's pressure distribution. Therefore:  P. =  9  £.o V t a n  C  + 2 XX sinn*. n=l oo  C  = Bfo^T  P,  + 2  ZX n=l  (3.20) 8inn  */  Camber lines are usually designed to give optimum performance in a particular situation; i.e., for given values of lift and, perhaps, moment. Since this ideal or design performance is achieved by minimizing adverse pressure gradients, B  j0  is usually chosen  to be zero. Then: B,  n  ='iy  C . An n0.de. p  (3.21) B  sinntfy d0j  and the lifts are simply: Ct=z+£j[»B.i+V-f)B ] fl  C  tt  = -nB  tl  ;  (3.22) C  tj  =  -*B . fl  Of course, each group in Equation 3.19 affects the forces on only one airfoil element since it affects the pressure on only one airfoil element. Again, these simple results and the simple construction of the camber design flow function are also applicable to multi-element in-line airfoils.  General Tandem Thin Airfoil Theory  /  48  3.2  A better way of handling the camber design problem is to use closed form functions since camber pressure distributions are usually specified by simple algebraic functions. In the following example, tandem airfoil camber lines which combine uniform chordwise loading with zero pitching moment about the aerodynamic centre are designed. The camber function: ...  =  iff,. IT  z — s  iffr.  In—  ,  z—1  f i n —  Z  Z  IT  3.23  — J  gives constant airfoil and flap pressures which are odd in 6, and 6j (ff, and ffj are constants): C . = T2ff,  ;  P  C  = T2ff .  pj  (3.24)  f  The airfoil boundary is: dy .  o-,,  c  -TT = "T dX  c  that y ,(s)  1~~  +  X  ffr, 1- x "T 7—Z • M  IT  J  3  (-  25  — X  =0:  c  Cl  s—x  is now obtained by integration. The constant of integration is arbitrarily chosen so  y ,(x)  y (x)  IT  m —  = ^f- [s In s — (s — x) ln(s — x) — x In x]  + — [ ( / - *) ln(/ - *) - (1 - x) ln(l - x) + {1-3) l n ( l - « ) - (/-a) l n ( / - a ) l . IT  (3.26) Similarly, y Ax)  is obtained by integration in such a way that y Af)  e  c  =0:  y ,(x) = — [ ( z - s ) l n ( x - s ) - x l n x - ( / - s ) l n ( / - s ) + / l n / j c  *  IT  + % ! - / ) l n ( l - / ) - (1 - x)ln(l - x) - (x - /)ln(x - /)].  (3.27)  It is interesting that:  =  ffj  Mil ff,  in(l-a) + / l n / - ( / - a ) l n ( / - 3 ) ] .  =  (3.28)  IT  The choice to zero the camber lines at x = 5 and x = / is arbitrary. Some stagger between the airfoil elements is acceptable as long as the linearizing requirement, that deviations from the average line be small, is kept in mind.  General Tandem Thin Airfoil Theory / 3.2  49  Now: C  <=7j^Z7k. + (W)<7]  (3.29)  and: c  - o  7-/)^sV«+^-f >f\2  (3+  =  ( - °) 3  3  If the flap chord to total chord ratio isfixed(so that only the slot size is variable), then it is convenient to define: = sTTJ  S = f-s.  ;  (3.31)  The lift is then independent of /—s: C = 4<x, + 4(ay - cr )R . t  t  (3.32)  f  From Equation 2.49, the aerodynamic centre is:  \ " (i "  ZflC =  5  (3  '  33)  and its associated pitching moment is: Cm„ = ~<r. + (cr/ - O c  If C '  mac  /  L  [2^/(1 + S) - 3 - 5] _ j i. 7 V  l  , (3.34)  can be designed to be zero, then: CTf  _  1-5  T.  =sl  +  R [2R {l + S)-3-S]' f  f  ( 3  '  3 5 )  If 5 is small and the ratio o~f/o~ is required to be positive, then this last equation only has t  solutions for Rf > |. As an example, consider a tandem airfoil configuration with Rf — 0.625 and S = 0.2. A number of current tandem wing light aircraft designs use this type of arrangement. It should be noted, however, that many other considerations are important in the design of an aircraft wing and that this example is chosen mainly to exhibit the way in which camber design can be accomplished with thin airfoil theory. The slot size in this example is probably too small for a tandem wing aircraft and it is certainly too large for a slotted slat arrangement. It is chosen so that K M K'.  General Tandem Thin Airfoil Theory / 3.2  50  If, in this example, a design lift coefficient of 0.6 is chosen and if C  mae  is  required to be zero, then a„ = 0.2833 and <Tj — 0.07. Figure 8a shows the camber lines that satisfy these requirements (Eqs. 3.26 and 3.27). If the NACA thickness distribution of the previous section (Eq. 3.14) is combined with these camber lines (according to Eqs. 1.7), then the airfoil shapes in Figure 8b are obtained. Figure 8a also compares the thin airfoil theory predictions for chordwise loading, Eqs. 3.24, with those which exact-numerical potential flow theory (Appendix A) would predict for the same geometry. The Becond order effects of thickness are shown by including exact-numerical computations for both 1% and 15% thick airfoils. Figure 8b shows the actual pressure distribution predictions over the 15% thick airfoils. The thin airfoil theory predictions of these pressures are obtained from Equations 3.24, 3.17, and 3.18. Table 1 compares the lifts and moments of the various computations for a = 0°. The linearized predictions in this table can be added to the a and n problem solutions of the previous chapter to solve the general force problem for this configuration.  AC  t  AC,,  AQ,  AC  m o  ACL  &C .  m 0  mac  0.6  1.133  0.28  -0.244  -0.567  -0.14  0  Ex-num 1% thick  0.585  1.087  0.284  -0.241  -0.542  -0.138  -0.003  Ex-num 15% thick  0.623  1.154  0.305  -0.253  -0.553  -0.144  0.000  T.A.T  TABLE 1  Potential Sow predictions of the forces on the tandem airfoils of Figure 8.  The overall thin airfoil theory force predictions agree well with the exactnumerical results, with the same trend as was seen in the one-element theory of the exact theory overpredicting the linearized theory for normal 15% thick airfoils. This accuracy is consistent with that of AC in Figure 1. The individual airfoil element forces also show e  good agreement. As in one-element thin airfoil theory, pressure distribution predictions can  General Tandem Thin Airtoil Theory / 3.2  1  •  •  •  i  i  -  51  1"  1  ,  •,  1  ,  • •  —  • •  •  •  Thin airfoil theory.  •  • Exact-numerical, 1% thick.  •  • Exact-numerical, 15% thick.  • 0.8  ••  •  • • 0.4  r  •  —* (LC.  0  .1  1  1  .2  .3  "  *  —*  • *  • •  • 1  .4  1  I  I  i  •  .5  .6  .7  .8  .9  1.0  x  a) The camber lines and chordwise load distributions for 1% and 15% thick airfoils, a = 0°. A C = C - C p  FIGURE 8  p  Potential Bow predictions of tandem airfoils designed by thin airfoil theory to give uniform chordwise load distributions at a = 0°. f—s = 0.2, ( ! - / ) / ( « + ! - / ) = 0.625, K = 1.07K'-, x = 0.325. ac  General Tandem Thin Airfoil Theory /  3.3  52  be fairly inaccurate at localized points but, nevertheless, do show good overall qualitative agreement with the exact potential theory. The leading edge singularity in the thin airfoil theory thickness pressure distributions of Figure 8b is an interesting irregularity. It is caused by a hidden discontinuity in ±dyjdx  through x = 0, the second term in Equation 3.15.  as x —• 0 in Equation 3.17 and since B is negative, C  pi  Since C  pt  ~ O (-Bln^)  —* —oo. This is, perhaps, sur-  prising because thickness functions usually give positive infinite pressure singularities at leading and trailing edges due to the presence of stagnation points. However, the rounded leading edge of an airfoil is like that of an ellipse and, as with an ellipse, the pressure singularity associated with the stagnation point is not present on the airfoil's surface.  3.3  C A M B E R LINE ANALYSIS: G E N E R A L FUNCTIONS  This section presents two different camber functions which can be used for calculating pressures, lifts, and moments generated by arbitrary camber line shapes, including the simple shapes of the a and n problems. The first function is the simplest in that it uses the elementary functions of the previous two sections of this chapter and provides easily understood qualitative information. The coefficients of its Laurent series reduce to the one-element theory coefficients when f—s = l. However, since one set of these coefficients fails to converge when / — s = 0, this function will be of limited practical use for small slot configurations. The second function is also constructed from elementary functions but ones that are of a more complicated nature. It has the advantage that the coefficients of its series converge like the one-element theory coefficients at both slot extremes. Both of these general flow field descriptions are unsatisfactory in that the coefficients of their series must be recalculated each time the slot size changes, even though every other aspect of the geometry may stay the same. This problem is compounded by the fact that all the coefficients are calculated via elliptic integrals, and all of them are needed for calculating pressures or localized forces. The problem is less acute when only  General  Tandem  Thin  Airfoil  Theory  53  3.3  /  overall forces are being considered since, as with the one-element theory camber function, the overall forces from both functions are described with only the first two or three of the coefficients. Even so, the direct methods of the next chapter are more efficient when evaluating either localized or overall forces. Interestingly, these direct methods were first motivated by the qualitative information provided by the camber functions of this section, as presented in the last section of this chapter. If pressure distributions are required and the camber line is specified as an analytical function, then the methods of the next section of this chapter are the most efficient ones to use. If the camber lines are numerically specified, then the two functions of this section are probably as efficient a way of obtaining the pressure distributions as any other. If the slot is small, then the second camber function should be used. The First Camber function is: B  ,  f0 1  where the  and  (3.36) n=l  V  planes are as shown in Figure 6 and as given in Equations 3.1-3.6.  The first group of functions in this equation is purely real on the flap since the square-root expression in it is purely imaginary on the flap. Similarly, the second group of functions is purely real on the airfoil. Thus:  (  /T=L  dVc  L=  dx  y f-x  X  '  dy,  X  L  and:  .  TT  J  dx  V1-  '  2 ^  ,  n=l  dx  IT J  (3.37)  cos nO.  —x  cos n9 V Z ( l - x)(s : — x)  g  dx.  (3.38)  Similarly: dx  c  B  fn /B  s  B fo B  +  COSTl9  J2 fn  (3.39)  f  n=l  and: 2 f dy I ^ = - -r^-J KJ dx\x-s  —  2 t d]ic I ~~x cosn9,dO =: - - 7 A / 7 w H T ; rcosn0 dx. irj dx V (x - s)(x - /)(1 - x) f  1  1  f  J  v  (3.40) '  General Tandem Thin Airfoil Theory / 3.3  54  When f — 3 —• 1, Equation 3.36 reduces to the conventional one-element theory representation of camber for two airfoils infinitely far apart. There is no problem with the coefficients in this geometrical extreme. However, since the Bf 'a blow up when / —• s, n  this representation of camber is of little practical use for small slot configurations. The overall lift is: C  t = Z ^ n B  t  o  + B  t  l  ) + (l-f)(B  + B )]  fo  (3.41)  fl  and the overall moment is:  C  - o  =  4 ( 3 +  I_/)2[3(3-2+2/)^o  +2«(»-l+/)B. + s B, 2  1  2  + (l-/)(3/+l-2«)fl + 2(l-f)(l+f-s)B  fl  /o  +  (l-f) B ]. 2  h  (3.42)  These expressions, through the B  and Bf coefficients, give the effects of the independent  tn  n  airfoil and flap boundary conditions on the overall forces. However, unlike Eqs. 3.22, the separate groupings of B  tn  and Bf coefficients in the expressions do not give the effects n  of the localized airfoil and flap forces. In general, all of the B  ln  and Bj coefficients from n  Equation 3.36 will have an effect on the localized forces. The Second Camber Function's Laurent series are defined in the f-plane which is arrived at through the elementary mappings illustrated in Figure 9. As f-s -* 0, this f-plane reduces to the f-plane of Figure 3 — the one-element theory's solution plane. When f—s —*• 1, the first Riemann surface of the f-plans becomes the f,-plane of Figure 6 , while the second surface becomes-the ( -plane. f  The intermediate z -plane is arrived at through a simple Joukowski trans2  formation with critical points at ZQ and S  Q  in the 2-plane, corresponding to Z2 and ^ Q  Q  in  the 2^-plane. Each side of the circle outlined in dots in the z-plane is transformed to the vertical line joining Z the z - p l 2  a n e a n  2 Q  and 2 2 ; the left side of the circle is in the first Riemann surface of 0  d the right side is in the second Riemann surface. One changes surfaces by  crossing the branch cut joining ^  and Z  2 Q  . The airfoil slit is also in the first Riemann sur-  face and stretches between —2 and 2 . It is superimposed exactly over the flap slit which  General Tandem Thin Airfoil Theory /  2Q  3.3  55  z - PLANE  z -PLANE 2  V2  XO —•—  -2  X2  2 •I  • ••• • •  2  2o  20  -Branch cut joining Z and z g via oo. 2  a)  Linearized  physical  plane.  b) Intermediate  $0  c)  Solution  plane  FIGURE  with  9  the  The  airfoil  in  elementary  the  first  Riemann  mappings  surface  for  the  2  Q  solution  plane.  f -PLANE  and  second  the  flap  camber  in  the  function.  second.  General Tandem Thin Airtoil Theory / 3.3  56  is in the second Riemann surface. -  z  z 2  2  ZQ  =  z  z  s  z  1  ( - f ) + ( ~ )( - ) z(z - f) - (z - s)(z  X  0  + iy = Q  2  - 1)  =  +  («-«)(«-!) Z-XQ  a + ,V«(l-/)(/-«) 3+ 1 - /  1  z — - (l + «-/)«2 + 2(1 +S + f ) -  8  2  Q  (3.43) (3.44)  (l +  S-f)yJ(z -Z )(z -Z ) 2  2Q  2  2Q  [4s - (1 + s + f)(l + s-f)] + tW«(l-/)(/-«)  % = *2 +m = Q  «(»-/)  a+1-/  («+l-/)  a  (3.45)  (3.46)  The f-plane is simply the 2 -plane with the slit opened up to the unit circle: 2  *2  or:  (3.47)  = f+ ~  /) + >/(ar-«)(ar-l)  (3.48)  Since Equation 3.47 has the same effect in both Riemann surfaces, the opened up airfoil andflapslits are also superimposed in the f-plane—but still in different Riemann surfaces. Note that infinity in the z and f planes corresponds to z = oo or z = x 2  Q)  depending on  which side of the branch cut and in which Riemann surface one is located. The second camber function is:  ? +1 n—i  s"  * - *o  (+1  ~  n=l  r  (3.49)  Here, the function t'/(f + 1) solves the tandem airfoil a-problem and is equivalent to Equation 2.32. Since, in general, the disturbance velocity will be non-zero at z = x and 0  since the functions 1/($ + !) and 1/f both have simple zeros here, the function 2yo/(2-Xo) helps ensure the generality of w' (z) without introducing a singularity into theflowfield. c  It also enables the second group of functions of Equation 3.49 to take on different values on the unit circle in the different Riemann surfaces of the f-plane for the same value of 9, whereas the first group has the same value in each Riemann surface for a given 9.  General Tandem  Thin  Airfoil  Theory  3.3  /  57  Thus, using the notation of Figure 9 a : 2yp Z  [z -  _  - z ) - (z- x )(z - x )  ZQ)(Z  Q  a  h  (3.50)  — XQ  Z —  XQ  so that: co  (x  3g(  It)  ^ = ^+fXc°sng+ ~ *~ x-x dx  0  n=l  2  + ^  Z? cos n0 n  n=l  X-Xn  Y + 52 » OO  ( X - Z  0  ) ( X - S  0  )  X-XQ  D  cosnd (3.51)  n=l  This separation is important since the function: (z-x )(z-x ) a  b  1+3-/  =  (2 ~ z  Z-XQ  A  (3.52)  *2 ) n  takes on the same values in each Riemann surface—i.e. it is even about the branch cut. On the other hand, the function: ( Z - Z  1+a-/  ) ( Z - Z Q )  0  (3.53)  •yJ{z2-z ){z -z ) 2o  Z -  2  2o  XQ  is odd about the branch cut. Thus, the first line on the right hand side of Equation 3.51 takes on exactly the same values on both the airfoil andflapfor given 0's, while the second line takes on values that are equal in magnitude but opposite in sign. Hence: d  y, dx c  =  _  2  (x - z )(x - 2 ) 0  Q  -^ + J2D cosn9  (3.54)  n  dx  X —  XQ  n=l  This equation is only true in the first Riemann surface (the RHS would have to change sign in the second surface) so that 0 < x < s. The i?„'s, then, are given by: X - X dx  0  (X-ZQ)(X-ZQ)  o  dy,  cos n9 dO  dx  X-XQ —r  (x-z )(x-z ) 0  cos  n6d9  (3.55)  0  o  where 0 < x < s for the first integral and / < x < 1 for the second integral. Since the second integral is evaluated on the flap, the sign of its integrand must be changed relative to what it is in Equation 3.54.  General Tandem Thin Airfoil Theory  3.3  /  58  Continuing on:  dx  -^ + ]Tfl cosn0 + 2 n  dx  (x-x )(x-x ) a  DQ  6  + ^  n=l  D  N  cos nO  n=l  oo  dy ,  dy  dx  dx  c  = 2 n=l  (x-x )(x-x ) a  C}  (X -  b  Z Q ) ( X  -  ZQ)  so that: dVc =  1 +  dx  n=l  (x-xj(x-x ) 0  >  +  (x-z )(x-5o) 0  <*y,  j _  (x - x )(x a  -  X ) fr  (x-2 )(x-2 ). (3.56)  dx  0  0  with 0 < x < s. Then: 2  M2/c.  o  (*-*o)  5  (ix ( x - Z o ) ( x - 5 )  cos n0 d9 +  2 ? rfy (x - x )2 cosnfldfl T 7 dx ( x - z ) ( x - 2 ) o e/  0  0  0  (3.57)  0  where, again, 0 < x < s for the first integral and / < x < 1 for the second. This requires a change of the sign in the middle of the last square-bracketed expression of Equation 3.56 for the integrand of the second integral. Changing variables, the coefficients can be rewritten as: 2  f dy  (x — x) cosnft 0  Ct  dx > / x ( 3 - x ) ( / - x ) ( i - x )  dx  +  1/ 2  f dy  (x — x )cosn0 o  Cf  dx  yJ (x-  S){X- / ) ( 1 - l )  X  dx (3.58)  cos n9 n  i r j d x  x/xTT^ Vx(s-x)(/-x)(l-x)  dx  1  /"<ty  cosn0  c  dx >/x(x-s)(x-/)(l - x )  dx. (3.59)  The behaviour of these coefficients as f—s —• 0 and f—s —• 1 is best examined for constant ify = (l-/)/(s+l-/). When / - a —• 0, the D 's become 0(ln(/-s)). However, n  yl = /fy(i — Rj}(f— ) go that the second group of functions in Equation 3.49 disappears. S  In addition, the f?„'s reduce to the one-element coefficients, the B 's of Equation 1.25, and, n  since the f-planes of the one-element and two-element mappings are equivalent for f—s = 0, the one-element camber function, Eq. 1.23, and those functions remaining in Equation 3.49 will also be equivalent.  General Tandem Thin Airfoil Theory /  For  3.3  59  t h e i n f i n i t e slot e x t r e m e , w h e r e a —• 0 a n d /  —• 1, a c h a n g e t o t h e  ( z ' , y ' ) v a r i a b l e s i s m a d e s o t h a t , i n t h e first i n t e g r a l s o f E q u a t i o n s 3 . 5 8 a n d 3 . 5 9 , x' =  x/s  a n d y ' = y/s. T h e s e c o n d i n t e g r a l s a r e d e a l t w i t h i n a s i m i l a r w a y , r e s u l t i n g i n :  B - (1 - R )B n  1  n ^ ^  D  where the  f  t  + RB  tn  f  as f-s  fn  1  1  n  - - B  f  (3.60) '  a s / - a - + l  n  B *s a n d i ? / ' s a r e t h e o n e - e l e m e n t t h e o r y c a m b e r f u n c t i o n c o e f f i c i e n t s f o r a n tn  n  i n f i n i t e l y s e p a r a t e d a i r f o i l a n d f l a p , r e s p e c t i v e l y ( E q s . 3.38 a n d 3 . 4 0 w i t h / = 1, a = 0 ) . On  t h e a i r f o i l , a s f-s —-°  -* 1 :  - -2R  0 < x < a  ,  f  (3.61)  Z-XQ  w h i l e o n t h e flap: ^  2(1  -R )  f<x<\  ,  f  (3.62)  Z-XQ  so t h a t o n t h e a i r f o i l :  —iB  B  00  <{*)-—r7-*E7? n=l  A  >  0 < z < a  (3.63)  .  / < * < ! •  (3-64)  a n d o n t h e flap: u  —iB  B  00  'cW-7Tf- 'E-X +  1  i  „ = i s/  T h u s , t h e s e r i e s o f t h e s e c o n d c a m b e r f u n c t i o n c o n v e r g e at b o t h f—s e x t r e m e s i n t h e s a m e way t h a t t h e one-element c a m b e r f u n c t i o n does. Therefore, there s h o u l d also b e reasonable convergence b e h a v i o u r f r o m these series f o r i n t e r m e d i a t e slot sizes. The C C  t  mn m  °  second c a m b e r f u n c t i o n ' s overall forces are easily o b t a i n e d b y residues:  = - 7 r ( 2 ? + B{)  (3.65)  0  =  . * _ 4(a+l-/) f  )  2  2  { [1 - ( / - a ) ( 3 / + a - 2 ) ] Z ? 2  + (l+ -/) 2? 5  2  -  0  2  +  2  2  2(l+a -/ )B,  8y (l+*-/)(A> +  £>!)}  (3.66)  or: (3.67)  s+l-J The  e x p r e s s i o n s f o r Ct a n d C  mae  camber function, with C  t  a r e m u c h s i m p l e r t h a n t h o s e f o r t h e first  c o m p a r i n g f a v o u r a b l y w i t h E q u a t i o n 1.26 i n s i m p l i c i t y .  General Tandem Thin Airfoil Theory /  3.4  60  The methods of the previous chapter can be used to obtain the localized forces from the first and second camber functions. However, all the coefficients in the series of these functions will contribute to the forces and, in order to "harvest" their contributions, one needs a recurrence relation to obtain all the terms in the expansion of the airfoil selector A[w) Eq. 2.39. This is provided in Appendix C. t  General functions are not often used for actual calculations and no numerical examples have been worked out for either the first or second camber function. Their practical worths have yet to be proven. Qualitatively, however, they are important as relatively simple general descriptions of what has historically been considered a complicated flow problem. Note that the extension of these camber functions to the problem of an arbitrary number of in-line airfoils is straightforward, although the algebra for the second type camber function could get messy.  3.4  C A M B E R LINE ANALYSIS: SPECIAL FUNCTIONS  When camber lines are specified as analytical equations, the procedures of this section can be used to solve for pressure distributions, lifts, and moments. If only lifts and moments are required, the procedures of the next chapter are more direct. The special functions presented here solve for w' (z) when dy /dx is given as c  c  a simple algebraic function of x. Thefirstfunction is simplest in that it solves a generalized o-problem using only elementary functions. It was first found by Richardson (1981) who used it in describing the aerodynamics of perforated bridge decks. The function takes dy /dx to be an n c  t h  order polynomial which is continuous over both the airfoil and flap.  Its coefficients of x are highly dependent on variations of f—s. n  The second function assumes dy Jdx and dy /dx each have their own indec  pendent n  t h  Cf  order polynomial representation. As f—s varies, the coefficients of x" in these  functions also vary, but in a simple manner. For example, if the airfoil's boundary condition  Genera/ Tandem Thin Airfoil Theory / 3.4  61  were specified as a polynomial in (x/s) , its coefficients would not vary at all with f—s and n  the coefficients multiplying its x terms would only vary as l/s . Note that the coefficients n  n  for these polynomials can be obtained if the one-element theory's Fourier series representation of the airfoil element's camber line is known, Eqs. 1.24 and 1.14. In this case, the order of the polynomial will equal the order of the truncated Fourier series. The problem with the second function as a general representation of tandem airfoil camber is that the coefficients of the polynomials are 0(l/a ) and 0(l/(l—/)"). Thus, if f—s is large or if the airfoil or n  flap is relatively small, special precautions must be taken to handle the very large high order coefficients, which probably have the least effect on the aerodynamic characteristics. Now, for the first function, if: ^ = a + a x + a x + a x*+ --- + a x ax 0  l  2  2  3  n  n  0 < z < s,  ;  /<z<l  (3.68)  then the following function is required:  where the 5,'s are chosen so that A (;z) —• 0 as z -* oo. That is: m  A  m  ( z ) ~ ^ + % - + % - + --Z JZ z°  ,  asz-oo.  (3.70)  To find the A . constants, the expansion for F(z) is needed: m  F(z) =  >-«)(*-1)1 i z{z - f)  F F 1 + -7+ z z + •••  ,  as z -> oo.  (3.71)  A recurrence relation for the F,'s is given in Appendix C. Then, the 5,'s are obtained via: 6i + Fi = 0 6 + 6^ + F = 0 2  2  6 + 6 F + 6 F + F =0 3  2  S + 6 _ F +6 _ F m  m  l  l  m  2  2  X  {  2  (3-72)  3  + --- + F  m  =0  and the A . ' s are: m  A  m i  = - [6 Fi + 6 ^F m  m  i+l  + ••• + «  +  m  _ , + F ]. i+m  (3.73)  General Tandem Thin Airfoil Theory  / 3.4  62  Thus, A ( z ) has its real part equal to x on the airfoil and flap, it goes to zero as z ~* co, m  m  and it satisfies the Kutta condition. The flow function for the boundary condition of Equation 3.68 is: = » X)  a  mA (2).  (3.74)  m  m=0  The forces are:  m=0  4TT ^f"  0  = ^  a+  1  "  _ ^ 2 YI  a  m^m  •  a  (3.76)  Richardson ( 1 9 8 1 ) found this flow function, Eq. 3 . 7 4 , for a parabolic camber line and noted that its extension to include higher order terms was straightforward, as it is. He also extended the case of tandem airfoil parabolic camber to N in-line airfoils in a parabolic camber configuration. This latter case can be further extended to include camber lines with higher order terms by using the square-bracketed expression of Equation 2.51 to replace F(z) in Equation 3 . 6 9 . The nature of Richardson's perforated bridge deck problem is such that Equation 3.68 is a practical description of the camber line. For multi-element airfoils, however, it is important to be able to use independent descriptions for each airfoil element's camber line. Thus, in the tandem airfoil case: % i = a + a x + a x -\ dx 0  dyc  -L  dx  =  {  2  2  +ax  1  2  2  ;  0<z<s (3.77)  d + d x + d x + --- + d z 0  n  n  n  ;  n  /<x<l  for which the following function is required  A,„M - « " { T If §M " h T ? '• • S } ' +  +  +  +  +  FW  which can be used to solve the flap boundary condition if: A  /  J z ) ~ ^  +  ^ + ^p- + ---  ,  as*->oo.  (3.79)  (3  General Tandem Thin Airfoil Theory / 3.4  63  The recurrence relation for the ^'s in the expansion: J(a;) = — + —  E + -^-+ - § + • . .  ,  3  asz-oo  (3.80)  is given in Appendix C, so that:  6  and the A r  h+ fJi+6fA  + V3F = T3  S  (3.81)  3  's are: A/  m <  =  Tm+i -  [8 Fi + Sj^F^ fm  + ••• + fyji^+m-i + E F , 3  + m  ].  (3.S2)  For the airfoil boundary condition, the required function is: A. W = A ( z ) - A m  m  and so:  / m  (z)  (3.83)  n  = ~ J2 [«m(Am(^) - A (z)) + d A Jz)\ l  /m  m  (3.84)  f  m=0  solves the camber line problem of Equation 3.77. The forces are:  4JT C  N  fl  A  A  ' = ^pE[ «( »..- / ) + A ] 4JT " X [ m(Am ~ A ) + d A J . ~~ (a+1-/) ^ ( i  mi  1  (3-85)  m=0  C m o  ^ — 2  v  7  '  a  2  n  3  / m a  m  /m  (3.86)  m=0  The pressures are obtained from the flow functions. For the case of independent airfoil and flap camber an elliptic function must be evaluated each time the pressure is evaluated. The procedures of Appendix B greatly facilitate this kind of calculation but the series presented in Whittaker and Watson (1927) could also be used. Other special functions dealing with other analytically specified boundary conditions could also be found. For example, the four and five digit series NACA cam-  Genera/ Tandem Thin Airfoil Theory /  3.5  64  ber lines are specified with two polynomials in x, one applying over the forward segment of the airfoil element and the other over the back segment. This discretization of the boundary condition can also be handled with analytical functions but it gets so complicated that the general Laurent expansions of the previous section may be more efficient. If an analytical description of the flowfieldsaround NACA type airfoil elements were to be attempted, then the airfoil and flap segment selectors of the next chapter would be needed. These functions use logarithmic singularities of the type used in the one-element theory's ^-problem, Eq. 1.19, to obtain functions with imaginary parts zero everywhere on the airfoil and flap except over the desired segment of interest, where they are constant. Analytical expressions for the lifts and moments generated by the NACA camber lines are obtained in the next chapter using these functions.  3.5  THICKNESS DISTRIBUTION DESIGN AND OTHER PROBLEMS  In the thickness design problem, desired pressure distributions are specified and the associated thickness distributions then worked out. The functions that allow this are the same as the camber functions of the previous two sections, except that they are multiplied by ». This allows the real part of the functions on the airfoil and flap to be specified as the boundary condition. The functions of this chapter can be combined in other ways that allow the designer the versatility of solving mixed problems. For example, suppose an airfoil's shape is given and the designer wants a flap with a specified pressure distribution. The camber problem for this situation could be solved with theflowfunction: . (3.87) which uses the  and ( planes of Figure 6. The first group of functions here provides a f  general camber line boundary condition for the airfoil and is purely imaginary on the flap. The real part of the second group of functions provides a general pressure distribution boundary condition for the flap and is purely real on the airfoil.  General Tandem Thin Air/oil Theory /  3.6  3.6  65  MUNK'S INTEGRALS EXTENDED  In the same way that Munk's integrals were obtained from the general one-element theory camber function, Eqs. 1.33 and 1.34, one can use the general camber functions of section 3.3 to solve directly for the tandem airfoil forces in terms of the boundary conditions. The general camber functions' overall lifts and moments are given in terms of the first two or three coefficients from their Laurent expansions, Eqs. 3.41 it 3.42 and Eqs. 3.65 & 3.66. If, in these equations, the coefficients are replaced with their respective integral expressions, Eqs. 3.38 & 3.40 and Eqs. 3.58 U 3.59, then, after some simplification, one gets:  C,=  -4 s+1-/  Ci  'mo  "  f  2  +  4 (s+1-/)  (3  «  2  /  x(/-x) dx — x)(l - x)  dy . c  dx  I  +  fdy^  I  x(x-f)  J dx V (*-«)(!- dx  x{f-x)  V(a-x)(l-x)  dx +  J dx  (x — s)(l — x)  (3.88)  dx (3.89)  By letting / —* 3, these expressions reduce in a very simple manner to Munk's original integrals. Obviously, if only lift and moment are desired, there is little reason to use the procedures of sections 3.3 or 3.4. Moreover, the simplicity of these expressions suggests that there must be some simpler way of deriving them — as, indeed, there is. This is presented in the next chapter along with efficient methods of evaluating the integrals.  CHAPTER 4 Tandem Airfoil Forces an Easy Way  At the end of the previous chapter, general integrals which give the overall lift and moment for any tandem airfoil camber line were presented. In this chapter, these integrals are rederived using a simpler, more direct method and the analogous integrals for localized lifts and moments are also obtained. Although all of these force integrals can be solved numerically, they can also be solved analytically by using selector functions. As an example, analytical expressions are found for the forces on tandem NACA airfoils. These thin airfoil theory predictions are then compared with exact-numerical potential flow theory predictions of the forces.  4.1  T H E OVERALL FORCES  The general expression for overall lift is: (4.1)  C^j^-jfwMdz  c where w' (z) is a general camber function that solves the a and v problems as well as c  problems involving conventional camber line shapes. The closed curve C in this integral can collapse onto the surfaces of the tandem slits since w' (z) has no simple poles in the c  z-plane:  66  Tandem Airfoil  Forces  an Easy Way- / 4.1  67  However, on both the airfoil and flap, dy /dx is even in 7, Fig. 4d, and dz is odd in 7 so that e  the imaginary part of this integral is zero (i.e. the drag is zero, as proved in section 1.3). On the other hand, C dz is even in 7 which not surprisingly leads to the result: p  l  a-  2  "  C„ dx  s+1-/  (4.3)  Thus, the problem with Equation 4.1 is that when C collapses onto the slit boundaries, the object of the integration becomes the unknown pressure distributions while the known boundary conditions are lost. This situation can be reversed in the following way. w' (z) —• 0 as z —• 00: c  ^ ( 2 ) " - ^ - + -!+ - f + ' • • z z  .  asz-00  (4.4)  so that: w' (z) dz = 27rtiy . c  (4.5)  1  c Now suppose there is an C(z) such that: j  w'c(z)dz  = j  C  .  w'c(z)C(z)dz  (4.6)  C  Then it must have the form: £(z)~l + ^  + ^  + ^  + ...  ,  asz-00  (4.7)  in order for the residue of the new integrand to be equal to wj, and since the u^'s are arbitrary while the £ ' s will be fixed (zn n  terms, n > 0, are possible in this expansion  of £(z) but, since the residue must remain equal to wlt the coefficients of zn would depend on the unknown tu„'s). Further, suppose £,(z) is purely imaginary on the slit surfaces and that  w'c(z)C(z)  has no simple poles in the z-plane. Then, the new integrand will collapse  in a simple manner to the surfaces of the slits where the imaginary £(z) will combine with the imaginary part of w' (z) c  to form the real part of the integrand, which will give the lift,  while the imaginary part of the integrand containing the unknown pressure terms must integrate to zero.  Tandem Airfoil Forces an Easy Way / 4.1  68  The conditions on £ ( z ) are very similar to those on the flow functions of the previous chapters: 1) the real part of C(z) is zero on the slit surfaces. 2)  C{z)  -*  1 as z  -*  oo  3) j£(z) must have its singularities at the trailing edges since the Kutta condition ensures that  w' (z) c  will not have any there. This prevents the singularities of  from possibly combining with those of  L(z)  to form simple poles in the z-plane.  w' (z) c  Since the conditions on C[z) are as strong as those on the flow functions, it will be unique:  m =  (4.8)  and: r  .  I  z  (z  -  v  f) J  .  (z-s)(z-l)  ,  (4.9)  J  >  dz.  Although it is not a requirement, it is not surprising that Z(z) is odd in 7 on the slit surfaces as this ensures that the imaginary part of Equation 4.9 is zero. On the slit surfaces, then, Equation 4.9 reduces to: -4 .+1-/  *</->  (s-x)(l-x) 0  dx + j ^ .  '  « < « - / )  .dx  (4.10)  /  which is identical to Equation 3.88. For the overall moment, let: j  zw' {z) dz = j c  C  zw' (z)C'(z) c  (4.11)  dz .  C  In order for the result of these integrals to be 2Kiw , 2  c  c  c  t'(z) must be of the form: as  z  —• 00.  (4.12)  Also, if L'{z) is linear in £ ( z ) : L\z)  =  g{z)C{z)  (4.13)  Tandem Airfoil Forces an Easy Way / 4.1  69  so that the boundary condition is selected when C collapses to the slit surfaces, then g(z) must be purely real on the slit surfaces. In order for £ \ to be zero: g { z ) ~ l -  S  +  * ~  f  +  ,  o(j£j  as 2 - o o .  (4.14)  Since any higher order terms in this expansion for g(z) create singularities which will either create unwanted simple poles in the z-plane or change the effect of £(z) on the integrand, it must be truncated after the second term. Thus:  (4.15)  C  t  ~T  2  / ,,,,,, kW  /  »(«-/)  • b  (7+wFJ" *V(.-.)(»-i)'  and on the slit surfaces:  _ Ct m o  ~ 2  +  4 (s+1-/)  Mft. 2  /  x ( / - x ~  ; rfx ]J ( - x ) ( l - x ) X  S  a X  +  J  dx  I  / x(x-/) . Y ( x - 3 ) ( l -- x )  dx  (4.16) which is identical to Equation 3.89. The imaginary part of the integrand again integrates to zero since it is odd in 7. Note that for in-line airfoil elements, the imaginary part of Equation 1.12 is also zero since w'(z) can be replaced with  VJ' (Z) C  and (dy /dx)dz c  is odd  across the slit surfaces. These integrals for lift and moment in terms of the boundary condition, Eqs. 4.10 and 4.16, can now be evaluated numerically using appropriate methods to take account of the integrable singularities at the trailing edges of the airfoil elements. If the airfoil elements' boundaries are specified numerically, then this is probably the easiest way to evaluate the camber line contributions to the forces, A C , and A C  m o  . However, if  the boundary condition is specified analytically, then it may be simplest to carry out the integration analytically, as done in the remainder of this section. If dyjdx  is given as an analytic function or as a group of analytic functions,  then it is convenient to find a complex function, -y^(z), which is single-valued in the flow dx  Tandem Airfoil Forces an Easy Way / 4.1  70  field and which: 1) on the slit surfaces has its real part even in 7 and equal to dy /dx and has its c  imaginary part (if it is nonzero) odd in 7. 2) does not create any simple poles in the z-plane when combined with £(z) in the product ^ ( z ) £ ( z ) . dx Then: C t  1  = -J L ; / s+l-fj  C  K  "° - f -  (4.17)  \: „dz 'y (z-s)(z-l)  2i  t  c  c  fo.)./  dx  tj+w?t  Z i  f )  fdy .  z(z - /)  I  Cf  v  J  («-.)(.-1)iz •  {z)2  li V  1  t4  -18)  C  Because of (1), these expressions are purely real and, because of (2), the curve C can collapse unhindered to the slit surfaces where the expressions are equivalent to Equations 4.9 and 4.15. Note that the conditions on -y^(z), while similar to those on u>c(z), are not d as strong. These conditions do not define a unique -r^(z) since there is no requirement dx ^ for it to meet as z —• 00. Of course, one obvious representation for -r^-(z) is to simply  ax set it equal to iw' (z), but this is a rather trivial use of the added flexibility the function c  provides. Consider the following examples. Suppose: ^Y± = a + a x + a x + a x + •••; dx 0  then:  l  2  2  3  0 < x < s, / < x < 1  3  ^(z) = a + a z + a z + a z + • • •. 0  dx  x  2  2  3  3  (4.19)  (4.20)  This allows the simple evaluation of Equations 4.17 and 4.18 by residues. A recurrence relation for the £„'s of Equation 4.7 is given in Appendix C, so that: C = ~**_ t  s  f  [apili + a, r + a l + • • •] 2  2  3  (4.21)  Tandem Airfoil Forces an Easy Way / 4.1  71  When o = - a , these last two equations solve the a-problem. Although they are similar 0  to the special function lift and moment results, Eqs. 3.75 and 3.76, they are simpler to evaluate. Now suppose the airfoil and flap have their own independent camber line representations: ^  = a + ax + ax Q  dy  ct  —  T-?-  dx  x  2  d + diX + d x 0  + •• •  2  ,  2  -\  ,  0 < x < s  (4.23)  / < x < 1.  ,  Then: ^(z)=  (ao + aiZ + a^ + 2  •••)[!-?(u)]  + {d + d z + d z + Q  where [l—7{UJ)]  i  2  (4.24)  ---)T{u)  2  and 7(u) are airfoil andflapselectors with a square-root pole at the flap  leading edge in the z-plane so they will not combine with £ ( z ) to form any simple poles: J(w) =  h — -^(w) = 0 + 1 odd(7) * = l + i odd(7)  *  K  Also:  2  7  7(UJ)  on the airfoil  7  7  ~ 7 +— ~ + --z + -^ z +z 2  Q  3  on the flap.  ,  as z  oo.  (4.25)  (4.26)  Note that 7(u) could have been constructed from the i?4.(w) function which has its zero at the airfoil leading edge. The process of evaluating overall lift and moment using the residue theorem is now straightforward. A recurrence relation for theflapselector 7 '8 is given in Appendix C. N  The lift is: C< = - ^ j ^ { « o [ - * + ^ ( l - ^ + d [7 + £J ]+d [7 0  i  0  l  +flJ  2  l  }  + C 7 ] + --2  0  (4.27)  which reduces to Equation 4.21 when the a 's equal the d 's. When d = —r), Equation 4.27 n  n  0  Tandem Airfoil Forces an Easy Way / 4.2  72  solves the rj-problem. A more general version of this equation and the results for overall moment are given in section 4.3, in Equations 4.60 and 4.61.  4.2  T H E LOCALIZED FORCES  The airfoil lift is: C  tt  =  -- jw' (z)dz. s  c  C Cauchy's theorem enabled the derivation of Equation 2.41 which is in terms of the airfoil selector A(u) (Eq. 2.38). It is convenient to replace F(z) in this equation with w' (z) so c  that: w' (z)A{u)dz  j  = j  c  C  w' (z)dz c  + j  C.  w' [z)[A{u)-l]dz  + j' vj' (z)A{u)dz.  c  c  C,  C  (4.28)  f  However, as z —• oo: ^(w)~l-E and  is O(lfz)  VJ' (Z) C  4  +  o  Q  so that, using the residue theorem: j  w' (z)A(u) d z = ( l - E ) j c  4  w' (z) dz c  (4.29)  The airfoil lift, then, is: C  - £ ) Q + l j w' {z) [A{u)-l] 4  tt  c  C,  dz+- j  w' {z)A{u>) dz . c  s  (4.30)  Cj  If these integrals are evaluated on the surfaces of the slits they enclose, where [A(u>) — l] and A(OJ) are purely imaginary and odd in 7, then the boundary conditions are selected for the real parts of the integrations while the imaginary parts containing the unknown pressure terms integrate to zero: s+l-f,  0 „  ^  .  ^ ITT  0<.1<  y  /  (4.31)  Tandem Airfoil Forces an Easy Way /  4.2  73  Here, the choice has been made to integrate over the upper surfaces of the slits. This expression can be evaluated in this form if necessary but it is simpler to integrate by parts. The first integral is:  / % W 'Tf)- ] * L  where:  ,  1  o  d  =  »e.(*)W»l)-l] ~ / yci')^{il)-l]dx  (4.32)  0  0  * y (x) = j-^dx.  (4.33)  Ct  Since the simple zero in y ,(x) at z = s overpowers the square-root pole in A(v) at z = s c  and since [A{u)-l]  l  at z = 0, the first term on the RHS of Equation 4.32 is zero.  2  \sO{z l )  Also, using the relations from section C.2 (with G and Gj defined in Equations 2.37): 2  G -G  <  1  H  d z '  i  2  ~  z  l  s  (~)  '/(*-*)  2 y/ (z - )( - f)( - 1) z  z  s  ( 4  .  3 4 )  z  where, on the slit surfaces: y/z(z - s)(z - f)[z - 1) = ±i|x(x - s)(x - /)(z - 1)|* .  (4.35)  The following diagram shows which sign to use on which surface of which airfoil element:  0  For the second integral, it is convenient that A(u>) has square-root zeros at both z = / and z — 1 since this then allows y (x) to also be defined with a simple zero Cj  at z = s: dy  c  (4.36)  Later, this consistency of definition (y {s) = 0, j = s,f) will facilitate the derivation Cj  of complex versions of the following two equations (note the singularities at x = s in the  Tandem Airfoil Forces an Easy Way / 4.2  74  —A(u) terms of all the integrands):  3J  S  i  Replacing tu£(z) with  c  2  ~* -f ^ )  .*/  y  -  l x  — X)(f  \/XlS  X](l-  —  x  G -G  1  3  (~)  in Equation 4 . 2 8 gives:  zw'c(z)  [G - 0,1=1] c  (i - * ) + ^  2  Vx( -x)(/-x)(i-x) a  d I +  5  2  ;  £  V«(*-*)(«-/)(i-«)  c / C  (4.38)  where:  *  :  .  .(x) = y x ^ L  j =  ,  S >  /.  (4.39)  Equations 4.37 and 4 . 3 8 are the localized force integrals analogous to the extended Munk integrals for the overall forces, Eqs. 3.88 and 3.89. They reduce to give the localized forces on a one-element airfoil as follows: G i —• 0 ,  2  1  G -* - V / ( l - / ) , E -> - c o s 2  4  _ 1  (2/-l) =  0  -i-  ;  as f-s  —0  (4.40)  where 0j is as shown in Figure 3b. Therefore, as f—& —• 0 :  ^^V-TjJT  + ^ — J + Jjd  —  J — , - , ^ .  (4.42)  These last two equations were also derived independently of Equations 4.37 and 4.38 using only elementary functions.  Tandem Airfoil Forces an Easy Way /  4.2  75  Now, if dyjdx is given by Equations 4.23 (independent airfoil andflapcamber), then it is convenient to construct the complex functions: 2  y (z) =  s)  (a (z -  c  0  +a  z -s  2  + •  2  x  i  2  z — s' + (d {z  -s)+d  0  2  z -s (oo  Y (z)= c  +  (4.43)  2  2 2  (do  -+  l  + « i — z — 4  z - s' Z + di  * - 8  ••)[i-n»)]  + 3  (4.44)  +  Then: z  (  't.  l  s  G —G (~)  ^(i-s^+i(/»,(.) . . ' W 8  3 ^  .  1V  >/z(z-s)(;Z-/)(*-!)  — 2ir» (residue of integrand at a = /)  (4.45)  and:  z(l-s) .  — 2;r»(residue at z = /)  ^  (4.46)  which can be evaluated by residues. Note that since 7{u) i s O f t z - s ) / ) as z —* 8, the 1  ?(U)-J-A(<JJ)  dz  2  product in these integrals has a term 0((z — s ) ) , a simple pole which is -1  eliminated by the zeros in the y (x) and Y .(x) functions. However, the 7(u)/\/z — f Cj  c  quotient does produce a simple pole which contributes a residue to the contour integrals that must be subtracted. In the next section, a more complicated analytically specified boundary condition problem is solved. The results, Eqs. 4.69 and 4.70, reduce in a simple manner to give the results for the above expressions for Ci and C t  for the polynomial coefficients.  mtQ  by choosing appropriate values  Tandem Airfoil Forces an Easy Way / 4.3  4.3  76  S E G M E N T E D BOUNDARIES: T H E SEGMENT SELECTORS  Suppose both the airfoil and flap are composed of two segments, each of which has its own independent boundary condition: dy,Ct d  *  dy i — = a +ax + o i + ax a  0  )*pL =  b  Q+  b  x  l  X +  0 < x< x  2  b2X  2  +  .  X\  dx 2  x  c  f < x <x  2  -p  dx  < x< s (4.47)  dyd = do + d x + d x + dx 1 I dy, o = eo + eix + e2* + ax  dy  x  2  x <x< 1 2  so that 0 < Xy < s and / < x < 1. Then, in order to construct the associated complex 2  function — ( 2 ) , both airfoil and flap segment selectors are required. dx The airfoil segment selector is based on the one-element theory's r/-problem solution,  w' (z) n  (Eq. 1.19). Using the f-plane of Figure 4f and if ft =  e  , 9 i  andftcorrespond  to z = Xj, then:  2ir  V| I = 1. f  f  -l-HoddW  (4.48)  Xi < X < s  for  This function has an unwanted effect on the flap, |f| = R, which is negated by using Milne-Thomson's circle theorem to add a series of image singularities. This leads to the elliptic function: a  2TT  a  2  fi(f-ft) (y-fi)\^-ft) (y-fi) ---  i  2  a  ft(<-fi) (f-ft) (*f-fO (f-ft) •  (4.49)  . c? -ft)n [ -  T  1 1 1  (ft)  Rin  R2n  &+7)+*]  (; - a) n [^*" - - R ' (^+7). ] 2  +i  Tandem Airfoil Forces an Easy Way / 4.3  77  Transforming to the w'-plane, Fig. 4e, using f = e  2,w  ' , R = l/q', and ft = e~ * «, this can 2  x  be rewritten as:  i *  ^  l n s i  («0)  -=Vfr Ki x  n(.' yI)n[»-^"-2(.' x' ) ^] +  +  1  +  *  ' " l (  + x l l r t  which is periodic in ir in the w'-plane as required. Whittaker and Watson (1927, 21.3) gives the infinite products representation for t?j(w). The real part of Equation 4.50 still has the characteristics of the real part of Equation 4.48 on the airfoil. On the flap, however, its real part, while constant, is not equal to zero as required. This can be rectified by multiplying the flap selector,  by the appropriate constant and then subtracting it  from Equation 4.50. As well, it is convenient to convert to the w-plane, Fig. 4d, so that the airfoil segment selector is:  * "r *J(=+^)--r^MM  t a  (4  -  51)  If 0 < t'7x < ^ , then: S(u) = 0 + todd(7)  ,  = 1 + iodd(7)  ,  = 0 + iodd(7)  ,  0<x<zi z, < x < s  (4.52)  /<z<l.  The expansion coefficients for S (u/) are derived in Appendix C:  cc  S(w) ~ $o + — ' + - £ + •• • , z z*  asz->oo.  (4.53)  1?'  Note that, since the -—-(UJ) function in Equation 4.51 is purely imaginary on both the airfoil and flap, the way it complements the logarithm function must be different from the way 7(u) was used to complement the logarithm function of Equation 4.50. When the latter combination was converted to Equation 4.51, terms linear in u/ from each of Equation 4.50 and /(w) cancelled. As a result, the logarithm function of Equation 4.51 has ail the required boundary condition characteristics on the airfoil and flap (Eq. 4.52)  Tandem Aittoil Forces an Easy Way /  4.3  78  0'  but not the required periodicity in JIT in the w-plane. This is corrected by the — (w) 1*2  function so that 5(w) is periodic in jrr. The flap segment selector can now be obtained directly from  by replac-  ing w-i'7! with u - (f-»72) * w+*7i with w - (f+»7 ). This preserves the relative orientation of the vectors w—1'71 and w-M'71 in the logarithm function. As decribed above, a n (  2  0'  the constant in front of -^(w) is chosen to ensure that T(u) is periodic in JTT: #7 7 »  = -ln  2  +-22-i(c).  1 2  t/2(W — I 7 2 J  JT  (4.54)  JT t/2  K 0 < »7 < ^ , then: 2  The  T(w) = 0 + »odd(7)  ,  0<z<s  = 1 + *odd(7)  ,  f < x < x  = 0 + iodd(7)  ,  x <x<l.  (4.55)  2  2  Tn's are derived in Appendix C:  T  T  7(w) ~ 7 J + — + - f + --•  ,  asz-»co.  (4.56)  5(w) and T(w) are similar to elliptic integrals of the third kind (Whittaker and Watson, 1927, 22.74). The -^(w) function, with its square-root pole at z — f, v  2  was chosen for the selector functions since they will be used to construct a -j^-(z) function dx  which should not combine with £ ( z ) , in Eqs. 4.17 and 4.18, to form simple poles at z = s or z = 1. S(<J) and T(w), then, will satisfy the Kutta condition requirement onflowfunctions and so can be used to construct A ( z ) functions, section 3.4, which could provide m  a closed form analytical description of the flow field for the segmented boundary problem and  permit the calculation of pressure distributions therefrom. Now,  the overall lift and moment for the segmented boundary problem are  given by Equations 4.17 and 4.18 with: ^(z) ax  =  {a + a Q  2  l Z  + az+  •••)[!-T{OJ)-S{U)]  2  + (6 + 6 z + 6 z + •••)$(") + (d + d z + d z 0  1  + (e + e z + e z Q  2  2  l  2  Q  2  + ---)[?(u)-T{u)]  x  .  2  2  +  ---)T(u)  (4.57)  Tandem Airfoil Forces an Easy Way  /  79  4.3  In order to solve for the lift and moment using the residue theorem, the following expansions are needed, all as z —• oo: C7(z)  = C(z)7(u>) ~ t 7 + ^ z  0  = C(z)S(u>) ~ £S + — + ^ z z  JLT(Z)  = £ ( * ) 7 » ~ £T + Q z  + •••  0  +^ z  0  £Mi = £ Mi  0  + •••  1,  CS(z)  where: CM = £ Mo\  + ^ z  0  + CiMo;  0  CM = £ M 2  0  +••• c  +11^  2  (4.58)  +£ M etc. 2  (4.59)  Q  Then: —4TT f  ^  =  7HZ/I  a  o ( A - ^ i - W W i W * 2 - £ s ) + --2  + b CSi + b CS Q  1  + e {£7 -llT ) Q  l  + -- + d £Ti  2  + & LT + • • •  Q  X  + e (L7 -LT )  l  l  2  2  + ---  2  }  (4.60)  and:  4JT (  C c  m  =  0  y +  (a+1  _  /)2  {  +  M£2-£?2-£S ) 2  +  6 £$2 + 0  a (£ -C73-£S3) l  2  2  •••  6 £ 5 + • • • + rf £T + d ^ T a + • • • X  3  0  + e (£7 -£T )+e {C7 -£T ) 0  +  3  l  3  2  + ---  3  }.  (4.61)  As usual, the localized forces are more difficult. The segmented boundary conditions, Eqs. 4.47, can be used directly in the expression for C( Eq. 4.31, but not l}  in the integrated by parts version of that expression, Eq. 4.37, since y (x) is no longer c  continuous over an airfoil element. If, for j = a, 6, d, and e: X X { x ) = J^dx and Yj(x) = f ^ d x * • then the first integral of Equation 4.31 becomes:  (4.62)  V j  0  with: x  *i  0  *i  G  G  X  '  Tandem Airfoil Forces an Easy Way /  4.3  80  at z — 0. Hence, the segmented boundary condition version of  since Equation 4.37 is: 3+1-/  (l-^)C  + j{[y {x )-y (x )][A{uJ )-l]  e  a  l  b  l  +  l  G -G  Go — G i  2  dx + V*(«-*)(/-*)(i-«)" *y  +  1  0  v^-x)(/-x)(i-x)  w  dx  * i  ^ xjl-s) G - G, *  J  n  2  sy  j »  7  G  da:  - - / &(*) ^ —  where  /(*-•)  1  , * > ~T  u  */x(x-s)(x-/)(l-x)  /  y e (  3 ^  x  2  ~ ~ x(l-3) - G , A - ? — f r  )  x-/)(l-x)  T/X(X-S)(:I  dx (4.63)  and w are on the upper surfaces of the airfoil and flap. However, if: 2  y {z) 5 y . W [ l - 7 ( u ) - 5 ( « ) ] + y 6 ( « ) 5 ( w ) + y ,(*)T(w) + y.W[7(w)-T(w)] <  c  (4.64)  then the integrals of Equation 4.63 become: t -< 3  y  Go — G i  y (z) —7= c  z(l-s)  w  / ^  w  ^  = ^  2TTI (residue  —  of integrand at z = /) •. (4.65)  Now, if:  G  i  "  G  l  / ( . - . )  as 2 —> oo  (4.66)  then: 1-3  (4.67) G.s  ,  n  n> 1  and it is convenient to define the »7 's and, for z —* oo, the Jmn's: m  4(z)=  £(*)$(*) ~  ?  J (z) = LT{z)Q(z) ~ 1 2  / (z) = LZS{z)9(z) ~ 1 3  Z  Z  Z ^  (4.68)  7 i ^ 4. ^2 , ^30 + — + - y +  z  z*  7 J . ^ 1 X ^ ± Mz) = rrwfl*) ~ 1 JiO H 1 T +  Tandem Airfoil Forces an Easy Way / 4.3  81  The recurrence relation for the J 's in terms of the previously defined expansion coeffimn  cients is similar to Equations 4.59. Thus, the airfoil lift is:  [ydM-y M]iAM}  + ^{[y«(*i)-y6(*i)][«^(wi)-i] +  e  _  - — |  _  ao(Vio ^20-«^3o) + y ( ^ n - ^ 2 i - ^ 3 i ) + y ( ^ i 2 ^ 2 2 - ^ 3 2 ) H u r ^1 T + O(W30 + Y  31  g  ^2 T "3" "*  +  j w *"  3 2  -  0  4 0+  _  di _ (^2 T "2" If ^  (4.69)  +  _  + o(^20 ^4o) + y(^21 ^4l) + y(^22 ^42) H +  T" (T"' T-)  - [*(/)-*(/)]£- [»(/>-*(/)]£]}.  Following similar procedures for the airfoil moment, one gets: 2  G G  - ^{[y.CxO-nMIWCwi)-!] + + ^2~|  4. °  + y ^ 7l 4. +1 7  4.  7 j. + ~£ M • 6  2  J  + y(^2i-^4i) + - ^ ( ^ 2 2 - ^ 4 2 ) +  d  -  6  4-  -  3  ~ T  [Y {x )-Y (x )]iAM)  y(^ii ^2i-^3i) + j t ^ - ^ 2 2 - ^ 3 2 ) b  i  1  r°y 7 v/  d  4.  2  e  + -^(^13-^23-^33)  +  (4.70)  l7 4. ^ 7 4. + — */ + "^-'43 H  rf  4 1  2  + -^(^23-^43)  2  42  H  ir ( T - " T " ) M - W ) -fc(/)-n(/)]£- [W)-W)]£]}  where, in these last two equations,  ut  = 171.,  w = § +1'7 , (and from Fig. 4c) A = 2  x  2  iu  u  and A2 = i f + tV2 are all on the upper surfaces of the airfoil and flap. It is important to appreciate the generality and versatility of these last two equations and the analogous expressions for overall lift and moment, Eqs. 4.60 and 4.61. They contain the force solutions for all the analytical boundary conditions discussed so far. The a and i] problems are solved by simply choosing the appropriate values for the polynomial coefficients in the boundary condition, Eq. 4.47. If o„ = b , d = e„, or if n  a  n  n  = b = d = e , then the force expressions simplify in a straightforward way. These n  n  n  Tandem Airfoil Forces an Easy Way / 4.3  82  coefficients are the "geometric amplitude parameters" mentioned in Chapter 1 and the coefficients from the selector function expansions, the "force solution constants", measure their relative effects on the forces. These expressions are exact analytical solutions to the linearized boundary condition problem. Note that the segment selector functions allow an airfoil element to be divided into any number of discrete segments. They allow, for example, the solution of a problem where both the airfoil and flap have attached flaps and/or slats at arbitrary deflections. Or, these selector functions could be used to describe numerically specified camber lines by fitting n  th  order polynomials to any number of segments, each containing n + 1 camber  line coordinates. In this case, a closed form analytical solution to a numerically specified problem would be achieved. It is questionable, however, as to whether the numerical problem is best handled by selector functions or by simply numerically integrating the force integrals of the previous section. It is convenient to have the a and n problem solutions in terms of the selector function expansion coefficients, as provided on the next page in Equations 4 . 7 1 . All these equations have been shown to be identical to those derived in Chapter 2, Eqs. 2 . 4 7 . Rather than calculate leading edge moment variations with a and n, it is probably more meaningful to calculate aerodynamic centers and their associated pitching moments. If:  defines j J j x  a  c  dC  which is associated with x /(s+l—f), ac  x /s Sac  }  and (x/ —/)/(l — f), then ac  the aerodynamic center characteristics are as shown in Equations 4 . 7 3 . A C refers to the shift in a coefficient's value at a = n = 0° due to the presence of a camber line. In the next section, all of the above aerodynamic characteristics for tandem NACA airfoils are evaluated and the results compared with exact-numerical calculations.  Tandem Airfoil Forces an Easy Way / 4.3  a)  b)  c)  d)  dC  An  t  da  s+l-f  dCt  Air  dr,  S+l-f  dC ,  S+l-f  t  da da  LL dr,  h)  s  i  _ s+l-f 1-f  dC  s+l-f  dC  tf  j  -  r  s  +  1  _  la  s  t  )  dC  e>  1-f  ^  ^  2  =  (4.71)  dr, 1 +  -2  4(/- )(l-/) S  {s+l-f?  Air  2 dr,  dr,  s  7T /  Zl2  {s+l-f)*  ^  M  ,  da  mp  S + l-f  _  2s  G  da  dr,  7o  da dC, m  fo  S+1  = (  _  dr,  f  ~)  N R  d  g  mo  -  (-^y)  .  2  ac  2  2K'  [2  2 d C m $ 0  f2K'  _  ir \ ir  21  f  '  (W)2^f-fe) da  X  2  ~da ~ a ~2~  r  2TT |V 2  2JT JH  1—3 ^ 2  / \ i - v 1 < - [G - — —  4  3  1)  2  s+1-.  (1 - E )  * 0  dC„  :  f  1-f da  dr,  ldC,  mo  m  (  2 E ' \ ,  dC  Alt =  da  da  s  t  da  1-f  0 ^ j)  or,  "»0  da  s 2K'f2K'  dr, da  g)  oa  -  da 3C  *)  Hi  s  da e)  83  2E'\ / ir  2  -3  2  J  d C t j  1-f  da  1-f  dr,  mo  da l™±  »ac dC,  ac dr,  da dC JO dr,  /  m  m  _j_  da  dC.xj Jac dr, c  . *  x -f  =  fac  1-/ AC  m m  dC  da •  3ac  /dC  mfo  /  tj  da  = ACL,- + A C / . m  Jo ^  ?  j  (4.73)  Tandem Airfoil Forces an Easy Way / 4.4  4.4  84  T A N D E M N A C A AIRFOILS, A N E X A M P L E  The NACA 4 and 5 digit airfoils have camber lines (Abbott and von Doenhoff, 1959, Chapter 6) that are described by two different polynomials: ^  = a' + oil* + o' x' 0  2  ,  2  0 < x' < x\ (4.74)  d x  = 6Q + b\x'  ,  xi < x' < 1 .  Although it is not a requirement of thin airfoil theory, these polynomials are matched so that the magnitude of dy' /dx' is continuous through x' = x\. The coefficients in the c  polynomial boundary conditions (Eqs. 4.47) for tandem NACA airfoils become: a = a' , 0  A a' 0 = 0 d  O i = o'./s ,  Q  a  YZ7a'l ++ (\JZ~f) \ f  f  a  a = a /s 2  2  ;  2  ~ Y—j,' l  \ ' A °2 » l d  1  a  b = b' , Q  ~  '  Q  ,' 0  (1=7)2 2 ' 2  f  6i = b\/s A 2 ~ d  1  Ja  (IZ7)2 2  (4.75)  o = d - jzy^i ' i JZ/ '  e  6  Xi=axi,  c  =  6 1  x = / + (l-/)xi 2  These coefficients illustrate what could become a serious problem when they are used in the expressions for the aerodynamic forces of the previous section. That is, when the slot is large or one of the airfoil element chord lengths is small, the coefficients become large. If cy is  is the smallest of s and 1 — / , then the largest coefficient in Equation 4.75  0(V y )* Thus, c  mln  if £ is the error in the selector function expansion coefficients, then  the maximum error in an aerodynamic coefficient will be:  •(f)  (4.76)  since the pitching moment coefficients are nondimensionalized by Cy. Further, for the NACA 23012 airfoil, A C  m o c  is O ( l 0 ) and so the percentage error in its value would be: -2  \  Jmin /  It is therefore important to keep £ as small as possible.  Tandem Airfoil Forces an Easy Way / 4.4  85  The segment selector expansion coefficients, S and T„, are different from the n  other expansion coefficients, A and 7 (which solve the a and n problems in Chapter 2), in n  n  requiring the evaluation of elliptic functions and integrals with imaginary arguments (such as when ui =  =  The procedures of section B.2, which were sufficiently accurate  when evaluating elliptic functions with real arguments, lose accuracy when the arguments are imaginary. For example, the first truncated term in the expression for X) , Eqs. 2.48k 3  and B26, is: 0( Replacing Xoo  g  1 6  sinl6  X o o  )-0(g  1 6  ).  :'7i makes this error:  and if x± is close to zero, ifi will be close to  and the error will be 0(q ). This loss of &  accuracy, combined with the previously mentioned need for greater accuracy, motivated the development of the highly accurate procedures in section B.3 for calculating Wj, u , 2  ^( i)> ^( 2)» SO^oo)* w  w  a n c  * 7"( oo) 3 w  8 w e  ^  *  3 8  n e  usual constants q, G G , and E . lt  2  3  The following is a step-by-step procedure by which the aerodynamic coefficients of the previous section can be evaluated for polynomial boundary conditions of the type given in Equation 4.47. At the end of each step, the order of magnitude of the error in that step, £,, is noted. Given s and / , calculate: 1) k and e from Eqs. 2.48a &b. £i = 0 0  2) /tj and q from Eqs. B27,28,<k29, using only the first two terms of Eq. B28. 2K  2K 2E 1  3) — , IT  4) 2c  , IT  Xoo  1  , and then G and G using Eqs. 2.48d, e, f, and i & j. £ = 0(Z2q ) x  2  Z2  z  IT  , dn(/i  O02  , k ), and Xoo fr° 2  terms of Eq. B32. £ = 0(q )  m  Eq - B23,30,31,32,«k33, using only thefirsttwo 8  Z2  4  5) - ^ ( 2 X 0 0 , q % - r - ( X o o ) , 3 & l - B37,38,&26. Only the term 0(q ) from v& t/ the series in Eq. B38 is included. £ = 0(l6q ) a  n  dE  u s i n  E c  s  iC  3  h  32  Tan Jem Airfoil Forces an Easy Way / 4.4  86  6) e , e , dn(A ,A: ), dn(A ,fc ), ( M, Wl  W2  l2  2  22  2  7i,and 7 from Eqs. B35,30,31,32,  2  2  &36, using only the first two terms of Eq. B36. If £  2 m a x  is the maximum of £ (wi) 2  and £ (w ), each of which is always positive, then £ = <2(32£ 2  2  6  2mix  <j ) 32  7) i -(2w ,g ) and -r -(2w ,g ) from Eqs. B38,39, & 25, truncating the series in Eq. B38 3  4  1  3  2  4  after the 0(<? ) term. Then, Eqs. B15 are used to get the expression: 16  .  „  ,  2w  K'#  . , __ __ _ i  M  ,  l f  M  ,  2w  JiT'i?,, ,  . , __ __ j  and using Eq. B25, one can show that  £ _ _  Ji  2  / _ . 5V x  j + \]~7u~ j x  — Xj  1/Zy + V - ^ l  ( 4  .  7 7 )  K' OL  1  K  x  ~> j)  ^3  "  (1-3)(1-/)  0(l6^ q ) mix  32  0  s  9) 5 X and T from Eqs. C12& 17. £ = 0 ( l 6 £ 9  x  2 m a x  mix  g ) 36  g ) 3 2  $nt and 7J, for n = 0 to 4 for the NACA camber lines, using the recurrence  relations of Appendix C. £ 11) £7 , LS n  _  2  8) 5 0 and T from Eqs. C8 & 16 and B40&41. £ = 0(^|  10) -Cm  _  — l] and i>t(w ) are equal to:  1 _ x  A , '  7  _  cnA  G  3  where; = 1,2. £ =  M  2tf'  ni  l 0  = 0  and LT from Eqs. 4.58&59. J n  mn  from Eqs. 4.67&63. m = 1 to 4 and  n — 0 to 3 for the NACA camber lines. £u = 0 12) all the aerodynamic characteristics from Eqs. 4.60,61,69,70,71, & 73. £  = 1 2  0  13) The maximum error in the above quantities will be: (4.79) where p is the order of the highest order polynomial in the boundary condition. Of course, this does not include the roundoff error associated with these computations, determined by the number of digits the computer or calculator uses. The above procedure has been programmed into Texas Instrument's Compact Computer 40 (CC-40), a textbook size portable microcomputer (34K of ROM and  Tandem Airfoil Forces an Easy Way / 4.4  87  6K of RAM) using TPs enhanced basic language. Calculations use a minimum of 13 digits per number. The program, which is listed in Appendix D, utilizes most of the available R A M memory and takes a little over half a minute to solve for all the aerodynamic characteristics. As a comparison, the procedures of Chapter 2, which solve only the a and TJ problems, use about IK of memory when programmed into this computer and are evaluated in only 6 or 7 seconds. Thus, in order to have the additional capability of calculating A C ^ and A C  m j a c  , all of which are relatively small forces in most situations, about five  times as much memory and computation time are needed. The thin airfoil theory predictions for the variation of the aerodynamic characteristics with geometry for tandem NACA 23012 airfoils are presented in the next six figures. The figures also make comparisons with exact-numerical potential flow theory (Appendix A). For the numerical theory, the variations of the coefficients with a and rj are approximated by assuming that they are constant over a 1° interval, beginning at a = T) = 0°. Also for the numerical theory, the aerodynamic center functions are calculated using the exact version of Equation 4.72: % ' C  mj  ac  where, again, dC j m  = mj0 C  +  (Qy  C  O  S  C  <* + dj  « H <*)  (4.80)  /da = 0 defines X y /cy. The exact-numerical predictions of the aero-  dynamic centers, Figure 13, are calculated in this way. For a proper comparison of C • m  in Figures 14 and 15, the exact-numerical predictions of Equation 4.80 are calculated using thin airfoil theory values for x- . Thus, thefigurescompare moment predictions about lac  the same point on the airfoil elements. In all cases, the error in the numerical calculations is estimated to be within the size of the symbol representing the point on the graph. In some instances, this required representing each airfoil element by as many as 100 straight line segments — as described in Appendix A. Parts (a) of Figures 10 to 15 show how the aerodynamic characteristics vary with slot size, with the flap chord to total chord ratio (Rf) held constant at a typical value  Tandem Airfoil Forces an Easy Way / 4.4  88  of 23%. The thin airfoil theory predictions, as indicated by the solid line curves, vary monotonically (except for 4 localized prediction curves) between the one-element theory's predictions at f — s = 0 (an airfoil with a simple flap) and f-s = 1 (two airfoils infinitely far apart). The small horizontal arrows in thefiguresindicate the predictions of the oneelement theory. Of course, all of the integral expressions for the tandem thin airfoil theory force characteristics have been shown to reduce to the one-element predictions at f—s — 0 and, with the TI CC-40 able to accurately calculate values for the tandem airfoil predictions at f—s values as small as 10~ where the predictions are virtually identical to those of the 6  one-element theory, there is now the necessary confirmation that the analytical solutions to the force integrals are also behaving properly at this important limit. For f — s values approaching 1, the calculations are limited in their accuracy by the error £  m a x  , Eq. 4.79.  For Figures 12 and 15, p = 2 and the calculations are only accurate up to / — s « 0.93, although there is little doubt that the curves are approaching the one-element prediction at f — s = 1. For the other figures, the a and rj problem results, | ^ 2 roundoff error limits the accuracy of the calculations before £  m a x  max  l —  P  =  ®> ^ an<  does.  As afinalindication that the present theory is indeed a correct representation of the flow about tandem thin airfoils, the exact-numerical calculations were also performed for 1% thick airfoils (23001 airfoils) and in all cases the numerical predictions collapsed to the linearized predictions (approximately 90% closer than the 23012 numerical predictions). This also happened when these "thin" exact-numerical calculations were performed for the curves infigures(b). Parts (b) of Figures 10 to 15 show the variations of the same aerodynamic characteristics with flap chord to total chord ratio, while keeping the ratio K'/K constant at 1.3. This is accomplished by simply keeping the modulus k constant:  This value was chosen from figures (a), where f — s = 0.05 is a slot size sufficiently far from the one-element extreme at f—s = 0 to make the comparisons worthwhile yet not so  Tandem Airfoil Forces an Easy Way  F I G U R E 10  / 4.4  Lift curve slope characteristics of tandem NACA 23012 airfoils. linearized theory; •, o, • exact-numerical theory; j-* one-element TAT prediction for j airfoil element. th  Tandem Airfoil Forces an Easy Way /  F I G U R E 11  4.4  Flap defection effects on lift for tandem NACA 23012 airfoils. linearized theory; •, o, • exact-numerical theory; j-+ one-element TAT prediction for j airfoil element. th  Tandem  Airfoil  Forces  an Easy  Way  4.4  /  0.16  0.14  -  0.12  0.08  0.06  0.04  a)  AC .  1  0.16  vs. f-s;  t  R = ( l - / ) / ( a + l - / ) = 0.2308 f  i  1  i  -  °  o „  I  \oo  ° m  o  6  i  i _ _ _ J  0  o  tc,.  -  ^  o  \ \  o  m m  •• ' *  •  V* • \ • \ •  J  .  0.12 ' — -  Aq  •  o  \ 0  0.08  \  •  -  o  1 --0.2 AC.  0.06  **/  •  o •  0.04  (> -0.6  -0.8  7  i  i  D  0.2  b)  F I G U R E 12  i  i  Camber  M  AC  tj  R,  vs. Rf]  one-element  theory; TAT  •  i  i  0.6  0.8  1.0  f—s = 0.05, K' = l.ZK  line lift characteristics linearized  j—*  i  of tandem  NACA  •, o, • exact-numerical  prediction  for;  t h  airfoil  23012  airfoils.  theory; element.  Tandem Airfoil Forces an Easy Way /  F I G U R E 13  4.4  The aerodynamic linearized j-*  centers of theory;  tandem NACA 23012 airfoils. °, • exact-numericai theory;  one-element TAT prediction  for 3  t h  airfoil element.  Tandem Airfoil Forces an Easy Way  F I G U R E 14  / 4.4  93  Flap deflection effects on moment for tandem NACA 23012 airfoils. linearized theory; •, °> • exact-numerical theory; j—* one-element TAT prediction for ; airfoil element. th  Tandem Airfoil Forces an Easy Way / 4.4  4.02  1  1  94  1 —  -i  r  — i  i -  1  1  —  - 1 . /  -0.01 I / O  If*  o •  •  •  a  9  S  •  •  •  o  •  *> 3  -  1.01' I1  1  1  1  0  0.2  1  1  1  0.6  0.4  • 0.8  i 1.0  /•J  vs.f-3;  •<f  *  r-  *  = (!-/)/(,+  !-/) = 0.2308  o««c  \©  S  © _1  i•  l_  o b _J  « 2  o o  I  r«  o 6  n  o o •  —  o o  o I  -0.02  AC, ml  cc  0.01  0.03  b)  F I G U R E 15  AC y m  a c  vs. R ; F  / - s = 0 . 0 5 , K' = 1.3iT  Camber line moment characteristics of tandem NACA 23012 airfoils. linearized theory; •, o, • exact-numerical theory; j—* one-element TAT prediction for j ^ airfoil element. 1  Tandem Airfoil Forces an Easy Way / 4.4  95  large that it is too far from typical slot sizes for slotted slat orflapairfoils. The functions shown in parts (a) and (b) of the figures, then, have common values at f — s = 0.05 and Rf = 0.2308. Note that for figures (b), f-s  is minimum at Rf = 0 and Rf = 1  (where it is zero) and maximum at Rf = 0.5. Figure 15b includes a scale showing the variation in f-s over the range of Rf values. Where appropriate, figures (b) also include the horizontal arrows indicating the one-element theory's prediction of the limit point at Rf = 0 or 1. When calculating the tandem thin airfoil theory predictions forfigures(b), the largest £  m a x  was O ( l 0  -31  ) . However, roundoff error became a problem on the CC-40  for extreme Rf values (Rf < 0.01 and Rf > 0.99).  Figures 10 to 15 show that tandem thin airfoil theory provides good quantitative predictions of the aerodynamic characteristics of the exact potential flow. Also, the linearized theory's predictions of the qualitative aspects of these characteristics is excellent, with the only discrepancies occurring in extreme geometries which are of little practical interest in aeronautics. It is particularly encouraging to see the consistency with which the overall forces are modelled, as this is a good indication that the linearized theory's predictions will correlate well with real flow forces. Forfigures(a), the quantitative relationship between the overall force predictions of the exact and linearized theories at the one-element extremes is maintained throughout the range of f—s values. This relationship is also maintained for the overall lift variations offigures(b). The overall moment characteristics offigures(b) do exhibit this consistent relationship to some extent, but with a great deal of sensitivity to the overall thickness of the configuration (the overall thickness is greater at Rf = 0 and Rf = 1 than it is at Rf = 0.5). Not surprisingly, the localized forces are less consistent than the overall forces when thin airfoil theory predictions are compared with exact potential theory, although this does not necessarily mean that the linearized predictions will be equally inconsistent when  Tandem Airfoil Forces an Easy Way / 4.4  96  compared with the real flow forces. Indeed, Foster et al (1970) include lift comparisons for their 31% slotted flap configuration which show that the experimental values of the airfoil's localized lift are overpredicted by exact potential theory by about the same amount as the configuration's overall lift is overpredicted (note the linearized and exact theories' relationships for C and C in Figures 10 and 11), whereas experimental values of the flap's t  tt  lift were very closely predicted by the exact potential theory (which is, again, compatible with the relationships of Figures 10 and 11). The camber line force prediction comparisons (Figs. 1 2 and 15) are apparently the least accurate of all the figures. However, these forces tend to be small in comparison to the total forces on the airfoil, as shown in Figures 16 and 17. Also, it must be remembered that the localized force coefficients have been nondimensionalized using the localized chord length so that, for example, if Rf —* 1 (so that s —• 0), then A C ^ , AC  m # a c  , and their errors are magnified relative to the actual forces. Figures 16 and 17 show the actual nondimensionalized forces and their vari-  ations with a and ij for Rf = 0.2308 and for f-s = 0.025 (figures (a), a realistic slot size) and /—s = 0.1333 (figures (b), an interesting midpoint between f—s = 0 and f—s = 1 where K = K'). For Figures 17,  is again calculated from Equation 4.80 using the  thin airfoil theory values for x- . Except for high values of a+rj, the linearized predictions Jac  of the overall forces agree well with the exact potentialflowforces, with the relationship between the exact and linearized theories being much the same as for the one-element 23012 airfoil of Figure 1. Although some of the localized linearized moment predictions show poor agreement with the exact predictions, it i3 in proportion to their decreasing magnitudes and, therefore, importance. Thus there is strong evidence that the consistent close relationship between the linearized and exact overall force predictions seen in one-element airfoil theory is maintained in the tandem two-element regime as well, so that one would expect the tandem thin airfoil theory overall force predictions to again be closer to the realflowoverall forces  F I G U R E 16  C  (  vs. a for different n for tandem NACA 23012 airfoils. R = 0.2308 linearized theory; •--<>--• exact-numerical theory f  J  1  '  5  "  1*  10  IS  1  I  -5  1  0  5  1  L*.  10  15  a°  a  a)  f-s  F I G U R E 17  I  -5  1 0  1 5  o°  = 0.025  (  10  !_». 15  |  -5  | 0  | 5  |  10  •  15  -5  I 0  • 5  •  •  10  15  -5  o°  a°  b)  f-s  0  • 5  I  I,  10  15  a°  = 0.133  C , vs. a for different TJ for tandem NACA 23012 airfoils. R, = 0.2308 linearized theory; «--o--• exact-numerical theory m  I  <P  Tandem Airfoil Forces an Easy Way /  4.4  99  than are exact potential flow theory predictions. Furthermore, the linearized predictions should correlate better with the real flow forces than they did in the one-element regime since properly designed multi-element airfoils reduce boundary layer effects on the overall flow. This completes the discussion of tandem thin airfoil theory. There is much else that could be done with it, such as finding the most efficient way of handling numerically specified airfoil profiles, or conducting a series of experiments on tandem airfoils in order to directly compare the linearized predictions and real flow results. However, since the tandem airfoil configuration is not the most common of configurations, there is little motivation for giving it more attention. The present theory has served its purpose by showing what can be accomplished by exploiting the theory of elliptic functions in conjunction with the use of modern hand calculators or microcomputers. In the next chapter, a beginning is made on the general linearized theory for staggered airfoils.  CHAPTER 5 Staggered Two-Element Thin Airfoil Theory  In this chapter, expressions for the general, two-element, linearized airfoil theory forces are derived. Overlap is accounted for so that the two airfoil elements of arbitrary chord length can be arbitrarily positioned relative to one another, as long as each airfoil element's chord line does not deviate appreciably from an orientation parallel to the direction of the flow at infinity — the linearization requirement. This flexibilty comes at a cost. As previous investigators have discovered (Munk, 1922b; Glauert, 1926; Millikan, 1930; Garrick, 1936), the staggered slit mapping is much more difficult to solve than the tandem slit mapping. Solving for its parameters requires a numerical approach that is not pursued to completion in this thesis. The force solutions are in terms of these mapping parameters.  5.1  T H E S T A G G E R E D P A R A L L E L SLIT M A P P I N G  This mapping is shown in Figure 18. The linearized physical plane is again the z-plane, Fig. 18a, with the airfoil slit on the real axis stretching between x  x  = 0 and x . n  The  flap slit is displaced downwards a distance h from the real axis, with its leading edge at z = X2 — ih and its trailing edge at z 2 2  2  =  X2 — ih. When the plane potential flow 2  about these slits is mapped into the infinite strip in the A-plane (0 < / i < K), Fig. 18c, the leading and trailing edges of the airfoil and flap no longer coincide with the corners of 100  Staggered  Two-Element  Thin  Airfoil  Theory  /  5.1  101  X -PLANE t -PLANE  z -PLANE  V A XQ  x  xii  U 22  h  hi  t  22  222  Cb)  (a)  F I G U R E 18  The staggered parallel slit mappings.  the rectangle as in the tandem airfoil mapping. The airfoil slit still opens to form the left side of the rectangle but its leading edge now maps to the point Aj and its trailing edge to A . The flap forms the right side of the rectangle, with its leading and trailing edges n  at A and A . 2  22  The f-plane, Fig. 18b, is a convenient secondary solution plane which helps with some of the mathematical formulations that follow. The leading and trailing edges of the airfoil and flap correspond to the points t  u  t  iu  t , and * . If h is positive, as shown 2  22  in Figure 18, then it turns out that t± and t are on the bottom surfaces of the slits in the 2  4-plane while tn and t  22  are on the top surfaces.  The mapping from the <-plane to the A-plane is the same tandem slit mapping used in Chapter 2, Eq. 2.4, with A = ±iK\ 0, K, and K ± iK' corresponding to t = 0, a, / , and 1 as before. In the z-plane, these points correspond to z = x , x , Zj — Xj — ih, 0  t  and z — x — ih. For the geometry shown, x and Zj will be on the top surfaces of the e  c  Q  airfoil and flap slits while x and z will be on the bottom surfaces. t  c  The mapping from the z-plane to the A-plane is obtained by considering the  Staggered Two-Element Thin Airfoil Theory / 5.1  FIGURE 19  102  Doublet How in the X-plane corresponding to zero incident uniform Bow in the z-plane.  exact formulation for the zero incident uniform flow, {7, about the slits in the z-plane: w(z) = U .  (5.1)  This is correct since the slits do not obstruct or redirect the flow in any way. Integrating tv(z) gives the exact complex potential, fl(z), which has the same value at corresponding points of a conformal mapping:  n(a) = Uz = n(z(A)).  (5.2)  In the A-plane, the corresponding flow is certainly not uniform as there must be a doublet at A = A^, corresponding to z — oo. This is the only singularity that can be present in theflowin the rectangle. The lines fi = 0 and p = K are solid boundaries across which there is no flow, a condition which is satisfied by placing an infinite series of image doublets about these lines, as shown in Figure 19. Since the flow must also be periodic in 2iK', this pattern of singularities continuously repeats itself up and down the imaginary axis. The generality of the flow is ensured by allowing the doublets to be oriented at any angle /?, and allowing for a superposed vertical flow over the plane.  Staggered Two-Element  Thin Airfoil  Theory  / 5.1  103  Consider the complex velocity: w(X) = C [e sn (A+/*«,) if}  sn  2  2  ( A - + tfc]  (5.3)  where C and b are unknown real constants. The singularities of the sn function are located along A = fi ± i K ' , v/here sn(jt ± i K ' ) = 1/k sn p: w(fi ±  iK')  = C  +  (5.4)  ib  The expansion of sn A for small A is: sn A ~ A + O (A )  ,  3  so that:  as A -* 0  ^  c  (5.5)  These are the complex velocities for four doublets of equal strength symmetrically placed about the imaginary axis in the A-plane and, since sn (A db 2K ± 2iK') = sn A, w(X) is 2  2  indeed the complex velocity representing the doubletflowof Figure 19. As a further check, one can use the addition theorem for the sn function to show that to (A) is purely imaginary along the lines ft = 0 and fi — K so that, as anticipated, theflowdoes not cross the surfaces of the opened airfoil and flap slits. Integrating to(A) givesfl(A)which must be equivalent to fl(z), Eq. 5.2, at corresponding points. Upon choosing an appropriate constant of integration, the staggered slit mapping is:  A *(A) =  +  \e sn (X+fL ) iff  - -^sn (A-  2  e  00  2  M o o  ) + ib] dX .  (5.6)  o As in Whittaker and Watson (1927, 22.73), it is convenient to choose:  A E{X) = Jdn XdX 2  (5.7)  Staggered Two-Element  Thin Airtoil Theory  /  5.1  104  as the fundamental elliptic integral of the second kind. Then, using the relation between the sn and dn functions (Eqs. 2.10) and absorbing terms into new real constants Cj and b  lt  z(X) becomes: z{\) = x, + d { c o s / ? [ £ ( A - / 0  -  ^(A+Moo)  ib \-i8\nP[E{\- i )  +  l  +  2£(/0]  + E{\+ii )]}.  (5.8)  00  l 00  The first line on the RHS of Equation 5.8 is periodic in 2tK': F(\  n  \  F(\±,  ^ _ 2cn/x dn/x 00  l  A; snVcom A 2  o0  2  while the second line can be made so by choosing an appropriate value for 6^ Since: % £ ( A - ^ )  +  -/»  v  £ ( A + ^ )  =  _^./vx 2  2cnAdnA  £(A)  m  k sn « „ sn A L  X  2  2  f  J  ,  2  .  ^  (  A  ' -  . 1  0  )  and: JE(A ± 2ii)C) = £(A) ± 2i(K' - E'),  (5.11)  then 6, must be:  6 = 2sin/?(l-|Q.  (5.12)  i  Also, on theflapwhere A = K + iv, the imaginary part of Equation 5.8 must equal — ih. On the flap, Equation 5.9 is purely real and Equation 5.10 is purely imaginary except for one term: E(K + iv) = E + imaginary terms so that, using Legendre's relation, Eq. B16: Cx = — • TT sin/3  (5.13) '  v  1  Thus, the mapping is: Z(AJ =  X  -i  i  cn  2K' h a — - < COS/7 i smjJI  dn fi^ sn Poo  • • * 17,  +  k sn 2  l-k ^  sn A T ~ T sn* y.^ sn A  2  TTj l  x  fi^  2  5  2  n/n  ,sm/?[(l--JA-S(A)  +  cnAdnA  ji-sn ^J^^J} A sn 2  2  2  A l)  (5.14)  Staggered Two-Element Thin Airfoil Theory  / 5.1  105  which, for /? ^ 0, is only quasi-doubly periodic. This is more easily seen when the doubly periodic components of this equation are written in terms of the ^-variable and the quasidoubly periodic term is written in terms of a -^-(w) function. One gets: H Z  =  XQ  +  sin/?  ^(t-s)(t-f)(t-l) t cos 0 - i sin /?  where:  -ih  ir  K  }  tV  2K'  H  (5.15)  (5.16)  Of course, z(t) is not infinite at f = 0. The square-root pole in the first square-bracketed term of Equation 5.15 is cancelled by that in the -^(w) function. When h and /? go to zero, this mapping reduces to the tandem slit mapping, although it is not yet obvious that iT/sin/? reduces to the correct value. Note that when h —* 0 with overlap (/? •/+ 0), none of the terms in Equation 5.15 disappear since K'  co. When /? = TT/2, the airfoil is in an  unstaggered biplane configuration. The eight unknowns in the z(X) mapping: x  »>  P> > /^ooi i » k  *2i 2 2  A  ( -  A  5  1 7  )  are determined by simultaneously solving eight equations: *  dz ^ ( A ) = 0;  z(\ ) = z , j  j  ;  (5.18)  ; = 1,11,2,22.  In the f-plane, the k, fi^, and Xj variables are replaced with / , s, and tj. When z is written entirely in terms of the f-variable, it must be written as an integral, which can be conveniently put in the form:  Z(t)  = XQ  H +  sin/?  sin/? *cos/? — i  « 2 * - ( l + s+/)* + /  /  2  1-<i-.>!  y/t{t-a){t-f)(t-l)  dt } . (5.19)  dz The -r-[tj) = 0 requirement becomes: dt cos{3^/tj(tj-s)(tj-f)(tj-l) where:  6s f 1 -  = »sin/? F'  il-») , Y  2  l+s+/  j  2  6 >' 2 +  (5.20)  (5.21)  Staggered Two-Element  Thin Airfoil  Theory  / 5.1  106  Squaring each side of Equation 5.20 produces the final relation for the iy's:  <J.-(l+a/)*J.+ +  (/+  3+3  /)cos /? + 2  8in /3 ( ^ 2  sfcos p + 2  1  +  s  + ^  +  f  »  8in p^-^-6 tj + s'm2 p— = 0 . 2  (5.22)  In the f-plane, the uniformflowof the z-plane, Eq. 5.1, is transformed into a flow incident to the tandem slits at an angle of attack a = p. The above equation for the tj's gives the stagnation points of that flow. For the general problem, this quartic polynomial does not have simple factors but can be solved analytically using a method given in most mathematical handbooks. However, this first requires knowledge of the first four variables of Equation 5.17. These variables can be solved only through a trial and error process involving the repetitive evaluation of elliptic functions (see the last term of Eq. 5.15). Although the methods of section B.3 would greatly facilitate this calculation, the calculation time needed on a small microcomputer would still be relatively long. The problem could be easily handled on a large mainframe computer but then it can be argued that one might just as easily solve the exact-numerical problem. In the remainder of this section, alternative ways of posing the staggered parallel slit mapping problem are considered in the hope that one might simplify the problem's solution.  Suppose the angles (3 of the doublets in Figure 19 were set to zero. The stagnation points of the flow would then move to the corners of the rectangle, A = 0, and A = K, and the airfoil and flap boundaries would be enclosed by the same streamline. In the *-plane, the angle of attack of this new flow relative to the tandem slits would be zero. In the z-plane, the angle of attack would be a = —/? and the stagnation points in  Staggered Two-Element Thin Airfoil Theory /  this flow would be at z = x  Q}  5.1  107  x , Xj-ih> and x -ih, t  as shown in Figure 18a. Using the  c  mapping, Eq. 5.15, one can show that:  X  '-  X o  =  ^  '  3  Then, with k given by:  X f  -  X  ° = T*Vp ' f  X c  -  X 0 =  T*h-  (5  23)  -  '  .2  =  /-«  =  (l-0)(/-3)  /(!_,)  (/_o)(l- )' s  the corresponding points to-0, s, / , and 1 in the z-plane also give: [x -x )(x -x,) f  0  c  and: s n  2  = l-s = ° ~ " . x - x X  X  c  (5.25)  0  Similar relationships for c n ^ and dn /x^ are obtained. J^Note that the last relation of 2  2  Equations 5.23 can be rewritten as: 2K'  h  K  tan/3  = >/(*/-*<>)(*«-*.) •  (5-26)  It is. now easy to see, in Equation 5.15, that this constant reduces to the correct value needed for the tandem slit mapping as /3, h —> 0 .j Thus, one way of obtaining the first four variables of Equation 5.17, without having to deal directly with the last four, is to solve for that unique circulationless flow in the z-plane which encloses both the airfoil and flap slits within the same streamline. Although this would still require a numerical approach, perhaps this could be more easily handled in the physical plane where all the geometry is known. A method for conformally mapping general airfoil shapes to circles, such as by Ives (1976), might prove useful, especially since the staggered parallel slit geometry is so simple. Another approach to solving for the mapping parameters is one which takes advantage of the properties of the theta functions.  Since E(X) is closely related to  Staggered Two-Element Thin Airfoil Theory / 5.1  108  the -r^(w,g) function, the mapping, Eq. 5.8, can be rewritten: #4  z(u) - x  t  27ihb  l n 9  {  cos/3  t?' ^ ( w - X o o . 9) - ^ ( w + X o o , 9) + 2^(Xoo>q)  + xsin/3  >9  }  t?  4  (5.27)  or, using the u'-plane and q' parameter (Fig. 4e): cos/3 t sin (3  fa-Woo,  fl  2  +  (5.28)  fa+i-LJ)  2  This last equation is the simplest form the mapping can take on, with the mq terms and the term linear in w from Equation 5.27 being absent. For large gaps, q' —* 0 and it should be possible to invert the q' series representation of Equation 5.28 and obtain u' as a function of z. However, since slotted flap airfoils have very small gaps, these q' series would be inefficient (q' —» 1 as the gap goes to zero). This problem is further compounded by the fact that, for a given separation distance of the airfoil elements, q' is closer to 1 in a staggered slit than in a tandem slit configuration. The corollary of this is that q will be closer to 0 — so the q series representations of Equation 5.27 should be very efficient indeed. Since Equation 5.27 is in terms of u, \nq, and series of sinw and q terms, a trial and error or double series solution is unavoidable (the familiar x = In i , solve for x problem). Nevertheless, a procedure for solving Equations 5.18 by taking advantage of the q series expansion of Equation 5.27 has been developed for the biplane configuration (/3 = 7r/2 => 4 equations in 4 unknowns, including a trial and error solution in 2 of the unknowns). The procedure uses Newton's method of successive approximations and occupies virtually all the available memory in a TI programmable 59 hand calculator, which takes from l\ to 2 minutes to solve for the four parameters k, fi^, A ( = —Ax), n  and A 2( = A ). Although it neglects terms O(g ), the procedure is accurate for gap sizes 2  2  as large as the airfoil elements themselves.  4  Staggered Two-Element Thin Airfoil Theory / 5 . 1  109  Besides showing the inadequacy of today's hand calculators for solving the general staggered slit thin airfoil theory problem, the biplane mapping solution produced some hard information on the values of the important parameters of the staggered slit mapping at this, the opposite extreme to the tandem slit mapping. Table 2 shows some of these values for a 23% "flap". Note the extremely small values of both f-s  and q for  typical slotted flap gap sizes of from 1% to 4%, although in a true slotted flap configuration P would likely be closer to 0 than JT/2. Nevertheless, f—s and q appear to be consistently the smallest of all the staggered slit mapping parameters and an efficient general solution of this mapping must surely take advantage of this.  h  f-s —  :  —  s+l-f  ;  O  Tandem slit mapping  50%  15.4%  4%  1%  f  0.846  0.800  0.778  0.772  f-s  0.333  0.133  0.038  0.010  q  0.102  0.043  0.013  0.003  q'  0.013  0.043  0.103  0.175  0.369  0.343  0.326  0.321  f  0.834  0.890  0.951  0.970  f-s  0.152  0.007  1.4(10 )  4.1(10" )  Q  0.053  0.004  1.8(10 )  8.6(10" )  q'  0.035  0.163  0.474  0.789  0.323  0.219  0.142  0.112  u /K'  0.706  0.875  0.966  0.991  V22/K'  0.560  0.625  0.765  0.896  E'/K'  0.772  0.389  0.151  0.048  R  Ci + C  2  parameters, /? = 0  Biplane slit mapping parameters, n  TABLE 2  -6  -6  19  19  Comparisons of staggered slit mapping parameter values at the tandem and biplane slit extremes, C2/(c +c ) = 0.2308. 1  2  Staggered Two-Element Thin Airfoil Theory /  5.1  UO  Although q appears to be the smallest parameter for the greatest range of P values and gap/slot sizes, f—s is sufficiently small for typical slotted flap configurations that one is tempted to work directly with Equations 5.19 and 5.22 to obtain solutions for the mapping parameters by expanding these equations in powers of f—s. However:  2  jgt  K'~Hf=s)  a3/  a  - -*°  and so the 5 term (Eq. 5.21) in these equations creates the same problem for a series expansion solution of z(t) in powers of f—s as the lncj term does in the g-series expansion solution of z(u), Eq. 5.27. E'/K' - » O a s / — s -* 0 but at a very slow rate, as shown in Table 2. Still, it is useful to at least consider the zeroth order term of such an expansion. The following relations are obtained from Equation 5.19 in the limit that h —• 0 with overlap. They are not completely general since they assume x z  22 >  a  2  > 0 and  typical slotted flap arrangement. As f—s —• 0:  *(0 -*  1  ~ l °  z  K ) ~*  S  ^  X  2 2  + x (tcoaP - isinpy/t^t-l)) 22  ia'mP\/z{z-x )) .  ^—- + —(zcosP+ 2  22  (5.29) (5.30)  3C22  These are just simple Joukowski transformations of the single slit that results from the staggered slits coming together. The surfaces of the airfoil and flap that come together in the overlap region (x  2  <  x  x  < n)  a r e  l * *° the external flow. On the single slit o s  in the <-plane, the airfoil surface ends at t = / , with the top of the slit corresponding to r = i n at this point and the bottom to x = x . Thus, P is determined by requiring 2  that t(x ) = t(x ): n  2  z  t^p-*-]— \/X (x 2-X ) 2  2  2  "-*2  as/-s->0.  (5.31)  +\yx (x -X ) n  22  n  For a small gap staggered slit mapping, the differences between the above predictions of the mapping, Eqs. 5.29 to 5.31, and the true predictions will be 0(s'mP(E'/K')).  Staggered  Two-Element  Thin  Airfoil  Theory  /  111  5.2  Finally, using the expressions for z(t) and z(A), it can be shown that: z~  He-*.  sin/? •t + Z + -j- + -jjt + '-' > as*->oo 0  (5.32)  where:  Z = x + »(jr^y-^-fcE ) 0  0  (5.33)  8  and: *, = - i f f  4  i  f  +  3  +  S  }  )  ~  4  S  .  (5.34)  8 E  3  is as defined in Chapter 2. In the following sections, this expansion is used for evaluating  the forces by residues.  A recurrence relation for higher order terms can be obtained  from z(t), Eq. 5.19, by considering dz/dt  5.2  zst-*oo.  T H E OVERALL FORCES  In the remainder of this chapter, general expressions for the aerodynamic forces are obtained by generalizing the selector functions of the previous chapter to apply to staggered slit configurations. No doubt theflowitself could be solved for by finding the "general" and "special" functions analogous to those described in Chapter 3. Indeed, some of these tandem slit functions could be used on the "tandemized" staggered slits in the i-plane, with appropriate adjustments to account for the change in position of the leading and trailing edges and the fact that camber and thickness are no longer purely even or odd in any parameter. Also, the initial complicated mapping from the z-plane to the i-plane cannot be avoided. Although a general solution to the linearized staggered two-element airfoilflowfieldis not pursued further, the selector functions presented in this chapter are easily rearranged to give theflowfieldsolutions for those simple geometries of the a and n problems. Now, if the overall lift is: C = t  ci+c j 2  I w'(z) dz =  c  I w'(z)C(z) dz , c +c y 1  2  (5.35)  c  where u)'(z) is a general flow function that solves the a, r/, camber, and thickness prob-  Staggered Two-Element Thin Airfoil Theory / 5.2  112  lems, then the conditions on £{z), the generalized staggered slit version of the boundary condition selector, are: 1) its real part is zero on all airfoil element surfaces. 2) £(z)  —*  1 as z —• oo  3) the poles in £(z) must be at the airfoil and flap trailing edges so that w'(z)£(z) has no simple poles in the z-plane. Although not a problem for aeronautical applications, it now appears that thickness distributions must be restricted to those without infinite slope trailing edges. Otherwise, w[(z) will have a square-root pole at the trailing edge (see Eq. 3.12 for tandem ellipses) which would combine with a square-root pole in £(z) to form a simple pole, preventing the necessary unhindered collapse of C to the slit surfaces. This is not a problem in the tandem theory since the residue of this pole is real and so does not affect lift. For staggered thin airfoil theory, it is not obvious this would still be true. Again, £(z) is unique: £ { z )  Vt(t-s)(t-f)(t-l) (t-t )(t-t )  =  n  +  1 T  t  22~ ll  22  \/<22(<22-*)(<22-/)(t22-l) t— t  ~ 1)  \All(*ll ~ *)(*!! '""^11  22  . (5.36)  This function is obtained through a linear superposition of the -^-(w-Wn) and -r^(w—w ) 22  functions which put simple poles at the airfoil and flap trailing edges in the w-plane. These singularities also occur as simple poles in the f-plane but become square-root poles in the z-plane since the trailing edges are critical points of the z(t) mapping. In the f-plane, it is important to distinguish between tj and tj, the upper and lower surfaces of the tandem slits at t — tj. When h > 0 (0 < /? < TT), t and t are more properly written n  as t  xl  22  and t . Using Equation 4.35, it can be shown that the second line on the RHS 22  of Equation 5.36 combines with the first line to keep £(z) finite at t = t^ and t = t . 22  Alternatively, when h < 0, tj = t^t^ft^ ,kt  22  and £{z) isfiniteat t & t . xl  22  Staggered  Two-Element  Thin Airtoil  Theory  5.2  /  H3  To calculate the forces, the large z expansion of £(z) is needed. If: +  +  ^  +  ,  +  as*-oo  (5.37)  defines the A„'s, then the z(t) expansion, Eq. 5.32, can be used to obtain the £ ' s : n  £{z)  n  f>  p  ~l + -± + ±l + Ji + ... , Z  as z-> oo  Z>  Z  (5.38)  and: Cl==  £l  i2 =ZoCl +  sw >>  £t  etc  -^T >> -  ^  Taking advantage of the relations (from the stagnation point equations, Eqs. 5.20 and 5.22):  l+*+/ = *i + *n + *11 ~ <1 + *22 ~  2 cos/?  h  =  ~i \ Z * 2 sin/?  2  (5.40)  *2 + *22  ~ V^llC^ll ~ s)(*U-f){t - 1)  (*22-s)(*22-/)(*22-l)  n  '22  —  (5.41)  * l l + * 2 2 - ( * l + *2)(*ll+<22) + g 2 cos/? •-»  <22\/<22(^2-s)(<22-/)(*22-l)  sin/3  ~ hi\A11  - / ) ( ' 11 ~ 1)  <22-<n (5.42)  it can be shown that: £ 1  *n-ti+ti2-hja 2 cos/?  Ht -t t -t sin/? 2 cos/?  t  n  1 +  2 2  2  '  v  and: 4(/+s+ /)-(l + +/) -45 8 3  =  3  2  ^ + e -(f +f )(f +f ) + g ^ W ' 2  1  +  2  1  2  u  22  6  ( 5  -  4 4 )  When /? - » 0, #/sin/? is finite (Eq. 5.26) as are Equations 5.41 and 5.42. The LHS's of these latter two equations are used in C  tl  and £, in order to eliminate the indeterminacy t2  of the RHS's at /? = 0. However, when P —* n/2, the LHS's become indeterminate and need to be replaced by the RHS's. Now, the contour integral giving the overall lift can be evaluated on the surfaces of the staggered slits, where:  Staggered Two-Element Thin Airfoil Theory / 5.2  114  The overall moment analogue to this equation is somewhat more difficult to get. Although it is true that:  j> zw'(z)dz  when  C  collapses to the slit surfaces  ('i+c2)2  tol 4(*-£Ju/(*K(*)<fe  on the flap and  z = x—ih  t(z)  does not select the  boundary condition in one of the resulting terms, namely:  dz  ihjiw\z)Z(z) flap  The following procedure avoids this problem. C  can be rewritten as:  m o  He-*  C_„ — t + Z,,) tv'(z) dz m ° ( c i + c 2 ) 2 " ) / V sin/?  (5.46)  •C  since the difference between z and the term replacing it in Equation 5.46 is 0(l/z) as z—>Q (Eq. 5.32) and, by the residue theorem, this will make no difference to the result of the integration. This substitution is useful because t is purely real on both the airfoil and flap. Thus:  'mo  (  C l  +C ) 2  5R 2  Z j w'(z) dz - iH j tw'(z) dz + H cot /? j tw'(z) dz 0  C  But, with H' =  tfe-^/sin/?,  C  (5.47)  C  and since dz/dt = H' +  0(l/t ): 2  j tw'(z) dz = j H'tw'(z) dt. C  C  Also:  ( l - ^ £ ( z ) ~ l + o(j^  ,  a s f ^ o o  (5.48)  dz  (5.49)  so that:  jtw\z)d* = f(f-jji ) S  Staggered Two-Element Thin Airfoil Theory / 5.2  115  Finally, using Equation 5.35, and the fact that the second integral in Equation 5.47 already selects the boundary condition when C collapses to the slit surfaces, one gets: Cl cos /? - x 2  °  m  o  ~  c  1  +  c  0 i  c  2  '-T^*kiz  M {  *  ) d  *  '  (5 50)  where M(z) is purely imaginary on the slit surfaces: M{z) =  Ht{cotl3C{z)-i)  ~*+(tfcot/?A -Z ) + ^ 1  0  z  + ^  z  +  (5.51)  as 2 —•co  and, like Ci, Mi is purely real:  2 sin'/?  f+S + sf -  ( l + S  + / ) 2 4  (5.52)  + < U + *22 - (<l+<2)(<ll+*22)  A recurrence relation for the M 's can be obtained from those for the £ ' s and Z 's. n  n  n  Although dy/dx can be easily broken up into components which are even and odd about the surfaces of the staggered slits: dy _ dy dx dx  c  dy dx ' t  C(z) and M(z) are not so readily separated. Equations 5.45 and 5.50 are probably the simplest forms of these two-element staggered slit versions of Munk's integrals for the overall forces. They could each be evaluated numerically as the sum of four line integrals but this would require as many evaluations (or inversions) of the complicated z(t) function as there were points in the numerical discretization of the integrand. Complex versions of these equations which enabled their evaluation by residues would be highly desirable. For this, it is convenient to again find a complex function, -7^(2), which is single-valued in the dx flowfieldand which: 1) has its real part equal to dy/dx on the slit surfaces. 2) preferably does not combine with £(z) or M{z) to form any simple poles in the z-plane. This condition is not strictly necessary because the residues of any simple poles that are created can be accounted for.  Staggered Two-Element Thin Airfoil Theory / 5.2  UG  Since the oddness or evenness of the boundary condition can no longer be taken advantage of (as with Eqs. 4.17 and 4.18), the overall lift and moment become: (5.53)  mo  (5.54)  dz  (ci+c2)2  The simplest problem is the o-problem: y (z) = -a dx  d  (5.55)  ai  and, using the above equations, one can easily obtain the linearized predictions for the lift curve slope and aerodynamic center of staggered two-element airfoils — results which do not appear to have been previously presented: Q  _ 1  4na  c +c x  d C m o / d a  dC /da  ji  —  (c  l +  2  1  2K' hi l +• ^2 _ 2ira — - n '-i -2 a h c +c TT ^ / ( l - s ) sin/? cos/? _  x  0  1  cos /? + 9± 2  Ci  t  = x + Q  {2K'/n)h  (5.56)  2  c ) = x -£ 2  —  (5.57)  (1+s+f)  2  f+s+sf  , ,2  + tn +1 t n — ti + f 2 2 — t  \/f{l-s) tan/?  22  {ti+t ){t +t ) 2  n  22  hi — h + '22 — h  C( and x ] both reduce to the correct limit when h —* 0, with or without overlap. It is ac  informative to look at the former case in more detail. When h —• 0 for /? ^ 0: 8 —• /, f -* s, and the stagnation point equation (Eq. 5.22) gives (for 0 < /? < TT/2): ^ _> (1 _ cos/?)/2 ,  hi-*f,  h~*f,  *  2 2  - ( l + cos/?)/2.  Also, from the last of Equations 5.23: 2K'  h  x sin/? 22  \  Thin airfoil theory makes no prediction of the vertical position of the aerodynamic center since it predicts each  airfoil element's drag to be zero.  Staggered Two-Element Thin Airfoil Theory / 5.2  since x —• z (l+cos/?)/2 22  c  a n c  117  * %o -* x (l-cos/?)/2 (Eq. 5-29). Thus: 22  C -*2irax /(c +c ) t  22  1  and  2  x  ac  -+ 1/4  as required. (For a single slit, of course, x /(ci+c ) is taken to be 1.) 22  2  When (3 -* TT/2, the indeterminate terms in C( and x  ae  can be replaced with  the RHS of Equation 5.41. Glauert (1926) solves for the exact overall lift on an unstaggered biplane arrangement of flat plates of equal chord length. The coefficient of a in Equation 5.56 reduces to give the coefficient of sin a in Glauert's solution at this extreme. For the r;-problem, the generalized staggered slit flap selector is: r f  ,  2w  K'v'2,  »  t -f  2K'/n  s  ff*2  lt(t-s)(t-l)  2  t-t  V ? ( l - S )  2  v  = 0 + imaginary terms  on the airfoil  = 1 4- imaginary terms  on the flap.  2u 7(oj) is just the  \/t (t -s)(t -f)(t -l)  +  2  2  t-f  2  2  t-t  2  (5.58)  K 0' 1  1- —-^-(w-u/ ) function with an imaginary constant subtracted so 2  that ?{woo) is purely real: ?o + ?Jr+ I T  ,  ~ /o H  ,  +  as*-oo (5-59)  1 ,  *i  •  *2  1—~ H z z  ,  as 2r —• oo.  Using Equation 2.36 and taking advantage of the stagnation point equation, Eq. 5.20, and the relation: r  one gets: , 2JT'  _ 1  «• G  =  (  1+  a +  /_4f ) 8  t -f 2  (  1+s+f  2  V  2  2  ^  1  6\  2  t  2  2  ( i + / -  a  )  2  l n  with #  n  &  (5.61)  / c ^  - / ( l  «  +  Ti and 7 are obtained by replacing £ 2  2K'  08  + / -  a  )  ,  4  T  2  l  i  '  in Equations 5.39. Although T and 7\  are purely real, in general the higher order T^s are not.  Q  Staggered  Two-Element Thin Airfoil Theory / 5.2  118  For the ^-problem, then:  ^W  =  (5-62)  -?W  and Equations 5.53 and 5.54 enable a straightforward solution of the flap effectiveness for a slotted flap airfoil. Again, this appears to be a new result:  Ci+C  2  H  J  Ci+c sin/3  j  4*1 2  t  n  -ti + t  22  -t  2  2cos/3  G -Gi  1  |  2K'(ti-f)(t -f)  2  2  2  cos/3  n  > .  y/fil^s) (5.63)  The last line of Equation 5.63 clearly shows how this expression reduces to the tandem airfoil configuration predictions for dC /dr], Eq. 2.47b. However, when /? —• TT/2, the first t  line must be used in conjunction with Equation 5.20 to remove the indeterminacy in the expression. Just as one would expect, this linearized prediction of dC^/dr] for staggered airfoil elements has the same value regardless of h being positive or negative. The expression for C C  m  «  m  =  mac  is relatively long:  *i » (ci+c r { 7 l g_ c o s ^,,)_^ e -, % j 4  (  ,  6 4 )  2  but it is independent of both E and x . 3  0  The method of solution of the a and r) problems via Equations 5 . 5 3 & 5 4 is more complicated than necessary. If C{z) is written as Z{z\ t ,t ) xx  22  t  then theflowfunctions  that solve these problems are: < ( * ) = ia[l-£(*;*i»<2)] (5.65)  v' (z) = n  ir [7{u)-7 C{z- t t )\ i  Q  )  u  2  and the large z expansion coefficients of these functions give the lift and moment directly. The main reason for developing Equations 5 . 5 3 and 5 . 5 4 is so the forces for the general camber and thickness problems can be easily solved. This requiresfinding-j^-{z) dx  and ^j^{z) functions capable of solving polynomial boundary conditions such as those in dx  Staggered Two-Element Thin Airtoil Theory / 5.3  119  Equations 3.14 and 4.47. No doubt staggered slit versions of the segment selectors of Chapter 4 would play a part in these formulations. However, finding the —~{z) dx  functions for  other than the zeroth order term of a polynomial boundary condition is relatively difficult and  is not pursued further.  5.3  T H E LOCALIZED FORCES  The  airfoil lift is again calculated using Equation 4.28, but with u)' (z) replaced with w'(z)  and  using the staggered slit version of the airfoil selector:  c  A( ) = 1 - _ - _ . - ! ( „ ) + U  hi-s  t{t-f)(t-l)  t-t \  t-3  n  , V'II('II-«)(«II-/)('U-1) = 1 + imaginary terms  on the airfoil  = 0 + imaginary terms  on the  2LJ  flap.  (5.66)  K't?'  This is the 1  j"r( K Vy  IT  u _ w  u ) function less an imaginary constant so that Afa^) is  purely real. Thus: j w'{z) dz = AQ j w'(z) dz- j w\z) C\  C  [X (w) - l]dz - j w'(z)A{u)) dz .  Cy  (5.67)  c%  On the slit surfaces, the unknown pressure terms in w'(z) are completely contained within the imaginary parts of the RHS of this equation (and they must sum to zero), so that the airfoil lift is: = h ± £ l  C t i  A  o  C  l  - *  Cy  If yi(x) and y2{x)  a r e  Cy  u  J  airfoil  upon integrating by parts, c  Cy  j  dx  flap  i L+ 2 A „ c  Q  yy(x),  at least, has a simple zero at z —  then,  Xyy,  becomes:  Ce = - ^-A C l  CiX  continuous functions over their respective domains (including around  the leading and trailing edges) and if  „  (5.68)  £^[A{ )-i]d*-££%lAMdz.  e  +  2  * X  -  , ^A(UJ)  j  ^yy(x)-^dz+-^y2{x)-j-J-dz  airfoil  2i  t flap  . dA{u) s  ,  ,  .  (5.69)  r  Staggered Two-Element Thin Airfoil Theory / 5.3  120  where: 2K'/n dz  n  n  u -  sin/?  dA(u)  ^t(t-3)(t -f)(t -l)  2  + t-t  ^/t[Jt - )(t-f)(t-l) n  2  S  n  2\/t(t — s)(t—f)(t — l) cos/? —tsin  H  (5.70) which is purely imaginary on the slit surfaces and reduces to Equation 4.34 when /? —> 0. If a y(z) function with the necessary characteristics can be found, C  tl  could  be rewritten: C  tl  =  —— A C + 3i j j p j y 0  dz - 27ri(residues of integrand) j  e  (5.71)  and then evaluated by residues. Note that dA{u)Jdz has a pole 0((z—x )~ / ) at the n  3  2  airfoil's trailing edge. Also, for the a and n problems, yi{x) and y {x) are polynomials 2  linear in x. The localized airfoil moment is a particularly difficult problem. u>'(z) cannot simply be replaced with zw'(z) in the above procedure since z is no longer purely real on both slit surfaces. One way of overcoming the problem is to introduce the complementary selector function A(u>) which is just a flap selector function with its singularity at the flap trailing edge. If A(OJ) is written as A{u}\Un), then:  A(u) = 1 - A(u;u ) 22  and  ^(Woo)  = 0 + imaginary terms  on the airfoil  = 1 + imaginary terms  on the flap  is also purely real. Note, from Eq. 5.36, that £(z) oc 1  —  (5.72)  [>l(w)+>l(w)]. Now,  writing A and A for A{u)) and A(UJ), consider:  dtlf[z{A+l-A)  +  ih{A-l)A]w'(z)dz\ J  K  c  (5.73)  = SR j 2 j zw'{z) dz - ih j w'(z) dz\-i £ [x{A - 1 - #) + ih(A - 1)A] ^ Ci  Ci  dz  airfoil  - i £[x{A+l-$) - ihA{l-A)}^}dz flap  .  Staggered Two-Element Thin Airfoil Theory / 5.3  121  The real part of the integral on the LHS of this equation can be evaluated (by residues) in terms of the overall lift and moment (the real part of the second term in the integrand integrates to zero). The real part of the first integral on the RHS gives the airfoil moment, while that of the next integral is zero. The last tv/o integrals on the RHS can be evaluated numerically, with all the problems that entails, but it would be more convenient to convert them to a form that is more easily solved by residues. Until now, this was accomplished by integrating by parts so that the resulting integrands consisted of a boundary condition multiplied by a function whose complex form was purely imaginary on both slit surfaces — a characteristic which greatly simplified the formulation of the complex version of the boundary condition. However, the A A products in the above integrands negate the effectiveness of this process. The problem can be seen in a different way, as suggested by the forms of the mapping function, Eq. 5.15, and the A(OJ) Z = x + 0  H i  r  n *  and A(OJ)  functions. If:  ,- : >/*('-')(«-7)('-l) + >/*ii(«n-«)(*n-/)(*n-l) n n fl  sin/3 i cos p — \ sin p  and, writing Z as Z(t;t ), n  t-hi  (5.74)  if: 2 = Z(t;t ),  (5.75)  22  then one can show that:  z- Z =  ih{A-\)  -  (5.76)  z - Z = -ihA .  Z and 2 are purely real on the surfaces of the slits. The integrands of the last two integrals of Equation 5.73 can now be rewritten: 7 9  2  x{A-l-$)  + ih{A-l)A  = x{A-l)  - ZA =  ~  X  2  x{A+l-A)  - ihA{l-A)  = xA + Z{l-A)  =  x--Z2  (5-77)  ih  The Z2 product has only square-root poles at the airfoil and flap trailing edges. The RHS's of Equations 5.77 probably should not be used for numerical calculations involving  Staggered Two-Element Thin Airtoil Theory /  5.3  122  small gap configurations since they become indeterminate as h —• 0. Nevertheless, they are informative. One way of writing the airfoil moment, then, is: ( c y + c m i  o  \  2  Cy  \  +  2  ( C i + C  2  J  M  °  V  airfoil  2  c\  \ )  Ay-fit  2  1  flap  The Z f - ^ integrands can be handled by converting to one complex integral: ax  -ji»yzzV(*)<fe.j  SR<J-»|> ZZw'{z)dz\.  (5.79)  The difficulty is with the remaining x — terms. They show that the problem of turning dx the integrals of Equation 5.78 into complex integrals i3 the same one encountered in the previous section. That is, since z  n  is not purely real on both slit surfaces, finding the  complex versions of polynomial boundary conditions requires the development of new procedures and/or functions. No doubt there is a solution to this problem, but it is left to be dealt with in some later work.  CHAPTER 6 Conclusions and Recommendations  The tandem thin airfoil theory of Chapters 2, 3, and 4 has shown how successful a linearization of the multi-element airfoil flow field can be. The theory's predictions of the aerodynamic characteristics are all analytical expressions, some of which are simple algebraic functions of the initial geometry while the rest are calculable on a small hand held computer. It has been shown that the consistent relationship between the exact potential and linearized theories' overall force predictions extends, at least, into this region of the multi-element regime. This, in conjunction with the arguments of Chapter 1, suggests that a thin airfoil theory analysis of a properly designed multi-element airfoil configuration would give overall force predictions which correlate with the real flow forces as well as and probably better than is the case in one-element theory. In addition, for a + r\ values less than 20° or 30°, these predictions should be better than those of exact potential flow theory.  i As mentioned in Chapter 1, however, the breakdown of the correlation be-  tween thin airfoil theory and the realflowcan be abrupt and come without warning since, of course, the linearized theory does not model the boundary layer effects that are responsible for this breakdown. In addition, the user must have other input in order to judge just what constitutes a "properly designed configuration". While not insignificant, these 123  Conclusions and Recommendations /  6  124  disadvantages are the same ones with which the successful exact potential and one-element linearized theories must also contend. The above expectations of the correlation between the multi-element linearized and real flow results have yet to be fully proved. This will require calculations of the overall force predictions for the staggered airfoil elements of Chapter 5 in conjunction with a solution of the staggered slit mapping. Comparisons with experiments on realistic slotted flap configurations, such as presented by Foster et al (1970), should then be possible. It is unfortunate that the simplicity of the tandem slit mapping is lost when the more realistic staggered slit geometry is accounted for. However, this author does not see this as a major problem and suggests that a form of Newton's method could be successfully applied to Equation 5.15 to obtain the mapping parameters. The assembly language programming capabilities of a small computer such as the TI CC-40 might be useful for achieving the necessary reduction in computation time. Given a solution to the mapping, the staggered two-element thin airfoil theory force predictions can be evaluated from the expressions given in Chapter 5, with the overall forces for the problems of incidence and flap deflection having simple analytical solutions. The general effects of mean line curvature and thickness on the overall forces are summarized in the integrals of Equations 5.45 and 5.50. These equations are probably as efficient a way as any of handling numerically specified boundary conditions. The analogous expression for the localized airfoil lift is given in Equation 5.69.  The airfoil  moment integral (Eq. 5.78 with the integrands replaced with the LHS's of Eqs. 5.77 for small gaps) can certainly be used as presented here, but it would be preferable to put it in a form where the non-boundary condition factor in the integrand was at least expressable as an algebraic function of the f-variable—for any gap size. It would also be desirable to obtain analytical solutions for the incidence and flap deflection problems from each of the localized force integrals.  Conclusions and Recommendations  /  6  125  A time constraint has prevented the completion of the staggered slit problem to the degree desired. However, it is felt that there are no major obstacles to this completion and that, if done, the results might very well suggest a simplification of the general theory as presented in Chapter 5.  With modern computational hardware, engineers today tend to concentrate on the numerical approach, using powerful computers to get quick numerical answers. The approach used in-this thesis has been analytical, making extensive use of conventional elliptic functions and their supporting documentation. This has led to direct expressions for a comprehensive set of doubly connectedflowfieldcharacteristics, as well as very efficient procedures for evaluating some of these characteristics. In this case, at least, using modern computational hardware together with powerful analytical techniques has provided readily accessible numerical answers, as well as a better overall understanding of the problem.  126  References Tests in the Variable-Density Wind Tunnel of the 23012 Airfoil with Plain and Split Flaps. N A C A Report No. 661  Abbott, I.H. and Greenberg, H . (1939)  N.A.C.A.  Abbott, I.H. and von Doenhoff, A . E . (1959)  Theory of Wing Sections. Dover Publications, New  York  Die Tragende Wirbelfache als Hilfsmittel zur Behandlung des ebenen Problems der Trag&ugeltheorie. ZAMM, vol. 3, p. 290  Birnbaum, W. (1923)  Foster, D.N., Irwin, H.P.A.H., and Williams, B.R. (1970)  Slotted Flap.  Garrick, I.E. (1936)  The Two-Dimensional Flow Around a  ARC RM 3681  Potential Flow About Arbitrary Biplane Wing Sections. N A C A Report  No. 542 Glauert, H . (1926)  The Elements of Aerofoil and Airscrew Theory.  Cambridge University  Press  Potential Flow Analysis of Multielement Airfoils Using Conformal Mapping. AIAA Journal, Vol. 17, No. 12, p. 1281-8  Halsey, N.D. (1979)  Hess, J.L. and Smith, A.M.O. (1967) Calculation of Potential Flow Progress in Aeronautical Sciences, Vol. 8, Pergamon Press  About Arbitrary Bodies.  A Modern Look at Conformal Mapping Including Multiply Connected Regions. AIAA Journal, Vol. 14, No. 8, p. 100G-11  Ives, D . C . (1976)  „  Kennedy, J.L. (1977) The Alberta  •  Design and Analysis of Airfoil Sections. Ph.D. Thesis, University of  An Extended Theory of Thin Airfoils and Its Application to the Biplane Problem. N A C A Report No. 362  Millikan, C B . (1930) Munk, M . M . (1922a)  General Theory of Thin Wing Sections. N A C A T.R. No. 142  Munk, M . M . (1922b)  General Biplane Theory. N A C A Report No. 151  Richardson, J.R. (1981)  Aerodynamic Forces on Perforated Bridge Decks. National Maritime  Institute, England, R 118 Stewart, H.J. (1942)  A Simplified Two-Dimensional Theory of Thin Airfoils.  Journal of  Aeronautical Sciences, Vol. 9, No. 12, p. 452-456 Theodorsen, T. (1931)  Theory of Wing Sections of Arbitrary Shape. N A C A Report No. 411  Whittaker, E . T . and Watson, G.N. (1927)  A Course of Modern Analysis.  4  t h  edition, Cambridge  University Press  An Exact Test Case for the Plane Potential Flow About Two Adjacent Lifting Aerofoils. ARC R M No. 3717  Williams, B.R. (1971) Woods, L . C . (1961)  The Theory of Subsonic Plane Flow.  Cambridge University Press  APPENDIX A Exact-Numerical Potential Flow Computations  This appendix describes the numerical procedure used in this thesis to obtain predictions of the exact, two-dimensional, incompressible, potential flow about multi-element airfoil configurations. The procedure is a slightly modified version of one of the more efficient surface singularity methods, by Kennedy (1977).  A.l  T H E FLOW FIELD F O R M U L A T I O N  In Kennedy's method, airfoil surfaces are represented by a series of straight line segments over which vorticity is distributed. The ends of these segments are the specified coordinates of the airfoil while the middle of each segment is a control point at which the tangential flow boundary condition is applied. This boundary condition can be accurately satisfied by using the stream function formulation for the flow and requiring the value of the stream function at each control point of the A; airfoil element to be the same (rp say). The th  k  Kutta condition is satisfied by requiring this V"jt streamline to also go through a point just off the trailing edge, on the bisector of the trailing edge angle. The general stream function formulation for an arbitrary point (x, y) in the flowfieldaround some multi-element airfoil configuration is: ^(x, y) = ycosa — xsina +  / l{S) ln(r(x, y\S)) dS .  27T J  (Al)  S  127  Exact-Numerical Potential Flow Computations /  A.1  128  Here, theflowat infinity is at an angle a to the z-axis; S is all the surface in theflowover which the dimensionless vorticity, 7(5), is distributed (the equation has units of length, each side having been divided by the magnitude of the velocity at infinity); and r(z, y;S) is the distance from a point on S to [x,y). This integral equation is solved numerically by dividing surfaces S into N straight line segments and then assuming that 7(5) is constant over each segment. Thus: 1  f  N  V'K y) » y cos ct - z sin a + —  7y / In (r(z, y; 5 )) dSj  (A2)  y  S  }  so that at the i  t h  control point on the fc airfoil element, the boundary condition becomes: th  N ip  k  = y,- cos a —  Xf  sin a + ^  IjKij  >  t = l,2, . . . , i V  (A3)  3=1  where the Jf.-y's are the influence coefficients and are calculated from the geometry of the airfoil configuration. They are given in Kennedy's thesis. In this problem, there are N+M  unknowns: N 7y's and M V'jt's, where there  are M airfoil elements in the configuration. The additional M equations are the Kutta conditions. They are the same as Equations A3 but with (z,-, y,) replaced with (z , y ), the tp  coordinates of the trailing points. The N + M equations are linear in the N+M  tp  unknowns  and so can be solved through a solution of an (N + M) x(N + M) matrix. This formulation has the added advantage that the matrix solution gives the 7y'8 directly and, since it can be shown that the velocity of the flow inside the boundaries of the airfoil elements is everywhere zero, 7(5) is also the velocity distribution over the surfaces of the airfoil elements. Pressures are then easily calculated and, as Kennedy shows, if C . = 1 - |7y| is taken to be the pressure at the control point (z -,yy), then it p  2  ;  will closely predict the true value of the pressure at that point if each airfoil element is represented by a sufficient number of segments (about 30 according to Kennedy). Also, since the sign of 7y gives the direction of the velocity over the j  th  segment, it is simple to  find the stagnation points and verify that the Kutta condition is working.  Exact-Numerical Potential Flow Computations /  A.2  129  The approach used in this thesis differs from the above approach only in the way in which the Kutta condition is applied. Since airfoil surfaces are discontinuous at the trailing edge, it is inevitable that the segmentation of the boundary results in the trailing edge being the junction of two straight line segments.  Hence, it is not a  control point and, therefore, it is not redundant to specify as the Kutta condition that the ipk streamline go exactly through the trailing edge. This simplification eliminates the need to calculate trailing point locations and gives essentially the same numerical answers as Kennedy's method. Care must be taken, however, in evaluating 2M of the Jf,-y influence coefficients since each of the M trailing edge control points is exactly on the ends of two of the Sj segments, resulting in OlnO and Otan (0/0) terms in the expression for -1  This trailing edge Kutta condition has been found to work well for slotted flap airfoil configurations at angles of attack up to a = 90°. In one instance at this high angle of attack, it was found that the flow around the deflected flap had actually reversed, with the ipf streamline first encountering the flap at the trailing edge and leaving it somewhere on its upper surface. Despite this, the circulation determining condition, that requires the ipli  streamline to branch into or out of the body contour at the trailing edge, still worked.  A.2  T H E FORCES  The remainder of this appendix deals with the overall and localized force calculations. Obtaining the overall forces is a reasonably straightforward procedure. The Kutta-Joukowsky law, for example, gives the overall lift in terms of the total circulation around the configuration, which is just the sum of all 7y£y products (£y being the length of segment j). Kennedy only needed to calculate the overall lift and did so in this manner. The localized forces are more difficult to calculate. They are usually obtained through a numerical integration of the surface pressures including, perhaps, a curve fit of the pressure distribution to improve accuracy.  Exact-Numerical  Potential Flow Computations  /  A.2  130  In this thesis, the forces are all calculated the same way, by using the Blasius equations to solve for the exact expressions for the forces in terms of the 7y's, so that no error is introduced beyond that which already exists in the approximation of Equation A2. To begin, the exact formulation for the flow, Eq. A l , is rewritten in its complex form. If the stream function is combined with its harmonic conjugate, the velocity potential, to form the complex potential fl(z) =  0(«) = ze~ + ± J  <f> + iip  t  then:  7(5) ln(* - z ) dS  ia  (A4)  s  s where z is on surface 5 in the complex plane z — x+iy. The segmentation approximation s  is:  j 7(5)ln(* - z )dS*Yl7y  Jln(z - z .)dSj  s  S  i=  (A5)  s  l  Sj  and, using the notation of Figure 20, the complex potential is:  0(«) « ze~ + JL £ Tj j Hz - z ) d£ . J o  (A6)  ia  (  Note that this approximation does not necessarily require 7,- to be constant over all of the j  t h  segment. The mean value theorem, for example, allows Equation A6 to be exact for a  given z. It is by dealing directly with Equation A6 and avoiding unnecessary assumptions about the physical meaning of 7y that the expression's accuracy is fully retained in the following force calculations. Performing the integration yields:  0(*) = ze~ + JL £ 7y[(*-*i>-''> ln(*-*i y ) ia  -  )«-''> l n ( * - « * , ) -  (,•].  (A7)  J The real part of the square-bracketed expression gives the matrix influence coefficients, the if,y's, when z = z,-. The dimensionless complex velocity is:  "(*) = ^  = e-'  a  +^ E V " y  J  "  (A8)  Exact-Numerical Potential Flow Computations /  A.2  131  z - PLANE  2  U  The control point z - in the middle of the j segment  FIGURE 20  th  }  The Blasius equations for lift and moment are: D-iL=^\w \ e 0o  2  fw {z)dz  ia  (A9)  2  M = £\tv \ Vl{|yj>zw zw (z)dz (z) dz^ 0  QO  2  (A10)  22  where the drag is taken to be parallel and the lift perpendicular to the flow at infinity. These expressions can be evaluated by residues using the expansion, as z —» oo: real terms  +  (All)  Thus: (A12)  c =o A  (A13)  C, = 1  (ci + -  Z  •+c ) r M  (ci+'-Hw) ^ ' ^ V - J — 2  •  (AH)  Calculating the localized forces on the k^ airfoil element is less straightforward. For the lift and drag: - j w {z) dz 2  Ck  (A15)  Exact-Numerical Potential Flow Computations / A.2  132  This expression is broken up into two integrals, as suggested by the form of Equation A 8 in the integrand. Thefirstintegral is evaluated by residues: 7  i  = / Ev""'M - ii) z  2  -  ln(*-J&,)]<fa = 2 i r t £ ;  7  /  y  .  (A16)  This is correct since, for j £ k there are no singularities within C/. and so no contribuf  tions to the integral; for j € k,  can expand out to infinity, allowing the expansion of  Equation A l l to be used in an application of the residue theorem. The remaining integral is:  / = /EEw2  i , w z-zl}. In  In  z-z  X}  (A17)  dz.  It is convenient to break up the double summation as follows:  £X>££+E£+££+££ • j  J  jek Jek  jekHk  &k Jek  (A18)  j& J$k  The last double summation term on the RHS of this equation does not contribute to I2 since it has no singularities within the Ck contour. Neither does the first term make a contribution since: Jreal terms \  E£~°( \  jek Jek  z*  as z —• oo  J  The middle two double summation terms are equal. Thus, I2 can be rewritten  J  2  =EE^/ 2  ftk Jek  G  WN*-*ii) -  M*-*2,)]  dz  (A19)  l  k  where G(z) is regular within and on Ck but has singularities elsewhere which prevent Ck from expanding out to infinity. The integration is carried out by joining z  lf  and z  2j  with  a branch cut and integrating around it. One gets: *2j  j  v7(s)[ln(*-* ,) - m{z-z )] dz = 1(  2]  2JT» j  G{z) dz .  (A20)  This expression can now be evaluated conventionally: *2,  2fft<T'(«/+M j [ln(z-z ) - ln{z-z .)]dz = 2iri(Zj - - Zjj) ly  2  3  (A21)  Exact-Numerical  Potential  Flow  Computations  /  A.2  133  where:  Zjj = Xjj - iYjj = e -i(Oj+°j) z  2j  In  -  In-  Z l  -±  (A22)  Thus, the forces are: C  *. = ^ I1IJ\EV[™«( JJ- JJ) X  -8ina(y  X  L  ^  { 2£j - ^ £ 7y [cos a (y K j&k  = f E Jek k  jy  j y  -F  - Y ) + sin a ( X }J  (A23)  )]}  y j  j y  - JT )]}. (A24) J yj  The localized moment is calculated in the same manner. If:  W j ^ U j j - i V j ^ e - ^ ' * )  2^  In  (A25)  - zf, In z  lj- 2 Z  }  then: C.,.,  (A26)  —  The localized force "influence coefficients", the Zjj's  and Wjj'a,  are most  easily calculated simultaneously with the matrix influence coefficients, the iif.y's, since they are all calculated from common geometrical parameters. The efficiency of this calculation is improved if one takes advantage of the fact that, if:  HJJ  =  In ^  and  (A27)  ujj = z In Xj  then: Zjj  ~ ZjJ =  ~ Pjj ~ A*y+i,J + h j ] ~ e  i [ 0 3 + 9 j )  (A28)  and: (A29) Note the J *-» j symmetry in these equations. This is why, when all the localized force coefficients are added together to form an overall coefficient, these terms all sum to zero.  Exact-Numerical Potential Flow Computations / A.3  A.3  134  COMPUTATIONAL ACCURACY  The numerical calculations are performed on UBC's Amdahl 470 mainframe computer. All calculations are made in double precision (numbers are represented with 16 significant figures). An indication of the accuracy of the present method (for which forces are calculated using Eqs. A12, 13, 14, 23, 24, and 26) is obtained by comparing its overall lift predictions to Williams' (1971) exact solutions for the potential flow around a slotted flap airfoil configuration. This is done in Table 3, where the numerical predictions are made by 1) using all the coordinates specified by Williams (61 segments per airfoil element) and 2) only using every other of the coordinates (31 segments per airfoil element).  Ci  a = 0°  a = 10°  TABLE 3  error  Williams, exact  2.7230  61 segments/element  2.7180  0.18%  31 segments/element  2.7240  0.04%  Williams, exact  3.7440  61 segments/element  3.7389  0.14%  31 segments/element  3.7434  0.02%  Comparisons of the exact and present numerical predictions of the lift on Williams' conBguration A slotted Bap arrangement. «  The accuracy of the present method is seen to be quite good. However, its predictions do not appear to be converging to the exact values as the number of segments per airfoil element (N say) increases. An explanation for this anomaly is found in the k  following investigation. Figures 21 to 25 show in greater detail how the numerical force predictions on a slotted flap airfoil configuration vary with N . The configuration consists of NACA 23012 k  airfoil elements in an arrangement similar to that of Williams' configuration A, with the  Exact-Numerical Potential Flow Computations / A.3  135  airfoil's chord line stretching between 0 and 1 on the z-axis, the leading edge of the flap located at (.99,-.02), the flap chord length equal to 0.373, and the flap deflected 30°. a = 0° for these calculations. The analytical expressions for the 23012 profiles given in Abbott and von Doenhoff (1959) make it easy to segment the airfoil boundaries at will. Throughout this thesis this segmentation is done in the same way, by dividing the chord line up according to equal increments in 6 (Eq. 1.14) and then finding the profile coordinates opposite those divisions. This provides a concentration of smaller segments around the leading and trailing edges of the airfoil element boundaries. Figures 21 to 25 show that the force predictions of the present theory do indeed converge to a constant value as N becomes very large but that, initially, the k  predictions have a tendency to oscillate. This is the probable explanation for what occurs in Table 3. That is, the good agreement between the numerical and exact lifts at N = 31 k  is only coincidental and, as can be seen from Figure 21a, cannot be counted on to repeat itself in other situations or with other profiles. Note that in the interpretation of Figures 21 to 25, the not unreasonable assumption is made that the asymptotes to which the forces are converging are identical to the values of the exact forces. Thesefiguresalso show the force predictions of a simpler but cruder analysis in which all the forces are calculated by integrating the pressure distributions and assuming C . = 1 — |7y| is constant over all of segment j. These force predictions appear to converge 2  p  to the values of the present theory as N becomes large but are generally less accurate for k  a given N , especially for the drag calculations. k  All of the "present method" force predictions for N > 60 in Figures 21 to 25 k  are within 0.2% of their values at N = 100. Throughout the thesis, this is the method k  and degree of segmentation used to calculate exact-numerical potential flow forces. The accuracy is needed for making accurate predictions of the variation of these forces with small increments of a and n.  ExactrNumeric&l Potential Flow Computations /  A.3  136  2.22  C,  0%  1  1  •  ••• 2.20  -1% 2.18  _1  2.16  20  L_ 60  40  80  100  a) Overall lift vs. Nf..  0.08 •  present t h e o r y  •  Cj P  a s s u m e d c o n s t a n t o v e r a l l o f s e g m e n t /.  0.06  0.04  0.02  •  20  i  T  40  60  •  *• 80  100  b) Overall drag vs. N . k  F I G U R E 21  The accuracy of the overall lift and drag exact potential flow numerical calculations as a function of the number of segments per airfoil element. NACA 23012 airfoil elements in a slotted Bap conGguration at a = 0°.  Exact-Numerical  Potential  Flow  Computations  /  A.S  137  2.44  0%  "1  •  1  2.42  2.40  20  40  80  60  a) Airfoil lift vs. N  1.64  100  k  +1%  -/2  S  •i  1  1  1  1.62  -1% 1.60  1.58  present t h e o r y Cj  1.56  P  assumed constant over all o f segment /  1.54  0  20  40  60  80  100  b) Flap lift vs. N  k  F I G U R E 22  The accuracy of the localized lift numerical calculations. NACA 23012 airfoil elements in a slotted flap configuration at a = 0°.  Exact-Numerical Potential Flow Computations /  A.3  138  -0.26 -1%  •  +1%  " - -  -  w  • •  -  w •  z  1  1  -0.24  •0.22 • present theory • C j assumed constant over all of segment /. P  -0.20  0  20  40  60  80  100  a) Airfoil drag vs. N . k  0.72  0.70  0.68+1% •  -I  m  1  1  •  •  •1% 0.66  -• 20  40  60  80  100  b) Flap drag vs. N . k  FIGURE 23  The accuracy of the localized drag numerical calculations. NACA 23012 airfoil elements in a slotted Sap configuration at a = 0°.  Exact-Numerical Potential Flow Computations / A.3  •  139  present theory  • C j assumed constant over all of segment / P  -0.92  -0.90 -1%  -*  •0.88  "I  »-  1"  +1% -0.86 20  F I G U R E 24  40  60  80  100  The accuracy of the overall leading edge moment numerical calculations. NACA 23012 airfoil elements in a slotted Bap conBguration ata = 0°.  Exact-Numerical Potential Flow Computations /  •1.02  A.3  140  •  present theory.  •  C j assumed constant over all of segment /. p  -1.00 -m 10  •1%  -0.98  1  1  1  *  +1% •0.96 20  40  60  80  100  80  100  Nk  a) Airfoil moment vs. iVj..  -0.58  •0.56  -  •  -m 20  -  -0.54  •  •  •  1  +1% _L 60  _L  -0.52 20  40 N  k  b) Flap moment vs. N^.  F I G U R E 25  The accuracy of the localized leading edge moment numerical calculations. NACA 23012 airfoil elements in a slotted flap conBguration at a = 0°.  APPENDIX B Elliptic Functions  In this appendix, the reader is presented with the various properties of elliptic functions needed for the present theory. In the first section, some of these properties are summarized from Chapters 21 and 22 of Whittaker and Watson (1927) (often referred to as "WW" followed by a section number) and they are then used to develop further simple relationships. In the remaining sections, more complicated relationships are developed for the evaluation of the two-element thin airfoil theory lift and pitching moment coefficients.  B.l  GENERAL  In the tandem airfoil mapping, Eq. 2.3, elliptic integrals of the first kind, which can be written in terms of an inverse Jacobian elliptic function, are encountered. The Jacobian elliptic functions, sn A, cn A, and dn A, are doubly periodic with respect to their argument A, leading to cells of periodicity in the A-plane. These periods are described in terms of K and K', the complete elliptic integrals of the first kind (;i.e. sn A = sn(A ±4K±  UK')).  K  and K' are determined via their moduli, k and k', for which numerical values are known in the tandem airfoil theory. The parameters of the rectangle, Fig. 4c, can all be calculated by using either conventional series in k or A;', tables, or Landen's transformation which is a repetitive way of reducing an elliptic integral to an elementary function. None of these procedures, however, are as efficient as the ones described in WW and those developed in this appendix. Ml  Elliptic Functions /  B.l  142  When numerical results are required, the elliptic functions to use are the four theta functions: t? (ur, g) = 2q« sin u — 2q* sin 3w + 2q~* sin 5w 1  i  \  1  9  21  V2[u,q) = 2q* cos a; + 2q* cos 3a; + 2g < cos5a; H  (Bl)  ^3( » q) = l + 2q cos 2a; + 2g cos 4a; + 2g cos 6a; H 4  w  9  t? (a;, g) = 1 - 2q cos 2a; - f 2q* cos 4a; - 2g cos 6a; H 9  4  where, in keeping with most of the notation used in WW: .K'  _*KL  TTA  It is interesting that: = i? (2oM ) + t? (2o;, ) 4  3  2  g  4  v {u q) = # {2uj,q*)-v {2u,q ). i  3  i  (B2)  4  2  The theta functions, while periodic in ir or 2TT in the ^-direction, are not perfectly periodic in the ^-direction and so are referred to as quasi-doubly periodic functions. If i?'(o;) is the derivative of t? with respect to a;, then for any of the theta functions these periodicities give (WW21.11): 0' 0' -(a;+*r) = -(a/)  tf' ;  -{U+KT)  =  where t?(a;) is written for 0(u q) and "^(w) for  -2i + -(w)  (B3)  . However, since theflowfunctions  t  must be periodic in JTT, Eq. 2.18, the following function, which has this required periodicity, is of interest:  7  M +5 i T  (W)  Alternatively, if the theta functions are defined relative to the a/ variable, Fig. 4e, and the q' parameter, then: j(u'+ir,q') =  (B5)  Elliptic  Functions  /  B.l  143  is also a possible flow function. In fact, Jacobi's imaginary transformation (WW21.51) gives:  ,  „  2u  K'fi'  ,  ( B 6 )  0' The -r^-(w) function has a simple pole at the origin of the w-plane which can 01  be transferred to the other corners of the rectangle using the half-period properties of the theta functions (WW21.11):  Now, evaluating the theta functions at their zeros (WW21.61), gives:  and:  /ivic\W  /oirxVa  V2  0i(O) = 0 (O)0 (O)«/ (O). 2  3  (B9)  4  Using Equations BI and B8, then, K can be related to the q parameter through an infinite series. Also, since y/k = t/ (0)/t7 (0) and using Equations B2, an expression for q in terms 1  4  3  of a known parameter, e , is developed (WW21.8): 0  9, = ~ ~' l  - * » ( o . g ) - * 4 ( Q , q ) _ 02(0,g )  y/k  4  " - l + v/it  U^)  7  +  U0 q) t  4  «/ (0,? )' 3  m  1  i  m  '  This infinite series for e in powers of q can be inverted to give: 0  q = c + 2e + lSe + 150e£ + O (lOOOeJ ) . 0  5  0  9  0  3  7  (Bll)  Thus, q is calculated using Equation B l l , K is obtained as above, and: K' = -—lnq.  (B12)  Elliptic Functions  /  B.l  144  For typically small values of f—s, these are all very rapidly converging series. As f — s —* 0, k —> 0, k' —> 1, and q —» 0; for most configurations q < 0.1. However, even if A: is as large as 0.999 (so that q = 0.33), these series will give four significant figures. Note that for small slots, y/k will be very close to 1 and so the number of significantfiguresin 1  a calculation will be reduced when evaluating Equation B10. This problem is avoided by using Equations 2.48a &b. The Jacobian elliptic functions can be written in terms of theta functions (WW22.11):  °^ =OT ^=yfS ™=^£R;  Also:  d  ;  d  (B13)  d  ,  — sn A = cn A dn A ; — cn A = — sn A dn A ; — dn A = —k sn A cn A oA dX dX  (B14)  so that, taking the derivative of the logarithm of Equations B13: cn A dn A sn A  2K  sn A dn A cn A  (B15)  IK IT —[to]  •  2K dnA A suitably chosen combination of quasi-doubly periodic functions, then, will yield a doubly  periodic function. Also, taking the derivative of a quasi-doubly periodic function can do the same. From WW22.73, Equations B15, and using Legendre's relation: (B16)  EK' + E'K - KK' = -  one can show that: d  2ui  2E'  —  +  —  +  TT sn A 2  dX d  2E'  dX  _ TC  d  '2ui  d dX  +  TT  - % )  '2ui  7T  2  2K'dn X 7r cn 2 A  2E'  —  dX  1  2K'  2  v 2E'  K^ \  n  dn A 2  k sn A  + U  2  2K'k cn X  TT  (B17)  Elliptic Functions / B.2  145  Here, E and E' are complete elliptic integrals of the second kind and, as with K and K\ they are the same function of their respective moduli k and k': l  E  o The easiest way to evaluate E is through one of its three q series representations. Using versions of the last three of Equations B17, Eqs. B l and B8, and evaluating the expressions at u = 0 gives:  I (2T) I= ( =  V2  7=k  (  2?1/4  ) " (1+  2  +  2 5 , 6 +  4  "  1 2 +  "  0  (  B  1  9  )  8  ^  +  2^)  V 2  (  *  + V  " + ^ +** ) '  +  9  1 6  C  E' is then obtained through Legendre's relation.  B.2  E L L I P T I C F U N C T I O N SERIES R E P R E S E N T A T I O N S IN P O W E R S OF q 4  In Chapter 2, Eq3. 2.47 & 48, the a and rj problems are solved in terms of the three constants Gi, G2, and E3. G and G are simple functions of K' and E' and are evaluated x  2  with the series representations of section B . l . E , however: 3  _2  Xco  ,K'0'  requires evaluating both an elliptic integral of the first kind (Eq. 2.11) to obtain x<x» ^ an  0' then the  -^(Xoo)  ' function, which is similar to an elliptic integral of the second kind. This  03 can be done conventionally, using Landen's transformation to get Xoo series for  -^-(Xoo)  a n c  * *  n e n  Fourier  (WW21, chapter problems), or the following more efficient procedures  can be used. Inverting a Jacobian elliptic function, which is equivalent to evaluating an elliptic integral of the first kind, can be accomplished by exploiting the theta function  Elliptic Functions f B.2  146  representation for dnA, Eq. B13, in conjunction with Equations B2. Thefirststep is to obtain dnA, by using Equations 2.10 if necessary. Then: _  2 e € w  =  dnX-Vk  1  _ 0 (2u;,g ) 4  2  dnA + V F " 0 (2u/,g ) 4  3  2q cos 2w + 2g cos 6w + 2g cos lOw H 9  25  1 + 2g cos 4u/ + 2g cos 8w + 4  16  or: cos 2w = — 1 + 2g (2 cos 2u>-1) - g (4 cos 2w-3) 9 4  2  8  ]  2  By continually resubstituting this equation into itself, a series for cos 2w in powers of g and in terms of:  £ = ^  4  (B21)  is obtained: cos2w = tf[l - g ( 2 - 4 £ ) + g ( 3 - 2 0 £ + 3 2 £ ) - g ( 6 - 7 6 £ + 2 7 2 £ - 3 2 0 £ ) + 0 ( g ) ] . 4  2  8  2  4  1 2  2  4  6  16  (B22) Note that u can be real or imaginary. For u real, it can be shown that |£| < 1 and, therefore, that the coefficients of powers of g in Equation B22 are no worse than the 4  coefficients of powers of CQ in the expression for g/c , from Eq. B l l . 0  To solve for u = x<x»  o  n  e  g * e  s:  1  B23  *«» =r=—  VF  ( > + 1  since dn fi^ — s/f. Using Equations B21&22 completes the procedure. 0' Evaluating the -^{v) functions is easily done using their Fourier series representations. For example: n=l  W  (B24)  which converges in powers of g. However, since the Jacobian elliptic functions of Xoo 0' are known (Eqs. 2.9),  -^(Xoo)  03  c  a  n D e  evaluated with a series which converges in powers  Elliptic Functions  /  B.3  147  I?'  of q . Although the following procedure is applied to the -r^(w) function, any of the other 4  v  3  —(u>) functions can be included by first converting them to the -r^(w) function using Equations B15. Equations B2 can be rewritten as : i/ (o;,q) + 0 (w, q) = 2d {2u, q ). 3  A  4  z  Taking the derivative of the logarithm of this expression and rearranging gives:  t? (w,g) 4  or:  2/f  , snAcnA Ar  =  d n X -  +  ^  2  ^  '  C > B2B  For E , then: 3  2  Xoo  2*'  ,  + tf ^  K'^{-Yq^smAnx  8  f T ^  •  (B26)  n=l  Again, u can be complex in these equations but this would reduce the rate of convergence of the series representation for -r -(2w,g ). This problem is overcome in the next section. 3  v  B.3  4  3  E L L I P T I C F U N C T I O N SERIES R E P R E S E N T A T I O N S IN P O W E R S O F q 16  In Chapter A and Appendix C,first,second, and third type elliptic functions with imaginary arguments are inverted and evaluated. While Equations B22 & 25 can be used for the first two of these computationsif f—s is small enough, they cannot be used over as large a range of f—s values as is required in this thesis. In this section, the procedures of section B.2 (that gave the q* series) are applied twice in succession to get series which converge in  Elliptic Functions /  powers of q . 16  B.3  148  This process begins with the derivation of a more efficient version of the  series representations for q Eq. B l l . t  In the same way that k determines q, Landen's transformation (WW22.42) shows that: 1ki =  k'  1-  k =  —  determines  q  k'  determines  g  JJ 1+ k  2  2  4  x  etc. In fact:  * —  ko  4  (l + A ^ l + v/* ) 7  Now, if e  02  4  ( 2 e  °  ) 4  '  (B27)  is defined: _l-V^ 0 (O,g  2,  =  2  1 6  )  then, from Eqs. BIO and B l l :  9  4  = e +2€g + 15eg + 150e 5 + .-. 02  so that:  2  2  1 , 1^4  1/4  1  +  2 %  0  117^8 , n(„M\] +0{q  £  -l2- 0  2  (B28)  where: 2jo  1/4 °2  ~  r  .  — x 2 2T l / 4  . . . .  (B29)  [2(1 + k' ) (1 + y / % ) ] 2  As before, to invert a Jacobian elliptic function one first obtains dnA; 2e is then given by Equation B20. However, since: w  0 {u,q) 2  t/ (0,<7)cn(A,A:) t/ (0,9)dn(A,A:) ' 2  3  2c can be rewritten as: w  2e = 2e w  0  cn(A ,fc ) 2  2  dn(A ,fc2) 2  which, when squared and rearranged gives: dn (A ,A: ) 2  2  L/2  2  1-  (4*o02  (B30)  Elliptic  /  Functions  B.3  149  Now, in the spirit of Equation B20: 2e  dn(A ,A; ) - \A£ 2  =  w a  2  "dn(A A ) + / ^ 2 >  2  >  so that, as in Equation B22: cos 4w = However, c  W2  1-? |2  L  w  16  2  +  and g both have terms in them O(CQ) which should be cancelled before their 4  ratio is calculated, g is given by Equations B28 & 29 and 2e 4  W2  is rearranged as follows:  dn (A , A ) - k' 2e,.,„ = " [dn(A ,k ) + y/lcl] [dn (A ,k ) + k' ] 4  2  2  2  2  2  2  2  2  2  2  2  with: dn (A , A ) 4  2  2  2  * = ^M(2e 2  0  )  4  \%  g^)' - 1 -  (2 J £  4  so that if: (1 + A )(1 2  +  V  ^) [2(^) -1-(26 2  2  W  )  4  c=  (B31)  2  1+  \ ,^)J ^)J L dn (A f ' dn(A dn(A ,,* )J 2  +  2 >2  22  2  2  then: 2 f? = £ 2 [ l - 2 g 1 6 + 5g 3 2 + 0(g 4 8 )] 2  9  and: cos4UJ = ( {l 2  +  4g [el-l] + 4 g 16  32  An investigation of the characteristics of £  [ 3 - l l £ + 8 £ ] + O (lOO^g )}. 2  2  m  2  4  48  (B32)  the rectangle of Figure 4d reveals that the  only very active term in it is the expression:  GO 2  When u = |,  1( 2 0  4  is zero. For w real, |e | < e and |£ | < 1. However, from Equation B20 it w  0  2  can be seen that c has a pole at u = | ± 27rr, so that in the rectangle (—^ < 17 < T^) it w  Elliptic Functions /  B.3  150  will take on its maximum magnitudes along u = x± from the pole. At w values of  although this is still some distance  | ± ^ , and § ± y-, |ej = 1 so that £ =  0{l~2)-  2  Thus the maximum error in Equation B32 for u> anywhere in the rectangle is O(l00g ). 36  To obtain  XOOJ  o n e  again uses the fact that dn/x^ = s/f to calculate e and w  then Equations B27&30 to get £ and Equation B32 for cos4w. Since e is odd in x 2  w  about u — J, one can use its sign to determine which value to choose from the arccosine operation that gives u>: € > 0  for  0 < Au < n  € < 0  for  TC  w  for w  To and  calculate u>i = t'71  < Au)  <  u real.  (B33)  lit  and w = | +i*7 , corresponding to x (0 < x± < s) 2  2  x  x (/ < x < 1), one makes use of the fact that: 2  2  s(l-x ) t  >' V(i^;  dnA  =  °^ ^  for  x  ( )  1  B34  so that:  _ y/Ajl-xJ-y/xj  Sf  For j = 1 or 2, cos4u>y = cosh 47y and, for 0 < »7y < ^ : 7y = ^ l n cosh 47y + yjcosh 47y — 1 2  (B36)  2  0*  When evaluating the •^•(w) functions, the procedure to follow is the same as in the previous section up to and including Equation B25. Then: 2K  2  2  k  sn(A ,A:2)cn(A , k ) 2  2  2  2  V^2  where, from WW22.42: K = 2  (1 +  Vk1)7  K .  Also:  dn(A ,A: ) 2  2  2e  0  2  '  V 1  (  4 e  o0  2  Elliptic Functions /  B.3  151  The sign of the sn function is determined by where w is in the rectangle. If: 2  ™S±V,_(&)  (B37)  then:  "  *J  '  IT  1+  dn(A , fc ) 2  When w is real, as when E  I  £3  2  is being evaluated, and when 0 < x 5; §> as  3  is always the case in the rectangle, then sn(A ,A; ) is real and positive and the + sign is 2  2  chosen in T(u;). When u = t'7 or u = | +1'7:  T M = n/(g -i. 2  ±  sn(A , k ) is odd in 7 for x — 0 2  a n  2  d X  §• For the top surface of the airfoil it is positive  =  imaginary and for the top of the flap it is negative imaginary. However, 2e /2e is always w  0  positive on the airfoil and always negative on the flap, sb that in Equation B38: , M . \ 2  ,  2  2e  K  J  Q  2e I n  V V 2e  n  "  J  This allows the evaluation of A{u{)  and  2  ~  ~ 2  in Chapter 4 and Appendix C.  A(OJ ) 2  The elliptic functions of the third kind that are evaluated in Chapter 4 contain the function: »  l n  04(Xoo-Wy)  This is of the form: - In so that:  04(Xoo-^y)  i, n  n- = -(2t6)  l  n  M x ^ ) ^ ^  2 T  , e  {  M  u x  -  +  U  j  )  }  '  Elliptic Functions  152  / B.3  Now, from Equations B2: MXoo+Uj)  = 0 (2(Xoo+Wy),9 ) "  ^(Xoo+Wy),?4)  4  3  ^(2(Xoo+^y),9 ) 4  x  2  =r2%  t> (2(Xoo+"y),<? 4 ) 3  ^( (Xoo+"y),? ) 4  /p^(2(Xoo+^y),g ) 4  dnC/x^+Ay) + \/A7 ' The numerator in this function is:  k  x  - Hi = 1 + 2$ cos4(Xoo+"y),+ 2g cos8(Xoo+Wy) + 2g cosl2(x +wy) + 4  16  36  00  and using the addition theorem for the dn function, WW22.21, the denominator is: R -tI 2  /—. dn /ioo dn Ay — k2 sn /iQQ cn [AQQ sn Ay cn Ay = Vk + — - ^ - — ^  2  a [ « / ( l - « ) ( W ) ] + Vxy(l-xy) " (d=t)|(xy-a)(xy-/)|> T  Thus, if: £ j = 1 + 2g cos 4Xoo cosh 47y + 2g cos 8Xoo cosh 87y + 0 ( ^ . g ) 4  16  Ji = 2g sin 4Xoo sinh 4 4  36  + 2q « sin 8 l  7 y  Xoo  sinh 8 y + O ^ f .g ) 7  36  + Vxy(l-Xy)  R = 2  7 = ± | (xy - s) (zy—/) | 2  5  ,  (B40)  + for airfoil upper surface — for flap upper surface  then:  i t?4(Xoo-"y) m  2  =  a r  (Rj-iJA (B41)  *  11  «? (Xoo+Wy)  *  4  where u>y = i~fi or u/y = | +t'7  2  a  \ R - il 2  u d 0 < :*7y <  2  j  APPENDIX C The Selector Function Expansion Coefficients  This appendix gives the coefficients for the large z expansions of the selector functions L(z)}  Cl £(z)  A(OJ), 7(u), 5(w), and T(w) as well as for the flow function F(z).  £(z)  AND  F(z)  isfirstused in section 4.1:  ^2) = y  ( z  _  a ) ( z  _  i  1 )  +  +  +  ~ +T ^ ^- "  If v = 1/z, then: /t . /0/1 > /i2^.2 = 1 1 +. o2CiV + ( 2 £ + C] )v + ( 2 £ + 2£ii:2)w + _  (1-5V)(1-V)  2  2  3  3  = 1 + [ l - ( / - ) ] v + [ l - ( / - ) ( l + ) ] v + [ l - ( / - s ) ( l + 5+a )]v + s  s  2  5  2  3  and equating coefficients of v : n  2 ^  =  1 - ( / - « )  2£  2  = l-(/-s)(l+s)-£  2£  3  = 1 - (/-s)(l+s+a ) - 2 t i t  2  (Cl)  2  ,_  2  n-l  2£» = 1 - 7 - ^ ( 1 - 0 - £ 4 » 4 . - m  , » >  m=l 153  The Selector Function Expansion Coefficients /  C.2  154  Although F[z) is used rather loosely in Chapter 2, it is given a unique definition in section 3.4: l(z-s)(z-1)  Of course, F[z) = l/C{z) of the £ 's. n  ,  F,  F  F  2  3  ,  as z —* oo.  which leads to a simple recurrence relation for the F 's in terms n  However, if F{z) is being used, C(z) is probably not being used, or vice versa,  so that a direct relationship for the F„'s is desirable:  2F = /-(l+a) 1  2 ^ = [/ -/(l+s) + s ] / 2  n-l  B  _  2  - £ FmK-m  , n>2.  ( C 2 )  m=l  C.2  AND 7{u)  The airfoil selector A{yj) is introduced as Equation 2.38 in section 2.3 and is used again in Chapter 4, beginning in section 4.2. „, , , Jlw =1  2UJ  .  K' \. A  - T T T  A  JT  1  "  . ~  >f  2  z  1/j  Using Equations B7 and since  A  4 +— +  as 2 —• oo.  + z  l  = Xoo ± y  It is convenient to define:  w  K  #,4  .o  (Xoo)  .o  (Xoo)  and using Equations B15, 2.5, and 2.9: £4  = £3 +  (C3)  y  where (?i (and the soon to be used G ) are defined in Equations 2.37. 2  The Selector Function Expansion Coefficients /  C.2  155  Now: dA{u) _ dAju) dX dv dX dv where, again, v = 1/z. From Equation 2.1: 1  dX _ y/fjl-a) 2  dv  (C4)  /(l-sv)(l-/t;)(l-v)  N  so that, using Equations B17: 2K' z(f-s) TT f(z — s)  dAju) dv  2E'  1-fv (1 -av)(l - v) 1-fv  JT  Ai + 2X v + 3A v  2  2  +  3  or: 1-s 2  G^l+sv+aV+sV-r-- • •) - G [l + £ i V + Z v  2  L/  2  2  3  + £v 3  + •••]  = Ai + (2A - fA{)v + (3X - 2fA )v  2  3  2  +  2  Thus: A = 1- E  4  0  l i s -J~ 1  i  G  2A  2  =  G  -  2  (C5)  ( / + s ) Ai + Ai Ci + sG  nA  n  = ( / + « ) ( « - l ) / l n _ i - af{n-2)A n  The flap selector  +  3  +  5  » > 3.  is introduced in section 2.2 and is used throughout  the thesis: ,,  ,  2u  K'$' . .  _  2  Ti  7  2  as z —• oo.  In much the same way as the A„B were obtained, one gets: 2  K'v",  Xoo  _ 2 ~ G  1  2^2 = M  nT = / ( n - 1 ) 7 -i + n  n  ~  2  + ?id  G  + ^ i{C . n  l  G  2  n > 3.  +  1J  + f£ n  + / £ _4 + • • • + 2  3  n  (C6)  Y~ I G  n  2  f~)  The Selector  Function  Expansion  Coefficients  /  156  C.3  The recurrence relations for the A '8 and ^,'s are certainly not unique. The N  relation for the A 's has another version similar to that of the ^,'s, and vice versa. n  The methods of Appendix B ensure the fast and accurate evaluation of the constants G j , G , and E 3 . 2  C.3  S(w) A N D T(w)  The airfoil segment selector is derived in section 4.3:  T  ^(W + H l )  /.  5i  ~5 H 0  z  » *J  5  2  T--5--1  z  ,  ,  as 2 —>  7  = Xoo ^  *  )  00.  7 l  <  2 L 2  ) in the w-plane, Fig. 4d, and  it corresponds to the point X\ which is on the airfoil in the z-plane (0 < 00  7  l  Here, j is on the upper surface of the airfoil (0 < t  Now, So = 5(w ),  C  a n c  *  « 4Uoo+*7l) 7  < s).  periodicity of the theta function  *t?  3  which can be evaluated using the methods of section B.3. Taking the derivative of S(UJ) with respect to z, or v — 1/z as is done here, leads to a doubly periodic function which can be represented by algebraic functions of z: dS(uj) dX  „  „„  .  o  ,\ ' — ~ 5i + 25 v + 3S tT + • • • aA av 2  ,  3  The following expression from dS{w)/dX  ass->oo.  (from the first term on the RHS of Eq. C7) can  be rewritten as follows (WW21, problem 18):  IK  ^•((j-u/J -  ^-(w+wj  2iy/{s-X )(f-Xl) l  TT =  "  t  4  2 T ^  \/l-x N  where:  2  2sn A cn A dn A  t?'  x  +  - y/Ax[  /x7+ /Z(l^x7) N  t  sn'A—n'A, N  /X (1-.T ) 1  Z —  1  Xt  x  ( C Q )  The Selector  Function  Expansion  Coefficients  /  157  C.3  The last line in Equation C9 is in its most desirable form since all of its terms take on only moderate values for any value of x between 0 and s. This, despite the square-root poles x  at Xy — 0 in each of the terms on the RHS of the previous line (they have cancelled out). The other expression from dS(uj)/d\ that needs to be considered is: 2E' dX  K  #  '  2  K' TT  1  2  TC cn A  K  2  TT  2E' K'  2K'dn X  1  f-s2K'z-l  1—3  JT  (Cll)  TT  K  Z—f  where I/J is associated with 7j and is the variable shown in Figure 4c. dX/dv is given in Equation C4, and so: 1  -s/{ -x ){f-x )- =^-— == 3  l  l  T  7  2 if'  G  2  -  (C12)  K  Now, if: (C13) and since, as z —• oo:  1  ~ , *, (i +  1 - sv){l - fv){l - v)  l-fv  +i: v + £ t; H — ) 2  3  2  3  K  1  then: (1 _ f )^M dv v  „ [Si - M v - Miv Q  2  3  - Mv 2  -...](! +Civ  + £v  2  + -•)  2  Si + (25 - fSi)v + (35 - 2/5 )v + 3  2  2  2  and: nS = / ( n - l ) 5 n _ i + n  n > 2.  SiCn-i  ~ M£0  n 2  - M^n-z  M _ £, - M _ n  3  n  2  (C14)  Tie Selector Function Expansion Coefficients / C.S  158  The flap segment selector T(w) is given in section 4.3:  T  _i_  (CIS)  j_ ^ 2 _i_  as 2 —> oo.  Note that this equation is written slightly differently from Equation 4.54 since it is in terms of w = § 4- 172, (0 < t 7 < ^ ) . -  2  2  w corresponds to x 2  o n 2  *  n e n a  P  m  2-plane  (/ < x < 1). Because of the similarity between T(w) and S(w), obtaining the 7J,'s is 2  straightforward. 'o = - in  (C16)  — r r r + ~x"Uoo)  2 If  2  ~  G l  F~~  (C17)  If: = \^x {x -s){x -f){\-x )x 2  then T and N n  n  2  2  2  n 2  - —  G  ^  replace S and M , respectively, in Equation C14. n  n  r  (C18)  APPENDIX D Computer Program: The Tandem Airfoil Forces  This appendix lists the program used to calculate the tandem thin airfoil theory force predictions presented in Chapter 4 in Figures 10 to 17. The program incorporates all of the 13 steps listed on pages 85 and 86. The programming language is Texas Instrument's "enhanced basic", as described in the user manual for the TI CC-40. In order to run the program in the available 5.7 K of RAM memory, stringent memory conservation measures have been used (the number of lines and variables in the program have been minimized). The table preceeding the program, on the next page, lists the main variables and gives the final values they are assigned by the program. The program accepts four quantities as input: 1) the NACA mean line designation for the slat (i.e., "230" for a 23012 airfoil or " 4 4 " for a 4412 airfoil), 2) the mean line designation for the flap, 3) the flap chord to total chord ratio Rf, and 4) the slot size f-s.  Only NACA four and five digit profiles can be processed. The program finishes  by displaying £  m a x  .  The effects of other mean lines (specified with appropriately sized polynomial boundaries) can also be accounted for by the program by manually defining values for the  TT  V  array and  "H","!!",  and  "X4''  variables immediately before line  185.  This allows, for  example, calculation of the effects of slats or flaps attached to the tandem airfoil elements.  159  Computer Program: The Tandem Airfoil Forces /  D  160  MA = dCmJda/(*/2)  AC  —  x  XS  =  X  N  = dC /dri/(2ir)  MN = dCmJdV/(*/2)  SA  = dCJda/(2n)  MSA = dCmJda/(*/2)  SN  dCJdr,/(2v)  MSN = dCmJdr,/(*/2)  e  ac »ac  XF = (*I*c -/)/(! M =  dC Jdr, m  FA  =  dC /da/(2n)  MFA = dCmJda/{*/2)  MS =  dC ,Jd  FN  =  dC /dr,/(2v)  MFN = dCmJdr,/(*/2)  MF =  dC J  DMO =  DM = A C m o c  tf  tf  D = ACt DS DF  AQ.  = AC  DMS =  DMSO = DMFO =  tf  S =s  SI = l-a  F =f  Fl =  l-f  FS =  m  m  V  mj  *c  dri  miac  DMF = A C m  /o  m  Iae  YI == uJK*  Gl = G  f-S  x  CT = 3 + 1- f  G2 =  R = Rf  Y4 == U /K'  G  2  2  KK = 2K/ir  q =Q ER = ^max  X == 2Xoo  100 RAD:DIM A(3,2),B(7.4),L(3.2) 105 PRINT "Enter NACA mean l i n e deBgnatna:":PAUSE:PRINT " S l a t : •; 110 ACCEPT SIZE(3)VALIDATE(DIGIT," ") ST:PRINT • ,Flap:";ST;:Q=ST:J=0 J  116 ACCEPT AT(18)SIZE(-3)VALIDATE(DIGIT.  B  "),FP:PRINT  120 INPUT "Enter R: ";R."Enter f - s : ";FS:PRINT ST;FP;R;FS; 125 I F Q=230 THEN M=.2025:K=1B.957:GOTO 165 130 I F q=210 THEN M=.068:K=361.4:G0T0 165 135 I F q=220 THEN M=.126:K=51.64:GOTO 165 140 I F q=240 THEN M=.29:K=6.643:G0T0 165 145 I F q=250 THEN M=.391:K=3.23:G0T0 165 ELSE I F H<1 THEN H=l 150 M=INT(q/10):P=(q-10*M)/10:M=M/50 155 L(J,0)-M/P:L(J,1)=-L(J,0)/P:J=J+1:L(J,1)=-M/(1-P)-2:L(J,0)=-L(J,1)*P 160 IF,J=3 THEN X4=P:G0T0 185 ELSE X1=P:J=2:IF q=FP THEN 165 ELSE q=FP:G0T0 125 165 L(J.1)=-K*M:L(J.2)=K/2:D=L(J,1)*M/6:L(J,0)=D*(M-3):L(J+1,0)=D*M:H=2 170 I F 3=2 THEN X4=M:G0T0 186 ELSE X1=M:J=2:IF q=FP THEN 165 ELSE q=FP:GOTO 125 175  BREAK:ER=0:X1=X1/S:X4=(X4-F)/Fl  180 INPUT "Enter f - s : ";FS:PRINT ST;FP;R;FS; 185 F1=R*(1-FS):F=1-F1:S1=FS+F1:S=1-S1:CT=1-FS:X1=S*X1:X4=F1*X4+F:F0R 190 N=S~J:A(0,J)=L(0,J)/N:A(1,J)=L(1,J)/N:NEXT:F0R 195  J=0 TO H  1=2 TO 3:A(I,2)=L(I,2)/Fl-2  N=F*A(I.2):D=L(I,l)/Fl-N:A(I,0)=L(I,0)-F*D:A(I,l)=D-N:NEXT:H=H+2  200 K=FS/F/Sl:N=SqR(l-K):D=(l+SqR(N))"2:M=K/(l+N)/D:SA=M"4:P=SqR(l-SA) 205 210  SN=(1+P)*(l+SqR(P))*2:FA=SA/2/SN:FA=(FA~4*2+1)*FA:q=FA".25:FN=q*FA*2 W=FA-4*2:KK=(W+(FN+FA+q)*2+l)-2:DS=-KK*LN(q)/PI:DF=SqR(F*Sl):MSA=KK/DS  216 G2=(K-(8*W+9*FN+4*FA+q)*8/KK*2.5)*DS+2/PI/KK:DMFO=DS*(DS-G2):G2=G2*DF  Computer Program:  The Tandem Airfoil Forces /  D  161  220 G1=FS*DS*DF/S1:MFA=S1/F*G1:B(4.0)=MFA-G2:G0SUB 600:IF C<0 THEN X=PI-X 225 MFN=W*SIN(4*X):DMS0=W*C0S(4*X):DMF=8*MFN/MSA+DS*T+G1/F/AC:B(3,0)=X/PI-DOT 230 B(3.1)=(G2-G1)/2:B(0,0)=1:T=S:B(0.1)=CT/2:D1=FS/S1:DMO=G1*F1/2:FOR 1=2 TO H 235 K,C=0:Z=I-1:F0R J=l TO Z:C=C+B(0,J)*B(0,I-J):K=K+B(0,J-1)*F~(Z-J):NEXT 240 B(3,I)=(DMO*K+B(3,1)*B(0,Z)+Z*F*B(3,Z))/I:B(I+3,0)=MFA*T+B(I+2.0)*F:T=T*S 245 B(0,I)=((T-1)*D1+1-C)/2:KEXT:Z=X1:DMS=(S*F/S1/F1)-.25:MA=SQR(S1/S/N) 250 MSN=DS*3/DF:N=SqR(N)*DF:FN=2*FA*XI:FA=2*FA*SIN(2*X) 255 K=B(3,1)-DF/KK/PI:J=1:G0SUB 605:Y1=Y4:D1=D4:Z=X4:J=2:C0SUB 605:Y4=-Y4:P=2*PI 260 FOR 1=1 TO 3:F0R J=H TO 1 STEP -1:F0R K=l TO J:B(I,J)=B(I,J)+B(O.K)*B(I.J-K) 265 NEXT:NEXT:NEXT:FOR 1=4 TO 7:N=I-4:F0R J=l TO H:T=J+3:C=0:FOR K=4 TO T 270 C=C+B(K,0)*B(N,T-K):NEXT:B(I,J)=C/J:NEXT:NEXT 275 AC.XS,XF,DM,MS,MF=0:SA=S:SN=FS:C=X1-S:T=X4-S:H=H-2:F0R 1=0 TO H 280 N=A(0.I)-A(l.I):D=A(2,I)-A(3,I):M=A(0.I)-A(3,I)-N*Y1-D*Y4:AC=AC+N*C 285 XS=XS+D*T:XF=XF+M*SN:K=I+2:SA=S*SA:SN=(F-K-SA)/K:C=(XrK-SA)/K:T=(X4*K-SA)/K 290 DM=DM+N*C:MS=MS+D*T:MF=MF+M*SN:NEXT:AC=(AC * D1+XS * D4)*4-DMF0*XF*P 295 XS=DMFO*MF*P-(DM*D1+MS*D4)*4:N=B(3.1)/B(0.1):XI=B(3,0):C=CT*(1-XI-G1/F) 300 SA=(C-B(4,1))/S:SN=(C*N-B(7,1)-DMF0*FS)/S:Z=B(0.1)-2:MA=(1-B(0.2)/Z)*2 305 MN=(N-B(3.2)/Z)*2:MSA=B(4,2):MSN=B(7,2):F0R 1=1 TO H+2 310 B(0.I)=B(0.I)-BC3.I)-BC1,I):B(4.I)=BC4.I)-B(7,I)-B(5,I):B(3.I)=BC3.I)-B(2.I) 315 B(7.I)=B(7.I)-B(6,I):NEXT:D,DS.DM0.DMS0=O:F0R J=0 TO H:T=J+1:M=J+2 320 FOR 1=0 TO 3:K=I+4:D=D-A(I,J)*B(I.T):DS=DS+A(I.J)*B(K.T) 326 DMO=DMO+A(I.J)*B(I,N):DMSO=DMSO+A(I.J)*B(K.M):NEXT:NEXT:D=2*P*D/CT 330 DS=(C*D+AC+P*DS)/S:C=CT*C:T=CT*B(4,0)*2:I=S*S:MSA=(C*MA-T+MSA*4)/I 335 MSN=(C*MN-T*N+(F+S)*FS*DMF0*2+MSN*4)/I:DM0=D/2+DM0*PI/Z 340 DMS0=(C*DM0-T*D/4+XS-P*DMS0)/I:I=CT/Fl:J=S/Fl:FA=I-J*SA:FN=I*N-J*SN:P=PI/2 345 DF=I*D-J*DS:K=I*I:T=J*J:Z=F/F1:MFA=K*MA-T*MSA+Z*FA*4:MFN=K*MN-T*MSN+Z*FN*4 350 DMFO=K*DMO-T*DMSO+Z*DF:AC=-MA*CT/4:XS=-MSA/SA/4:XF=-MFA/FA/4:M=(MN-N*MA)*P 355 MS=(MSN+SN*XS*4)*P:MF=(MFN+FN*XF*4)*P:DM=DMO+D*AC/CT:DMS=DMSO+DS*XS 360 DMF=DMF0+DF*XF:ER=(ER*2*W)"2*8:Z=H+2:IF FKS THEN ER=ER/F1~Z ELSE ER=ER/S"Z 365 PRINT:PRINT "Max err=";ER:PAUSE:GOTO 175 500 I=-l:AC=SqR(S/F/N)+l:XS=l-2/AC:G0T0 510 605 I=0:XF=SqR(Z):MS=SqR(l-Z):AC=DMS*MS+XF:XS=1-2*XF/AC 610 MF=(l-(M*XS)-2)/P:DM=l+SqR(MF):C=XS/M 515 XI=SN*(C-2*2-l-XS-4)/(l+MF)/(DM*P)-2:XI=((XI-2-l)*2*W+l)*XI 520 T=ABS(XI):IF I AND T>1 OR NOT I AND T<1 THEN XI=SGN(XI) 625 T=D*C*SqR(ABS(C-2-l)/P)/DM/2*SA:IF I THEN X=ACS(XI)/2:RETURN 630 MF=XF*MS:DM=MF+N:I=SqR(XI"2-l):C=SGN(C):IF ER<XI THEN ER=XI 635 XS=DMS0*(2*XI-2-l)+FN*XI+l:D4=LN(XI+I)/PI/2:Y4=D4*MSA*C:T=16*W*I*XI*C/KK+T 540 I=(MFN*2*XI+FA)*I*C:MN=Z-S:XI=SqR(MN*(Z-F)):MF=MF*XI/PI 645 B(J.0)=ATN((XS*XI-DM*I)/(XS*DM+XI*I))*2/PI+Y4*DMF 650 B(J,1)=(XI*(DMS*XF-MS)/AC-DF*T)/PI-Y4*K 655 D4=D4-MSN*XI/ABS(MN)*(MS+MA*XF)/AC-DS*T*C:AC=Y4*DM0:F0R 1=2 TO H:XI=I-2 660 XS=0:FOR T=0 TO XI:XS=XS+(MF*Z-T+AC*F~T)*B(O.XI-T):NEXT:XI=I-1 665 B(J,I)=(F*XI*B(J,XI)+B(J.1)*B(0,XI)-XS)/I:NEXT:RETURN  

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0096761/manifest

Comment

Related Items