A hybrid method to solve a fractional-order Newell–Whitehead–Segel equation

This paper solves fractional diﬀerential equations using the Shehu transform in combination with the q-homotopy analysis transform method (q-HATM). As the Shehu transform is only applicable to linear equations, q-HATM is an eﬃcient technique for approximating solutions to nonlinear diﬀerential equations. In nonlinear systems that explain the emergence of stripes in 2D systems, the Newell–Whitehead–Segel equation plays a signiﬁcant role. The ﬁndings indicate that the outcomes derived from the tables yield superior results compared to the existing LTDM in the literature. Maple is utilized to depict three-dimensional surfaces and ﬁnd numerical values that are displayed in a table.


Introduction
Leibnitz created arbitrary-order derivatives shortly after integer-order derivatives, and this topic has attracted the interest of numerous academics recently.As soon as its potential was known, fractional calculus (FC) quickly replaced classical calculus as the preferred modeling tool for practical issues.FC theory provides a superbly methodical framework for understanding the physical world.The specific explanation it offers for nonlinear complex systems has recently attracted the attention of numerous experts.Using fractional models of differential equations has the benefit of being nonlocal, which can be used in any model.Locality is preserved by derivatives of integer order, but nonlocality is introduced by those of fractional order.This demonstrates that the future state of a physical system is dependent on not only its current state but also its past states.Models built with the use of fractional-order derivatives are thus more accurate representations of reality.Classical derivatives and fractional derivatives are both mathematical concepts related to the rate of change of a function, but they differ in their definition and application.While classical derivatives deal with integer orders and provide information about instantaneous rates of change, fractional derivatives extend this concept to noninteger orders, offering a way to model systems with memory or long-range dependence.Fractional calculus is a field that has found applications in diverse scientific and engineering disciplines [1].
However, the fractional order varies based on time and space.The scenario leads to a rapidly expanding field of FPDEs with fractional operators of variable order.Several potent numerical approaches were established in the scientific literature, and numerous eminent scholars contributed to this topic.These techniques include the adomian decomposition method (ADM), the homotopy perturbation method (HPM), the homotopy analysis method (HAM), the collocation method, the Sumudu transform method (STM), the conformable Shehu homotopy perturbation method (CSHPM), the conformable q-Shehu homotopy analysis transform method (Cq-SHATM), the conformable sumudu decomposition method, the conformable Laplace decomposition method, and the differential transform method (DTM) .
In recent years, there have been studies on ST in the literature [33][34][35][36][37][38][39][40].Using the qhomotopy Shehu analysis transform method (q-HSATM), we analyze and find an approximated analytical solution for the time-fractional NWS equation in the current framework.The methodology under consideration has been tweaked, and it is now a sophisticated hybrid of q-HATM and the Shehu transform (ST).Since q-HATM is a modified method of HAM, it avoids the need for linearization, discretization, and perturbation, and it also requires less memory and fewer calculations than traditional HAM and no complicated polynomials, integrations, or physical parameters.Recently, many researchers have found it to be highly helpful in interpreting results for many types of nonlinear problems due to its consistency and effectiveness.Srivastava et al. determine the q-HATM solution for the model of the vibration equation and analyze its behavior [41].Veeresha et al. examine the approximation of an analytical solution for a smoking epidemic model of fractional order [42].Singh et al. suggest that the efficiency of q-HATM may be shown in the coupled system [43].Numerical simulations of the Gardner and Cahn-Hilliard equations are presented by Prakasha et al., who verify that the proposed method is more accurate than alternative classical techniques [44].Veeresha et al. provide an illustrative example of the novel numerical surfaces used in the mathematical model of cancer chemotherapy [45].The solution for the combined modified Boussinesq and approximate long wave was reported by Veeresha et al. [46].Using the energy-dependent Schrödinger potential as an example, Veeresha et al. demonstrated the superior performance of the suggested strategy over the coupled fractional reduced differential transform method (CFRDTM) for solving the associated coupled equation [47].The proposed method is also being considered by many researchers as a means to obtain an approximate analytical solution to a wide range of issues [48,49].
The fluctuations in the sand, the lines on seashells, and many other striped patterns like these can be modeled by amplitude equations.The NWSE is one of the most important amplitude equations in applied sciences.It shows how stripes appear in two-dimensional systems [50][51][52][53].The NWSE has the form [50,51,54] where m is a positive integer, a, h and c real numbers.The Shehu transform is used to solve linear equations such as the Laplace transform.However, to solve nonlinear equations, q-homotopy analysis must be used together with a numerical method such as the transform method.Therefore, in this study, the Shehu transform and the q-homotopy analysis transformation method were combined.This paper aims to present a new method, the q-SHATM and to use it to obtain new numerical solutions for a NWS equation.
In light of developments, this paper proposes to investigate the use of the q-HSATM to solve the time-fractional Newell-Whitehead-Segel equation numerically, and to compare the results with those obtained using a more recent technique.
The first aim of this study is to introduce the q-homotopy Shehu analysis transform method.The second aim is to solve a time-fractional Newell-Whitehead-Segel equation for the first time with the newly produced hybrid method.Several previously unmentioned solutions in the literature are retrieved and their graphic characteristics are plotted in their entirety.The main motivation of writing this paper is to obtain the new numerical solutions of a time-fractional Newell-Whitehead-Segel equation by using a new hybrid method, namely the q-homotopy Shehu analysis transform method (q-HSATM).Also, it has been observed that the results obtained from the table give better results than the LTDM [50] in the literature.
The structure of this paper is as follows.In Sect.2, we provide a quick recap of some of the preliminary definitions of the Caputo fractional derivative and other conclusions that can be helpful in the study of fractional differential equations.Considering the issue at hand, a potential q-HSATM solution procedure is presented in Sect.3. Convergence of the suggested approach is given in Sect. 4. The approximate solutions are found using q-HSATM on several numerical test situations presented in Sect. 5. Section 6 contains the discussion and results.Finally, the conclusion is given.

Preliminaries
Several fundamental definitions are provided in this section.

Definition 3 ([57])
The Mittag-Leffler function E a is defined as Definition 4 ( [33]) The Shehu transform (ST) of the function f (t) is given by Definition 5 ([33]) If V (s, u) is the ST of the function f (t), then the ST of CFD is defined by 3 The methodology of the q-homotopy Shehu analysis transform method In this section, the q-HSATM for nonlinear FPDEs is presented.In order to illustrate the technique for the suggested method, the nonlinear FPDEs are written in standard operator form with the initial condition where A is a linear operator, N is a nonlinear operator, g(x, t) is a source term, and D a t is a time-fractional derivative operator of order α.Now, by performing a Shehu transform on Eq. ( 7) and using the initial condition, it is acquired as The nonlinear operator by the assistance of HAM for a real function ϕ(x, t; q) is defined as where q [0, 1  n ].This establishes a homotopy as follows: where, h = 0 is an auxiliary parameter and S represents the Shehu transform.For q = 0 and q = 1 n , the results in Eq. ( 12) are, respectively, provided as ϕ(x, t; 0) = u 0 (x, t), ϕ x, t; Therefore, by amplifying q from 0 to 1 n , the solution ϕ(x, t; q) converges from u 0 (x, t) to the solution u(x, t).Employing the Taylor theorem around q and expanding ϕ(x, t; q), it is obtained as where Equation ( 14) converges at q = 1 n for the appropriate u 0 (x, t), n, and h.Then, we have one of the solutions of the original nonlinear equation of the form If we differentiate the zeroth-order deformation Eq. ( 13) m times with respect to q and we divide by m!, respectively, then for q = 0, it is obtained as where the vectors are defined by When the inverse Shehu transform to Eq. ( 17) is applied, then it is obtained as where and Here, H m is a homotopy polynomial and presented as and By utilizing q-HSATM, the series solution is defined by
Proof The solution of nonlinear FPDEs Eq. ( 7) is presented as where Assume that u and w are two different solutions of Eq. ( 7), then with the help of the aforementioned equation, we obtain Now, using the convolution theorem for the Shehu transform, it is obtained as Then, using the integral mean-value theorem, it yields As a result, we have (1τ )|w -u| ≤ 0. As τ ∈ (0, 1), |w -u| = 0, which implies that w = u.Hence, the solution is unique.
Theorem 2 (Convergence theorem [58]) Assume that X is a Banach space and F : X→X is a nonlinear mapping.If the inequality exists, then G has a fixed point in view of Banach fixed point theory [59].Furthermore, for the arbitrary choice of a 0 , b 0 ∈X, the sequence created by the q-HSATM converges to a fixed point of G and Proof Let us take a Banach space (C[J], .) of all continuous functions on J with the norm expressed as g(t) = max t∈J |g(t)|.Now, we demonstrate that the sequence {w n } is a Cauchy sequence in the Banach space: Now, utilizing the convolution theorem for the Shehu transform, it is obtained as Next, by the application of the integral mean-value theorem [59], it is obtained as Let m = n + 1, then we obtain By using the triangular inequality, we have Since τ ∈ (0, 1), 1τ m-n-1 < 1, then we obtain For w 1w 0 < ∞, so as m → ∞ then w mw n → 0. Thus, the sequence {w n } is a Cauchy sequence in C[J], and so the sequence is convergent.

The numerical solutions of the time-fractional Newell-Whitehead-Segel equation
Let us assune a = 1, h = 2, c = 3, and m = 2 for the nonlinear time-fractional Newell-Whitehead-Segel equation (NTFNWSE) [52,60]: with initial condition where β is an arbitrary constant.Now, by applying the Shehu transform to Eq. ( 39) and by using Eq. ( 40), then it is obtained as The nonlinear operator by using Eq. ( 41) is described by By applying the proposed algorithm, the mth-order deformation equation is defined by where On applying the inverse Shehu transform to Eq. ( 43), it is obtained as By the use of the initial condition, it is acquired as To find the values of u 1 (x, t), putting m = 1 in Eq. ( 45), it is obtained as Similarly, to find values of u 2 (x, t), putting m = 2 in Eq. ( 45), it is found as In this way, the other terms are generated.Thus, the q-HSATM solution of Eq. ( 39) has the form Substituting the values of α = 1, n = 1, h = -1 in Eq. ( 49), the obtained results Figure 1 demonstrates a 3D graph of the q-HSATM solution, the exact solution, and the absolute error in Eq. ( 49) for h Figure 2 shows a 2D graph of Eq. ( 49) for different values of α.
Table 1 shows the comparison of the absolute errors of the methods ATHPM, LTDM, and q-HSATM.

Result and discussion
For Eq. ( 39), the graphs of the temperature u(x, t) for different values of α = 0.9, α = 0.7, α = 0.6, and α = 1 are drawn.Figure 1 shows three-dimensional graphs of the q-HSATM solution, the exact solution, and the absolute error at h = -1, n = 1, α = 1, β = 1.For increasing values of t, the two-dimensional graphs of the numerical solutions to this equation for various values α and the value β = 1 are depicted in Fig. 2. Figure 2 shows that the Table 1 Comparative study between ATHPM [61], LTDM [50], and q-HSATM for the numerical solutions u(x, t) at α = 1 and β = 0.01 temperature u(x, t) rises as any values of the space variable x and the time variable t increases.It can be seen in Fig. 2 that as the α value approaches one, the temperature u(x, t) converges.It is shown in Table 1 that the absolute errors of the third-order q-HSATM solution were obtained.Table 1 indicates that the absolute error increases significantly when any values of the space variable x and the value of time t increases.Table 1 demonstrates that while the q-HSATM yields the same results as ATHPM, it yields far more robust results than LTDM.

Conclusion
This work investigates the performance of NTFNWSE using q-HSATM.It is imperative to illustrate the influence of the fractional operator incorporated into the model being examined.Furthermore, the MAPLE software has been utilized to construct 2D and 3D graphs that depict the solutions to this equation for different values of α.The variability of the general structure of surface graphs created by the Maple software for Eq. ( 39) is evident.In addition, MAPLE software was used to obtain the graphs of the numerical solutions of this equation for the various α and β = 1 values.For NTFNWSE, it is observed that the general structure of surface graphs plotted in Maple software differs.The numerical solutions for NTFNWSE have been quickly and successfully obtained.Therefore, it may be extrapolated that q-HSATM is overly effective and robust for obtaining numerical solutions for various fractional nonlinear partial differential equations.