Chapter 13 Compressible Thin Airfoil Theory
Total Page:16
File Type:pdf, Size:1020Kb
CHAPTER 13 COMPRESSIBLE THIN AIRFOIL THEORY 13.1 COMPRESSIBLE POTENTIAL FLOW 13.1.1 THE FULL POTENTIAL EQUATION In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy from an inviscid, irrotational model of the flow. Recall the equations developed in Chapter 6 governing steady, irrotational, homentropic (¢s0= ) flow in the absence of body forces. ¢l()U = 0 UU ¢P ¢£¥--------------- +0-------- = ¤¦2 l . (13.1) P l a ------ = £¥------ ¤¦ P0 l0 The gradient of the isentropic relation is 2 ¢Pa= ¢l. (13.2) Recall from the development in Chapter 6 that P a – 1 ¢P ¢£¥--- = £¥------------ -------- (13.3) ¤¦l ¤¦a l Using (13.3) the momentum equation becomes a P UU ¢£¥£¥------------ --- + --------------- = 0 (13.4) ¤¦¤¦a – 1 l 2 Substitute (13.2) into the continuity equation and use (13.3). The continuity equa- tion becomes 2 2 U ¢a + ()a – 1 a ¢ U = 0 . (13.5) Equate the Bernoulli integral to free stream conditions. bjc 13.1 5/31/13 Compressible potential flow 2 2 2 2 a UU a' U' a' a – 1 2 ------------ + --------------- ==------------ + -------- ------------£¥1 + ------------M =C T (13.6) a – 1 2 a – 1 2 a – 1¤¦2 ' p t Note that the momentum equation is essentially equivalent to the statement that the stagnation temperature T t is constant throughout the flow. Using (13.6) we can write 2 a UU £¥------------ = h – --------------- (13.7) ¤¦a – 1 t 2 The continuity equation finally becomes UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- = 0 (13.8) ¤¦t 2 ¤¦2 The equations governing compressible, steady, inviscid, irrotational motion reduce to a single equation for the velocity vector U . The irrotationality condition ¢ × U = 0 permits the introduction of a velocity potential. U = ¢\ (13.9) and (13.8) becomes ¢\ ¢\ 2 ¢\ ¢\ ()a – 1 £¥h – ------------------------ ¢ \¢\¢– £¥------------------------ = 0 . (13.10) ¤¦t 2 ¤¦2 For complex body shapes numerical methods are normally used to solve for \ . However the equation is of relatively limited applicability. If the flow is over a thick airfoil or a bluff body for instance then the equation only applies to the sub- sonic Mach number regime at Mach numbers below the range where shocks begin to appear on the body. At high subsonic and supersonic Mach numbers where there are shocks then the homentropic assumption (13.2) breaks down. Equation (13.8) also applies to internal flows without shocks such as fully expanded nozzle flow. 13.1.2 THE NONLINEAR SMALL DISTURBANCE APPROXIMATION In the case of a thin airfoil that only slightly disturbs the flow, equation (13.8) can be simplified using small disturbance theory. 5/31/13 13.2 bjc Compressible potential flow Consider the flow past a thin 3-D airfoil shown below. v y=f(x,z) U' + u U' y x z The velocity field consists of a freestream flow plus a small disturbance UU= ' + u Vv= (13.11) Ww= where . (13.12) uU ' « 1 , v U' « 1 , w U' « 1 Similarly the state variables deviate only slightly from freestream values. PP= ' + P' TT= ' + T' (13.13) ll= ' + l' and . (13.14) aa= ' + a' This decomposition of variables is substituted into equation (13.8). Various terms are 2 2 2 2 UU U' u v w --------------- = --------uU++++----- ----- ------ 2 2 ' 2 2 2 (13.15) ¢ U = ux ++vy wz and bjc 13.3 5/31/13 Compressible potential flow UU ¢£¥--------------- = (u U +++uu vv ww , ¤¦2 x ' x x x (13.16) uyU' +++uuy vvy wwy, uzU' +++uuz vvz wwz ) as well as UU ()a – 1 £¥h – --------------- ¢ U = ¤¦t 2 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ ux + ¤¦¤¦2 2 2 2 (13.17) 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ vy + ¤¦¤¦2 2 2 2 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ wz ¤¦¤¦2 2 2 2 and, finally UU 2 U ¢£¥--------------- = u U +++uu U vv U ww U + ¤¦2 x ' x ' x ' x ' 2 uu U +++u u uvv uww + x ' x x x . (13.18) 2 vuyU' +++vuuy v vy vwwy + 2 wuzU' +++wuuz wvvz w wz Neglect terms in (13.17) and (13.18) that are of third order in the disturbance velocities. Now 5/31/13 13.4 bjc Compressible potential flow UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 £¥£¥U' ()a – 1 ²´ht – ²´-------- + uU' ux + ¤¦¤¦2 2 £¥£¥U' (13.19) ()a – 1 ²´ht – ²´-------- + uU' vy + ¤¦¤¦2 2 £¥£¥U' 2 ()a – 1 ²´ht – ²´-------- + uU' wz – ()uxU' +++uuxU' vvxU' wwxU' – ¤¦¤¦2 uuxU' – vuyU' – wuzU' or UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 (13.20) ()a – 1 ()h' – uU' ()ux ++vy wz – ()uxU' +++uuxU' vvxU' wwxU' –uuxU' – vuyU' – wuzU' 2 Recall that ()a – 1 h' = a' . Equation (13.20) can be rearranged to read UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 2 (13.21) a'()ux ++vy wz – U'ux – ()a + 1 uuxU' – ()vvxU' + wwxU' ()a – 1 ()uvyU' + uwzU' – vuyU' – wuzU' 2 Divide through by a' . bjc 13.5 5/31/13 Compressible potential flow UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 ()a + 1 M' ()1M– ' ux ++vy wz – --------------------------- uux – (13.22) a' M' ---------()()a – 1 ()uvy + uwz ++vuy wuz ++vvx wwx a' This equation contains both linear and quadratic terms in the velocity disturbances and one might expect to be able to neglect the quadratic terms. But note that the first term becomes very small near M' = 1 . Thus in order to maintain the small disturbance approximation at transonic Mach numbers the uux term must be retained. The remaining quadratic terms are small at all Mach numbers and can be dropped. Finally the small disturbance equation is 2 ()a + 1 M' . (13.23) ()1M– ' ux ++vy wz – --------------------------- uux = 0 a' The velocity potential is written in terms of a freestream potential and a distur- bance potential . (13.24) \ = U'x + q()xyz,, The small disturbance equation in terms of the disturbance potential becomes 2 M' (13.25) ()1M– ' qxx ++qyy qzz = ()a + 1 ---------qxqxx a' Equation (13.25) is valid over the whole range of subsonic, transonic and super- sonic Mach numbers. 13.1.3 LINEARIZED POTENTIAL FLOW If we restrict our attention to subsonic and supersonic flow, staying away from Mach numbers close to one, the nonlinear term on the right side of (13.25) can be dropped and the small disturbance potential equation reduces to the linear wave equation. 5/31/13 13.6 bjc Compressible potential flow 2 .` qxx – ()qyy + qzz = 0 (13.26) 2 where ` = M' – 1 . In two dimensions 2 ` qxx – qyy = 0 . (13.27) The general solution of (13.27) can be expressed as a sum of two arbitrary functions q()xy, = Fx()– `y + Gx()+ `y . (13.28) 2 Note that if M' < 1 the 2-D linearized potential equation (13.27) is an elliptic equation that can be rescaled to form Laplace’s equation and (13.28) expresses the solution in terms of conjugate complex variables. In this case the subsonic flow can be analyzed using the methods of complex analysis. Presently we will restrict our attention to the supersonic case. The subsonic case is treated later in the chapter. 2 If M' > 1 then (13.27) is the 2-D wave equation and has solutions of hyperbolic type. Supersonic flow is analyzed using the fact that the properties of the flow are constant along the characteristic lines x ± `y = constan t. The figure below illustrates supersonic flow past a thin airfoil with several characteristics shown. Notice that in the linear approximation the characteristics are all parallel to one another and lie at the Mach angle µ' of the free stream. Information about the flow is carried in the value of the potential assigned to a given characteristic and in the spacing between characteristics for a given flow change. Right-leaning characteristics carry the information about the flow on the upper surface of the wing and left-leaning characteristics carry information about the flow on the lower surface. bjc 13.7 5/31/13 Compressible potential flow x – `y = constant y µ' x M' > 1 µ -1 ' = Sin (1/M') x +const`y = tan All properties of the flow, velocity, pressure, temperature, etc. are constant along the characteristics. Since disturbances only propagate along downstream running characteristics we can write the velocity potential for the upper and lower surfaces as q()xy, = Fx()– `y y>0 . (13.29) q()xy, = Gx()+ `y y<0 Let yfx= () define the coordinates of the upper surface of the wing and ygx= () define the lower surface. The full nonlinear boundary condition on the upper surface is v df ---- = ------ . (13.30) U dx yf= In the spirit of the thin airfoil approximation this boundary condition can be approximated by the linearized form v df -------- = ------ (13.31) U' dx y0= which we can write as 5/31/13 13.8 bjc Compressible potential flow ,q()xy, df -------------------- = U £¥------ (13.32) ,y '¤¦dx y0= or U' df F'()x = –-------- £¥------ . (13.33) ` ¤¦dx On the lower surface the boundary condition is U' dg G'()x = -------- £¥------ . (13.34) ` ¤¦dx In the thin airfoil approximation the airfoil itself is, in effect, collapsed to a line along the x-axis, the velocity potential is extended to the line y0= and the sur- face boundary condition is applied at y0= instead of at the physical airfoil surface. The entire effect of the airfoil on the flow is accounted for by the vertical velocity perturbation generated by the local slope of the wing. The linearized boundary condition is valid on 2-D thin wings and on 3-D wings of that are of “thin planar form”.