Next Article in Journal
Adaptive Fractional-Order Multi-Scale Optimization TV-L1 Optical Flow Algorithm
Next Article in Special Issue
Impressive Exact Solitons to the Space-Time Fractional Mathematical Physics Model via an Effective Method
Previous Article in Journal
The Soliton Solutions for Nonlocal Multi-Component Higher-Order Gerdjikov–Ivanov Equation via Riemann–Hilbert Problem
Previous Article in Special Issue
The Finite Difference Method and Analysis for Simulating the Unsteady Generalized Maxwell Fluid with a Multi-Term Time Fractional Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System

by
Obaid J. Algahtani
Department of Mathematics, College of Sciences, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
Fractal Fract. 2024, 8(3), 178; https://doi.org/10.3390/fractalfract8030178
Submission received: 26 February 2024 / Revised: 14 March 2024 / Accepted: 16 March 2024 / Published: 21 March 2024

Abstract

:
This article aims to examine the nonlinear excitations in a coupled Hirota system described by the fractal fractional order derivative. By using the Laplace transform with Adomian decomposition (LADM), the numerical solution for the considered system is derived. It has been shown that the suggested technique offers a systematic and effective method to solve complex nonlinear systems. Employing the Banach contraction theorem, it is confirmed that the LADM leads to a convergent solution. The numerical analysis of the solutions demonstrates the confinement of the carrier wave and the presence of confined wave packets. The dispersion nonlinear parameter reduction equally influences the wave amplitude and spatial width. The localized internal oscillations in the solitary waves decreased the wave collapsing effect at comparatively small dispersion. Furthermore, it is also shown that the amplitude of the solitary wave solution increases by reducing the fractal derivative. It is evident that decreasing the order α modifies the nature of the solitary wave solutions and marginally decreases the amplitude. The numerical and approximation solutions correspond effectively for specific values of time (t). However, when the fractal or fractional derivative is set to one by increasing time, the wave amplitude increases. The absolute error analysis between the obtained series solutions and the accurate solutions are also presented.

1. Introduction

In recent decades, due to the broad applications of soliton theory in physics, mathematics, and other areas of engineering and applied sciences, the analysis of explicit accurate solutions in the form of solitary wave solutions of evolution equations has played a significant role [1,2,3]. In soliton theory, several analytical methods can be applied to calculate the approximate solution to nonlinear partial differential equations [4,5,6]. To investigate fractal fractional nonlinear coupled Hirota systems, we consider the Hirota equation [7]
ψ ( x , t ) t + 3 δ | ψ ( x , t ) | 2 ψ ( x , t ) x + γ 3 ψ ( x , t ) x 3 = 0 , t > 0 , x [ , ] ,
where the function ψ ( x , t ) is a complex-valued function, and δ , γ R + . The Hirota equation is a modified nonlinear Schrödinger equation (NLSE) that considers higher-order dispersion and time-delay corrections to the cubic nonlinearity [8]. It can be observed as a generalization for the NLS equation when the higher dispersion and time-delay changes are considered, as well as complex generalization of the mKdV equation [9]. It is also known as soliton equation because this usually has an exact solution where localized moving excitations that preserve their shape in the evolution, similarly to the NLS soliton [10,11]. The Hirota equation was initially studied to illustrate the ultra-short pulses experienced from the self-steepening effect and high-order dispersion [8]. It has numerous applications, such as localized wave structures, transmission of optical pulses in optical waveguide arrays, rogue waves, solitary waves, and breathers [12,13]. It has also been widely studied for rational solitary and rogue wave solutions, interaction properties of complex modified solitons, multi-solitons, dark-bright-rogue wave phenomena, and rogue wave pairs [14,15,16,17]. The Hirota equation has been analytically investigated using different methods containing the Hamiltonian formalism, perturbation theory, Darboux transformation, sine–cosine, and tangent methods [18,19,20,21]. The Hirota Equation (1) possess more than one parameter solutions in the form solitons, analogous to NLS equation subjected to velocity and amplitude as
ψ ( x , t ) = A sech η ( x μ t ) e i k x ω t ,
where A = k 2 γ / δ is the amplitude and μ = γ ( η 2 3 k 2 ) represents the wave speed. The frequency of the phase is given by ω = γ ( 3 η 2 k 2 ) , where γ represents nonlinearity dispersion coefficients, η is the wave envelope’s width, and k is the propagation number. By using assumption ψ = ϕ ( x , t ) + i φ ( x , t ) , the Hirota equation can be transformed to coupled Hirota system [22,23]
ϕ t + 3 δ ϕ 2 + φ 2 ϕ x + γ 3 ϕ x 3 = 0 ,
φ t + 3 δ ϕ 2 + φ 2 φ x + γ 3 φ x 3 = 0 .
It should be noted that converting the Hirota equation to a coupled Hirota system allows for the study of growth or decay in the system (the real and imaginary parts of the function ψ ( x , t ) ) separately, leading to deeper insights into the solitary wave solutions. Similar to Equation (1), the coupled Hirota system has also been investigated to describe propagation of electromagnetic pulses in optical fibres by employing the inverse scattering transform [24]. Several other interesting results have been examined for coupled Hirota system, such as relation amongst dark–bright solitons, Lax pair, Painlevé analyses and rogue waves solutions [25,26,27]. Here, we use the concept of fractional fractal derivative of order ( α , β ) to study the coupled Hirota system
D t α , β F F P ϕ + 3 δ ϕ 2 + φ 2 ϕ x + γ 3 ϕ x 3 = 0 ,
D t α , β F F P φ + 3 δ ϕ 2 + φ 2 φ x + γ 3 φ x 3 = 0 ,
where the classical time partial derivatives are replaced with fractional derivatives of the functions of order α ( 0 < α 1 ) in the Caputo sense with power law kernel and fractal derivative β ( 0 < β 1 ). The initial conditions considered herein are
ϕ ( x , 0 ) = f ( x ) , φ ( x , 0 ) = g ( x ) .
The concept of fractional calculus (FC) originated from the challenge that even the concept of a derivative of a classical order could potentially be expanded to stay authentic when the order was not an integer. Following this exceptional theory, the subject of FC caught great attention in mathematics and other areas of scientific research. Since then, FC has been extensively used to examine various physical facts, such as electromagnetic, viscoelastic, and damping theories, artificial intelligence, fluid dynamics, wave propagation theories, chaotic dynamics, heat transfer analysis, appliances, control systems, and various inherent algorithms [28,29,30,31,32,33,34].
In FC, the classical integrals and derivatives are expanded to fractional order as the standard operators are inadequate for understanding many complex systems, because fractional operators produce more realistic and precise results. Numerous fractional derivatives have been derived with distinct kinds of kernels, such as Hilfer, Riemann–Liouville (RL), Caputo, Grünwald, Caputo–Fabrizio, and Atangana–Baleanu derivatives [35,36]. Similarly, several novel kinds of fractional operators have proposed to link the ideas of Caputo and other fractional operators [37,38].
In addition to the fractional-order operators, another new theory of differentiation has been suggested, where the derivative is represented as the fractional order with fractal dimension [39]. The fundamental concept is called fractal fractional integrals and derivatives [40]. If the system under investigation is derivable, the fractal derivative becomes β t β 1 , and the classical operator is regained when the fractal order reaches one [41]. In a fractal derivative, the variable is ascended, corresponding to t α . This modified derivative was proposed to model certain real-world physical systems for which conventional physical phenomena are not appropriate and cannot be applied to the non-integral means of fractal dimension [42,43]. The theory of fractal fractional integrals and derivatives has been extensively used in engineering applications (fluid motion, aquifer, and crystallization), chaotic processes, and electronic networks [44,45]. Here, we will investigate the localized excitations in a coupled Hirota system with time fractional derive in Caputo’s sense with fractal dimension. The Laplace transform with Adomian decomposition is a very effective analytical technique, and has been successfully applied to investigate numerous mathematical problems in classical as well as in fractional calculus. We will derive the series solution of the considered system in a systematic manner using LADM.
The article is organized as follows: In Section 2, various fundamental definitions of fractional calculus are included. In Section 3, the general series solution of the coupled Hirota system with order ( α , β ) is calculated using LADM. The convergence analysis of LADM is examined in Section 4. For confirmation of the obtained findings, a particular example is studied in Section 5 with absolute error analysis. The systematic approximate solutions are compared with numerical results with comprehensive discussions in Section 6.

2. Preliminaries

Here, we state fundamental definitions associated with fractal fractional calculus and the Laplace and inverse Laplace transforms that will be applied in the next section.
Definition 1. 
Let ψ H 1 , and α ( 0 , 1 ] ; we define the Caputo’s derivative as
D t α a C ψ ( t ) = I t ( n α ) a d n d t n ψ ( t ) = 1 Γ ( p α ) a t ( t s ) p α 1 ψ ( s ) d s , p 1 < α p .
Definition 2. 
Let ψ ( t ) be an open differentiable interval ( a , b ) . If ψ ( t ) is fractal fractional differentiable on ( a , b ) , with fractal order β, then the fractal fractional derivative of order α in Caputo sense with power law kernel is given by [40]
D t α , β a F F P ψ ( t ) = 1 Γ ( p α ) a t ( t s ) p α 1 d d s β ψ ( s ) d s , < p 1 < α p , < p 1 < β p ,
where d ψ ( s ) d s β = lim t s ψ ( t ) ψ ( s ) t β s β . The above operator can also be represented as
D t α , β a F F P ψ ( t ) = 1 Γ ( p α ) a t ( t s ) p α 1 d γ d s β ψ ( s ) d s , p 1 < α , β , p ,
where d γ ψ ( s ) d s β = lim t s ψ γ ( t ) ψ γ ( s ) t β s β .
Definition 3. 
Let ψ ( t ) be a continuous function on the interval I, then the fractional integral of ψ ( t ) with fractional order α can be described [40]
I t α a F ψ ( t ) = α 0 t s α 1 ψ ( s ) d s .
From the above, the fractal fractional integral can be obtained as [40]
I t α t F F P ψ ( t ) = β ( α 1 ) ! 0 t s α 1 ( t s ) α 1 ψ ( s ) d s .
Definition 4. 
For a function ψ ( t ) , the Laplace transform is symbolized by L and is defined as [46]
L [ ψ ( t ) ] = F ( s ) = 0 e s t ψ ( s ) d s , t > 0 .
Definition 5. 
If ψ ( t ) = L 1 F ( s ) , then the inverse Laplace transform L 1 is defined as [46]
L 1 F ( s ) s = 0 t ψ ( t ) d t .
Let ψ ( t ) be a continuous function; the fractal Laplace transform of order α can be written as [40]
L a α F ψ ( t ) = 0 t α 1 e a t ψ ( t ) d t .
Definition 6. 
For the function ψ ( x , t ) , applying Laplace transform to Caputo fractional derivative, we obtain
L a C D t α ψ ( x , t ) = s α L ψ ( x , t ) k = 0 n 1 s α k 1 ψ k t ( x , 0 ) , n = [ α ] + 1 .

3. General Solution of Fractal Fractional Coupled Hirota System

Here, we discuss properties of Laplace transform with order ( α , β ) and calculate coupled the Hirota system with fractal fractional derivative of order α using the Laplace Adomian decomposition method. In 1980, George Adomian established a novel approach to solving nonlinear functional equations. Since then, this technique has been called the Adomian decomposition method (ADM), and has been applied in numerous studies [47,48]. The ADM involves breaking down the equation under consideration into linear and nonlinear parts and generates the series solution, whose terms are determined by a recursive relation using the Adomian polynomials. Numerous fundamental studies on various aspects of the extended version of Adomian’s decomposition method have been performed [49,50,51]. Similarly, the extended version of the Laplace decomposition method has been used for the solution of different PDEs [52].
Suppose ψ ( t ) H 1 is a continuous function ∀ t [ 0 , T ] , such that
D t α , β F F P ψ ( t ) = F ( x ) .
Using the definition presented in [45], we can write
D t α C ψ ( t ) = β t β 1 F ( x ) .
Implementing the Laplace transformation to Equation (9) yields
L [ ψ ( t ) ] = ψ ( 0 ) s + 1 s α L [ β t β 1 F ( x ) ] , ψ ( t ) = L 1 ψ ( 0 ) s + 1 s α L ( β t β 1 F ( x ) ) , applying L 1 ψ ( t ) = ψ ( 0 ) + L 1 Γ ( β + 1 ) s α + β F ( x ) = ψ ( 0 ) + Γ ( β + 1 ) t α + β 1 Γ ( α + β ) F ( x ) .
This is the Laplace transformation of fractal fractional derivative with power law kernel.
To obtain the general solution of system (5) and (6), one can apply definition of fractal derivative, which gives
D t α C ϕ = β t β 1 3 δ ϕ 2 ϕ x 3 δ φ 2 ϕ x γ ϕ x x x ,
D t α C φ = β t β 1 3 δ ϕ 2 φ x 3 δ φ 2 φ x γ φ x x x .
Applying the Laplace transformation
L [ C D t α ϕ ] = L β t β 1 ( 3 δ ϕ 2 ϕ x 3 δ φ 2 ϕ x γ ϕ x x x ) , L [ C D t α φ ] = L β t β 1 ( 3 δ ϕ 2 φ x 3 δ φ 2 φ x γ φ x x x ) .
Using the definition above gives
L [ ϕ ] = f ( x ) s + 1 s α L β t β 1 ( 3 δ ϕ 2 ϕ x 3 δ φ 2 ϕ x γ ϕ x x x ) ,
L [ φ ] = g ( x ) s + 1 s α L β t β 1 ( 3 δ ϕ 2 φ x 3 δ φ 2 φ x γ φ x x x ) .
Now, consider ϕ and φ in the series form
ϕ = n = 0 ϕ n , φ = n = 0 φ n ,
the non-linear terms are decomposed as
ϕ 2 ϕ x = n = 0 A n , φ 2 ϕ x = n = 0 B n a n d ϕ 2 φ x = n = 0 C n , φ 2 φ x = n = 0 D n .
The values of A n , B n C n , and D n are known as the Adomian polynomials [53], obtained in the form
A n = 1 n ! d n d λ n k = 0 n λ k ϕ k 2 k = 0 n λ k ϕ k x λ = 0 , B n = 1 n ! d n d λ n k = 0 n λ k φ k 2 k = 0 n λ k ϕ k x λ = 0 , C n = 1 n ! d n d λ n k = 0 n λ k ϕ k 2 k = 0 n λ k φ k x λ = 0 , D n = 1 n ! d n d λ n k = 0 n λ k φ k 2 k = 0 n λ k φ k x λ = 0 .
which gives
A 0 = ϕ 0 2 ϕ 0 x , A 1 = ϕ 0 2 ϕ 1 x + 2 ϕ 0 ϕ 1 ϕ 0 x , A 2 = ϕ 1 2 ϕ 0 x + 2 ϕ 0 ϕ 2 ϕ 0 x + 2 ϕ 0 ϕ 1 ϕ 1 x + ϕ 0 2 ϕ 2 x , A n , B 0 = φ 0 2 ϕ 0 x , B 1 = φ 0 2 ϕ 1 x + 2 φ 0 φ 1 ϕ 0 x , B 2 = φ 1 2 ϕ 0 x + 2 φ 0 φ 2 ϕ 0 x + 2 φ 0 φ 1 ϕ 1 x + φ 0 2 ϕ 2 x , B n , C 0 = ϕ 0 2 φ 0 x , C 1 = ϕ 0 2 φ 1 x + 2 ϕ 0 ϕ 1 φ 0 x , C 2 = ϕ 1 2 φ 0 x + 2 ϕ 0 ϕ 2 φ 0 x + 2 ϕ 0 ϕ 1 φ 1 x + ϕ 0 2 φ 2 x , C n , D 0 = φ 0 2 φ 0 x , D 1 = φ 0 2 φ 1 x + 2 φ 0 φ 1 φ 0 x , D 2 = φ 1 2 φ 0 x + 2 φ 0 φ 2 φ 0 x + 2 φ 0 φ 1 φ 1 x + φ 0 2 φ 2 x , D n .
Applying L 1 to Equations (13) and (14), we obtain
n = 0 ϕ n = f ( x ) + L 1 1 s α L β t β 1 3 δ n = 0 A n 3 δ n = 0 B n γ n = 0 ϕ n x x x ,
n = 0 φ n = g ( x ) + L 1 1 s α L β t β 1 3 δ n = 0 C n 3 δ n = 0 D n γ n = 0 φ n x x x .
which give
ϕ 0 = f ( x ) , ϕ 1 = L 1 1 s α L β t β 1 3 δ A 0 3 δ B 0 γ ϕ 0 x x x , ϕ 2 = L 1 1 s α L β t β 1 3 δ A 1 3 δ B 1 γ ϕ 1 x x x , ϕ 3 = L 1 1 s α L β t β 1 3 δ A 2 3 δ B 2 γ ϕ 2 x x x , φ 0 = g ( x ) , φ 1 = L 1 1 s α L β t β 1 3 δ C 0 3 δ D 0 γ φ 0 x x x , φ 2 = L 1 1 s α L β t β 1 3 δ C 1 3 δ D 1 γ φ 1 x x x , φ 3 = L 1 1 s α L β t β 1 3 δ C 2 3 δ D 2 γ φ 2 x x x ,
The other terms can be calculated in the similar fashion.

4. Convergence Analysis

Here, the convergence of the approximate solution obtained through LADM is studied employing the following theorem.
Theorem 1. 
Let X and Y be Banach spaces, and T : X Y is a contraction. If c ( 0 , 1 ] such that x , y X ,
| | T ( x ) T ( y ) | | c | | x y | | ,
then by Banach contraction theorem, the approximate solutions (17) and (18) are convergent.
Proof. 
Let the series ϕ ( x , t ) = n = 0 ϕ n obtained by LADM in the form
ξ ( n ) = T ( ξ n 1 ) , ξ n 1 = i = 1 n 1 x i , n N ,
further, we consider that
T ( ϰ i ) = x i B r ( x ) where B r ( x ) = x 0 : | | x x 0 | | < r ,
then we have
  • ξ ( n ) B r ( x ) ;
  • lim n ξ ( n ) = x .
Proof. 
Mathematical induction can be applied to prove (1) and (2).
(1).
for n = 1 , we have
| | T ( ϰ 1 ) T ( x ) | | = | | ϰ 1 x | | c | | x 0 x | | .
Let the result be true for n 1 , then
| | ξ n 1 x | | c n 1 | | x 0 x | | , | | ξ ( n ) x | | = | | T ( ξ n 1 ) T ( x ) | | c | | ξ n 1 x | | c n | | x 0 x | | , | | ξ ( n ) x | | = | | T ( ξ n 1 ) T ( x ) | | c | | x 0 x | | c n | | x 0 x | | , | | ξ ( n ) x | | c n | | x 0 x | | c n r < r ,
which implies that ξ ( n ) B r ( x ) .
(2).
Since we have | | ξ ( n ) x | | c n | | x 0 x | | , and
lim n c n = 0 , therefore lim n | | ξ ( n ) x | | = 0 , lim n ξ ( n ) = x .
Similarly, the result can be proven for φ ( x , t ) = n = 0 φ n . □

5. Applications

In this section, the numerical examples on the coupled Hirota system is studied with detailed error analysis.
Example 1. 
Here, an example is provided that emphasizes the findings from the preceding section. Consider Equations (5) and (6) with initial conditions
ϕ ( x , 0 ) = A η sech ( η x ) cos ( ω x ) , φ ( x , 0 ) = A η sech ( η x ) sin ( ω x ) .
Plugging the estimation of the obtained Adomian polynomials and apply the proposed method gives
ϕ 0 = A η sech ( η x ) cos ( ω x ) , ϕ 1 = β ! t α + β 1 ( α + β 1 ) ! 3 δ ϕ 0 2 + φ 0 2 ϕ 0 x γ 3 ϕ 0 x 3 , ϕ 2 = β ! t α + β 1 ( α + β 1 ) ! 3 δ ϕ 0 2 + φ 0 2 ϕ 1 x 3 δ 2 ϕ 0 ϕ 1 + 2 φ 0 φ 1 ϕ 0 x γ 3 ϕ 1 x 3 , φ 0 = A η sech ( η x ) sin ( ω x ) , φ 1 = β ! t α + β 1 ( α + β 1 ) ! 3 δ ϕ 0 2 + φ 0 2 φ 0 x γ t 3 φ 0 x 3 , φ 2 = β ! t α + β 1 ( α + β 1 ) ! 3 δ ϕ 0 2 + φ 0 2 φ 1 x 3 δ 2 ϕ 0 ϕ 1 + 2 φ 0 φ 1 φ 0 x γ t 2 3 φ 1 x 3 , .
Further, inserting the values gives
ϕ 1 = β ! t α + β 1 A η ( α + β 1 ) ! [ γ ω 3 η 2 ω 2 cosh 3 ( η x ) sin ( ω x ) γ η 3 ω 2 η 2 sinh ( η x ) cos ( ω x ) cosh 2 ( η x ) + 3 ω η 2 A 2 δ 2 γ sin ( ω x ) cosh ( η x ) + 3 η 3 A 2 δ 2 γ sinh ( η x ) cos ( ω x ) ] sech 4 ( η x ) , φ 1 = β ! t α + β 1 A η ( α + β 1 ) ! [ γ ω w 2 3 η 2 cos ( ω x ) cosh 3 ( η x ) η γ 3 w 2 η 2 sin ( ω x ) sinh ( η x ) cosh 2 ( η x ) 3 ω η 2 A 2 δ cos ( ω x ) sin 2 ( ω x ) A 2 δ sin 3 ( ω x ) 2 γ cos ( ω x ) cosh ( η x ) + 3 η 3 sinh ( η x ) sin ( ω x ) A 2 δ cos ( ω x ) sin ( ω x ) + A 2 δ sin 2 ω x 2 γ ] sech 4 ( η x ) .
In a similar way, other values can be obtained. The complete result can be expressed as
ϕ ( x , t ) = ϕ 0 + ϕ 1 + ϕ 2 + ,
φ ( x , t ) = φ 0 + φ 1 + φ 2 + .
Equations (22) and (23) combined gives the final solution
ψ ( x , t ) = n = 0 ϕ n ( x , t ) + i φ n ( x , t ) = ϕ ( x , t ) + i φ ( x , t ) .
Here, the coupled Hirota system with order ( α , β ) in the Caputo’s sense with power law kernel is calculated via the combination of Adomian decomposition and Laplace transform. To demonstrate the solutions, a specific illustration is considered and obtained the solutions in the form of a series. The convergence of the suggested method also is examined.

Absolute Error Analysis

The error analysis among the exact solutions (2) against solutions (22) and (23) is provided in the table below.
From Table 1, it follows that the absolute error between accurate and approximately obtained solutions of the Hirota equation reduces when spatial variable increases for particular value of time (t). It is observed that combining in iterations reduces the absolute error.

6. Results and Discussion

For numerical demonstration of the consistent nonlinear configurations connected by Equations (2) and (24), the nondimensional variables are considered as: γ = 1 ,   δ = 1 ,   η = 0.2 and ω = 0.5 . It is worth mentioning that the dispersion (nonlinearity) effect γ ( δ ) originates from the oscillation steepening (spreading); η represents the wavenumber and the frequency; and ω illustrates the nonlinear variations (oscillations) of the wave packets, specified by Equations (5) and (6).
The numerical solution (24), together with the absolute value of the exact solution (2), are displayed versus the spatial variable x by the dashed and solid curves, as shown in Figure 1a. A pulse-shaped wave profile results from the wave packet’s internal oscillations dispersing. The comparison of both curves shows the confirmation of numerical solution. The three-dimensional description to numerical solution (24) is illustrated in Figure 1b against the temporal and spatial variables at δ = 0.25 , γ = 1 , and η = 1 . It reveals that solitary potentials are spatially localized, keeping a stable amplitude for a short time.
The solid curve in left panel of Figure 2 illustrates the real part of Equation (2) versus 10 x 10 at time (t) = 0.1 . The dashed curved line signifies the real function of the solution (22) for Equation (24) attained by the proposed method (LADM). One can see that the numerical solution accurately meets the obtained systematic approximation. In the right panel of Figure 2, the imaginary value of Equation (2) is compared with the numerical result obtained in Equation (23) versus 10 x 10 . It shows the excitations of perturbations, keeping the carrier wave, and as a result establishes a wave packet. The surface plots of analytical solutions (22) and (23) are depicted in Figure 3. One can see that the internal fluctuations of the waves become constant and sequentially localize the solitary wave. It is important to note that the soliton excitations effect caused by a stability of steepening and the spreading impacts of the nonlinear terms. The nonlinear wave steepening influences the wave breaking for relatively insignificant wave dispersion.
Figure 4 illustrates the influence of the fractal variable β by fixing fractional variable α for time (t = 0.1) for obtained solutions ϕ and φ . Similarly, the influence of the variable α with fixed value of β for obtained solutions ϕ and φ is shown in Figure 5. In conclusion, it is found that decreasing the fractal derivative β causes the amplitude to increase. Similarly, reducing the fractional order α modifies the shape of the solitonic solution and reduces its amplitude.
Figure 6 and Figure 7 show the behaviour of the solutions (22) and (23) with various values of α and β for particular value of the spatial variables ( x = 1 ) versus time ( t ) . It is noted that the solitary waves are in quite good agreement for comparatively particular rate of time ( t ) . It is also found that, for α = β = 1, increasing time ( t ) rapidly increases the wave propagation. Figure 8 constitutes the surface plots of the absolute error of the solutions obtained by LADM and accurate solutions for ϕ and φ for α = β = 1 , respectively.

7. Conclusions

The coupled nonlinear system was investigated analytically and numerically with fractal fractional derivative with power law kernel. The suggested method (LADM) is applied and calculated the general series solution. The convergence analysis is also examined. For validation, a particular example is considered, and the obtained results are competed with the numerical simulation with physical interpretations. The confinement of the carrier wave and the presence of confined wave packets are observed from the numerical analysis. It is observed that the dispersion nonlinear parameter reduction equivalently influenced the wave amplitude and spatial width. It is revealed that the amplitude of the solitary wave solution increases by decreasing the fractal dimension. It is also observed that decreasing the fractional order α reduces the nature of the solitary wave solutions and slightly decreases the amplitude. The numerical and approximation solutions resemble effectively for specific values of time (t). However, when the fractal or fractional derivative reaches its maximum value by increasing time, the wave amplitude increases. An error analysis is performed between the obtained series solutions versus the accurate solution. It is found that the error between the exact and approximation solutions is minimized for a sufficiently small value of time ( t ) .
As a future work, it will be interesting to study the fractal fractional Hirota equation and coupled Hirota equations with variable coefficients for high-order solitary wave solution, rogue wave solution and dark–bright solutions using the Mellin transform, Chebyshev collocation scheme and Finite difference method reported in [54,55].

Funding

Researchers Supporting Project number (RSP2024R447), King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Lu, D.; Hong, B.; Tian, L. New explicit exact solutions for the generalized coupled Hirota–Satsuma KdV system. Comput. Math. Appl. 2007, 53, 1181–1190. [Google Scholar] [CrossRef]
  2. Khan, K.; Ali, A.; Irfan, M.; Algahtani, O. Time-fractional electronacoustic shocks in magnetoplasma with superthermal electrons. Alex. Eng. J. 2022, 65, 531–542. [Google Scholar] [CrossRef]
  3. Irfan, M.; Alam, I.; Ali, A.; Shah, K.; Abdeljawad, T. Electron-acoustic solitons in dense electron-positron-ion plasma: Degenerate relativistic enthalpy function. Results Phys. 2022, 38, 105625. [Google Scholar] [CrossRef]
  4. Liu, G.-L. New research directions in singular perturbation theory: Artificial parameter approach and inverse-perturbation technique. In Proceedings of the 7th Conference of the Modern Mathematics and Mechanics, Shanghai, China, September 1997; Volume 61, pp. 47–53. [Google Scholar]
  5. He, J.H. A coupling method of homotopy technique and perturbation technique for nonlinear problems. Int. J. Non-Linear Mech. 2000, 35, 37–43. [Google Scholar] [CrossRef]
  6. Ali, A.; Gul, Z.; Khan, W.A.; Ahmad, S.; Zeb, S. Investigation of Fractional Order sine-Gordon Equation Using Laplace Adomian Decomposition Method. Fractals 2021, 29, 2150121. [Google Scholar] [CrossRef]
  7. Hirota, R. Exact envelope-soliton solutions of a nonlinear wave equation. J. Math. Phys. 1973, 14, 805–809. [Google Scholar] [CrossRef]
  8. Ankiewicz, A.; Soto-Crespo, J.M.; Akhmediev, N. Rogue waves and rational solutions of the Hirota equation. Phys. Rev. E 2010, 81, 046602. [Google Scholar] [CrossRef]
  9. Li, L.; Wu, Z.; Wang, L.; He, J. High-order rogue waves for the Hirota equation. Ann. Phys. 2013, 334, 198–211. [Google Scholar] [CrossRef]
  10. Anco, S.C.; Ngatat, N.T.; Willoughby, M. Interaction properties of complex modified Korteweg–de Vries (mKdV) solitons. Phys. D Nonlinear Phenom. 2011, 240, 1378–1394. [Google Scholar] [CrossRef]
  11. Guo, B.; Ling, L.; Liu, Q.P. Nonlinear schrödinger equation: Generalized darboux transformation and rogue wave solutions. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2012, 85, 026607. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Dong, K.; Jin, R. The Darboux transformation for the coupled Hirota equation. AIP Conf. Proc. 2013, 1562, 249–256. [Google Scholar]
  13. Huang, X. Rational solitary wave and rogue wave solutions in coupled defocusing Hirota equation. Phys. Lett. A 2016, 380, 2136–2141. [Google Scholar] [CrossRef]
  14. Jia, T.T.; Chai, Y.Z.; Hao, H.Q. Multi-soliton solutions and Breathers for the generalized coupled nonlinear Hirota equations via the Hirota method. Superlattices Microstruct. 2017, 105, 172–182. [Google Scholar] [CrossRef]
  15. Xin, W.; Yong, C. Rogue-wave pair and dark-bright-rogue wave solutions of the coupled Hirota equations. Chin. Phys. B 2014, 23, 070203. [Google Scholar]
  16. Chen, S. Dark and composite rogue waves in the coupled Hirota equations. Phys. Lett. A 2014, 378, 2851–2856. [Google Scholar] [CrossRef]
  17. Tao, Y.; He, J. Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2012, 85, 026601. [Google Scholar] [CrossRef]
  18. Dai, C.; Zhang, J. New solitons for the Hirota equation and generalized higher-order nonlinear Schrödinger equation with variable coefcients. J. Phys. Math. Gen. 2013, 39, 723–737. [Google Scholar] [CrossRef]
  19. Bhrawy, A.H.; Alshaery, A.A.; Hilal, E.M.; Manrakhan, W.N.; Savescu, M.; Biswas, A. Dispersive optical solitons with Schrödinger-Hirota equation. J. Nonlinear Opt. Phys. Mater. 2014, 23, 1450014. [Google Scholar] [CrossRef]
  20. Faddeev, L.; Volkov, A.Y. Hirota equation as an example of an integrable symplectic map. Lett. Math. Phys. 1994, 32, 125–135. [Google Scholar] [CrossRef]
  21. Eslami, M.; Mirzazadeh, M.A.; Neirameh, A. New exact wave solutions for hirota equation. Pramana J. Phys. 2015, 84, 3–8. [Google Scholar] [CrossRef]
  22. Saeed, R.K.; Muhammad, R.S. Solving Coupled Hirota System by Using Homotopy Perturbation and Homotopy Analysis Methods. J. Zankoi Sulaimani-Part Pure Appl. Sci. 2015, 201, 201–2017. [Google Scholar] [CrossRef]
  23. Saeed, R.K.; Muhammad, R.S. Solving Coupled Hirota System by Using Variational Iteration Method. Zanco J. Pure Appl. Sci. 2015, 201, 69–80. [Google Scholar]
  24. Tasgal, R.S.; Potasek, M.J. Soliton solutions to coupled higher-order nonlinear Schrödinger equations. J. Math. Phys. 1992, 33, 1208–1215. [Google Scholar] [CrossRef]
  25. Bindu, S.G.; Mahalingam, A.; Porsezian, K. Dark soliton solutions of the coupled Hirota equation in nonlinear fiber. Phys. Lett. A 2001, 286, 321–331. [Google Scholar] [CrossRef]
  26. Porsezian, K.; Nakkeeran, K. Optical solitons in birefringent fibre-Backlund transformation approach. Pure Appl. Opt. 1997, 6, L7–L11. [Google Scholar] [CrossRef]
  27. Chen, S.H.; Song, L.Y. Rogue waves in coupled Hirota systems. Phys. Rev. E 2013, 87, 032910. [Google Scholar] [CrossRef]
  28. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic: New York, NY, USA, 1974. [Google Scholar]
  29. Podlubny, I. Fractional Differential Equations; Academic: New York, NY, USA, 1999. [Google Scholar]
  30. Hilfer, R. Application of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  31. Zaslavsky, G.M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2005. [Google Scholar]
  32. Li, L. The cardiovascular promotion of college taekwondo based on fractional differential equation. Appl. Math. Nonlinear Sci. 2023, 8, 1577–1582. [Google Scholar] [CrossRef]
  33. Jiang, F. Regional Logistics Routing Optimization System Based on Fractional Differential Equation Modeling. Appl. Math. Nonlinear Sci. 2023, 8, 2865–2872. [Google Scholar] [CrossRef]
  34. Nishant; Bhatter, S.; Purohit, S.D.; Nisar, K.S.; Munjam, S.R. Some fractional calculus findings associated with the product of incomplete N-function and Srivastava polynomials. Int. J. Math. Comput. Eng. 2024, 2, 97–116. [Google Scholar] [CrossRef]
  35. Atangana, A.; Gómez-Aguilar, J.F. Numerical approximation of Riemann–Liouville definition of fractional derivative: From Riemann–Liouville to Atangana–Baleanu. Numer. Methods Partial. Differ. Equ. 2018, 34, 1502–1523. [Google Scholar] [CrossRef]
  36. Khader, M.M.; Saad, K.M.; Hammouch, Z.; Baleanu, D. A spectral collocation method for solving fractional KdV and KdV-Burgers equations with non-singular kernel derivatives. Appl. Numer. Math. 2021, 161, 137–146. [Google Scholar] [CrossRef]
  37. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel; Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef]
  38. Baleanu, D.; Fernandez, A.; Akgül, A. On a Fractional Operator Combining Proportional and Classical Differintegrals. Mathematics 2020, 9, 360. [Google Scholar] [CrossRef]
  39. Atangana, A.; Akgül, A.; Owolabi, K.M. Analysis of fractal fractional differential equations. Alex. Eng. J. 2020, 59, 1117–1134. [Google Scholar] [CrossRef]
  40. Atangana, A. Fractal-fractional differentiation and integration: Connecting fractal calculus and fractional calculus to predict complex system. Chaos Solit. Fractals. 2017, 102, 396–406. [Google Scholar] [CrossRef]
  41. Chen, W.; Sun, H.; Zhang, X.; Korošak, D. Anomalous diffusion modeling by fractal and fractional derivatives. Comput. Math. Appl. 2010, 59, 1754–1758. [Google Scholar] [CrossRef]
  42. Sun, H.G.; Meerschaert, M.M.; Zhang, Y.; Zhu, J.; Chen, W. A fractal Richards’ equation to capture the non-Boltzmann scaling of water transport in unsaturated media. Adv. Water Resour. 2013, 52, 292–295. [Google Scholar] [CrossRef]
  43. Saifullah, S.; Ali, A.; Goufo, E.F.D. Investigation of complex behaviour of fractal fractional chaotic attractor with mittag-leffler Kernel. Chaos Solitons Fract. 2021, 152, 111332. [Google Scholar] [CrossRef]
  44. Akgül, A.; Siddique, I. Analysis of MHD Couette flow by fractal-fractional differential operators. Chaos Solitons Fract. 2021, 146, 110893. [Google Scholar] [CrossRef]
  45. Atangana, A. Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos Solitons Fract. 2020, 136, 109860. [Google Scholar] [CrossRef]
  46. Algahtani, O.; Saifullah, S.; Ali, A. Semi-analytical and numerical study of fractal fractional nonlinear system under Caputo fractional derivative. AIMS Math. 2022, 7, 16760–16774. [Google Scholar] [CrossRef]
  47. Adomian, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers: Boston, MA, USA, 1995; p. 19. [Google Scholar]
  48. Tatari, M.; Dehghan, M. Numerical solution of Laplace equation in a disk using the Adomian decomposition method. Phys. Scr. 2005, 72, 345–348. [Google Scholar] [CrossRef]
  49. Andrianov, I.V.; Olevskii, V.I.; Tokarzewski, S. A modified Adomian’s decomposition method. Appl. Math. Mech. 1998, 62, 309–314. [Google Scholar] [CrossRef]
  50. Venkatarangan, S.N.; Rajalakshmi, K. A modification of Adomian’s solution for nonlinear oscillatory systems. Comput. Math. Appl. 1995, 29, 67–73. [Google Scholar] [CrossRef]
  51. Venkatarangan, S.N.; Rajalakshmi, K. Modification of Adomian’s decomposition method to solve equations containing radicals. Comput. Math. Appl. 1995, 29, 75–80. [Google Scholar] [CrossRef]
  52. Hussain, M.; Khan, M. Modified Laplace decomposition method. Appl. Math. Sci. 2010, 4, 1769–1783. [Google Scholar]
  53. Adomian, G. Modification of the decomposition approach to heat equation. J. Math. Anal. Appl. 1987, 124, 290–291. [Google Scholar] [CrossRef]
  54. Ata, E.; Onur Kıymaz, I. New generalized Mellin transform and applications to partial and fractional differential equations. Int. J. Math. Comput. Eng. 2023, 1, 45–66. [Google Scholar] [CrossRef]
  55. Demir, A.; Bayrak, M.A.; Bulut, A.; Ozbilge, E.; Çetinkaya, S. On new aspects of Chebyshev polynomials for space-time fractional diffusion process. Appl. Math. Nonlinear Sci. 2023, 8, 1051–1062. [Google Scholar] [CrossRef]
Figure 1. (a) Comparison between exact versus approximate solutions (2)/(24) for 6 x 6 . (b) The surface plot of approximate solution (24).
Figure 1. (a) Comparison between exact versus approximate solutions (2)/(24) for 6 x 6 . (b) The surface plot of approximate solution (24).
Fractalfract 08 00178 g001
Figure 2. (a) The real value representing accurate solution Re ( ψ ) versus numerical solution ϕ from Equations (2)/(22), are depicted versus x (b) The imaginary value representing accurate solution Im ( ψ ) /numerical solution φ from Equations (2)/(23) are depicted versus x.
Figure 2. (a) The real value representing accurate solution Re ( ψ ) versus numerical solution ϕ from Equations (2)/(22), are depicted versus x (b) The imaginary value representing accurate solution Im ( ψ ) /numerical solution φ from Equations (2)/(23) are depicted versus x.
Fractalfract 08 00178 g002
Figure 3. The wave profiles of analytical solutions (22) and (23) versus spatial and temporal variables.
Figure 3. The wave profiles of analytical solutions (22) and (23) versus spatial and temporal variables.
Fractalfract 08 00178 g003
Figure 4. The influence of stable fractional order α with variation in fractal variable β with time ( t = 0.1 ) for obtained solutions ϕ and φ .
Figure 4. The influence of stable fractional order α with variation in fractal variable β with time ( t = 0.1 ) for obtained solutions ϕ and φ .
Fractalfract 08 00178 g004
Figure 5. The influence of stable fractal order β with variation in fractional order α for the obtained solutions ϕ and φ at time (t).
Figure 5. The influence of stable fractal order β with variation in fractional order α for the obtained solutions ϕ and φ at time (t).
Fractalfract 08 00178 g005
Figure 6. The conduct of the functions ϕ ( x , t ) and φ ( x , t ) for distinct estimates of β keeping α stable versus time ( t ) .
Figure 6. The conduct of the functions ϕ ( x , t ) and φ ( x , t ) for distinct estimates of β keeping α stable versus time ( t ) .
Fractalfract 08 00178 g006
Figure 7. The behaviour of ϕ ( x , t ) (Equation (22)) and φ ( x , t ) (Equation (23)) for uneven values of α by keeping β stable versus time ( t ) .
Figure 7. The behaviour of ϕ ( x , t ) (Equation (22)) and φ ( x , t ) (Equation (23)) for uneven values of α by keeping β stable versus time ( t ) .
Fractalfract 08 00178 g007
Figure 8. Plot for error analysis reported in Table 1.
Figure 8. Plot for error analysis reported in Table 1.
Fractalfract 08 00178 g008
Table 1. The analyses are performed for α = β = η = γ = 1 , δ = 0.25 , a = 1.3 for Example 1.
Table 1. The analyses are performed for α = β = η = γ = 1 , δ = 0.25 , a = 1.3 for Example 1.
xTime (t)Exact ( Re ( ψ ) ) ϕ (Equation (22))∣ Exact− ϕ Exact ( Im ( ψ ) ) φ (Equation (23))∣ Exact− φ
−40.010.05220.05580.00370.09440.09230.0021
−3 −0.2089−0.20290.00600.20480.21080.0060
−2 −0.6767−0.68420.0075−0.3916−0.37860.0130
−1 −0.2089−0.20290.00600.20480.21080.0060
0 −0.6767−0.68420.0075−0.3916−0.37860.0130
1 0.47450.45780.0167−1.8295−1.83490.0054
2′ 2.82572.82840.0028−0.0481−0.04760.0006
3 0.50430.48750.01681.70331.70710.0037
4 −0.6130−0.61950.00650.38310.37060.0125
−40.050.06880.07050.00170.10670.10370.0030
−3 −0.2286−0.22220.00650.25690.25700.0001
−2 −0.8197−0.81370.0061−0.4025−0.39470.0078
−1 0.39130.39840.0071−2.0833−2.09050.0072
0 2.76082.82840.0676−0.2357−0.23790.0022
1 0.54280.54690.00411.45991.45150.0084
2 −0.4998−0.49000.00980.36190.35450.0074
3 −0.1793−0.17390.0054−0.1430−0.14150.0016
4 0.03310.03380.0007−0.0778−0.07530.0024
−40.10.09510.08880.00630.12310.11780.0053
−3 −0.2520−0.24630.00570.33700.31480.0221
−2 −1.0292−0.97550.0537−0.4007−0.41480.0140
−1 0.24030.32410.0838−2.3828−2.41000.0272
0 2.57162.82840.2569−0.4422−0.47570.0335
1 0.55790.62120.06321.18181.13200.0498
2 −0.3827−0.32820.05460.33020.33440.0042
3 −0.1558−0.14980.0060−0.1039−0.08370.0202
4 0.02150.01550.0060−0.0655−0.06120.0044
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Algahtani, O.J. Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System. Fractal Fract. 2024, 8, 178. https://doi.org/10.3390/fractalfract8030178

AMA Style

Algahtani OJ. Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System. Fractal and Fractional. 2024; 8(3):178. https://doi.org/10.3390/fractalfract8030178

Chicago/Turabian Style

Algahtani, Obaid J. 2024. "Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System" Fractal and Fractional 8, no. 3: 178. https://doi.org/10.3390/fractalfract8030178

APA Style

Algahtani, O. J. (2024). Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System. Fractal and Fractional, 8(3), 178. https://doi.org/10.3390/fractalfract8030178

Article Metrics

Back to TopTop