“…A two dimensional, finite volume solver [56,57] was used to discretize the 2D flow domain through a first-order order upwind scheme. Several quadrilateral grids were generated for the purpose of ensuring that the solution is grid-independent, as detailed in Section 5.…”