Solution of viscous transonic flow over wings

Abstract

Knowledge of the viscous flow about wings is very important in 3-D wing design. In transonic flow about a typical supercritical wing, the viscous effect results in a sizable reduction of the lift-to-drag ratio. The Reynolds number dependence of the flow is not clearly defined, and no known similitude exists that can be used to scale the experimental data for a particular design. Recent advances in computer technology and numerical technique have relieved the difficulty of obtaining a theoretical solution somewhat, but the lack of a proper reliable method of treating the turbulence in a time-averaged Navier-Stokes solution remains the major stumbling block. For this paper, a “zonal” approach has been used for a viscid-inviscid interaction analysis to yield an iterative solution for the viscous flow about wings in the transonic flow regime. The chord Reynolds number considered was of the order of 10 6 and above so that the flow was predominantly turbulent. The inviscid flow field was obtained by solving the 3-D potential flow equation. A parabolic coordinate mapping was used in the computation, in conjunction with a finite volume formulation. A new approximate factorization scheme has been developed for the iterative solution of the inviscid flow. A special far field asymptotic boundary condition that improves the accuracy and convergence of the method was derived. For the 3-D boundary layer calculation, the integral method of Myring-Smith-Stock was extensively modified to make it suitable for the interaction calculation. The effect of wing thickness was taken into account and the 3-D viscous wake was computed. The interaction calculation was formulated with a set of coupling conditions that includes the source flux distribution due to the surface boundary layer on the wing, the flux jump distribution due to the viscous wake, and the effect of the viscous wake curvature. The transpiration boundary conditions have been used for the inviscid flow in the coupled calculation. In addition, a method was devised so that the results of an analysis of the trailing edge strong interaction solution for a 2-D viscous airfoil could be adapted for the normal pressure correction near the trailing edge. The theory has been applied to supercritical wing geometries of practical interest. The converged viscous flow results compare favorably with experimental pressure data.