A numerical method using Laplace-like transform and variational theory for solving time-fractional nonlinear partial differential equations with proportional delay

Time-fractional nonlinear partial differential equations (TFNPDEs) with proportional delay are commonly used for modeling real-world phenomena like earthquake, volcanic eruption, and brain tumor dynamics. These problems are quite challenging, and the transcendental nature of the delay makes them even more difficult. Hence, the development of efficient numerical methods is open for research. In this paper, we use the concepts of Laplace-like transform and variational theory to develop a new numerical method for solving TFNPDEs with proportional delay. The stability and convergence of the method are analyzed in the Banach sense. The efficiency of the proposed method is demonstrated by solving some test problems. The numerical results show that the proposed method performs much better than some recently developed methods and enables us to obtain more accurate solutions.


Introduction
Fractional calculus is a valuable concept in applied mathematics dedicated to the study of integrals and derivatives of arbitrary positive order [1,2]. The fact of the matter is that fractional differential equations (FDEs), which are considered to be generalization of differential equations, enable us to describe functional values at continuous points not just at discrete points. For this reason, they are a popular choice for modeling complicated phenomena marked by memory and hereditary behaviors. Hence, fractional derivatives and integrals facilitate the analysis of real-world situations. It is also worth mentioning that nonlinear FDEs are preferable when it comes to dealing with various real-world events such as earthquake propagation, volcanic eruption, population growth, and the likes [3,4].
In the case when one is interested in incorporating the present time as well as history in the mathematical formulation of a certain physical situation, one needs to consider delay differential equations. These kinds of differential equations (DEs), which are also referred to as time-delay differential equations, are commonly used for modeling various phenomena in science and engineering. Combining the concepts of fractional derivatives and time-delay differential equations gives rise to fractional delay differential equations which are appropriate for capturing and modeling real-world events that involve memory and history. In this paper, we are interested in time-fractional nonlinear partial differential equations (TFNPDEs) with proportional delay. TFNPDEs are popular in different fields such as environmental sciences, economics, chemistry, physics, and life-sciences [5][6][7].
TFNPDE equation (1.1) can be used to represent models of real-world phenomena by selecting f wisely. For example, the time-fractional differential equation given by is known as Korteweg-de Vries equation arising in various studies related to shallow water waves. The Klein-Gordon equation which is a TFNPDE with proportional delay is given by where h(x, t) is a known analytic function, and F is a nonlinear function of u(x, t). This equation is usually used in the description of nonlinear wave interaction which arises in quantum field theory. Since TFNPDEs, and in general nonlinear differential equations, are difficult to solve and it is very hard to obtain exact solutions, many attempts have been made to develop numerical methods that provide good approximate solutions. For this reason, several techniques have been proposed to solve these kinds of DEs. Abazari and Ganji [8] proposed a two-dimensional differential transform method and its reduced form to solve nonlinear partial differential equations with proportional delay. Variational iteration method (VIM), initially proposed by He [9], has been proved to be a powerful technique for solving nonlinear DEs [10]. However, the success of this method mainly depends on accurate identification of Lagrange multipliers by using variational theory and integration by parts [9,11]. Moreover, application of VIM for obtaining series solutions of nonlinear problems requires repeated calculations and computations of a large number of unnecessary terms which reduce the efficiency of the method [12]. To avoid some of the shortcomings of VIM, Abassy et al. [12,13] proposed the modified variational iteration method (MVIM) and used it to find an approximate power series solutions of nonlinear problems. It is reported that MVIM facilitates and minimizes the computational work significantly thereby leading to fast convergence. Sakar and Saldir [14] also proposed an improved version of VIM based on incorporating an auxiliary parameter. The authors obtained the optimal value of this parameter by using residual error function and applied their method to solve various TFNPDEs. Saldir et al. [15] successfully applied reproducing kernel method (RKM) for solving time-fractional Kawahara equation with variable coefficients. The authors managed to construct orthonormal bases functions on specific Hilbert spaces that facilitate the solution process. Sakar et al. [16] developed an iterative method based on RKM to solve the time-fractional Cahn-Allen differential equation with Caputo derivative. Sakar et al. [17] designed a homotopy perturbation method (HPM) to solve various TFNPDEs with proportional delays in which the fractional derivative was taken in the Caputo sense. The authors compared their approach with traditional methods and reported that HPM is cost efficient. Singh and Kumar [18] combined HPM with Laplace transform and proposed a hybrid method called homotopy perturbation transform method (HPTM) to solve TFNPDEs with proportional delay. They reported that HPTM provides solutions in a series form and is a powerful technique for treating various physical phenomena in science and engineering. Singh and Kumar [19] developed an analytical method called alternative variational iteration method (AVIM) which is capable of proving approximate power series solutions to TFNPDE problems. The authors reported that AVIM competes favorably with DTM, HPM, and HPTM. Wang et al. [20] proposed two analytical methods called residual power series method (RPSM) and homotopy analysis transform method (HATM) to solve TFNPDEs with proportional delay. It is reported that the above methods provide approximate solutions which are in good agreement with HPM, HPTM, and exact solutions.
In this paper, we propose a new numerical method which benefits from the nice properties of Laplace-like transform and variational theory. We have observed that the new method facilitates computing generalized Lagrange multipliers which are crucial for solving nonlinear DEs. This paper is organized as follows: In Sect. 2, we discuss some basic concepts from fractional calculus and provide brief descriptions of a new Laplace-like transform and variational theory. Section 3 is dedicated to developing a new numerical method. In Sects. 3.1 and 3.2, we analyze the stability and convergence of the proposed method in the Banach sense. In Sect. 4, we present and analyze various numerical results. Moreover, we compare the efficiency of some methods which are used to solve TFNPDEs with proportional delay. Finally, in Sect. 5, we provide the conclusive remarks.

Preliminaries
In this section we present some basic concepts from fractional calculus which are crucial for understanding the remaining part of the paper. We also briefly describe Laplace-like transform and variational theory which will be used to design the new method.

Laplace-like transform
Laplace-like transform (LlT), also known as "Shehu transform", is a new integral transform derived to facilitate solving DEs [21]. LlT is defined over the set where the constant M is finite and k 1 , k 2 are either finite or infinite. Then LlT of the func- , is defined as follows: provided the integral on the right-hand side exists [21]. In particular, where is the gamma function. Moreover, LlT of the fractional order derivative C 0 D α t v(x, t) for α ≥ 0 is given by [22]

Variational theory
Modified variational iteration method (MVIM) is an efficient technique which has been implemented successfully to solve a large class of nonlinear problems [12,13,23]. Efficient implementation of MVIM rests on establishing the correction functional and identifying optimal Lagrange multipliers.
To illustrate the basic idea of MVIM, we consider the following TFNPDE problem: where C 0 D α t is the Caputo fractional derivative operator, R is a linear operator, Nv(x, t) is the nonlinear term, and g(x, t) is a given continuous function [11,23,24].
One of the most important steps of MVIM is to construct the correction functional for (2.2) which is defined as follows [12]: where λ is a general Lagrange multiplier, v n is the approximate solution at the nth iteration, andṽ n denotes a restricted variation, i.e., δṽ n = 0 with variational operator δ [25].
Making the correction functional in equation (2.3) stationary, i.e., δv n+1 (x, t) = 0, we obtain The other crucial step in MVIM is computing the Lagrange multiplier λ by solving the system of stationary conditions which can be derived from (2.4). In ordinary calculus, the optimal value of λ can be obtained by using variational theory and integration by parts. However, in fractional calculus, it is not possible to perform integration by parts on the second term of the above equation [24].

Banach's fixed point theorem
the image Tx coincides with x. [26]) Let X = (X, d) be a metric space. A mapping T : X → X is called a contraction on X if there is a nonnegative real number γ < 1 such that, for all x, y ∈ X,

Definition 2.5 (Contraction
(2.5) Suppose that X is complete, and let T : X → X be a contraction on X. Then T has a unique fixed point.
and posterior estimate

Theorem 2.2 ([27, 28])
Let (X, d) be a Banach space and T : X → X be a self-map of X satisfying

A new numerical method for solving TFNPDEs with proportional delay
Even if there are methods which approximate the solution of TFNPDE problems, for instance HPM [17], HPTM [18], RPSM, and HATM [20], the challenging nature of these DEs inspires the search for new and efficient methods. In the sequel, we derive a new method which makes use of Laplace-like (Shehu) transform and variational theory. We name this new method Laplace transform variational iteration method (LTVIM). Consider the general TFNPDE problem given in equation (1.1). Applying LlT as given in equation (2.1) to transform the fractional derivative in (1.1) gives (3.1) Applying MVIM as given in equation (2.3) to construct a correction functional of the transformed equation (3.1), we obtain , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) . Hence, the value of the Lagrange multiplier which satisfies the stationary condition 1 + λ( s u ) α = 0 is Inserting the Lagrange multiplier found above into equation (3.2) gives , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) . (3.4) Taking the inverse LlT of both sides of (3.4) gives , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) .
In this paper, we are interested in α ∈ (0, 1]. Therefore, m = 1 and k = 0 in (3.5) and the first term on the right-hand side of this equation becomes = v n (x, 0). , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) , (3.6) which is the (n + 1)th approximate solution of the TFNPDE problem in (1.1) for 0 < α ≤ 1.

Stability analysis
In the sequel, we state and prove an important result regarding the stability of LTVIM. To show Picard stability, it suffices to show that the mapping associated with LTVIM satisfies the conditions of Theorem 2.2. = v n (x, 0) , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) Proof For the sake of convenience, we drop (x, t) from v n (x, t) and v m (x, t). Thus, for n, m ∈ N, we have , q 1 t), . . . , ∂ n ∂x n v m (p n x, q n t) .  , q 1 t), . . . , ∂ n ∂x n v n (p n x, q n t) , q 1 t), . . . , ∂ n ∂x n v m (p n x, q n t) .
Using the linearity of S, S -1 and the partial differential operator f , we obtain , q 1 t) , . . . .

Convergence analysis and error estimate
To analyze the convergence of LTVIM and to estimate its error, we state and prove two important results. (3.12) By the triangle inequality and equation (3.12), for any m, n ∈ N such that n > m, we have Since 0 ≤ γ < 1, the sum 1 + γ + γ 2 + · · · + γ n-m-1 represents a finite geometric progression whose total sum is 1-γ n-m 1-γ . Therefore, Since 0 ≤ γ < 1 and v 1v 0 is fixed, the right-hand side of equation ( This implies that v -Tv = 0. Since · is a metric, we have Tv = v, i.e., v is the fixed point of T which is unique.

Corollary 3.1 The maximum absolute error of the approximate solution generated by LTVIM is estimated to be
Proof Taking n → ∞ in equation (3.13) gives The rest follows from Corollary 2.1.

Numerical results and discussion
In this section, we show the efficiency of LTVIM by solving some TFNPDE problems with proportional delay. We also compare the performance of LTVIM with recently proposed methods, namely HPM [17], HPTM [18], RPSM, and HATM [20].
The (n + 1)th approximate solution of (4.1) is obtained by applying the LTVIM iteration formula given in equation (3.6): First iteration Since v 0 (x, 0) = x, the first iteration solution is given by .

Third iteration
The third iteration solution, for v 2 (x, 0) = x, is given by
It can be seen from Fig. 1 that the approximate solution by LTVIM gets closer to the exact solution as α → 1. We have also compared some numerical methods on the basis of absolute error and recorded the results in Table 1. The intention is to compare LTVIM with HPM [17], HPTM [18], RPSM, and HATM [20]. The results in this table show that the solutions of the proposed method are in better agreement with the exact solution than those of the other methods. The maximum error obtained with LTVIM is 2.25387 × 10 -3 , but the maximum error obtained with the others is 7.46137 × 10 -3 , which is almost three times bigger than that of LTVIM. In general, the numerical results recorded in this   17,18,20]) Consider the following fractional partial differential equation with proportional delay: 3) x, t ∈ [0, 1] and 0 < α ≤ 1, with initial conditions v(x, 0) = x 2 . The results in Fig. 3 indicate that the approximate solution by LTVIM gets closer to the exact solution when the value of α approaches 1. In addition, it is obvious from the numerical results recorded in Table 2 that LTVIM provides more accurate solutions, which are very close to the exact solution, than HPM [17], HPTM [18], RPSM, and HATM [20]. One can also observe that the maximum error obtained with LTVIM is 4.91832×10 -4 , whereas the maximum error obtained with the others is 5.59603 × 10 -3 . This means LTVIM is at least ten times more accurate than other existing methods. Once again, the experimental results assure that LTVIM outperforms its counterparts and is well suited for solving nonlinear DEs of fractional order.

Conclusion
In this paper, we have studied TFNPDEs with proportional delay. A new numerical method called LTVIM is designed using the concepts of Laplace-like transform and variational theory. Very good conditions for the stability and convergence of LTVIM are constructed and analyzed in the Banach sense. Moreover, the efficiency of the new method is illustrated by solving some test problems. The numerical results show that LTVIM is very efficient and provides more accurate solutions than some recently developed methods. The promising experimental results signal that LTVIM could be applied successfully for other similar nonlinear problems.