A direct method of calculating through separated regions in boundary layer flow

A direct, finite-difference method is presented for calculating through regions of recirculation in laminar boundary layer flow. Analysis is based primarily on steady boundary layer equations, but the concept of incompressible 'Bernoulli interaction' is introduced to let the pressure gradient depend on the normal component of velocity at the edge of the boundary layer. This appears to eliminate the separation point singularity and results in a shift of the predicted separation point downstream more in line with recent solutions of the Navier-Stokes equations.