Hybrid Shifted Gegenbauer Integral–Pseudospectral Method for Solving Time-Fractional Benjamin–Bona–Mahony–Burgers Equation

This paper introduces a novel hybrid shifted Gegenbauer integral–pseudospectral (HSG-IPS) method to solve the time-fractional Benjamin–Bona–Mahony–Burgers (FBBMB) equation with high accuracy. The approach transforms the equation into a form with only a first-order derivative, which is approximated using a stable shifted Gegenbauer differentiation matrix (SGDM), while other terms are computed with precise quadrature rules. By integrating advanced techniques such as the shifted Gegenbauer pseudospectral method (SGPS), fractional derivative and integral approximations, and barycentric integration matrices, the HSG-IPS method achieves spectral accuracy. Numerical results show it reduces average absolute errors (AAEs) by up to 99.99% compared to methods like Crank–Nicolson linearized difference scheme (CNLDS) and finite integration method using Chebyshev polynomial (FIM-CBS), with computational times as low as 0.04–0.05 s. The method’s stability is improved by avoiding ill-conditioned high-order derivative approximations, and its efficiency is boosted by precomputed matrices and Kronecker product structures. Robust across various fractional orders, the HSG-IPS method offers a powerful tool for modeling wave propagation and nonlinear phenomena in fractional calculus applications.