A $C^0$-Weak Galerkin Finite Element Method for the Two-Dimensional Navier-Stokes Equations in Stream-Function Formulation. JCM. 2020;38(2):310-336. doi:10.4208/jcm.1806-m2017-0287