Three dimensional flows around airfoils with shocks

The present work describes a mathematical model and numerical scheme for the computer-aided calculation of two- and three-dimensional transonic flow over an isolated yawed wing with oblique shock waves and a trailing vortex sheet. The flow is modeled by the potential equation for irrotational flow which is hyperbolic at supersonic points and elliptic at subsonic points. A coordinate-invariant difference scheme is used in which retarded difference formulas are constructed to conform with the local flow direction. The resulting 'rotated' difference scheme allows complete flexibility in the choice of a coordinate system. Shock waves are located automatically in the form of compression bands spread over a few mesh widths. The scheme has proven to be stable and convergent throughout the transonic range. Calculations have been performed for Mach number up to 1.2 and yaw angles up to 60 deg, the likely operating range of a yawed-wing transport designed to fly at supersonic speeds. Calculations becomes less accurate towards upper end of range, because the difference scheme is first-order accurate in the supersonic range.