Next Article in Journal
Design of the Algorithm for Packaging of Water Molecules in a Fixed Volume
Previous Article in Journal
Theoretical Study of Complexes of Tetravalent Actinides with DOTA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Reliable Way to Deal with the Coupled Fractional Korteweg-De Vries Equations within the Caputo Operator

by
Thongchai Botmart
1,
Badriah M. Alotaibi
2,*,
Rasool Shah
3,*,
Lamiaa S. El-Sherif
4,5 and
Samir A. El-Tantawy
6,7
1
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Physics, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
3
Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23200, Pakistan
4
Department of Physics, College of Arts and Science in Wadi Al-Dawaser, Prince Sattam bin Addulaziz University, Wadi-Dawaser 11991, Saudi Arabia
5
Department of Physics, Faculty of Science, Ain Shams University, Cairo 11566, Egypt
6
Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt
7
Research Center for Physics (RCP), Department of Physics, Faculty of Science and Arts, Al-Baha University, Al-Mikhwah 65431, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Symmetry 2022, 14(11), 2452; https://doi.org/10.3390/sym14112452
Submission received: 3 November 2022 / Revised: 8 November 2022 / Accepted: 15 November 2022 / Published: 18 November 2022
(This article belongs to the Section Mathematics)

Abstract

:
The development of numeric-analytic solutions and the construction of fractional order mathematical models for practical issues are of the highest concern in a variety of physics, applied mathematics, and engineering applications. The nonlinear Kersten–Krasil’shchik-coupled Korteweg–de Vries-modified Korteweg–de Vries (KdV-mKdV) system is treated analytically in this paper using a unique method, known as the Laplace residual power series (LRPS) approach to find some approximate solutions. The RPS methodology and the Laplace transform operator are combined in the LRPS method. We provide a detailed introduction to the proposed method for dealing with fractional Kersten–Krasil’shchik-linked KdV-mKdV models. When compared to exact solutions, the approach provides analytical solutions with good accuracy. We demonstrate the effectiveness of the current strategy compared to alternative methods for solving nonlinear equations using an illustrative example. The LRPS technique’s results show and highlight that the method may be used for a variety of time-fractional models of physical processes with simplicity and computing effectiveness.

1. Introduction

Fractional calculus (FC) has attracted increased interest from scientists and academics due to its applications in a variety of scientific disciplines, including engineering, chemistry, and social science [1,2,3]. Iterative methodology [4], Laplace transforms [5], and the operational tool [6] are some of the methods used to explore fractional equations. Different approaches were employed by various researchers to examine various fractional nonlinear equations found in the field of nanotechnology, including the shock wave equation, the KdV Burgers–Kuramoto equation, and the differential-difference equation [7,8,9]. Singh et al. studied the Tricomi equation in relation to the local fractional derivative of the fractal transonic flow, which led to the discovery of the system’s non-differentiable solution [10]. Using the integral transform method, Choudhary et al. [11] investigated a fractional system of heat flux in the semi-infinite solid. Liouville–Caputo operators were employed by Gomez Aguilar et al. in their study of a fractional derivative for an electrical RLC circuit [12]. In [13], a fractional-order model of the exponential kernel in the Euler–Lagrange and Hamilton equations was investigated. As a result, symmetry analysis is an excellent tool for understanding partial differential equations, particularly when looking at equations derived from accounting-related mathematical ideas [14,15]. The majority of natural observations lack symmetry, despite the fact that symmetry is the cornerstones of nature. Unexpected symmetry-breaking occurrences are a sophisticated method of hiding symmetry. Finite and infinitesimal symmetry are the two types. Finite symmetries can either be discrete or continuous. Parity and temporal inversion are discrete natural symmetries, whereas space is a continuous transformation. Patterns have always captivated mathematicians. In the seventeenth century, the classification of spatial and planar patterns really got going. Regrettably, precise solutions to fractional nonlinear differential equations have shown to be exceedingly challenging.
In recent years, people have been able to simulate numerous processes with more freedom by employing fractional differential equations (FDE). This concept has been realized in numerous disciplines, including engineering, economics, control theory, and finance, with impressive achievements. People are becoming more and more interested in using FDE to simulate challenging and real-world issues; however, integer characteristics do not make differential equations more realistic (IDE) [16,17,18]. FDE is non-local, has a storage effect, and performs better when compared to local IDE. For the investigation of various situations, the model’s future circumstances will vary depending on both current events and historical data. These features make it possible for FDE to successfully model non-Markov events in addition to non-Gaussian phenomena realistically. Additionally, standard IDEs are unable to explain FDE’s assistance or offer information between two different integer values. Numerous non-integer order derivative operators have been put out in various studies occurring to get over the restriction of only differentiating integer values. There are several disciplines where fractional-order differential operators are used [19,20,21,22,23,24].
The degree of freedom of its differential operator in modern calculus (including classical calculus) in a specific circumstance is more significant than that of the local ordinary differential operator. The main application of calculation can be found in [25,26,27]. As a result, the study of non-integer order differentiation and integrations is highly valued by scholars. From a geometric standpoint, the entire function’s accumulation or the full global integration range is explained by the arbitrary order derivatives, which are mostly definite integrals. The optimization of differential equations, and numerical and qualitative research has benefited tremendously from the work of researchers. It is important to remember that recent derivative operators were developed using definite integral techniques. It is a well-known reality that there is no fundamental formula for solving this. As a result, different definitions are involved with both types of kernels. The Atangana, Baleanu (ABC) fractional-order derivative [28] is the fundamental formulation that is the most attractive. Nonlinear equations are typically difficult to solve analytically or exactly. As a result, numerous numerical methods have been developed to evaluate the aforementioned equations. In order to examine FDE under ABC derivatives, many academics have recently looked into numerical approaches [29,30,31,32,33].
There are many published papers about the interpretation of nonlinear phenomena that can be generated and propagated in different plasma models. Most researchers have focused on reducing the basic equations of different plasma models to evolution equations in the form of partial differential equations such as the Korteweg–de Vries (KdV) equation, modified KdV (mKdV) equation, extended KdV equation, Kawahara-type equation, and so on [34,35,36]. A few studies reduced the fluid equations of different plasma models to fractional differential equations (FDEs) to gain a lot of information about how these waves propagate. Fractional coupled systems are frequently employed to study the complicated behavior of many nonlinear structures in different plasma models. Many experts have made an effort to assess this behavior. Recent research on the KdV equation and mKdV equation was conducted by Paul Kersten and Joseph Krasil’shchik. They proposed absolute complexity between coupled KdV-mKdV nonlinear systems for the analysis of nonlinear system behavior [37,38,39]. This Kersten–Krasil’shchik-coupled KdV-mKdV nonlinear system has been the subject of various variations, as noted in [40,41,42,43]. The nonlinear fractional Kersten–Krasil’shchik-coupled KdV-mKdV system is one of these variations and offers a mathematical explanation for the behavior of multi-component plasma for waves traveling down the positive κ axis:
D ρ u + u κ κ κ 6 u u κ + 3 v v κ κ κ + 3 v κ v κ κ κ 3 u κ v 2 + 6 u v v κ = 0 , > 0 , κ R , 0 < ρ 1 , D ρ v + v κ κ κ 3 v 2 v κ 3 u v κ + 3 u κ v = 0 , > 0 , κ R , 0 < ρ 1 ,
where κ denotes a spatial coordinate and denotes a time coordinate. Factor ρ denotes the order of the fractional operator. This operator is studied using the Caputo form. The fractional coupled system turns into a classical system when ρ = 1 , as shown below:
u + u κ κ κ 6 u u κ + 3 v v κ κ κ + 3 v κ v 2 κ 3 u κ v 2 + 6 u v v κ = 0 , > 0 , κ R , v + v κ κ κ 3 v 2 v κ 3 u v κ + 3 u κ v = 0 , > 0 , κ R .
The well-known KdV system is created from the Kersten–Krasil’shchik-linked KdV-mKdV system as follows if we put v = 0 .
u + u κ κ κ 6 u u κ = 0 , > 0 , κ R .
The well-known modified KdV system is created as follows when the Kersten–Krasil’shchik-connected KdV-mKdV system at u = 0 .
v + v κ κ κ 3 v 2 v κ = 0 , > 0 , κ R .
The Kersten–Krasil’shchik-linked KdV-mKdV system can, therefore, be thought of as a combination of the KdV system and the mKdV system, which are described by (2) to (4). In this work, we also look at the third-order KdV system with two components that are fractionally nonlinear and homogeneously coupled in time:
D ρ u u κ κ κ u u κ v v κ = 0 , > 0 , κ R , 0 < ρ 1 , D ρ v + 2 v κ κ κ u v κ = 0 , > 0 , κ R , 0 < ρ 1 ,
where ρ is the order factor of the fractional operator, is the temporal coordinate amd κ is the spatial coordinate. This operator is studied using Caputo form. The fractional coupled system turns into a classical system when ρ = 1 , as shown below:
u u κ κ κ u u κ v v κ = 0 , > 0 , κ R , v + 2 v κ κ κ u v κ = 0 , > 0 , κ R .
Finding the exact solution to nonlinear partial differential equations is still a major problem in physics and applied mathematics, necessitating the use of various techniques to obtain innovative approximate or exact solutions. Many approximation and numerical techniques have been used to solve fractional differential equations [44,45,46]. In this work, we propose a new easy and effective semi-analytical method to solve fractional PDE systems with variable coefficients. Our recommended method, LRPSM, which combines RPSM and the Laplace transform, was put into practice. In comparison to the RPS technique, which requires fractional differentiation in each phase, the main advantage of the present technique is that it can identify the unknown components of the suggested solutions by utilizing limits in the Laplace space, which, in turn, reduces the number of calculations needed and saves time [47,48,49]. For a variety of FDEs and time-fractional PDEs, the LRPS approach has been effectively used to develop approximative series solutions in closed forms [50,51,52].
The remaining parts of the ongoing work are structured as follows: The Laplace transform, Laplace fractional expansion, and certain fundamental definitions and theorems relating to fractional calculus are reviewed in Section 2. The concept of the suggested method for constructing the approximate solution of the fractional model under consideration (5) is provided in Section 3. To demonstrate the applicability and effectiveness of examining the solutions of time-PDEs of fractional order, the LRPS technique is used in Section 4 to solve fractional Kersten–Krasil’shchik-coupled KdV-mKdV systems. Finally, Section 5 presents possible interpretations of our findings.

2. Preliminaries

Here we presented some essential definitions related to our present work.
Definition 1. 
The Caputo fractional derivative is stated as [53,54,55]
C D ρ u ( κ , ) = J m ρ u m ( κ , ) , m 1 < ρ m , > 0 .
where m N and J ρ is the fractional Riemann-Liouville (RL) integral stated as
J ρ u ( κ , ) = 1 Γ ( ρ ) 0 ( t ) ρ 1 u ( κ , t ) d t
considering that the given integral exists.
Definition 2. 
For a function u ( κ , ) , the Laplace transform is stated as [53]
u ( κ , ϖ ) = L [ u ( κ , ) ] = 0 e ϖ u ( κ , ) d , ϖ > ρ ,
where the inverse LT is stated as
u ( κ , ) = L 1 [ u ( κ , ϖ ) ] = l i l + i e ϖ u ( κ , ϖ ) d ϖ , l = R e ( ϖ ) > l 0 ,
where l 0 is in the right half-plane of the Laplace integral’s absolute convergence.
Lemma 1. 
Let us consider that u ( κ , ) is a continuous piecewise function having exponential order ζ and U ( κ , ϖ ) = L [ u ( κ , ) ] , we get
1.
L [ J ρ u ( κ , ) ] = U ( κ , ϖ ) ϖ ρ , ρ > 0 .
2.
L [ D ρ u ( κ , ) ] = ϖ ρ U ( κ , ϖ ) k = 0 m 1 ϖ ρ k 1 u k ( κ , 0 ) , m 1 < ρ m .
3.
L [ D n ρ u ( κ , ) ] = ϖ n ρ U ( κ , ϖ ) k = 0 n 1 ϖ ( n k ) ρ 1 D k ρ u ( κ , 0 ) , 0 < ρ 1 .
Proof. 
Check Ref. [56]. □
Theorem 1. 
Let u ( κ , ) be a piecewise continuous function on I × [ 0 , ) with exponential order ζ. Assume the U ( κ , ϖ ) = L [ u ( κ , ) ] function has the below fractional expansion:
U ( κ , ϖ ) = n = 0 f n ( κ ) ϖ 1 + n ρ , 0 < ρ 1 , κ I , ϖ > ζ .
Thus, f n ( κ ) = D n ρ u ( κ , 0 ) .
Proof. 
Check Ref. [53]. □
Remark 1. 
The inverse LT of Equation (11) is stated as [53]:
u ( κ , ) = i = 0 D ρ u ( κ , 0 ) Γ ( 1 + i ρ ) i ( ζ ) , 0 < ζ 1 , 0 .
which is equal to the fractional Taylor’s formula as given in [57].
The FPS convergence of Theorem (1) is determined in the below theorem.
Theorem 2. 
Consider u ( κ , ) is piecewise continuous on I × [ 0 , ) having order ξ and is as given in Theorem (1), u ( κ , ϖ ) = L [ u ( κ , ) ] can be stated to be in the form of fractional Taylor’s formula. If ϖ L [ D i ρ + 1 u ( κ , ) ] M ( κ ) , on I × ( ξ , γ ] where 0 < ρ 1 , then R i ( κ , ϖ ) is the remainder of the form of fractional Taylor’s formula in Theorem (1), proving the below inequality.
| R i ( κ , ϖ ) | M ( κ ) S 1 + ( i + 1 ) ρ , κ I , ξ < ϖ γ .
Proof. 
Check Ref. [53]. □

3. Idea of LRPS

In this part, we will discuss the general methodology of LRPS to solve the system of fractional partial differential equations.
On taking LT of Equation (5), we have
U ( κ , ϖ ) f 0 ( κ , ξ ) ϖ + 1 ϖ ρ L L 1 [ U κ κ κ ] + L 1 [ U ] L 1 [ U κ ] + L 1 [ V ] L 1 [ V κ ] = 0 , V ( κ , ϖ ) g 0 ( κ , ξ ) ϖ + 1 ϖ ρ L 2 L 1 [ V κ κ κ ] L 1 [ U ] L 1 [ V κ ] = 0 .
Considering that Equation (14)’s solution has the below expansion
U ( κ , ϖ ) = n = 0 f n ( κ , ϖ ) ϖ n ρ + 1 , V ( κ , ϖ ) = n = 0 g n ( κ , ϖ ) ϖ n ρ + 1 .
The series for the k t h -truncated term is
U ( κ , ϖ ) = f 0 ( κ , ϖ ) ϖ + n = 1 k f n ( κ , ϖ ) ϖ n ρ + 1 , V ( κ , ϖ ) = g 0 ( κ , ϖ ) ϖ + n = 1 k g n ( κ , ϖ ) ϖ n ρ + 1 . k = 1 , 2 , 3 , 4
The Laplace residual functions (LRFs) [58] are
L R e s u ( κ , ϖ ) = U ( κ , ϖ ) f 0 ( κ , ϖ ) ϖ + 1 ϖ ρ L L 1 [ U κ κ κ ] + L 1 [ U ] L 1 [ U κ ] + L 1 [ V ] L 1 [ V κ ] , L R e s v ( κ , ϖ ) = V ( κ , ϖ ) g 0 ( κ , ϖ ) ϖ + 1 ϖ ρ L 2 L 1 [ V κ κ κ ] L 1 [ U ] L 1 [ V κ ] .
The k t h -LRFs as:
L R e s u , k ( κ , ϖ ) = U k ( κ , ϖ ) f 0 ( κ , ϖ ) ϖ + 1 ϖ ρ L L 1 [ U κ κ κ , k ] + L 1 [ U k ] L 1 [ U κ , k ] + L 1 [ V k ] L 1 [ V κ , k ] , L R e s v , k ( κ , ϖ ) = V k ( κ , ϖ ) g 0 ( κ , ϖ ) ϖ + 1 ϖ ρ L 2 L 1 [ V κ κ κ , k ] L 1 [ U k ] L 1 [ V κ , k ] .
There are some properties that arise in the LRPSM [58] to point out some facts:
  • L R e s ( κ , ϖ ) = 0 and lim j L R e s u , k ( κ , ϖ ) = L R e s u ( κ , ϖ ) for each ϖ > 0 .
  • lim ϖ ϖ L R e s u ( κ , ϖ ) = 0 lim ϖ ϖ L R e s u , k ( κ , ϖ ) = 0 .
  • lim ϖ ϖ k ρ + 1 L R e s u , k ( κ , ϖ ) = lim ϖ ϖ k ρ + 1 L R e s u , k ( κ , ϖ ) = 0 , 0 < ρ 1 , k = 1 , 2 , 3 , .
Now to get coefficients f n ( κ , ϖ ) and g n ( κ , ϖ ) , we recursively solve the below system
lim ϖ ϖ k ρ + 1 L R e s u , k ( κ , ϖ ) = 0 , k = 1 , 2 , , lim ϖ ϖ k ρ + 1 L R e s v , k ( κ , ϖ ) = 0 , k = 1 , 2 , ,
Lastly, we employ inverse LT to Equation (16) to obtain the k t h approximate solutions of u k ( κ , ) and v k ( κ , ) .

4. Numerical Problem

In this part, we implemented the suggested scheme for solving fractional Kersten–Krasil’shchik-linked KdV-mKdV systems.

4.1. Problem

Assume a time-fractional homogeneous two component-coupled KdV system of third-order as:
D ρ u u κ κ κ u u κ v v κ = 0 , D ρ v + 2 v κ κ κ u v κ = 0 .
subjected to initial conditions:
u ( κ , 0 ) = 3 6 tanh 2 κ 2 , v ( κ , 0 ) = 3 c 2 tanh κ 2 .
On taking LT from Equation (20) and by means of Equation (21), we obtain
U ( κ , ϖ ) 3 6 tanh 2 κ 2 ϖ 1 ϖ ρ L L 1 [ U κ κ κ ] + L 1 [ U ] L 1 [ U κ ] + L 1 [ V ] L 1 [ V κ ] = 0 , V ( κ , ϖ ) 3 c 2 tanh κ 2 ϖ 1 ϖ ρ L 2 L 1 [ V κ κ κ ] + L 1 [ U ] L 1 [ V κ ] = 0 ,
The series for k t h -truncated term is
U ( κ , ϖ ) = 3 6 tanh 2 κ 2 ϖ + n = 1 k f n ( κ , ϖ ) ϖ n ρ + 1 , V ( κ , ϖ ) = 3 c 2 tanh κ 2 ϖ + n = 1 k g n ( κ , ϖ ) ϖ n ρ + 1 , k = 1 , 2 , 3 , 4
and the k t h -LRFs as:
L R e s u , k ( κ , ϖ ) = U k ( κ , ϖ ) 3 6 tanh 2 κ 2 ϖ 1 ϖ ρ L L 1 [ U κ κ κ , k ] + L 1 [ U k ] L 1 [ U κ , k ] + L 1 [ V k ] L 1 [ V κ , k ] , L R e s v , k ( κ , ϖ ) = V k ( κ , ϖ ) 3 c 2 tanh κ 2 ϖ 1 ϖ ρ L 2 L 1 [ V κ κ κ , k ] + L 1 [ U k ] L 1 [ V κ , k ] .
Now, to find f k ( κ , ϖ ) and g k ( κ , ϖ ) , k = 1 , 2 , 3 , , we put the k t h -truncated series Equation (23) into the k t h -Laplace residual function Equation (24), multiply the resulting equation by ϖ k ρ + 1 , and now solve the relation recursively lim ϖ ( ϖ k ρ + 1 L R e s u , k ( κ , ϖ ) ) = 0 , and lim ϖ ( ϖ k ρ + 1 L R e s v , k ( κ , ϖ ) ) = 0 , k = 1 , 2 , 3 , .
We get some components as:
f 1 ( κ , ϖ ) = 3 6 tanh 2 κ 2 , g 1 ( κ , ϖ ) = 3 c 2 tanh κ 2 , f 2 ( κ , ϖ ) = 6 sech 2 κ 2 tanh κ 2 , g 2 ( κ , ϖ ) = 3 c 2 sech 2 κ 2 tanh κ 2 , f 3 ( κ , ϖ ) = 3 [ 2 + 7 sech 2 κ 2 15 sech 4 κ 2 ] sech 2 κ 2 g 3 ( κ , ϖ ) = 3 c 2 2 2 + 21 sech 2 κ 2 24 sech 4 κ 2 sech 2 κ 2 .
and so on.
Now, by substituting f k ( κ , ϖ ) and g k ( κ , ϖ ) , k = 1 , 2 , 3 , , in Equation (23), we obtain
U ( κ , ϖ ) = 3 6 tanh 2 κ 2 ϖ + 6 sech 2 κ 2 tanh κ 2 ϖ ρ + 1 + 3 [ 2 + 7 sech 2 κ 2 15 sech 4 κ 2 ] sech 2 κ 2 ϖ 2 ρ + 1 + , V ( κ , ϖ ) = 3 c 2 tanh κ 2 ϖ + 3 c 2 sech 2 κ 2 tanh κ 2 ϖ ρ + 1 + 3 c 2 2 2 + 21 sech 2 κ 2 24 sech 4 κ 2 sech 2 κ 2 ϖ 2 ρ + 1 + .
On taking the inverse LT, we have
u ( κ , ) = ( 3 6 tanh 2 κ 2 + 6 sech 2 κ 2 tanh κ 2 ρ Γ ( ρ + 1 ) + 3 [ 2 + 7 sech 2 κ 2 15 sech 4 κ 2 ] sech 2 κ 2 2 ρ Γ ( 2 ρ + 1 ) + ) , v ( κ , ) = ( 3 c 2 tanh κ 2 + 3 c 2 sech 2 κ 2 tanh κ 2 ρ Γ ( ρ + 1 ) + 3 c 2 2 2 + 21 sech 2 κ 2 24 sech 4 κ 2 sech 2 κ 2 2 ρ Γ ( 2 ρ + 1 ) + ) .
we obtain exact solutions by taking ρ = 1 as
u ( κ , ) = 3 6 tanh 2 κ + 2 , v ( κ , ) = 3 c 2 tanh κ + 2 .

4.2. Problem

Assume the time-fractional Kersten–Krasil’shchik-coupled KdV-mKdV system as:
D ρ u u κ κ κ 6 u u κ + 3 v v κ κ κ + 3 v κ v κ κ 3 u κ v 2 + 6 u v v κ = 0 , D ρ v + v κ κ κ 3 v 2 v κ 3 u v κ + 3 u κ v = 0 .
subjected to initial conditions:
u ( κ , 0 ) = c 2 c sech 2 ( c κ ) , v ( κ , 0 ) = 2 c sech ( c κ ) .
On taking the LT of Equation (28) and by means of Equation (29), we obtain
U ( κ , ϖ ) c 2 c sech 2 ( c κ ) ϖ 1 ϖ ρ L [ L 1 [ U κ κ κ ] + 6 L 1 [ U ] L 1 [ U κ ] 3 L 1 [ V ] L 1 [ V κ κ κ ] 3 L 1 [ V κ ] L 1 [ V κ κ ] + 3 L 1 [ U κ ] L 1 [ V 2 ] 6 L 1 [ U ] L 1 [ V ] L 1 [ V κ ] ] = 0 , V ( κ , ϖ ) 2 c sech ( c κ ) ϖ 1 ϖ ρ L [ L 1 [ V κ κ κ ] + 3 L 1 [ V 2 ] L 1 [ V κ ] + 3 L 1 [ U ] L 1 [ V κ ] 3 L 1 [ U κ ] L 1 [ V ] ] = 0 ,
The series for the k t h -truncated term is
U ( κ , ϖ ) = c 2 c sech 2 ( c κ ) ϖ + n = 1 k f n ( κ , ϖ ) ϖ n ρ + 1 , V ( κ , ϖ ) = 2 c sech ( c κ ) ϖ + n = 1 k g n ( κ , ϖ ) ϖ n ρ + 1 , k = 1 , 2 , 3 , 4
and the k t h -LRFs as:
L R e s u , k ( κ , ϖ ) = U k ( κ , ϖ ) c 2 c sech 2 ( c κ ) ϖ + 1 ϖ ρ L [ L 1 [ U κ κ κ , k ] + 6 L 1 [ U k ] L 1 [ U κ , k ] 3 L 1 [ V k ] L 1 [ V κ κ κ , k ] 3 L 1 [ V κ , k ] L 1 [ V κ κ , k ] + 3 L 1 [ U κ , k ] L 1 [ V k 2 ] 6 L 1 [ U k ] L 1 [ V k ] L 1 [ V κ , k ] ] , L R e s v , k ( κ , ϖ ) = V k ( κ , ϖ ) 2 c sech ( c κ ) ϖ + 1 ϖ ρ L [ L 1 [ V κ κ κ , k ] + 3 L 1 [ V k 2 ] L 1 [ V κ , k ] + 3 L 1 [ U k ] L 1 [ V κ , k ] 3 L 1 [ U κ , k ] L 1 [ V k ] ] .
Now, to find f k ( κ , ϖ ) and g k ( κ , ϖ ) , k = 1 , 2 , 3 , , we put the k t h -truncated series Equation (31) into the k t h -Laplace residual function Equation (32), multiply the resulting equation by s k ρ + 1 , and now solve the relation recursively lim ϖ ( ϖ k ρ + 1 L R e s u , k ( κ , ϖ ) ) = 0 , and lim ϖ ( ϖ k ρ + 1 L R e s v , k ( κ , ϖ ) ) = 0 , k = 1 , 2 , 3 , .
We get some components as:
f 1 ( κ , ϖ ) = c 2 c sech 2 ( c κ ) , g 1 ( κ , ϖ ) = 2 c sech ( c κ ) , f 2 ( κ , ϖ ) = 8 c 5 2 sinh ( c κ ) sech 3 ( c κ ) , g 2 ( κ , ϖ ) = 4 c 2 sinh ( c κ ) sech 2 ( c κ ) , f 3 ( κ , ϖ ) = 16 c 4 [ 2 cosh 2 ( c κ ) 3 ] sech 4 ( c κ ) g 3 ( κ , ϖ ) = 8 c 7 2 [ cosh 2 ( c κ ) 2 ] sech 3 ( c κ ) .
and so on.
Now, by substituting f k ( κ , ϖ ) and g k ( κ , ϖ ) , k = 1 , 2 , 3 , , into Equation (31), we obtain
U ( κ , ϖ ) = c 2 c sech 2 ( c κ ) ϖ + 8 c 5 2 sinh ( c κ ) sech 3 ( c κ ) ϖ ρ + 1 + 16 c 4 [ 2 cosh 2 ( c κ ) 3 ] sech 4 ( c κ ) ϖ 2 ρ + 1 + , V ( κ , ϖ ) = 2 c sech ( c κ ) ϖ + 4 c 2 sinh ( c κ ) sech 2 ( c κ ) ϖ ρ + 1 + 8 c 7 2 [ cosh 2 ( c κ ) 2 ] sech 3 ( c κ ) ϖ 2 ρ + 1 + .
On taking inverse LT, we have
u ( κ , ) = ( c 2 c sech 2 ( c κ ) + 8 c 5 2 sinh ( c κ ) sech 3 ( c κ ) ρ Γ ( ρ + 1 ) 16 c 4 [ 2 cosh 2 ( c κ ) 3 ] sech 4 ( c κ ) 2 ρ Γ ( 2 ρ + 1 ) + ) , v ( κ , ) = ( 2 c sech ( c κ ) 4 c 2 sinh ( c κ ) sech 2 ( c κ ) ρ Γ ( ρ + 1 ) + 8 c 7 2 [ cosh 2 ( c κ ) 2 ] sech 3 ( c κ ) 2 ρ Γ ( 2 ρ + 1 ) + ) .
we obtain exact solutions by taking ρ = 1 as
u ( κ , ) = c 2 c sech 2 ( c ( κ + 2 c ) ) , v ( κ , ) = 2 c sech ( c ( κ + 2 c ) ) .

5. Results and Discussion

The solutions to u ( κ , ) using the exact and suggested approach are calculated in Figure 1 with ρ = 1 . The graphical representations of u ( κ , ) for ρ = 0.8 and 0.6 are shown in Figure 2. Figure 3 illustrates the 2D and 3D behavior of u ( κ , ) for various fractional orders. The solutions to v ( κ , ) using the actual and suggested approach are calculated in Figure 4 with ρ = 1 . The graphical representations of v ( κ , ) for ρ = 0.8 and 0.6 are shown in Figure 5. Figure 6 illustrates the 2D and 3D behavior of v ( κ , ) for various fractional orders. The graphical representation of System 1 demonstrates great agreement between our solutions and the correct answer. Similar to Figure 6, Figure 7 presents the actual and suggested methods solutions for u ( κ , ) for ρ = 1 .  Figure 8 exhibits the graphical representations of u ( κ , ) for ρ = 0.8 and 0.6 , respectively, while Figure 9 depicts the 2D and 3D behaviors of u ( κ , ) for various fractional orders. Additionally, Figure 10 shows both the actual and proposed technique’s solutions for v ( κ , ) for ρ = 1 . Figure 11 presents the graphical representations of v ( κ , ) for ρ = 0.8 and 0.6 , respectively, while Figure 12 depicts the 2D and 3D behaviors of v ( κ , ) for various fractional orders. In a similar manner, System 2’s graphical appearance demonstrates the fact that our solutions are in close agreement with the accurate solutions. Each and every figure was generated at c = 0.1 and [ 0 , 0.1 ] within the range 3 κ 3 . Additionally, Table 1 and Table 2 display the approximate solution to System 1 for various values of κ and , while the approximate solution of System 2 for various values of κ and is shown in Table 3 and Table 4, respectively. The figures and tables show that the exact solution and our approach solution are quite similar to one another and have a greater level of precision. Further, from the figures and tables, it is observed that the proposed method solution gets closer to the exact solution as the value of ρ tends from the fractional order toward the integer order. Moreover, from the comparison, it is clear that the LRPS technique solution and the exact solutions are very close. Thus, the LRPS technique is a dependable new study that requires less computation, is adaptable, and is simple to use.

6. Conclusions

The fractional Kersten–Krasil’shchik-linked KdV-mKdV analytical approximate solution, together with the necessary initial data, was successfully developed in this work using the LRPS technique. The major goal of the suggested method is to use the limit concept to find the unknown LFSE coefficients for the new equation in Laplace space. Without perturbation, discretization, or physical hypotheses, the analytical approximations for the solved fractional Kersten–Krasil’shchik-coupled KdV-mKdV systems’ starting value equations are obtained in rapidly convergent MFPS formulas. Two illustrated examples were used to study the LRPS technique’s performance and reliability. As a result, the LRPS technique is a straightforward, simple, and useful tool for treating a variety of nonlinear time-fractional PDEs that occur in engineering and science problems.

Author Contributions

Conceptualization, T.B. and B.M.A.; methodology, R.S.; software, L.S.E.-S.; validation, S.A.E.-T., T.B. and R.S.; formal analysis, B.M.A.; investigation, L.S.E.-S.; resources, S.A.E.-T.; data curation, T.B.; writing—original draft preparation, R.S.; writing—review and editing, S.A.E.-T.; visualization, B.M.A.; supervision, S.A.E.-T.; project administration, B.M.A.; funding acquisition, T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project (Grant No. PNURSP2022R32), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This research received funding support from the NSRF via the Program Management Unit for Human Resources & Institutional Development, Research and Innovation, (grant number B05F650018).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  2. Odibat, Z.M.; Momani, S. Application of Variational Iteration Method to Nonlinear Differential Equations of Fractional Order. Int. J. Nonlinear Sci. Numer. Simul. 2006, 7, 27–34. [Google Scholar] [CrossRef]
  3. Razminia, K.; Razminia, A.; Machado, J.A.T. Analytical Solution of Fractional Order Diffusivity Equation With Wellbore Storage and Skin Effects. J. Comput. Nonlinear Dyn. 2016, 11, 011006. [Google Scholar] [CrossRef]
  4. Diethelm, K.; Ford, N. Analysis of Fractional Differential Equations. J. Math. Anal. Appl. 2002, 265, 229–248. [Google Scholar] [CrossRef] [Green Version]
  5. Diethelm, K. An algorithm for the numerical solution of differential equations of fractional order. Electron. Trans. Numer. Anal. 1997, 5, 1–6. [Google Scholar]
  6. Luchko, Y.; Srivastava, H. The exact solution of certain differential equations of fractional order by using operational calculus. Comput. Math. Appl. 1995, 29, 73–85. [Google Scholar] [CrossRef] [Green Version]
  7. Kumar, D.; Singh, J.; Kumar, S.; Singh, B.P. Numerical computation of nonlinear shock wave equation of fractional order. Ain Shams Eng. J. 2015, 6, 605–611. [Google Scholar] [CrossRef] [Green Version]
  8. Bulut, H.; Belgacem, F.B.M.; Baskonus, H.M. Some new analytical solutions for the nonlinear time-fractional KdV-Burgers-Kuramoto equation. Adv. Math. Stat. Sci. 2015, 2015, 118–129. [Google Scholar]
  9. Kumar, D.; Singh, J.; Baleanu, D. Numerical Computation of a Fractional Model of Differential-Difference Equation. J. Comput. Nonlinear Dyn. 2016, 11, 061004. [Google Scholar] [CrossRef]
  10. Singh, J.; Kumar, D.; Nieto, J.J. A Reliable Algorithm for a Local Fractional Tricomi Equation Arising in Fractal Transonic Flow. Entropy 2016, 18, 206. [Google Scholar] [CrossRef] [Green Version]
  11. Choudhary, A.; Kumar, D.; Singh, J. Numerical simulation of a fractional model of temperature distribution and heat flux in the semi infinite solid. Alex. Eng. J. 2016, 55, 87–91. [Google Scholar] [CrossRef] [Green Version]
  12. Gómez-Aguilar, J.F.; Morales-Delgado, V.F.; Taneco-Hernández, M.A.; Baleanu, D.; Escobar-Jiménez, R.F.; Al Qurashi, M.M. Analytical solutions of the electrical RLC circuit via Liouville-Caputo operators with local and non-local kernels. Entrop 2016, 18, 402. [Google Scholar] [CrossRef]
  13. Coronel-Escamilla, A.; Gómez-Aguilar, J.F.; Baleanu, D.; Escobar-Jiménez, R.F.; Olivares-Peregrino, V.H.; Abúndez-Pliego, A. Formulation of Euler-Lagrange and Hamilton equations involving fractional operators with regular kernel. Adv. Differ. Equat. 2016, 2016, 283. [Google Scholar] [CrossRef]
  14. Dai, B.; Zhang, B.; Niu, Z.; Feng, Y.; Liu, Y.; Fan, Y. A Novel Ultrawideband Branch Waveguide Coupler With Low Amplitude Imbalance. IEEE Trans. Microw. Theory Tech. 2022, 70, 3838–3846. [Google Scholar] [CrossRef]
  15. Feng, Y.; Zhang, B.; Liu, Y.; Niu, Z.; Fan, Y.; Chen, X. A D-Band Manifold Triplexer With High Isolation Utilizing Novel Waveguide Dual-Mode Filters. IEEE Trans. Terahertz Sci. Technol. 2022, 12, 678–681. [Google Scholar] [CrossRef]
  16. Kovalnogov, V.N.; Kornilova, M.I.; Khakhalev, Y.A.; Generalov, D.A.; Simos, T.; Tsitouras, C. New family for Runge-Kutta-Nyström pairs of orders 6(4) with coefficients trained to address oscillatory problems. Math. Methods Appl. Sci. 2022, 45, 7715–7727. [Google Scholar] [CrossRef]
  17. Hao, R.-B.; Lu, Z.-Q.; Ding, H.; Chen, L.-Q. A nonlinear vibration isolator supported on a flexible plate: Analysis and experiment. Nonlinear Dyn. 2022, 108, 941–958. [Google Scholar] [CrossRef]
  18. Luo, G.; Yuan, Q.; Li, J.; Wang, S.; Yang, F. Artificial Intelligence Powered Mobile Networks: From Cognition to Decision. IEEE Netw. 2022, 36, 136–144. [Google Scholar] [CrossRef]
  19. Wang, J.; Tian, J.; Zhang, X.; Yang, B.; Liu, S.; Yin, L.; Zheng, W. Control of Time Delay Force Feedback Teleoperation System with Finite Time Convergence. Front. Neurorobot. 2022, 16, 877069. [Google Scholar] [CrossRef]
  20. Lu, S.; Ban, Y.; Zhang, X.; Yang, B.; Liu, S.; Yin, L.; Zheng, W. Adaptive control of time delay teleoperation system with uncertain dynamics. Front. Neurorobot. 2022, 16, 928863. [Google Scholar] [CrossRef]
  21. Abdilraze, A. Adomian Decomposition Method: Convergence Analysis and Numerical Approximations. Master’s Dissertation, McMaster University Hamilton, Hamilton, ON, Canada, 2008. [Google Scholar]
  22. Khan, A.; Khan, T.S.; Syam, M.I.; Khan, H. Analytical solutions of time-fractional wave equation by double Laplace transform method. Eur. Phys. J. Plus 2019, 134, 163. [Google Scholar] [CrossRef]
  23. Kumar, A.; Kumar, S.; Yan, S.-P. Residual Power Series Method for Fractional Diffusion Equations. Fundam. Informaticae 2017, 151, 213–230. [Google Scholar]
  24. Baleanu, D.; Ghanbari, B.; Asad, J.H.; Jajarmi, A.; Pirouz, H.M. Planar System-Masses in an Equilateral Triangle: Numerical Study within Fractional Calculus. Comput. Model. Eng. Sci. 2020, 124, 953–968. [Google Scholar] [CrossRef]
  25. Zhang, Y.-Z.; Yang, A.-M.; Long, Y. Initial boundary value problem for fractal heat equation in the semi-infinite region by Yang-Laplace transform. Therm. Sci. 2014, 18, 677–681. [Google Scholar] [CrossRef]
  26. Spiga, G.; Spiga, M. Two-Dimensional Transient Solutions for Crossflow Heat Exchangers With Neither Gas Mixed. J. Heat Transf. 1987, 109, 281–286. [Google Scholar] [CrossRef]
  27. Hassan, H.E.; Kilicman, A. A note on solutions of wave, Laplace’s and heat equations with convolution terms by using a double Laplace transform. Appl. Math. Lett. 2008, 21, 1324–1329. [Google Scholar]
  28. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm Sci. 2016, 20, 763–785. [Google Scholar] [CrossRef] [Green Version]
  29. Alshehry, A.S.; Imran, M.; Khan, A.; Shah, R.; Weera, W. Fractional View Analysis of Kuramoto–Sivashinsky Equations with Non-Singular Kernel Operators. Symmetry 2022, 14, 1463. [Google Scholar] [CrossRef]
  30. Botmart, T.; Agarwal, R.P.; Naeem, M.; Khan, A.; Shah, R. On the solution of fractional modified Boussinesq and approximate long wave equations with non-singular kernel operators. AIMS Math. 2022, 7, 12483–12513. [Google Scholar] [CrossRef]
  31. Alaoui, M.K.; Nonlaopon, K.; Zidan, A.M.; Khan, A.; Shah, R. Analytical Investigation of Fractional-Order Cahn–Hilliard and Gardner Equations Using Two Novel Techniques. Mathematics 2022, 10, 1643. [Google Scholar] [CrossRef]
  32. Shah, N.A.; Hamed, Y.S.; Abualnaja, K.M.; Chung, J.-D.; Shah, R.; Khan, A. A Comparative Analysis of Fractional-Order Kaup–Kupershmidt Equation within Different Operators. Symmetry 2022, 14, 986. [Google Scholar] [CrossRef]
  33. Shah, N.A.; El-Zahar, E.R.; Akgül, A.; Khan, A.; Kafle, J. Analysis of Fractional-Order Regularized Long-Wave Models via a Novel Transform. J. Funct. Spaces 2022, 2022, 2754507. [Google Scholar] [CrossRef]
  34. Alharthi, M.R.; Alharbey, R.A.; El-Tantawy, S.A. Novel analytical approximations to the nonplanar Kawahara equation and its plasma applications. Eur. Phys. J. Plus 2022, 137, 1172. [Google Scholar] [CrossRef]
  35. Alyousef, H.A.; Salas, A.H.; Matoog, R.T.; El-Tantawy, S.A. On the analytical and numerical approximations to the forced damped Gardner Kawahara equation and modeling the nonlinear structures in a collisional plasma. Phys. Fluids 2022, 34, 103105. [Google Scholar] [CrossRef]
  36. Alyousef, H.A.; Salas, A.H.; Alharthi, M.R.; El-Tantawy, S.A. New Periodic and Localized Traveling Wave Solutions to a Kawahara-Type Equation: Applications to Plasma Physics. Complexity 2022, 2022, 9942267. [Google Scholar] [CrossRef]
  37. Naeem, M.; Zidan, A.; Nonlaopon, K.; Syam, M.; Al-Zhour, Z.; Shah, R. A New Analysis of Fractional-Order Equal-Width Equations via Novel Techniques. Symmetry 2021, 13, 886. [Google Scholar] [CrossRef]
  38. Iqbal, N.; Akgül, A.; Shah, R.; Bariq, A.; Al-Sawalha, M.M.; Ali, A. On Solutions of Fractional-Order Gas Dynamics Equation by Effective Techniques. J. Funct. Spaces 2022, 2022, 3341754. [Google Scholar] [CrossRef]
  39. Sunthrayuth, P.; Zidan, A.M.; Yao, S.-W.; Shah, R.; Inc, M. The comparative study for solving fractional-order Fornberg-Whitham equation via ρ-Laplace transform. Symmetry 2021, 13, 784. [Google Scholar] [CrossRef]
  40. Qasim, A.F.; Al-Amr, M.O. Approximate solution of the Kersten-Krasil’shchik coupled Kdv-mKdV system via reduced differential transform method. Eurasian J. Sci. Eng. 2018, 4, 1–9. [Google Scholar]
  41. Kalkanli, A.K.; Sakovich, S.Y.; Yurdusen, I. Integrability of Kersten-Krasil’shchik coupled KdV-mKdV equations: Singularity analysis and Lax pair. J. Math. Phys. 2003, 44, 1703–1708. [Google Scholar] [CrossRef] [Green Version]
  42. Hon, Y.C.; Fan, E.G. Solitary wave and doubly periodic wave solutions for the Kersten-Krasil’shchik coupled KdV-mKdV system. Chaos Solitons Fract. 2004, 19, 1141–1146. [Google Scholar] [CrossRef]
  43. Al-Sawalha, M.M.; Khan, A.; Ababneh, O.Y.; Botmart, T. Fractional view analysis of Kersten-Krasil’shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Math. 2022, 7, 18334–18359. [Google Scholar] [CrossRef]
  44. Hietarinta, J. Introduction to the Hirota bilinear method. In Integrability of Nonlinear Systems; Springer: Berlin/Heidelberg, Germany, 1997; pp. 95–103. [Google Scholar]
  45. Shah, N.A.; Agarwal, P.; Chung, J.D.; El-Zahar, E.R.; Hamed, Y.S. Analysis of optical solitons for nonlinear Schrödinger equation with detuning term by iterative transform method. Symmetry 2020, 12, 1850. [Google Scholar] [CrossRef]
  46. Sunthrayuth, P.; Alyousef, H.A.; El-Tantawy, S.A.; Khan, A.; Wyal, N. Solving Fractional-Order Diffusion Equations in a Plasma and Fluids via a Novel Transform. J. Funct. Spaces 2022, 2022, 1899130. [Google Scholar] [CrossRef]
  47. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, U.K.S. Computational Optimization of Residual Power Series Algorithm for Certain Classes of Fuzzy Fractional Differential Equations. Int. J. Differ. Equations 2018, 2018, 8686502. [Google Scholar] [CrossRef]
  48. Prakasha, D.G.; Veeresha, P.; Baskonus, H.M. Residual Power Series Method for Fractional Swift–Hohenberg Equation. Fractal Fract. 2019, 3, 9. [Google Scholar] [CrossRef]
  49. Bataineh, M.; Alaroud, M.; Al-Omari, S.; Agarwal, P. Series Representations for Uncertain Fractional IVPs in the Fuzzy Conformable Fractional Sense. Entropy 2021, 23, 1646. [Google Scholar] [CrossRef]
  50. Alaroud, M. Application of Laplace residual power series method for approximate solutions of fractional IVP’s. Alex. Eng. J. 2022, 61, 1585–1595. [Google Scholar] [CrossRef]
  51. Burqan, A.; El-Ajou, A.; Saadeh, R.; Al-Smadi, M. A new efficient technique using Laplace transforms and smooth expansions to construct a series solution to the time-fractional Navier-Stokes equations. Alex. Eng. J. 2022, 61, 1069–1077. [Google Scholar] [CrossRef]
  52. Alquran, M.; Ali, M.; Alsukhour, M.; Jaradat, I. Promoted residual power series technique with Laplace transform to solve some time-fractional problems arising in physics. Results Phys. 2020, 19, 103667. [Google Scholar] [CrossRef]
  53. Alderremy, A.A.; Aly, S.; Fayyaz, R.; Khan, A.; Shah, R.; Wyal, N. The Analysis of Fractional-Order Nonlinear Systems of Third Order KdV and Burgers Equations via a Novel Transform. Complexity 2022, 2022, 4935809. [Google Scholar] [CrossRef]
  54. Zidan, A.M.; Khan, A.; Shah, R.; Alaoui, M.K.; Weera, W. Evaluation of time-fractional Fishe’s equations with the help of analytical methods. AIMS Math. 2022, 7, 18746–18766. [Google Scholar] [CrossRef]
  55. Al-Sawalha, M.M.; Ababneh, O.Y.; Shah, R.; Khan, A.; Nonlaopon, K. Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators. AIMS Math. 2023, 8, 2308–2336. [Google Scholar] [CrossRef]
  56. Areshi, M.; Khan, A.; Shah, R.; Nonlaopon, K. Analytical investigation of fractional-order Newell-Whitehead-Segel equations via a novel transform. AIMS Math. 2022, 7, 6936–6958. [Google Scholar] [CrossRef]
  57. Abu Arqub, O.; El-Ajou, A.; Momani, S. Constructing and predicting solitary pattern solutions for nonlinear time-fractional dispersive partial differential equations. J. Comput. Phys. 2015, 293, 385–399. [Google Scholar] [CrossRef]
  58. Alshehry, A.S.; Shah, R.; Shah, N.A.; Dassios, I. A Reliable Technique for Solving Fractional Partial Differential Equation. Axioms 2022, 11, 574. [Google Scholar] [CrossRef]
Figure 1. The accurate and suggested approach solution for u ( κ , ) for System 1 at ρ = 1 .
Figure 1. The accurate and suggested approach solution for u ( κ , ) for System 1 at ρ = 1 .
Symmetry 14 02452 g001
Figure 2. The suggested approach solution for System 1 at ρ = 0.8 , 0.6 of u ( κ , ) .
Figure 2. The suggested approach solution for System 1 at ρ = 0.8 , 0.6 of u ( κ , ) .
Symmetry 14 02452 g002
Figure 3. The solution of the suggested approach for u ( κ , ) of System 1 at various ρ values.
Figure 3. The solution of the suggested approach for u ( κ , ) of System 1 at various ρ values.
Symmetry 14 02452 g003
Figure 4. The accurate and suggested approach solution for v ( κ , ) for System 1 at ρ = 1 .
Figure 4. The accurate and suggested approach solution for v ( κ , ) for System 1 at ρ = 1 .
Symmetry 14 02452 g004
Figure 5. The suggested approach solution for System 1 at ρ = 0.8 , 0.6 of v ( κ , ) .
Figure 5. The suggested approach solution for System 1 at ρ = 0.8 , 0.6 of v ( κ , ) .
Symmetry 14 02452 g005
Figure 6. The solution of the suggested approach for v ( κ , ) of System 1 at various ρ values.
Figure 6. The solution of the suggested approach for v ( κ , ) of System 1 at various ρ values.
Symmetry 14 02452 g006
Figure 7. The accurate and suggested approach solution for u ( κ , ) for System 2 at ρ = 1 .
Figure 7. The accurate and suggested approach solution for u ( κ , ) for System 2 at ρ = 1 .
Symmetry 14 02452 g007
Figure 8. The suggested approach solution for System 2 at ρ = 0.8 , 0.6 of u ( κ , ) .
Figure 8. The suggested approach solution for System 2 at ρ = 0.8 , 0.6 of u ( κ , ) .
Symmetry 14 02452 g008
Figure 9. The solution of the suggested approach for u ( κ , ) of System 2 at various ρ values.
Figure 9. The solution of the suggested approach for u ( κ , ) of System 2 at various ρ values.
Symmetry 14 02452 g009
Figure 10. The accurate and suggested approach solution for v ( κ , ) for System 2 at ρ = 1 .
Figure 10. The accurate and suggested approach solution for v ( κ , ) for System 2 at ρ = 1 .
Symmetry 14 02452 g010
Figure 11. The suggested approach solution for System 2 at ρ = 0.8 , 0.6 of v ( κ , ) .
Figure 11. The suggested approach solution for System 2 at ρ = 0.8 , 0.6 of v ( κ , ) .
Symmetry 14 02452 g011
Figure 12. The solution of the suggested approach for v ( κ , ) of System 2 at various ρ values.
Figure 12. The solution of the suggested approach for v ( κ , ) of System 2 at various ρ values.
Symmetry 14 02452 g012
Table 1. Solution to u ( κ , ) at various fractional orders for System 1.
Table 1. Solution to u ( κ , ) at various fractional orders for System 1.
( κ , ) u ( κ , ) at ρ = 0.5 u ( κ , ) at ρ = 0.75 ( LRPSM ) at ρ = 1Exact Result
(0.2, 0.01)2.95229512.94630312.94039772.9403977
(0.4, 0.01)2.78920152.77722632.76625782.7662578
(0.6, 0.01)2.52279232.50680212.49082172.4908217
(0.2, 0.02)2.95232422.94635222.94039772.9403977
(0.4, 0.02)2.78920412.77621902.76625782.7662578
(0.6, 0.02)2.52278632.50680152.49082172.4908217
(0.2, 0.03)2.95232452.94636212.94039772.9403977
(0.4, 0.03)2.78919872.77621242.76625782.7662578
(0.6, 0.03)2.52279212.50681052.49082172.4908217
(0.2, 0.04)2.95231832.94635612.94039772.9403977
(0.4, 0.04)2.78921242.77623212.76625782.7662578
(0.6, 0.04)2.52278332.50681032.49082172.4908217
(0.2, 0.05)2.95231332.94635622.94039772.9403977
(0.4, 0.05)2.78921892.77624882.76625782.7662578
(0.6, 0.05)2.52280212.50681112.49082162.4908216
Table 2. Solution to v ( κ , ) at various fractional orders for System 1.
Table 2. Solution to v ( κ , ) at various fractional orders for System 1.
( κ , ) v ( κ , ) at ρ = 0.5 v ( κ , ) at ρ = 0.75 ( LRPSM ) at ρ = 1Exact Result
(0.2, 0.01)−0.0004074−0.0004127−0.0004228−0.0004228
(0.4, 0.01)−0.0008051−0.0008202−0.0008373−0.0008373
(0.6, 0.01)−0.0011923−0.0012126−0.0012359−0.0012359
(0.2, 0.02)−0.0004024−0.0004087−0.0004228−0.0004228
(0.4, 0.02)−0.0008051−0.0008212−0.0008373−0.0008373
(0.6, 0.02)−0.0012003−0.0012217−0.0012359−0.0012359
(0.2, 0.03)−0.0004034−0.0004126−0.0004228−0.0004228
(0.4, 0.03)−0.0008111−0.0008202−0.0008373−0.0008373
(0.6, 0.03)−0.0012004−0.0012127−0.0012359−0.0012359
(0.2, 0.04)−0.0004024−0.0004106−0.0004228−0.0004228
(0.4, 0.04)−0.0008061−0.0008203−0.0008374−0.0008374
(0.6, 0.04)−0.0012005−0.0012127−0.0012359−0.0012359
(0.2, 0.05)−0.0004033−0.0004126−0.0004228−0.0004228
(0.4, 0.05)−0.0008072−0.0008213−0.0008374−0.0008374
(0.6, 0.05)−0.0012016−0.0012137−0.0012359−0.0012359
Table 3. Solution to u ( κ , ) at various fractional orders for System 2.
Table 3. Solution to u ( κ , ) at various fractional orders for System 2.
( κ , ) u ( κ , ) at ρ = 0.5 u ( κ , ) at ρ = 0.75 ( LRPSM ) at ρ = 1Exact Result
(0.2, 0.01)−0.0991522−0.0991843−0.0992005−0.0992005
(0.4, 0.02)−0.0967362−0.0967994−0.0968306−0.0968306
(0.6, 0.03)−0.0928263−0.0929185−0.0929647−0.0929647
(0.2, 0.01)−0.0991503−0.0991826−0.0991989−0.0991989
(0.4, 0.02)−0.0967322−0.0967953−0.0968275−0.0968275
(0.6, 0.03)−0.0928212−0.0929134−0.0929601−0.0929601
(0.2, 0.01)−0.0991480−0.0991812−0.0991973−0.0991973
(0.4, 0.02)−0.0967281−0.0967922−0.0968244−0.0968244
(0.6, 0.03)−0.0928152−0.0929094−0.0929556−0.0929556
(0.2, 0.01)−0.0991473−0.0991795−0.0991957−0.0991957
(0.4, 0.02)−0.0967250−0.0967891−0.0968212−0.0968212
(0.6, 0.03)−0.0928102−0.0929044−0.0929510−0.0929510
(0.2, 0.01)−0.0991454−0.0991776−0.0991941−0.0991941
(0.4, 0.02)−0.0967214−0.0967866−0.0968181−0.0968181
(0.6, 0.03)−0.0928051−0.0928992−0.0929464−0.0929464
Table 4. Solution to v ( κ , ) at various fractional orders for System 2.
Table 4. Solution to v ( κ , ) at various fractional orders for System 2.
( κ , ) v ( κ , ) at ρ = 0.5 v ( κ , ) at ρ = 0.75 ( LRPSM ) at ρ = 1Exact Result
(0.2, 0.01)0.63111400.63116410.63119020.6311902
(0.4, 0.02)0.62727310.62737430.62742440.6274244
(0.6, 0.03)0.62100900.62115810.62123220.6212322
(0.2, 0.01)0.63111120.63116240.63118760.6311876
(0.4, 0.02)0.62726720.62736830.62741940.6274194
(0.6, 0.03)0.62100130.62115060.62122480.6212248
(0.2, 0.01)0.63110830.63115960.63118510.6311851
(0.4, 0.02)0.62726200.62736320.62741440.6274144
(0.6, 0.03)0.62099210.62114230.62121750.6212175
(0.2, 0.01)0.63110520.63115640.63118260.6311826
(0.4, 0.02)0.62725610.62735830.62740940.6274094
(0.6, 0.03)0.62098440.62113580.62121010.6212101
(0.2, 0.01)0.63110230.63114570.63118000.6311800
(0.4, 0.02)0.62725000.62735320.62740440.6274044
(0.6, 0.03)0.62097640.62112750.62120270.6212027
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Botmart, T.; Alotaibi, B.M.; Shah, R.; El-Sherif, L.S.; El-Tantawy, S.A. A Reliable Way to Deal with the Coupled Fractional Korteweg-De Vries Equations within the Caputo Operator. Symmetry 2022, 14, 2452. https://doi.org/10.3390/sym14112452

AMA Style

Botmart T, Alotaibi BM, Shah R, El-Sherif LS, El-Tantawy SA. A Reliable Way to Deal with the Coupled Fractional Korteweg-De Vries Equations within the Caputo Operator. Symmetry. 2022; 14(11):2452. https://doi.org/10.3390/sym14112452

Chicago/Turabian Style

Botmart, Thongchai, Badriah M. Alotaibi, Rasool Shah, Lamiaa S. El-Sherif, and Samir A. El-Tantawy. 2022. "A Reliable Way to Deal with the Coupled Fractional Korteweg-De Vries Equations within the Caputo Operator" Symmetry 14, no. 11: 2452. https://doi.org/10.3390/sym14112452

APA Style

Botmart, T., Alotaibi, B. M., Shah, R., El-Sherif, L. S., & El-Tantawy, S. A. (2022). A Reliable Way to Deal with the Coupled Fractional Korteweg-De Vries Equations within the Caputo Operator. Symmetry, 14(11), 2452. https://doi.org/10.3390/sym14112452

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop