The laminar flow over a slender delta wing at incidence has been extensively studied both experimentally and theoretically using vortex sheet methods. These vortex sheet methods have generally been successful apart from the prediction of the secondary boundary-layer separation induced by the primary vortex. This paper revisits the problem using computational fluid dynamics (CFD) and focusses on the effects of the secondary flow separation. The modelling approach is briefly summarised, and the results are compared with flow measurements and results from vortex sheet methods. The computations show very good agreement with measurements for the surface pressures and total head contours. The results help to understand the complex structure of the leading edge vortex flow, and the associated secondary separation of the boundary layer. They indicate that inviscid mechanisms dominate the larger scale features, and highlight a possible mechanism for the development of an instability in the leading edge vortex sheet.