Numerical calculation of the transonic flow past a swept wing

A numerical method is presented for analyzing the transonic potential flow past a lifting, swept wing. A finite difference approximation to the full potential equation is solved in a coordinate system which is nearly conformally mapped from the physical space in planes parallel to the symmetry plane, and reduces the wing surface to a portion of one boundary of the computational grid. A coordinate invariant, rotated difference scheme is used, and the difference equations are solved by relaxation. The method is capable of treating wings of arbitrary planform and dihedral, although approximations in treating the tips and vortex sheet make its accuracy suspect for wings of small aspect ratio. Comparisons of calculated results with experimental data are shown for examples of both conventional and supercritical transport wings. Agreement is good for both types, but it was found necessary to account for the displacement effect of the boundary layer for the supercritical wing, presumably because of its greater sensitivity to changes in effective geometry.