SUMMARYA combined analytical-numerical method based on a matching asymptotic algorithm is proposed for treating angular (sharp corner or wedge) singularities in the numerical solution of the Navier-Stokes equations. We adopt an asymptotic solution for the local ow around the angular points based on the Stokes ow approximation and a numerical solution for the global ow outside the singular regions using a ÿnite-volume method. The coe cients involved in the analytical solution are iteratively updated by matching both solutions in a small region where the Stokes ow approximation holds. Moreover, an error analysis is derived for this method, which serves as a guideline for the practical implementation. The present method is applied to treat the leading-edge singularity of a semi-inÿnite plate. The e ect of various in uencing factors related to the implementation are evaluated with the help of numerical experiments. The investigation showed that the accuracy of the numerical solution for the ow around the leading edge can be signiÿcantly improved with the present method. The results of the numerical experiments support the error analysis and show the desired properties of the new algorithm, i.e. accuracy, robustness and e ciency. Based on the numerical results for the leading-edge singularity, the validity of various classical approximate models for the ow, such as the Stokes approximation, the inviscid ow model and the boundary layer theory of varying orders are examined. Although the methodology proposed was evaluated for the leading-edge problem, it is generally applicable to all kinds of angular singularities and all kinds of ÿnite-discretization methods.