“…The momentum equation, Equation (6), is often referred to as the stress-divergence form [28], which is the dominantly adopted form in FE-based flow solvers [5,26,27], due to the convenience of enforcing the natural boundary conditions in the weak form of Equation (6). Alternatively, the incompressible momentum equations can also be written as…”