Higher-order streamfunction-vorticity formulation of 2D steady-state Navier-Stokes equations

Based on regular boundary element method and a kind of linearity invariance under homotopy, a kind of numerical scheme of 2D steady-state Navier-Stokes equation in streamfunction-vorticity formulation is described. The flow inside a square cavity is used to illustrate this numerical scheme.