New Analytical Technique for Solving a System of Nonlinear Fractional Partial Differential Equations

This paper introduces a new analytical technique (NAT) for solving a system of nonlinear fractional partial differential equations (NFPDEs) in full general set. Moreover, the convergence and error analysis of the proposed technique is shown. The approximate solutions for a system of NFPDEs are easily obtained by means of Caputo fractional partial derivatives based on the properties of fractional calculus. However, analytical and numerical traveling wave solutions for some systems of nonlinear wave equations are successfully obtained to confirm the accuracy and efficiency of the proposed technique. Several numerical results are presented in the format of tables and graphs to make a comparison with results previously obtained by other well-known methods


Introduction
Over the last few decades, fractional partial differential equations (FPDEs) have been proposed and investigated in many research fields, such as fluid mechanics, the mechanics of materials, biology, plasma physics, finance, and chemistry, and they have played an important role in modeling the so-called anomalous transport phenomena as well as in theory of complex systems, see [1][2][3][4][5][6][7][8].In study of FPDEs, one should note that finding an analytical or approximate solution is a challenging problem, therefore, accurate methods for finding the solutions of FPDEs are still under investigation.Several analytical and numerical methods for solving FPDEs exist in the literature, for example; the fractional complex transformation [9], homotopy perturbation method [10], a homotopy perturbation technique [11], variational iteration method [12], decomposition method [12], and so on.There are, however, a few solution methods for only traveling wave solutions, for example; the transformed rational function method [13], the multiple exp-function algorithm [14]), and some references cited therein.
The system of NFPDEs have been increasingly used to represent physical and control systems (see for instant, [15][16][17] and references cited therein).The systems of nonlinear wave equations play an important role in a variety of oceanographic phenomena, for example, in the change in mean sea level due to storm waves, the interaction of waves with steady currents, and the steepening of short gravity waves on the crests of longer waves (see for example, [18][19][20][21][22]).In this paper, two systems of nonlinear wave equations with a fractional order are studied; one is the nonlinear KdV system (see [23,24]) and another one is the system of dispersive long wave equations (see [24][25][26]).Some numerical or analytical methods have been investigated for solving a system of NFPDEs, such as an iterative Laplace transform method [27], homotopy analysis method [28], and adaptive observer [29].Moreover, very few algorithms for the analytical solution of a system of NFPDEs have been suggested, and some of these methods are essentially used for particular types of systems, often just linear ones or even smaller classes.Therefore, it should be noted that most of these methods cannot be generalized to nonlinear cases.
In the present work, we introduce a new analytical technique (NAT) to solve a full general system of NFPDEs of the following form: where L i and N i are linear and nonlinear operators, respectively, of ū = ū( x, t) and its partial derivatives, which might include other fractional partial derivatives of orders less than q i ; f i ( x, t) are known analytic functions; and D q i t are the Caputo partial derivatives of fractional orders q i , where we The goal of this paper is to demonstrate that a full general system of NFPDEs can be solved easily by using a NAT without any assumption and that it gives good results in analytical and numerical experiments.The rest of the paper is organized in as follows.In Section 2, we present basic definitions and preliminaries which are needed in the sequel.In Section 3, we introduce a NAT for solving a full general system of NFPDEs.Approximate analytical and numerical solutions for the systems of nonlinear wave equations are obtained in Section 4.

Basic Definitions and Preliminaries
There are various definitions and properties of fractional integrals and derivatives.In this section, we present modifications of some basic definitions and preliminaries of the fractional calculus theory, which are used in this paper and can be found in [10,[30][31][32][33][34][35].Definition 1.A real function u(x, t), x, t ∈ R, t > 0, is said to be in the space C µ , µ ∈ R if there exists a real number p(> µ), such that u(x, t) = t p u 1 (x, t), where u 1 (x, t) ∈ C(R × [0, ∞)), and it is said to be in the space C m µ if and only if Definition 2. Let q ∈ R \ N and q ≥ 0. The Riemann-Liouville fractional partial integral denoted by I q t of order q for a function u(x, t) ∈ C µ , µ > −1 is defined as: where Γ is the well-known Gamma function.
Theorem 1.Let q 1 , q 2 ∈ R \ N, q 1 , q 2 ≥ 0 and p > −1.For a function u(x, t) ∈ C µ , µ > −1, the operator I q t satisfies the following properties: Definition 3. Let q, t ∈ R, t > 0 and u(x, t) ∈ C m µ .Then is called the Caputo fractional partial derivative of order q for a function u(x, t).

NAT for Solving a System of NFPDEs
This section discusses a NAT to solve a system of NFPDEs.This NAT has much more computational power in obtaining piecewise analytical solutions.
To establish our technique, first we need to introduce the following results.
Then the system (1) admits at least a solution given by where L i(k−1)t denote the fractional partial integral of order q i for L i(k−1) and E i(k−1) respectively with respect to t.
Proof.Let the solution function u i ( x, t) of the system (6) to be as in the following analytical expansion: To solve system (1), we consider with initial conditions given by Next, we assume that, system (12) has a solution given by Performing Riemann-Liouville fractional partial integral of order q i with respect to t to both sides of system (12) and using Theorem 1, we obtain for i = 1, 2, . . ., n.By using the initial condition from the system (1), the system (15) can be rewritten as for i = 1, 2, . . ., n. Inserting ( 14) into ( 16), we obtain By using Lemma 1 and Theorem 3, the system (17) becomes Next, we use Definition 4 in the system (18), we obtain By equating the terms in system (17) with identical powers of λ, we obtain a series of the following systems Substituting the series (20) in the system (14) gives the solution of the system (12).Now, from the systems (11) and (14), we obtain By using the first equations of ( 21), we see that Inserting (20) into ( 21) completes the proof.

Convergence and Error Analysis
Theorem 5. Let B be a Banach space.Then the series solution of the system (20) Proof.Define the sequences S in , i = 1, 2, . . ., n of partial sums of the series given by the system (20) as . . .
and we need to show that {S in } are a Cauchy sequences in Banach space B. For this purpose, we consider For every n, m ∈ N, n ≥ m, by using the system (23) and triangle inequality successively, we have, Since 0 Therefore, the sequences {S in } are Cauchy sequences in the Banach space B, so the series solution defined in the system (21) converges.This completes the proof.Theorem 6.The maximum absolute truncation error of the series solution (11) of the nonlinear fractional partial differential system (1) is estimated to be sup where the region Ω ⊂ R n+1 .
Proof.From Theorem 5, we have But we assume that S in = ∑ n k=0 u ik ( x, t) for i = 1, 2, . . ., n, and since n → ∞, we obtain S in → u i ( x, t), so the system (28) can be rewritten as So, the maximum absolute truncation error in the region Ω is sup and this completes the proof.

Applications to the Systems of Nonlinear Wave Equations
In this section, we present examples of some systems of nonlinear wave equations.These examples are chosen because their closed form solutions are available, or they have been solved previously by some other well-known methods.
Example 1.Consider the nonlinear KdV system of time-fractional order of the form [24] For q = 1, the exact solitary wave solutions of the KdV system (31) is given by where the constant α is a wave velocity and β, γ are arbitrary constants.
To solve the system (31), we compare (31) with the system (1), we obtain where we assume Next, we assume the system (31) has a solution given by To obtain the approximate solution of the system (31), we consider the following system.
subject to the initial conditions given by and we assume that the system (36) has a solution of the form By operating Riemann-Liouville fractional partial integral of order q with respect to t for both sides of the system (36) and by using Theorem 2 and the system (37), we obtain By using Remark 1 and system (38), in the system (39), we obtain By equating the terms in the system (40) with identical powers of λ, we obtain a series of the following systems.
Hence the third-order term approximate solution for the system ( 31) is given by In Table 1, the numerical values of the approximate and exact solutions for Example 1 show the accuracy and efficiency of our technique at different values of x, t.The absolute error is listed for different values of x, t.In Figure 1a, we consider fixed values α = β = 0.5, γ = 1 and fixed order q = 1 for piecewise approximation values of x, t in the domain −20 ≤ x ≤ 20 and 0.20 ≤ t ≤ 1.In Figure 1b, we plot the exact solution with fixed values α = β = 0.5 and γ = 1 in the domain −20 ≤ x ≤ 20 and 0.20 ≤ t ≤ 1.
x t q = 0.5  Example 2. Consider the nonlinear dispersive long wave system of time fractional order [24-26] for 0 < q 1 , q 1 < 1, with initial condition given by u For q 1 = q 2 = 1, the system (43) has the following exact solitary wave solutions: where α, β are arbitrary constants.
By comparing the system (43) with the system (1), the system (43) can be rewritten as where N 1 (u, v) = −uu x and N 2 (u, v) = −(uv x + vu x ).To solve the system (46) by NAT discussed in Section 3, we assume that the system (46) has a solution given by Forgetting the approximate solution of the system (43), we consider the following system.
subject to the initial conditions given by Assume that the system (48) has a solution given by + − 4Γ (q 1 + q 2 + 1) e β+αx e β+αx + 1 e β+αx e β+αx − 3 + 1 Γ (q 1 + 1) Γ (q 2 + 1) Γ (q 1 + 2q Hence the third-order term approximate solution for the system (43) is given by In Table 2, we obtain similar approximation values as in the exact solution for Example 2 at different values of x, t.The absolute error is listed against different values of x, t.In Figure 2a, we plot the approximate solution at fixed values of the constants α = β = 0.5 and fixed order q = 1.In Figure 2b, we plot the exact solution with fixed values of the constants α = β = 0.5.The graphs of the approximate and exact solutions are drawing in the domain −20 ≤ x ≤ 20 and 0.20 ≤ t ≤ 1.

Discussion and Conclusions
In this paper, a nonlinear time-fractional partial differential system via NAT was studied.Moreover, the convergence and error analysis was also shown.From the computational point of view, the solutions obtained by our technique were in excellent agreement with those obtained via previous works and also conformed with the exact solution to confirm the effectiveness and accuracy of this technique.Moreover, approximate traveling wave solutions for some systems of nonlinear wave equations were successfully obtained.We used Mathematica software to obtain the approximate and numerical results as well as drawing the graphs.