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

Saved in:
Bibliographic Details
Title: Fourier–Gegenbauer Integral Galerkin Method for Solving the Advection–Diffusion Equation with Periodic Boundary Conditions.
Authors: Elgindy, Kareem T.
Source: Computation; Sep2025, Vol. 13 Issue 9, p219, 37p
Subject Terms: ADVECTION-diffusion equations, GALERKIN methods, ITERATIVE methods (Mathematics), BOUNDARY value problems, MATHEMATICAL optimization, SIMULATION methods & models, INTEGRAL equations
Abstract: 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. [ABSTRACT FROM AUTHOR]
Copyright of Computation is the property of MDPI and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Complementary Index
Description
Abstract: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. [ABSTRACT FROM AUTHOR]
ISSN:20793197
DOI:10.3390/computation13090219