A numerical solution of the flow field over a transonic airfoil including strong-shock-induced flow separation

An approximate numerical method has been developed for computing the entire transonic flow field over an airfoil, including the boundary layer, the shock-wave/boundary-layer interaction, and the shock-induced-separation. The inviscid flow, the turbulent boundary layer in the weak interaction zone, and the strong interaction elliptic region are solved separately. The inviscid flow is represented by the solution of the transonic small-disturbance equation, while the boundary layer is solved by the generalized Galerkin method, and the strong interaction region (including the wake) is computed from the time-averaged Navier-Stokes equations. A second-order interaction is computed iteratively using the inviscid pressure distribution as a boundary condition for the boundary layer and Navier-Stokes codes, and the viscous displacement thickness to modify the airfoil geometry for the inviscid computation. The computed pressure distributions and shock-wave profiles are compared with experimental data. Agreement is very good except that the pressure recovery on the trailing edge is greater than the experimental data.