Next Article in Journal
Extinction Dynamics and Equilibrium Patterns in Stochastic Epidemic Model for Norovirus: Role of Temporal Immunity and Generalized Incidence Rates
Next Article in Special Issue
The Investigation of Nonlinear Time-Fractional Models in Optical Fibers and the Impact Analysis of Fractional-Order Derivatives on Solitary Waves
Previous Article in Journal
Properties and Applications of Complex Fractal–Fractional Operators in the Open Unit Disk
Previous Article in Special Issue
Analytical Scheme for Time Fractional Kawahara and Modified Kawahara Problems in Shallow Water Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solutions of the Fractional Hirota–Satsuma Coupled KdV Equation along with Analysis of Bifurcation, Sensitivity and Chaotic Behaviors

by
Yongyi Gu
1,
Chunling Jiang
2,* and
Yongkang Lai
1
1
Department of Statistics and Mathematics, Guangdong University of Finance and Economics, Guangzhou 510320, China
2
Department of Basic Courses Teaching, Software Engineering Institute of Guangzhou, Guangzhou 510990, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(10), 585; https://doi.org/10.3390/fractalfract8100585
Submission received: 1 August 2024 / Revised: 12 September 2024 / Accepted: 1 October 2024 / Published: 3 October 2024

Abstract

:
This paper explores the exact solutions of the fractional Hirota–Satsuma coupled KdV (fHScKdV) equation in the Beta fractional derivative. The logistic method is first proposed to construct analytical solutions for the fHScKdV equation. In order to better comprehend the physical structure of the solutions, three-dimensional visualizations and line graphs of the exponent function solutions are depicted with the aid of Matlab. Furthermore, the phase portraits and bifurcation behaviors of the fHScKdV model under transformation are studied. Sensitivity and chaotic behaviors are analyzed in specific conditions. The phase plots and time series map are exhibited through sensitivity analysis and perturbation factors. These investigations enhance our understanding of practical phenomena governed by the fHScKdV model, and are crucial for examining the dynamic behaviors and phase portraits of the fHScKdV system. The strategies utilized here are more direct and effective, and can be applied effortlessly to other fractional order differential equations.

1. Introduction

Due to the inherent complexity of various phenomena, linear systems are often theoretical approximating of simpler nonlinear systems. However, nonlinear systems more accurately capture the essence of the objective world. Therefore, comprehending and studying nonlinear phenomena is crucial for modern science and technology. From a mathematical physics standpoint, numerous nonlinear phenomena can be reduced to solve nonlinear differential equations [1,2,3]. Among these, nonlinear fractional differential equations (NLFDEs) have garnered considerable attention. NLFDEs represent an extension of partial differential equations in which the order of derivatives can be fractional rather than integer-valued. These equations have significant applications in physics, engineering, biology, finance, and numerous other fields, particularly in describing systems with memory effects, nonlocal phenomena, or complex dynamics. NLFDEs typically involve fractional derivatives with respect to spatial or temporal variables. Fractional derivatives can be defined in several ways, the most common of which are the Caputo [4], Grünwald–Letnikov [5], Riemann–Liouville [6] and Beta [7] definitions. NLFDEs symbolise a vibrant research field that provides powerful tools in describing the dynamics of complex systems. As computational methods and theories continue to evolve, we can see that even more novel applications and solution techniques are emerging for NLFDEs.
In the literature, many effective approaches have been constructed to search analytical solutions of NLFDEs, such as the F-expansion method [8], the fractional sub-equation method [9,10,11], the generalized Riccati equation mapping method [12], the first integral method [13,14], the Kudryashov methods [15,16], the exp ( φ ( z ) ) -expansion method [17,18], the ( G / G ) -expansion method [19,20,21,22], the tanh-function method [23,24], the modified simple equation method [25,26,27], the truncated Painlevé expansion method [23], the Sine–Gordon expansion method [28], the complex method [29,30,31] and so on.
The HScKdV equation [32] typically describes the interaction between two long waves and different dispersion relations. By extending it to the fractional form, the fHScKdV equation captures additional layers of complexity that cannot be fully captured by integer-order derivatives. The fHScKdV equation [33] is given as follows:
D t β u = 1 4 u x x x + 3 u u x + 3 ( v 2 + w ) x , D t β v = 1 2 v x x x 3 u v x , D t β w = 1 2 w x x x 3 u w x , 0 < β 1 ,
where D β denotes the Beta differential operator.
The fHScKdV equation has potential applications in various fields, including physics, engineering, and finance. For example, in fluid mechanics, it can be used to model wave and turbulence phenomena with fractional properties. In materials science, it can describe the behavior of materials with memory effects. In finance, it can analyze market dynamics, exhibiting long-memory effects. The fHScKdV equation is an interesting system that attracts researchers to analyze it in depth. Alam et al. [33] acquired the closed-form solutions using the ( G / G ) -expansion method as well as the Riemann–Liouville derivative. Yin et al. [34] utilized the homogeneous balance method, improved F-expansion method, and unified method, and extracted various analytical solutions. Ganji et al. [35] exploited the homotopy perturbation method to obtain exact solutions. Yan [36] constructed exact doubly periodic solutions by applying the extended Jacobian elliptic function expansion technique. Kaplan and Bekir [37] applied the exp ( φ ( z ) ) -expansion method and local fractional derivatives to derived exact solutions for the fHScKdV equation. They were all committed to using effective methods to quickly obtain solutions to nonlinear differential equations. Inspired by these studies, we put forward the logistic approach, which is founded on the notion that the precise solutions of a nonlinear differential equation can be expounded by a polynomial in N ( η ) , and N ( η ) complies with the logistic ordinary differential equation. Our results indicate that the logistic method is direct and effective. This enriches the study of the fHScKdV equation and may provide reference for researchers in related fields. In this paper, we propose a new approach, named the logistic method, to study the fHScKdV equation, and investigate its bifurcation, sensitivity and chaotic behaviors.
The paper is organized as follows. The logistic method is introduced in Section 2 to execute analytical solutions of the fHScKdV equation. Section 3 extracts solutions to the fHScKdV equation, and shows the dynamic behaviors of these solutions in 3D, contour plot and line plots. Phase portraits and bifurcation behavior analysis of the fHScKdV system are presented in Section 4. In Section 5, sensitivity analysis to the initial value is explained. The chaotic behavior regarding the fHScKdV equation is interpreted in Section 6. The conclusions are given in Section 7.

2. Proposal of the Logistic Method

Recently, the Beta fractional derivative (FD) concept was put forward by some researchers [38], signifying significant progress in the study of mathematical derivatives. In particular, Beta FD has a stronger ability to precisely predict real-time phenomena compared to the standard derivative. The key advantage of FD is its non-locality, which reveals the effect of distant elements on the behavior of a system and brings fundamental value. A wide range of fields, like dielectric polarization, viscoelasticity, electrochemistry, and image processing, extensively employ Beta FD. Most of these applications are in engineering and physics, for instance wave theory, biotechnology, mass transfer, heat, robotics, and so on.
First, some definitions of the Beta fractional derivative are introduced and a detailed explanation of the logistic method is given to better elucidate the obtained results.
Let f ( z ) be specified as a function of all non-negative z with the β derivative [39,40], hence
D z β ( f ( z ) ) = d β f ( z ) d z β = lim δ 0 f z + δ z + 1 Γ ( β ) 1 β f ( z ) δ , 0 < β 1 ,
in which Γ ( β ) is the Gamma function. Some effective properties [39,40,41,42] of the above definition are included:
D z β ( f ( z ) ) = z + 1 Γ ( β ) 1 β d f ( z ) d z ,
D z β ( f g ( z ) ) = z + 1 Γ ( β ) 1 β g ( z ) f ( g ( z ) ) .
Using this knowledge about Beta derivatives, analytical solutions of the fHScKdV equation are sought by exerting the logistic method, which is presented according to the following steps.
Consider a space-time NLFDE as follows:
F ( u , D x γ u , D t β u , D x γ D t β u , D x 2 γ u , ) = 0 , 0 < γ 1 , 0 < β 1 ,
where F is the polynomial of unknown function u with its fractional derivatives and D x 2 γ u means D x γ ( D x γ u ) .
Mechanism I. Inserting the transformation
η = μ γ x + 1 Γ γ γ λ β t + 1 Γ β β .
into space-time NLFDE (5), we yield
P ( u , u , u , ) = 0 ,
where u = d u d η , and P is the polynomial of u along with its derivatives.
Mechanism II. Assume that Equation (7) has the following solution:
U ( η ) = c 0 + i = 1 n c i N ( η ) i , c n 0 ,
where c i for i = 0 , 1 , . . . , n are constants, and N ( η ) satisfies the following logistic equation:
N ( η ) = s 0 N ( η ) 1 N ( η ) s 1 ,
Mechanism III. At this step, impose the homogeneous balancing method on the highest order derivative of U and highest order nonlinear term in Equation (7), and we can gain the integer n of (8). Apply (8) and Equation (9) to Equation (7), and then we are able to acquire an algebra system through the gathering N with the same order. Utilizing the results of the algebra system, we can obtain the values of c i and related restriction conditions.
Mechanism IV. By substituting the results of previous steps into Equation (8) and applying the general solutions of Equation (9) as follows,
N ( η ) = s 1 1 + a s 1 exp ( s 0 η ) ,
provided a , s 0 , s 1 are arbitrary constants, we can achieve closed-form solutions of Equation (5).

3. Applications to the fHScKdV Equation

We start with the application of the following transform to (1):
u ( x , t ) = 1 λ U ( η ) 2 , v ( x , t ) = λ + U ( η ) , w ( x , t ) = 2 λ 2 2 λ U ( η ) , η = x λ β t + 1 Γ β β .
By substituting the transform (11) into (1), we can obtain the following system:
4 U U = 1 2 λ ( 3 U U + U U ) + 6 λ 2 U 3 U , λ U = 1 2 U 3 λ U 2 U , 2 λ 2 U = λ U + 6 U 2 U .
By combining the first and second equations of (12), and eliminating U , we arrive at
U ( λ U 2 λ 2 U + 2 U 3 ) = 0
To seek nontrivial solutions, we proceed by reducing the fHScKdV equation to the following ordinary differential equation:
λ d 2 U d η 2 2 λ 2 U + 2 U 3 = 0 .
By balancing between U and U 3 in Equation (14), then n = 1 can be determined for (8), which yields the solution of U, as follows:
U = c 0 + c 1 N ,
where N = N ( η ) is satisfied (9). To gain meaningful solutions, displacing (15) along with Equation (9) into Equation (14), we can collect coefficients of the same order of N j , ( j = 0 , . . . , 3 ) to zero. The algebraic system is attained as below:
N 3 : 2 c 1 c 1 2 s 1 2 + λ s 0 2 = 0 , N 2 : 3 c 1 λ s 0 2 2 c 0 c 1 s 1 = 0 , N 1 : λ c 1 s 0 2 2 λ 2 c 1 + 6 c 0 2 c 1 = 0 , N 0 : 2 λ 2 c 0 + 2 c 0 3 = 0 ,
and the results produced for coefficients are as below.
Set 1:
c 0 = λ , c 1 = 2 λ s 1 , s 0 = 2 λ .
Set 2:
c 0 = λ , c 1 = 2 λ s 1 , s 0 = 2 λ .
By inserting these results into (15) and using the definition of N, the following solutions are extracted for Equation (1).
Case 1:
U 1 ( η ) = λ a s 1 exp ( 2 λ η ) 1 a s 1 exp ( 2 λ η ) + 1 , u 1 ( η ) = λ a s 1 exp ( 2 λ η ) 1 a s 1 exp ( 2 λ η ) + 1 2 , v 1 ( η ) = λ + λ a s 1 exp ( 2 λ η ) 1 a s 1 exp ( 2 λ η ) + 1 , w 1 ( η ) = 2 λ 2 2 λ 2 a s 1 exp ( 2 λ η ) 1 a s 1 exp ( 2 λ η ) + 1 ,
with η = x λ β t + 1 Γ β β .
Case 2:
U 2 ( η ) = λ 1 2 1 + a s 1 exp η , u 2 ( η ) = λ 1 2 1 + a s 1 exp η 2 , v 2 ( η ) = λ 1 + 2 1 + a s 1 exp η , w 2 ( η ) = 4 λ 2 2 λ 2 2 1 + a s 1 exp η ,
where η = x λ β t + 1 Γ β β .
Matlab is a powerful software developed by MathWorks, designed for numerical computing, data analysis, algorithm development, and visualization. With the help of Matlab (Version 9.14.0.2306882 (R2023a)), we can reveal the physics performance by presenting the simulations. Since u i , v i , and w i are based on the function U i and are similar to each other, we have only shown the simulation of v i . We are particularly interested in the impact of the fractional derivative β on the function.
To describe the effect of β , we set the parameters λ = 2 , a = 2 , and s 1 = 2 , while varying the value of β from 1 to 0.1 . The graphs of v 1 , v 2 are depicted in Figure 1 and Figure 2, respectively. In Figure 1a–e, the behaviors of v 1 are illustrated for β values of 1 , 0.99 , 0.9 , 0.5 , and 0.1 , respectively. The characteristics at x = 0 for various values of β are shown in Figure 1f. Similarly, v 2 uses the same parameters as v 1 , and the corresponding plots are shown in Figure 2a–f.
Figure 1a–e and Figure 2a–e clearly demonstrate that as β decreases, the region for t < 0 is compressed into a constant value, and a jump region gradually forms near t = 0 . Additionally, based on Figure 1a,b and Figure 2a,b, even a slight deviation of β from 1 disrupts the original periodicity. These findings play a significant role in the study and analysis of the dynamic behavior of the fHScKdV system. The choice of β affects the nonlinear properties of the function, altering its smoothness near t = 0 . Therefore, it is crucial to carefully consider the impact of the fractional order β when studying related systems. Furthermore, these properties may potentially characterize tsunami; however, further research is still needed.

4. Phase Portraits and Bifurcation Behaviors

If
d U d η = V ,
we can rewrite (14) with Hamiltonian dynamics:
d U d η = V , d V d η = 2 λ U 2 λ U 3 .
The Hamiltonian of system (21) is
H = 1 2 V 2 λ U 2 + 1 2 λ U 4 .
The system (21) has three equilibrium points:
E 1 ( 0 , 0 ) , E 2 ( λ , 0 ) , E 3 ( λ , 0 ) .
The Jacobian of (21) is
J ( U , V ) = 0 1 2 λ 6 λ U 2 0 = 6 λ U 2 2 λ .
Then, we have
J | E 1 = 2 λ , J | E 2 , 3 = 4 λ .
For the equilibrium point E of a Hamiltonian system, E can be analyzed in three situations: when J | E > 0 , E is a center; when J | E = 0 , E is a cuspid; and when J | E < 0 , E is a saddle. Hence, we conclude that the bifurcation behavior of system (21) is as follows:
Case I: If λ > 0 , E 1 is a saddle and E 2 , 3 are center.
Case II: If λ < 0 , E 1 is a center and E 2 , 3 are saddle.
To analyze system (21), we list the effect of bifurcation by plotting the phase portraits of the system in Figure 3. Case I is shown in Figure 3a by letting the value of the system parameter λ = 1 , while Figure 3b depicts case II by setting the value of system parameter λ = 1 . From Figure 3, we can study the bifurcation behaviors when the values of parameters alter. Therefore, the dynamics behaviors are remarkably affected by the bifurcations.

5. Sensitive Analysis

Sensitivity analysis pertains to the examination of how a system’s response is sensitive to alterations in parameters or initial conditions. It represents a mathematical method for evaluating the influence of variations within a variable framework on its output. Comprehending the capability and dependability of dynamic structures is of utmost significance. This type of analysis is frequently utilized to explore the changes in configurations or variables that impact the performance of systems in various disciplines like dynamical frameworks and energy.
We examined the sensitivity behaviors of Equation (14) concerning the initial values of cases I and II, as previously discussed. The results for the parameter λ according to case I are illustrated in Figure 4, while Figure 5 shows the outcomes for case II. In Figure 4 and Figure 5, the red curve represents the initiation at ( U , V ) = ( 0.15 , 0.1 ) , and the blue curve represents the initiation at ( U , V ) = ( 0.75 , 0.1 ) . Figure 4a and Figure 5a exhibit the time series for the model (21). In addition, Figure 4b and Figure 5b display phase portraits of system (21) initiated from two distinct initial conditions, highlighting sensitivity impacts. To evaluate the impact of sensitivity variations, we undertook a comparative analysis of the pathways emanating from diverse initial conditions, specifically focusing on the distinguishing patterns between the blue and red curves presented in Figure 4a,b and Figure 5a,b. The results indicate that significant changes in initial conditions have a considerable impact on the dynamic behavior of the system when the remaining parameters remain constant. Moreover, we investigated the system’s sensitivity to parameter variations by comparing Figure 4a,b with Figure 5a,b, demonstrating a significant sensitivity of the system to alterations in its parameters.

6. Chaotic Behavior

Chaos, in the realm of science and mathematics, refers specifically to the highly complex and completely unpredictable behavior that can be observed within nonlinear systems. This concept implies a state in which the patterns and outcomes are not straightforward or easily foreseen, presenting a level of intricacy that makes it extremely challenging to accurately predict or understand the subsequent developments in such systems. Through the study of chaotic phenomena, we can acquire comprehension of the nonlinear characteristics of the fHScKdV equation and disclose their intricate temporal evolution courses. It is of great significance to consider the dynamic behaviors of system (14) in this paper.
In this segment, our primary concern is the chaotic dynamics exhibited by Equation (14) when subjected to noise-induced perturbations. We undertook an analysis of the chaotic nature of Equation (14) by contemplating a periodic perturbation model as a framework for our discussion. The model is read as:
d U d η = V , d V d η = 2 λ U 2 λ U 3 + f 0 cos ( ω η ) ,
where ω and f 0 denote the frequency and amplitude of noise, respectively. Equation (26) is an exact forced non-linear pendulum. There have been many similar studies, such as Neishtadt et al.’s [43] study on the dynamics of nonlinear pendulums under periodic forces with small amplitudes and slowly decreasing frequencies. Bazzani et al. [44] applied the perturbative techniques to Hamiltonians. We analyzed the sensitivity of the noise as follows. When either the amplitude or frequency was fixed and the other one varied, we studied the time series as well as phase portraits of the system (21). To study the effect of amplitude f 0 , frequency ω was fixed; the time series and phase portraits are shown in Figure 6 and Figure 7. Similarly, we discuss the impacts of frequency ω when amplitude f 0 is fixed. The results are shown in Figure 8 and Figure 9. In addition, parameter λ is fixed as λ = 1 in Figure 6 and Figure 8, while λ = 1 is set in Figure 7 and Figure 9.
We valued the parameters f 0 and ω as follows. In Figure 6 and Figure 7, we kept ω = 0.1 and let f 0 = 0 . The red curve represents the simulation results for the time series and phase orbits. Then, we adjusted f 0 to 0.1 and 0.5 ; the graphical depiction of the simulation outcomes is distinctively presented via the blue and green lines, respectively. The simulation highlights a pivotal observation: the divergence between the red line and the aforementioned curves underscores the emergence of chaotic system dynamics contingent upon the presence of noise. Notably, in both scenarios under consideration, the system manifests a pronounced sensitivity to variations in the amplitude parameter f 0 .
In Figure 8 and Figure 9, we set f 0 to 0. The red curve illustrates the simulation results for both the time series and phase orbits. We then adjusted f 0 to 0.5 and ω to 0.1 and 0.5 , with the results represented by the blue and green curves, respectively. The distinction between the red curve and the others indicates that the system exhibits chaotic behavior while the noise term is present. In both cases, the system’s sensitivity to the frequency ω is significant.

7. Conclusions

In this paper, a transformation using the Γ ( β ) function is applied to simplify fractional differential equations into ordinary differential equations. The exact solutions of the fHScKdV equation are constructed by proposing and utilizing the logistic method. These derived solutions enable valuable understanding of the inherent physical phenomena in such systems. Furthermore, three-dimensional visualizations and line graphs are presented to depict the influence of different β -order derivatives on solutions. By analyzing the plots, we conclude that the choice of β has a significant impact on the function. Therefore, the careful selection of β is essential in related research. Additionally, we studied the phase portraits and bifurcation behaviors of the fHScKdV model under wave transformation, examining the related sensitivity as well as chaotic behaviors. These investigations not only augment existing knowledge but also shed invaluable light on the intricate dynamics that govern the system, offering fresh insights for further exploration and understanding.

Author Contributions

Conceptualization, Y.G.; Methodology, C.J.; Software, Y.L.; Validation, Y.L.; Writing—original draft, Y.G. and Y.L.; Writing—review and editing, C.J.; Supervision, Y.G. and C.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the NSFC (11901111), Science Research Group Project of SEIG (ST202101) and Science Research Project of SEIG (ky202211).

Data Availability Statement

All data were contained in the main text.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ullah, M.S.; Ali, M.Z.; Roshid, H.-O.; Seadawy, A.; Baleanu, D. Collision phenomena among lump, periodic and soliton solutions to a (2+1)-dimensional bogoyavlenskii’s breaking soliton model. Phys. Lett. A 2021, 397, 127263. [Google Scholar] [CrossRef]
  2. Wu, X.-S.; Liu, J.-G. Solving the variable coefficient nonlinear partial differential equations based on the bilinear residual network method. Nonlinear Dyn. 2024, 112, 8329–8340. [Google Scholar] [CrossRef]
  3. Kumar, D.; Seadawy, A.R.; Haque, M.R. Multiple soliton solutions of the nonlinear partial differential equations describing the wave propagation in nonlinear low-pass electrical transmission lines. Chaos Solitons Fractals 2018, 115, 62–76. [Google Scholar] [CrossRef]
  4. Tuan, N.H.; Mohammadi, H.; Rezapour, S. A mathematical model for COVID-19 transmission by using the Caputo fractional derivative. Chaos Solitons Fractals 2020, 140, 110107. [Google Scholar] [CrossRef]
  5. Machado, J.A.T. The bouncing ball and the Grünwald-Letnikov definition of fractional derivative. Fract. Calc. Appl. Anal. 2021, 24, 1003–1014. [Google Scholar] [CrossRef]
  6. Hilfer, R. Fractional diffusion based on Riemann-Liouville fractional derivatives. J. Phys. Chem. B 2000, 104, 3914–3917. [Google Scholar] [CrossRef]
  7. Bagheri, M.; Khani, A. Analytical method for solving the fractional order generalized KdV equation by a Beta-fractional derivative. Adv. Math. Phys. 2020, 2020, 8819183. [Google Scholar] [CrossRef]
  8. Ozkan, E.M. New exact solutions of some important nonlinear fractional partial differential equations with beta derivative. Fractal Fract. 2022, 6, 173. [Google Scholar] [CrossRef]
  9. Tang, B.; He, Y.; Wei, L.; Zhang, X. A generalized fractional sub-equation method for fractional differential equations with variable coefficients. Phys. Lett. A 2012, 376, 2588–2590. [Google Scholar] [CrossRef]
  10. Bekir, A.; Aksoy, E.; Cevikel, A. Exact solutions of nonlinear time fractional partial differential equations by sub-equation method. Math. Methods Appl. Sci. 2014, 38, 2779–2784. [Google Scholar] [CrossRef]
  11. Nisar, K.S.; Inan, I.E.; Inc, M.; Rezazadeh, H. Properties of some higher-dimensional nonlinear schrödinger equations. Results Phys. 2021, 31, 105073. [Google Scholar] [CrossRef]
  12. Khater, M.M.A.; Jhangeer, A.; Rezazadeh, H.; Akinyemi, L.; Akbar, M.A.; Inc, M. Propagation of new dynamics of longitudinal bud equation among a magneto-electro-elastic round rod. Mod. Phys. Lett. B 2021, 35, 2150381. [Google Scholar] [CrossRef]
  13. Mirzazadeh, M.; Eslami, M.; Biswas, A. Solitons and periodic solutions to a couple of fractional nonlinear evolution equations. Pramana 2014, 82, 465–476. [Google Scholar] [CrossRef]
  14. Eslami, M.; Fathi-Vajargah, B.; Mirzazadeh, M.; Biswas, A. Application of first integral method to fractional partial differential equations. Indian J. Phys. 2014, 88, 177–184. [Google Scholar] [CrossRef]
  15. Hosseini, K.; Mirzazadeh, M.; Baleanu, D.; Raza, N.; Park, C.; Ahmadian, A.; Salahshour, S. The generalized complex Ginzburg-Landau model and its dark and bright soliton solutions. Eur. Phys. J. Plus 2021, 136, 709. [Google Scholar] [CrossRef]
  16. Hosseini, K.; Mirzazadeh, M.; Salahshour, S.; Baleanu, D.; Zafar, A. Specific wave structures of a fifth-order nonlinear water wave equation. J. Ocean Eng. Sci. 2022, 7, 462–466. [Google Scholar] [CrossRef]
  17. Kaplan, M.; Bekir, A. A novel analytical method for time-fractional differential equations. Optik 2016, 127, 8209–8214. [Google Scholar] [CrossRef]
  18. Hosseini, K.; Bekir, A.; Ansari, R. Exact solutions of nonlinear conformable time-fractional boussinesq equations using the exp(-φ(ϵ))-expansion method. Opt. Quantum Electron. 2017, 49, 131. [Google Scholar] [CrossRef]
  19. Khan, K.; Akbar, M.A.; Koppelaar, H. Study of coupled nonlinear partial differential equations for finding exact analytical solutions. R. Soc. Open Sci. 2015, 2, 140406. [Google Scholar] [CrossRef]
  20. Zafar, A.; Raheel, M.; Asif, M.; Hosseini, K.; Mirzazadeh, M.; Akinyemi, L. Some novel integration techniques to explore the conformable m-fractional schrödinger-hirota equation. J. Ocean Eng. Sci. 2022, 7, 337–344. [Google Scholar] [CrossRef]
  21. Nisar, K.S.; Ali, K.K.; Inc, M.; Mehanna, M.S.; Rezazadeh, H.; Akinyemi, L. New solutions for the generalized resonant nonlinear schrödinger equation. Results Phys. 2022, 33, 105153. [Google Scholar] [CrossRef]
  22. Nisar, K.S.; Inan, I.E.; Yepez-Martinez, H.; Inc, M. Some new type optical and the other soliton solutions of coupled nonlinear hirota equation. Results Phys. 2022, 35, 105388. [Google Scholar] [CrossRef]
  23. Hu, H.; Li, X. New interaction solutions of the similarity reduction for the integrable (2+1)-dimensional boussinesq equation. Int. J. Mod. Phys. B 2022, 36, 2250001. [Google Scholar] [CrossRef]
  24. Raslan, K.R.; Ali, K.K.; Shallal, M.A. The modified extended tanh method with the riccati equation for solving the space-time fractional EW and MEW equations. Chaos Solitons Fractals 2017, 103, 404–409. [Google Scholar] [CrossRef]
  25. Kaplan, M.; Bekir, A.; Akbulut, A.; Aksoy, E. The modified simple equation method for nonlinear fractional differential equations. Rom. J. Phys. 2015, 60, 1374. [Google Scholar]
  26. Kaplan, M.; Koparan, M.; Bekir, A. Regarding on the exact solutions for the nonlinear fractional differential equations. Open Phys. 2016, 14, 478–482. [Google Scholar] [CrossRef]
  27. Rahman, Z.; Abdeljabbar, A.; Harun-Or-Roshid; Ali, M.Z. Novel precise solitary wave solutions of two time fractional nonlinear evolution models via the MSE scheme. Fractal Fract. 2022, 6, 444. [Google Scholar] [CrossRef]
  28. Kumar, D.; Hosseini, K.; Kaabar, M.K.A.; Kaplan, M.; Salahshour, S. On some novel solution solutions to the generalized Schrödinger-Boussinesq equations for the interaction between complex short wave and real long wave envelope. J. Ocean Eng. Sci. 2022, 7, 353–362. [Google Scholar] [CrossRef]
  29. Gu, Y.; Yuan, W.; Aminakbari, N.; Lin, J. Meromorphic solutions of some algebraic differential equations related Painlevé equation IV and its applications. Math. Methods Appl. Sci. 2018, 41, 3832–3840. [Google Scholar] [CrossRef]
  30. Gu, Y.; Wu, C.; Yao, X.; Yuan, W. Characterizations of all real solutions for the KdV equation and WR. Appl. Math. Lett. 2020, 107, 106446. [Google Scholar] [CrossRef]
  31. Gu, Y.; Liao, L. Closed form solutions of Gerdjikov-Ivanov equation in nonlinear fiber optics involving the beta derivatives. Int. J. Mod. Phys. B 2022, 36, 2250116. [Google Scholar] [CrossRef]
  32. Hirota, R.; Satsuma, J. Soliton solutions of a coupled Korteweg-de Vries equation. Phys. Lett. A 1981, 85, 407–408. [Google Scholar] [CrossRef]
  33. Alam, M.N.; Seadawy, A.R.; Baleanu, D. Closed-form wave structures of the space-time fractional Hirota-Satsuma coupled KdV equation with nonlinear physical phenomena. Open Phys. 2020, 18, 555–565. [Google Scholar] [CrossRef]
  34. Yin, Q.; Gao, B. New solutions of the time-fractional Hirota-Satsuma coupled KdV equation by three distinct methods. Int. J. Geom. Methods Mod. Phys. 2023, 20, 2350170. [Google Scholar] [CrossRef]
  35. Kurt, A.; Rezazadeh, H.; Senol, M.; Neirameh, A.; Tasbozan, O.; Eslami, M.; Mirzazadeh, M. Two effective approaches for solving fractional generalized Hirota-Satsuma coupled KdV system arising in interaction of long waves. J. Ocean Eng. Sci. 2019, 4, 24–32. [Google Scholar] [CrossRef]
  36. Yan, Z. The extended jacobian elliptic function expansion method and its application in the generalized Hirota-Satsuma coupled KdV system. Chaos Solitons Fractals 2003, 15, 575–583. [Google Scholar] [CrossRef]
  37. Kaplan, M.; Bekir, A. Construction of exact solutions to the space-time fractional differential equations via new approach. Optik 2017, 132, 1–8. [Google Scholar] [CrossRef]
  38. Atangana, A.; Baleanu, D.; Alsaedi, A. Analysis of time-fractional hunter-saxton equation: A model of neumatic liquid crystal. Open Phys. 2016, 14, 145–149. [Google Scholar] [CrossRef]
  39. Atangana, A.; Goufo, E.F.D. Extension of matched asymptotic method to fractional boundary layers problems. Math. Probl. Eng. 2014, 2014, e107535. [Google Scholar] [CrossRef]
  40. Atangana, A.; Alqahtani, R.T. Modelling the spread of river blindness disease via the Caputo fractional derivative and the beta-derivative. Entropy 2016, 18, 40. [Google Scholar] [CrossRef]
  41. Hosseini, K.; Mirzazadeh, M.; Gómez-Aguilar, J.F. Soliton solutions of the Sasa-Satsuma equation in the monomode optical fibers including the beta-derivatives. Optik 2020, 224, 165425. [Google Scholar] [CrossRef]
  42. Hosseini, K.; Kaur, L.; Mirzazadeh, M.; Baskonus, H.M. 1-Soliton solutions of the (2+1)-dimensional Heisenberg ferromagnetic spin chain model with the beta time derivative. Opt. Quantum Electron. 2021, 53, 125. [Google Scholar] [CrossRef]
  43. Neishtadt, A.I.; Vasiliev, A.A.; Artemyev, A.V. Capture into resonance and escape from it in a forced nonlinear pendulum. Regul. Chaotic Dyn. 2013, 18, 686–696. [Google Scholar] [CrossRef]
  44. Bazzani, A.; Capoani, F.; Giovannozzi, M.; Tomás, R. Nonlinear cooling of an annular beam distribution. Phys. Rev. Accel. Beams 2024, 26, 024001. [Google Scholar] [CrossRef]
Figure 1. (af) | v 1 ( x , t ) | at λ = 2 , a = 2 , s 1 = 2 ; (a) β = 1 ; (b) β = 0.99 ; (c) β = 0.9 ; (d) β = 0.5 ; (e) β = 0.1 ; (f) x = 0 .
Figure 1. (af) | v 1 ( x , t ) | at λ = 2 , a = 2 , s 1 = 2 ; (a) β = 1 ; (b) β = 0.99 ; (c) β = 0.9 ; (d) β = 0.5 ; (e) β = 0.1 ; (f) x = 0 .
Fractalfract 08 00585 g001
Figure 2. (af) | v 2 ( x , t ) | at λ = 2 , a = 2 , s 1 = 2 ; (a) β = 1 ; (b) β = 0.99 ; (c) β = 0.9 ; (d) β = 0.5 ; (e) β = 0.1 ; (f) x = 0 .
Figure 2. (af) | v 2 ( x , t ) | at λ = 2 , a = 2 , s 1 = 2 ; (a) β = 1 ; (b) β = 0.99 ; (c) β = 0.9 ; (d) β = 0.5 ; (e) β = 0.1 ; (f) x = 0 .
Fractalfract 08 00585 g002
Figure 3. Phase portraits of case I (a) and case II (b). The stars represent the equilibriums.
Figure 3. Phase portraits of case I (a) and case II (b). The stars represent the equilibriums.
Fractalfract 08 00585 g003
Figure 4. Sensitivity analysis of case I.
Figure 4. Sensitivity analysis of case I.
Fractalfract 08 00585 g004
Figure 5. Sensitivity analysis of case II.
Figure 5. Sensitivity analysis of case II.
Fractalfract 08 00585 g005
Figure 6. Chaotic behavior analysis and the amplitude for case I.
Figure 6. Chaotic behavior analysis and the amplitude for case I.
Fractalfract 08 00585 g006
Figure 7. Chaotic behavior analysis and amplitude for case II.
Figure 7. Chaotic behavior analysis and amplitude for case II.
Fractalfract 08 00585 g007
Figure 8. Chaotic behavior analysis and frequency for case I.
Figure 8. Chaotic behavior analysis and frequency for case I.
Fractalfract 08 00585 g008
Figure 9. Chaotic behavior analysis and frequency for case II.
Figure 9. Chaotic behavior analysis and frequency for case II.
Fractalfract 08 00585 g009
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

Gu, Y.; Jiang, C.; Lai, Y. Analytical Solutions of the Fractional Hirota–Satsuma Coupled KdV Equation along with Analysis of Bifurcation, Sensitivity and Chaotic Behaviors. Fractal Fract. 2024, 8, 585. https://doi.org/10.3390/fractalfract8100585

AMA Style

Gu Y, Jiang C, Lai Y. Analytical Solutions of the Fractional Hirota–Satsuma Coupled KdV Equation along with Analysis of Bifurcation, Sensitivity and Chaotic Behaviors. Fractal and Fractional. 2024; 8(10):585. https://doi.org/10.3390/fractalfract8100585

Chicago/Turabian Style

Gu, Yongyi, Chunling Jiang, and Yongkang Lai. 2024. "Analytical Solutions of the Fractional Hirota–Satsuma Coupled KdV Equation along with Analysis of Bifurcation, Sensitivity and Chaotic Behaviors" Fractal and Fractional 8, no. 10: 585. https://doi.org/10.3390/fractalfract8100585

APA Style

Gu, Y., Jiang, C., & Lai, Y. (2024). Analytical Solutions of the Fractional Hirota–Satsuma Coupled KdV Equation along with Analysis of Bifurcation, Sensitivity and Chaotic Behaviors. Fractal and Fractional, 8(10), 585. https://doi.org/10.3390/fractalfract8100585

Article Metrics

Back to TopTop