Comments
Description
Transcript
N73 - ,,_sa ,/5 NASA CR-2228
N?3 25045 N73 - ,,_sa ,/5 NASA NASA CONTRACTOR CR-2228 REPORT PART I m C_I " I Z AN IMPROVED ANALYSIS IN OF SUBSONIC Part I- by F. A. Langley NATIONAL THE AERODYNAMIC CONFIGURATIONS AND FLOW and SUPERSONIC Application Woodward by AEROPHYSICS Bellevue, FOR WING-BODY-TAIL Theory Prepared for METHOD RESEARCH Wash. CORPORATION 98009 Research Center AERONAUTICS AND SPACE ADMINISTRATION ° WASHINGTON, D. C. ° MAY 1973 AN IMPROVED OF METHOD FOR THE WING-BODY-TAIL IN PART SUBSONIC I By Analytical AERODYNAMIC ANALYSIS CONFIGURATIONS AND THEORY SUPERSONIC AND FLOW APPLICATION F. A. Woodward Methods, Incorporated SUMMARY A new method sure distribution tail combinations A computer program calculations. has been developed for calculating and aerodynamic characteristics of in subsonic and supersonic potential has been developed to perform the the preswing-bodyflow. numerical The configuration surface is subdivided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. The normal components of velocity induced at specified control points by each singularity distribution are calculated and make up the coefficients of a system of linear equations relating the strengths of the singularities to the magnitude of the normal velocities. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations using an iterative procedure. Once the singularity strengths are known, the pressure coefficients are calculated, and the forces and moments acting on the configuration determined by numerical integration. Several this program data. Good achieved. examples of pressure are presented, and correlation between distributions calculated compared with experimental theory and experiment has by been TABLE OF CONTENTS Page SUMMARY.......................... 1 INTRODUCTION ....................... 2 LIST OF SYMBOLS ...................... 3 AERODYNAMICTHEORY .................... 7 Description Derivation Derivation of Method ................. of the Incompressible Velocity of the Compressible Velocity Aerodynamic Representation The Boundary Calculation Condition of Pressures, COMPUTER PROGRAM Program Program Operating Program Program Components Components . .............. Equations Forces, ............ and Moments ..... ..................... 7 7 . . . 34 43 47 57 59 Description Structure Instructions Input Data Output Data ................. ................... ................. .................. .................. 59 59 59 61 76 VERIFICATION ................. 78 EXPERIMENTAL Isolated Bodies Isolated Wing-Body Wings .................... Combinations ................. 82 87 ........................ 91 CONCLUSIONS .................... APPENDIX I: Integration APPENDIX II: Panel APPENDIX III: Sample REFERENCES Procedures Geometry Case Calculation ................. ........................ 78 ........... Procedure 93 .... 94 i00 125 iii 1. Report NASA No. 4. Title AN 2. CR-2228, and Pt. _THOD CONFIGURATIONS I Accession No. Subtitle IMPROVED PART Government - FOR IN THEORY _ AERODYNAMIC SUBSONIC AND AND ANALYSIS SUPERSONIC OF 9. A. Report Nameand Research Addre_ Date May 197._ FLOW mUnder Corporation Subcontract Analytical 187 9320 Bellevue, Washington Sponsoring Agency 6. Performing Organization 8. Performing Organization 10. Work No. 11. Contract 98009 Name and by: Aeronautics Washington, D. 15. Supplementary 16. Abstract Report A new computer E. Inc. Bellevue, Shoreland Drive Washington or Grant No. NASI-I0408 98004 13. Type of Report 14. Sponsoring and method has of been Administration developed has been a wing surface vortex and to is panels. points by each singularity linear equations relating calculating The perform having a normal components linear are strengths a constant large Report Agency Code number the aerodynamic potential of panels, distribution variation in the velocity calculated of and supersonic flow. A calculations. source of distribution and numerical into A pressure subsonic the subdivided distribution the the in distribution. distribution tall for combinations developed singularity and the Space wing-body-tail configuration panels, and is make up to the of which used on the streamwise induced singularities each direction at the is specified used control coefficients magnitude contains body of of the a system normal velocities. The singularity control points system of for given using Several numerical satisfy number interative are the and boundary angle of procedure. calculated, and condition attack Once the forces are of tangential determined the singularity and moments by flow strengths acting on are the Key Words of pressure experimental Potential by data. distributions Good calculated correlation by between 18. Author(s)) Distribution Lifting Surface this theory program and are presented, experiment has Distribution Statement Unclassified Theory Representation Solution of 19. Classif. Security Linear known, configuration Flow Pressure Vortex (Suggested the this achieved. 17. at solving integration. examples with which Mach an coefficients by compared strengths a equations pressure determined Covered 20546 program The Period Notes aerodynamic the No. Addr_s and C. characteristics on Code 501-06-01-06 Methods, S. Unit Contractor National an No. APPLICATION Organization Aerophysics 12. 5. Catalog Woodward Performing Box Recipient's WING-BODY-TAIL 7. Author(s) F. 3. I Equations (of this report) 20. Unclassified Security Classif. (of this page) 21. No. Unclassified *For sale by the National Technical Information of Pages 128 Service, Soringfield, Virclinia 7_151 been and of INTRODUCTION A unified approach body-tail configurations originally presented in to the aerodynamic in subsonic and references 1 and analysis supersonic 2. This of wingflow was method has been extended by the introduction of several new aerodynamic singularity distributions which substantially improve its capability to represent arbitrary shapes. For example, the new method permits the analysis of non-circular bodies, provides a more accurate representation of rounded wing leading edges, and allows the determination of wing interference effects in the presence of body closure. A computer program has been developed to perform the numerical calculations. The program accepts the standard geometry input format currently in use at the Langley Research Center, and described in reference 3. The graphics capability of the program of reference 3 may be used to obtain a visual display of the configuration input geometry. In addition, the new program has two boundary condition options available for determining the pressure distribution on lifting surfaces. In the first option, the aerodynamic singularities are located on the mean plane of the surface, and approximate planar boundary conditions applied to determine the singularity strengths. In the second option, the aerodynamic singularities are located on the upper and lower surfaces of the lifting component, and exact surface boundary conditions applied. This results in a more accurate pressure distribution, more computer time. Surface boundary applied in the determination of the Part describes compares isolated contains cluding but requires conditions body pressure I of this report outlines the aerodynamic theory, the input requirements of the computer program, and the program output with experimental data for several wings, bodies, and wing-body combinations. Part II a a detailed complete description of program listing the and The by Mr. tance author wishes to acknowledge E. W. Geller to the aerodynamic given by Dr. T. S. Chow in the solution ment of 2 considerably are always distribution. techniques, the computer and by program. Mr. D. N. computer sample program, case. in- the contributions made theory, and the assisformulation of the matrix Bergman in the develop- LIST OF SYMBOLS A consistent a set of units is assumed throughout Aerodynamic influence coefficient, panel inclination angle 8, or wing parameter (_2 - Matrix of aerodynamic cross-sectional area b Wing span, or c Panel chord C Aerodynamic d Distance or body influence influence major axis length, body slope coefficients, coefficient, of ellipse or reference or chord or wing panel length of control point from singularity point from wing origin, diameter Diagonal e Distance of intersection E Off-diagonal F, G, H Velocity distribution I Integral expression k Supersonic K Kernel function Length of block matrix control block factor, source edge or Body panel slope M Mach number, n Direction cosine of panel component normal to panel or force, or tip functions scaling line panel matrix m Normal of edge coefficient D N tangent panel report. _i ) A thickness this or vortex, iteration or body number length dy/dx pitching number moment normal of vector, aerodynamic or velocity singularities 3 NW Number of wing and NB Number of body singularities q Magnitude r Radial R Reynolds s Auxiliary S Wing t of tail velocity singularities at control point distance number variable reference Auxiliary thickness area velocity T Tangential force U, Components of distribution induced function, or wing velocity Vw W V Induced X, Cartesian velocity at control coordinates of point points Y, z Greek Angle of 6 Mach number Y Ratio of singularity attack parameter, specific heats strengths Inclination angle & Incremental value E Minor 0 Inclination Tangent cosine axis (i of of - for M2) ½ air, or panel with x panel with x,y aerodynamic axis ellipse angle of of panel sweepback angle of coordinate transformation plane dx/dy, or direction Sweepback angle Direction cosines Integration of coordinate variables along transformation x and y axes n Ratio of circumference to diameter of Radial distance of control point from line through wing panel tip intersection Velocity potential, and x axis × Integration or angle between variable Velocity component normal to Subscripts B Body base Body base c Wing camber D Drag i Index of panel J Index of influencing k Index of panel L Lift M Pitching max Maximum N Normal P Pressure t Wing moment force thickness control corner point panel point panel circle streamwise velocity vector T Tangential W Wing Xf Refer Y, z 6 to force x, y, z axes AERODYNAMICTHEORY Description of Method The configuration surface is divided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. A typical configuration panel subdivision is shown on Figure i. Analytical expressions are derived for the perturbation velocity field induced by each panel singularity distribution. These expressions are used to calculate the coefficients of a system of linear equations relating the magnitude of the normal velocities at the panel control points to the unknown singularity strengths. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations by an iterative procedure. The pressure coefficients at panel control points are then calculated in terms of the perturbation velocity components, and the forces and moments acting on the configuration obtained by numerical integration. The following perturbation singularities, tion equations, coefficients, non-standard are given in Derivation paragraphs describe the derivation of the induced by the aerodynamic solution of the boundary condiused to calculate the pressure on the configuration. Two repeatedly in these derivations velocity components the formation and and the procedure forces, and moments integrals appearing Appendix I. of the Incompressible Velocity Components Formulas for the perturbation velocity components u, v, and w induced by the aerodynamic singularity distributions in incompressible flow are derived by superposition of elementary line sources or vortices located in the plane of the panel. The resulting expressions are subsequently transformed by Gothert's rule to subsonic obtain the compressible and supersonic flow. Elementary induced by a distance a _ line point from source.- source of the origin v = velocity The component velocity unit strength is given by: 1 at a located formulas point on for P(x, the y, x z) axis (1) + y2 + z 8 0 -4 4_ -4 4_ -4 O_ Figure 1 - Aerodynamic Representation 0 -_1 _-_1 XE_ _0 The source velocity is and the field directed point P. along the line joining the point The u, v, and w components of velocity at the point P induced by a unit strength line source coincident with the axis and having a length _ is obtained by resolving V into x, y, and z components and integrating with respect to E. geometry is illustrated on the following sketch: x its The _y £ _V X sin Z w__ v Y u =0J£V_ cos 1 _ [i dE (x- = _0J_x- _) d_ (2) _;2+ y2_ z2],_2 v=/ v s,ncos 0 0 = (x - _)2 + y2 _ sin 8 dE _)2 d_ + y2 + z2] + z21%_ 4_ r 2 3/2 2 I" W 0J_ V sin £ = z / 41T [(X- _ z 4nr 2 x__ d2 (4) dz 0 9 where r = /y2 x I = dl =_ z2 x x2 = The + + tan-* three r2 x r - x2 = d2 = /(x 6 8 components of velocity x = - £ - tan -I satisfy £)2 + r2 zY Laplace's equation, velocity potential since _u _v _w = ° and of u the line The a _--_x'v = 8__y, w (x The that respect to sin _ i 4_r There i ) used in the as the basis this report. The velocity coincident sin - = [(x _ 6) 2 is - at with of the the point P the x axis Biot-Savart's more induced and law to d_ + y2 + z2 normal _)2 r + y2 to + the plane z z] ½ and containing the integrating with xl 1 dl is no axial component The v and w components and z components. Thus, u= is 6, V vortex. its y source is derived velocity vector and the point P. Noting _ £ is obtained by applying the vortex and integrating. v ¼ 0/ axis _-_z' where Elementary line vortex.unit strength line vortex having a length each element of x = source. elementary line source distributions complex by = 0 (5) of are velocity obtained induced by by resolving the V line into (6) The components equation. v = - V sin e - w 8 = V cos notation is of velocity The elementary of the more complex Care must be taken vortex lines form vortex theorem. z 4zr 2 =-Y 4zr2 defined can [xz d 2 xl]I d (7) [ x2 d2 xl] dl (8) following equation (4). The three also be shown to satisfy Laplace's line vortex solution is used as the basis vortex distributions derived in this report. during these derivations to ensure that all closed rings and thus satisfy Helmholtz's Rotation of coordinates.In the following applications, the line source or vortex coordinate system is in general rotated with respect to the reference coordinate system of the panel. Using primed coordinates to refer to the rotated line source or vortex, and defining I = tan A to be the tangent of the sweep angle of the rotated system, the following coordinate transformations apply: x' = Ix + y (i + 12)9 (9) y, = ly - x (i + 12)½ (i0) z (ii) z ' = in The geometry the following of the sketch: rotated coordinate A system is illustrated x, P(x,y,z) X ii The distance d from the field point to the origin is unchanged in this transformation, but the perpendicular distance of the point from the line source or vortex is given by =f(x r' The coordinate - velocity system = (i - + = w = (i are transformed into the reference v' 12)% Iv I + v +z 2 components as follows: lu' u ly) 2 + (12) u l 12)½ (13) w' (14) Constant source distribution on unswept panel with streamwise taper.The velocity components induced at a point P by a constant source distribution in the plane of the panel are derived by summing the influences of a series of elementary line sources extending across the panel parallel to the leading edge. The geometry of the elementary line source located a distance _ from the leading edge and having a strength dE is illustrated in the following sketch: b 1 (_, _/line sou_c y m _ _, b + m2_) i 2 c 4 P(x, y, z) x lies In in left end 12 the the of following derivation, it x, y plane. The distance the line source is d I = is assumed that of the point P [(y - m1_) 2 + (x the from - panel the _)2 + z2]½ and the distance from the right end of the line source is d 2 = [(y - b - m2_) 2 + (x - _)2 + z2]½ . The panel edge slopes m = dy/dx may be arbitrary. The velocity components are obtained by applying a 90 degree coordinate rotation to the line source velocity formulas given by equations (2) - (4), and integrating across the panel chord as follows: C u = -v' -i _ = / (X _x - _) _)2 d_+ z 2 [__i__ dl _ y - b d2 - m2_] (15) 0 C v= u' i/[i = 1 d2 ] d_ (16) 3 c w = w, = -z _w Only the as the second translation. / 0 (x d_ _)2 - first integral integral may For the same only at the lower limit. respond to the influence ing edge. Denoting these ul = __I 4n _ [ + m;m12)½ (i sinh -I Y -i l 4_(i wl sinh ---- = m12) 1 -, 4-_ tan The corners origin + z2 [_[- dlm_- [- bd2 The resulting velocity components of the inboard corner of the panel results by the subscript one, l [(y mix) x 2 + + mly (I - + corlead- (18) J x -I z(x 2 + y2 -x(y - mix) (17) m12)z2]½ ] ½ m2_] in each formula need be evaluated, be obtained by a simple coordinate reason the integrals are evaluated . sinh_ (x2 + z2)½ V + [(y + + - mlx) z2)½ mlzZ + 2 m,y + (l (19) + m12)z2]½ (20) velocity components induced by the remaining three are obtained by applying the above formulas with the shifted to the corner under consideration, and using the appropriate edge slope. 13 The influence of the influences of the to the corner numbers the complete four corners, shown on the panel where sketch. is obtained by the subscripts summing refer u = uI - u2 - u 3 + u_ (21) v = v I - v 2 - v 3 + v_ (22) w = w I - w 2 - w 3 + w_ (23) The velocity components given by equations (18) (20) are expressed in terms of a coordinate system lying in the plane of the panel. One additional rotation of coordinates about the y axis is required to obtain the formulas used in the computer program. Referring to the following sketch, the panel coordinate system now denoted by primes, is rotated through an angle 6 with respect to the unprimed reference coordinate system. The reference system also has its origin at the inboard corner of the panel leading edge, but the x axis is parallel to the body reference axis. Z ! X ! _x Defining a = tan_, x' = x (i y' = y z' 14 the = coordinate + + az a2)½ transformations are (24) (25) z - ax (i + a2)½ (26) m m' = (I + a2)½ Similarly, the u' u - velocity (27) components become: aw' = (i + a2)½ [ mG 4_(i v = w = v' + = w' (I -G(I + + - H - aF (28) ] aZ)½ + a2)½ 47 (29) au' a2)½ 1 = + 4_(i where F = [F a2)½ + a(mG tan -1(z - ax)(x - z mx) + y2 - z(ay + = sinh(I = x + 2 + (30) z2)½ - mz) 1 H H)] : -x(y G - + sinh- a2 i + m2)½ my + az i [(y - mx) (ay - mz) 2 + (z - ax) 2]½ y (x2 + z2)½ Constant taper.stant derived series parallel source distribution on The velocity components induced source distribution in the plane in a similar manner by summing of elementary line sources to the leading edge. In swept panel with spanwise at a point P by a conof a swept panel are the influences of a extending this case, across the panel the line sources 15 are swept back by the angle A. The geometry of an elementary line source located a distance _ from the leading edge, and having strength dE, is illustrated on the following sketch: y=b , (_, 0) _y --line source c (_ x The of [(x tance d z = panel is assumed the point P from - _)2 + y2 + z2]½ of the line where I is source is the tangent The velocity of the line integrating lu v u P(x, d y, + Ib, b) z) to lie in the x, y plane. The disthe left end of the line source is and the distance from the right end - _x _f-t-_ components are source velocity across the panel - _ leading obtained formulas chord Ib) 2 edge + (y - b) 2 sweepback by rotating through the as follows: + z2]½, angle A. the coordinates angle A, and v' = (1 + 12)½ 1 = 4Tr(l + 12)½ I o + x - _ - ly r2 /el[ [ 1(xL 1 1 d2 4z _) + y dz l 1(x d2 d_ (31) 16 IV V I + U I = (i + 12)½ C 1 4_(i + / 0 12)% I 1 1 d2 l(x- W ----W _- Ib)d2 - + _ r 2 y- - I (x - _) + y ly) b l dl I (32) dE I z(l where grals 1(X d z + 12)½ /Cd__ 4_ 0 r2 r2 = (x - In order to are divided _ - [ l(x- ly) 2 obtain by (i + _)+ dl (i + y- 12 ) the results + 12)½ prior l(x- _- Ib)+ d2 Y - b] (33) z2 in standard to their form the evaluation. inteAs before, only those integrals associated with the inboard edge of the panel require evaluation, and then only at their lower limit. The resulting velocity components correspond to the influence of the inboard corner of the panel leading edge. Denoting these results ul = 4_(I Vl = z _Z by -i + 12)½ [ (i + the subscript sinh-1 x12)½ [(x sinh-1 _ - one, ly) Ix 2 + + [(x sinh- - y(i Ix ly)+ 2 ¥ + i (y2 x + + 12)z (I + ] z2)½ J (34) 2]½ 12)z 2]½ (35) 17 w I i[ - z x2+ y2+ tan -I 47 -xy + l(y 2 + - tan -I z 2) z] -Y (36) The velocity components induced by the remaining three corners are obtained by applying the above formulas with the origin shifted, and using the value of I corresponding to the leading or trailing edge. The influence of the complete panel is obtained by summing the influences of the four corners as indicated by equations (21) (23). Linearly var_in@ source distribution on swept panel with spanwise taper.The velocity components induced by a source distribution having a linear variation in the chordwise direction are derived in the same manner as described in the preceding section for the constant source distribution. In this case, however, the expressions under the integral signs in equations (31) (33) are multiplied by _ prior to integration. The velocity components induced by the inboard corner of the panel leading edge are given below: ul _ -I 4_ + I (i x + - Y _y 12)% sinh-1 x + (y2 z[tan_ I L v I 18 = 1 -4_ _ sinh- - Iz z(x 2 -xy I (x I - I sinh-1 ly) + + y2 [ tan -I t + 1(y [ x + ly) 2 _x + +(i y + 12)z 2]% z2)½ 2 z2)½ + z2)½ + sinh-1 12)½ ] + (37) z 2) 1 (i (y2 [(x - x + [x 2 z(x2 + y2 + z2)½ -xy + l(y 2 + z 2) [(x + - y2 Ix ly) 2 + + + y (I + 12)z 2] ½ z2]½ (38) W, = iI _ (X - + z [(I _ i sinh-1 ly) + It an -l 12)½ sinh x + (y2 z(x2 -xy -I z 2) ½ + + [(x ] y2 + z2)½ 1(y 2 + z 2) - Ix + y ly) 2 + (i - tan- iz] y + 12)z2]½ (39) 1 The velocity components induced by the remaining three corners are obtained by applying the above formulas with the origin shifted, and using the appropriate value of I. The influence of the complete panel is obtained by summing the influences of the four corners as indicated by equations (21) (23). Constant vortex distribution on swept panel with spanwise taper.The velocity components induced at a point P by a constant vortex distribution in the plane of a swept panel are derived by summing the influences of elementary line vortices extending across the panel parallel to the leading edge, and concentrated edge vortices extending back to infinity from the panel side edges. The geometry of an elementary line vortex located a distance _ from the leading edge, and having strength d_, is illustrated on the following sketch: y=b , 0 (_, 0) --y bound _ line vortex P trailing y / X vortices / _ 19 The influence of the bound vortices are considered first. The distance of the point P from the left end of the vortex is dl = [(x - _) + y2 + z2]½, and the distance from the right end of the vortex is d 2 = [(x - _ - Ib) 2 + (y - b) 2 + z_ %, where is the tangent of the leading edge sweepback angle as before. The velocity components are obtained by rotating the coordinates of the line vortex velocity formulas through the angle A, and integrating from the leading edge to infinity as follows: u =;_ lu' 0J - v' ; z 47 K d_ (40) oo 0 v = - lu w = -i 4-_ (41) oo (x - _ - ly)K _) + y d_ (42) 0 where K and = r 2 = i__ [l(x r2 L (x - - l(x - _ d I _ - - Ib) + [ - b] J d2 ly) 2 + (i + 12)z 2 Only those integrals corresponding to the inboard the panel require evaluation, since the outboard edge obtained by a coordinate translation. In this case, both upper and lower limits of the integrals must be edge of can be however, evaluated to obtain the correct results. The resulting velocity components give the influence of a semi-infinite region bounded by the leading edge and the x axis, with origin at the inboard leading edge corner of the panel. They are identified by the subscript one. u I 2O = _1 [ tan -z z(x2 -xy + + y2 _(y2 + + z2)½ z 2) tan-I z ] -Y (43) v I = - lu wl [ 1 4_ - I (44) (i I sinh It should be noted ed by considering edges of the panel integral approaches + 12)% sinh -I (y2 +x -I [(x z2)½ - that the last the influence simultaneously infinity. - I Ix 2 + + y (i ly) log (y2 + + 12)z 1 z2)½ = _0 = _[x AW I _ Therefore, -Y[ 4n x the y2 + [(x - x _) 2 z2 (45) term of equation (45) is obtainof both inboard and outboard as the upper limit of the The edge vortex contributes only to the v and of velocity. The velocity components are obtained grating equations (7) and (8) for a line vortex of length with respect to _, as follows: Avl 2] ½ _ + y2 + z w components by inteinfinite 2]½] J + (x2 + y2 + y2 z 2 + z2)½] (46) + (x2 y2 (47) velocity + + y2 z2 + z2)½] components induced by the inboard lead- ing edge corner of the panel are given by equation (43), the sum of equations (44) and (46), and the sum of equations (45) and (47). The velocity components induced by the remaining three corners are obtained by applying these equations with the origin fluence fluences (23). shifted, and using the appropriate of the complete panel is obtained of the four corners as indicated value of I. The inby summing the inby equations (21) - 21 Linearly varying vortex distribution on swept panel with spanwise taper.A vortex distribution is considered which has a linear variation in the chordwise direction, and lies within the triangular region bounded by the panel leading and trailing edges extended to intersection, and the panel inboard edge. The velocity components induced at a point P by this vortex distribution are derived in three steps. In the first step, the velocities induced by a horseshoe vortex of strength _ d_ having its bound segment located along a radial line from the intersection of the leading and trailing edges are evaluated and integrated across segments sketch. the panel of the chord. horseshoe The geometry vortex are of the bound shown on the y=b , 0 and trailing following _Y (_, O) bound vortex c " d1_ /2 //_ X The origin. bound The __//_ /_ vortices ¸ vortex point is P is located located a a distance distance + z 2] ½ from the inboard end of the2vortex, d2 = [(x - _ - Ib) 2 + (y - b) 2 + z ]½ from In thfs derivation, the slope of the vortex of _, I = 11 + are the slopes respectively. angle ponents 22 A, trailing _ d I and the is from = the [(x - panel _)2 a distance outboard a linear + y2 end. function a_/c, where a = 12 - 17, b = c/a, and 11 and _2 of the leading and trazling edges of the panel, The line vortex formulas are rotated through the as before, of the bound to obtain expressions for the vortex prior to integration. velocity com- The velocity components are given below in integral form: c z u = v =- (48) K_ 7r / lu (49) c w where I = K = 4_ci = 11 + l(x / o (x - _ r2- l(x - I[)K_ These substitution - _) + y _ - = (x - _ + y - b d2 - integrals in terms ly) 2 are of - + - (i ay) 2 ay) + After a lengthy integration ponents induced by the inboard In the following formulas, the leading edge c 4rip 2 z[ 7 12)z - a2z 2 making use variable allz2]. of the X- - x following (51) 2 procedure, edge of the parameter I the velocity companel are obtained. is redefined as the slope. z U + evaluated by the integration 11y)(c (c + lb) m c[(x = L panel (50) a_/c d I r2 d_ --+ lad C (Cl - ( t ax)as - (c - ay)e 2 ] G 1 c i[ p2 (ci - ax)az 2 + (c - ay) s ] (52) F 0 23 V = - (CA W = - (C t _. - ax)(c - ay)t - ay)u/p + (cA - ax)azu/p _) where 2 - azt -iI [ cx+ 4,rrp 2 d + [_ (cA 2 - ay)(cA - ay)s + 7 - ay)e 2 - (C + (cA - a_{,)] J p2 - ax)s _ e__ p_2 [ (c (54) - ax) (c 2c [(c (53) + ae2z2] - ax)az (cX 2] - y + ar--_2 G22c ] G x iic - ax)as (55) F, o and d = [(x r 2 = y2 p2 = e2 = (c - _)2 + (c = tan r2]½ z 2 - ay) 2 + (cA = + - a2z 2 ax) 2 + - ay)(x - G x = I 1 _ sinh-1 e Ay) 2 (56) (I - a_/c)r (cA c[[x- + alz zd -I F p2 Ay + 2 - y(x ax)(x - _) (C - ay)_/C] - _) + y(c 2 + - z2[l ay) + - az 2 (A - a_/c2]]½ (57) 24 G = --1 sinh X -- (58) r It should from those be noted defined The functions axial that the following distribution of can be determined velocity u for z functions equation FI, GI (30). and vorticity corresponding by examining the = 0. From - c(x equation G2 differ to these behaviour of the (51), - ly) U 4 (C Along Along the the panel trailing x vorticity inversely edge, = c + x = 12y, ly and therefore / z y2 + 0. ay) distribution as the local [ _ d_ 4_C = (59) - is seen to vary chord spanwise. linearly The contribution of the trailing vortex originating the inboard edge of the panel is considered next. This contributes only v and w components of velocity, which tained by multiplying equations (7) and (8) for a line of infinite length by _, and integrating. The results follows: _v = u therefore -c = 4(c the and ay) 2 leading edge u Thus, wise, - 1 + Z2 [ [(X - -2 _) + y + z2]½ -z/x t chord- along vortex are obvortex are as ] 1c = 8-Jc" _-_ x - _ + [cx- _)_ + r_]_+ G_o z[ 4zr2 x - c + [(x x - _ - c)2 + r2]½ ] (60) Similarly, Aw= -y 8_C + I x- r2 4nrY 2 + [(x - 612 + rZ] ½ I 1c + Gz 0 Ix - c + [(x - c) 2 + r 23½] (61) 25 The first term in the braces gives the velocities induced by a pair of line vortices of quadratic strength along the x axis, and the last term gives the velocities induced by a linearly varying vortex from the panel trailing edge. The combination gives the contribution of a line vortex of quadratic strength to the trailing edge, followed by a constant vortex of strength c/2 extending downstream in the wake. A constant vortex of equal but opposite strength trails downstream from the outboard tip of the triangular panel. In the second step, the velocities induced by a vortex distribution having a linear variation in both chordwise and spanwise directions is derived and subtracted from those given above to obtain the velocity components corresponding to a vortex distribution having a linear variation chordwise, but remaining constant spanwise. In this step, the bound vortex located along the radial line from the intersection of the panel leading and trailing edges is given a linear variation in the spanwise direction prior to performing the chordwise integration. The linearly varying bound vortex is made up by superimposing a series of horseshoe vortices of strength _d_dq with inboard edge located at q, and outboard edge located at b. The geometry is illustrated below: y=q 0 y=b _y + c , ( P X 26 I + lb, b) The contribution of the bound segment of this elementary horseshoe vortex is obtained from the line vortex formulas, with the origin shifted to the point (_ + IH,_), and the coordinates rotated through the angle d I = [(x - _ - I_) _ + the bound vortex, and from the outboard end. in integral form: (y d2 A. The c u = v = w - 4_c - P is r2 0 -i distance dn (62) 0 (63) c [ _ (x- J 0 11 + a_/c a = I - 1 K = 2 1(x _ - ly) d_ - 1 _ - In) + y 1(x - - _ d1 (x (64) dD r2 0 = = a inboard end of - b)2 + z2]½ are given below lu I r2 located from the Ib) 2 + (y components b K 4zc where point _ _)2 + z2]½ = [ (x - _ The velocity - _ - - Ib) + y - b d2 ly)2 + (i + 12)z 2 Only the first term in the K integral requires the second term cancels in the superposition grating this with respect to n, evaluation, as process. Inte- b I = 1(x d : [(x d and d 2 _) + 0 = where - is - y dl - (i dn + 12)_ (65) d 2 the _)2 same + y2 as + z2]½ previously defined. 27 The integrals (60) - z u (d = = - d 2)_ d_ (66) c 4_c v (62) may now be written o r2 -lU (67) c W _- These by equation board edge the panel -i 4_c / o (d - d2)(x r 2 _ - I[)_ d_ integrals are evaluated, using (51). The velocity components of the panel are given below, leading edge slope. -c z[x- U 4_p 2 z[ p2 (cl 2(CI - ax)cs - - ax)(y(c e2(y(c p2 - (68) the substitution given induced by the inwhere I is redefined as - ay) ay) - - az2) az 2) ] G2 - zd G z c i[ + V (cl - ax)cz + s(y(c- ay)- az2)] (69) FII 0 where 28 v = - (cl w = - (c c t- - - ax) (c ay)t 4_p 2 c + -i [(Cl + - ay)u/p (cl - 2 - azt ax)azu/p 2 ax)(y(c - (71) az' p2 - (70) ay) - az2)s + ce2z 2 2] + B 2c G2 + + e2[ V (y(c -p2 (cl - - ay) - ax)cs az2)s + c(cl e2 - (y(c - ] - ay) ax) z 2 - az 2) GI Iic (72) F i 0 and the ing equation remaining The functions distribution velocity equation -cy The is (x - axial the u are defined follow- vorticity given by corresponding the value of to u for z these new = From 0. ly) If ay) 2 velocity trailing is edge, = the c + leading edge, and 12y, the ay) new axial velocity function ay = 4(c - the = ay) 4(c combined - = sF 2 I by a/c + z [ e2G -- 4 functions 1 - ay) on the panel, which is along the trailing edge, direction. The velocity distribution are given 41Tp multiplied 1 + Thus, is from the original, the value of u along the trailbe constant. This can be seen by multiplying by a/c and subtracting from equation (59). The -c u along x (73) - and subtracted ing edge will equation (73) result is: U zero where -cy = 4(c bution stant wise vortex variables _--- 4(c along of functions (68) U and (55). give the zero along and varies components below: I - (cA - desired the leading linearly corresponding ax) G 2 II c vortex distri- edge, in the to conchordthis (74) 0 v = - (cl w = - (c - - ax)(c ay)t - + ay)u/p (cl - 2- (75) azt ax)azu/p 2 (76) 29 where : cls[ t 4_P 2 V ( Y + e2Gl arc 2 )G 2 - (Cl + cA - F ax)G2] c ax d 1 (77) Ic 0 and the remaining ing equation ity functions simpler than The functions (55). given either and variables are defined It by of should be noted that the final velocequations (74) (77) are considerably the preceding sets. of the derivation velocity component formulas vortex distribution is completed by adding the the wake. Returning to the sketch on page 26, that the elementary horseshoe vortices generate vortex sheet of constant strength. This vortex utes nent (7) only to the v and w of velocity will be for a line vortex of components derived infinite c AV _ d, and I The = [(x = _ inner AV - this of velocity. The first by integrating length, as follows: v compoequation (Y - _ _)2 (Y - [1 + z2 + z2] ½ + x - _ 1 - )'q,] (78) 0 - + f d_ 0 where for contribution of it can be seen a trailing sheet contrib- b f --a---{z 4nc2 = follow- _ - ln) 2 + n) 2 a_/c integral is tanc[ [ -a 4_c 2 evaluated first, giving [(x -- _)+ 2 + r2]% -y(x _r 2 i _tan-_ z] _ d_ 0 -tan-1 = -a 8-_ I [tan-Z zt z -y(x (x- (x c where r 3O = u (y2 and + t z2)½, - ly)(c (c are given and _ - by is - _)2 _) - + + ay) ay) 2 + + alz 2 a2z 2 equations redefined z ] J rIr2 2 ½ (74) as the ul c (79) 0 and (77) leading respectively, edge slope _,. The w component of velocity by integrating equation (8) for Here, c b 4zc2 where d I is The _ d_ 0 0 defined inner and is I = 11 evaluated + manner, length. 1 + x - _d, - In ] (y _ D)2 + z 2 above, integral is derived in a similar a line vortex of infinite (80) a_/c. first, giving c Aw = sinh 4_c 2 - (i o sinh + the last Aw= [(x ½ - _ y + l(x - ly) 2 + _) (i + 12 )z 2] ½ x- -I + r Only Thus 12) -I two log r] integrals a 4Trc 2 [I1-I _ can +I] 2 (81) d_ be evaluated in closed form. (82) 3 where II c : (i + I_ 12)½ sinh- i [(x- _- Y + l(xlyi _ + _) (i + (83) d_ 12)z21½ 0 I2 'I = 4 (3x + _) [ d - (x - _)G2] + d2G2 1c (84) o c2 I - log a (85) r 2 where I equation computer It = 11 + (55). a_/c, and Equation d, r, and G 2 are (83) is integrated defined following numerically in the program. should been multiplied rectly account be by for noted -a/c the that Av prior to contribution and Aw as integration of the derived in wake. above order have to cor- 31 The velocity components induced by a vortex distribution which has a linear variation in the chordwise direction, and remains constant in the spanwise direction have now been derived for a triangular region bounded by the panel leading and trailing edges, and the inboard side edge. In the third step of this analysis, these velocity component formulas are combined to give the influence of a swept, tapered panel of arbitrary span. This is accomplished by superimposing two of these triangular regions having common outboard intersections and equal values of the leading and trailing edge slopes. The superposition process is illustrated by the following sketch: y=b y=c/a y a C 2 C 32 = 12 - 11 The upper triangular panel has a concentrated vortex of strength c I trailing from the inboard edge, and a vortex sheet of strength a behind the trailing edge. There is no concentrated vortex shed from the outboard tip, since the circulation around the trailing vortex sheet is equal and opposite to that of the concentrated edge vortex. A similar vortex pattern is shed by the second triangular panel, except that the concentrated vortex has a strength c 2. The influence of a swept, can be obtained by superimposing indicated. It should be noted trailing if the vortex from panel sheet tapered panel of finite span b the two triangular panels as that the concentrated vortices the edges of this panel is tapered, the difference in the wake. The vortex are of unequal being made distribution up on strength by the the panel is zero along the leading edge, and varies linearly in the chordwise direction to a constant value along the trailing edge. The axial component of velocity u is given by equation (74), the v component of velocity is given by the sum of equations (60), (75) and (81), and the w component of velocity is given by the sum of equations (61), (76), and (82). care If the influence must be taken in of a triangular the evaluation panel is required, of equations (74) special and (77). In this case, the chord of the outboard panel subtracted in the superposition process is zero, and two terms in the equations become indeterminate. The limiting values of these terms are given below. First, the function [GI] c = 0 _ 2 are the GI becomes: 2 lim c÷0 where and 1 2 log (x (x - slopes fly) 2 X2y)2 of + + the (i (i panel + + 2 l_ )z X2)2z2 leading (86) and trailing edges. Second, the last two terms in the expression for t become: C lim c+0 The [i c remaining (r2Ga + xd)] = (x 2 + r2)½ (87) 0 terms in the equations are unchanged. 33 Derivation of the Compressible Velocity Components The compressible velocity components induced by the source and vortex distributions are obtained by applying Gothert's rule to the incompressible velocity components derived in the previous section. The original derivation of Gothert's rule presented in reference 4 considered only compressible subsonic flows; here the rule is extended to include supersonic flows as well. The extended rule states that the velocity components u, v, and w at a point P(x, y, z) in a compressible flow are equal to the real parts of u i, 8vi and 8wi, where ui, vi and wi are the incompressible velocity components evaluated at a point P(x, BY, 8z), and 8 = (i - M2)½. In subsonic flow, this rule agrees exactly with that given by Gothert if each of the compressible velocity components are divided by the constant B2. In supersonic flow, the compressible velocity components become complex functions, and care must be taken to extract the real parts of these functions in order to obtain the correct results. However, this procedure evaluating the velocity a straightforward method fields corresponding to tion. A forming simple example the velocity along the equations that d I = terms are is generally much simpler than formally components by integration, and provides for obtaining the supersonic velocity any existing incompressible flow solu- of the components extended induced x axis. The velocity (2) (4) are unchanged (x 2 + B2r2)½ and d2 = real in subsonic flow; rule by is obtained by transa line source located component formulas given by this transformation, ((x - _)2 + B2r2)½. Both but in supersonic flow, by except these both are imaginary ahead of the Mach cone from the origin, d I is real but d 2 is imaginary between the Mach cone from the origin and the rear Mach cone, and both are real behind the rear Mach cone. Thus the velocity components are zero ahead of the Mach cone from the origin, and the finite length of the source has no influence on the velocity field except within the rear Mach cone. Considerable advantage is taken of this ability to correctly define ity component basic sented the regions formulas in of influence the following of each term applications. The compressible velocity components singularity distributions used in in the following sections. Constant source distribution on for each this method unswept panel in veloc- of the five are pre- with wise taper.The incompressible velocity components for source distribution are given by equations (28) (30). corresponding compressible velocity components are: 34 the streamthis The u -- -k + B2a2-_ 4_(i /aF_ (B2mG _ (88) v = _a 2)% 4_ ----- G (89) W= k 4_(i + S2aZ)_ IF + a(B2mG _ H)] (90) where F = tan-l -x(y _ (z rex) a×) d o2_. S z (ay - mz) (91) G - 1 e sinh-1 x___' Br' (92) H = _ Sinh-1 By (×2 + B2 z2)% (93) and B2 = 1 - M2 k = f d = e -- [i d' = r' i for 2 for M _ 1 M > I (x2 + S2r2)½ + 82(a 2 + m2)]% de = [(y _ rex)2 + (z - a×)2 + x' = sonic Velocity The COnstant x + 82 (my k Components. gives the + B2 (ay - mz)2]½ az) COrrect Scaling factor for the SUper- 35 In supersonic flow, the real parts of the functions F, G, and H must be determined. The function F is zero everywhere ahead of the Mach cone from the origin, except for panels having supersonic side edges, when F = ± _ within the "twodimensional" region bounded by the Mach waves from the side edge and the Mach cone from the origin. The boundaries of the two-dimensional region are given in the following sketch, which shows the traces of the Mach waves and Mach cone from the origin on a plane perpendicular to the x axis. _ zT 1 = - e 2 ae)x+ my X az t z = + a2 Imz + - ayle' m 2 ax _ y \ \ \ (m + 1 the tion The function G takes relative sweepback on in logarithmic form, x' G or G = = _ e + log several the side different edges. ae')x - e2 forms depending Expressing the d' for e2 > 0 for e2 = 0 for e2 < 0 and - 8'r' --d X I 1 or 36 G = --COS e' --I X' B'r' 8'r' < x' < B'r' on func- or G = + --_ e' for e 2 and x < > 0, my x' + _< -B'r' az + a2 or G where B' = e' = Finally, = in H = 0 lay + [-i - 1)% - B 2 (a 2 supersonic 8' + flow, tan -Id mzle' m 2 elsewhere (M 2 - (94) m 2)] ½ the for function x > H becomes: 8r B'y or H = 0 elsewhere Constant source distribution _.The incompressible velocity distribution are given by equations ing compressible velocity components (95) on swept panel with spanwise components for this source (34) (36). The correspondare: -kG I u - v = -4_ FI IG I k[ w- where (96) 4n 4_ = F 1 - - tan_1 GI F 2 z -xy F 2 (97) = tan_ I = 1 -- sinh-1 e + ] (98) d (99) Ir 2 z -Y (100) Ix + 82y ly) 2 + (82 (i01) 8[(x- + 12)Z2]½ 37 G2 = sinh_ I 8rx and (102) 1 for M _< 1 2 for M > k d = (x 2 + 82r2) e 2 = (32 + 1. 2 r 2 = y2 + Z2 1 ½ In supersonic flow, the real parts of the functions Fz, F 2 , G z and G 2 must be determined. The function F 2 is always real, and can be dropped from equation (98) without affecting the results since the contributions from the four corners of the panel ahead of always the Mach cancel. cone The from the function origin, F z is except zero everywhere for panels having supersonic leading edges, when F I = ± w within the "twodimensional" region bounded by the Mach waves from the leading edge and the Mach cone from the origin. The boundaries of the two-dimensional region for this case are shown on the following sketch: f f y = ix/8 ,2 / x = / x Xy / / \ J 38 / / --, _-y + e'lzl The function on the relative the function in 1 e GI or GI or G I = or G = = log + where 8' = (M 2 e' = (-82 d' = or function different panel leading for e2 > 0 for e2 = 0 for e2 < 0 and -8'r' for e2 and x forms depending edge. Expressing x w 8'r' -e' G The x' + d' 8'r' 1 --cos e' or 1 several of the form, d _ i = G I takes sweepback logarithmic 0 < < > 0, ly x' < 8'r' x' < -8'r' + e'Iz I elsewhere - i)½ - 12) ½ (103) x' = Ix r' = [ (x + 82y - ly) 2 + e2z2]½ ed G2 G2 = log G2 = 0 becomes: x + d -_ 8'r for x > 8r elsewhere (104) Linearly varyin_ source distribution on swept panel spanwise taper.The incompressible velocity components source distribution are given by equations (37) (39). corresponding compressible velocity components are: u- v -k 4_ - [ (x- k[ 4_ (x ly)G - I ly)(IG, + yG 2 - - G 2) z(F I + x -F + d with for this The 2 )] - Iz(F (105) 1 - F 2 ] ) (106) 39 W where = the _ in the - functions by equations of functions and the real (104). The by d from (x ly) "te o- (F 1 FI, F2, GI, G2, d, e, k, (107) and r are defined (99) (102). In supersonic flow, the behaviour F I and F 2 is described following equation (102) parts of G I and G 2 are given by equations (103)and sum (x + d) appearing in equation (106) is replaced supersonic origin. Constant flow, vortex and distribution _.The incompressible vortex distribution are corresponding compressible U = k 4-"_ v = -lu w =- IF is I _ real on only swept within the panel Mach with cone spanwise velocity components for the given by equations (43) - (45). velocity components are: bound The Fa ] (108) (109) GI 4_ - IG2 ] (ii0) where the functions F I, F 2, G,, G z, d, e, k, and r are defined by equations (99) (102). In supersonic flow, the behaviour of the functions F I and F 2 is described following equation (102), and the real parts of G I and G 2 are given by equations (103) and (104). The is given contribution of Au = 0 edge vortices in compressible flow Av = kz 4_ Aw = - (iii) x ky 4_ In supersonic flow, equations is replaced cone from the origin. 4O the by + d r2 (112) x (113) + d r 2 the sum by d, (x + d) appearing and is real only in the above within the Mach Linearly varyin@ vortex distribution on swept panel with spanwise taper.The incompressible velocity components for the bound vortex distribution are given by equations (74) (77). The corresponding compressible velocity components are: -k p2 I sF + 1 u - V = - (Cl W = - (C t - ck _ s 4_P z ( V z [e2G - 1 (Cl - ax)G2] 0 - - ax)(c ay)t - + ay)u/p (cl - 2 - ax)azu/p azt (115) 2 (ii6) i where (114) Ic 2 ze [e2G - i (cl - ax)G 2] p F 2 1 C ar2)G 2 C where k + 1 for M ( 1 2 for M > 1 - 1 + 82r2]½ + a2z (cl - ax) d (117) ____ a = I d = [(x- r 2 = p2 = (c e2 = (cl s = (c y2 2 _)2 + Z 2 - ay) 2 - - ax) 2 ay)(x 2 + 82p 2 - ly) + alz 2 41 and zd -I F = -y(x- _) + (I- (cl GI (118) tan 1 _ 1 - a_/c)r ax)(x 2 - _) + 82[y(c - ay) - az 2] sinh-1 8c[[x e - ly + (C - ay)_/c] 2 + z2182 + (I + a_/c) 2]]½ (119) G = --I sinh X -- (120) 8r In described and in G 2 supersonic following are flow, the equation given x' = (cl r' = c[[x d' = e e' : [- - by equations ax)(x - ly behaviour of function (102), and the real - + _) (c (103) + - 82[y(c ay)_/c] 2 and (104), - ay) - + z2182 F I parts is of G I where az 2] + (I + a_/c) 2]] ½ d (cl Finally, compressible Au the = ax) 2 _ 82p2]½ contribution flow is given of by the vortex sheet in the wake (121) 0 I Av -a _ = I k(F I - F 2) zt c C (x - ly) (c - ay) + alz 2 ) (122) + where u and Aw where Ii, the Gothert 42 t = 12 (c are 2k-_a 4_c - given [II and 13 are transformation - ay) 2 + a2z 2 by equations I2 + u (114) and (117) respectively. I3 ] given by equations applied. (123) (83) - (85), with Aerodynamic Representation The source and vortex distributions derived in the preceding sections provide the basis for the aerodynamic representation of the configuration. The strengths of these singularities are determined by satisfying the boundary condition of tangential flow at panel control points for given Mach number and angle of attack. In general, the control points are located at the panel centroids, except where noted below. The body is represented by constant source distributions on surface panels, but two optional methods are available to represent the wing and tail surfaces.(Here, tail surface implies any horizontal or vertical tail or canard surface.) Planar boundary condition option.In this option, the panels are located in the mean plane of the wing or tail sur, faces. Linearly varying source distributions are used to simulate the airfoil thickness, and linearly varying vortex distributions are used to simulate the effects of camber, twist, and incidence. The slope of the airfoil thickness distribution is approximated by linear segments between the panel leading and trailing edges. This linear distribution is constructed by superimposing a series of triangular source distributions extending over two adjacent panels. The strength of the triangular source distribution is determined by the slope of the thickness distribution at the intermediate panel edge, as illustrated below: X / dz t dx x Chordwise slope thickness distributions and x x i-i X i Triangular distribution i+l source 43 The same method is used to approximate the chordwise vortex distribution on the wing. In this case, the strengths of the vortex distributions are determined by satisfying the boundary condition that the resultant normal velocity is zero at panel control points. A typical chordwise vortex distribution is shown below: x Chordwise vortex distribution In subsonic flow, or if the trailing edge is swept behind the Mach line in supersonic flow, the Kutta condition implies that the vorticity goes to zero along the trailing edge. In this case, the control points are located at the panel centroids. If the trailing edge lies ahead of the Mach line in supersonic flow, an additional vortex singularity is added at the trailing edge, as indicated by the dashed line in the above sketch. In this case, an additional control point is added on the trailing edge of the wing, and the intermediate control points adjusted and the trailing on or is swept point is located ing edge control linearly between the leading edge control point edge. If the leading edge of the wing lies behind the Mach line, the leading edge control at the centroid as before, otherwise the leadpoint is located on the wing leading edge. In the non-planar boundary condition option, the panels are located on the upper and lower surfaces of the wing and tail, and linear vortex distributions alone are used to simulate both lift and thickness effects. The upper and lower surface vortex distributions are similar to those described above, 44 and the two vortex sheets are joined together at the leading edge by equating the vortex strengths of the leading edge panels. The resulting continuous distribution of vorticity around the chord is illustrated below. T T In subsonic flow, the non-planar boundary condition option presents the problem that one more control point exists than the number of vortex distributions if the Kutta condition is enforced at the trailing edge. An additional source or vortex distribution must be included to make the resulting system of equations determinate. One way to resolve this problem is to introduce an additional pair of trailing edge vortices having equal and opposite strength, as indicated by the dashed line on the above sketch. Another method is to add an internal line source at some point in the interior of the airfoil. In either method, the strength of the additional line source or vortex approaches zero and has small effect on the final solution. The second method is recommended, however, since the first tends to generate an ill-conditioned system of equations for airfoils with small trailing edge angles. In supersonic flow, a similar problem exists if the trailing edge lies on or is swept behind the Mach line. If the trailing edge lies ahead of the Mach line, additional trailing edge vortex singularities must be added on the upper and lower 45 surfaces, and the strengths of these determined by satisfying the boundary conditions at additional control points on the trailing edge. The remaining vortex strengths are determined as described above. If the leading edge lies ahead of the Mach line, the vortex distributions on the upper and lower surfaces of the airfoil are determined independently using control points located on the panel leading and trailing edges. The influence of the trailing vortices in the wake is included in the velocity component formulas derived for the constant and linearly varying vortex distributions. Since the wake is assumed to lie in the plane of the panel in these derivations, the wake vortices must be rotated at the leading edge of each downstream panel to follow the contour of the upper or lower surface of the wing to the trailing edge. The paths of the trailing edge vortices are illustrated on the following sketch. 46 The Boundary Condition Equations A system of linear equations is established which relates the magnitude of the velocity normal to the surface at each panel control point to the aerodynamic singularity strengths. The geometrical relationship between each influencing panel and control point is required to evaluate the coefficients of this system of equations. Wing and body panel geometry.A typical panel subdivision of a configuration which includes a wing, body, and tail is illustrated on figure 1 (page 8). A reference coordinate system is established with origin at or near the nose of the body, having its x axis on the center line and parallel to the body axis, and a vertical z axis. Since symmetry about the xz plane is assumed throughout this analysis, panels are located only on the positive y (right hand) side of the configuration. The body panel corners are defined by the intersections a series of planes normal to the x axis, and longitudinal ian lines. A maximum of 30 rings of panels may be used, containing up to 20 rows of panels around the circumference. of merideach The top to body panels are numbered of each ring, starting in with The wing and tail surface intersections of a series of x axis, and lines of constant columns of panels all other horizontal may sequence from the the forward ring. panel vertical percent be used, or vertical corners planes chord. including tail or each column may contain up to 30 rows of tail panels are numbered in sequence from the trailing edge of each column, starting column of the wing. bottom the are defined by the parallel to the A maximum of 20 those canard on the wing surfaces, and and panels. The wing the leading edge with the inboard and to For each panel, the corner point coordinates, centroid coordinates, inclination angles, area, and chord length through the centroid are calculated, using the method outlined in Appendix II. It should be noted that the panel inclination angles $ and 8 are related to the direction cosines of the normal as follows: nx = - sin ny = - cos nz = cos 6 6 cos sin 8 (124) 8 47 A primed system of coordinates is introduced, at corner point k of panel j, and inclined at the originating angle 8.3 with respect to the xy plane. For body panels, the x' axis is parallel to the reference x axis, and the y' axis lies in the plane of the panel through the leading edge. The panel is inclined at the angle 6j to the x'y' plane. For wing panels, the x' axis lies in the plane of the panel along the inboard side edge, and is perpendicular to the y axis. The z' axis is normal to the panel. In this case, the x' axis is inclined at the angle _j to the x axis. The geometry of the wing and body panels, trated Z and panel corner point numbering convention, is illus- below: ! Z Z v j' L v X ! x y! y! 2 4 4 1 X v 2 Body panel Wing panel The control point of a panel is defined as that point on the panel where the boundary conditions are satisfied, and each panel has a unique control point associated with it. The control point of panel i is normally located at the panel centroid. Exceptions to this rule exist for wing or tail surfaces using 48 planar boundary conditions, as described in the previous section. The coordinates of the control point are given in terms of the primed system originating at corner k of panel j as follows: For body panels, condition option, and wing panels x' = AX y' = Ay cos e z_ = Az cos O. 3 Ax = xi - xk Ay = Yi - Yk Az = zi - zk using the planar boundary l i where For wing panels using the J + Az sin e - Ay sin O. ] non-planar IIAx + 12AY + 13Az ' Yi = _ + _ + _ z_i = _IAx + _2AY + _3 Az Az are Ax 2 Ay (125) boundary x i' = I J 3 condition option, Az (126) where and Ax, Ay and defined above, -n nz l = I _1 k (nx2 = + nz2)½ nx = 0 3 (nx2 2 _2 = 3 ny (n x 2 ) ½ + nz2)½ n z z 2 + = + -nyn -nxny _I x nz2 _2 : (nx nz2)½ 3 (nx 2 + nz 2 ) ½ (127) 49 The direction with noted n n x, and Y n z are given by equation (124) subscript j applied to the angles 8 and 6. It should that equations (126) reduce to (125) for 6. = 0. 3 The opposite (126) cosines coordinates side of the with by calculate = - Yi panel of xz the image plane are - Yk" symmetry The of control point given by equations image control point be i on the (125) or is used to effects. Calculation of the normal velocity at the control points.The resultant velocity normal to panel i at the control point is the sum of the normal component of the free stream velocity vector and the normal velocities induced by the panel singularity distributions. In the following analysis, the free stream velocity vector is assumed to have unit magnitude, and lie in the xy plane at an angle _ to the x axis. The component of the velocity vector normal to panel i is = sin _ cos i axes 8 cos _ i The three at control - cos _ sin i 6 (128) i components of velocity point i are given by parallel to the following the reference equations: N Aui = _ [ (11 u_. _3 + _ i v'. l] + _ iw_.) _3 +(Ix u_lj + + _lwij -' j=l (129) -' _xvij 'JIIYj N :X +.v j÷ j=l (130) U2 ij j j N Aw i = I:[< 13uij + _3vij + u3wij) j=l (131) 5O where 1], is the total is the strength number of of the are the three components at control point i by v_. 13, primed coordinate singularities .th 3 singularity of panel system velocity j, given induced in the of panel j of j of velocity given in panel j induced the I, are W_. l] U w , ° 13, are the three components at image point i by panel primed coordinate system V w 13, _'. 13 and the coefficients by equations with panel of (127), the where transformation the direction v, are given associated j. The normal component of velocity the panel singularity distributions above velocity components as follows: A_. = i _ l The coefficients (127), except panel U, cosines Au. l + _ 2 _i, that Av. l + _ is at panel i induced given in terms of (132) Aw 3 v2, and _s the direction by the i are also cosines given by equation are associated with i. Combining equations velocity at control point n. = 1 _. + (128) i is and (132), the resultant normal A_. 1 1 N = _. l + _ (133) a..y. l] ] j=l where from the aerodynamic equations (129) For body panels, condition option, the axis, and 6. = 0. In ] - influence (132). coefficient a can be obtained ij and wing panels using the planar boundary x' axis is parallel to the reference x this case, the normal and tangential 51 velocity at control point i can be written N w: 1 lj cos (8.3 -8.) 1 + v'. 13 sin (8 j -8 i) w'jl cos (8.3 + + V_. 13 sin (8 j + 8i)]y.j 3 [v_ lj cos (8.-8 3 13 sin (8.3 - 8i ) cos (8 w_ lj sin (8 + O j=l + 8.) 1 (134) N v"l = _ i )-w_. j=l -V_. + 13 and axes the at three control 8 j components point of are: i ) + i velocity j parallel to )ly i J j the (135) reference N Au i = _ (u'..13+ uqj)Yj j=l AV i = v[ cos 8i - w_ sin 8i _w. = w[' i cos 8 + v[' 1 sin 8. 1 1 Then, from _. = 1 i equation w? 1 cos _. (132) - du. 1 Formation of condition equations (136) sin _. 1 (137) 1 the boundary are obtained condition equations.by setting n. = 0 in The boundary equation 1 (133). The complete set N of equations N = i=lj=l in matrix [Aij] 52 be written: N _aijYj Alternatively, may - _ _i (138) i=l notation, {Yj} =- I_i} (139) where A.. is the matrix of aerodynamic influence 13 and _.1 is given by equation (128). In general, coefficients, the matrix is subdivided into four partitions in order to simplify the calculation procedures. The first partition, ABB, gives the influence of the body panels on the body control points, the second, ABW, gives the influence of the wing panels on the body control points, the third, AWB, gives the influence of the body panels on the wing control points, and the fourth, AWW, gives the influence of the wing panels on the wing control points. Equation (139) is rewritten in terms of these four partitions below. [w If 1 The subscripts present program, W and the B refer maximum The right side of modified if the planar In this case, the slope tan where _ (14o) to wing or body panels. order of each partition In the is 600. the boundary condition equations is boundary condition option is selected. of the wing surface is given by: = ( I_ dzc dx i + dzt dx (141) /i dz c is dx dz__it is the slope of the wing camber surface. the slope of the wing thickness distribution. dx The positive sign to the In addition, sign lower the applies to the surface. normal velocity N A_ i = _a bij NW = cos is at surface, control the point i negative is given by NW i3 .y j + cos e _ j =i where upper bi j (dzt dx (142) 1 /j j =i _i ( wijc°s the number vij sin@ Oi - of wing i - u ij tan _i ) panels, 53 and are the point panel W0 velocity i j by components the source induced at distribution on control wing 0 13 The second the projection of the wing. term in equation of the free stream Combining n i equations = cos + _ L[sin _ ij y j +cos i _aN (128) cos (142) is velocity and e Setting n. 1 = 0, cos _b NW j =i in by cos _, the plane (142) - i multiplied vector _ tan _ i ] (143) i j (dzt) d-_- j j =i the new N boundary condition N equations are: N (144) _aijYj = i=lj=l where _ i = cos - _ cos i _ mi i=l [cos _ tan _ - i sin _ cos 8 i ] NW.(dzt 1 _ bi3\a-_-lj _ (145) j=l On for the wing, tan panels lying 6. l is in the dzt dx So by plane 6 cos i equation of = the (141). Furthermore, wing, NW(dzti___ bij dx j=l lj that i = For non-coplanar used. 54 given cos 6 i [cos wing _ or dzc dx tail li - sin segments, _ cos (146) 8i] equation (145) must be Solution of the boundary condition equations.Several methods could be employed to solve the boundary condition equations for the unknown source and vortex strengths. For example, equation (139) could be solved by direct inversion, even though this is generally impractical for dense matrices of orders up to 1200. On the other hand, the partitioned matrix of equation (140) can be solved using the method described in reference i, which requires the inversion of only the diagonal partitions, having a maximum order of 600, together with matrix multiplications of the off-diagonal partitions. A rapidly convergent iteration scheme for solving large order systems of equations has been reported in reference 5. In this method, as applied in this report, the partitions are further subdivided into smaller blocks, with no block exceeding order 60. The matrix elements in each block are carefully chosen to represent some well defined feature of the original configuration. For example, a body block represents the influence of one ring of panels around the body, while a wing block represents the influence of one chordwise column of wing panels. For wings using the non-planar boundary condition option, the block size corresponds to the total number of panels on the upper and lower surface of the column. The initial iteration calculates the source and vortex strengths corresponding to each block in isolation. For this step, only the diagonal blocks are present in the aerodynamic matrix. Once the initial approximation to the source and vortex strengths is determined, the interference effect of each block on all the others is calculated by matrix multiplication. The incremental normal velocities obtained are subtracted from the normal velocities specified by the boundary conditions. This process is repeated 15 to 20 times, or until the residual interference velocities are small enough to ensure that convergence has occured. At present the computer program repeats the iteration a fixed number of times, namely 15. The consisting are To procedure of nine designated solve where yj is illustrated blocks. The , the specified Aij A.. z3 yj = = below unknown normal for an aerodynamic matrix singularity strengths velocities _ . i _.l AIi Al 2 AI 3 A21 A22 A23 A31 A32 A33 55 Put A = D + E D = Therefore [D + or First All 0 0 A22 0 0 E] A + A A2z 0 { Y} = {_ } {Y} = D-z[ ,_ A_ I = E - E{y Second Similarly, approximation k th {yl I : approximation: {._}k= D-_{_ - A_k-_ } Note'that 56 D -z = 0 32 }} A -z ii 0 0 0 A -z 22 0 0 0 A_3 1 1 A23 A A31 33 approximation: Calculate A 12 Calculation of Pressures, Forces, and Moments Once the strengths of the aerodynamic singularities have been determined, the three components of velocity at a point i can be determined as follows: where _u 1 , Av. 1 u.1 = _u. v.1 = _v.1 w. 1 = _w. 1 + are given and 6w. 1 l + cos _ (147) (148) sin _ (149) by equations (129) - (131) If the planar boundary condition option has been selected, the incremental velocity components induced by the wing thickness distribution must also be calculated and added to the above equations. the exact where The pressure coefficient isentropic formula Cp i _ YM 2-2 q2 i = u 2 i = 0, = 1 For M IE + 1 + Y-12 v 2 i + w2 i M2 (i is - ql) then 1 - calculated 1 using 1 (150) 2 CPi The forces - (151) qi and moments acting on the then be calculated by numerical integration. tangential force, and pitching moment about coordinates of panel i are given by: N.l =- Ti = AiCPiCOS AiCPi sin 8.cos l $i 6i configuration can The normal force, the origin of (152) (153) 57 M. = N.x. - T.z. l 1 i i i where A i 8 x is the panel (154) area l 8. l are the panel inclination by equation (124) l z. l are the point coordinates angles, of the panel defined control The total force and moment coefficients acting on the configuration are obtained by summing the panel forces and moments on both sides of the plane of symmetry. N 1 CN S 2Nl (155) 2T. i (156) 2M i (157) i=l N S1 CT _ i=l N CM 1 S_ - Z i=l Finally, the lift and C drag = C L C = The computer program acting on the body, the complete configuration. 58 be calculated output. cos _ - C N D moment may an optional coefficients C are: sin (158) cos (159) T sin _ + N C T computes the forces wing and tail surfaces, In addition, section for the wing and tail and moment and the forces and surfaces as COMPUTER Program PROGRAM Description A computer program has been developed to calculate the pressure distribution and aerodynamic characteristics of wingbody-tail combinations in subsonic and supersonic flow. The program is written in CDC FORTRAN IV, version 2.3 for a SCOPE 3.0 operating system and library file. It is designed for the CDC and eral 6000 series operates in disc files of computers, occupies 70,000 (octal) words, OVERLAY mode. The program requires five periphin addition to the input and output files. Program Structure The Figure calls SOLVE. CONFIG, while LINVEL, overlay structure of the program is illustrated on 2. The main overlay program is designated USSAERO, and the three primary overlay programs GEOM, VELCMP, and In turn, GEOM calls seven secondary overlay programs NEWORD, WNGPAN, NEWRAD, BODPAN, NUTORD, and TALPAN, VELCMP calls three secondary overlay programs BODVEL, and WNGVEL. The complete 19 subroutines. routine are given program Detailed in Part tions give the purpose the method used, and and constants. of list Operating The program deck and sequence: job card, system program deck, end-of-record file card. The input data section. consists of 14 overlay programs and descriptions of each program and subII of this report. These descripthe the program names of or subroutine, the principal outline variables Instructions data deck are loaded in the following control cards, end-of-record card, card, input data deck, and end-ofdeck is described in the following 59 OVERLAY(0,0) USSAERO OVERLAY(3,0 ) SOLVE OVERLAY(2,0) VELCMP OVERLAY( 1,0 ) GEOM INVERT PARTIN DIAGIN ITRATE PRESS FORMOM TRAP PANEL DERIV DERIVl DERIV2 CUBIC2 COMCU SCAMP4 OVERLAY(2,1) BODVEL SORPAN OVERLAY(2,2 ) LINVEL CONFIG OVERLAY (i, 1 ) SORVEL VORVEL NEWORD OVERLAY(I,2) WNGPAN OVERLAY(I,3) I I OVERLAY( 2,3 ) WNGVEL VORPAN TRANS NEWRAD OVERLAY(I,4) I BODPAN OVERLAY (i, 5 ) NUTORD OVERLAY (i, 6) 1 TALPAN OVERLAY(I,7) I Figure 6O 2 - Program Overlay Structure Program Input Data The input to this program consists of two basic parts, namely, the numerical description of the configuration geometry as described in reference 3, and an auxiliary data set specifying the singularity paneling scheme, program options, Mach number, and angle of attack. The program input is illustrated by the sample case presented in Appendix III. Description of configuration _eometry configuration is defined to be symmetrical therefore only one side of the configuration The convention used in this program is to the configuration located on the positive plane. The number of input cards depends ponents used to describe the configuration, detail used to describe each component. Card tifying 1 1-3 Identification.- information Card 2 each punched may be used lowing: Columns - in columns Card 1 input cards.The about the xz plane, need be described. present that half of y side of the xz on the number of comand the amount of contains any desired iden- 1-80. - Control integers.Card 2 contains 24 integers, right justified in a 3-column field. Columns 73-80 in any desired manner. Card 2 contains the fol- Variable J0 Value 0 1 4-6 Jl 0 1 -i 7-9 J2 0 1 -i Description No reference area Reference area to No wing Cambered Uncambered be read data wing data to be read wing data to be read No fuselage data Data for arbitrarily shaped fuselage to be read Data for circular fuselage read (With J6=0, fuselage be cambered. With J6=-l, lage will be symmetrical plane. uration With will J6=l, entire be symmetrical to be will fusewith xyconfigwith xy-plane) 10-12 J3 No pod (nacelle) data Pod (nacelle) data to be read 61 Columns Variable 13-15 J4 16-18 J5 19-21 J6 Value 0 1 NWAF No fin (vertical tail) data Fin (vertical tail) data to be read No canard (horizontal Canard (horizontal be read A cambered fuselage 2-20 Number to 25-27 NWAFOR if circular J2 is tail) data tail) data to or arbitrary nonzero Complete configuration is symmetrical with respect to xy-plane, which implies an uncambered circular fuselage if there is a fuselage Uncambered circular fuselage with J2 nonzero -i 22-24 Description 3-30 of describe Number of airfoil the sections used wing ordinates used to define each wing airfoil section. If the value of NWAFOR is input with a negative sign, the program will expect to read lower surface ordinates also 28-30 NFUS 1-4 31-33 NRADX(i) 3-30 34-36 NFORX(1 ) 2-30 Number NRADX(2 ) 3-30 40-42 NFORX(2 ) 2-30 43-45 NRADX( 3) 3-30 62 fuselage segments Number of half-section segment. the program number of points used to represent of first fuselage If fuselage is circular, computes the indicated yand z-ordinates Number stations lage 37-39 of of for first fuse- segment Same as fuselage NRADX(1), segment but for second Same as fuselage NFORX(1), segment but for second Same as fuselage NRADX(1), segment but for third Columns Variable Value Description 46-48 NFORX(3) 2-30 Same as NFORX(1), fuselage segment but for third 49-51 NRADX(4) 3-30 Same as NRADX(1), fuselage segment but for fourth 52-54 NFORX(4) 2-30 Same as NFORX(1), fuselage segment but for fourth 55-57 NP 0-9 Number of pods described 58-60 NPODOR 4-30 Number of stations at which radii are to be specified 61-63 NF 0-6 Number of fins to be described 64-66 NFINOR 3-10 Number of ordinates used to describe each fin (vertical tail) airfoil section 67-69 NCAN 0-2 Number of canards (horizontal tails) to be described 70-72 NCANOR 3-10 Number or ordinates used to define each canard (horizontal tail) airfoil section. If the value of NCANORis input with a negative sign, the program will expect to read lower surface ordinates also, otherwise the airfoil is assumed to be symmetrical (vertical pod tails) Cards 3, 4, . . . - remainin@ input data cards.The remaining input data cards contain a detailed description of each component of the configuration. Each card contains up to I0 values, each value punched in a 7-column field with a decimal point and may be identified in columns 73-80. The cards are arranged in the following order: reference area, wing data cards, fuselage data cards, pod data cards, fin (vertical tail) data cards, and canard (horizontal tail) data cards. in Reference columns 1-7 tains Wing the area and card: may be data cards: locations in The reference identified as The first wing percent chord area REFA value is in columns punched 73-80. data card (or cards) at which the ordinates conof 63 all the wing airfoils ly NWAFORlocations identified the last are to be specified. in percent chord given. in columns location in 73-80 percent by the chord symbol given There will Each card be exact- XAFJ where J on that card. may be denotes The next wing data cards (there will be NWAF cards) each contain four numbers which give the origin and chord length of each of the wing airfoils that is to be specified. The card representing the most inboard airfoil is given first, followed by the cards for successive airfoils. These cards contain the following: Columns Contents 1-7 x-ordinate of airfoil leading edge 8-14 y-ordinate of airfoil leading edge 15-21 z-ordinate of airfoil leading edge 22-28 airfoil 73-80 card identification, WAFORGJ where J denotes the paricular airfoil, thus WAFORGI denotes the most inboard airfoil data values If a cards of cambered is the delta z streamwise chord length wing has been specified, the next set of wing mean camber line cards. There will be NWAFOR referenced to the z-ordinate of the airfoil leading edge, each value corresponding to a specified percent chord location on the airfoil. These cards are arranged in the order which begins with the most inboard airfoil and proceeds outboard. Each card may be identified in columns 73-80 as TZORDJ where J denotes the particular airfoil. Note that the z-ordinates are dimensional. Next are the wing ordinate cards. There will be NWAFOR values of half-thickness specified for each airfoil expressed as percent chord. These cards are arranged in the order which begins with the most inboard airfoil and proceeds outboard. Each card may be identified in columns 73-80 as WAFORDJ where J denotes the particular airfoil. Fuselage the x There in of 64 data cards: values of the fuselage will be NFORX(1) values columns the last 73-80 by fuselage The first stations and the card (or cards) of the first cards may be the symbol XFUSJ where station given on that J denotes card. specifies segment. identified the number If the fuselage is circular, the next card (or cards) gives the fuselage cross sectional areas, and may be identified in columns 73-80 by the symbol FUSARDJ where J denotes the number of the last fuselage station given on that card. If the fuselage is of arbitrary shape, NRADX(1) values of the y-ordinates for a half-section are given and identified in columns 73-80 as YJ where J is the station number. Following the y-ordinates are the NRADX(1) values of the corresponding z-ordinates for the half-section identified in columns 73-80 as ZJ where J is the station number. Each station will have a set of y and z, and the convention of ordering the ordinates from bottom to top is observed. For each fuselage segment a new set of cards as described must be provided. The segment descriptions should be given in order of increasing values of x. Pod data cards: The first the location of the origin fies contains the pod (nacelle) of the first data pod. card speciThe card following: Columns Contents 1-7 x-ordinate of origin of first pod 8-14 y-ordinate of origin of first pod 15-21 z-ordinate of origin of first pod 73-80 card where The ordinates, of be the zero next pod input data referenced to the pod radii and the cards may where J For must gram cate (or cards) origin, at to be specified. x value is the identified in denotes the pod each additional columns of Fin data describe a contains the zero The length 73-80 PODORGJ pod number contains the xwhich NPODOR values first x value of the pod. by the symbol must These XPODJ number. pod, new be provided. Only single pods assumes that if the y-ordinate is located symmetrically with y-ordinate to be are last card pod identification, J denotes the implies cards: Exactly fin (vertical a PODORG, XPOD, and PODR cards are described but the prois not zero an exact duplirespect to the xz-plane, a single pod. three tail). data The input first cards are used fin data card following: 65 Columns Contents 1-7 x-ordinate on inboard leading edge airfoil 8-14 y-ordinate leading 15-21 z-ordinate leading chord 29-35 x-ordinate 43-49 airfoil of inboard airfoil length leading edge y-ordinate leading edge z-ordinate leading x inboard edge 22-28 36-42 of edge of inboard of outboard airfoil of outboard airfoil of outboard airfoil edge 50-56 chord 73-80 card identification, J denotes the fin The expressed second in length of outboard J The card denotes The third fin input data the fin airfoil half-thickness Since the fin airfoil must be on the positive y side of the The card identification FINORDJ where J denotes the fin number. may the be identified fin number. airfoil FINORGJ number fin input data card contains NFINOR percent chord at which the fin airfoil are to be specified. 73-80 as XFINJ where airfoil in where values of ordinates columns card contains NFINOR values of expressed in percent chord. symmetrical, only the ordinates fin chord plane are specified. may be given in columns 73-80 For each fin, new FINORG, XFIN, and FINORD cards must be provided. Only single fins are described but the program assumes that if the y-ordinate is not zero an exact duplicate is located symmetrically with respect to the xz-plane, a yordinate of zero implies a single fin. Canard data cards: airfoil is symmetrical, a canard, and the input fin. If, however, the 66 If the canard (or exactly three cards is given in the same canard airfoil is not horizontal tail) are used to describe manner as for a symmetrical (indicated by a negative value of NCANOR), a fourth canard input data card will be required to give the lower ordinates. The information presented on the first canard input data card is as follows: Columns Contents 1-7 x-ordinate of inboard leading edge airfoil 8-14 y-ordinate of inboard leading edge airfoil 15-21 z-ordinate of inboard leading edge airfoil 22-28 chord 29-35 x-ordinate of outboard leading edge airfoil 36-42 y-ordinate of outboard leading edge airfoil 43-49 z-ordinate of outboard leading edge airfoil 50-56 chord 73-80 card where The second canard input x expressed in percent chord nates are to be specified. columns 73-80 as XCANJ where the card. The be punched must may the the be identified canard number. lower ordinates length of identification, J canard, new airfoil of outboard denotes the in card contains expressed in columns 73-80 If the canard are presented airfoil CANORGJ canard CANORG, XCAN, NCANOR percent number values chord. as CANORDJ airfoil is on a second program expects both upper and lower as positive values in percent chord. For another be provided. inboard data card contains NCANOR values at which the canard airfoil ordiThe card may be identified in J denotes the canard number. The third canard input data canard airfoil half-thickness This card J denotes metrical, length and of where not symCANORD ordinates CANORD of to cards 67 Description Card identifying Card planar Non of Auxiliary i.i - Identification.information in columns Input Cards Card i.i 1-80. 1.2 - Boundary condition boundary conditions are contains and control point always applied on however card 1.2 permits the selection of boundary apply on a wing, fin (vertical tail), or canard This card also selects the output print options. tains the following: Columns Variable 1-3 LINBC Value 0 any desired definition.a body, conditions (horizontal This card to tail). con- Description Control points on surface fin (vertical tail), and (horizontal tail). This ferred to as ary condition of wing, canard is re- the nonplaner option. bound- Control points in plane of wing, fin (vertical tail), and canard (horizontal tail). This is referred to ary condition 4-6 THICK 7-9 PRINT 0 Do not matrix 1 Calculate if LINBC as the planar option. calculate 0 Print forces 1 Print out wise loads canards = wing 1 wing out the pressures and moments option on the 0 negative output LINBC, integers. 68 value indicated of and the the spanfins, and and the velocsource and the steps Print the axial print adds for options THICK, and PRINT THICK is not used and wing, matrix vortex strengths Print out option 2 and the iterative solution and A the thickness thickness Print out option 1 ity components and to bound- out normal the 1-4. are punched if LINBC = option 3 velocity panel as 0. and matrices geometry right in print justified out Card control right Columns 1-3 4-6 2.1 - integers.justified Variable K0 K1 Revised configuration The integers K2 10-12 K3 13-15 K4 card 2.1 Value No reference Reference 0 1 No wing data Wing data to sharp leading Wing data to round leading 0 1 3 punched as lengths length data to be read, edge be read, edge be read wing has a wing has a No body data Body data follows Not 0 1 description are Description 0 1 3 7-9 paneling contents of as follows: used No fin (vertical tail) data Fin (vertical tail) data to read, fin has a sharp leading be edge Fin (vertical read, fin has be tail) data to a round leading edge 16 -18 K5 0 1 3 No canard (horizontal tail) data Canard (horizontal tail) data to be read, canard has a sharp leading edge Canard (horizontal tail) data to be read, leading 19-21 K6 Not 22-24 KWAF Number define canard edge KWAFOR a round used of wing sections the inboard and panel edges. panel edges the geometry 25-27 has Number of define the edges KWAFOR defined of = If KWAF are defined input ordinates leading used to outboard = used and 0, by the NWAF in to trailing the wing panels. If 0, the panel edges are by NWAFOR in the geometry input 69 Columns Variable Value Description 28-30 KFUS 31-33 KRADX(i) 0, 3-20 Number of meridian lines used to define panel edges on first body segment. There are three options for defining the panel edges. If KRADX(1) = 0, the meridian lines are defined by NRADX(1) in the geometry input. If KRADX(1) is positive, the meridian lines are calculated at KRADX(1) equally spaced PHIKs. If KRADX(1) is negative, the meridian lines are calculated at specified values of PHIK 34-36 KFORX(i) 0, 2-30 Number of axial stations used to define leading and trailing edges of panels on first body segment. If KFORX(1) = 0, the panel edges are defined by NFORX(1) in the geometry input 37-39 KRADX(2) 0, 3-20 Same as KRADX(1), body segment but for second 40-42 KFORX(2) 0, 2-30 Same as KFORX(1), body segment but for second 43-45 KRADX(3) 0, 3-20 Same as KRADX(1), body segment but for third 46-48 KFORX(3) 0, 2-30 Same as KFORX(1), body segment but for third 49-51 KRADX(4) 0, 3-20 Same as KRADX(1), body segment but for fourth 52-54 KFORX(4) 0, 2-30 Same as KFORX(1), body segment but for fourth The number of fuselage segments. The program sets KFUS = NFUS The program is restricted to 600 body singularity panels. For this program there is an additional restriction that the total number of singularity panels in the axial direction on the body (fuselage) cannot exceed 30. The arbitrary body (fuselage) capability of this program is limited to those shapes for which the radius is a single-valued function of PHIK for each cross section of the body. 7O Card description punched as Columns 1-3 2.2 - Additional revised configuration control integers.The right justified integers Variable 0, 2-20 card 2.2 are Description Value KF(1) paneling contents of as follows: Number define of the fin sections inboard and used to outboard panel edges on the first fin. If KF(1) = 0, the root and tip chords define the panel edges 4-6 7-9 i0-12 13-15 16-18 19-21 22-24 25-27 28-30 31-33 34-36 KFINOR (i) KF (2) KFINOR 0, 2-20 (2) KF(3) KFINOR KF (3) (4) 0, 3-30 0, 2-20 (5) KF(6) KFINOR 0, 3-30 0, 2-20 KF (5) KFINOR 0, 3-30 0, 2-20 (4) KFINOR 0, 3-30 0, 3-30 0, 2-20 (6) 0, 3-30 Number define of ordinates the leading used to and trailing edges first of the fin. fin panels If KFINOR(1) on = the 0, the panel edges are by NFINOR for second but for for third but for for fourth but for for fifth but for for sixth but for Same as defined for KF(1), for KFINOR(1), but fin Same as second Same fin as for KF(1), as for KFINOR(1), but fin Same third Same fin as for KF(1), as for KFINOR(1), but fin Same fourth Same fin as for KF(1), as for KFINOR(1), but fin Same fifth Same fin as for KF(1), as for KFINOR(1), but fin Same sixth fin 71 Columns Variable 37-39 KCAN( 1 ) Value Description Number define of canard the inboard sections used and outboard to panel edges on the first canard. If KCAN(1) = 0, the root tip chords define the panel edges. If KCAN(N) negative, no vortex sheets carry through the body and concentrated vortices are shed from the inboard canard or tail 40-42 KCANOR(1 ) Number the the of leading first the panel NCANOR 43-45 KCAN(2 ) 46-48 KCANOR(2 ) KCAN(3) 52-54 KCANOR(3) , , 3-30 55-57 KCAN(4) 0 , 2-20 58-60 KCANOR( 4) 61-63 KCAN(5 ) 0 , KCANOR(5) , 3-30 67-69 KCAN(6) 0 ! 2-20 70-72 KCANOR(6) 0! 3-30 72 of the used to define and canard. trailing edges of If KCANOR(1)=0, edges are defined but Same as for KCANOR(1), second canard Same third as for canard KCAN(1), Same third as for canard KCANOR(1), Same as for KCAN(1), fourth canard Same as for KCANOR(1), fourth canard 2-20 64-66 ordinates Same as for KCAN(1), second canard 3-30 49-51 edge surface Same fifth as for canard KCAN(1), Same fifth as for canard KCANOR(1), Same sixth as for canard KCAN(1), Same sixth as for canard KCANOR(1), by for but but for for but but for for but but for for but but for for but for panels The program is restricted on the wing-fin-canard to a total combination. of 600 singularity For this program there is an additional restriction the total number of singularity panels in the spanwise on the wing-fin-canard combination cannot exceed 20. that direction Cards 3, 4, . . . - remainin_ input data cards.The remaining input data cards contain a detailed description of the singularity paneling of each component of the configuration. Each card contains up to i0 values, each value punched in a 7-column field with a decimal point and may be identified in columns 73-80. The cards are arranged in the following order: reference lengths, wing data cards, fin (vertical tail) data cards, canard (horizontal tail) data cards, fuselage (body) data cards, and finally Mach number and angle of attack case cards. Note that the present program will not handle a pod and therefore there are no pod panel inputs. However, if the geometry input contains a pod description it will be read and ignored. in Reference length columns 73-80 and Columns card: contains This the card may following: Variable be identified as REFL Description i- 7 REFA Wing reference the reference the value of input 8-14 REFB Wing value semispan. of 1.0 reference area. If REFA = 0, area is defined by REFA in the geometry If REFB = 0, used for the is a semispan 15-21 REFC Wing reference a value of 1.0 reference chord 22-28 REFD Body (fuselage) reference diameter. If REFD = 0, a value of 1.0 is used for the reference diameter 29-35 REFL Body (fuselage) reference length. If REFL = 0, a value of 1.0 is used for the reference length 36-42 REFX X coordinate of moment center 43-49 REFZ Z coordinate of moment center chord. is used If REFC for the = 0, 73 Wing data cards: The first wing data card is the wing leading edge radius card and is required only when K1 = 3. This card contains NWAFvalues of leading edge radius expressed in percent chord. It may be identified in columns 73-80 as RHOJ where J denotes the number of the last radius given on that card. Next is the wing panel leading edge card. This card contains KWAFORvalues of wing panel leading edge locations expressed in percent chord. This card may be identified in columns 73-80 as XAFKJ where J denotes the last location in percent chord given on that card. Omit if KWAFOR= 0. The last wing data card gives the wing panel side edge data. This card contains KWAF values of the y ordinate of the panel inboard edges. This card may be identified in columns 73-80 as YKJ where J denotes the last y ordinate on that card. These values are arranged in the order which begins with the most inboard panel edge and proceeds outboard. Omit if KWAF = 0. Fin (vertical tail) data cards: The first fin data card is the fin leading edge radius card and is required only when K4 = 3. This card contains NF values of leading edge radius expressed in percent chord, one value for each fin. It may be identified in columns 73-80 as RHOFIN. Next is the fin panel leading edge card for the first fin. This card contains KFINOR(1) values of fin panel leading edge locations expressed in percent chord. This card may be identified in columns 73-80 as XFINKJ where J denotes the fin number. Repeat this card for each fin. The last fin data card gives the fin panel side edge data for the first fin. This card contains KF(1) values of the z ordinate of the panel inboard edges. This card may be identified in columns 73-80 as ZFINKJ where J denotes the fin number. These values are arranged in the order that begins with the most inboard panel edge and proceeds outboard. Repeat this card for each fin. Canard (horizontal tail) data cards: The first canard data card is the canard leading edge radius card and is required only when K5 = 3. This card contains NCAN values of leading edge radius expressed in percent chord, one value for each canard. It may be identified in columns 73-80 as RHOCAN. Next is the canard panel leading edge card for the canard. This card contains KCANOR(1) values of canard leading may be edge locations identified in canard number. 74 Repeat expressed in percent columns 73-80 as XCANKJ this card for each chord. where canard. J first panel This card denotes the The data for of the y identified number. the most card for last canard data card gives the canard panel side edge the first canard. This card contains KCAN(1) values ordinate of the panel inboard edges. This card may be in columns 73-80 as YCANKJ where J denotes the canard These values are arranged in the order that begins with inboard panel edge and proceeds outboard. Repeat this each canard. Fuselage (body) data cards: The first body card is the body meridian angle card. This card contains KRADX(1) values of body meridian angle expressed in degrees and may be identified in columns 73-80 as PHIKJ where J denotes the body segment number. The convention is observed that PHIK = 0. at the bottom of the body and PHIK = 180. at the top of the body. Omit unless KRADX(1) is negative. Repeat this card for each fuselage segment. The second body card is the body axial station card. This card contains KFORX(1) values of the x ordinate of the body axial stations and may be identified in columns 73-80 as XFUSKJ where J denotes the body segment number. Omit if KFORX(1) = 0. Repeat this Mach identified card for number and in columns Columns each angle 73-80 fuselage segment. of attack as MALPHA card: and Variable i- 7 This contains card may be the following: Description MACH The subsonic Mach number (includ- ing the value MACH = 0.) or the supersonic Mach number at which it is desired to calculate the aerodynamic 8-14 ALPHA A for with series of Mach The angle of attack expressed degrees at which it is desired calculate the aerodynamic data number the same geometry may the desired values. A of the such a value present terminal of MACH case. card. = data -i. and be on Geometry angle of calculated this cards attack by card for combinations repeating signifies a new in to this the case card termination can follow 75 Program Output Data All output is processed by a standard 132 characters-perline printer. The output from each run is always preceded by a complete list of the input data cards. The amount and type of the remaining output depend on the PRINT option selected, the number of panels used, and whether the configuration being analyzed is an isolated wing, an isolated body, or a complete wing-body-tail combination. The program output options are described below: PRINT = 0 The program prints the case description, Mach number and angle of attack, followed by a table listing the panel number, control point coordinates (both dimensional and non-dimensional), pressure coefficient, normal force, axial force, and pitching moment. Separate tables are printed for the body and wing panels, noting that any tail, fin or canard panels are included with the wing output. If the planar boundary condition option has been selected, the results for the wing upper surface are given in one table, followed by a separate table giving the results for the wing lower surface. Additional tables giving the total coefficients on the body, the wing and the complete configuration follow the pressure coefficient tables. These include the reference area, reference span and reference chord, the normal force, axial force, pitching moment, lift, and drag coefficients, and the center of pressure of the component. PRINT = 1 In addition to the output described for PRINT = 0, the program prints out additional tables giving the normal force, axial force, pitching moment, lift and drag coefficients, and the center of pressure of each column of panels on the wing and tail surfaces. In addition, the indices of the first and last panel in the column are listed, together with the span, chord and origin of the column. PRINT = 2 In addition to the output described for PRINT = i, the program prints out tables listing the panel number, the source or vortex strength of that panel, and the axial velocity u, lateral velocity v, and vertical velocity w at the panel control point. The normal velocity is also calculated for 76 body panels. Separate tables are printed for the body and wing panels, noting again that any tail, fin, or canard panels are included with the wing output. If the planar boundary condition option has been selected, separate tables are given for the wing upper and lower surfaces. PRINT = 3 In addition the the at PRINT = 4 to the program source each step In addition the program normal ments prints number output described of the to the prints iterative solution by number printing PRINT = 2, and number, obtained procedure. output described for PRINT = out tables of the axial and velocity components which of the aerodynamic matrices. out the matrix row number, of elements in that row. four matrix option is for prints out the iteration and vortex strength arrays partitions selected, each will of 3, make up the eleThe program and gives the A maximum of be printed which is and its influence description the velocity component tables. if this identified prior to If a negative value of PRINT is selected, the program prints all the information described above for the positive values, together with the complete panel geometry description of the configuration following the list of input cards. This consists of tables giving the wing panel corner points, control points, inclination angles, areas, and chords. If the configuration has a horizontal tail, fin or canard, additional tables are printed giving the same information as listed above for the wing. panel angles Finally, if corner points, are listed. The sented in program Appendix the configuration control points, output III. is illustrated includes areas, by a and the body, the inclination sample case body pre- 77 EXPERIMENTAL VERIFICATION Several examples of pressure distributions calculated by the program are presented in this section, and compared with experimental data. The examples include isolated bodies, isolated wings, and wing-body combinations in both subsonic and supersonic flow. Isolated Bodies O_ive-cylinder body with boattail in subsonic flow.The theoretical pressure distribution calculated for this body at M = .40 and _ = 0 degrees is presented on Figure 3. The experimental data has been obtained from reference 6, which also contains experiment the blunt 12 degree additional for this base and cone aft comparisons between the present theory and body at M = .61 and .83. In this example, sting were replaced by an arbitrarily chosen of the boattail region in an attempt to sim- ulate the flow separation ment between the theory of the body. retical having Haack-Adams pressure A body with distribution /A base and = .532 region behind experiment base and in supersonic calculated for _/d max the body. is achieved = i0 is a Good over agreemost flow.The Haack-Adams presented on theobody Figure max 4, for M = 2.01 and e = 0 degrees. The experimental data for this body is obtained from reference 7, which also gives the pressure distribution calculated by characteristics theory. The present method agrees closely with the experimental data and the characteristics theory for this body. Elliptic cone in supersonic flow.The theoretical pressure distribution on an elliptic cone is compared with experimental data on Figure 5, for M = 1.89 and e = 0 and 6 degrees. The experimental data was obtained from reference 8. Again, the theory agrees well with experiment except near the leading edge on the lower over-estimated. 78 surface, where the positive pressure is slightly I,M E' E', O m I-,4 Z H o 4-1 Ill 1 • I -,4 I • O i • II 1i • U I 0 • _-t \ 0 I • I ! • I "0 0 D -;-I m -,-4 N 0 0 -,-4 -,--4 4-1 m -,-..t 1.4 m m ill -,-i u? 0 m ..Q I (9 _4 -,-4 79 0 _r_ ,< I 0 80 o o II o '_ II ,< N r_ o 1,4 0 I .,.4 1,4 I o I 0 o L; o 0 o I I 0 0 -r-I 0 -,--t N ,< 0 0 -,.-t -,-'t I/1 .,-I 0 m m 0 _4 U -,-I 0 m 14 0 ! "1"t r_ Z 0 U H o_ II _J \ \_.___//I \V i o I % II 1.4 0 © © -D I r-t U -H X © ,-I o D _5 II o o I _ ._ D .J o o o 0o o ,-I o _I" r-I 0 0 0 ,--I o o o cxl 0 U -,-I ._I 0 0 .¢-I .H ._I b_ 1.4 C) -_I 0 I -_I I 81 Isolated Two-dimensional distribution on a degrees is compared airfoil NACA 64A010 with the Wings in subsonic airfoil at experimental flow.The pressure M = .167 and e = 8 data from reference on Figure option was 6. In this example, the surface boundary utilized in the theoretical calculations. agreement surface with the of experiment airfoil, is excellent indicating that on both viscous condition The upper and lower effects are small. The potential flow solution obtained by the present method also agrees closely with that given by the viscous flow solution presented in reference i0 for this airfoil, except for a small region near the trailing edge. In general, potential flow theory tends to over-estimate the negative pressure peaks in two-dimensional flow. Figure 7 compares the results of the present program with the exact incompressible pressure distribution around a 10 percent thick Karman-Trefftz airfoil. Here, the program results agree closely with the exact solution, and give considerable confidence in the capability of the present method to reproduce theoretical two-dimensional flows. Variable sweep butions calculated section, 72 degrees angles, are compared win@ in subsonic flow.on a variable sweep wing inboard sweep, and two with experimental data The pressure having a NACA outboard wing from reference distri64A006 sweep ii at M = .23 and e = 10.5 degrees on Figure 8. In this example, the boundary conditions are applied in the plane of the wing. The theory agrees reasonably well with experiment at the root and at the mid-span break point. The pressure distribution is less accurate near the wing tip, although the net loading appears to be approximately correct. The agreement between and experiment is considered to be acceptable, considering relatively high angle of attack chosen for this comparison. theory the Cambered arrow win_ in bu_ersonic flow.The pressure distributions calculated on a cambered and twisted arrow wing having a 3 percent circular arc section and 70 degrees sweepback are compared with experimental data from reference 12 at M = 2.01 and e = 4 degrees on Figure conditions are applied in the plane can be seen to agree reasonably well entire wing, except in the immediate edge. 82 9. Here, the boundary of the wing; and the theory with experiment over the vicinity of the leading 9 < NACA 64A010 AIRFOIL 6.0 M = .167 e = 8° R = 4.1 × I0 _ 5.0 Theory Q Experiment 4.0 - Cp 3.0 2.0 1.0 0.0 - 1.0 Figilre 6 - Pressure Pistribution on Two-dimensional Airfoil 83 KARMAN-TREFFTZAIRFOIL 2.0 t/c = 10% M = 0. 2.5% Camber 1.6 -Cp _ = 4° Exact 1.2 O Present Theory Method .8 .4 --°4 8 Figure 84 7 - Incompressible Pressure Distribution on Two-dimensional Airfoil )4 0 E_ 0 0 0 0 u_ $4 C 0 $4 0 0 0 4-) -,4 C O C_ C O C ,--4 C_ O M U_ O C O X [] O 11 o uD II ! C_ 7 "0 0 0 o o 0 ;) r_ ml 0 0 o ,.-I C 0-,._ ,-Io ._< < -,-I kO >u gZ C_ O C C O 4J • • C •,40 D m ._ cq CrO _._ I co 85 86 o II ,-I o II 44 X 0 0 @ O 0 o\ \\ \\ ! \ I i t::::)_ (D G) ( E (D (D -,-I ra_ _r4 o O_ 0 c_ rd I-4 0 -,.-I •,..-Ir_ 4..1 -,.-I .14 -,-I _-_ m m_ -,--I O_ mr..) _4 q) .1::_ I _4 .,-t Wing-Body Combinations O@ive-cylinder body with swept win 9 in supersonic flow.The planform of this simple wing-body combination, and the paneling scheme used in the aerodynamic representation are shown on Figure i0. The wing has a NACA 65A004 section, is centrally mounted on the body, and the quarter-chord line is swept back 45 degrees. The ogival nose is 3.5 body diameters in length. The pressure distribution on the wing is compared with experimental data from reference 13 for M = 2.01 and _ = 5 degrees on Figure ii. The theoretical curve was calculated using the planar boundary condition option. The agreement between theory and experiment is reasonably good, except near the wing leading edge. Pa:t of this discrepancy is probably due to shock wave detachment ahead of the round leading edge of the airfoil for this supersonic leading edge circulation to develop around the leading calculations assume an attached Mach wave leading edges, prohibiting the development flow. The lines 12. iment pressure distribution on of the body are compared with In this example the agreement is extremely good. the wing, allowing a small edge. The theoretical along supersonic of any circulatory upper and experimental between theory lower data and meridian on Figure exper- 87 88 / / 0 I -,-I f-I .;-I r_ 0 0 .r-I C_ -,-I I 0 o t_ -_1 0 I • 0 I I ! I I 0 Q Q |e ,--I I I• f_l I !• ("7. I I 'lJ O -,-I (D CL q-q O -H _D ,.el Z_ O O -,--I C_ © (D 0-, {.) O (D I (D Zl 89 o II L_ ,-t o II Z 0 H H H ! m 0 9O 0 I I o • ,a r-i i" 0 ! I o_ N Q / | I U I • o,-I I I 0"_ I • _1 I ,--I I I I • o_ -r-I I+ 0 ,-a o i ,,-i ,-i -,-t u,-t 0 0 0 0 -,-t m -,--I m _4 0 -,-4 0 m I1) I r-t -,.-t CONCLUSIONS The aerodynamic analysis method described in this report has been developed to succeed the earlier methods reported in references 1 and 2. Considerable progress has been made in the achievement of this goal, but additional work remains if the full potential of the new method is to be realized. Several promising described areas briefly for the below. future development of this method are Program refinements.Increased geometrical capability is required to permit the analysis of engine pods, nacelles, or fairings. In addition, improved programming techniques, including far field approximations to the aerodynamic singularity distributions would be desirable to reduce the time required to calculate the matrix of aerodynamic influence coefficients. Finally, the extension of the force and moment subroutine to include the calculation of additional aerodynamic stability derivatives would be valuable. Development of computer program is sure distribution, The development of increase its range problem of determining which will generate the presence of an drag and interference based on the design the pro@ram as a design tool.The present restricted to the determination of the presforces and moments on given configurations. the program as a design tool would greatly of application. For example, the important the wing camber and twist distribution favorable surface pressure distributions in arbitrary body, or for minimizing pressure effects could be included in this program procedures described in reference i. Development of leading ed@e vortex model.The use of linearly varying vortex distributions to represent the circulatory flow around lifting wings permits the Kutta condition to be imposed along leading edges, as well as trailing edges. Using this flow model, the vortex sheet from leading edge panels can be modified to trail downstream from the leading edge or wing tip to simulate a separated flow, and provide a first approximation to the lift distribution on wings at high angles of attack. Application of the method to the analysis of transonic flows.The non-linear effects of transonic flow can be approximated by using the local Mach nun_er calculated from the potential flow solution to redefine the regions of influence and the magnitude of the velocity field induced by the aerodynamic singularities. An iterative solution of the boundary condition 91 equations can then be established, in which the coefficients of the equations depend on the local Mach number distribution of the preceding step. The iterative procedure would be continued until a convergent pressure distribution is obtained. Analytical Bellevue, December 92 Methods, Washington 31, 1972 Incorporated APPENDIX I Integration The velocity may all tables, component be reduced to forms and two non-standard J (v2 = The + e2)[av [ J2 method i. integrals appearing integrals used to results dv 2 (V 2 + e2)[_v 2 evaluate these are summarized J2 Y_ - y_ + b2e 2 y_ + b2e 2 a real - c)y 2 + 2bv + c]½ = tan tanh- - b2e 2 2 i e[av sinh - be root y[av -I = y2F non-zero and G c]_J is given in refer- _2 = (ae 2 F + = z is 2bv integrals below: Y Y + vdv J J where appearing in the text in standard integral given below: = 1 The ence Procedures 2 of = + G ] 0, 2by 72 - bv vy + e_ + 2bv (ae 2 -y2)v2 equation _ be/y = + C]½ + C]½ vY -x the + + e_ (c - _a)ea]½ 93 APPENDIX II Panel Geometry Calculation Procedure The analytical procedure presented here follows closely the method first developed in reference 14. A quadrilateral surface element is described by four corner points, not necessarily lying in the same plane, as shown in the sketch. Note that the numbering convention of the corner points differs from that used in the preceding text. The quadrilateral element is approximated by a planar panel as follows: 2 I \J / / The fied coordinates by their vectors We may cross T in the subscripts. and 1 T 2 reference The of system are the diagonal identi- are T ix = x 3 - x i T ly = Y 3 - Yl T iz = z3 - z T2x = x _ - x2 T 2y = Y 4 - Y2 T 2z = z_ - z2 by taking of the now product obtain a vector N (and diagonal its components) vectors. N 94 coordinate components = T 2 × T 1 the N = T T - T T x 2y Iz ly 2z N = T T y Ix N = T as N unit unit normal divided normal). by - 2x ly vector, n, to its length own T T 2x T z The 2z T T Ix the lZ _y plane N of the (direction element cosines is of taken outward Nx nx - N ny N NZ where nz - N N = [ N 2 x + N 2 y + N2z]½ The plane of the element is now completely determined if a point in this plane is specified. This point is taken as the point whose coordinates, _, y, z are the averages of the coordinates of the four input points. 1 Y = 41 [Yl l[z Now the element input along points will the normal + x2 + x + x ] + Y2 + Y3 + y ] +z +z be projected vector. The +z] into the resulting plane points of the are the 95 corner are points of equidistant d = = quadrilateral the plane, Inx(X The coordinates ordinate system x' k the from - x I) + ny(y of the corner are given by x + (I) - k+* n k and - element. The this distance yl) points + in nz(Z the - input is points zl)I reference co- d x k+l | Now the requires vectors, Yk = Yk + (-1) z' k = z + (i) - k nYd k+1 n k = i, 2, 3, 4 d z element coordinate the components of one of which points system must be constructed. three mutually perpendicular along each of the coordinate This unit axis of the system, and also the coordinates of the origin of the coordinate system. All these quantities must be given in terms of the reference coordinate system. The unit normal vector is taken as one of the unit vectors, so two perpendicular unit vectors unit in vectors the plane t of and t 1 by its own length the element The 2" T 1 t , i.e., ix TIX T t 1 TIY ly T 1 TIZ lZ T 1 96 t needed. is 1 - t vector are taken Denote as T these divided 1 Ti where The are vector t2 is = [T ix2 +T ,y2 +T_]% lZ defined t by = n 2x = n 2y = The vector t is × - n - that its components t n lZ t y vector so ly x ly t I, t z t unit n ix x the - t n 2z n Iz z t = t y t t ix parallel to the x or _ axis 1 of the element coordinate system, while t is parallel to the 2 y or _ ordinate The axis, and system. corner ordinate points n points is parallel are now system based have coordinates system. with this on to plane of the element the element, coordinate which defines is a multiple coordinate n in the by _ they into a of axis the 0. of this element co- co- as origin. reference These coordinate coordinate Because system they lie zero Also, the x or _ axis of the "diagonal" and the coordinate _ point in the have system. or element , _ , k k z or because _ the in coordinate vector the element coordinate vector from point 1 to n are equal. In the I coordinate z transformed the average x', y', z' k k k Their coordinates origin are denoted the the in t I , system, 3, the (_,n) 3 system, the corner points of the element are: ! _k = n = t k t ix (x - x.') K + t ly (_- x') + t 2x k (y - yk ) + t _z (F- y') + t 2y k (z - z_) (_- z') 2z k 97 These corner quadrilateral I points are taken as the corners as illustrated in the following of a plane sketch. (_3 ' 1 n3) A (g_,, n 4) The origin of the element coordinate system is to the centroid of the area of the quadrilateral. average point as origin the coordinates of the element coordinate system are: _ o - 1 3 1 - T] T] [_ ,+ (q , - q 2 ) + _ 2 (D now transferred With the centroid in the _ - q 1 ) ] 2 1 T] 0 3 __ These are subtracted from the coordinates of the corner points in the element coordinate system based on the average point as origin to obtain the coordinates of the corner points in the element coordinate system based on the centroid as origin. Accordingly, these latter coordinates are _k = _k - _o k nk 98 = nk - n o = I, 2, 3, 4 Since the centroid is element, its coordinates are required. These the in =x+t Y0 = y + tl y_ 0 = z + t 0 area A IX of = the lZ the 1 -- (_ 2 3 _ _ 0 coordinate 2 0 + t2y_ 0 + t2z_ of the system 0 quadrilateral - point are +tx_ 0 control reference coordinates X0 z Finally, to be used as the _ )(_ l is 2 _ ) 4 99 APPENDIX III SAMPLE CASE i00 UN|FJEO SUbSUN|L-SUPE_SUNIC LLS| AckODYhAM[CS UF INPU[ P_UG&A_ VEgSIUN AO0 _AkDS _U_11L_x1IZZ_2_Z_3J_333_3_k_4_55_5_5566_06_77_77?_d _1_9_3_5678_J_9_L_5_9_3_56_89_34_6789_|_3456789_3_56_89_ UGAVE -1 CYLJNOER O, _. .5 _0. ,15 _5. I*_5 6u. 2*5 _5. 5= 50. 75* BO. 85. 90. 95. lO0o 27.65 1£, O, Z, O, 1.1925 ._Ol_ L,89_5 ,3_ 1,966 ,_1_ L.997> ,85_ 1,994 ,8745 L,9675 I,I?_ .969 ,7k5 ._SU ,2465 °009 O. 1.1925 ,_O?b 1.8955 .373 1.966 ._155 1o9975 .6_L5 1.996 ,8765 1,_475 L.176 .9_9 .715 ._SU .266_ ,009 O. b.8333 .50J3 o,4ibZ L, 1867 1.G Lo15 7o5B33 2o33_3 _.1667 11.86611_._ 18.75 2_._ 25.85 O, .08805 6o9583_5°595olo, 8.72l 8.7_1 SJNGULA&]I¥ .3_573 ,69Z_1 l,[8095|,109202,31.5622.980803,627136,ZgZ59 i91178ol_Z69?.gb2_tT.bB9558.o_555u.365818.>5o|08o886_i 8.727 8.727 8.121 8ol27 8.721 PANELING FOR SAHPLE CASk 0 1 -L 3 -L 8UUY 2 I wiTH _6 L 45 5 UEGkEE Z8 SwEeP 7*5 55. IO, bO. HIO-WING 15, 65° 20. 70. XAFI XAFZ _AFO_G2 6 lI I 5 I,06 1.B57 /o21O 1.728 1,_63 1,5615 I,o50_ 1,3815 _AFORDI _AFURO2 _AFORD3 l. O6 1.857 [.ZI6 1.128 Lo¢63 1,5615 1o6_05 1.3815 Z.9167 8.75 3.5 9.3333 _o0_33 9.9L61 6.8887 10.5 5._5 11.0833 _AFOkUI _AFURO2 29._ 32.95 36,5 XFUS3 8o127 FOSAKUI FUSAKUZ FUSAK03 xAFOgD_ XFUS1 XFUSZ 18 12. .20_ b,89 3._ _6,5 ZU°813 O° O, lO, _0, _0° _0° 50. 60° 100. l,60Z 2,9Z 5,31 7,1_ iO, O. L._ 4.5 7.5 10.5 11.667 _9°6 33° 38.5 _°01 -I° 85A006 XAF3 1_6° °_0_ _2,10_8_,_U_SZO,_8 _o01 O° NACA l KEFL KN02 70. 80° 90° XAFK_O XAFKII ¥Kb lg,O /5,596811.37_619,1_0320°928 XFUSKIO XFUSKI6 NALPHA _. NALPHA NALPHA 000000000H1111111122222222223333_33336666646666555_55_5560686o06867177177777_ _Z3_5_?89_g3_56?89_3_5_89_365_89_ig3_5_89_6_7_9_3_5_8_3_5_?89_ OGLVE CVLINOEg BOU¥ WITH _5 UEGREE 3_EEP NACA bSAO04 MID-WiNG i01 M|NG 1ANU PANEL COkNE_ J ZNUICAIk PANEL POINI WING X COORG[NATE$ PANEL LEADING-EDGE ¥ [ L L POiNT$_ X 1 2 AND Y 2 6 Z 2 _RAILING-EDGE |NOICATE Y X Z PUINT$ Z 3 3 X 3 Y 6 L 6 6 16.68370 17.37257 1*66700 1.6b?00 0.00000 0.00000 17.[_500 A?.91700 2.9/000 2.97000 0,00000 0.00000 17.91700 18.71900 2.9?000 2.97000 0°00000 0.00000 0.00000 O.OOO00 10.2_163 19.15030 1.66700 1.66700 0.00000 0.00000 18.71900 19.52100 2.9?000 2.9/000 0.00000 0.00000 [9.52100 20.32300 2.97000 2.9?000 0,00000 0o00000 O.00000 Z0.03917 1.66700 O.0OO00 20,32300 2.97000 0.00000 21.12500 Z.97000 0.00000 1,6o700 1o66100 0.00000 0.00000 20.92803 21.81690 [.b6700 1.b6700 0.00000 0.00000 21.12500 21*92700 2.97000 2.91000 0.00000 0.00000 21.92700 LZ,?Z900 2.97000 2,9/000 0.00000 0.00000 21*dlbgO 22°10571 l.bb?O0 1.06700 0.00000 0.00000 22.70577 23.59663 l.bbTO0 lobblO0 0.00000 O.O000U 22.?2900 23.53100 2.97000 0.00000 23,53100 2.97000 0.OOOOO 10 L[ 23.576b3 17.11_00 l.b6?O0 Z.97000 OeUO000 0.00000 26.68350 17.91100 lebb?O0 2.97000 0.00000 0.00000 26.33300 19.91500 2.9_000 2.97000 0o00000 0.00000 26°33300 25.13500 2.97000 2.97000 0.000U0 0.00000 18 I_ 17.9110U 18.71900 2.97U00 _.97000 0.00000 0.00000 18.71g00 19._Z100 2.91000 2.97000 0.00000 0.00000 20.55700 _1.19900 5.3?000 5.3?000 0.00000 0.00000 20.55100 21.19900 5.37000 5.37000 0.00000 0.00000 16 [9.52100 2.91000 0.00000 20.32300 2.91000 0.00000 21.8_100 5.3?000 5.3?000 0.00000 0.00000 21,86100 22.48300 _o31000 5.31000 0.00000 0.00000 0.00000 0.00000 2_.|2500 2_.76100 5.37000 5.37000 0.00000 0.00000 1 Z 15.59683 lb.683?O 1.66700 1.66700 0.00000 0o00000 J ¢ _7.J?257 _8._b16_ l.b6?O0 1.bb?00 5 19.15030 1.66700 b ? 20.03911 20.9280J 8 9 [5 16 20.32300 21.18500 2.97000 2.97000 0.00000 0.00000 21.1_500 21.92700 Z.97000 2.9?000 0.00000 0.00000 22.68300 _3.12500 5.31000 5.37000 L? 2[.92700 Z.91000 0.00000 22.12900 2,97000 0.00000 23.76700 5.37000 0.00000 26°60900 5.37000 0.00000 18 L9 22.72900 23.53100 Z.97000 2.97000 0o00000 0.00000 25.53100 26.33300 2.97000 2.9?000 0.00000 0.00000 26.60900 _5.05[00 5.37000 5.3?000 0.00000 0.00000 25.05100 25,69300 5.37000 5.37000 0.00000 0.00000 ZO 21 _6.3_00 19.91500 2.97000 5.37000 0.00000 0.00000 25.13500 20.55700 2.9?000 5.37000 0.00000 0.00000 25*09300 22.66833 5.3?000 0.00000 20.33500 5.37000 0.00000 22 2_ Z0.§5100 21.19900 5.37000 5.3?000 0.00000 0.00000 21.19900 21.86[00 5.3?000 5.37000 0.00000 0.00000 2_.15300 23.63767 ?.73000 7.73000 0.00000 0.00000 23.[5300 23o63767 /.?3000 7,?$000 0o00000 0,00000 26 21.06100 5.37000 0.00000 22.68300 5.3?000 0.00000 26.12233 7.?3000 ?.?3000 0,000OO 0.00000 26°12233 26ob0700 ?.73000 7.?3000 0°00000 0.00000 25 2b 27 22.68_00 83.12500 23o10100 5.37000 5.3?000 5.3?000 0.00000 0.00000 0.00000 23.12500 23.70100 26.60900 5.3?000 5.37000 5.37000 0.00000 0.00000 0.00000 26.60?00 25.09167 25._1633 ?.73000 ?.?3000 0,00000 0.00000 25.09167 25.5?633 ?.?3000 ?.?3000 0,00000 0.00000 28 26°60900 5.3?000 0.00000 Z5.05100 5.37000 0.00000 26.06100 ?.73000 7.?3000 0.00000 0.00000 2b.06100 26.56567 7*?3000 7,?3000 0.00000 0.00000 ?./3000 0.00000 2?.03033 /.?3000 0.00000 7.73000 L0.10000 0.00000 0o00000 27.51500 25.76000 /.73000 10.10000 0.00000 0.00000 10.10000 10.10000 0.00000 0.00000 20.08667 Zb.61333 L0o10000 10.10000 0o00000 0o00000 [0.10000 10.10000 0.00000 0.00000 26°?6000 2?.0666? 10.10000 10.L0000 0o00000 0.00000 0.00000 0.00000 2?.39333 2?.72000 |0.10000 LO.LO000 0.00000 0.00000 29 _O 25.05100 25.09_00 5.37000 5.37000 0.00000 0.00000 25,69300 26.33500 5.3?000 5.37000 0.00000 0°00000 26.5656? 2?.03033 Jl 32 22.bb_33 23.L§_00 7.?_000 /.1J000 0.00000 0.00000 23.15300 23.63767 ?.13000 7.73000 0.00000 0.00000 25.63333 25.76000 33 36 23.63761 26.1223_ ?.73000 7.73000 0.00000 0.00000 26.12233 26.00100 ?.13000 7,73000 0.00000 0.00000 26.08667 Z6,61333 35 J6 26.60700 25,09L67 7.73000 7.13000 0.00000 0.00000 25.09167 25,57633 7.73000 7.73000 0.00000 0.00000 20.16000 27*06667 37 25.51033 7.13000 0.00000 26.00100 _.73000 0.00000 27.39333 [0.10000 LO.IO000 0.00000 28.06067 10.10000 0.00000 J8 2o.ObEOU /.1_000 0o00000 2b,56567 7.73000 0.00000 27,12000 10.10000 J9 60 26._567 27.03033 ?.?3000 7.73000 0.00000 0.00000 27.03033 27.51500 7.73000 ?.73000 0o00000 0.00000 28.0666? 28.37333 10.10000 0.00000 28.37333 [0oL0000 0.00000 _1 62 6J 25._3_3_ 25.76000 2b,OSbb? [0.10000 10.10000 10.100OO 0.00000 0.00000 0.00000 ZS.?bO00 Zb,OBbb7 26.61333 10.10000 L0.1 0000 10.10000 0.00000 0.00000 0.00000 2?.65000 2 ?.85000 28.05000 10.10000 [2.00000 12.00000 0.00000 0.00000 0,00000 28.70000 2?.85000 20.05000 I0.10000 12.00000 12.00000 0.00000 0.00000 0.00000 66 65 26.6[_33 26.76000 10.10000 10,10000 0.00000 0.00000 2b.7_000 Zl. Obbb? IO. L0000 10.10000 0.00000 0.00000 28.25000 28.65000 12,00000 12.00000 0.00000 0.00000 28.25000 28°65000 12.00000 [2.00000 0.00000 0.00000 60 27.06667 10,10000 0.00000 27.39333 10.|0000 0.00000 28.65000 12.00000 [2.00000 0.00000 0,00000 28.65000 28. B5000 12.00000 12.00000 0.00000 0.00000 _? 27.39333 IO. 10000 0.O0000 27.72000 [0.10000 0.0OOO0 28.85000 12.00000 0.00000 29.05000 12.00000 0.00000 68 69 27.72000 28,06661 10.10000 I0. I0000 0°00000 0.00000 28.O6667 28,3?333 [0.[0000 10.10000 0.00000 0.00000 29.05000 29.25000 12.00000 12.00000 0.00000 0.00000 29°25000 29°65000 12.00000 I2o00000 0.00000 0o00000 50 2_.J?_ 10. 0.00000 28.70000 10.10000 0.00000 29.65000 L2.00000 0.00000 29o65000 [2°00000 0,00000 102 tCO00 _ING PANEL CJNEKOL PUJN| PU|NT$ X AND ANbLES Z CP X tNCL|NA||ON y CP |HE[A CP CAHBEK $LUPE tHiCKNeSS SLOPE le,3_Lvu _.3073_ O,OUOOO O,O000U 0,00000 ,LTB?9 l?.L_dU8 _*_0734 O*OUO00 0.00000 0.00000 *05896 1_,0_ LO.BBU_j _.3073_ Z*301_ O,OOOOU U,OUOOU 0,00000 UoOOUOU OoOUO00 0.00000 ,0_265 *01709 5 19. _,3073_ OoOOOOU U, U,UOOUO o 7 _0,51_ 2i*_LB_b £,301_ Z,30734 0,00000 U*GOOUO 0,00000 U.O00UO O,UO00U O*O000U -°01379 -*02918 9 TL6_L ,00315 Z_.Zo_£_ 23, L£LJZ _,301J_ Z*307J_ U, U, 0,00000 0°00000 O, OOOO0 O*UO00O -,0J946 -,04623 10 23.95749 2o30134 0*00000 OoUOOOO 0*00000 -o04706 LL LZ _4.a03o7 L_,4o_9 _,JUl_ _.1Z5_ O,OOOOO U,OUOOO U°OOOOu U°UO000 U.O0000 0°00000 -,0_710 ,17879 i_ L4 Lg.L_8_b 19.9£_10 4oIZ_bB _.£_508 U,OOUOU 0.00000 0*00000 0°00000 0,00000 U.UO000 *05696 °UJ2b5 £5 ZO.C_L6 4. U.OOUOO 0.00000 0.00000 16 L7 ZL._o_ii 2Z,Od_U? 4°IZSb_ _°1_5o_ 0.00000 0,00000 U. O0000 U.UOOO0 0.00000 0,00000 .O0_l§ -*01379 IB L9 _Z,SL_OZ 23,5_198 4,1Z568 _o1_568 0,00000 0°00000 0,00000 U, OUUOO 0,00000 O,OOUOO -,0£918 -*03g_.6 ZO _i _._b_93 24.9e7a8 _°LZS8_ 4. i_Sba 0.00000 U,_OUO0 O.O0000 O,O00UO 0.00000 O, OOO00 -.0462_ -*0_70_. _z _5.IIz_ 6. O,UUOUO O,O000U 0,00000 -,04710 _ Ze ZL.ZZ759 21,19_58 b*49507 o,4950_ 0*00000 O*OOOOO U*UO000 O*UOOOO 0.00000 O*OOOOO ,L1879 ,05696 _ Z_°_6L_B 0°49507 O°OOO00 0*00000 O.OOO00 *O_b5 Z6 _7 Z_.9Z8_7 _3.49_57 b,69507 6o_9_07 u,O000U U.OOOOO O,O00UO O°UO00U 0,00000 U.00000 ZU Z9 Z4,UbZ_6 2_._Z9>6 0.49507 Oo_9_0_ O,OOOOU 0.00000 0,00000 O.OOOO0 0.00000 O.O000O -.01379 -.0Z918 _0 _i _5,£9o_5 25.76355 O°_9_O_ 6._9507 0.00000 O,OUO00 U. O000U O,UOUOO OeO0000 O, O00UO -.01966 -,0_42_ _ 2o.3_05_ 6.4950_ 0.00000 0.00000 0.00000 -.O_?O_ 33 3_ 26,_91>4 Z3.98_09 0*_9601 a*83_UU 0o00000" U.OOOOO 0°00000 0°00000 O,O00UO O°UO00U -,04*710 *L_819 _5 J_ Z¢,31188 z6.18Zo8 o,o_aoo _.83_0_ 0.00000 U,00000 Oo00000 o.O0000 0°00000 0.00000 _Z _ _5°£93_1 Z5,_0_£? 8*83808 8. B_UO_ 0°00000 O°OOOO0 0,00000 O,UOOOU 0°00000 O,O0000 39 20*01506 U*d3bO_ O*O0000 0*00000 0°00000 -*01379 -,02918 -*039_6 LZSb_ iZ568 UUUUU O00UU UOOO0 .0L?09 .01709 oUU3L5 .0569b .03265 ,Ol/O9 ,003E5 _U 41 26.4_8o 26°B_605 8,03_0_ 8°e3808 O,UO000 0.00000 0,00000 U*O0000 0.00000 0°00000 4_ _7. 8._800 O.OOO00 0.00000 0.00000 -*O_b2J 4J 4_ ZT.O58Z4 Z_°ObgO_ 8, 8. o,0000o O. O0000 0,0000o 0.00000 0,00oo0 0.00000 -.0_704 -.U_?LU _5 _6 2_._5Z81 Zb.?Zt2Z LU.91384 i0,973_ 0o00000 O. G0000 U.O0000 0.00000 0.00000 0.00000 .L7879 *05696 _ 4_ 26.9_9o_ Z_,_SaO5 L0.9138_ 10.9738_ 0°00000 0,00000 O.O0000 0,00000 U.O0000 0,00000 .0_265 ,01709 _9 50 21*5Z6_O _7°79_7 10.g1_84 10.9_8_ 0.00000 OoO0000 0*00000 O,0OOOU 0.00000 U,O0000 .00315 -.OLd/9 _1 28°00_28 10°97384 U*OO00U 0.00000 O*OOO00 -°02918 52 5_ ZB.33169 28°60010 10,91384 _0, g7384 0°00000 O*OOO00 U,OOUO0 0*00000 OoUO000 0.00000 -.039'bb -*O4bZ_ 5_ 55 ZB*Oo_SL Z°I_O_ 10°913_ L0*91_84 O, OOOOO O,OOUOO 0°00000 0,00000 0,00000 0,00000 -.0_706 -*06710 Z_7_5 a3_oa B_08 103 NING PANEL PANEL AKEA A_A$ AND _NUAD 1 1.101o0 .8_b18 Z loiOlbO i.i01oO loLOlbO .8_b18 ._4blB .8_618 5 b lolOloO loiOlbO °8_b18 o84bI_ Z I.LOIoU o84b18 8 9 l.lOlO0 1.lOlbO .8_b18 o84618 10 i./O/oO ._4bl_ II 12 iJ L.I_Z_O 1.7_ZeO 1°7_280 .1Z_95 .7_95 .72_95 L_ L.TJ_80 .72_95 15 i6 L./3ZSO 1.13ZBO °72_95 .12_95 L1 L.13Z_O °7Z_95 18 19 Lo73Z80 1.1J2UO .72_95 .72_95 gO 21 1.13Z80 1.32947 .12495 °50700 _Z Z3 l,JZg_7 i. dZg_7 °56700 °50700 _ 1._Z9_7 .56700 22 _6 1.3_9_1 1,._9_7 oSblOO .SbTOQ Z7 l. .5o700 Z8 Z9 30 [.329_1 L._Z947 L.._9_I Bdg_l .SblO0 °50700 .5_700 31 .9b[_3 .4|079 32 3_ 34 .9b143 .96L_3 .9bi,3 °4|019 .41079 .41019 35 .9bI_3 .41079 3b 31 .961_3 .9b1_3 .41019 °41019 38 _9 .961_3 .9bl_ ._1079 ._I019 _0 .9b143 .41079 _1 _Z . 50033 .5003_ .Zb_l ._b841 _3 _ °50033 .50033 °26841 ._b841 _§ _b .500J3 .500_3 ._b84l .268_l _I .50033 .2684L _ _9 °50033 .50033 .2684L .268_1 _0 .50033 .2o841 104 _H_KO$ BOUY I &NJ PANEL CUNN_ _ LNu|C_IE _A_EL PU[NT _UU_UJNATh_ OUOY PANEL LEAD|N_-EOGE x y L l Z L l Z POINTer x AND y Z _ |NO|CATE Z 2 2 UoOUOuo O.OOUUU LoSO00O °00000 -o_06Z2 UoUOUO0 O.OOOOO 0°00000 OoOOOOO UoOU00O 0°00000 OoOOOO0 Lo_OOOO lo_OOOO o£_Z_ o4062Z -.Zd72_ .00000 OoOOOO0 O.OOOO0 5 0o0000o L._OOUO O.UGUO0 .00000 OoO000O -o_06_ ko_UOOO _._0000 .2_174 .00000 o Z 1.50000 L,SOOUO o£O1_4 ,_OZ£ -oZ8/Z4 .OOOCO _.SOOOO 4.50000 .7_Ud3 l,O_B/ 9 Lo)O000 _.o000o .Z_IZ_ .O000o .LUTZ_ -1.0_487 6.50000 7.50000 .738_3 .O000D LU _._UO00 .7_d83 -._dU3 7.50000 &,O_O_7 II LZ 4,_OOOO 4,_0000 _.0_8/ .7_803 1.50000 7.50000 L,_51_0 1.03U_7 L_ 7._0000 i_ i_ Z.5000O 7o50O0o .00000 ,l_8_J OoUO00O l._O00O -./388_ .00000 i.50000 1.50000 .TJBd_ -1._57J0 l._OOO0 6.50000 -_.0J047 _.50000 oOOOOO 1.0_0_7 4.50000 _,50000 .Z_724 -°00000 -.2812_ .00000 _,50000 _.5000U .73883 1,0_87 -.73803 .00000 6.50000 ¢.50000 .73883 -.OuO00 .7_B83 1.044_7 7.50000 i.03047 1.50000 /.50000 1.4_/30 Io03047 • 2_724 • 40_£Z • 2877¢ -.00000 °7388_ 1.04_1 • 7_8u_ -L.650_0 1.50000 -.00000 t.O_O¢! -[.£bbg_ .00000 1o50000 7.50000 X.45730 I.L669_ -1.60670 -L.loo_4 .00000 LL.66700 LL._6700 L.17854 _._6670 [.L6O94 LL.66700 L*1785_ -.00000 L.[78_4 1.05030 -L.LTd54 IZ.66700 L5.59¢80 -.00000 1.1785¢ l.o6b/O 1.1/854 .00000 I.LTd54 15.59_80 L_.594_0 L.6_6/0 1.17854 -.00000 1.6bbTO 15.59480 -.00000 1/.37260 17.372o0 L.17854 L.Ob_TU LT. 372bO 17._726U 1.L/856 -.00000 -1.L785¢ .00000 19.150_0 19.150]0 L.L7854 1.66670 1.1/_54 L9.1_030 L.LT_54 I.ob_70 -I.17d54 Lg,LS030 Z0.92800 -.00000 1.1/854 /.1oo_ -i.lo69_ 11.6o700 I,I785_ IO.50OOo lO.50uO0 L.obO_O 1._b09¢ t_.o6700 //.OoTOO L._oo?O L.17854 2i _Z Li.60100 LL.ooIOU .OOOOo i. Lldb4 15._4_0 15.59480 o00000 1o1185_ Z_ Z¢ i_oOO70U LI.oo70U /.6bo?O 1.17_4 15o59480 1_._94_0 I.b6o/O 1.17_5_ _5 i5.59_00 -_.o_70 I7._7ZO0 -i.bb610 I_._9_80 Zb Z1 L_.)9_80 [5.59_d0 I./7_b¢ l°O_/O -1.I785_ *O0000 LI*_7Z00 I7.a7ZoO L.1105¢ L.bOOlO -I./785_ .OOOO0 I5.59_80 LS.)9¢OU 18 _9 15._9_0 17._7_b0 /.17854 .OOOOO I.L?_5_ -t.bbblO I7.3]760 19.LbO_O I.I7d54 .OOOOO 1.1705_ -L.O_blO 15._80 I7._1160 _0 17*_7Zo0 I.I7db4 -i. 1_.I50_0 _,I1_54 -L.11054 I7.372b0 L.1785_ 1.6b070 Jl JZ I7.$77b0 L7o_77O0 L.OC_70 I.11_5¢ L9.15030 I9.,5U$0 1.b6670 1.1785¢ L1._7260 17.37200 1.L785¢ J3 J4 I9./505U 19.L5030 .OOOOO i. L7854 20.92800 ZOoNZOUO .OOOOO L.1785¢ -.00000 1.L7854 II_5_ -I.00670 -L.11854 .OOO00 1.17854 -L°Ob670 -L.L185_ I1.0o700 I/.Ob700 i9./5030 Lg. LS030 _ L9°/5030 i.bb_70 20.9_U0 LoOooIO zg.l_O_O 20._OOO 1°118_4 °OOOOO l,Lld54 -1.oo070 /U. gZO00 22./0580 [.I785¢ .OOoOO L.LI85_ -L.66o70 19oL5030 20°92800 _B a_ ZOo92dUU 20.97_U0 L.II_ /.b_o70 -1o11854 °00000 72.70)B0 22°105_0 I./7054 t.bbo70 -L.I185_ .00000 Z0.92_00 20.92800 _O ZU.gZoO0 i.iloo4 2L.70500 L.17854 1.L1854 £OogZSOO _I _ ZZ./05_O 22,105BU .OOOOO 1.11_4 2¢,_8J50 24.4_50 .O00oO i.1785_ _ 4_ Z2,IU_O 22.105_0 I.b60/O i. L7_5_ Z¢._B_SO Z4._50 l,boblO 1.L7_5_ _5 2¢,_U -/.hollo 20.2_000 .00000 -L.ObblO Z_.4_J50 • 0 _1 _8 Z4.485_0 Z¢°_50 2_._d050 1.17_ l.b6_/O 1.17d5_ -/.I7_5_ .O00OO 1.118_4 Z0.28000 Zb./_OOO 2o.28000 I./7854 £,06070 L.i7854 -l.i185% .OOO00 1.17856 24.68350 2¢°4_50 Z_o¢d350 • _ 50 Zo.ZoOOO 2o. ZdO00 .O000O i. I7854 -*.oO0?O -Lo_185_ Z_o¢0000 zg._oOOO .OOOOO L.I7_5_ -I.b6070 -I.17854 26.78000 26.28000 5I IO.Z_O0_ i.ooo/O _2 _ 26.2_000 29._0000 i. 5_ 55 zg._uOOO 29._0000 L.I1_54 L.OoO70 50 79.40000 L.II_54 57 _8 _°OOOOO 33.00U00 .OOOO0 L.L/85_ _9 O0 _°OOOO0 _3.OOOOO I.bb07O I./7_5_ .00000 L185¢ .OOOOO L.LId54 -l.6b6lO -L.L7_5_ .00000 I.ildb_ .O0000 .O00OO /1.0o700 IL.bO700 Jb _7 .00000 -l.ObblO -i.17U54 .00000 L.LI_54 1_o1_030 ZZ. 10580 72.70580 2Z,70580 27.70580 79.40U00 1°_6_70 1o17656 -1.06670 79.40000 33.00000 L.Lld54 °00000 /.17854 -l.bbb70 2_.2_00U Z_.40000 -I.17854 °OOOOO _a. OOO00 _.00000 L./7854 1.b6070 -1o17854 .00000 79._0000 29.60000 3_.00000 1.1/854 _b._O00O Sb.50000 *OOOO0 1.17_54 _oo50000 36.50000 1.66bi0 _.lldS_ 1.1185_ -l.bOOlO -£.I1_54 °00000 L.;185_ .UOOO0 L.LI85_ -l,bbblO -i.iT_5¢ o00000 i. I185_ LOO5030 L.lb094 I0o_0000 I_ _O oOO00O l,ll_ -o00000 1olb694 LO.50000 L.Lbb94 L.o50_O i8 .OOO00 ?._0000 10.50000 10,50000 I. .00000 o00000 /0.50000 10,_0000 10.50000 IloOOlO0 .00000 I.LlSb¢ 1o0_¢8/ -1.030_7 L.lob94 -.00000 1.030_7 -L._5030 .00000 i.il_5_ .00000 °/3885 IU.50000 LO.50000 1°0_047 .00000 -L.bO610 -L.I7854 -o73_ i.0_0_7 1.45730 7.)0000 IO.)ooOo -I°_bolO -i.i185_ .1U72_ .¢OOZZ .Z0124 ._0622 L.O_O_7 -.OUO00 Lo LI .00000 1.I7854 ZO°ZdOOU 29.40000 33.00000 _3.00000 33.00000 _J.OOO00 -o287_4 .00000 Xo50000 1.50000 1.50000 1.50000 IoIOO_4 1.6b030 -L,IIB54 4 OoOOOOO O.OOO00 IO.)O000 LO.SOO00 .OOO00 A.I66_ 4 O.OOOO0 0.00000 L0.50000 7._0000 /0.50000 Z y 4 0.00000 0.00000 -I.OJO_l .O00GO i6694 oOOOOO X 3 0.00000 0o00000 -L.457_0 .00000 l 3 1.03067 _._7_0 .O00OO POINTS Y 3 U.UUUOO .dBIZ4 -L.0_87 TKA_LIN_-ED_E X L.L;854 L.b6o/O IoL/BS_ -o00000 -l.L/d54 .O000U 1oL7854 |.bbo#O 20.92800 1.b0670 L.L785¢ -.00000 1.1/85¢ 1.66670 20.92U00 Z0.92800 1.17_5_ -.00000 l.LT_b¢ 1.6oo/0 -L.L7d5_ .00000 22.70580 ZZ./0580 L.L7854 L.h66/O 22.70580 _Z.;0580 t.1/854 -.00000 26.68350 L.L785_ Z_.48350 26,68350 1.06670 L.L7854 74.48_50 ZO. Z8OUO -.00000 1.t/854 Zb.28000 Z6.28000 L.b6blO l. L785¢ 2o.28000 Z9.40000 -.00000 1.17854 29.60000 1o066/0 L.o66/0 L.L7854 -.00000 L.L7_54 L.oboTO L.17854 -o00000 LoL78)_ L.6bo/O L.L;dS¢ -.00000 1.17_54 Lo60070 .00000 L.1785_ 1.660/0 -l.17ub4 .UO000 t,I7854 _o06b_O -1o17854 .00000 l.lldS_ 1.66o70 -1.17854 o00000 1.L785_ -.00000 1.[785¢ 1.66o70 29.60000 29._0000 1.L7854 -.00000 1.17_5_ L.660/0 -1.L7_5_ .00000 3J.O0000 3_.00000 1.L/854 L.66670 Lo1785_ -,00000 1o17_5_ [o66_70 _.00000 _3,00000 L.1185_ -,00000 L.L7854 -1.L7854 ,Z872_ .406_2 -I°03047 .00000 1.03047 1._7_0 -Lo1669_ .00000 1.16694 1.65030 -1oL1_54 .00000 X*LI_5¢ 1.66b/d -I.L785_ .00000 1.17854 1.6667U -L._78_4 .00000 I.L785_ l.bOO;U -1.1785¢ .00000 1.1785¢ I.bbbTO -1.I785¢ .OUO00 1.17854 1.06670 -1.L785_ .00000 1.L785_ I.bbo]O -L. LT_5_ .00000 L.1_854 _.b66/0 -i. LTdS_ o00000 &.L785_ 1.6bb70 -l.l/dSe o00000 1.17d5_ L.66070 -1.1785_ .00000 L,lld§_ Lob6blO 3o.50000 1.L7854 t.66670 1.17854 .00000 1.17854 3o.50000 _b.50000 _.666/0 1.17854 -I. .00000 1.17_5_ L785¢ -.UOUO0 1.66_70 _o.50000 -.00000 1.66670 105 00UY PANEL PU|N_ CONTKUL X CUORDINATES Y GP 1 2 PUINT Z LP CP 1.00000 t.OOO00 °095?5 .23L16 I.OOUO0 ._3[[o 4 5 o L.O0000 3.22006 J.dSOOb .09575 ,_1308 °65928 7 3.880Ub .0_928 .21308 8 9 J,ZZOOb 6.08_e2 .21_08 .*_b33 .65928 -[°07754 j -.2JILb -.09575 .09575 ,23Lib -,b5928 -.27308 10 0.08_48 1.07754 -,64633 11 12 6°08242 0.08248 _,G7754 .4_b33 .446_3 1.01754 [j 14 9°03[05 9.0)105 .55006 1.3219o [_ 9,03105 1.32796 [b 17 9.03[05 IE. 08446 18 19 11.0844b 1[.084_b 20 21 Ll. OU44b 1Jo63090 22 13.630U0 J,68262 83 24 1).63090 13.63090 k.48262 .58927 25 20 L6,48_70 L6.48310 =5892? 1.422b2 27 28 16.48_70 10.48_70 L._2_68 .58927 29 30 L8°26145 L8.26145 .58927 1._2268 JL 38 L8,26L_5 L8°261¢5 L._2262 .58921 3_ 20,03915 34 20.03915 1,42262 35 30 _7 20.039[5 20.039L_ 21.81690 1.62268 .5B921 .58927 38 39 2L.8_090 21.8169U 1.42262 1,_2268 40 21°81690 .58927 4[ 48 _3.594o5 23.59_65 .58927 1.62262 4_ _6 23.59465 23,594o5 1.42Z62 ,_8927 45 46 25.38115 25,3811_ .58927 k.62262 _1 25.38115 1.42268 48 49 25.38175 21.8_000 ,58921 .58927 1.42262 -[.62262 50 27°84000 1.42262 -,58927 51 52 27.84000 21.84000 [.48868 .58927 _3 54 31.20000 31.20000 ,58927 1.62_b2 55 50 31.20000 3L.20000 L,62262 ,5892? 57 3_,15000 .58927 58 59 60 34.75000 3_o15000 34,75000 1.42262 1.62_b8 °5898? 106 .550U6 ._8637 1,61563 1.41503 ,5Ub31 .58921 .58927 -1.32796 -.55006 .5_006 1.3279b -1.41563 --.58037 ,58637 1,41563 -1.62_02 -.58_21 ,58_21 1,42262 -1o42262 -,58927 °5892? 1.42262 -I°62262 -,58921 ,58927 1,q2_62 -1._2202 -,58927 ,58927 1o42_62 -1°48_02 -,58927 .5U927 1.62262 -1.62262 -.58921 .58_27 1.42202 -[._2262 -,58921 .58927 .58921 1.42262 -1._2262 -.58927 ,58927 1,42262 -1.62262 -.58927 ._8921 1o62262 BODY PANEL AREA:> PANEL AREA ANU INCL|NATIUN D='-L IA ANbLES THETA i .2_031 .26037 .24517 .zSbl? -2.76889 -1.96350 J 6 .24037 ._60_1 .24511 .265IT -l.17_lO -.39270 5 t.691_ ,19_20 -2.14d89 6 ? L.697U4 1,09786 .29_0 ,I962U -L,96350 -L.t?8tO B 9 1.597d6 _,89570 ,19680 ,IZ536 -.39270 -2°74689 IO L1 2°89570 2.8957U .I2036 .I_03_ -1.96350 -_.ll_LO EZ i3 _,895?U 3._139_ .L_03_ ,05931 -,3W210 -2.16889 I4 3°57398 .G5931 -&.96350 X5 3.5739_ .U_93T -t. 1o I7 3.5739_ L._8161 .05937 ._129d -._9270 -2.14889 18 19 Io68L_1 1,68*67 .g1298 .0tZ98 -L.96350 -£oL78i0 20 ZL 1.4816T 5.0i045 .G1_98 O. O00OO -t39ZTO -2.74889 22 5.0104_ U._UUO0 -1.90350 _3 5.010_5 U. -L.l?StO 24 Z5 _6 5.0/U_5 _.2b78J 2.26783 O. GO000 O._UUUU OoCOOO0 -.J9270 -2.76889 -1.96350 -Z,ITBiO COOUO 27 2.25783 O. 28 29 2,Z5783 2.26770 O. O0000 O.GOO00 -.39210 -2.76889 30 2.goTlU U.UUUUU -i.9_350 3i 32 Z.2bTTU 2.20710 O. OOOOO 0.00000 -I.I/SlO -.39270 _ 34 Z.Z_?_O _.2bTIo U,UOOOU O.O000U -2.14889 -1.90350 35 36 2.207?0 2.26710 O. O000O 0.00000 -t./78LO -.392T0 37 _.26183 O. 38 _9 Z.26783 2.2o783 O._OOOO O.O000O -/.9635U -l.l?810 40 _l Z._6783 2.86770 0.00000 0.00000 -.3921U -2.14889 4Z 43 2.261?0 2.26770 O.OOOOO o. oOOOO -1.9o350 -_.l?8tO 44 2,2071U 0.00000 65 4o 2.29168 2,29168 0,00000 OoO0000 -2. T6889 -_o9635U 47 2.2916d O. -l.llBiO 68 _9 50 2.29108 _.91999 3.97999 0.00000 O. OOO00 O. OOOOO -°39270 -2.14889 -1.96350 51 3.97999 O.OOOO0 -&,£7810 52 b_ 3.91999 _,59229 O.OOOOO u. UOUOU -.39210 -2.74889 5_ 4. O. -1.96350 55 55 _1 4.59229 4.592Z9 4.66473 O. COOOO O.OOOOO O. GO000 -I._?810 -._927U -2.7_889 58 59 4.6661_ 4.66_7_ U.OOOO0 0,00000 -i.9635U -1.11810 60 4._6613 0,00000 b92_9 OOOOD L1BLU COOO0 OOOOO GOUUU PAgTITION INFLUENCE = UF i TIME 80UY ON = 69o82100 80OY PARTITIUN INFLUENCE = OF 2 TIME WING ON • 61.21900 80bY PARTITION INFLUENCE = UF 3 TIME _OUY ON = 72.2itUO _ING PAWTIIION INFLUENCE • OF 6 I|NE WING O_ = 82o8510U WING -Z.76889 -.39270 -.39270 94.18100 T_ME = rINE = 95.90300 TIME = LO3,35?UO i07 VELuCIT[E$ ON 60_Yt MACH=2,OIO ALPHA= PANEL Nu. 3uUR_E bTgENG[ 1 2 .1964_ ,19662 -.09678 -.0907d ,0864B .20879 .1904Z o190_Z -.09_78 -. 09078 ,20879 .08668 5 .i126i -.0701b .00998 -.ibBb ,19668 b 7 o17_61 .I?Zbi -.070i6 -°07016 .L8896 .1_896 -°06998 .08998 .19668 .Lb68 8 9 ,17_61 olLL9o -°07016 -.0_539 .06998 ,06689 o/bBb -.11319 .19668 ,I2701 iO 1k .I129b oi1296 -.035_9 -.0_539 .I1319 .I1319 -,06689 .06689 .12101 ,12101 12 ld °it290 oOSL9B -.03_39 -.00181 .046d9 °02270 oll_19 -°05481 ,12701 ,05966 1_ iS °0_198 .U_198 -°00181 -.U0181 ,05481 .06681 -,02270 .02270 .05944 .05966 ib .05198 -.0018l .02270 H AXIAL VELOC|[V O.UO0 LATLRAL VELOC| IY V_k[iCAL VELUC| r¥ NORHAL VELOCITY -,Z0879 -.0864_ .25020 °25020 .08b_8 .Z0879 ,25020 ,25020 .05481 -,o_oo -.o2Lu6 .020_0 .02050 .00501 .01226 19 _0 -.0L28O -.022_b .02050 .02050 .012Z4 .00507 .00507 .01224 21 2L -.OOLll -.00217 .01801 .UIBO/ .00000 .00000 .00000 .00000 -.00000 --.00000 2_ -.00277 .01807 -.00000 .00000 -.00000 24 2_ -.00271 .0_891 .01801 .00714 -.00000 .O00UO .0U000 .00000 -.00000 -.00000 2b 21 .02u91 .0_91 .U0714 .00714 -,00000 .00000 -.00000 -.00000 .00000 -.00000 28 29 .02891 -.00208 *00714 .00751 -.00000 -.00000 .00000 -.00000 -.00000 -.00000 30 _1 .02_71 .0_31L -.00807 -.00807 -.019_9 .019_9 -.00000 -.00000 ,00151 .00037 .00000 -.01960 -.00000 -.00800 .00055 -.00635 -.015_6 .00055 -.00635 32 33 -.UULO_ -.02S2_ 34 .01089 3_ .01089 3b 37 -.02323 -.040L_ 38 39 -.01999 -,01999 60 41 -.06016 .00854 42 63 -.006_? -.00627 64 45 .008_ .0_622 4b 47 -.013_6 -.0t326 -.00814 -.00814 .0003_ -.0010| .01634 .01434 -.00101 -. 00006 .0191B .019L_ -.00006 .00587 -.0196b -.01001 .00025 .00025 -.01001 .00246 -.012_4 -.00507 .05946 iI 18 .0153_ .0129B .01298 .0129B .01298 ,00000 .00000 -.00000 -.00000 .00806 -.00415 .00000 .00000 .00060 -.00060 .00000 .00000 ,OOZE5 .00101 .00711 .00711 .Olllb -.OiTlb °00246 .01490 -.00101 .006L7 .00000 -.00000 .00000 .00000 .00000 -.00000 .01811 .01811 .0122_ .0122_ .02952 -.02952 68 .0_b22 .00587 .01690 --.00611 49 50 .01696 --.ObOb6 .01650 .O01II .020_b .00965 5[ -.ObO00 .00717 .O09b_ -.02_9 °0_450 .01317 *0203b .02397 -,00843 .00993 --.00000 .00000 52 5_ ,01694 .0_094 .00863 .02329 .00000 .00000 -.00000 -.00000 .00000 .00000 _4 -.0_9_ .011_ ,0_1_8 .02698 .00000 55 5b 51 -.0_393 .03096 .0/268 .01141 .01317 .00842 .01118 .0z397 .02074 -.0209B -.009_ .00859 .00000 .00000 .00000 58 -.01801 .00925 .012_0 .02993 .00000 59 60 -.01_01 .01206 .00925 .00842 .01260 .02074 i08 -.0299_ -.00859 .00000 .00000 d_I¥c CYLINULK S[nWGULAKIIY INIE_kATION UN IHc MACH= ALPHA WITH UF IHE FUR P_LSSUKE _,50_-GgEE SWE_-P SaMPLh NACA 65A00_, MIU-WING CASE UISTK|6UI|UN ouUY = X Pu_NI 6GOY PAN_L|NG 2.0100 0.0Q00 Y L XlC 2Y/d LIC _P CN CT CM POINI 1 k,O0000 ,095?5 -.23110 .0_1@0 .0_875 -.069_2 .L5199 ,QJZ7_ .00887 .64669 I d l.O000U L,O0000 .gJ116 ._3110 -.0g_75 .09575 .0Z140 .02760 .069q2 °06962 -,UZb?_ .0Z675 .L_L_ .l_gg .0135_ -.OLJSO .00687 .00867 .26101 -.267B7 2 3 .02760 o0Z675 ,06942 .15199 -.03216 .00887 -.66609 .08622 °08822 .0_201 .19798 .03699 .03699 J.O32ZI L.Z5598 .0_822 .088_ oL9l_8 .08201 l.OOOOo .0957_ ) o 3.2Z000 _.220_6 .27_06 .65928 l 6 J._2006 J.2ZO00 ,o_92_ .Z1_06 O,d82_Z ,6_t33 .Z7306 .655Z8 -l.ullS_ .0_101 2.Z26o5 .057_9 .05759 ,0o_30 -,06330 .0_101 .O_iOl .9Z31S -.9_316 I.dI15_ -1.3_796 ,16664 ,2_74_ .13603 .I_518 ,32359 -._9079 .05/59 .00010 -.L5283 ,0003_ .0_101 .00002 -._5006 .Z6763 .39879 -.16_18 .00010 .00014 .55006 1.3279o ,_763 ,24743 °39879 .16516 .30368 ._0_68 .17609 .*2_II ,JU_68 ._0368 .626L1 .17609 .J7365 .L7696 -.427Z1 -.03515 -o162?0 0.00000 -Lo1685_ .3736_ .627_1 -.17o96 -.03515 --.06_39 0.00000 -.6840Z ._?345 ._7365 .65161 .6Z?ZI .11096 .1769o .17696 .6_1_1 -.62721 -,U_15 -.0351_ -.01613 .00139 .1o210 --,0_960 0.00000 0.00000 O.OOOdU .4860_ 1.1685_ --.12817 -.589i/ .589Z7 ._5161 .65101 .621ZI ._Z121 -.i1696 .17_96 -.0i61_ -.0141_ -.01226 .O/ZZ6 0.00000 0.00000 -.05309 .05309 _.4Z262 .65Lbl ./7090 -°0161_ .OZ_oU .50031 ,50031 .17696 ._27_1 -.0i_86 .0ib02 --.0_113 °01390 0.00000 0.00000 -.079_2 .035_7 ._UU31 .50031 ._Z?ZI .11696 .0160_ -.01680 -.01390 .03Li2 0.00000 0.00000 -.035#1 .01942 32 .5_902 °11696 -.6Zing -.00lL8 -.00246 0.00000 -.00192 -.56921 .5690_ .62_Z1 -,11690 -.00137 -,00119 0.00000 -.0009_ 33 3_ .58921 1.422o2 -l.6ggoZ .5690Z .56902 ,69172 .6Z?_l .l?bb .11696 .17696 .42721 -.6272L -°OOi_7 -.001/6 .00191 0.00000 0.00000 0.00000 .00092 .00192 -.00_01 -.58921 .5_997 .5977Z ._11Z .421_/ .697Zl -.11696 .17696 -,0Z805 -.02805 0.00000 0.00000 .02443 -._Z663 38 39 l.62Zbg ,5917Z .11696 .00401 60 --1.6_26_ -.56927 .06663 .0_643 ./709_ ._27Z1 .66063 o6_?Z[ ,11o96 _.42_6Z -1.62262 .64643 °69539 .L_696 .L7696 ._Z?21 -.627z1 -.58921 .5_9_7 .09539 .b9539 .68721 .627Z1 -.17696 .l?09b .07539 .76Z76 .|7696 .11696 -,58927 .16976 .e2?Zl .589Z? 1.6ZZoZ .?bZ?e .16214 .e_121 .11696 .17696 ._Z721 -.0i596 -.OgS_l .65679 .65679 o17696 ._21_1 -.42121 -.1769o -.0Z665 -.02_2_ .0_67_ .85619 .6Z121 .1169o ,i?6gb ._27Zl .95205 .95205 .1/696 ,_2721 .95205 .95205 .6g/2L .1/696 9,03105 9.03105 1.32790 .5_06 i? 18 1[,06666 11.08_6 .58637 I.61_6J -1._15o3 -.68_37 i9 20 11.08_6 II.0_6,6 /._lbb3 .5B631 ._631 1.615o3 21 1_.6_090 _Z 13o63090 i,_g6_ -._8_2/ i3 Z_ Z5 13.6_090 13.030_U 1o.66310 l.qZ_bZ ._6_Z? ._69d7 .589_7 1._226Z -1.6226z _6 il 10.46370 16._6370 i.6_ZbZ 1._2_oZ 26 16._831U .569/7 Z9 30 i_.2oi_) 18.20145 .589_7 1._2_o2 31 3g iS.gbl_5 i6.261s5 l.6gg6_ .58927 _3 Z0.0_91> 36 20.0_9t) l._Zzo_ 3_ _o 37 _0,03915 _U.0391_ 21._1690 1.4226z .5_927 .56927 3_ 39 Zl.81690 Z1,61690 i.6ZZbZ 1.6_/oZ _0 ZL.BIbO .569Zl 61 _ 23,59_o5 z_.59_05 .569Z7 1.4_zb_ _3 Z_,Sg40_ /.4Zgb2 66 _5 _3.59_6_ Z_.36175 60 _1 Z},3_115 25.38L75 _6 69 gS._i15 27.86000 50 _7.66000 1.6g_6_ 5t 5_ 27,86000 17.8_000 l._Z6Z .56927 53 5_ 31.g0000 31.2000o .58921 1._2262 55 _6 _1._0000 _i.g0000 1.6ggb_ .58927 57 58 36.15000 36._5000 .589Z/ 1.42Z62 59 6d 3_.75000 3_.75000 1._b_ .5892? ON I_.UC00 20.8130 CNCT- .OUUU .003? CL= -I._226_ -.569_1 .>6921 1.4L_62 -1._226g ._8921 1._20_ -1.42262 -1.62262 -.58927 ._Z? 1._2262 -I.4ZZO2 -.58_Z/ .5892? 1.62262 .1_609 .S2511 ._27Z/ -._2_ZI -olTbb .i1096 .62_gl .62721 .00010 .00010 -.0_988 -.0_966 -.03986 -.03988 -.05_58 -.02261 .OZZOI .05458 .00119 .002_6 .00600 -.02636 .0_3_ .00002 .00002 -.00071 -.00011 -.00077 -.00071 0.00000 12 13 ,00163 -.00163 -.00372 15 16 -.52986 -.21968 17 .21948 .52986 .1_817 18 19 20 2Z 2_ 2_ 25 Z6 21 28 Z9 3O 35 36 31 -.00400 0.00000 .00022 -.0_259 0.00000 O. O0000 -.00061 .09066 .03Z59 0.00000 -.09066 *000_0 -.01[UU -.OOOZ_ -.0Z515 0.00000 0.00000 -.03616 -.03o/0 -.0_111 .03171 0.00000 0.00000 .14690 -.14F_.90 66 .6Z1ZZ -.4Z/Z1 -.0118_ -.02881 .0Z51_ -.10593 0.00000 0.00000 -.1169L .16638 _8 -,17096 -.01596 -.0263t 0.00000 -._Z121 -.17696 -.621Z1 -.17696 o17696 .6_721 .00191 -.0001_ -.0003_ -2.22865 .00J92 .O00kO -,0J750 -.0375o -,O_Z_ -.020_5 -.01112 -.0192_ -.UIZ6 -.017L_ ,02_31 .L0593 -.1122_ -.0_08_ .0408_ .11222 -.070bl -.03288 .0_Z86 .07061 .00061 ,11691 .17085 0.00000 0.00000 -.|7085 -.74438 0.00000 0.00000 1.16565 .62408 0.00000 0.00000 --.62608 -l.16565 0.00000 0_00000 0.00000 0.00000 .98_06 .6_826 -.45826 -.98606 _2 _5 50 51 52 53 5_ 55 56 5? 58 59 60 BUuY R_FX= CM= -l._ZzbZ .1o518 .39_/9 -.6_511 -.17o09 .00002 10 11 _ucF_I_I_NIS fH¢ REFA= CU= XCP- 9 .05759 15 lo ._6927 ._6_Z7 ? B -.13_03 .1_603 1.J2196 ._UZ? .58_27 -L.25598 -3.03222 -.323_9 9,03105 .LSZ83 .03699 .03699 .38359 .3_359 6.08g&2 _.03105 l.eZZO2 l._lZO_ -.0719? -.1737_ 5 6 .13603 16 .569z? .1126B .11200 ,17376 .Ollg? .16666 .10064 12 1_ ,b_927 .08201 ,1979B .1126¢ .ll2_d ,1o666 6,08g_Z b.0626g .66¢_ .SbOOb -.19198 -.06201 -.6_633 ._633 IU 11 IUTAL 1.01#56 1.07756 ,Z3116 -.6592_ -.Z1308 _EFO= J.330_ k_FL= 0.0000 kEFL = 36.5000 -.OOUO .0000 .003? -L.5_IC 109 VELC_||IE$ PANEL NO, A 2 UN WLNG UPPEK $UKFACLt MACH-2,0AU ALPHA- 0.000 VU_TEX STreNgTH AX|AL VELUCITY LAIEKAL VEL_C|T¥ VEKTIGAL VELOCITY -,00000 -,00000 -,I_191 o023Z5 .1q872 -.06J09 ,17879 °05696 J ,00000 .00695 -,03668 ,03265 b • 00000 oOUOUU ,00518 ,00868 -,02358 -°01760 ,018Z1 ,0_0b -.02061 -o0L522 -,01379 -,02918 6 1 -,OdOOO -o00000 .0|109 °003|§ 8 ,00000 ,0_614 -,0091T -,039_6 9 10 ,00000 ,00000 ,0_86_ ,0£682 -,01025 -,00_3_ -,0_623 -,0670_. Ln LL -o00U00 ,O00UO .022_7 -,1_687 ,0016_ ,15178 -,0_710 ,11879 1_ 1_ -.00000 -.00000 -.U_o_ -.00618 .0_060 .00164 .05696 °03265 E_ -,00000 ,02_67 -,05002 ,01109 16 17 -,00000 ,00000 °01651 o0228_ -°03580 -o0_163 ,00_15 -,01379 18 £9 ,00000 -.00000 *0_992 *03_71 -.042_6 -.06606 -*02918 -*039_6 _0 _I -,OuOOO -.OOOOO ,U3T_ .03_39 -,0619_ -.03595 -,0¢623 -.0_70_ _2 2_ -,00000 --,00000 ,0_269 -.LZ591 -,03250 -.0_710 _ Z5 -,00000 ,0000U -,0350Z -.018_6 26 27 -*00000 -,0U00U -,0093¢ ,0005o 28 -,00000 ,0229_ -,03670 -,01J79 29 JO -,00000 -,OUO00 ,0611_ ,0_3_9 -,07500 -,06_82 -,0291_ -,039_6 _1 _ -,00UOU -,00000 ,U_270 ,O_ZT7 -,0535_ -,05200 -.0_623 -,0_706 33 -.0U000 .06111 -.0_856 -.0_710 _6 Jb -,00000 ,00000 -,12906 -,03766 ,15068 .0_531 ,17879 ,05696 _6 31 _8 -.O000U -,OOOOO -,00000 -,01850 -,00791 -,00023 _01967 +005_0 -,00397 .03_65 ,01709 ,00315 _9 -.00000 .01061 -.OL68H -.0[379 %o 6_ _00000 -,00000 .019_9 _026J_ -.02_9_ -.033_1 -oOZgL_ -,0_966 _ -_00000 ,02815 -,03980 -,0_623 _3 -.00000 o02_09 -.06336 -.0_10_ ¢6 65 60 .00000 -*00000 -.00000 -o0¢263 .1713_ ,06159 -,0¢710 *17879 .USbgb 67 .00000 -,01966 ,02050 68 *OUUOU -,0_057 .00896 _9 50 -,00000 .0000u 51 52 5_ *00000 -,00000 -,00000 5_ 65 110 ,03005 -.1_1_9 -.06152 • 16680 ,0_206 .021Z3 ,0C966 -,00363 .17879 eo_bgb ,03205 ,01709 ,00315 ,032_ ,01709 -.00036 -.012_9 .00315 -.01379 ,0106B o0Z3O6 .02826 -.02250 °°03263 -o0_052 -,02918 --.03966 -.0_6Z3 -.00000 o0_756 -,06298 -,06106 -.00000 ,0_755 -,0_707 -,06710 -.00_87 .00763 O_IVL CYLINDER S|NGULAR|IY INTEGRA[|U/_ ON THE HACH= ALPHA P_INI OF NJNO BUU¥ MIIH PANEL|NG THE UPPER FUR PgE$SURE _5 O_G_EE SAMPL_ SWEEP NACA 65A00_ M[U-WING GASE O|$T_[bU||ON SURFACE 2.0100 = X O.OUO0 Y Z X/C ZY/B Z/C CP CN POINt 2,3073_ 2.3013_ O. O0000 0.00000 .05000 .1_000 ,19228 .19228 0,00000 0.00000 .07831 -.0_365 J [8._5734 Z.30_3_ O.O0000 .25000 .19228 0.00000 -.01353 .01470 -.0003/ .03511 3 19o30352 20.14910 Z.30734 2.30734 0.00000 0.00000 ._5000 .45000 .17228 .19228 0.00000 O.UO000 -.0_427 -.02670 .01512 .0294| -,O001b .00016 °02313 .0[951 4 5 80,995H7 2.30734 0.00000 .55000 .19828 0.00000 -.03900 .04386 .00094 -.00802 6 1 21.8_Z05 22o60823 8.30134 2.30734 0.00000 O,OOUO0 =65000 °75000 .19228 .19228 0.00000 0o00000 -.04500 -.05d37 .05046 .05169 .00173 °00241 -.05172 -.10819 1 8 9 Z3.5_64t 24.3o0_8 Z.30134 8.30734 0.00000 0.00000 .8_000 .95000 .19228 .19228 U, OOOO0 0.00000 -.05418 -,0_851 .05969 .053_3 .00278 .00252 -.16263 -.|9063 9 10 18.82571 19o55072 4°12_68 _.12>68 U.O0000 0.00000 .05000 .15000 .34381 .3_30t 0.00000 0,00000 °02878 .00289 -.48523 -.081_1 11 IZ 20.275O8 _.1Z560 0.00000 ,25000 .34381 0.00000 -.018J9 .03187 -.00079 21.00063 _1.725_9 4o12568 6.L2560 0.00000 O.O0000 °35000 ._5000 .34_81 ._4381 0.00000 0.00000 -.06071 -.03961 .07065 .00864 -.00071 .00037 -.0132_ -.0626_ 14 15 22,_5054 4./_560 O, O0000 .5_000 .34_81 O,OOO00 -.05263 .09_19 .00196 -.14_33 16 11 2_.17_50 23o900*5 24,625_1 4.12560 4.12508 4.12568 0°00000 0.00000 0.00000 .O5000 °75000 .85000 ,_43_1 .34381 o34381 0.00000 0.00000 0.00000 -.06497 -.01210 -.07093 .11258 ,12494 .12292 .00386 .00535 .00573 -.26_96 -.J8515 -.46061 11 18 19 20 ZL 25.35036 4.12_68 0.00000 .9_000 .3_381 0.00000 -.066_7 .11534 .00543 -.5_336 20 ZZ 21.51108 ZZ.07000 6.49507 6.49507 0.00000 O.dO000 .05000 .15000 .54126 .56126 0.00000 0.00000 °14331 .05821 -.19052 -°06950 .02246 .00311 .|3300 .08192 21 22 ZJ Z_ 22.64507 Z3.ZIZO? 6o_9507 6°49507 O.O00OO O°OUOOG .25000 .35000 .541Z6 .5412O O. O. °02754 .00870 -.03_2 -.01157 °00091 .0001_ .0 6108 °02776 23 26 Z_ g6 23.77906 6°_9507 0.00000 .65000 °54E26 0.00000 -°0Z337 °03107 .00017 -.0921_ 25 2_.36606 84.71305 6.49507 6.49507 0.00000 0.00000 .55000 .o5000 .54126 °54126 0.00000 0.00000 -.06915 -.08921 .09194 .11860 .00198 °00401 -.32682 -.48628 26 27 .11279 °11L66 .0048_ .005_1 -,52638 -.58_35 28 29 .10958 -.14_ .00516 °01690 -.63568 .4807[ 30 31 12 tJ 1_ A_ 16 2? Z8 OOUO0 OOO00 .1409l .03/17 -.Z_18 -.06441 .0101/ -.00166 CM i6.16_99 17.611t/ Io II -.0U621 .03/07 CI I -,369_| .t1869 .01712 1 2 13 25.48005 2©.06704 6°49507 6.49507 0°00000 0.00000 .7_000 .05000 .54226 °54126 0,00000 0.00000 -.08_84 -.08398 26.61404 24.16649 6,49507 0.d3_00 0.00000 0.00000 .95000 °05000 .54126 .73651 0.00000 0.00000 -°082_2 .16910 84.51720 8.8_808 0.00000 .15000 .7_651 0.00000 .05497 -.05286 .00237 .19893 32 26.98808 85o_9807 0.03000 0.03808 0.00000 0.00000 ,25000 ._5000 .13051 .73651 0.00000 0.00000 .Ozb08 .00806 -.08507 -.00715 .00062 .00008 .10460 .03552 33 36 25.80967 26.22066 8.03808 8.83008 0°00000 0o00000 .45000 ._000 .7_6§1 */3651 0.00000 0.00000 -.01044 -°03020 .0100_ .02906 °00005 .00062 -.05018 _.15701 35 36 31 _8 3_ 26,6_126 8.83008 0.00000 .65000 .13651 0.00000 -o04427 .04247 .00146 -.24711 37 21.04205 21.45285 8.83808 8.83808 0.00000 0.00000 .1500G .85000 °73652 .7365i 0.00000 0.00000 -.05322 -.05/28 .0§117 .05501 .00219 .00257 -.31873 -.36568 38 39 _0 27.86304 26.58102 8.83808 10.97_8_ 0.00000 0.00000 °95000 .05000 o73_51 .91449 0°00000 0.00000 -°05967 .16990 °05737 -.08500 .00210 .01002 -.40669 ,49082 40 61 26.0554_ 10o97_84 0°00000 .15000 .91_49 0.00000 -°032?9 .00t_8 65 21.12384 27.39225 27.66006 10.97384 10.91386 10.97384 0.00000 0.00000 0.00000 .25000 .35000 .45000 .91649 o91_49 ,91449 0,00000 0°00000 0°00000 °02978 °01342 -.00_83 -.01_90 -°00672 °00262 .00031 .00007 °0000l .09402 .044|9 -°01656 43 44 45 _6 27.9890/ 10.91384 0.00000 .55000 °91649 0°00000 -°02441 .01281 .00026 -.08692 46 28.1_748 28°46589 10o9738_ 10.91386 0.00000 0.00000 .o5000 .75000 .91449 °91449 0.00000 0.00000 -.0_061 -°05268 .02035 .02636 .O00lO °00113 -°15026 -.Z0112 47 60 20.73430 29.00811 10,97384 10°91_86 0°00000 0°00000 .85000 .95000 °9|449 ,91449 0o00000 0.00000 -o050_9 -°05653 .02867 .02828 o0013_ .001_3 -°22549 -°23164 49 50 g9 30 31 3Z _3 3# 35 3o _t _g #? 5O .06594 .19936 42 iii VELUC|IiE$ PANtL ON W/NG LUNLK VOKTEX N_. SUKFAC_t MACH-Z.UIO AXIAL 5TAEN_[H VELOL|IY I -°00000 -.A219X 2 -.00000 °00000 ALPHA= 0.000 LATEkAL VEKI|CAL VELOCITY VELUCLIY ,L6872 -.L7879 ,02_25 .00695 -,06309 -.0366_ -o05O96 -.OJ2b_ °00000 .00_1_ -.023_8 -,01709 5 .00000 ,00_b8 -,01760 -,00J15 6 7 8 -.OdUO0 -o00000 .00000 °01_L ,022U0 .02416 -.0206L -.01_Z -.00917 ,013#9 ,0291_ ,03946 .02864 --*01025 ,04023 .0_502 ,02_47 -.00434 °00_44 °00000 £0 L? .00000 -,00000 1_ 1_ .000U0 -°00000 -._687 -.U3266 14 15 -.0000U -,00000 -.00018 .023_7 lo -.00000 °01651 .02284 .029_2 -,04163 -.04226 .01379 .02918 .0394_ .0462J 17 18 .00000 ,00000 .04104 °04710 °15118 .04060 -.17879 -.05_9b .00164 -.0_00_ -.0_265 -.01709 -,03580 -,00J1_ 19 20 -.00000 -°00000 .03571 .037_3 -.06606 -.04196 21 -.00000 .03_139 -°03595 22 23 -°00000 -.00000 .OJZ09 -,12591 2_ 25 -.00000 .00000 -.03502 -.01u56 20 Z7 -,00000 -.00000 -.00934 .00056 28 29 -.00000 -.000U0 *U2Z9£ ,04716 --.03670 -.07500 .01379 ,02918 JO 31 --*00000 -,00000 .04J49 .06270 -.OOZBZ -.O_SJ .0396b ,04b23 J_ J3 -,00000 -.00000 .0427/ .04L11 -°05200 -.0_850 34 -,00000 J5 _o _7 -.03250 .14680 .04706 .04710 -.17879 .0_20o .02123 -.0_o96 -,OJ265 .0_966 --,003b_ --.01109 --,0031_ °06704 .04710 -,1_906 ,150_0 -°17879 .00000 I.OOOOU -°03766 -.01850 .06531 .01967 -.0509b -.0_205 -.00000 -.00191 °00560 30 _9 -.00000 -.00000 -.0002_ °01001 40 61 .00000 -,OUUUU .01939 *026_4 .028_5 -°01709 -.00397 -.01688 -.00_15 ,01_79 -.02695 -,03_1 .02918 .U3946 -.03980 .0402_ 62 -.00000 4J 46 45 -.00000 .000UU -.00000 .02009 .03005 -.16159 _b -.00000 -.04752 .06159 -.05090 -.0196O -*01057 .02050 ,00896 -.0_26_ -.01709 47 68 .00000 °00000 -.00036 -.01269 .06106 .04710 -.17079 49 50 -.00000 .00000 51 52 .00000 -.00000 ,01048 .02_6_ -*02250 -.0326_ .0291_ .03946 53 -.00000 .OZBZb -.04052 .04o2J 54 55 -°00000 -,00000 .02754 ,0Z755 -,04Z98 -.0_707 .06704 .04710 112 -,00287 .0076_ -.06334 -.06_4_ .L7134 -o00_15 .0L379 OG|¥E GYLINDER BO0¥ $|N6ULAR||Y PANELING INTEGRAT/UN UF ON LOWE_ THE MACH_ THE PRESSURE 45 DEGREE SME_P SAMPLE CASE NACA 05A004 M|O-W|NG D|$TR|SUT|ON SURFACE Z*0100 ALPHAPO|NT WING W|TH FOR 0.0000 X Y L X/C ZY/B Z/C GP CN 1 2 _6.76499 17.61_L1 2o30_34 2.30734 O. UOO00 0.00000 .05000 .15000 .19228 .192Z8 0.00000 0.00000 ,07831 -.03365 .08627 -.0310 3 4 18.65736 19.30352 d._0736 2.30734 0.00000 0.00000 .25000 .35000 .19228 .192Z8 0.00000 0.00000 -.0|353 -.01421 -.01490 -.015TZ 5 6 30736 1 CT CM PO]NI .01017 -.O0|bb .34921 -.[|869 t 2 -.00031 -.00010 -.0351| -.02313 3 6 20.16970 Z. 0.00000 .45000 .19228 0.00000 -.02670 -*0Z961 oO00[O 1 20.99587 21.86205 _.30 134 2°30136 0.00000 0.00000 .55000 .05000 .t9ZZ8 *19228 0.00000 0*00000 -.03980 -o04580 -.04384 -.05046 .00094 .00173 -o0[951 .0080_ .05192 5 6 7 8 9 22o68_23 23.53441 2°30136 2.30336 0.00000 0.00000 .15000 .85000 .19_Z8 .19228 0o00000 0.00000 -°05_J7 -.05618 -.05169 -.05969 .00261 .00218 .[0819 .16_43 10 11 24.38058 18,82511 2.30136 4.12568 0.00000 0.00000 .95000 .05000 .19228 .36381 0.00000 0.00000 -.0_851 °16091 -.05343 .24610 .OOZSz .02810 .19003 .685Z3 10 II 19.55012 4o12568 0o00000 .15000 .34301 0.00000 E3 t_ Z0.215O8 21.00063 6.12568 6.1Z 568 0.00000 0.00000 .25000 .35000 .34381 .36381 0.00000 0°00000 -.01839 -°04011 -°03187 -.07065 _5 _1.7Z559 4.1Z508 O.OOOUO .65000 .34381 0.00000 -°03961 -.06866 16 11 22.65056 23.17550 6.12508 4.|2568 O. UUO00 _.OOOO0 .§50UU .05000 .34381 .34301 O. OOOUO O. OOUO0 -.05263 -.0_497 18 19 23.90045 26°6Z54| 4.12568 4.12568 0o00000 O. OUO00 .#_000 °85000 .34381 .34381 0.00000 O. OOUO0 20 Z[ Z5.35036 21._1108 4.1Z568 6.49507 U. OOUO0 0.00000 .95000 .0_000 °34381 .561Z6 O*OOUO0 0.00000 22 22.07808 6.69507 0.00000 .15000 .561Z6 13 26 2Z.64501 6.69507 0.00000 .25000 .54126 23.Z1¢07 Z3°71900 bo49507 b.49507 O°O00UO O. O000U .35000 .45000 .54126 °5_126 0.00000 0.00000 °00810 -.02331 27 24.3600o 24°91305 6.69501 6.49501 0.00000 O. UO000 .55000 °65000 .56126 .56126 O. OOUO0 0.00000 28 29 30 25.68005 26o06104 6.69507 6.69501 0.00000 0.00000 .15000 .85000 .56126 .56126 26.61406 6,49501 0.00000 .95000 .56_Z6 26.16649 26o5T728 0.83008 0°83808 0.00000 0.00000 .05000 .15000 .73651 .73651 0°00000 0°00000 _149|0 .05497 °|6_35 *05285 _6.9_808 25.39081 8.83008 0.83800 0°00000 0.00000 °25000 .35000 .13651 °73051 0°00000 0.00000 .02608 .00806 o0_501 .00715 25.809ol Z6.22066 0.83808 8.83808 0.00000 0.00000 ._5000 .55000 .13651 .13651 0.00000 0.00000 -.01046 -,03020 -°01004 -.02904 .00005 ,00062 .05018 .|5701 35 36 Z6.03126 0°83800 0.00000 .o5000 .13651 0.00000 -°04611 -.06247 .00146 .Z4711 37 27.06Z05 27.6§Z05 8.83_08 8.03808 0.00000 0.00000 .15000 .85000 .7365| .73651 0o00000 O.OUO00 -°05322 -.057Z8 -.05111 -.05507 .00219 .0025Z .31813 .36568 38 3_ 41 Z7°86366 20.5870_ 8.83808 10.97384 O. UO000 U. OOOOO .95000 .05000 °73651 .91449 0.00000 O. OOOOU -o05901 .10990 -.05137 .08500 .00270 .O|OU_ .60469 -.49082 60 41 6Z 43 26.85563 27.12306 10o91384 10.9738_ O° 00000 0.00000 .15000 .25000 °91469 .91449 0.00000 0.00000 .06596 .0_918 .03299 .01690 ° 00_48 .00037 -.19936 -.09_0Z 62 43 64 65 21.39225 Z?obbOb6 10.97304 £0.97386 0.00000 0.00000 .35000 ._OOG .91449 .91669 0.00000 0.00000 .0136Z --.00483 .00672 -.00Z6Z .00001 .0000| -,06619 .01056 6_ 45 40 47 27°92907 L0°97386 0.00000 .55000 .91669 0o00000 -.02641 -°01221 .00020 .0869Z 66 28.|9748 _8.46589 10.97386 10.97384 0.00000 0°00000 o05000 ,75000 .91469 °91669 0°00000 0°00000 -.06001 --°05268 -°0Z035 -°0Z630 .00010 °0011_ .15026 .Z0112 67 48 28.13630 29.00ZIL 10o97386 10.91384 0.00000 0.00000 .85000 .9500U .91469 °91449 0.00000 0.00000 -.05689 -°05653 -.02847 -.02B2d °00133 .00L33 .22569 .2_L6_ 49 50 Z5 2o 31 3Z 3J 34 35 36 3? 38 39 40 68 49 _U .03711 .00461 .00289 8 9 .08131 |2 -.0[1|2 .013Z6 13 |6 .00031 o06_64 |5 -.09119 -.11258 .00196 .00300 .|6933 .26596 16 17 -°07210 -.01093 -.1Z494 -°12292 .00535 °00573 °38575 .4686| L8 19 -,06051 .14331 -.1_534 .1905_ .00543 .02266 .5Z336 -.13300 ZO 21 0.00000 .05221 .06950 .00311 -.08792 22 0.00000 .02156 o03002 .00091 -.06108 23 .0115T -.0310_ .00012 .00011 -.OZlIb ,09215 24 25 -,06915 -.08921 -.0919_ -.11860 .0019d .00607 .32682 °68628 26 21 0o00000 0°00000 -.08686 -.08398 -.11279 -°11166 .00483 .00521 .52o30 .58435 28 29 0.00000 -.08262 -.1095_ °00510 .03568 30 °01690 ,00237 -°68071 -.19893 3[ 32 .OOU6Z .O000U -°10_68 -.0355Z 33 34 -.00019 -.00071 113 TOTAL cJN THE REFA" REFX= (.|ENTS N|NG ZU.0130 -,OOO0 .O04b CN_ CL" .OOOO -.OUO0 THE REFAREFXa KEFG= 6.8900 kbF¢- b*Bguo 0.0000 GOEFFIG|ENT$ G_MPLETE CONFXGURAT|UN I_6,0000 2UoSL30 CNCT- o0000 .0083 GNCL- -.0000 o0000 CD" L_.O000 s -.LOSo TUTAL XCP- AEFZ ,0046 XCP- ON KEFS[ 1_4,0000 CN" CT- CO- 114 GUEFF| .0U83 -_.S/LG AEFB., REFL= TNE ON IviNU _E_L= CN= _T= -.OuUO o0U34 CM= -.U_UU CL = CO= -.UUOU ,U036 X_P= = _.4UOU CN= -,UdO0 CT= ,uOol LM= LL= LO= .UUUU -,UUO0 .O0_l XCP lboJ4X9 kEeL = b._90U XLE = I_.4633 6.d9UO XLE = 21.2216 = ZJ.9611 -.l_l = UELY = Z.3_OU LN= REFL= -,OU00 CT= CM= °0072 .UUOO LL = bD= -oUU0O °UOIL XCP: -o_ll SbCTIUN UN TME UELY= LUEkFICIENTS _&Nb Z.3/OO KEFL = 0.890G XLE kEkL = b,sgoo XLE- -.OOUO .UObi CN= _T= oUUUU CM= LL= CD= -.OUO0 ,OUbl XCP= -°b32_ UEL¥= 1,9UUO CN= ,UUUU CT= CM- .0065 .0000 CL= CD= XLP= = XLE= ._33_ UELV CP_[Au _,890C Zb._528 oOUOU ,UOb5 _,OJ31 2.4_o5U TAME = IZZ.42_00 TIME = IL4,14_00 TIME = L_l.bO_OO LPCR1T = L,i}OgL LPVAC = -,3_JO0 115 VELUCITIE_ PANEL NUo 1 ON BUDYp HACH-2.010 ALPHA= 30UKGE AX|AL STk_N_TH ._9821 VELOGiIV -.12592 5.000 LATEkAL VEkT|CAL NORMAL VELOG|rY .12762 VELUCIIY -.2700b VELOG|TY .32971 2 • ,23814 ,L5320 -,10863 --,08418 ,24986 ,16672 --,06568 ,10663 ,28260 ,21590 6 ,0931S .296_1 -°06689 -oU9585 ,06488 ,11639 o1_592 -,2JOb1 o168T3 °27665 6 7 ,2_353 ,1_038 -,OBOb5 -.05916 ,81698 ,12165 -,0386B ,10076 ,22928 ,1625B 8 ,0616& -°063?3 ,0_305 9 10 ._5703 ,112J8 -.0518k -.06502 .09712 °16318 11 12 .05207 -.03198 -,025_9 -°01161 ,06234 -.00371 13 °228?2 -,01979 ,0795b |6 1_ °12501 -.02A_0 16 -,18515 t? 18 19 20 21 _2 23 84 25 26 ,1003_ ,0_307 -.098o2 -.20589 .19295 ,078_1 -.08383 -,L98_8 ,2•512 ,1168b -,00985 ,005_6 °01618 .00860 .01566 .02560 .03845 .01318 ,01601 .08000 .02288 ,007_1 .007_U 21 28 --,05666 -,177_2 29 -.1_689 30 31 -,07005 ,12_88 -,03603 .01915 • _ _3 ,15216 ,00168 J6 °00691 3_ 36 ,01o19 -.0537u 37 38 ,0_961 -.0021_ 39 60 -.03708 -.1_940 41 -,10_69 -,01962 _2 6_ -,09_8_ .006_• _4 45 66 ,11910 .09712 -.035_ 41 68 69 ,00901 -,025_6 .0_..06t ,11155 -.00236 -°0_33 ,06631 .07365 -.06907 -.05621 .06679 ,06419 -.06679 -,06679 ,06566 .06566 ,10602 -.|1519 -,00891 ,08651 ,06973 -°11702 ,02886 .01409 -,00781 -°07381 ,0558b .0o590 -,06962 -,06032 ,00925 .06925 -,06032 -,05996 ,01137 ,t1561 .20105 ,15988 ,09318 .06601 ,E3973 ,09856 ,02586 -,02L31 ,09346 .06629 -,020_2 -.06759 ,08052 ,03335 -.03335 -,08058 ,08052 ,03335 ,00703 °00692 -,0656b -,06566 *011_1 -,0599o .00058 -.0020_ -°08800 .00506 -,02128 -.0168b -o03592 °03330 -,03340 ,016_0 -,03_65 ,00204 -.02645 -,00800 -,09813 -,08058 ,08054 -.0J168 -.02_00 -,16499 ,01123 -,01260 -,i16_9 -,08200 -°01938 -.0066b -°09522 -.10331 .00723 -.0005b -,10654 -.0_b89 .03603 -°07226 °08053 -,01116 .04962 ,01966 -.00515 -,0_030 -.0_459 ,03360 -.03330 ,01950 -,03091 -.03bb2 -,03113 .01107 -,010t0 -,07625 -.08255 -.11170 -,08051 .08050 .03342 ,03_55 .03616 -,06515 -.00961 .030|0 .06306 -o03335 -,08052 .08052 .0333| -,03340 -,08050 ,00056 .03338 -.03333 -.00069 ,07161 ,03_b5 -,17065 -,03329 ,06211 ,01866 -,09_91 -,08056 .04121 -,0700_ .08069 -,02•59 50 51 -.0o008 -,0607_ -.00391 ,01962 .03051 -,0|112 -.01_69 -,06012 .03363 --°03388 5_ b# 5_ -°006_1 .lbl6b -°0J317 .05_55 -.U06_6 -,0026b -,00061 .06800 ,06641 -.08696 -.05875 ,02651 --,08055 .08053 ,03349 55 -,016Z9 ,025_1 -,02383 -,02.926 -,03321 50 57 -.09919 ,1_181 ,0_50 ,0025B -,02013 .06929 -,07855 -°05867 -,08051 ,08053 58 59 -,0/301 -,162_1 ,0113b ,00707 ,06502 -°0399b O0 -,1_19_ °01620 -°02790 116 .06939 .00914 -.07559 °03351 -.03319 -.08051 dbiV_ bYL|NBcK _JN_ULAR[TY LNIEbRAT|_N ON TH,F HACHALPHA= 8UUY W|[H PANELING UF I,c FOK P_EbbUKE 45 OEGKEE SAMPLE SWLEP N_A o5A004 M|D-W|NG CASE _ISIR|BUIION _UDY _,OIUO 5.0000 PuINI X Y L XlC 2YIB LIC CP i L.OOUUU .ogbl> -.&311o .02740 ,02875 -°06?42 .23352 Z 3 1o000o0 L.O00OU .zJlib ._llo -.09515 .u9515 .02740 .02140 .0096Z .06942 -.02b75 .02_15 .L7365 °11402 l.oooou 3.2200o .u9575 ,Z?JOd ._31z6 -._59Z_ .0Z740 ,088_ .0Z87_ ,OCZOI .06962 -.19198 ,086?4 .18116 b _.Z20UO o6592_ -.Z7308 .08822 ,1979_ -oUBZOI oiZb74 ? 3.ZzOOo 3.ZZOOb .o_92_ .Z130d °08822 .08822 °19198 .08Z0/ .08Z01 .i9798 .01717 .058_0 9 OoO8_Z .6_633 °16_64 °L3403 -,J2_9 .iL548 ,/bob4 .16b06 °32359 °32359 --°13603 .13403 .06657 .OZ49U 1.01754 -/.3Z796 .lbOb4 .Z47_3 °13403 .L6518 ._359 -.39_79 -.5500b .55006 .Z4743 °Z4743 °_9079 .39879 -.1b_10 .LbSl8 iO II O,ObZ4_ O.O_i4Z IZ LJ o.OBg_/ 9.03105 16 L5 9.0_105 9.03105 .Z1308 .b5928 -1.07754 i°077_s 1.071_6 -._4b_3 °446_3 .4_0_ .55hUb 1._2196 1._279b °55000 .38b_1 CN °05031 CT CN POINI .Ol]b2 .99362 I .UL_U -._LOi? °0_01J .000_ .JOb04 -.20096 2 3 -o01869 .Z7883 .O050b .0593b -o36905 4.86624 4 5 °06153 L.410L9 & °08080 -°069Z0 -.09004 .025Z9 .01911 -.85865 -1.57165 7 8 ,04Z13 4°66908 9 °07_18 -.dZ14b .02429 .009|1 1.06719 -.400_Z 10 l| .01211 .04100 -.0321_ .13514 .00462 .00869 -,6687_ 1.58062 12 13 .00005 -.0_080 .00006 .03933 ._064_ .O000L -.O06LL .00073 °4599U 16 L5 l_ LT 9.d_105 11.0_66b /°3zTgz -L°41563 ,Z474_ .30_08 .k6518 .i7609 °398/9 -.4Z51l -°0_13_ -e01359 .103Z4 -.01860 -.00064 -.O002b 1oZ0759 -.1805/ 16 17 1_ 1i.084_0 i°41_b3 -°58631 o_03b_ .4ZSiL -.11609 -.06685 -.OZbS6 -.00090 -°Z5781 L8 19 Zd /[.0044b ll.0_4_b /.41563 .5_031 *5863? 1.4L503 °_036_ .30368 .40511 .l?o09 ZL /_.b3OgU .37345 .1769b -.sZlzi -.023Z0 -.L0138 0.00000 -o77118 2L _Z Z_ Z# 13._090 l_.o_O_O 1_,63090 .373_5 .37345 ._7345 °4Z/ZL .407Z1 .l?b9b -o/?bQb .l?b_b .SZ lZi -.05041 -.05?hi -.06140 -°09667 .IL046 °19190 0°00000 0.00000 0.00000 -.bQ42b .79_3b 1.3782_ 22 Z3 26 z_ kb.6_U _o ZI io._/U Lo.su3Tu i¢ Io.4_310 Z_ 30 i_._01_ id._b145 .5_9Z7 i.62ZbZ _i 3Z /_.ZO14_ iU.ib165 /.4Zgo2 ._89Z7 33 3_ 20.039|5 ZU.0391_ i. 3> Z0.03915 LogZ_OZ 30 _7 /0.03915 Z1.81090 J8 ZL.Ui¢gU 39 _0 ._927 1°_2_6_ i._gzoZ .549Z -1.422oz -°58927 ._89_? 1.6220Z l .170U9 .SZ_11 -.06477 -.0588Z °O3b?Z .O_05d -.OOLZ5 -.00113 .35647 °78153 _9 ZO .451bi .I7696 -._27Z1 -°0119_ -.02500 0°00000 -°108Z4 25 -.589Z7 .5_9Z1 .4_1bl .SSIbL .4212| ,6272i -.17o9b .1709b -°0_493 -.0J461 -.03031 ,03005 O. O00UO 0.00000 -.13126 .13009 26 2? 1.6ZZOZ ,451bL .1709b SZZb_ -._8_21 ._UO3i .5003/ .1709b .4Z12L *5_Z? /.?Z2bZ .50031 ._0031 .4_IZI ./?096 -1.62262 -._8SZ? ._4902 .b4902 °[?ogb .427Z1 ._Z7 .54902 ,42721 /.s2Zb2 -t.42ZoZ °>490/ °59772 °11690 .17o96 .42721 -.42721 l.*£Zog -._89Z7 ._9/?Z .40721 -.l?bgb 2/°dibgO Z/.dL690 1.4ZZ6_ ._gZ? *_89Z7 lo4ggbZ .b97?Z ._9?72 .42721 °llbgb 6£ 4g Z3.59SO5 Z3._94O_ .50927 i,SigOZ °b464_ °b4663 .1709_ .4Z721 -.627_1 -.ilb9b 43 4_ Z3.59q05 2_.59_ i. .04643 °o4043 ,4Z721 .1/b9o °17696 ._27_1 qb Z5.3_i75 .589Zl ,09539 °17690 -.4_?ZL 4o 41 Z_°3817_ gb.3817_ L._2_Oi L.?ZZbZ ,o9539 ,095)9 .4Z121 °42721 -.11096 .17696 48 69 _5°30115 Z?°d40du .589Z1 .Sd9Z? .o9939 .?bZ?4 °17696 .17696 5u 5i _?.84000 Z?.8_O00 i._ZlbZ Lo4g_bz -._8SZ? °5_9Z? .lb2?q °?0216 5_ 53 ZT.8400d 31.ZOOdO .5_Z? °589z? 1.4Z2OZ -i._ZZo2 .7627_ .o5679 _ _ _L._0OOO _£.LUUUO /.4_b_ t.4ZgbZ 5o _1 31.Z0000 _.75000 5d 34.?5000 1.62ZOZ b9 Od 34.750d0 34.?5000 1.4Z_bZ .5_927 .)bgZ7 TUTAL UN THE K_FA= gEFX_ .589Z7 .50927 SZ2bZ .589_7 °SueZ? glgOZ ._0_27 °58927 o_9_7 -l. ._ZIZL -._27Zl -°1709b .17096 *421Z1 -.011i6 .00049 °0842Z --°0_39U -°02084 *023_8 0.00000 ,I0121 28 .013_9 .07309 0.00000 0,00000 °03468 °18069 29 30 °029_9 °04361 0.00000 O.O00OO .0/_24 .1114Z 3i 3Z °16288 °0b001 0,00000 0.00000 .12b04 .04b90 33 36 -*05572 .04830 U.0000U °03?62 35 -.05800 .I0955 *12152 .2190b 0°00000 0.00000 .09_04 -.2L991 3b 37 ,OZJ4_ 0.00000 -.02357 38 .05b?b .154bl 0°00000 0o00000 --°05698 -°155Z2 39 60 .Ugdio .02420 0.00000 0.00000 -.2130_ -.0613_ 41 4Z .01371 .0b043 0.00000 O.O000U -.20_20 -.18478 63 66 .154bL O. O00OO -.10637 45 .082_1 -.Lg?5_ .OIZ_ .|118b O.O000O O,UOOUO -.33040 -._LiO> 6b _? .427_ -°427Zl -°0735Z °0_484 .L556b .20L05 0°00000 0.00000 o42721 o6Z?Zl -.L?o96 ,L/bgb .009L2 -.031Z8 .01_90 °06?66 0°00000 0°00000 .Llbgo °1769b .42721 -°4Z7_1 -.090i9 °01468 .331o4 .06230 0.00000 0.00000 ._5_7_ .05419 ._7ZI ,4Z?Z/ -.17o96 .£769b -°00172 -.0453b -°00303 .07970 0.00000 0.00000 .03142 -.82189 .85479 ,9_ZUb .l?09b .1/b96 ._Z?_L -,627_1 -.05531 -,00_11 .23_b? -,01308 0,00000 0.00000 -2.4375b .18229 -.58_Z? o95Z05 ,4ZIZI -,I?09b -.0_199 -.0717_ 0.00000 .5d927 1o6ZZ62 ,95205 .95Z05 .427Z1 .L/bgb .0Z943 *0880Z U.UO000 0.00000 -1.4Z26_ -°_gZl °589Z? /,_2ZbZ -L.4ZZ6Z -._09Z? °_SSZ? L.42Zo2 -1.62262 -.50_Z? ._dg_? 1.4_ZbZ -t.4Z_oZ -._Z?2i -ol7b9b .17696 °17696 .6ZTZI .17696 .4Z?ZL °07?74 .00984 .00705 -.065_1 -.01379 .0608_ .0Z789 -.08500 -o0_171 .U?3UZ -.017Z3 -o02L34 -.?iLl4 --1°41697 -.09166 -°3347b -2.33045 -°64/09 .999?9 -.4LOZ2 -LoZ26/9 68 69 50 5L 52 53 54 55 56 57 58 59 bO CucFF|_iCNI_ OU_¥ i_.UCOd ZU.81JO L_= ,052_ LT= LM= .003b o0U_3 bL" bO= oUZeL .UUSL XLP= -I,4_262 l._ZZbg L.4izOi olUZi KEkO= 3o330G R_kZ- 0.0000 kEFL= _O.SUO0 117 @ELLE. iT|E WiNG PANEL NU. VUMIcX SIRcN_IH S CN UPPE_ SUkF_LEp MACh=Z.OIO AXIAL VELuCiIY -o0099J ALPHA= LAIEkAL VELUCII¥ °01805 5.000 VERI|CAL VELOCIIY 1 o_L_Oo Z oZUOz2 • Io9U_ oiZ_Z8 o091_ -.180Z0 -.135_4 -,O_OZU -°05450 ,09163 4 5 .IZOJU .OlI_b .ObbSJ .0_@55 --.09981 -.07099 -.07007 -.0_601 ? .0_636 • UoZU? .O_OZJ .O_£b -.05336 --.0567Z -.llb3J -.|Zbb_ 9 .U1L_Ù ,0b_69 -.ObUg? -,L3_3B iO iI ,0d669 .LLUI6 .0091£ ,0Z769 -.05953 -,06001 -,136[9 -.136Zb IZ I3 ,1o_7 ,106u9 -.0_39d oUSO_b .0_506 -.05625 .091bJ -.oJO20 L4 A5 °17415 .1o9_4 .08080 .L0_55 -.099_1 -°16900 -°05650 -.01007 Io 1/ ,159Z_ .i_b_? .0961Z ,ugub7 -.12989 -.IZb63 -.08601 -,IO09b L8 .lU36b .08LbJ -,ll_Z6 -,1_633 _9 LU .0o9i7 .O_19 ,U7059 .06391 -.10Z_9 -.09_04 -.IZboI -.13338 2L _Z .06boi .062_3 .0_¥08 .05_3 -,0_101 -.08265 -.136_9 -.1_6_6 _ .169Zb Z_ Z5 ,lS90L .1507L ,0_965 °0_676 -,06683 -o0_057 _o ZI °L_Z ,15o18 .Ub7g9 ,07_91 -.01965 -.09635 -,07007 -.0840L Z_ 29 ,Io070 .iOU_U .I032Z .ILl30 -.12910 -.16727 -.I0095 -.[_633 _0 °l_o60 oLZL67 -olb3b7 -.120bl 31 JZ .14868 .13_£5 .iib9_ ,10B89 -.16I?l -.13537 -.1_338 -,13619 3_ 36 .ILZLL ,i63Jl ,097Lb -.05763 -.lZ65O .067L2 -.136Zb o09163 _5 3b .1434£ .I43_0 .03606 ,0_ -.0_833 -,066|Z -°03020 --*05650 _? b8 ,16624 *iS614 .Obs/! *OlZiL -.01806 -*088d5 -.0?001 -*08601 _9 °14340 .d8331 -.I0167 -.10095 _0 4L *166bb .1_815 .09Zb9 .09839 --*|I20Z -,1i893 --.I1635 --.l_bbl 6Z ./6990 .10_07 -.I_600 -.133_8 6_ _6 .I_LgB .15_98 .I0406 .iO/bl -.13016 -.IZ978 -o1_619 -.1_42b 65 40 *16UZZ .160_ .08951 -.0203_ .09L6] -.O_OZO 41 6_ .140_6 *160_b *050B0 *U59_6 -*06147 -*07316 -.05650 -*07007 49 .16119 .06769 -.08Zbl -*08601 50 bl .16Ibl .16I_2 .07836 .0813T -.096_8 -.10503 -.10095 -.1163_ _ b_ .I_Zi_ .16Z5_ *09669 .09951 -*11528 -.12328 -.IZbbl -.13338 54 .16Z09 .09895 -.1Z584 -.i3619 >5 ,16_Z5 .09913 -°_3002 -.£362b i18 -,05133 -.01161 .02Z68 .0_979 .09163 -.U3020 -.0_650 OGIVE _YL|NOE& SINGULA_IIY UF INTkGKAT|UN ON ThE MACH WINb = P_ESSURE 65 DE_&EE SNEEP SAMPLE CASE NACA bSAOO4 MIU-NING OI$lk|BUTI_N SURFACE 5.00OU X 1 Z N[TH FOR ¢.OlOO ALPHAPOINT ThE UPPEk BUOY PANELING Y L X/C 2Y/B Z/C CP CN CT CM PUIN! 10.7b_99 I7.61/i/ [8._5734 Z._07_ 2.J07_4 _.30_36 0.00000 O. O000O O.OOOOO .USOOO .iSO00 .25000 .19_2d .19228 .[9228 O.OOO00 O.OOOO0 O.O00OO -.1095B -.189ZZ -.[4[_0 .|ZO?L .ZO8_6 .[Sb[O -.Ot_Z3 -°00934 -.00368 5 _9._035Z _O.L_910 Z._O?_ Z°30136 O. O000O O.OOOO0 ._50OO .65000 .[9228 .IgZP_ O. O000O O.OOUUO -.0997[ -.01159 °[09_4 °07881 -BOO111 .00042 b 20.99581 _.30734 O. OOO00 .55000 .L9228 O.OOO00 -.Ob76Z .016_7 .00160 -.0[J58 6 1 B 2[.86Z05 22.@882_ _._073_ _.30734 0.00000 O. OOOUU .65000 .75000 .19228 ._92ZB O. OOOO0 O.OOOOO -.085[1 -.[0645 .09_83 .iLia? .0032_ .00502 -.09655 -,21991 8 ? lO 23.§34_[ _4°_805B 2.30734 _._0_6 0.00000 O°O00OO ._SOOO .9_OOO °19228 °igZ_ O. O00UO O°O000O -.1[8i4 -.[_JO .130_6 .i_[_ .OObOl °0066_ -.354i1 -.5062_ 11 lB*OL5?l _.IZ50_ 0°00000 .05000 .36_BI 0.00000 -.01494 .0258_ -.00305 .05143 il 19o5501_ 20.27508 4.1_5o8 4.[Z568 0.00000 O. OOO00 .15000 .ZSO00 .3_81 .3438l O.OOOOO O.O00OO -.t1151 -.1662b .20313 .28810 -*009i3 -.00716 .25711 °|5480 12 [3 _i°O00b_ Z1.72559 _./ZSoB 4.[ZSb8 O. O000O O. OOOOU .35000 .65000 *36381 .3438_ 0.00000 0.00000 -.1791[ -.16501 *31036 .Z859J -.00314 .00152 2Z._5056 _.iZSO8 O°OOOOO °55000 .3638l O.OUOUO -.[5J67 .ZbbZ7 2_°[755U 23.90065 6. i256d 6.1_568 0o00000 O.OOOOO .o5000 .7§000 .3638[ ._638L O.O000O O. OOOO0 -.13?bO -.[_3i5 .2J84# .2L339 Z4°bZ56[ Z5°35036 6. lZ568 6.lZ568 0.00000 0°00000 .85000 .95000 .36381 °3638l O.O00OO O.OOO00 -.l120& -.[030_ °[9416 .i_B53 _l._[lOU @._950_ O,O000O .05OOO .54l_6 O°O0000 2_ 22.07808 ZP°OB_Ol 6.69507 0.69507 O.OOOOO O.O0000 .[9000 .ZSOOU .56LZ6 .5412b O. OOO00 O.OOOOO -.08677 -.IIOSU 25 Z_.PLZO7 ZJ.I/9U_ 6.49501 o.69_U7 U. OOOO0 U. UOUOU .35000 .45000 .56126 .5_lZb O.OOOOO u.o0000 Z6 27 24°34600 Z6.9[_0_ 6._950_ b°69501 O. O000O O.O00OO °5SOU0 .6_000 .54126 .561Z6 29 _5.48005 26.06706 b°69bO1 o.6950# O. OOOOO O.OOOO0 .75000 .8500U 30 Z_,6L606 0.4950_ O.O000O JZ 26.[0069 ZB.SIIZ_ 8.a_Boa 8.83808 O. O0000 O.OOO00 15 ll 18 19 20 2¢ 22 .68866 .bb/40 .3@7_Z 1 2 3 .16580 .0523/ 5 9 IO -.058Z3 -*Zb094 [_ _5 °0057Z -.43603 16 .OOOl8 .00914 -°563_0 -.o5883 17 18 .00905 .008_0 -°740[Z -.81006 19 2O .00040 .00239 Z[ .LI535 .|470[ -.00517 -.00366 -.16593 -.26936 ZZ Z3 -.12933 -.15894 .11196 .2[131 -.00174 ,OOllz -.41251 -.6261b 26 25 O.OOOOO O. OO00O -.19700 -._[[Oi .Z6191 .2805_ .00503 .00903 -.9253_ -|.150[? 26 _1 °561_b .56[26 O.O0000 O.O00OO -.20309 -.19389 .27001 .25171 .0II57 .OIZOZ -t. ZbOi3 -1.349[b 28 29 .95000 .56[Zb 0.00000 -._79_9 .OtLZb -1.38137 30 .OSOOO .[5000 .7305[ .l_bSl O.OO00O 0°00000 .0[370 -°078_4 -.013Z3 .O_5?O .OOlSo -.00339 .0443b -.28496 31 32 °0025_ -.0036_ °2_916 _4°9B_0_ ZS._YS_I _._UO_ 8.8_80a 0°00000 O.OOOOO .2_000 °_5000 °?3b_ .73051 0.00000 O.O00OO -°_0_07 -°[_i21 .[0[02 .[[bb6 -.002_| -.O0_L8 -°6217_ -.53_Z 33 36 J5 _6 ZS.80961 _O.ZZO46 8.a_808 8.B_OU O°OO00O O°00000 .45000 ._UOO .7_65L .?_b5[ 0.00000 O. OOOOO -.|373b -.15_4_ .[_20b .1685Z .00070 .00319 -.6598_ -.80_13 35 36 J1 ZB°6_I_O 8.83608 O°OOO00 °bSOOO °?JbSI O. O0000 -°[oOB6 .[6043 .00_5l -.93341 37 B.U3UU8 8.83_0_ O.OOOOO 0°00000 .75000 oUSO00 .13651 .73651 0.00000 0°00000 -.[1560 -.IBOZ5 ./0863 .17330 .OOlZ2 .ooBo8 -1.05061 -1,15065 38 39 .lZbBb -.Oigi_ .O0831 .OUZZ5 -1.2_413 .[IOS_ 60 4L }9 27.U_20_ Z1.452_5 _0 _L Z1.86364 2O.SB?OZ 0.83_0_ [U.973_4 O.OOOOO O.UO000 .?_OOO .05000 .73651 .91669 O.O00OO 0.00000 -.L8353 .O_BZ3 4Z _J Zb. I0.97384 O°OOOOO .lbOOO .9[669 0°00000 -.06559 .0328Z -.OOk41 -.1983[ 42 4_ 4b Z?°IP_H_ ZT.39ZZ5 27.o0060 [0°91384 io.973a4 iU.97386 O.OO00O O.OOOOO O.OOOOO ._5000 .35BOO .4500G °91469 .91649 .91669 O°O000O O.O0000 O.O0000 -.09_9H -.1[357 -.12942 °U69_Z .05682 .06675 -°OOL_3 -.00051 .00034 -°3iZSZ -.313_6 -.44341 43 66 _5 4! zT°gggo/ ZB.IgZga 10.913_4 I0.97_84 O.OOOOO 0.00000 ._OOO .bSO00 .91449 .91669 Ù.O000O O. O00UO -.14bO? -.1_918 .07308 .07994 .001_? .OOZ14 -.5200b -._903Z 4b 67 Z0._6589 10.97384 O. OOOOO .15000 .91_69 0.00000 -.LTUOO .08505 .U0364 -.OSO9i 48 28.73430 Z9.00271 [0.91386 £0.9Z384 0.00000 0°00000 ._000 °95000 .91449 .91469 O.OOOOO O.O000O -.1139Z -°[76Zl .O_102 °0811b .0060o .00610 -.68928 -.11385 49 50 _d _9 50 Bg_43 119 VELOCIT|E_ PANEL NO, 120 _ WING LOWE& VOkTcX STK_NGTH SURFACEr AXIAL ¥ELOC|[¥ MACH=Z.010 ALPHA= LATERAL VELGC[iY 5.000 VERT]CAL VELOC[[Y k .ZZ608 -.ZJ_01 .27948 -.Zb595 3 4 .16902 .12030 -°01158 -.05491 .06587 o0525B -.11981 -oL0624 5 6 .07135 .03408 -°02?00 o00115 ,03570 °01651 -.09031 -.0733b 7 a .OJb3b ,06207 °003_7 -°00691 ,0Z280 °03830 -.05_90 -°04710 9 10 *OIAUO eOSb69 -°00131 -*01751 °040_9 °05091 -.04093 --.040|2 11 1_ .1101_ olbSOF -°03265 -,Z0986 o06299 °24851 -.0400b -°26595 13 .Lboo9 -.11573 .1375z -.1¢41z 16 15 .174£5 ,t6904 -,09_29 -°06128 ,10285 ,04905 -.1198L -.LO_2_ lb ,159Z2 -,06_11 *05825 -,09031 17 L8 °13b_1 ,10_43 -o04_90 -.0Z119 .04213 .0Z868 -.07330 -.05798 19 ZU oU6977 ,06319 .00082 .01012 ,01455 *0ILL5 -°04770 -.04093 21 22 oO_bbl .04233 ,OLEO/ °01151 .01517 .01750 -,0_01/ -.U400_ 2_ ,16925 -.20058 .Z$392 -.26595 2_ 25 .1_901 .15011 -.1095b -.09_95 .12903 .109L0 -,14412 -.11981 26 27 o15332 .l_b7_ -.08b03 -.0718b ,09883 °0871b -,10424 -.09031 28 29 .L_070 *16040 -,05748 -.U3310 .05580 ,01739 -.07J36 -.0519U 30 31 .15640 ,|_8 -,03_73 -.0315_ .0280b °03_o_ -*06770 -.06093 _2 .13225 -.023_b °03135 -°06012 33 3_ ,11211 .14331 -.01695 -*20075 °02941 °23432 -°04006 -,Z6595 35 3b .16342 *1_370 -.1093_ -.09038 *12899 .|0350 -*16_12 -.11981 31 38 .14424 o1_476 -.08007 -°07263 .08952 .08038 -.10624 -*09031 39 .14546 -.06215 .06778 -.07_36 60 41 ,t_b6b .1_8i5 -,05397 -.0497b ,05819 .05237 -,05798 -.04770 42 .14990 -.04683 ,0404] -°Uq093 43 • _ _5 e15198 ,15J98 ,14022 -°06792 -,04636 -,21168 °04356 .04500 ,25309 -°04012 -,O_UOb -,26595 46 • 7 • 8 .1_038 .14064 .160_5 -,11169 -.00974 -,08102 ,143_ .10268 .09111 -.14412 -.11981 -.10_24 • 9 50 .14119 .14151 -.073_9 -°06315 .08196 ,06996 -°090_1 -.07336 51 .14102 -.05445 .ObO07 -,U5798 b2 53 ,1_215 ,14Z54 -.041_6 -.04303 .05006 °042.3l -.04770 --,04093 54 .L4289 -.04_93 *03996 -.0401_ 55 ,1_325 -.04412 ,03598 -,04006 UGIVc CYL|NOLK _INbULAKAI¥ OF INILGKATIulS ON TH£ _ANd It6 LOdC_ _5 FO_ P_tSbOgk _R_ S_EEP SAMPC_ NACA 05A004 MIO-WING CASE dISIK[60IlUN bUgFACC Z.OIUO 5._000 HACH= ALPHA= Y g PuINT sllff _UOY YAN_LING Z X/C ZY/O L/C CP CN POINT L,3LOZb L .173B_ *lgiS7 .00858 .6[305 2 ,L49_9 .09_5_ .gb47_ .[U_O_ .OO410 .OOIO_ .388ll .1555_ 3 4 ,U_39_ .ddi6_ .O_?_d .001_5 -.OOO_O -.00004 .UZ_80 -,000_6 5 6 .00_1 ,O[865 °00946 .0205_ -,OOO3Z -.00088 -.00971 -.0385i I 8 O. OOdO0 O.GO000 .028dO .05_3_ .OJ[1Z ,05990 -.00168 -.0028/ -,0_633 -,ZlJI[ .3438L 0.00000 .316bO .545L$ ,Ob4ZO .i_O00 ._5000 .34_81 .34JSL O.OOO00 0.00000 ,ZZ05_ .l?_5 .JqZS/ °Z98OZ ,0[759 .O074J U. UOOO0 U°OUUO O ,35000 .6_000 .34_81 ._4_8L O. OO000 0.00000 ,i_iO7 .LZZ_I .Z4_4_ .ZI22L ,OUZ41 -.UOiL_ -.0_7 -.ig_bb L_ L5 _,£_o_ O. .55000 ,34_8L 0.00000 .Ol_(_ .L_Zg/ -.002_0 -.LLII4 L6 _,IZ_O_ 4.IZ_OO O. OOOOO U,OO000 .o_OOd .1500U ._4J8L ,3638L U.0OOOO O. OOO00 .02777 -.00_4 .0481Z -.OzOLl -.00L65 .000_3 -.LL367 .d3iZ2 L? L8 Z._ul_ O.UOO00 .05000 ,19ZZ_ O,OUO00 ,Z95_1 Z l?.Ollg/ Z._U/36 O. OOOOO .15000 .19Z/U Q,UUUUO 1_,q5136 I_,_U35_ Z,_O{_4 Z,JUi36 U,OUOO0 O, O0000 ._)OOO ._5000 ,L9_2_ .19ZZ8 O. UOOO0 d. OO000 b _U.L_9?U ZU,9_5_l _._0136 Z.JulJ6 O.OOOdO O.UOOOU ._000 ._UOO ,192_ .19Z28 O. O. I u Z[._420) _£.o_2J _._U?J_ Z,JO7_ 0.00000 O,UOOUO .b_O00 .15000 .[9Z28 ,[9_8 O,UOUUO O.O00OO 9 Z_.53S41 z4.Jo05_ Z,_U?_ _.JU/34 O.0000G O.O00OO ._bO00 .9_000 .i9228 .19Z_8 i_.dZb_? _.i_SbO U. ,05000 iZ lg.5507Z 2U,Z_SO_ 4,1_508 _.tZ_oo U. O0000 G°OOOUO 14 ZL.O00O3 ZLolZ_5_ _.iZbu_ _,[g)O0 _Z°450_9 Z_.lIS_O Z_.9006) iO CT ,0383_ L6,1o_99 ZU II CN .32516 I OOOdO O00UO OdOOO OOdO0 l. OB33L ,4955_ .i605b 9 tO il IZ L] i_.O_Ssl s, U,UOOUO .O500U ,363_L 0.00000 -.OLo3L -.02_25 .OO/JZ _J ZL ZS.J503b _.SilOd 4.1i_o8 0._9507 U. OOOUU 0,00000 .9_000 .05OOb ,_6Jdl ,5_t2b U. OOUO0 O.UOdO0 -,01120 ._0191 -.0298D ,_01_0 .OOL6U .041_2 .L_SZO -.Z8025 ZO Zl ZZ ZZ.OIBO_ Z_.b4501 _.4950? o,_gbO? O.OOOO0 U,UUUOU ,l_OOU ,Z_O00 ,54LZ6 ,_LZb O.UO000 O.UUUOO ,2ZOi9 ,L_/_L .29214 .Z625_ .OLJLZ ,OObbJ -.370_3 -,48L0_ 22 Z3 tJ.iLZUI _,_9)U_ O.OJO00 ,_SOOO ,)_i_b J. ._81Z7 .Z_099 ,UU2_4 -,97815 2# Z) ZO ZJ.?igo_ Z4,_46d0 Z_,9L_U5 C.49)0? 0._95UI _._9501 U,OOOJO O.OOOOO O, OOO00 .S)UQO .55000 ,bSO00 ._iido ._4t_6 ,5_LZ6 O._UOVO U.OUdOO O,UOO00 ,LSLZb .LO3LI .07770 ,20109 ,i_/O_ ,id33d -,OOtO? -,00Z95 -,00354 -._9645 -.4843_ -,4Z354 25 26 27 Ld _5._Ud9 o°_9_0_ O, UOO00 .75000 .54L26 O.OO000 .OI_?Z .0993_ -,O04Zb -.40_bZ Z8 2V JO Zb,O_Z04 Z6.bL404 _._9_07 0._9507 O.UUOOU U. OOUO_ ,o_OOO .9_OUU ._LZ6 .5_gb O.O0000 OaOOOO0 .OoLSg .0_4 .OOZL_ .Dg_bd -.003_3 --*OOZe6 -,43OLd -,3404U 29 _0 _Z Z_,IbU_9 o._300_ O.OOOOO .0_000 o?JO_i U._UOUO .bOl_z ,ZgUL_ .U$42U -,97310 3L 33 g4,_llZd _4.9_=0o 0.0_000 _._3U08 O. UOUO0 O.O0000 .L_O00 .zSOOO .7_051 .7365L O*O00OO O,OOdO0 .2ibJ3 .l_??l *Z0799 .IO0_/ .OU9_Z .0064_ -.78Z93 -,75367 3Z 33 J6 Z_°_9Od/ _*SJoO_ O.OO000 ._OOU .P3bSL 0.00000 .i09_O .ioggO ,OOLO5 -,T4733 J4 Z5,809o7 £o,ZZU4b 8.0_0_ U.O3_O_ u,ddO00 O.OOOO0 .s_O00 ._)0OU .?JoSl .?_OSL O, O00UO O,OOOOO ,I_O_ .LZ95Z .t_Z .LZ_SZ -.OOOll -.OOZO_ -./_Lb_ -.b?_30 $5 36 Zo,o_l_b zl,O4ZO_ O._3uO_ 0.83000 O,uO00O O, udOOO .o_OOO .15000 .?_OSL .1365[ O.OOOO0 O,OOOO0 ,£L53/ ,LOILL .LIOgZ .10298 -.003_i -.00441 -.6_530 -.641_o 37 38 _.8_U_ o.8_uOo -,00471 -.O04lJ _J z4 31 _o 39 _U _1.45Z_5 _l.do_4 LgbbU OOUO0 ,LOIIZ L9 O. OOOOO U.UUOUU ._SOOO ,950OU .7_oS1 .730_L 0,00000 O, OOOO0 .t0512 ,LO_o_ ,iOIO0 *ZOO59 -.o/iO_ -.?0920 39 40 41 2O°201OZ Zb.055_$ ZU._/_8_ LU,9?3U_ O. O000O O,dO_OO .oSO00 .LSO00 .9/_49 .9£_9 0.00000 O*O0000 ,$iO_ ,_349 .l_3Z .LLL_Z ,0L866 ,0050/ -.gi4L4 -.b7565 4L 4Z 6J _?,LZ_d_ ZI,39ZL_ I0.97_04 [U._lJ8_ O.O00OO U. uOOO0 .Z_O00 ,35000 .9L6_9 .91_4g O,OOOO0 0._0000 ,iOBO_ *tlib_ .0940U .O_SlJ .OOZ_ .OOO8? -.59_1_ -,50603 43 4# _0 41 ZI,bbUOO zI,gz9oI /O.91J84 10._138_ U.OOUO0 0.00000 ,_5000 ,SHOO0 .914_9 ,91_#9 O.O00UO O*O0000 ,/5ZO_ ,LJ09_ .Olb09 .00551 -.O004L -,OOL41 -.SZLO6 -.4bbil q5 qb _8.191_ lu,913_ O. UOOUO .bSO00 .91649 O.0OOO0 .LLJZ9 .05bb_ -,00195 -,4[859 67 sd ZO,40589 I0.97_0_ O. OOOOU .l_OOO .gi_g U.O0000 ,IOOSi .050/9 -.OOiI> -,38490 48 5U ZO,?_)O Zg,OOZ?i i0.9130_ iO°97_d_ U.OOOUO U°UOOOU .d5000 .95000 .91449 ,9L449 O. OUUO0 O.OUOUO ,0966b .0981W .0483b ,OWglO -.O02Zb -*OOZ_i -.383_0 -.WO2LI 49 50 121 _'U IAL ON THE _EFA= CI_L F F l,(., 1EN _¢_.0000 KEFX= 20oSt3U CN= .Agb9 C1= CH= CL= CO s THE KEFA= b.8900 COCFFAG|ENIS _OMPLc_E CUNFL_URA[I_N l_4oO000 ZU.SL30 GN= °Z_95 CT= CM= °U081 -.OOSI CO= XCP= 0.0000 REFG- .L95T .0217 AEFX= CL = A2.000C REPL= -.3bOO TOTAL ON REFU= .00¢6 -.0105 XCp= 122 i'S iCING .2¢79 .0298 -°2628 REF_= [Z°O000 KEFZ= 0.0000 REFC= 6.8900 fUN $ECI UN tHE _ELY= GN= o1V7_ CT= CM= °0037 ,035b CL- .L963 2,_GO0 CN= GT= XLE- 16.3419 REFL= 6.890C XLE= 1B._633 kEFL= 6.8900 XLE" 2[.2276 REFL= 6.890G XLE= 23.961[ 6.B900 XLE .2_05 ,U063 CL= CO- .2Z91 ,0_03 XCP= -°0590 OELY= Z.3600 CN= .28_3 C[= CH= ,0069 -,1_99 CL= CO- .2_6 .OJLU X_p= -,4563 SECTION ON THE OELY= COEFFICIENTS N_N_ _.3700 GN" ,28_1 CT= CM= ,0058 -._L4G CL = CO= XCP= .2825 ,0305 -°7577 _ELY = [.9000 _N= .Z73Z El= _M= ,006Z -.Z702 CL= CO= KEFL = = .271_ .U300 XCP= = 6.8900 -.013_ CM- T_ME REFL" ,0209 .1812 O_LY= " ENT$ 1.30_0 CO= XCP= CPSIAG LI CUEFF! _ING -1.0111 2,@5050 _PCkiI - 1.13092 _PVAC = -.35360 150.55900 123 REFERENCES • • • • • • • • Woodward, F. A., Tinoco, and Design of Supersonic Flow Properties in the August, 1967. Woodward, tions, at Aircraft, F. A.; Subsonic Vol. 5, E. N., and Larsen, J. W.; Analysis Wing-Body Combinations, Including Near Field• NASA CR-73106, Analysis and and Supersonic No. 6, Nov.- Craidon, C. B.) Description for Airplane Configuration September, 1970. Gothert, B.; Subsonic Speeds. Plane and NACA Labrujere, imate Method T. on Wing-Body Conference Combinations Proceedings Fox, for NASA a C. H. Jr.; Family of TM X-2439, S. H.; Peterson, istics of 1950. No. Experimental Axisymmetric December, 1105, TM Flow Slooff, J. of Pressure Surface Bodies 1971. at Distribution Mach The Boundary-Layer 64A010 Airfoil 10. Stevens, W. A., Goradia, S. H., matical Model for Two-Dimensional in Viscous Flow. NASA CR-1843, Ii. Lamar, Tunnel Sweep J. E., and McKinney, Investigation of a Pressure Wing Model• Program at High 1946. at Subcritical 71, September, at Combinaof Computer X-2074, W.; An ApproxDistribution Speeds• 1970. AGARD Pressure Subsonic Landrum, E. J.; Drag Bodies of Revolution TN D-3163, December, Pressure R. F.; the NACA Digital NASA Three-Dimensional TM of Elliptic Cross Section E8K05, January, 1949. • of a Plots• E., Loeve, W., and for the Calculation Harris, R. V. Jr., and of a Series of Low-Drag from 0.6 to 4.0. NASA Maslen, Design of Wing-Body Speeds• Journal Dec., 1968. on Distributions Speeds• Characteristics at Mach Numbers 1965. Thin Conical Number 1.89. and Section. Stalling NACA and Braden, J. Multi-Component July, 1971. L. W.; Half-Span NASA TN NACA Body RM CharacterTN 2235, A.; MatheAirfoils Low Speed Static Wind Fuselage and Variable D-6215, August, 1971. 125 12. 13. Carlson, H. W.; of Highly Twist and Gapcynski, a Pressure J. P., and Distribution of a 45 degree Angles of Attack November, 1958. 14. 12 6 Pressure on a Series Degrees of Distributions Swept Camber. Arrow NASA Landrum, E. Investigation Sweptback-Wing and Sideslip. at Mach Number Wings Employing TN D-1264, May, J.; Airplane NASA Tabulated Data at Mach Number from 2.01 Model at Combined MEMO 10-15-58L, Hess, J. L., and Smith, A. M. 0.; Calculation Potential Flow about Arbitrary Three-Dimensional Douglas Aircraft Company Report, No. ES 40622, _.S. 2.05 Various 1962. GOVERNMENT PRINTING of Nonlifting Bodies. March, 1962. OFFICE: 1973-739-027/31