Next Article in Journal
Riemann–Hilbert Approach for Constructing Analytical Solutions and Conservation Laws of a Local Time-Fractional Nonlinear Schrödinger Type Equation
Next Article in Special Issue
On Conditions for L2-Dissipativity of an Explicit Finite-Difference Scheme for Linearized 2D and 3D Barotropic Gas Dynamics System of Equations with Regularizations
Previous Article in Journal
The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment
Previous Article in Special Issue
Numerical Investigation of Fractional-Order Swift–Hohenberg Equations via a Novel Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform

1
Department of Mechanical Engineering, Sejong University, Seoul 05006, Korea
2
Mechanical Engineering Department, College of Engineering, King Saud University, P.O. Box 800, Al-Riyadh 11421, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(9), 1592; https://doi.org/10.3390/sym13091592
Submission received: 30 June 2021 / Revised: 10 August 2021 / Accepted: 18 August 2021 / Published: 30 August 2021

Abstract

:
In this article, we use the homotopy perturbation transform method to find the fractional Kersten–Krasil’shchik coupled Korteweg–de Vries (KdV) non-linear system. This coupled non-linear system is typically used to describe electric circuits, traffic flow, shallow water waves, elastic media, electrodynamics, etc. The homotopy perturbation method is modified with the help of the ρ -Laplace transformation to investigate the solution of the given examples to show the accuracy of the current technique. The solution of the given technique and the actual results are shown and analyzed with figures.

1. Introduction

Fractional calculus has a long history that began at the end of the 17 century and continued until the 20th century. Many scientists and mathematicians have spent the last few decades researching numerical and analytic results to nonlinear fractional partial differential equations in engineering and science. Due to their use in many implementations, such as viscoelasticity, fluid mechanic, physics, biology, dynamic schemes of control theories, electric networks, optics, chemical physics, and signal processing, fractional differential equations have gained special attention. Many books contain various definitions and basic concepts of fractional calculus [1,2,3,4].
A differential equation symmetry is a transformation that makes the differential equation invariant. The existence of such symmetries may aid in the solution of the differential equation. A scheme of differential equations line symmetry is a continous symmetry of a scheme of differential equations. Solving a linked set of ordinary differential equations can reveal symmetries. It is sometimes easier to solve these equations than it is to solve the original differential equations. The symmetry structure of the system consists of integer partial differential equations and fractional-order partial differential equations with the fractional Caputo derivative.
Using the symmetry construction, we consider two particular cases: pure fractional-order partial differential equations, whose symmetry condition is divided into two parts of integer-order and fractional, and the linear scheme of fractional partial differential equations, which acknowledges an infinite dimensional insignificant generator continuously. Second, we built a theoretical structure of potential symmetry and constructed three potential schemes to analyze potential symmetries of time-fractional partial differential equations with a divergence form using the composition rule of fractional derivatives. Still, when dealing with linear equations, it is impossible to find their exact results.
Therefore, many approximate and numerical methods have been introduced for the results of fractional-order linear and non-linear differential equations, such as variation iteration technique [5,6], Adomian’s decomposition technique [7], homotopy analysis technique [8], homotopy perturbation technique [9], differential transform technique [10], homotopy asymptotic technique [11], and Galerkin technique [12].
The fractional coupled Kersten–Krasil’shchik Korteweg–de Vries (KdV) nonlinear scheme and homogeneous two component fractional coupled third-order KdV schemes are significant fractional-order nonlinear systems for identifying wave behavior elaborating numerous nonlinear phenomena in physics. The fractional-order coupled system is widely used to investigate complex behavior containing multi-components, such as atoms, ions, and free electrons. Many scholars have attempted to study this behavior numerically.
Recently, Paul Kersten and Joseph Krasil’shchik analyzed and modified the KdV equation, proposing absolute complexity between coupled KdV nonlinear systems to analyze the behavior of nonlinear systems. Numerous variations of this Kersten–Krasil’shchik coupled KdV-mKdV nonlinear system have been introduced by many researchers [13,14,15,16]. Many researchers applied different analytical and numerical methods to investigate KdV equations, such as the homotopy analysis method [17], expansion method [18], shifted Legendre polynomials [19], and natural decomposition method [20].
The homotopy perturbation method (HPM) was first proposed by He [21]. For solving differential and integral equations, linear and nonlinear has been the subject abstract and applied analysis of extensive analytical and numerical studies [21]. The HPM is efficient and effective and eliminates an unconditioned matrix, complicated integrals, and infinite series. This method does not need a specific parameter of the model. The ρ -Laplace transformation is a modified transformation of the Laplace transform.
It should be remembered that, with the use of ρ -Laplace transform, absolute differential equations with variable coefficients cannot be solved by Laplace and Sumudu transforms [22,23]. The homotopy perturbation transformation method (HPTM) combines the ρ -Laplace transformation and the homotopy perturbation method. Numerous researchers have utilized HPTM to solve various equations, such as Navier–Stokes problems [24], heat-like problems [25], gas dynamic models [26], and Fisher’s and hyperbolic equation [27].

2. Basic Definitions

In this section, the generalized fractional integral, the generalized fractional derivative, the Mittag–Lefller function, and the ρ -Laplace transformation are defined.
Definition 1.
The fractional-order generalized integral ϱ of a continuous function g : [ 0 , + ] R is given as [22]
( I ϱ , ρ g ) ( χ ) = 1 Γ ( ϱ ) 0 χ χ ρ s ρ ρ ϱ 1 g ( s ) d s s 1 ρ ,
the gamma function denote by Γ, ρ > 0 , χ > 0 and 0 < ϱ < 1 .
Definition 2.
The fractional-order generalized derivative of ϱ of a continuous function g : [ 0 , + ] R is given as [22]
( D ϱ , ρ g ) ( χ ) = ( I 1 ϱ , ρ g ) ( χ ) = 1 Γ ( 1 ϱ ) d d χ 0 χ χ ρ s ρ ρ ϱ g ( s ) d s s 1 ρ .
where we define the gamma function Γ, ρ > 0 , χ > 0 and 0 < ϱ < 1 .
Definition 3.
The fractional-order Caputo derivative ϱ of a continuous function g : [ 0 , + ] R is expressed as [22]
( D ϱ , ρ g ) ( χ ) = 1 Γ ( 1 ϱ ) 0 χ χ ρ s ρ ρ ϱ β g ( s ) d s s 1 ρ .
where ρ > 0 , χ > 0 , β = χ 1 ρ d d χ and 0 < ϱ < 1 .
Definition 4.
The ρ-Laplace transform of a continuous function g : [ 0 , + ] R is defined as [22]
L ρ { g ( χ ) } = 0 e s χ ρ ρ g ( χ ) d χ χ 1 ρ .
A continuous function g as the fractional-order Caputo generalized rho-Laplace transform derivative is defined by
L ρ { D ϱ , ρ g ( χ ) } = s ϱ L ρ { g ( χ ) } k = 0 n 1 s ϱ k 1 ( I ϱ , ρ β n g ) ( 0 ) .

3. The General Methodology of HPTM

In this section, the HPTM for the general form of FPDEs
D ϱ φ ( χ , ) + M φ ( χ , ) + N φ ( χ , ) = h ( χ , ) , > 0 , 0 < ϱ 1 ,
with the initial condition
φ ( χ , 0 ) = g ( χ ) .
where D ϱ φ ( χ , ) = ϱ ϱ is the Caputo fractional derivative of order ϱ ; M and N, are linear and non-linear functions, respectively; and h is the source operator. Using the ρ -Laplace transform of Equation (1),
L ρ [ D ϱ φ ( χ , ) + M φ ( χ , ) + N φ ( χ , ) ] = L ρ [ h ( χ , ) ] , > 0 , 0 < ϱ 1 , φ ( χ , ) = 1 s g ( χ ) + 1 s ϱ L ρ [ h ( χ , ) ] 1 s ϱ L ρ [ M φ ( χ , ) + N φ ( χ , ) ] .
Now, by taking inverse ρ -Laplace transform, we obtain
φ ( χ , ) = L ρ 1 1 s g ( χ ) + 1 s ϱ L ρ [ h ( χ , ) ] L ρ 1 1 s ϱ L ρ { M φ ( χ , ) + N φ ( χ , ) } ,
where
φ ( χ , ) = g ( χ ) + L ρ 1 1 s ϱ L ρ [ h ( χ , ) ] L ρ 1 1 s ϱ L ρ { M φ ( χ , ) + N φ ( χ , ) } ,
Now, the perturbation method parameter p is defined as
φ ( χ , ) = 𝚥 = 0 p 𝚥 φ 𝚥 ( χ , ) ,
where the perturbation term p [ 0 , 1 ] .
The non-linear functions can be defined as
N φ ( χ , ) = 𝚥 = 0 p 𝚥 H 𝚥 ( φ 𝚥 ) ,
where H m are He polynomials of φ 0 , φ 1 , φ 2 , , φ m , and can be determined as
H m ( φ 0 , φ 1 , , φ m ) = 1 m ! m p m N 𝚥 = 0 p 𝚥 φ 𝚥 p = 0 , m = 0 , 1 , 2 .
putting Equations (7) and (8) in Equation (5), we have
𝚥 = 0 p 𝚥 φ 𝚥 ( χ , ) = g ( χ ) + L ρ 1 1 s ϱ L ρ [ h ( χ , ) ] p × L ρ 1 1 s ϱ L ρ { M 𝚥 = 0 p 𝚥 φ 𝚥 ( χ , ) + 𝚥 = 0 p 𝚥 H 𝚥 ( φ 𝚥 ) } .
With both sides comparison coefficient of p, we have
p 0 : φ 0 ( χ , ) = g ( χ ) + L ρ 1 1 s ϱ L ρ [ h ( χ , ) ] , p 1 : φ 1 ( χ , ) = L ρ 1 1 s ϱ L ρ ( M φ 0 ( χ , ) + H 0 ( φ ) ) ,
p 2 : φ 2 ( χ , ) = L ρ 1 1 s ϱ L ρ ( M φ 1 ( χ , ) + H 1 ( φ ) ) , p 𝚥 : φ 𝚥 ( χ , ) = L ρ 1 1 s ϱ L ρ ( M φ 𝚥 1 ( χ , ) + H 𝚥 1 ( φ ) ) , 𝚥 > 0 , 𝚥 N .
φ ( χ , ) = lim M 𝚥 = 1 M φ 𝚥 ( χ , ) .

4. Numerical Experiments

Example 1.
Assume time fractional Kersten–Krasil’shchik coupled KdV-mKdV nonlinear system as:
D ϱ μ + μ 3 χ 6 μ μ χ + 3 φ φ 3 χ + 3 φ χ φ 2 χ 3 μ χ φ 2 + 6 μ φ φ χ = 0 , > 0 , χ R , 0 < ϱ 1 , D ϱ φ + φ 3 χ 3 φ 2 φ χ 3 μ φ χ + 3 μ χ φ = 0 ,
with the initial conditions
μ ( χ , 0 ) = c 2 c sech 2 ( c χ ) , c > 0 , φ ( χ , 0 ) = 2 c sech ( c χ ) .
Using ρ-Laplace transform on Equation (13) by the application of the initial condition given by Equation (14), we find
L ρ [ μ ( χ , ) ] = c 2 c sech 2 ( c χ ) 1 s ϱ L ρ [ μ 3 χ 6 μ μ χ + 3 φ φ 3 χ + 3 φ χ φ 2 χ 3 μ χ φ 2 + 6 μ φ φ χ ] , L ρ [ φ ( χ , ) ] = 2 c sech ( c χ ) 1 s ϱ L ρ [ φ 3 χ 3 φ 2 φ χ 3 μ φ χ + 3 μ χ φ ]
Apply the inverse ρ-Laplace transform, and we obtain
μ ( χ , ) = c 2 c sech 2 ( c χ ) L ρ 1 1 s ϱ L ρ [ μ 3 χ 6 μ μ χ + 3 φ φ 3 χ + 3 φ χ φ 2 χ 3 μ χ φ 2 + 6 μ φ φ χ ] , φ ( χ , ) = 2 c sech ( c χ ) L ρ 1 1 s ϱ L ρ [ φ 3 χ 3 φ 2 φ χ 3 μ φ χ + 3 μ χ φ ]
Use HPM on Equation (16), and we obtain
n = 0 p n μ n ( χ , ) = c 2 c sech 2 ( c χ ) p L ρ 1 1 s ϱ S n = 0 p n μ n ( χ , ) 3 χ + n = 0 p n H n ( μ ) , n = 0 p n φ n ( χ , ) = 2 c sech ( c χ ) p L ρ 1 1 s ϱ S n = 0 p n φ n ( χ , ) 3 χ + n = 0 p n H n ( φ )
Nonlinear steps given by He’s polynomials H n ( μ ) and H n ( φ ) , which are given as
n = 0 p n H n ( μ ) = 6 μ μ χ + 3 φ φ 3 χ + 3 φ χ φ 2 χ 3 μ χ φ 2 + 6 μ φ φ χ , n = 0 p n H n ( φ ) = 3 φ 2 φ χ 3 μ φ χ + 3 μ χ φ .
Values of components of He’s polynomials are given by
H 0 ( μ ) = 6 μ 0 ( μ 0 ) χ + 3 φ 0 ( φ 0 ) 3 χ + 3 ( φ 0 ) χ ( φ 0 ) 2 χ 3 ( μ 0 ) χ ( φ 0 ) 2 + 6 ( μ 0 ) φ 0 ( φ 0 ) χ , H 1 ( μ ) = 6 μ 1 ( μ 0 ) χ 6 μ 0 ( μ 1 ) χ + 3 φ 1 ( φ 0 ) 3 χ + 3 φ 0 ( φ 1 ) 3 χ + 3 ( φ 0 ) χ ( φ 1 ) 2 χ + 3 ( φ 0 ) 2 χ ( φ 1 ) χ 3 ( μ 1 ) χ ( φ 0 ) 2 + 6 ( μ 0 ) χ φ 0 φ 1 + 6 ( μ 0 ) φ 1 ( φ 0 ) χ + 6 ( μ 0 ) φ 0 ( φ 1 ) χ 6 ( μ 1 ) φ 0 ( φ 0 ) χ , H 2 ( μ ) = 6 μ 2 ( μ 0 ) x i 6 μ 1 ( μ 1 ) χ 6 μ 0 ( μ 2 ) χ + 3 φ 2 ( φ 0 ) 3 χ + 3 φ 1 ( φ 1 ) 3 χ + 3 φ 0 ( φ 2 ) 3 χ + 3 ( φ 0 ) χ ( φ 2 ) 2 χ + 3 ( φ 1 ) χ ( φ 1 ) 2 χ + 3 ( φ 2 ) χ ( φ 0 ) 2 χ 3 ( μ 2 ) χ ( φ 0 ) 2 6 ( μ 1 ) χ φ 0 φ 1 6 ( μ 1 ) χ φ 0 φ 2 + 6 ( μ 2 ) φ 0 ( φ 0 ) χ + 6 ( μ 1 ) φ ( φ 0 ) χ + 6 ( μ 2 ) φ 0 ( φ 0 ) χ + 6 ( μ 0 ) φ 1 ( φ 1 ) χ + 6 ( μ 1 ) φ 0 ( φ 1 ) χ + 6 ( μ 0 ) φ 0 ( φ 2 ) χ H 3 ( μ ) = 6 μ 3 ( μ 0 ) χ 6 μ 2 ( μ 1 ) χ 6 μ 1 ( μ 2 ) χ 6 μ 0 ( μ 3 ) χ + 3 φ 3 ( φ 0 ) 3 χ + 3 φ 2 ( φ 1 ) 3 χ + 3 φ 1 ( φ 2 ) 3 χ + 3 φ 0 ( φ 3 ) 3 χ + 3 ( φ 0 ) χ ( φ 3 ) 2 χ + 3 ( φ 1 ) χ ( φ 2 ) 2 χ + 3 ( φ 2 ) χ ( φ 1 ) 2 χ 3 ( φ 3 ) χ ( φ 0 ) 2 χ 3 ( μ 3 ) χ φ 0 2 6 ( μ 2 ) φ 0 φ 1 6 ( μ 1 ) φ 0 φ 2 3 ( μ 1 ) χ φ 1 2 + μ 0 φ 3 ( φ 0 ) χ + 6 ( μ 1 ) φ 2 ( φ 0 ) χ + 6 ( μ 1 ) φ 2 ( φ 0 ) χ + 6 ( μ 3 ) φ 0 ( φ 0 ) χ + 6 ( μ 0 ) φ 2 ( φ 1 ) χ + 6 ( μ 1 ) φ 1 ( φ 1 ) χ + 6 ( μ 2 ) φ 0 ( φ 1 ) χ + 6 ( μ 0 ) φ 1 ( φ 2 ) χ + 6 ( μ 1 ) φ 0 ( φ 2 ) χ + 6 ( μ 0 ) φ 0 ( φ 3 ) χ
and
H 0 ( φ ) = 3 ( φ 0 ) 2 ( φ 0 ) χ 3 μ 0 ( φ 0 ) χ + 3 ( μ 0 ) χ φ 0 . H 1 ( φ ) = 3 ( φ 0 ) 2 ( φ 1 ) χ 6 φ 0 φ 1 ( φ 0 ) χ 3 μ 1 ( φ 0 ) χ 3 μ 0 ( φ 1 ) χ + 3 ( μ 1 ) χ φ 0 + 3 ( μ 0 ) χ φ 1 H 2 ( φ ) = 3 ( φ 0 ) 2 ( φ 2 ) χ 6 φ 0 φ 1 ( φ 1 ) χ 6 φ 0 φ 2 ( φ 1 ) χ 3 ( φ 2 ) 2 ( φ 0 ) χ 3 μ 2 ( φ 0 ) χ 3 μ 1 ( φ 1 ) χ 3 μ 0 ( φ 2 ) χ + 3 ( μ 0 ) χ φ 2 3 ( μ 1 ) χ φ 1 3 ( μ 2 ) χ φ 0 H 3 ( φ ) = 3 ( φ 0 ) 2 ( φ 3 ) χ 6 φ 0 φ 1 ( φ 2 ) χ 6 φ 0 φ 2 ( φ 1 ) χ 6 φ 0 φ 3 ( φ 0 ) χ 6 φ 1 φ 2 ( φ 0 ) χ 3 ( φ 3 ) 2 ( φ 3 ) χ 3 μ 3 ( φ 0 ) χ 3 μ 2 ( φ 1 ) χ 3 μ 1 ( φ 2 ) χ 3 μ 0 ( φ 3 ) χ + 3 ( μ 0 ) χ φ 3 + 3 ( μ 1 ) χ φ 2 + 3 ( μ 2 ) χ φ 1 + 3 ( μ 3 ) χ φ 0
Comparing the coefficients of same powers of p, we have
p 0 : μ 0 ( χ , ) = c 2 c sech 2 ( c χ ) , p 1 : μ 1 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 0 ) 3 χ + H 0 ( μ ) ] = 8 c 5 2 sinh ( c χ ) sech 3 ( c χ ) ρ ρ ϱ Γ ( ϱ + 1 ) ,
p 2 : μ 2 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 1 ) 3 χ + H 1 ( μ ) ] = 16 c 4 [ 2 cosh 2 ( c χ ) 3 ] sech 4 ( c χ ) ρ ρ 2 ϱ Γ ( 2 ϱ + 1 ) , p 3 : μ 3 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 2 ) 3 χ + H 2 ( μ ) ] = 128 c 11 / 2 [ cosh 2 ( c χ ) 3 ] sinh ( c χ ) sech 5 ( c χ ) ρ ρ 3 ϱ Γ ( 3 ϱ + 1 ) , p 4 : μ 4 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 3 ) 3 χ + H 3 ( μ ) ] = 256 c 7 [ 2 cosh 4 ( c χ ) 15 cosh 2 ( c χ ) + 15 ] sech 6 ( c χ ) ρ ρ 4 ϱ Γ ( 4 ϱ + 1 ) , .
and
p 0 : φ 0 ( χ , ) = 2 c sech ( c χ ) , p 1 : φ 1 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( φ 0 ) 3 χ + H 0 ( φ ) ] , = 4 c 2 sinh ( c χ ) sech 2 ( c χ ) ρ ρ ϱ Γ ( ϱ + 1 ) , p 2 : φ 2 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( φ 1 ) 3 χ + H 1 ( φ ) ] , = 8 c 7 2 [ cosh 2 ( c χ ) 2 ] sech 3 ( c χ ) ρ ρ 2 ϱ Γ ( 2 ϱ + 1 ) , p 3 : φ 3 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( φ 2 ) 3 χ + H 2 ( φ ) ] , = 16 c 2 [ cosh 5 ( c χ ) 6 ] sinh ( c χ ) sech 4 ( c χ ) ρ ρ 3 ϱ Γ ( 3 ϱ + 1 ) , p 4 : φ 4 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( φ 3 ) 3 χ + H 3 ( φ ) ] , = 32 c 13 / 2 [ cosh 4 ( c χ ) 20 cosh 2 ( c χ ) + 24 ] sech 5 ( c χ ) ρ ρ 4 ϱ Γ ( 4 ϱ + 1 ) ,
Hence, the series solution is given by
μ ( χ , ) = n = 0 μ i ( χ , ) = 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 ) + 128 c 11 / 2 [ cosh 2 ( c χ ) 3 ] sinh ( c χ ) sech 5 ( c χ ) ρ ρ 3 ϱ Γ ( 3 ϱ + 1 ) 256 c 7 [ 2 cosh 4 ( c χ ) 15 cosh 2 ( c χ ) + 15 ] sech 6 ( c χ ) ρ ρ 4 ϱ Γ ( 4 ϱ + 1 ) + ,
and
φ ( χ , ) = n = 0 φ i ( χ , ) = 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 ) 16 c 2 [ cosh 5 ( c χ ) 6 ] sinh ( c χ ) sech 4 ( c χ ) ρ ρ 3 ϱ Γ ( 3 ϱ + 1 ) + 32 c 13 / 2 [ cosh 4 ( c χ ) 20 cosh 2 ( c χ ) + 24 ] sech 5 ( c χ ) ρ ρ 4 ϱ Γ ( 4 ϱ + 1 ) ,
Putting ϱ a n d ρ = 1 in Equations (24) and (25), we obtain the solution of the problem as:
μ ( χ , ) = c 2 c sech 2 ( c χ ) + 8 c 5 2 sinh ( c χ ) sech 3 ( c χ ) 8 2 c 4 [ 2 cosh 2 ( c χ ) 3 ] sech 4 ( c χ ) + 64 3 3 c 11 / 2 [ cosh 2 ( c χ ) 3 ] sinh ( c χ ) sech 5 ( c χ ) 32 3 4 c 7 [ 2 cosh 4 ( c χ ) 15 cosh 2 ( c χ ) + 15 ] sech 6 ( c χ ) + ,
and
φ ( χ , ) = 2 c sech ( c χ ) 4 c 2 sinh ( c χ ) sech 2 ( c χ ) + 4 2 c 7 2 [ cosh 2 ( c χ ) 2 ] sech 3 ( c χ ) 8 3 3 c 2 [ cosh 5 ( c χ ) 6 ] sinh ( c χ ) sech 4 ( c χ ) + 4 3 4 c 13 / 2 [ cosh 4 ( c χ ) 20 cosh 2 ( c χ ) + 24 ] sech 5 ( c χ ) ,
The solution represented by Equations (26) and (27) is similar to exact solution in closed form as:
μ ( χ , ) = c 2 c sech 2 ( c ( χ + 2 c ) ) , φ ( χ , ) = 2 c sech ( c ( χ + 2 c ) ) .
In Figure 1, the actual and HPTM solutions of μ ( χ , ) is calculated at ϱ = 1 . In Figure 2, the 3D graphs for μ ( χ , ) for different fractional-order shows that the HPTM approximated solutions derived are in a strong agreement with the actual and the approximate solution. Similarly Figure 3, the actual and HPTM solutions of φ ( χ , ) is calculated at ϱ = 1 . In Figure 4, the φ ( χ , ) for different fractional-order shows that the HPTM approximated solutions derived are in a strong agreement with the actual and the approximate solution. This comparison shows that the HPTM and the actual solutions are very close. As a result, the HPTM is a dependable new study that requires less computation of computations, is adaptable, and simple to use.
Example 2.
Assume a homogeneous two-component time fractional coupled third order KdV system as:
D ϱ μ μ 3 χ μ μ χ φ φ χ = 0 , > 0 , χ R , 0 < ϱ 1 , D ϱ φ + 2 φ 3 χ μ φ χ = 0 ,
with the initial condition
μ ( χ , 0 ) = 3 6 tanh 2 χ 2 , φ ( χ , 0 ) = 3 c 2 tanh χ 2 .
Using the ρ-Laplace transform on Equation (29) by the application of initial conditions given by Equation (30), we obtain
L ρ [ μ ( χ , ) ] = 3 6 tanh 2 χ 2 + 1 s ϱ L ρ [ μ 3 χ μ μ χ φ φ χ ] , L ρ [ φ ( χ , ) ] = 3 c 2 tanh χ 2 1 s ϱ L ρ [ 2 φ 3 χ μ φ χ ] .
Applying the inverse ρ-Laplace transform, we obtain
μ ( χ , ) = 3 6 tanh 2 χ 2 + L ρ 1 1 s ϱ L ρ [ μ 3 χ μ μ χ φ φ χ , φ ( χ , ) = 3 c 2 tanh χ 2 L ρ 1 1 s ϱ L ρ [ 2 φ 3 χ μ φ χ ] .
Using HPM on Equation (32), we obtain
n = 0 p n μ n ( χ , ) = 3 6 tanh 2 χ 2 + p L ρ 1 1 s ϱ L ρ ( n = 0 p n μ n ( χ , ) ) 3 χ + ( n = 0 p n H n ( μ ) ) ,
n = 0 p n φ n ( χ , ) = 3 c 2 tanh χ 2 p L ρ 1 1 s ϱ L ρ 2 n = 0 p n φ n ( χ , ) 3 χ n = 0 p n H n ( φ ) .
Nonlinear steps given by He’s polynomials H n ( μ ) and H n ( φ ) , are given as:
n = 0 p n H n ( μ ) = μ μ χ + φ φ χ , n = 0 p n H n ( φ ) = μ φ χ .
Values of factors of He’s polynomials are given as
H 0 ( μ ) = μ 0 ( μ 0 ) χ + φ 0 ( φ 0 ) χ , H 1 ( μ ) = μ 1 ( μ 0 ) χ + μ 0 ( μ 1 ) χ + φ 1 ( φ 0 ) χ + φ 0 ( φ 1 ) χ , H 2 ( μ ) = μ 2 ( μ 0 ) χ + μ 1 ( μ 1 ) χ + μ 0 ( μ 2 ) χ + φ 2 ( φ 0 ) χ + φ 1 ( φ 1 ) χ + φ 0 ( φ 2 ) χ , H 3 ( μ ) = μ 3 ( μ 0 ) χ + μ 2 ( μ 1 ) χ + μ 1 ( μ 2 ) χ + μ 0 ( μ 3 ) χ + φ 3 ( φ 0 ) χ + φ 2 ( φ 1 ) χ + φ 1 ( φ 2 ) χ + φ 0 ( φ 3 ) χ ,
and
H 0 ( φ ) = μ 0 ( φ 0 ) χ , H 1 ( φ ) = μ 1 ( φ 0 ) χ μ 0 ( φ 1 ) χ , H 2 ( φ ) = μ 2 ( φ 0 ) χ μ 1 ( φ 1 ) χ μ 0 ( φ 2 ) χ , H 3 ( φ ) = μ 3 ( φ 0 ) χ μ 2 ( φ 1 ) χ μ 1 ( φ 2 ) χ μ 0 ( φ 3 ) χ ,
Comparing coefficients of the same powers of p, we have
p 0 : μ 0 ( χ , ) = 3 6 tanh 2 χ 2 , p 1 : μ 1 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 0 ) 3 χ + H 0 ( μ ) ] , = 6 sech 2 χ 2 tanh χ 2 ρ ρ ϱ Γ ( ϱ + 1 ) , p 2 : μ 2 ( χ , ) = L ρ 1 1 s ϱ L ρ ( μ 1 ) 3 χ + H 1 ( μ ) , = 3 [ 2 + 7 sech 2 χ 2 15 sech 4 χ 2 ] sech 2 χ 2 ρ ρ 2 ϱ Γ ( 2 ϱ + 1 ) , p 3 : μ 3 ( χ , ) = L ρ 1 1 s ϱ L ρ [ ( μ 2 ) 3 χ + H 2 ( μ ) ] ,
and
p 0 : φ 0 ( χ , ) = 3 c 2 tanh χ 2 , p 1 : φ 1 ( χ , ) = L ρ 1 1 s ϱ L ρ [ 2 ( φ 0 ) 3 χ H 0 ( φ ) ] , = 3 c 2 sech 2 χ 2 tanh χ 2 ρ ρ ϱ Γ ( ϱ + 1 ) , p 2 : φ 2 ( χ , ) = L ρ 1 1 s ϱ L ρ [ 2 ( φ 1 ) 3 χ H 1 ( φ ) ] , = 3 c 2 2 2 + 21 sech 2 χ 2 24 sech 4 χ 2 sech 2 χ 2 ρ ρ 2 ϱ Γ ( 2 ϱ + 1 ) , p 3 : φ 3 ( χ , ) = L ρ 1 1 s ϱ L ρ [ 2 ( φ 2 ) 3 χ H 2 ( φ ) ] ,
Hence, the series solution is given by
μ ( χ , ) = n = 0 μ n ( χ , ) , = 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 ) , φ ( χ , ) = n φ n ( χ , ) = 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 ) + .
Placing ϱ a n d ρ = 1 in Equation (40), we obtain the solution of the problem as:
μ ( χ , ) = 3 6 tanh 2 χ 2 + 6 sech 2 χ 2 tanh χ 2 + 3 2 2 [ 2 + 7 sech 2 χ 2 15 sech 4 χ 2 ] sech 2 χ 2 , φ ( χ , ) = 3 c 2 tanh χ 2 + 3 c 2 sech 2 χ 2 tanh χ 2 + 3 c 2 4 2 2 + 21 sech 2 χ 2 24 sech 4 χ 2 sech 2 χ 2 + .
The solution given by Equation (41) is similar to the closed form solution as:
μ ( χ , ) = 3 6 tanh 2 ( χ + 2 ) , φ ( χ , ) = 3 c 2 tanh ( χ + 2 ) .
In Figure 5 the actual and HPTM solutions of μ ( χ , ) is calculated at ϱ = 1 . In Figure 6, the 3D graphs for μ ( χ , ) for different fractional-order shows that the HPTM approximated solutions derived are in a strong agreement with the actual and the approximate solution. Similarly Figure 7, the actual and HPTM solutions of φ ( χ , ) is calculated at ϱ = 1 . In Figure 8, the φ ( χ , ) for different fractional-order shows that the HPTM approximated solutions derived are in a strong agreement with the actual and the approximate solution. This comparison shows that the HPTM and the actual solutions are very close. As a result, the HPTM is a dependable new study that requires less computation of computations, is adaptable, and simple to use.

5. Conclusions

In this paper, we calculated the fractional-order Kersten–Krasil’shchik coupled KdV-mKdV nonlinear system, using an ρ -Laplace transform. The suggested method was applied to obtain the solution of the given two problems. The HPTM solution is in close contact with the exact result of the given problems. We also calculated the results of the given problems with the fractional-order derivatives. The figures of the fractional-order results achieved demonstrated convergence toward the results of the integer-order. Furthermore, the present method is simple, straightforward, and requires less computational cost; the current technique can be modified to solve other fractional-order partial differential equations.

Author Contributions

Conceptualization, N.A.S. and A.H.S.; methodology, N.A.S.; software, A.H.S.; validation, J.D.C., N.A.S. and A.H.S.; formal analysis, J.D.C.; investigation, N.A.S.; resources, A.H.S.; data curation, J.D.C.; writing—original draft preparation, N.A.S.; writing—review and editing, A.H.S.; visualization, J.D.C.; supervision, J.D.C.; project administration, A.H.S.; funding acquisition, J.D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Researchers Supporting Project number (RSP-2021/373), King Saud University, Riyadh, Saudi Arabia. This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2017R1D1A1B05030422).

Acknowledgments

The authors would like to acknowledge the Researchers Supporting Project number (RSP-2021/373), King Saud University, Riyadh, Saudi Arabia and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2017R1D1A1B05030422).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zaslavsky, G. Book Review: “Theory and Applications of Fractional Differential Equations” by Anatoly A. Kilbas, Hari M. Srivastava and Juan J. Trujillo. Fractals 2007, 15, 101–102. [Google Scholar] [CrossRef]
  2. Kilbas, A. Partial fractional differential equations and some of their applications. Analysis 2010, 30. [Google Scholar] [CrossRef]
  3. Mainardi, F. Fractional Calculus: Theory and Applications. Mathematics 2018, 6, 145. [Google Scholar] [CrossRef] [Green Version]
  4. Mu’lla, M. Fractional Calculus, Fractional Differential Equations and Applications. Oalib 2020, 7, 1–9. [Google Scholar] [CrossRef]
  5. Shah, R.; Khan, H.; Baleanu, D.; Kumam, P.; Arif, M. A semi-analytical method to solve family of Kuramoto–Sivashinsky equations. J. Taibah Univ. Sci. 2020, 14, 402–411. [Google Scholar] [CrossRef]
  6. Sakar, M.; Erdogan, F.; Yildirim, A. Variational iteration method for the time-fractional Fornberg-Whitham equation. Comput. Math. Appl. 2012, 63, 1382–1388. [Google Scholar] [CrossRef] [Green Version]
  7. El-Sayed, A.; Behiry, S.; Raslan, W. Adomian’s decomposition method for solving an intermediate fractional advection-dispersion equation. Comput. Math. Appl. 2010, 59, 1759–1765. [Google Scholar] [CrossRef] [Green Version]
  8. Nonlaopon, K.; Alsharif, A.; Zidan, A.; Khan, A.; Hamed, Y.; Shah, R. Numerical Investigation of Fractional-Order Swift-Hohenberg Equations via a Novel Transform. Symmetry 2021, 13, 1263. [Google Scholar] [CrossRef]
  9. Zhang, X.; Zhao, J.; Liu, J.; Tang, B. Homotopy perturbation method for two dimensional time-fractional wave equation. Appl. Math. Model. 2014, 38, 5545–5552. [Google Scholar] [CrossRef]
  10. Srivastava, V.; Awasthi, M.; Kumar, S. Analytical approximations of two and three dimensional time-fractional telegraphic equation by reduced differential transform method. Egypt. J. Basic Appl. Sci. 2014, 1, 60–66. [Google Scholar] [CrossRef] [Green Version]
  11. Pandey, R.; Singh, O.; Baranwal, V. An analytic algorithm for the space-time fractional advection-dispersion equation. Comput. Phys. Commun. 2011, 182, 1134–1144. [Google Scholar] [CrossRef]
  12. Wei, L.; He, Y. Analysis of a fully discrete local discontinuous Galerkin method for time-fractional fourth-order problems. Appl. Math. Model. 2014, 38, 1511–1522. [Google Scholar] [CrossRef]
  13. Rui, W.; Qi, X. Bilinear approach to quasi-periodic wave solutions of the Kersten–Krasil’shchik coupled KdV-mKdV system. Bound. Value Probl. 2016, 2016, 1–13. [Google Scholar] [CrossRef] [Green Version]
  14. Keskin, Y.; Oturanc, G. Reduced Differential Transform Method for Partial Differential Equations. Int. J. Nonlinear Sci. Numer. Simul. 2009, 10, 741–750. [Google Scholar] [CrossRef]
  15. Kalkanli, A.; Sakovich, S.; 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]
  16. Goswami, A.; Sushila; Singh, J.; Kumar, D. Numerical computation of fractional Kersten–Krasil’shchik coupled KdV-mKdV system occurring in multi-component plasmas. AIMS Math. 2020, 5, 2346–2368. [Google Scholar]
  17. Jafari, H.; Prasad, J.; Goswami, P.; Dubey, R. Solution of the Local Fractional Generalized KDV Equation Using Homotopy Analysis Method. Fractals 2021, 29, 2140014. [Google Scholar] [CrossRef]
  18. Yang, Y.; Qi, J.; Tang, X.; Gu, Y. Further Results about Traveling Wave Exact Solutions of the (2 + 1)-Dimensional Modified KdV Equation. Adv. Math. Phys. 2019, 2019, 1–10. [Google Scholar] [CrossRef]
  19. Bhrawy, A.; Doha, E.; Ezz-Eldien, S.; Abdelkawy, M. A numerical technique based on the shifted Legendre polynomials for solving the time-fractional coupled KdV equations. Calcolo 2015, 53, 1–17. [Google Scholar] [CrossRef]
  20. Elbadri, M.; Ahmed, S.; Abdalla, Y.; Hdidi, W. A New Solution of Time-Fractional Coupled KdV Equation by Using Natural Decomposition Method. Abstr. Appl. Anal. 2020, 2020, 1–9. [Google Scholar] [CrossRef]
  21. He, J.H. Homotopy perturbation method: A new nonlinear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar] [CrossRef]
  22. Bhangale, N.; Kachhia, K.B.; Gomez-Aguilar, J.F. A new iterative method with ρ-Laplace transform for solving fractional differential equations with Caputo generalized fractional derivative. Eng. Comput. 2020, 1–14. [Google Scholar] [CrossRef]
  23. Wang, C. Hyers-Ulam-Rassias Stability of the Generalized Fractional Systems and the ρ-Laplace Transform Method. Mediterr. J. Math. 2021, 18, 1–22. [Google Scholar] [CrossRef]
  24. Jena, R.; Chakraverty, S. Solving time-fractional Navier-Stokes equations using homotopy perturbation Elzaki transform. SN Appl. Sci. 2018, 1, 16. [Google Scholar] [CrossRef] [Green Version]
  25. Mahgoub, M.; Sedeeg, A. A Comparative Study for Solving Nonlinear Fractional Heat-Like Equations via Elzaki Transform. Br. J. Math. Comput. Sci. 2016, 19, 1–12. [Google Scholar] [CrossRef]
  26. Das, S.; Gupta, P. An Approximate Analytical Solution of the Fractional Diffusion Equation with Absorbent Term and External Force by Homotopy Perturbation Method. Z. Fur Naturforschung A 2010, 65, 182–190. [Google Scholar] [CrossRef]
  27. Jan, R.; Khan, H.; Kumam, P.; Tchier, F.; Shah, R.; Bin Jebreen, H. The Investigation of the Fractional-View Dynamics of Helmholtz Equations within Caputo Operator. Comput. Mater. Contin. 2021, 68, 3185–3201. [Google Scholar] [CrossRef]
Figure 1. The graphs of actual and HPTM results at μ ( χ , ) of Example 1.
Figure 1. The graphs of actual and HPTM results at μ ( χ , ) of Example 1.
Symmetry 13 01592 g001
Figure 2. The various fractional-order graphs of μ ( χ , ) for Example 1.
Figure 2. The various fractional-order graphs of μ ( χ , ) for Example 1.
Symmetry 13 01592 g002
Figure 3. The actual and HPTM results graphs at φ ( χ , ) of Example 1.
Figure 3. The actual and HPTM results graphs at φ ( χ , ) of Example 1.
Symmetry 13 01592 g003
Figure 4. The various fractional-order of φ ( χ , ) for Example 1.
Figure 4. The various fractional-order of φ ( χ , ) for Example 1.
Symmetry 13 01592 g004
Figure 5. The actual and HPTM results graphs at μ ( χ , ) of Example 2 for ϱ = 1 .
Figure 5. The actual and HPTM results graphs at μ ( χ , ) of Example 2 for ϱ = 1 .
Symmetry 13 01592 g005
Figure 6. The various fractional-order graphs of μ ( χ , ) for Example 2.
Figure 6. The various fractional-order graphs of μ ( χ , ) for Example 2.
Symmetry 13 01592 g006
Figure 7. The actual and HPTM results graphs at φ ( χ , ) and φ ( χ , ) of Example 2 at ϱ = 1 .
Figure 7. The actual and HPTM results graphs at φ ( χ , ) and φ ( χ , ) of Example 2 at ϱ = 1 .
Symmetry 13 01592 g007
Figure 8. The various fractional-order graphs of φ ( χ , ) for Example 2.
Figure 8. The various fractional-order graphs of φ ( χ , ) for Example 2.
Symmetry 13 01592 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shah, N.A.; Seikh, A.H.; Chung, J.D. The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform. Symmetry 2021, 13, 1592. https://doi.org/10.3390/sym13091592

AMA Style

Shah NA, Seikh AH, Chung JD. The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform. Symmetry. 2021; 13(9):1592. https://doi.org/10.3390/sym13091592

Chicago/Turabian Style

Shah, Nehad Ali, Asiful H. Seikh, and Jae Dong Chung. 2021. "The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform" Symmetry 13, no. 9: 1592. https://doi.org/10.3390/sym13091592

APA Style

Shah, N. A., Seikh, A. H., & Chung, J. D. (2021). The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform. Symmetry, 13(9), 1592. https://doi.org/10.3390/sym13091592

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