Numerical experiments in the stability of leading edge boundary layer flow. A two-dimensional study

A numerical study is performed in order to gain insight to the stability of the infinite swept attachment line boundary layer. The basic flow is taken to be of the Hiemenz class with an added cross-flow giving rise to a constant thickness boundary layer along the attachment line. The full Navier-Stokes equations are solved using an initial value problem approach after two-dimensional perturbations of varying amplitude are introduced into the basic flow. A second-order-accurate finite difference scheme is used in the normal-to-the-wall direction, while a pseudospectral approach is employed in the other directions; temporally, an implicit Crank-Nicolson scheme is used. Extensive use of the efficient fast Fourier transform (FFT) algorithm has been made, resulting in substantial savings in computing cost. Results for the two-dimensional linear regime of perturbations are in very good agreement with past numerical and theoretical investigations, without the need for specific assumptions used by the latter, thus establishing the generality of our method.