Numerical integration of Navier–Stokes equations by time series expansion and stabilized FEM

This manuscript introduces an advanced numerical approach for the integration of incompressible Navier-Stokes (NS) equations using a Time Series Expansion (TSE) method within a Finite Element Method (FEM) framework. The technique is enhanced by a novel stabilization strategy, incorporating a Divergent Series Resummation (DSR) technique, which significantly augments the computational efficiency of the algorithm. The stabilization mechanism is meticulously designed to improve the stability and validity of computed series terms, enabling the application of the Factorial Series (FS) algorithm for series resummation. This approach is pivotal in addressing the challenges associated with the accurate and stable numerical solution of NS equations, which are critical in Computational Fluid Dynamics (CFD) applications. The manuscript elaborates on the variational formulation of Stokes problem and present convergence analysis of the method using the Ladyzhenskaya–Babuvska–Brezzi (LBB) condition. It is followed by the NS equations and the implementation details of the stabilization technique, underscored by numerical tests on laminar flow past a cylinder, showcasing the method’s efficacy and potential for broad applicability in fluid dynamics simulations. The results of the stabilization indicate a substantial enhancement in computational stability and accuracy, offering a promising avenue for future research in the field. • Convergence analysis of Stokes problem solution in the TSE-FEM framework. • Factorial Series algorithm for the incompressible Navier-Stokes equation in FEM. • Stabilization technique to improve space modes computation of the TSE. • The formulation is linked to the Navier Stokes alpha model. • The methodology is validated on the flow past a cylinder test case.