A numerical scheme to simulate the distributed-order time 2D Benjamin Bona Mahony Burgers equation with fractional-order space

In this study, a new class of the Benjamin Bona Mahony Burgers equation is introduced, which considers the distributedorder in the time variable and fractional-order space in the Caputo form in the 2D case. The 2D-modified orthonormal normalized shifted Ultraspherical polynomials are derived from 1Dmodified orthonormal normalized shifted Ultraspherical polynomials and 2D-modified orthonormal normalized shifted Ultraspherical polynomials and the orthonormal normalized shifted Ultraspherical polynomials are applied to approximate of the space and time variables, respectively. Moreover, the convergence analysis of these basis functions is investigated. Due to the time variable being in the distributed-order mode and the space variable being in the fractional-order case, to apply the desired numerical algorithm for this type of equation, operational matrices of ordinary, fractional and distributed-order derivatives are computed. In the proposed method, once the unknown function is approximated using the mentioned polynomial, the matrix form of the residual function is derived and then a system of algebraic equations is adopted by applying the collocation approach. An approximate solution is extracted for the original problem by solving constructed equation system. Several examples are examined to demonstrate the accuracy and capability of the method.