Abstract
The aim of the present work is to devote a friendly approach based on Adomian decomposition method (ADM) to find the numerical solution of the time-fractional Newell-Whitehead-Segel equation. Newell-Whitehead-Segel equation plays an efficient role in non-linear systems which describe the appearance of the stripe patterns in two dimensional systems. The numerical results obtained by proposed method are compared with exact solution for different values of fractional order α. Plotted graph illustrate the efficiency and accuracy of the proposed technique.
AMS Mathematics Subject Classification (2010): 44A99, 35Q99.
Introduction
Fractional calculus is a field of applied mathematics, three centuries old as the conventional calculus. Fractional calculus deals with derivatives and integrals of arbitrary orders. During the last decade, superb improvements have been visualized in the field of fractional calculus, very popular amongst science and engineering community. In recent year, differential equation containing fractional order derivatives has been contributed in various fields of science and engineering [1–4] such as diffusion equation, polarization, electro-magnetic waves, visco elasticity, electrode-electrolyte heat conduction, finance [5], control theory, biomedical engineering, biology [6] etc. In order to achieve the goal of highly accurate solution, many authors illustrate various techniques such as Adomian decomposition method [7], Finite difference method [8], Generalized differential transform method [9], Finite element method [10], Fractional differential transform method [11], Homotopy perturbation method [12, 13], Iterative methods [14], Variational iteration method [15], Homotopy analysis method [16], Differential quadrature method [17], Homotopy perturbation Sumudu transform method [18], Homotopy analysis transform method [19], Local fractional homotopy perturbation Sumudu transform method and Local fractional reduced differential transform method [20], Homotopy analysis Sumudu transform method [21] etc.
Recently various author used a new fractional derivative with Mittag-Leffler type kernel by different numerical method like Laplace decomposition method [22] and iterative method [23] etc.
The Newell-Whitehead-Segel equation model is the interaction of the effect of the diffusion term with the non-linear effect of the reaction term. Fractional Newell-Whitehead-Segel equation is written as
where a, b and k > 0 are real numbers and q is a positive integers. First term on the left hand side in Equation (1.1) represent the variation of u(x, t) with time at a fixed location, first term on the right hand side uxx represent the variation of u(x, t) with spatial variable at a specific time and term au − buq takes into account the effect of the source term. The function u(x, t) may be non-linear distribution of temperature in an infinitely thin and long rod or fluid flow as a velocity in an infinitely long pipe with narrow diameter.
Mostly two types of patterns are observed. First is the roll pattern in which cylinders form by fluid stream lines. These cylinders may be bend and form spiral like patterns. Second pattern is the hexagonal in which liquid flow is divided into honey comb cells. The same patterns, stripes and hexagons appear in different physical system. For example, stripes patterns are notice in human fingerprints, on zebra skin and in a visual cortex. Hexagonal patterns are obtained from the propagation of laser beams through a non-linear medium and in systems with chemical reaction and diffusion species [24].
Recently Newell-Whitehead-Segel equations were solved by S. S. Nourazar, M. Soori, and A. Nazari-Golshan by homotopy perturbation method [25], A. Prakash and M. Kumar [26] by He's variational iteration method. Also fractional model of Newell-Whitehead-Segel were solved by Kumar et al. [27] and Prakash et al. [28] by homotopy analysis Sumudu transform method and fractional variational iteration method, respectively. But fractional model of Newell-Whitehead-Segel has not been solved by Adomian decomposition method. Adomian decomposition method is very powerful and efficient numerical method for handling non-linear fractional model. Adomian decomposition method (ADM) demonstrates fast convergence of the solution and therefore provides several significant advantages. This method attacks directly on non-linear term, in a straightforward fashion without using linearization, discretization, perturbation or any other restrictive assumption. Many studies have shown that few terms of decomposition series provide numerical result of high degree of accuracy which makes the method powerful when compared with other existing numerical techniques.
The outline of this paper is as follow. First section is introductory, in the Basic Definition of Fractional Calculus the basic definition of fractional calculus is discussed, in Proposed Adomian Decomposition Method solution process of non-linear Newell-Whitehead-Segel equation by Adomian decomposition method is discussed, in Error Analysis of The Proposed Method error analysis of proposed technique is discussed, in Application of ADM to Fractional Newell-Whitehead-Segel Equation five test examples of fractional Newell-Whitehead-Segel equation are given to elucidate the proposed method ADM and in last Conclusion of the work is drawn.
Basic Definition of Fractional Calculus
In this section, we will introduce the basic definitions and properties of fractional calculus used to describe the proposed schemes.
Definition 2.1. A real function f(t), t > 0, is said to be in the space Cα, ϵ αR, if there exists a real number p, (p > α), such that where f1(t) ϵ C[0, ∞) and it is said to be in the space
Definition 2.2. The Riemann-Liouville fractional integral of order α ≥ 0, of a function f(t) ϵ Cβ, β ≥ −1 is defined as [29–31]:
For the Riemann-Liouville fractional integral, we have
where Γ. is the well-known Gamma Function.
Definition 2.3. The Caputo fractional derivative of is defined as [29–31]:
Where m − 1 < α ≤ m.
Proposed Adomian Decomposition Method
In this section, we illustrate the basic idea of the Adomian Decomposition method (ADM) for the time-fractional Newell-Whitehead-Segel equation.
Consider time-fractional Newell-Whitehead-Segel equation as
where a, b and k > 0 are real numbers and q is a positive integers with initial condition
Applying the operator on both sides of (3.1), we have
Next, we decompose the unknown function u(x, t) into sum of an infinite number of components given by the series
and the non-linear term can be decomposed as
where An are Adomian polynomial, given by
where n = 0, 1, 2, 3, …….
Components u0, u1, u2, u3, u4, …. are determined by substituting (3.3), (3.4), and (3.5) into (3.2) leading to
This can be written as
Adomian method uses the formal recursive relations as:
Error Analysis of the Proposed Method
Theorem 4.1. If we can find a constant 0 < ε < 1 such that ||um+1(x, t)|| ≤ ε||um(x, t)|| for each value of m. Moreover, if the truncated series is employed as a numerical solution u(x, t), then the maximum absolute truncated error is determined as
Proof. We have
Which proves the theorem.
Application of ADM to Fractional Newell-Whitehead-Segel Equation
In this section, five test examples of fractional Newell-Whitehead-Segel equation demonstrate the efficiency of proposed ADM.
Ex. 5.1. We study the linear time-fractional Newell-Whitehead-Segel equation
with initial condition
Applying the operator on both side of above defined problem, we have
This gives the following recursive relation:
Now, for the standard case when α = 1, we get u(x, t) = ex−t, which is the exact solution of the classical Newell-Whitehead-Segel equation as obtained by HPM [25] and VIM [26]. Here the numerical results obtained by ADM upto eight terms of approximation and exact solution as shown in Figures 1, 2 are almost identical. It can be observed that as the value of t increases, u decreases, and as x increases, u also increases. Hence, the accuracy of ADM can be enhanced by increasing the number of iterations.
Figure 1
Figure 2
Ex. 5.2. We study the non-linear time-fractional Newell-Whitehead-Segel equation
with initial condition
Applying the operator on both side of above defined problem, we have
This gives the following recursive relation:
In particular when α = 1, we get the solution in the form
Which converge to the exact solution of the classical Newell-Whitehead-Segel equation very fastly [25, 26].
Figure 3 shows the comparison of approximate solution for different value of fractional order α = 0.25, 0.50, 0.75, 1 and exact solution at α = 1, when η = 1. It is observed from the Figure 3 that there is a good agreement between exact solution and approximate solution at α = 1. It is also noticed that solution depends on the time-fractional derivative. Accuracy and efficiency can be enhanced by increasing the number of iterations.
Figure 3
Ex. 5.3. We study the non-linear time-fractional Newell-Whitehead-Segel equation.
With initial condition,
Applying the operator on both side of above equation, we get
This gives the following recursive relation:
In particular when α = 1, we get the solution in the form
Which converge to the exact solution of the classical Newell-Whitehead-Segel equation very fastly [25].
Figure 4 shows the comparison of third order approximate solution for different value of fractional order α = 0.25, 0.50, 0.75, 1 and exact solution at α = 1, for x = 1. It is observed from the Figure 4 that there is a good agreement between exact solution and approximate solution at α = 1. It is also noticed that solution depends on the time-fractional derivative. Accuracy and efficiency can be enhanced by increasing the number of iterations.
Figure 4
Ex. 5.4. We study the non-linear time-fractional Newell-Whitehead-Segel equation
with initial condition
Applying the operator on both side of above equation, we have
This gives the following recursive relation:
Taking α = 1, we get the solution in the form
Which converge to the exact solution of the classical Newell-Whitehead-Segel equation very fastly [25, 26].
Figure 5 shows the comparison of third order approximate solution for different value of fractional order α = 0.25, 0.50, 0.75, 1 and exact solution at α = 1 for x = 1. It is observed from the Figure 5 that there is a good agreement between exact solution and approximate solution at α = 1. It is also noticed that solution depends on the time-fractional derivative. Accuracy and efficiency can be enhanced by increasing the number of iterations.
Figure 5
Ex. 5.5. We study the nonlinear time-fractional Newell-Whitehead-Segel equation of the form
with initial condition
Applying the operator on both side of above defined problem, we have
This gives the following recursive relation:
Taking α = 1, we get the solution in the form
which converge to the exact solution of the classical Newell-Whitehead-Segel equation very fastly [25, 26].
Figure 6 shows the comparison of third order approximate solution for different value of fractional order α = 0.25, 0.50, 0.75, 1 and exact solution at α = 1, for x = 1. It is observed from the Figure 6 that there is a good agreement between exact solution and approximate solution at α = 1. It is also noticed that solution depends on the time-fractional derivative. Accuracy and efficiency can be enhanced by increasing the number of iterations.
Figure 6
Conclusion
In this article, we have successfully applied the ADM to obtain the approximate analytic solutions of fractional model of Newell-Whitehead-Segel equation. The plotted graph and numerical result shows the accuracy of proposed method. We observed an excellent agreement between ADM and the exact solution. The results reveal that ADM is an efficient and computationally very attractive approach to investigate non-linear fractional model. Therefore, ADM can be further applied to solve various types of linear and non-linear fractional model arising in the field of science and engineering.
Statements
Author contributions
AP and VV designed the study, collected the data, performed the analysis, and wrote the 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.
References
1.
OldhamKBSpanierJ. The Fractional Calculus: Theory and applications of Differentiation and Integration of Arbitrary Order. New York, NY: Academic Press (1974).
2.
PodlubnyI. Fractional Differential Equation. New York, NY: Academic Press (1999).
3.
JainS. Numerical analysis for the fractional diffusion and fractional Buckmaster's equation by two step Adam-Bashforth method. Eur Phy J Plus (2018) 133:19. 10.1140/epjp/i2018-11854-x
4.
MaginRL. Fractional Calculus in Bioengineering. New York, NY: Begell House Publishers (2006).
5.
RabertoMScalasMMainardiF. Waiting times and returns in high frequency financial data. An empirical study. Phys A (2002) 314:749–55. 10.1016/S0378-4371(02)01048-8
6.
MomaniSOdibatZ. Numerical approach to differential equations of fractional orders. J Comput Appl Math. (2007) 207:96–110. 10.1016/j.cam.2006.07.015
7.
RaySSBeraRK. Analytical solution of Bagley–Torvik equation by Adomian decomposition method. Appl Math Comput. (2005) 168:398–410. 10.1016/j.amc.2004.09.006
8.
MeerschaertMTadjeranC. Finite difference approximations for two sided space fractional partial differential equations. Appl Numer Math. (2006) 56:80–90. 10.1016/j.apnum.2005.02.008
9.
OdibatZMomaniSErturkVS. Generalized differential transform method: application to differential equations of fractional order. Appl Math Comput. (2008) 197:467–77. 10.1016/j.amc.2007.07.068
10.
JiangYMaJ. Higher order finite element methods for time-fractional partial differential equations. J Comput Appl Math. (2011) 235:3285–90. 10.1016/j.cam.2011.01.011
11.
ArikogluAOzkolI. Solution of a fractional differential equations by using differential transform method. Chaos Solitons Fract. (2007) 34:1473–81. 10.1016/j.chaos.2006.09.004
12.
ZhangX. Homotopy perturbation method for two dimensional time-fractional wave equation. Appl Math Model. (2014) 38:5545–52. 10.1016/j.apm.2014.04.018
13.
PrakashA. Analytical method for space-fractional telegraph equation by homotopy perturbation transform method. Nonlin Eng. (2016) 5:123–8. 10.1515/nleng-2016-0008
14.
DhaigudeCDNikamVR. Solution of fractional partial differential equations using iterative method. Fract Calc Appl Anal. (2012) 15:684–99. 10.2478/S13540-012-0046-8
15.
SafariMGanjiDDMoslemiM. Application of He's Variational iteration method and Adomain decomposition method to the fractional KdV-Burger-Kuramoto equation. Comput Appl. (2009) 58:2091–97. 10.1016/j.camwa.2009.03.043
16.
LiaoS. On the homotopy analysis method for nonlinear problem. Appl Math Comput. (2004) 147:499–513. 10.1016/S0096-3003(02)00790-7
17.
JiwariRPanditSMittalRC. Numerical solution of two dimensional Sine-Gordon Solitons by Differential Quadrature method. Comput Phys Commun. (2012) 183:600–16. 10.1016/j.cpc.2011.12.004
18.
SinghJKumarD. Homotopy perturbation Sumudu transform method for nonlinear equation. Adv Appl Mech. (2011) 14:165–75.
19.
KumarDSinghJBaleanuDRathoreS. Analysis of fractional model of Ambartsumian equation. Eur Phy J Plus (2018) 133:259. 10.1140/epjp/i2018-12081-3
20.
KumarDTichierFSinghJBaleanuD. An efficient computational technique for fractal vehicular traffic flow. Entropy (2018) 20:259. 10.3390/e20040259
21.
SinghJKumarDBaleanuDRathoreS. An efficient numerical algorithm for the fractional Drinfeld-Sokolov-Wilson equation. Appl Math Comput. (2019) 335:12–24. 10.1016/j.amc.2018.04.025
22.
KumarDSinghJBaleanuD. A new analysis of Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler type kernel. Eur Phy J Plus (2018) 133:70. 10.1140/epjp/i2018-11934-y
23.
KumarDSinghJBaleanuD. The analysis of chemical kinetics system pertaining to a fractional derivative with Mittag-Leffler type kernel. Chaos (2017) 27:103113. 10.1063/1.4995032
24.
GolovinbNA. General Aspect of pattern formation, pattern formation and growth phenomena in Nano-System. Alxaander (2007) 1–54.
25.
NourazarSSSooriM. On the exact solution of Newell-Whitehead-Segel equation using the Homotopy perturbation method. Aust J of Bas Appl Sci. (2011) 5:1400–11.
26.
PrakashAKumarM. He's variation iteration method for the solution of nonlinear Newell-Whitehead-Segel equation. J Appl Anal Comput. (2016) 5:123–8. 10.11948/2016048
27.
KumarDPrakashR. Numerical approximation of Newell-Whitehead-Segel equation of fractional order. Nonlin Eng. (2016) 5:81–6. 10.1515/nleng-2015-0032
28.
PrakashAGoyalMGuptaS. Fractional variational iteration method for solving time-fractional Newell-Whitehead-Segel equation. Nonlin Eng. (2018). [Epub ahead of print]. 10.1515/nleng-2018-0001
29.
HeJH. Homotopy perturbation technique. Appl Mech Eng. (1999) 178:257–62.
30.
HeJH. Homotopy perturbation method: a new nonlinear analytical technique. Appl Math Comput. (2003) 135:73–9. 10.1016/S0096-3003(01)00312-5
31.
HeJH. Homotopy perturbation method for bifurcation of nonlinear problem. J Nonlin Sci Numer Simul. (2005) 6:207–8. 10.1515/IJSNS.2005.6.2.207
Summary
Keywords
caputo fractional derivative, fractional newell-whitehead-segel equation, adomian decomposition method, fractional calculus, numerical method
Citation
Prakash A and Verma V (2019) Numerical Method for Fractional Model of Newell-Whitehead-Segel Equation. Front. Phys. 7:15. doi: 10.3389/fphy.2019.00015
Received
30 July 2018
Accepted
25 January 2019
Published
22 February 2019
Volume
7 - 2019
Edited by
Dumitru Baleanu, University of Craiova, Romania
Reviewed by
Mustafa Inc, Firat University, Turkey; Devendra Kumar, University of Rajasthan, India; Carlo Cattani, Università degli Studi della Tuscia, Italy
Updates
Copyright
© 2019 Prakash and Verma.
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.
*Correspondence: Amit Prakash amitmath@nitkkr.ac.in; amitmath0185@gmail.comVijay Verma vijay_mtech21@rediffmail.com
This article was submitted to Mathematical 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.