In this paper, the local singular behavior of Stokes flow is solved near the salient and re-entrant corners by the matching eigenfunction method. The flow in a rectangular and an L-shaped cavity are considered as a model for the flow generated by the motion of the upper lid. The solutions of the Stokes equation in polar coordinates are matched with a velocity vector components obtained by analytic or numerical solution for the streamfunction developed for any values of the heights of the rectangular and an L-shaped cavity. Streamline patterns near the corner are simulated for a different aspect ratio A. The techniques are tested on a flow problem undergoing Stokes or Navier-Stokes equations in a square cavity. It is seen that the method appears to be cheaper and more accurate than the numerical and analytical methods. It is expected that the study will lead to useful insights into the understanding of the flow topology near a re-entrant corner from a combined analytical-numerical method. Attention is then focused on the topological behavior near the re-entrant corner of the L-shaped cavity. Careful analysis of the streamlines of streamfunction near the re-entrant corner by using wall shear stress allows us to give a possible flow bifurcation of dividing streamline.