Numerical study of Newtonian laminar Flow around circular and square cylinders

In this study, we present a novel algorithm designed to address the computational challenges posed by the governing equations of incompressible Newtonian fluid flow, particularly of bluff bodies of circular and square in a channel. This a significant innovation algorithm is developed based on the Taylor-Galerkin/Pressure-Correction finite element method. Notably, a pivotal element of this approach involves the implementation of the two-step Lax-Wendroff scheme to resolve the momentum equation. To assess the efficacy of our novel algorithm, we conducted a comprehensive analysis using a standard benchmark scenario; the initiation of Poiseuille flow through an axisymmetric rectangular channel encompassing bluff bodies. In this context, our investigation centered on the influence of the Reynolds number (Re) Re ) on various pertinent variables. The results obtained from our study underscore the remarkable congruence between the outcomes generated by the new algorithm and established physical principles, as well as findings from prior research, as evidenced through thorough comparisons with the existing body of literature. This research represents a significant stride toward advancing our understanding of fluid dynamics and computational methodologies in this domain.