ORIGINAL RESEARCH article

Front. Phys., 20 June 2022

Sec. Statistical and Computational Physics

Volume 10 - 2022 | https://doi.org/10.3389/fphy.2022.900502

A New Modified Analytical Approach for the Solution of Time-Fractional Convection–Diffusion Equations With Variable Coefficients

  • 1. Department of Mathematics, Abdul Wali Khan University Mardan, Mardan, Pakistan

  • 2. Department of Mathematics, Near East University TRNC, Mersin, Turkey

  • 3. Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan

  • 4. Department of Mathematics, Taiwan Theoretical and Computational Science (TaCS) Center, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), Bangkok, Thailand

  • 5. Department of Mathematics, College of Science, King Saud University, Riyadh, Saudi Arabia

  • 6. Department of Mathematics, Intelligent and Nonlinear Dynamic Innovations Research Center, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok (KMUTNB), Bangkok, Thailand

  • 7. University College Dublin, Dublin, Ireland

Abstract

In this article, a new modification of the Adomian decomposition method is performed for the solution fractional order convection–diffusion equation with variable coefficient and initial–boundary conditions. The solutions of the suggested problems are calculated for both fractional and integer orders of the problems. The series of solutions of the problems with variable coefficients have been provided for the first time. To verify and illustrate our new technique, four numerical examples are presented and solved by using the proposed technique. The derived results are plotted, and the dynamics are shown for both fractional and integer orders of the problems. An excellent variation among the solutions at various fractional orders is observed. It is analyzed that the new technique based on the Adomian decomposition method is accurate and effective. The present method fits both the initial and boundary conditions with double approximations simultaneously, which increases the accuracy of the present method. For the first time, the present technique is used for the solutions of the problems with variable coefficients along with initial and boundary conditions. It is therefore suggested to apply the present procedure for the solutions of other problems with variable order and coefficients along with initial and boundary conditions.

1 Introduction

Fractional Calculus (FC) is the branch of mathematics that focuses on the study of fractional-order operators and their applications in mathematical theory, which is the expansion of the basic classical concept in calculus and has attracted the attention of many researchers in the last few decades. FC has useful applications in numerous fields of physical and numerical sciences, such as relaxation processes, diffusion, visco-elasticity, hydrology, biomedical engineering, electro-chemistry, relaxation vibration, finance, probability, seismology, and so on []. Other FC applications include earthquake’s oscillations modeling [], statistical and mechanic continuum [], control theory [], engineering and physics [, ], entropy [], image processing [], field theories, optimal control, classical and quantum mechanics [, ], chaos theory [], fractional diabetes model [], and human diseases []. The main contribution to the development of (FC) includes Podlubny et al. []; Miller and Ross []; Oldham and Spanier []; and Kilbas and Trujillo [].

A well-known fact about fractional-order partial differential equations (FPDESs) is that they are involved in many natural and physical phenomena. Many researchers have attracted towards them due to the wonderful applications in various fields especially anomalous diffusion, dielectric polarization, relaxation vibration, scattering, system identification, fluid flow, mathematical physics [], psychology, and acoustics [].

Over the last 10 years, attempts were made to find robust analytical and numerical techniques to solve (FPDES). Some of them are the homotopy analysis method [], the wavelet operational method [, ], the Legendre-based method [], the fractional recti sub-equation method [], the generalized fractional Taylor series method [], the Laplace Adomian decomposition method [], the fractional wavelet method [], the Adomian decomposition method [], the discrete Adomian decomposition method [], the Mohand decomposition method [], the Wlash function method [], the Chebyshev spectral approximation [], the wavelet optimization method [], the solitary ansatz method [], the multiple exp-function method [], the Hirotas direct method [], the Laplace transform decomposition method [], the natural transform decomposition method [], and so on. The time-fractional convection–diffusion equations (CDESs) have been conventional used in the mathematical models of computations and simulations, such as the transport of energy and mass, oil reservoir simulations, flow of heat, particles, dispersion of chemicals in reactors, and global weather prediction [, ]. Time-fractional CDEs are operated in the abnormal diffusion processes related to time.

Several research work have been done on the Adomian decomposition method (ADM), such as Wazwaz A.M. [] used it to solve linear and non-linear operators, initial-value problems solved by Abdelrazec and Pelinovsky [] and Ray and Bera [] applied it to the differential equations. Hashim [] applied it to the boundary-value problems, fractional differential equations solved by Hu and Luo [], Duan J.S. applied to the boundary-value problems [], Song and Wang [] solved fractional differential equations, initial value problems solved by Wu [] and Almazmumy et al. []. Patel et al. applied the ADM sumudu transform method in [, ] and the fractional reduced differential transform method in [, ]. For the convergence of ADM, see [, ]. Fatoorehchi et al. solved explicit Frost–Kalkwarf type equations by ADM []. The one thing about all the above appreciable work is that ADM is used either with initial or only boundary conditions. Furthermore, a new modification is generated in ADM by Ali in [, ] for initial–boundary problems. He used the same new technique in [] with the variation iteration method, and the results are good. We extend his work by applying it to fractional order problems. This procedure gives accuracy because we construct a new initial approximate for each iteration. Abdelrazec et al. have provided a comprehensive study on the convergence of ADM for initial value problems in []. Similarly, Aminataei et al. have done the stability of ADM with other methods in [].

In the present work, time-fractional CDEs with variable coefficients are considered as follows:with the initial conditionsand the boundary conditionswhere c(ζ), d(ζ) ≠ 0 are continuous functions and 0 < γ < 1. The fractional derivative is in the Caputo sense.

Many techniques have been used to obtain the solutions of time-fractional CDEs such as the Gegenbauer spectral method [], the Sinc–Legendre collocation method [], the Bernstein polynomials method [], the flatlet oblique multi wavelets [], and the improved differential transform method []. Marjan et al. applied the radial basis interpolation method to compute the solution of this equation with constant coefficients []. Luo et al. computed the equation by the quadratic spline collocation method with constant coefficient and without convection []. In [], Pirkhedri and Javadi computed the equations with variable coefficient by the collocation approach.

In this article, we have investigated the analytical solutions of various time-fractional CDEs by using a new technique of the Adomian decomposition method (ADM). In the year 1980, Adomian had developed an effective technique, called ADM to solve differential equations involved in physical phenomena []. Later on, some short coming related to ADM, the researchers have modified ADM and construct some effective tools for the solutions of various fractional partial differential equations and their systems. Generally, many techniques are not efficient for the solutions of fractional initial and boundary conditions together. However, in this work, the solutions of initial and boundary value problems related to time-fractional CDEs are obtained by using a new technique of ADM. The obtained solutions are represented through graphs and tables. The 2D and 3D graphs have confirmed a closed contact between the exact and new version of ADM results. The fractional-order solutions represent the valuable dynamics of the suggested problems. The solution convergence of the fractional solutions towards integer-order solutions confirmed the validity of the suggested method. The present method is simple and required fewer calculations and therefore can be extended for the solutions of other higher non-linear fractional PDEs and their systems.

2 Preliminaries

In this section, a few definitions related to our work are taken into consideration.

2.1 Definition

The expression for the Caputo derivative of fractional order γ is as follows:where , ζ > 0, , and t ≥ 1.

2.2 Definition

The Reimann–Liouville integral operator having order γ is given by []where Γ denotes the gamma function and can be written as

2.3 Lemma

For p − 1 < γp with and with t ≥ −1, then [].where ζ > 0 and p − 1 < γp.

2.4 Definition

The Mittag-Leffler function Eγ(ρ) for γ > 0 is []

3 Adomian Decomposition Method

This method was discovered by Adomian for the solution of linear and nonlinear differential and integro-differential equations. To understand about the method, consider the following equation []:where g is the known function and F (u(x)) is the differential operator, which can be decomposed as follows:where L is the invertible operator of highest order derivative, R represents the linear operator, and N is the non-linear term. Then, Eq. (1) has the following representation:Apply the inverse operator L−1 on both sides of Eq. (2),where φ is the constant of integration and = 0.The ADM solution is represented as an infinite series,

The non-linear term Nu is denoted by An and is defined as follows:With the help of the following formula, we can calculate An asThe other method to calculate Adomian polynomials can be seen in []. The series has the following relation to represent the solution of Eq. (1):

4 Modification of Adomian Decomposition Method

To understand the main idea of the proposed technique, consider Eq. (1), Eq.(1) and Eq.(1). The new initial solution is calculated by new iteration for Eq. (1) with the help of the proposed technique [],Using ADM, the operator form of Eq. (1) iswhere the differential operator L is defined asHence, L−1 is defined as

Applying L−1 to Eq. (3),and using the new technique, the initial approximation becomesand the iteration formula becomeswhere n = 0, 1, ….It is obvious that initial solutions of Eq. (1) satisfied both the initial and boundary conditions as follows:The proposed technique works effectively for the two-dimensional problems.

5 Numerical Results

In this section, we will present the solution of some illustrative examples by using the new technique based on ADM.

5.1 Example 1

Consider TFCDE of the form []with the following initial and boundary conditions:

The problem has the exact solution at γ = 1,

Apply the new technique based on ADM to Eq. (4), we get the following result:where n = 0, 1, ….Applying L to Eq. (4), we havewhere and L−1 isOperating Eq. (6) by L−1, we haveUsing the ADM solution, the initial approximation becomesand using the new technique of initial approximation , the iteration formula becomes

By putting initial and boundary conditions in Eq.(5), for n = 0,

From Eq.(7), we have

For n = 1, Eq.(5) becomesFrom Eq. (7), we have

For n = 2, Eq. (5) becomesFrom Eq. (7), we have

Thus, the ADM solution as a series is

5.2 Example 2

Consider the TFCDE of the following form []:With the following initial and boundary conditions,whereThe problem has the exact solution at γ = 1 is,Applying the new technique based on ADM to Eq.(15), we get the following result:where n = 0, 1, …. Applying L to Eq.(15), we havewhere and L−1 isOperating by L−1, we haveUsing the ADM solution, the initial approximation becomesUsing the new technique of initial approximation , the iteration formula becomesBy putting the initial and boundary condition in Eq.(16), for n = 0,

From Eq. (9), we haveFor n = 1, Eq.(9), we haveFrom Eq. (9), we haveFor n = 2, Eq.(11) becomes

From Eq. (9), we have⋮Thus, the ADM solution in a series form is

5.3 Example 3

Consider TFCDE of the form []With the following initial and boundary conditions,whereThe problem has the exact solution at γ = 1 is,Applying the new technique based on ADM to Eq. (12), we get the following result:where n = 0, 1, …. Applying L to Eq. (12), we havewhere and L−1 isOperating Eq. (12) by L−1, we haveUsing the ADM solution, the initial approximation becomesand using the new technique of initial approximation , the iteration formula becomes

By putting the initial and boundary conditions in Eq. (13), for n = 0,From Eq. (15), we haveFor n = 1, Eq. (13), we haveFrom Eq. (15), we haveFor n = 2, Eq. (13) becomes

From Eq. (15), we have⋮Thus, the ADM solution as a series is

5.4 Example 4

Consider the TFCDE of the form []With the following initial and boundary conditions,where

The problem has the exact solution at γ = 1 is,

Applying the new technique based on ADM to Eq. (17), we get the following result:where n = 0, 1, ….

Applying L to Eq.(17), we havewhere and L−1 is

Operating Eq. (19) by L−1, we have

Using the ADM solution, the initial approximation becomesand using the new technique of initial approximation , the iteration formula becomes

By putting the initial and boundary condition in Eq.(18), for n = 0,

From Eq. (20), we have

For n = 1 Eq. (18), we have

From Eq. (20), we have

For n = 2, Eq. (18) becomes

From Eq.(20), we have

Thus, the ADM solution as a series is

6 Results and Discussion

Figures 1, 2 are the 2D and 3D representations of the Exact and approximate solutions of Example 1 while Figure 3 is the comparison plot of exact and approximate solutions. Figure 4 is the 2D and 3D plots for different fractional orders of Example 1. Figures 5, 6 represent 2D and 3D plots of exact and approximate solutions for Example 2 while Figure 7 is the comparison plot of the exact and approximate solutions. In Figure 8, the 2D and 3D plots at different fractional orders of the derivatives are discussed for Example 2. Figures 9, 10 have shown the 2D and 3D plots of exact and approximate solutions of Example 3 are presented while Figure 11 has displayed the comparison plot of exact and approximate solutions. Figure 12 is the 2D and 3D plots for different fractional orders of the derivative of example 3. Figures 13, 14 represent the 2D and 3D plots of exact and approximate solutions of Example 4 while Figure 15 is the comparison plot of the exact and approximate solutions. In Figure 16, the 2D and 3D plots are displayed for different fractional orders of the derivative of Example 4.

FIGURE 1

FIGURE 2

FIGURE 3

FIGURE 4

FIGURE 5

FIGURE 6

FIGURE 7

FIGURE 8

FIGURE 9

FIGURE 10

FIGURE 11

FIGURE 12

FIGURE 13

FIGURE 14

FIGURE 15

FIGURE 16

7 Conclusion

In many situations, the analytical solutions of fractional partial differential equations with boundary conditions are very difficult to investigate. Here, a very useful and successful attempt is made to determine the solution of fractional convection–diffusion equations with variable coefficients along with the initial–boundary conditions. The simulations for the given problems are derived and applied it for the solutions of some numerical examples. The fractional-order solutions are obtained in a very sophisticated manner. The solution graphs are plotted with higher accuracy for both fractional and integer orders of the problems. It is shown that the solution sequence of fractional order problems is convergent towards the integer-order solutions. Moreover, it is investigated that the present method has a unique capability to solve fractional partial differential equations with variable coefficients with initial and boundary conditions. In the light of the above varieties and features, the proposed technique can be extended to solve other problems with higher accuracy.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, and further inquiries can be directed to the corresponding author.

Author contributions

HK supervision, PK funding, Hajira methodology, QK Investigation, FT Project Administrator, KS Funding, ID draft writing.

Funding

The authors acknowledge the financial support provided by the Center of Excellence in Theoretical and Computational Science (TaCS-CoE), KMUTT. This research project is supported by Thailand Science Research and Innovation (TSRI) Basic Research Fund: Fiscal year 2022 under project number FRB650048/0164.

Acknowledgments

Researchers Supporting Project number (RSP-2021/401), King Saud University, Riyadh, Saudi Arabia.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

Adomian decomposition method, initial–boundary value problems, Caputo derivative, fractional convection–diffusion equations, analytical method

Citation

Khan H, Kumam P, Hajira, Khan Q, Tchier F, Sitthithakerngkiet K and Dassios I (2022) A New Modified Analytical Approach for the Solution of Time-Fractional Convection–Diffusion Equations With Variable Coefficients. Front. Phys. 10:900502. doi: 10.3389/fphy.2022.900502

Received

20 March 2022

Accepted

20 April 2022

Published

20 June 2022

Volume

10 - 2022

Edited by

Guillermo Fernandez-Anaya, Ibero American University, Mexico

Reviewed by

Ramakanta Meher, Sardar Vallabhbhai National Institute of Technology Surat, India

Trushitkumar Patel, University of the People, United States

Hooman Fatoorehchi, University of Tehran, Iran

Mustafa Turkyilmazoglu, Hacettepe University, Turkey

Updates

Copyright

*Correspondence: Poom Kumam,

This article was submitted to Statistical and Computational Physics, a section of the journal Frontiers in Physics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics