scholarly journals Finite Integration Method with Shifted Chebyshev Polynomials for Solving Time-Fractional Burgers’ Equations

Mathematics ◽  
2019 ◽  
Vol 7 (12) ◽  
pp. 1201
Author(s):  
Ampol Duangpan ◽  
Ratinan Boonklurb ◽  
Tawikan Treeyaprasert

The Burgers’ equation is one of the nonlinear partial differential equations that has been studied by many researchers, especially, in terms of the fractional derivatives. In this article, the numerical algorithms are invented to obtain the approximate solutions of time-fractional Burgers’ equations both in one and two dimensions as well as time-fractional coupled Burgers’ equations which their fractional derivatives are described in the Caputo sense. These proposed algorithms are constructed by applying the finite integration method combined with the shifted Chebyshev polynomials to deal the spatial discretizations and further using the forward difference quotient to handle the temporal discretizations. Moreover, numerical examples demonstrate the ability of the proposed method to produce the decent approximate solutions in terms of accuracy. The rate of convergence and computational cost for each example are also presented.

Symmetry ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 497 ◽  
Author(s):  
Ratinan Boonklurb ◽  
Ampol Duangpan ◽  
Phansphitcha Gugaew

In this article, the direct and inverse problems for the one-dimensional time-dependent Volterra integro-differential equation involving two integration terms of the unknown function (i.e., with respect to time and space) are considered. In order to acquire accurate numerical results, we apply the finite integration method based on shifted Chebyshev polynomials (FIM-SCP) to handle the spatial variable. These shifted Chebyshev polynomials are symmetric (either with respect to the point x = L 2 or the vertical line x = L 2 depending on their degree) over [ 0 , L ] , and their zeros in the interval are distributed symmetrically. We use these zeros to construct the main tool of FIM-SCP: the Chebyshev integration matrix. The forward difference quotient is used to deal with the temporal variable. Then, we obtain efficient numerical algorithms for solving both the direct and inverse problems. However, the ill-posedness of the inverse problem causes instability in the solution and, so, the Tikhonov regularization method is utilized to stabilize the solution. Furthermore, several direct and inverse numerical experiments are illustrated. Evidently, our proposed algorithms for both the direct and inverse problems give a highly accurate result with low computational cost, due to the small number of iterations and discretization.


Author(s):  
Waleed M. Abd-Elhameed ◽  
José A. Tenreiro Machado ◽  
Youssri H. Youssri

Abstract This paper presents an explicit formula that approximates the fractional derivatives of Chebyshev polynomials of the first-kind in the Caputo sense. The new expression is given in terms of a terminating hypergeometric function of the type 4 F 3(1). The integer derivatives of Chebyshev polynomials of the first-kind are deduced as a special case of the fractional ones. The formula will be applied for obtaining a spectral solution of a certain type of fractional delay differential equations with the aid of an explicit Chebyshev tau method. The shifted Chebyshev polynomials of the first-kind are selected as basis functions and the spectral tau method is employed for obtaining the desired approximate solutions. The convergence and error analysis are discussed. Numerical results are presented illustrating the efficiency and accuracy of the proposed algorithm.


2021 ◽  
Vol 5 (3) ◽  
pp. 103
Author(s):  
Ampol Duangpan ◽  
Ratinan Boonklurb ◽  
Matinee Juytai

In this paper, the finite integration method and the operational matrix of fractional integration are implemented based on the shifted Chebyshev polynomial. They are utilized to devise two numerical procedures for solving the systems of fractional and classical integro-differential equations. The fractional derivatives are described in the Caputo sense. The devised procedure can be successfully applied to solve the stiff system of ODEs. To demonstrate the efficiency, accuracy and numerical convergence order of these procedures, several experimental examples are given. As a consequence, the numerical computations illustrate that our presented procedures achieve significant improvement in terms of accuracy with less computational cost.


2019 ◽  
Vol 24 (2) ◽  
pp. 176-188 ◽  
Author(s):  
Eid H. Doha ◽  
Mohamed A. Abdelkawy ◽  
Ahmed Z.M. Amin ◽  
Dumitru Baleanu

In this manuscript, we introduce a spectral technique for approximating the variable-order fractional Riccati equation (VO-FRDEs). Firstly, the solution and its space fractional derivatives is expanded as shifted Chebyshev polynomials series. Then we determine the expansion coefficients by reducing the VO-FRDEs and its conditions to a system of algebraic equations. We show the accuracy and applicability of our numerical approach through four numerical examples.


2016 ◽  
Vol 12 (4) ◽  
pp. 6133-6138
Author(s):  
Khadijah Abu Alnaja

In this paper, the homotopy perturbation method is implemented to derive the explicit approximate solutions for the time-fractional coupled Burger's equations. The including fractional derivative is in the Caputo sense. Special attention is given to prove the convergence of the method. The results are compared with those obtained by the exact at special cases of the fractional derivatives. The results reveal that the proposed method is very effective and simple.


2017 ◽  
Vol 22 (3) ◽  
pp. 643-682 ◽  
Author(s):  
Jian Zhao ◽  
Huazhong Tang

AbstractThis paper develops Runge-Kutta PK-based central discontinuous Galerkin (CDG) methods with WENO limiter for the one- and two-dimensional special relativistic hydrodynamical (RHD) equations, K = 1,2,3. Different from the non-central DG methods, the Runge-Kutta CDG methods have to find two approximate solutions defined on mutually dual meshes. For each mesh, the CDG approximate solutions on its dual mesh are used to calculate the flux values in the cell and on the cell boundary so that the approximate solutions on mutually dual meshes are coupled with each other, and the use of numerical flux will be avoided. The WENO limiter is adaptively implemented via two steps: the “troubled” cells are first identified by using a modified TVB minmod function, and then the WENO technique is used to locally reconstruct new polynomials of degree (2K+1) replacing the CDG solutions inside the “troubled” cells by the cell average values of the CDG solutions in the neighboring cells as well as the original cell averages of the “troubled” cells. Because the WENO limiter is only employed for finite “troubled” cells, the computational cost can be as little as possible. The accuracy of the CDG without the numerical dissipation is analyzed and calculation of the flux integrals over the cells is also addressed. Several test problems in one and two dimensions are solved by using our Runge-Kutta CDG methods with WENO limiter. The computations demonstrate that our methods are stable, accurate, and robust in solving complex RHD problems.


Acta Numerica ◽  
2016 ◽  
Vol 25 ◽  
pp. 567-679 ◽  
Author(s):  
Ulrik S. Fjordholm ◽  
Siddhartha Mishra ◽  
Eitan Tadmor

A standard paradigm for the existence of solutions in fluid dynamics is based on the construction of sequences of approximate solutions or approximate minimizers. This approach faces serious obstacles, most notably in multi-dimensional problems, where the persistence of oscillations at ever finer scales prevents compactness. Indeed, these oscillations are an indication, consistent with recent theoretical results, of the possible lack of existence/uniqueness of solutions within the standard framework of integrable functions. It is in this context that Young measures – parametrized probability measures which can describe the limits of such oscillatory sequences – offer the more general paradigm of measure-valued solutions for these problems.We present viable numerical algorithms to compute approximate measure-valued solutions, based on the realization of approximate measures as laws of Monte Carlo sampled random fields. We prove convergence of these algorithms to measure-valued solutions for the equations of compressible and incompressible inviscid fluid dynamics, and present a large number of numerical experiments which provide convincing evidence for the viability of the new paradigm. We also discuss the use of these algorithms, and their extensions, in uncertainty quantification and contexts other than fluid dynamics, such as non-convex variational problems in materials science.


Author(s):  
B. V. Rathish Kumar ◽  
Gopal Priyadarshi

We describe a wavelet Galerkin method for numerical solutions of fourth-order linear and nonlinear partial differential equations (PDEs) in 2D and 3D based on the use of Daubechies compactly supported wavelets. Two-term connection coefficients have been used to compute higher-order derivatives accurately and economically. Localization and orthogonality properties of wavelets make the global matrix sparse. In particular, these properties reduce the computational cost significantly. Linear system of equations obtained from discretized equations have been solved using GMRES iterative solver. Quasi-linearization technique has been effectively used to handle nonlinear terms arising in nonlinear biharmonic equation. To reduce the computational cost of our method, we have proposed an efficient compression algorithm. Error and stability estimates have been derived. Accuracy of the proposed method is demonstrated through various examples.


2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Shaheed N. Huseen ◽  
Said R. Grace

A modifiedq-homotopy analysis method (mq-HAM) was proposed for solvingnth-order nonlinear differential equations. This method improves the convergence of the series solution in thenHAM which was proposed in (see Hassan and El-Tawil 2011, 2012). The proposed method provides an approximate solution by rewriting thenth-order nonlinear differential equation in the form ofnfirst-order differential equations. The solution of thesendifferential equations is obtained as a power series solution. This scheme is tested on two nonlinear exactly solvable differential equations. The results demonstrate the reliability and efficiency of the algorithm developed.


2014 ◽  
Vol 2014 ◽  
pp. 1-7 ◽  
Author(s):  
Li Chen ◽  
Yang Zhao ◽  
Hossein Jafari ◽  
J. A. Tenreiro Machado ◽  
Xiao-Jun Yang

The local fractional Poisson equations in two independent variables that appear in mathematical physics involving the local fractional derivatives are investigated in this paper. The approximate solutions with the nondifferentiable functions are obtained by using the local fractional variational iteration method.


Sign in / Sign up

Export Citation Format

Share Document