Fourier–Gegenbauer Integral Galerkin Method for Solving the Advection–Diffusion Equation with Periodic Boundary Conditions

This study presents the Fourier–Gegenbauer integral Galerkin (FGIG) method, a new numerical framework that uniquely integrates Fourier series and Gegenbauer polynomials to solve the one-dimensional advection–diffusion (AD) equation with spatially symmetric periodic boundary conditions, achieving exponential convergence and reduced computational cost compared to traditional methods. The FGIG method uniquely combines Fourier series for spatial periodicity and Gegenbauer polynomials for temporal integration within a Galerkin framework, resulting in highly accurate numerical and semi-analytical solutions. Unlike traditional approaches, this method eliminates the need for time-stepping procedures by reformulating the problem as a system of integral equations, reducing error accumulation over long-time simulations and improving computational efficiency. Key contributions include exponential convergence rates for smooth solutions, robustness under oscillatory conditions, and an inherently parallelizable structure, enabling scalable computation for large-scale problems. Additionally, the method introduces a barycentric formulation of the shifted Gegenbauer–Gauss (SGG) quadrature to ensure high accuracy and stability for relatively low Péclet numbers. This approach simplifies calculations of integrals, making the method faster and more reliable for diverse problems. Numerical experiments presented validate the method’s superior performance over traditional techniques, such as finite difference, finite element, and spline-based methods, achieving near-machine precision with significantly fewer mesh points. These results demonstrate its potential for extending to higher-dimensional problems and diverse applications in computational mathematics and engineering. The method’s fusion of spectral precision and integral reformulation marks a significant advancement in numerical PDE solvers, offering a scalable, high-fidelity alternative to conventional time-stepping techniques.