“…The computational domain consists of 28 Â 28 bi-quadratic elements representing 57 Â 57 grid points in cartesian coordinate system and grid generation with iso-parametric mapping has been discussed in earlier work [46,47]. The Galerkin finite element method has been used with bi-quadratic elements corresponding to lesser number of nodes, as that smoothly capture the non-linear variations of the field variables, which are in contrast with finite difference/finite volume solutions available in the literature [33]. Bi-quadratic elements are used for the finite element method and the detailed description of element basis functions can be found elsewhere [48,51].…”