Front. Phys.Frontiers in PhysicsFront. Phys.2296-424XFrontiers Media S.A.10.3389/fphy.2019.00193PhysicsOriginal ResearchFractional Approach for Equation Describing the Water Transport in Unsaturated Porous Media With Mittag-Leffler KernelPrakashaD. G.1VeereshaP.2SinghJagdev3*1Department of Mathematics, Faculty of Science, Davangere University, Davangere, India2Department of Mathematics, Karnatak University, Dharwad, India3Department of Mathematics, JECRC University, Jaipur, India
Edited by: Devendra Kumar, University of Rajasthan, India
Reviewed by: Haci Mehmet Baskonus, Harran University, Turkey; Francisco Gomez, Centro Nacional de Investigación y Desarrollo Tecnológico, Mexico
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
In this paper, we find the solution for a fractional Richards equation describing the water transport in unsaturated porous media using the q-homotopy analysis transform method (q-HATM). The proposed technique is to use graceful amalgamations of the Laplace transform technique with the q-homotopy analysis scheme as well as the fractional derivative that is defined with the Atangana-Baleanu (AB) operator. The fixed point hypothesis is considered in order to demonstrate the existence and uniqueness of the obtained solution for the proposed fractional order model. In order to validate and illustrate the efficiency of the future technique, we analyze the projected model in terms of fractional order. Meanwhile, the physical behavior of the q-HATM solutions are captured in terms of plots for diverse fractional order and the numerical simulation is also demonstrated. The achieved results illuminate that the future algorithm is easy to implement, highly methodical, effective, and very accurate in its analysis of the behavior of non-linear differential equations of fractional order that arise in the connected areas of science and engineering.
Laplace transformAtangana-Baleanu derivativeRichards equationq-homotopy analysis methodfixed point theoremIntroduction
Fractional calculus (FC) was originated in Newton's time, but, lately, it has fascinated and captured the attention of many scholars. For the last 30 years, the most intriguing leaps in scientific and engineering applications have been found within the framework of FC. The concept of the fractional derivative has been industrialized due to the complexities associated with a heterogeneous phenomenon. The fractional differential operators are capable of capturing the behavior of multifaceted media as they have diffusion processes. It has been a very essential tool, and many problems can be illustrated more conveniently and more accurately with differential equations having an arbitrary order. Due to the swift development of mathematical techniques that use computer software, many researchers started to work on generalized calculus to present their viewpoints while analyzing many complex phenomena.
Numerous pioneering directions are prescribed for the diverse definitions of fractional calculus by many senior researchers, and these have prearranged the foundation [1–6]. Calculus with fractional order is associated with practical ventures and is extensively employed within nanotechnology [7], optics [8], human diseases [9], chaos theory [10], and other areas [11–39]. The numerical as well as analytical solutions for these equations illustrate that these models have an important role in portraying the nature of non-linear problems within connected areas of science.
In order to illustrate the importance of the novel fractional order derivative and future scheme, we, in the present framework, consider the Richards equation, which plays a vital role in describing the nature of the porous medium as well as the penetration of unsaturated regions in the soil. In 1931, Lorenzo A. Richards was the first person to pioneer work on the unsaturated porous material in order to model water movement. Later, he derived an equation based on continuum mechanics, which govern the water flow in the soil [40]. In the proposed model for the momentum equation, the continuity equation is an amalgam with Darcy's law, and is defined in a one-dimensional form as follows, with soil water diffusivity symbolized by ρ and hydraulic conductivity by σ for unsaturated soil moisture content u
∂u∂t=∂∂z(ρ∂u∂z−σ),
where z designates the elevation above a vertical datum. Recently, many authors employed numerical as well as analytical techniques in order to analyze and predict the suitable models for parameters in the equation and solve the governing equation of unsaturated flow in soils. Meanwhile, three models are generally applied, namely (i) the exponential model, (ii) the van Genuchten model, and (iii) the Brook-Coreysmodel (BCM). Among these models, BCM is extensively applied due to its well-defined configuration and because it is associated with the largest pore size. The following equations describe the complete wet ability of the BC model [41, 42]:
σ(u)=σ0uk,ρ(u)=ρ0(n+1)un,
where σ0, k, ρ0, and n are constants denoting particle shape, pore-size distribution and many other soil properties. For n = 0 and k = 2, Equation (2) simplified it to the classic Burgers equation [43, 44], and some particular values signify the generalized Burgers equation, which is essential to describing the important physical phenomena. In the present study, we consider that BCM employed the RC equation. In this case, for the (n, 1) order, the RC equation coincides with the Burgers equation, and this is presented here [45, 46]:
ut+a(un)x+buxx=0,a,b≠1,n≥1.
The analytical solution for the above equation is presented:
u(x,t)=(c2a(1+tanh(c(n−1)2b(x−ct))))1n−1.
In the present scenario, many important and non-linear models are methodically and effectively analyzed with the help of fractional calculus. There have been diverse definitions that have been suggested by many senior research scholars like, Riemann, Liouville, Caputo, and Fabrizio. However, these definitions have their own limitations. The Riemann-Liouville derivative is unable to explain the importance of the initial conditions; the Caputo derivative has overcome this shortcoming but cannot explain the singular kernel of the phenomena. In 2015, Caputo and Fabrizio solved the above issues [47], and many researchers consult this derivative in order to analyze and find the solution for diverse classes of non-linear complex problems. Some issues, however, were pointed out in the CF derivative; non-singular kernel and non-local properties are very essential in describing the physical behavior and nature of the non-linear problems. In 2016, Atangana and Baleanu introduced and natured a novel fractional derivative, namely the AB derivative. This novel derivative was defined with the aid of Mittag-Leffler functions [48]. This fractional derivative buried all the above-cited issues and helps us to understand the natural phenomena in the systematic and effective way.
In this framework, we consider the fractional RC equation of the form
aABCDtαu(x,t)+a(un)x+buxx=0,
where α is fractional order of the system and defined with AB fractional operator, u is the water content with depth x. The fractional order is introduced in order to incorporate the memory effects and hereditary consequence in the system, and these properties aid us in capturing the essential physical properties of the complex problems.
Recently, many mathematicians and physicists have developed very effective and more accurate methods in order to find and analyze solutions for complex and non-linear problems that have arisen in science and technology. In connection with this is the homotopy analysis method (HAM) proposed by Chinese Mathematician Liao Shijun [49, 50]. HAM has been profitably and effectively applied to study the behavior of non-linear problems without perturbation or linearization. But, for computational work, HAM requires significant time and computer memory. To overcome this, there is a possibility of using an amalgamation of the considered method and well-known transformation techniques.
In the present investigation, we analyzed the nature of the q-homotopy analysis transform method (q-HATM) solution for the FCDG equation by applying q-HATM. The future algorithm is the combination of q-HAM with LT [51]. The method of the considered scheme is merging two strong methods to solve linear and non-linear fractional differential equations both analytically as well as numerically. The future technique has many sturdy properties, including a non-local effect, straight forward solution procedure, and a promising large convergence region; moreover, it is free from any assumptions, discretization, and perturbation. Recently, due to its reliability and efficacy, the considered method has been exceptionally applied by many researchers to understand physical behavior in diverse classes of complex problems [52–60]. The novelty of the future method is that it aids a modest algorithm to evaluate the solution, and it is natured by the homotopy and axillary parameters, which provide the rapid convergence of the obtained solution for a non-linear portion of the given problem. Meanwhile, it has prodigious generality because it plausibly contains the results obtained by many algorithms like q-HAM, HPM, ADM and some other traditional techniques. The considered method can preserve great accuracy while decreasing the computational time and work in comparison with other methods.
The considered non-linear model recently caught the attention of researchers from different areas of science. Since RC equation plays a significant role in portraying several complex phenomena, many authors have found and analyzed the solution using analytical as well as numerical schemes; for instance, authors in [61] considered analytical techniques and found solutions for the considered model with arbitrary surface boundary conditions, and authors in [62] presented the compression approximation and infiltration of the RC equation with an analytical solution, authors in [45] applied the Adomian decomposition scheme, and authors in [46] applied HAM in order to find the approximated analytical solution. In this paper, we made an attempt to find the solution for the FRC equation using q-HATM.
Preliminaries
Recently, many authors considered these derivatives to analyze a diverse class of models in comparison with classical order as well as other fractional derivatives, and they prove that the AB derivative is more effective while analyzing the nature and physical behavior of the models [63, 64]. Here, we define the basic notion of Atangana-Baleanu derivatives and integrals [48].
Definition 1. The fractional Atangana-Baleanu-Caputo derivative for a function f ∈ H1(a, b) (b > a, α ∈ [0, 1]) is presented:
aABCDtα(f(t))=ℬ[α]1−α∫atf′(ϑ)Eα[α(t−ϑ)αα−1]dϑ.
Definition 2. The AB derivative of fractional order for a function f ∈ H1(a, b), b > a, α ∈ [0, 1] in the Riemann-Liouville sense is presented:
aABRDtα(f(t))=ℬ[α]1−αddt∫atf(ϑ)Eα[α(t−ϑ)αα−1]dϑ.
Definition 3. The fractional AB integral related to the non-local kernel is defined by
Theorem 1. The following Lipschitz conditions, respectively, hold true for both Riemann-Liouville and AB derivatives defined in Equations (6) and (7) [48],
‖aABCDtαf1(t)−aABCDtαf2(t)‖<K1‖f2(x)−f2(x)‖,
and
‖aABCDtαf1(t)−aABCDtαf2(t)‖<K2‖f1(x)−f2(x)‖.
Theorem 2. The time-fractional differential equation aABCDtαf1(t)=s(t) has a unique solution, which is defined as [48]
f(t)=1−αℬ[α]s(t)+μℬ[α]Γ(α)∫ats(ς)(t−ς)α−1dς.Fundamental Idea of the Proposed Scheme
Here, we consider the arbitrary order differential equation in order to demonstrate the basic solution procedure [65, 66]
aABCDtαv(x,t)+ℛv(x,t)+Nv(x,t)=f(x,t),n−1<α≤n,
with the initial condition
v(x,0)=g(x),
where aABCDtαv(x,t) symbolize the AB derivative of v(x, t). On using the LT on Equation (13), we have after simplification
Here, φ(x, t; q) is the real valued function with respect to x, t and (q∈[0,1n]). Now, we define a homotopy as follows
(1−nq)ℒ[φ(x,t;q)−v0(x,t)]=ℏqN[φ(x,t;q)],
where L is signifies LT, q∈[0,1n](n≥1) is the embedding parameter and ℏ ≠ 0 is an auxiliary parameter. For q = 0 and q=1n, the results given below are hold true
φ(x,t;0)=v0(x,t),φ(x,t;1n)=v(x,t).
Now, by intensifying q from 0 to 1n, then φ(x, t; q) varies from v0(x, t) to v(x, t). By using the Taylor theorem near to q, we define φ(x, t; q) in series form and then we get
φ(x,t;q)=v0(x,t)+∑m=1∞vm(x,t)qm,
where
vm(x,t)=1m!∂mφ(x,t;q)∂qm|q=0.
The series (16) converges at q=1n for the proper chaise of v0(x, t), n and ℏ. Then
v(x,t)=v0(x,t)+∑m=1∞vm(x,t)(1n)m.
On m-times differentiating Equation (17) with q and lately dividing by m! and then substituting q = 0, we get
ℒ[vm(x,t)−kmvm−1(x,t)]=ℏℜm(v→m−1),
where the vectors are defined as
v→m={v0(x,t),v1(x,t),…,vm(x,t)}.
On employing the inverse LT on Equation (22), we have
Then, the terms of vm(x, t) we can obtain using the Equation (28). The q-HATM series solution is presented as
v(x,t)=∑m=0∞vm(x,t).Solution for FRC Equation
In order to present the solution procedure and efficiency of the future scheme, in this segment we consider the DSW equation of fractional order with two distinct cases. Further, by the help of obtained results we made an attempt to capture the behavior of q-HATM solution for different fractional order. By the help of Equation (5) for the function of cubic water content and constant, we have
aABCDtαu(x,t)+u2ux−uxx=0,0<α≤1,
with initial conditions
u(x,0)=u0(x,t).
Taking LT on Equation (29) and then using Equation (30), we get
L[u(x,t)]=1s(u0(x,t))+1ℬ[α](1−α+αsα)L{u2ux−uxx}.
The non-linear operator N is presented with the help of future algorithm as below
On applying inverse LT on Equation (34), it reduces to
um(x,t)=kmum−1(x,t)+ℏL−1{ℜm[u→m−1]}.
On simplifying the above equation systematically by using u0 (x, t), we can evaluate the terms of the series solution
u(x,t)=u0(x,t)+∑m=1∞um(x,t)(1n)m.Existence of Solutions for the Future Model
Here, we considered the fixed-point theorem in order to demonstrate the existence of the solution for the proposed model. Since the considered model cited in Equation (30) is non-local as well as complex, there are no particular algorithms or methods that exist to evaluate the exact solutions. However, under some particular conditions, the existence of the solution is assured. Now, Equation (30) is considered:
0ABCDtα[u(x,t)]=G(x,t,u).
The foregoing system is transformed to the Volterra integral equation using the Theorem 2 as follows
where a = ||u|| and b = ||u1|| (since u and u1 are the bounded functions). Putting η = δ (a2 + b2 + ab) − δ2 in the above inequality, then we have
‖g(x,t,u)−g(x,t,u1)‖≤η‖u(x,t)−u(x,t1)‖.
The Lipschitz condition is thus obtained for g. Further, we can see that if 0 ≤ (δ (a2 + b2 + ab)− δ2) < 1, then it implies the contraction. The recursive form of Equation (36) is defined as
We prove the following theorem by using the above result.
Theorem 4. The solution for the Equation (30) will exist, and if we have specific t0, then
(1−α)ℬ(α)η+αℬ(α)Γ(α)η<1.
Proof. Let us consider the bounded function u(x, t) satisfying the Lipschitz condition. Then, by Equation (43), we have
‖ϕi(x,t)‖≤‖un(x,0)‖[(1−α)ℬ(α)η+αℬ(α)Γ(α)η]n.
Therefore, the continuity as well as existence of the obtained solution is proved. Subsequently, in order to show the Equation (46) is a solution for the Equation (29), we consider
From the above condition, it is clear that u(x, t) = u*(x, t), if
(1−(1−α)ℬ(α)η−αℬ(α)Γ(α)ηt)≥0.
Hence, Equation (53) proves our essential result.
Theorem 5. Suppose un (x, t) and u (x, t) are defined in the Banach space (𝔅[0, T], ||·||). The series solution defined in Equation (29) converges to the solution of the Equation (13), if 0 < λ < 1.
Proof: Consider the sequence {Sn}, which is the partial sum of the Equation (29), and we have to prove {Sn} is the Cauchy sequence in (𝔅[0, T], ||·||). Now consider
But 0 < λ < 1, therefore ||Sn − Sm|| = 0. Hence, {Sn} is the Cauchy sequence. This proves the required result.
Numerical Results and Discussion
In the present investigation, we have found the solution for equation describing the water transport in unsaturated porous media using q-HATM with the help of Mittag-Leffler law. Here, we consider two distinct cases to present the effectiveness of the proposed method.
Case 1: In this case, we consider the conductivity term as a function of cubic water content and constant σ=u33cm/h and ρ = 1 cm2/h. At a=c=13,n=3and b = −1, Equation (4) becomes
ABCaDtau(x,t)+u2ux−uxx=0,0<α≤1,
with initial condition
u(x,0)=12(1+tanh(−x3)).
Case 2: In this segment, we consider the conductivity term as a function of quadric water content and constant σ=u44cm/h and ρ = 1 cm2/h. At a=c=11,n=1and b = −1, the Equation (4) becomes
aABCDtαu(x,t)+u3ux−uxx=0,0<α≤1,
with initial condition
u(x,0)=12(1+tanh(−3x8))3.
Here, we demonstrate the numerical simulation for the considered non-linear. In Tables 1 and 2, the error analysis has been validated. From the tables we can see that the proposed scheme is more accurate, and we confirm that the iterations increase the q-HATM solutions so that they get closer to the analytical solution.
Numerical simulation presented for u(x, t) of FR equation consider in Case 1 at n = 1, ℏ = −1 and α = 1.
x
t
|uExact−uq−HATM(3)|
|uExact−uq−HATM(4)|
|uExact−uq−HATM(5)|
2.5
0.25
6.50636 × 10−7
1.84782 × 10−8
1.01181 × 10−10
0.50
5.35446 × 10−6
2.97193 × 10−7
3.15969 × 10−9
0.75
1.85802 × 10−5
1.51192 × 10−6
2.33755 × 10−8
1
4.52584 × 10−5
4.80032 × 10−6
9.57906 × 10−8
5
0.25
3.86055 × 10−7
8.87727 × 10−10
5.42849 × 10−11
0.50
3.08044 × 10−6
1.50958 × 10−8
1.76075 × 10−9
0.75
1.03664 × 10−5
8.10613 × 10−8
1.35526 × 10−8
1
2.44931 × 10−5
2.71248 × 10−7
5.78869 × 10−8
7.5
0.25
2.68674 × 10−7
1.50440 × 10−9
2.57157 × 10−12
0.50
2.16147 × 10−6
2.41095 × 10−8
8.01731 × 10−11
0.75
7.33582 × 10−6
1.22240 × 10−7
5.92186 × 10−10
1
1.74857 × 10−5
3.86892 × 10−7
2.42301 × 10−9
10
0.25
1.26101 × 10−7
8.43873 × 10−10
4.16169 × 10−12
0.50
1.01562 × 10−6
1.35692 × 10−8
1.33829 × 10−10
0.75
3.45096 × 10−6
6.90365 × 10−8
1.01990 × 10−9
1
8.23569 × 10−6
2.19278 × 10−7
4.31168 × 10−9
Numerical simulation presented for u(x, t) of FR equation consider in Case 2 at n = 1, ℏ = −1 and α = 1.
x
t
|uExact−uq−HATM(3)|
|uExact−uq−HATM(4)|
|uExact−uq−HATM(5)|
2.5
0.25
3.32962 × 10−7
5.38964 × 10−9
1.86945 × 10−11
0.50
2.70710 × 10−6
8.65185 × 10−8
5.83456 × 10−10
0.75
9.28381 × 10−6
4.39361 × 10−7
4.31477 × 10−9
1
2.23573 × 10−5
1.39264 × 10−6
1.76784 × 10−8
5
0.25
9.75190 × 10−8
5.33186 × 10−10
1.16072 × 10−11
0.50
7.75698 × 10−7
8.72022 × 10−9
3.74948 × 10−10
0.75
2.60229 × 10−6
4.51222 × 10−8
2.87429 × 10−9
1
6.12959 × 10−6
1.45752 × 10−7
1.22274 × 10−8
7.5
0.25
9.03607 × 10−8
2.66049 × 10−10
2.09943 × 10−13
0.50
7.25010 × 10−7
4.25260 × 10−9
7.54358 × 10−12
0.75
2.45406 × 10−6
2.15062 × 10−8
6.08009 × 10−11
1
5.83395 × 10−6
6.78929 × 10−8
2.69468 × 10−10
10
0.25
5.17685 × 10−8
1.88526 × 10−10
2.13869 × 10−10
0.50
4.15714 × 10−7
3.07483 × 10−9
3.36348 × 10−9
0.75
1.40835 × 10−6
1.56930 × 10−8
1.69010 × 10−8
1
3.35098 × 10−6
4.98643 × 10−8
5.31486 × 10−8
The surfaces of the obtained solution and the exact solution in comparison with absolute error have been captured, respectively, in Figures 1 and 2 for Case 1 and Case 2. The behavior of the obtained solution for different orders is presented in Figure 3 for both the cases in terms of 2D plots. In order to analyze the variations of the obtained solution for the FRC equation cited in Case 1 and Case 2 with respect to the homotopy parameter (ℏ), and the (ℏ) curves are drawn for diverse μ and presented in Figure 4 with distinct n. In the plots, the horizontal line signifies the convergence region of the q-HATM solution and these curves aid us to adjust and handle the convergence province of the solution. For an appropriate value of ℏ, the achieved solution quickly tends to the exact solution. The small deviation in the physical behavior of the complex models stimulates the enormous new results to analyze and understand the nature in a better and systematic manner. Moreover, from all the plots we can see that the proposed method is more accurate and very effective in its analysis of the considered non-linear fractional order equations.
Surfaces of (A)uq−HATM, (B)uExact(C)uAbs. Err. = |uExact − uq−HATM| for FR equation considered in Case 1 at ℏ = −1, n = 1 and α = 1.
Surfaces of (A)uq−HATM, (B)uExact(C)uAbs. Err. = |uExact − uq−HATM| for FR equation considered in Case 2 at ℏ = −1, n = 1 and α = 1.
Nature of theq-HATM solution for (A) Case 1 and (B) Case 2 with distinct α at ℏ = −1, n = 1 and x = 1.
ℏ-curves for (A) Case 1 (B) Case 2 with distinct α at x = 1 and t = 0.01 with n = 1 and 2.
Since every non-linear differential equation does not have an exact solution we look for an approximated analytical solution thorugh which we can prove the exactness or accuracy of the proposed scheme, as opposed to an exact solution. As we mentioned earlier, the q-HATM is a modified algorithm of HAM, and it thus does not require perturbation, dissertation, linearization, or any assumptions. More importantly, the future method generalizes many traditional techniques, such as HAM, HPM, FRDTM, and others, because these are a special case of q-HATM (n = 1, ℏ = 1). In connection with this, we capture the physical behavior of q-HATM solution to illustrate the accuracy. Further, we noticed that the considered non-linear phenomenon is highly dependent on a fractional operator. In order to illustrate the computational level and computational cost, the numerical simulation has been presented. From the table, it shows that as a number of series terms increases the solution converges to an analytical solution.
Conclusion
In this paper, the q-HATM is applied profitably to find the solution for an arbitrary order RC equation describing the water transport in the unsaturated porous media. Since AB derivatives and integrals having fractional order are defined with the help of generalized Mittag-Leffler function as the non-local kernel and non-singular, the present investigation illuminates the effectiveness of the considered derivative. The existence and uniqueness of the obtained solution is demonstrated by the fixed point hypothesis. The results obtained by the future scheme are more stimulating as compared to results available in the literature. Further, the proposed algorithm finds the solution of the non-linear problem without considering any discretization, perturbation or transformations.
The behavior of the obtained series solution has been captured in terms of 2D and 3D plots for distinct fractional order. These plots show that the q-HATM solution is more accurate and also conformed with the help of numerical simulation, and this is cited in the tables. Further, we confirm that, as the order of the solution increases, the obtained solutions converge to the exact solution. The present investigation illuminates how the considered complex non-linear phenomena noticeably depend on the time history and the time instant, which can be proficiently analyzed by applying the concept of calculus to fractional order. The present investigation helps the researchers to study the behavior of non-linear problems, and this gives very interesting and useful consequences. The proposed derivative provides non-singular kernel and non-local properties; these properties are very essential in describing the physical behavior and nature of the non-linear problems, and hence researchers can consider the AB derivative to solve many non-linear complex problems. Lastly, we can conclude the projected method is extremely methodical, effective and very accurate, and that it can be applied to the analysis of the diverse classes of non-linear problems that exist in science and technology.
Data Availability Statement
All datasets generated for this study are included in the article/supplementary material.
Author Contributions
All the authors have worked equally in this manuscript. All the authors have read and approved the final manuscript.
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.
ReferencesLiouvilleJ. Memoire surquelques questions de geometrieet de mecanique, etsur un nouveau genre de calcul pour resoudreces questions. (1832) 13:1–69.RiemannGFB. . Leipzig: Gesammelte Mathematische Werke (1896).CaputoM. . Bologna: Zanichelli (1969).MillerKSRossB. . New York, NY: Wiley (1993).PodlubnyI. . New York, NY: Academic Press (1999).KilbasAASrivastavaHMTrujilloJJ. . Amsterdam: Elsevier (2006).BaleanuDGuvencZBTenreiro MachadoJA. . London; New York, NY: Springer Dordrecht Heidelberg (2010). 10.1007/978-90-481-3293-5EsenSSulaimanTABulutHBaskonusHM. Optical solitons and other solutions to the conformable space-time fractional Fokas-Lenells equation. . (2018) 167:150–6. 10.1016/j.ijleo.2018.04.015VeereshaPPrakashaDGBaskonusHM. Solving smoking epidemic model of fractional order using a modified homotopy analysis transform method. (2019) 13:115–28. 10.1007/s40096-019-0284-6BaleanuDWuGCZengSD. Chaos analysis and asymptotic stability of generalized Caputo fractional differential equations. (2017) 102:99–105. 10.1016/j.chaos.2017.02.007VeereshaPPrakashaDGBaskonusHM. New numerical surfaces to the mathematical model of cancer chemotherapy effect in Caputo fractional derivatives Chaos. (2019) 29:013119. 10.1063/1.5074099BaskonusHMSulaimanTABulutH. On the new wave behavior to the Klein-Gordon-Zakharov equations in plasma physics. (2019) 93:393–9. 10.1007/s12648-018-1262-9PrakashaDGVeereshaPBaskonusHM. Analysis of the dynamics of hepatitis E virus using the Atangana-Baleanu fractional derivative, . (2019) 134:1–11. 10.1140/epjp/i2019-12590-5VeereshaPPrakashaDGBaleanuD. An efficient numerical technique for the nonlinear fractional Kolmogorov-Petrovskii-Piskunov equation. . (2019) 7:1–18. 10.3390/math7030265YokuşAGülbaharS. Numerical solutions with linearization techniques of the fractional Harry Dym equation. . (2019) 4:35–42. 10.2478/AMNS.2019.1.00004CattaniC. Haar wavelet-based technique for sharp jumps classification. . (2004) 39:255–78. 10.1016/S0895-7177(04)90010-6YangA-MZhangY-ZCattaniCXieG-NRashidiMMZhouY-J. Application of local fractional series expansion method to solve Klein-Gordon equations on Cantor sets. . (2014) 2014:372741. 10.1155/2014/372741ZhangYCattaniCYangX-J. Local fractional homotopy perturbation method for solving non-homogeneous heat conduction equations in fractal domains. . (2015) 17:6753–64. 10.3390/e17106753PrakashaDGVeereshaPBaskonusHM. Residual power series method for fractional Swift–Hohenberg equation. . (2019) 3:1–9. 10.3390/fractalfract3010009SinghJKumarDHammouchZAtanganaA. A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. . (2018) 316:504–15. 10.1016/j.amc.2017.08.048AsifNAHammouchZRiazMBBulutH. Analytical solution of a Maxwell fluid with slip effects in view of the Caputo-Fabrizio derivative. . (2018) 133:272. 10.1140/epjp/i2018-12098-6VeereshaPPrakashaDGMageshNNandeppanavarMMChristopherAJ. Numerical simulation for fractional Jaulent–Miodek equation associated with energy-dependent Schrödinger potential using two novel techniques. . (2019). 10.1080/17455030.2019.1651461. [Epub ahead of print].BaskonusHMMekkaouiTHammouchZBulutH. Active control of a chaotic fractional order economic system. . (2015) 17:5771–83. 10.3390/e17085771RavichandranCJothimaniKBaskonusHMValliammalN. New results on nondensely characterized integrodifferential equations with fractional order, . (2018) 133:1–10. 10.1140/epjp/i2018-11966-3GaoWGhanbariBBaskonusHM. New numerical simulations for some real world problems with Atangana-Baleanu fractional derivative. . (2019) 128:34–43. 10.1016/j.chaos.2019.07.037AtanganaAKocaI. Chaos in a simple nonlinear system with Atangana–Baleanu derivatives with fractional order. . (2016) 89:447–54. 10.1016/j.chaos.2016.02.012AbroKAMirbharMNGómez-AguilarJF. Functional application of Fourier sine transform in radiating gas flow with non-singular and non-local kernel. . (2019) 41:400. 10.1007/s40430-019-1899-0AbroKAKhanIGómez-AguilarJF. Thermal effects of magnetohydrodynamic micropolar fluid embedded in porous medium with Fourier sine transform technique. . (2019) 41:174. 10.1007/s40430-019-1671-5AbroKAGómez-AguilarJF. A comparison of heat and mass transfer on a Walter's B fluid via Caputo-Fabrizio versus Atangana-Baleanu fractional derivatives using the Fox-H function. . (2019) 134:101. 10.1140/epjp/i2019-12507-4Gómez-Aguilar1JFYépez-MartínezHTorres-JiménezJCórdova-FragaTEscobar-JiméneRFOlivares-PeregrinoVH. Homotopy perturbation transform method for nonlinear differential equations involving to fractional operator with exponential kernel. . (2017) 2017:68. 10.1186/s13662-017-1120-7Yepez-MartinezaHGomez-AguilarJFSosaaIOReyesaJMTorres-JimenezJ. The Feng's first integral method applied to the nonlinear mKdV space-time fractional partial differential equation. (2016) 62:310–6. 10.1155/2016/7047126Gómez-AguilarJFTorresLYépez-MartínezHBaleanuDReyesJMSosaIO. Fractional Liénard type model of a pipeline within the fractional derivative without singular kernel. . (2016) 2016:173. 10.1186/s13662-016-0908-1Morales-DelgadoVFTaneco-HernándezMAGómez-AguilarJF. On the solutions of fractional order of evolution equations. . (2017) 132:47. 10.1140/epjp/i2017-11341-0SinghJKumarDBaleanuD. New aspects of fractional Biswas-Milovic model with Mittag-Leffler law. . (2019) 14:1–23. 10.1051/mmnp/2018068BhatterSMathurAKumarDSinghJ. A new analysis of fractional Drinfeld–Sokolov–Wilson model with exponential memory. . (2020) 537:122578. 10.1016/j.physa.2019.122578KumarDSinghJBaleanuD. On the analysis of vibration equation involving a fractional derivative with Mittag-Leffler law. . (2019). 10.1002/mma.5903. [Epub ahead of print].GoswamiASinghJKumarDSushila. An efficient analytical approach for fractional equal width equations describing hydro-magnetic waves in cold plasma. . (2020) 524:563–75. 10.1016/j.physa.2019.04.058KumarDSinghJQurashiMABaleanuD. A new fractional SIRS-SI malaria disease model with application of vaccines, anti-malarial drugs, and spraying. . (2019) 2019:278. 10.1186/s13662-019-2199-9VeereshaPPrakashaDGKumarD. An efficient technique for nonlinear time-fractional Klein–Fock–Gordon equation. . (2020) 364:124637. 10.1016/j.amc.2019.124637RichardsLA. Capillary conduction of liquids through porous medium. . (1931) 1:318–33. 10.1063/1.1745010BrooksRHCoreyAT. . Hydrol Paper 3. Fort Collins, CO: Colorado State University (1964).CoreyAT. . Highlands Ranch, CO: Water Resource Publications (1994).WhitamGB. . New York, NY: Wiley. (1974).BashaHA. Burgers' equation: a general nonlinear solution of infiltration and redistribution. . (2002) 38:1–9. 10.1029/2001WR000954NasseriMShaghaghianMRDaneshbodYSeyyedianH. An analytic solution of water transport in unsaturated porous media. (2008) 11:591–601. 10.1615/JPorMedia.v11.i6.60JafariHFiroozjaeeMA. Application of homotopy analysis method for water transport in unsaturated porous media. (2010) 1:8–13.CaputoM. Fabrizio M. A new definition of fractional derivative without singular kernel. (2015) 1:73–85. 10.12785/pfda/010201AtanganaABaleanuD. New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model. (2016) 20:763–9. 10.2298/TSCI160111018ALiaoSJ. Homotopy analysis method and its applications in mathematics. . (1997) 5:111–25.LiaoSJ. Homotopy analysis method: a new analytic method for nonlinear problems. (1998) 19:957–62. 10.1007/BF02457955SinghJKumarDSwroopR. Numerical solution of time- and space-fractional coupled Burgers' equations via homotopy algorithm. (2016) 55:1753–63. 10.1016/j.aej.2016.03.028SrivastavaHMKumarDSinghJ. An efficient analytical technique for fractional model of vibration equation. (2017) 45:192–204. 10.1016/j.apm.2016.12.008PrakashaDGVeereshaPBaskonusHM. Two novel computational techniques for fractional Gardner and Cahn-Hilliard equations, . (2019) 1:1–19. 10.1002/cmm4.1021BulutHKumarDSinghJSwroopRBaskonusHM. Analytic study for a fractional model of HIV infection of CD4+T lymphocyte cells. . (2018) 2:33–43. 10.22436/mns.02.01.04VeereshaPPrakashaDG. Solution for fractional Zakharov-Kuznetsov equations by using two reliable techniques. . (2019) 60:313–30. 10.1016/j.cjph.2019.05.009KumarDAgarwalRPSinghJ. A modified numerical scheme and convergence analysis for fractional model of Lienard's equation. . (2018) 399:405–13. 10.1016/j.cam.2017.03.011PrakashAPrakashaDGVeereshaP. A reliable algorithm for time-fractional Navier-Stokes equations via Laplace transform. . (2019). 8:695–701. 10.1515/nleng-2018-0080VeereshaPPrakashaDGBaskonusHM. Novel simulations to the time-fractional Fisher's equation. . (2019) 13:33–42. 10.1007/s40096-019-0276-6PrakashAVeereshaPPrakashaDGGoyalM. A homotopy technique for fractional order multi-dimensional telegraph equation via Laplace transform. . (2019) 134:19. 10.1140/epjp/i2019-12411-ySinghJKumarDBaleanuDRathoreS. An efficient numerical algorithm for the fractional Drinfeld–Sokolov–Wilson equation. . (2018) 335:12–24. 10.1016/j.amc.2018.04.025ParlangJ-YBarryDAParlangeMBHogarthWLHaverkampRRossPJ. New approximate analytical technique to solve Richards equation for arbitrary surface boundary conditions. . (1997) 33:903–6. 10.1029/96WR03846ParlangJ-YHogarthWLBarryDAParlangeMBHaverkampRRossPJ. Analytical approximation to the solutions of Richards' equation with applications to infiltration, ponding, and time compression approximation. . (1999) 23:189–94. 10.1016/S0309-1708(99)00022-6AtanganaAAlkahtaniBT. Analysis of the Keller-Segel model with a fractional derivative without singular kernel. (2015) 17:4439–53. 10.3390/e17064439AlkahtaniBTAtanganaA. Analysis of non-homogenous heat model with new trend of derivative with fractional order. (2016) 89:566–71. 10.1016/j.chaos.2016.03.027SinghJKumarDSwroopRKumarS. An Efficient computational approach for time-fractional Rosenau-Hyman equation. (2018) 30:3063–70. 10.1007/s00521-017-2909-8VeereshaPPrakashaDG. A novel technique for (2+1)-dimensional time-fractional coupled Burgers equations. . (2019) 166:324–45. 10.1016/j.matcom.2019.06.005