Next Article in Journal
Delay Effects on Plant Stability and Symmetry-Breaking Pattern Formation in a Klausmeier-Gray-Scott Model of Semiarid Vegetation
Previous Article in Journal
Analogies between Lattice QCD and the Truncated Nambu–Jona-Lasinio Model
Previous Article in Special Issue
Non-Contact Impact Source Localization in Composite Symmetry Panels Based on A0 Mode of Lamb Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing Dynamics: Lie Symmetry Approach to Bifurcation, Chaos, Multistability, and Solitons in Extended (3 + 1)-Dimensional Wave Equation

1
IT4Innovations, VSB—Technical University of Ostrava, Ostrava, Czech Republic
2
Department of Computer Science and Mathematics, Lebanese American University, Byblos, Lebanon
3
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(5), 608; https://doi.org/10.3390/sym16050608
Submission received: 15 April 2024 / Revised: 2 May 2024 / Accepted: 6 May 2024 / Published: 14 May 2024
(This article belongs to the Special Issue Symmetry in the Soliton Theory)

Abstract

:
The examination of new (3 + 1)-dimensional wave equations is undertaken in this study. Initially, the identification of the Lie symmetries of the model is carried out through the utilization of the Lie symmetry approach. The commutator and adjoint table of the symmetries are presented. Subsequently, the model under discussion is transformed into an ordinary differential equation using these symmetries. The construction of several bright, kink, and dark solitons for the suggested equation is then achieved through the utilization of the new auxiliary equation method. Subsequently, an analysis of the dynamical nature of the model is conducted, encompassing various angles such as bifurcation, chaos, and sensitivity. Bifurcation occurs at critical points within a dynamical system, accompanied by the application of an outward force, which unveils the emergence of chaotic phenomena. Two-dimensional plots, time plots, multistability, and Lyapunov exponents are presented to illustrate these chaotic behaviors. Furthermore, the sensitivity of the investigated model is executed utilizing the Runge–Kutta method. This analysis confirms that the stability of the solution is minimally affected by small changes in initial conditions. The attained outcomes show the effectiveness of the presented methods in evaluating solitons of multiple nonlinear models.

1. Introduction

In the present age of technological advancements, significant progress has been analyzed on the subject of solitons. They play a key role in transferring data over extended distances through fiber optics. Solitons also matter significantly as they are pivotal for the advancement of computer systems’ computing power, while they present a wide variety of applications. Such applications cover image processing, data analysis, neurology, and fluids, among other fields [1,2,3]. As a result, experts have actively explored the evolution of soliton solutions for nonlinear partial differential equations (NLPDEs). Various effective methods have been employed to obtain soliton solutions for NLPDEs, including the Darboux transform [4], Hirota’s approach [5], the Painlevé test [6], the sub-equation method [7], the unified approach [8], the bilinear approach [9], and multiple techniques discussed in [10].
Among these approaches, Lie analysis stands out as one of the systematic approaches for deriving closed-form solutions to NLPDEs. In recent decades, the application of Lie symmetries has been widespread, extending to numerous instances in both mathematics and engineering [11]. Upon conducting Lie analysis on nonlinear models, symmetry reductions and group-invariant solutions become accessible. These resultant symmetries frequently prove instrumental in diminishing both the number of independent variables in PDEs and the order of an ordinary differential equation (ODE). The academic community has demonstrated significant interest in the exploration of symmetries, leading to a substantial body of noteworthy research in this field. In their work, Kumar et al. [11] applied the Lie approach to examine the generalized Kadomtsev Petviashvili (KP) equation, while the authors in reference [12] scrutinized the Oskolkov model using Lie analysis and extracted several soliton structures.
Substantial investigations have previously explored the behavior of solitons, utilizing various nonlinear frameworks like the Sakovich equation [13], Schrödinger equation [14], and Zakharov–Kuznetsov modified equal-width equation [15]. The integrable equations in (3 + 1) dimensions have been the subject of increased research attention lately. These equations play a key role in elucidating the physics behind numerous scientific and engineering phenomena. Various nonlinear extended equations in (3 + 1) dimensions, including the Korteweg–de Vries (KdV) equation, modified KP equation, and others, have been formulated in response to this growing interest.
Akinyemi [16] examined the equation in (2 + 1) dimensions:
A Υ x t + a Υ x x + b ( Υ 2 ) x x + c Υ x x x x + d Υ y y = 0 ,
and its extended version, comprising two additional linear terms:
A Υ x t + a Υ x x + b ( Υ 2 ) x x + c Υ x x x x + d Υ y y + γ 1 Υ t y + γ 2 Υ t t = 0 ,
Here, the constants A , a , b , c , d , γ 1 , and γ 2 in the given context are unrestricted real values. It is worth highlighting that when A = γ 1 = d = 0 , Equation (2) simplifies to the Boussinesq equation:
Υ t t + a Υ x x + b ( Υ 2 ) x x + c Υ x x x x = 0 .
Equations (1) and (2) have been demonstrated to exhibit Painlevé integrability, and their multiple solitons have also been obtained [16]. In the ongoing study, our intention is to address an extended version of Equation (1) [17]:
Υ x t + a Υ x x + b ( Υ 2 ) x x + c Υ x x x x + d Υ y y + e Υ x y + h Υ x z + k Υ y z + m Υ z z = 0 ,
where Υ = Υ ( x , y , z , t ) . It is evident that Equation (4) is constructed by introducing four additional linear terms, specifically, e Υ x y , h Υ x z , k Υ y z , and m Υ z z , to Equation (1). Furthermore, the coefficients a , b , c , d , e , h , k , and m are random real values. In their study [17], the authors investigated the integrability criteria for the discussed models (4) through the use of the Painlevé test. The study concluded by analyzing a set of lump solutions for the suggested equation.
This study deals with the use of the new auxiliary equation (NAE) method [12,18] in order to analyze as well as to generate soliton structures for the considered model. It is worth mentioning here that the methodology discussed in this work has not been applied in prior research to address Equation (4). Through the adoption of this novel approach, the obtained results manifest into different forms, such as exponential, trigonometric, hyperbolic, and rational functions. Additionally, the findings are depicted graphically using three-dimensional, two-dimensional, and density plots, thereby improving the clarity of the results. The applications of dynamical systems are diverse, covering fields such as biology, engineering, and economics [19].
In general, scholars have demonstrated a growing interest in investigating many aspects of dynamical systems, including bifurcation analysis, chaos, and sensitivity. Bifurcation studies how an orbit varies with respect to different parameters. The majority of continuous dynamical systems using differential equations (DEs) have parameters. It is possible for a little change in a parameter to have a big effect on the response. It is a qualitative examination of the planar dynamical system that takes place for distinct values of the physical parameters that are present in the system. Equilibria, periodic orbits, and other invariant sets typically experience changes in their local stability characteristics during a bifurcation.(Please revise this part to reduce duplication.) It has various applications in different fields such as climate science, ecology, economics, and finance. In deterministic nonlinear systems that behave irregularly and erratically, chaos can also arise. Numerous researchers have studied the chaos and bifurcation phenomena within planar systems. For instance, Kazmi et al. [20] reported the soliton solutions and examined the dynamics of bifurcation and chaos in the q-deformed Sinh–Gordon model. Similarly, Rafiq et al. [21] investigated the qualitative nature of shallow water equations, uncovering multiwave solitons. Additionally, Hosseini et al. [22] analyzed the dynamic characteristics of the Schrödinger equation. The chaotic phenomena of the suggested equation is analyzed through the presentation of two-dimensional plots, time plots, and Lyapunov exponents analysis. Furthermore, a sensitivity [23] of the system has been conducted employing the Runge–Kutta (RK) method.
In this work, we an explored extended integrable wave equation (EIWE) (4) from three distinct viewpoints:
  • Initially, the Lie approach is utilized to identify translational symmetries. Subsequently, leveraging these symmetries, the model undergoes a transformation into an ODE. The solution to the resulting ODE is obtained through the application of the NAE approach, facilitating the examination of soliton structures.
  • Secondly, it is essential to explore the parameter-dependent analysis of the studied system through the lenses of bifurcation and chaos. A detailed explanation of bifurcation for the undisturbed dynamic system is provided, illustrated via phase plots. Additionally, it explores chaotic tendencies in the perturbed dynamical system employing the phase plot, time plot, and Lyapunov exponents tool.
  • Lastly, we examine the sensitivity of the considered model under various initial conditions. Slight modifications to the initial values can lead to a subtle change in the outcome of the system. Therefore, our observation suggests that the suggested system exhibits sensitivity, though not to an extreme extent.
The structure of the manuscript is as follows: In Section 2, we identify the Lie symmetries governing the equation and transform the model into an ODE using the extracted symmetries. The analytical solution of the model and results obtained from the analysis are discussed in Section 3. A comprehensive investigation of the dynamic characteristics of the proposed equation, utilizing phase portraits of bifurcation, is conducted in Section 4. Chaotic phenomena are explored in Section 5. Section 6 is dedicated to a discussion on the multistability analysis of the system. Sensitivity analysis is covered in Section 7. Finally, the conclusion is provided in the last section.

2. Lie Symmetries

Lie analysis focuses on discerning transformations that maintain the integrity of a given DE. These symmetries are often derived from Lie operators, which collectively form a Lie algebra [24,25]. The application of these operators to the DE allows for the identification of symmetries that uphold its inherent structure. Extracting symmetries and reducing the dimensions of the PDE involves the following key steps:
Step 1: Identify classical Lie point symmetries for the suggested model.
Step 2: Construct an algebra (in our case, an abelian algebra) based on the identified symmetries.
Step 3: Determine similarity variables corresponding to each symmetry.
Step 4: Utilize the obtained symmetries to reduce the PDE to an ODE.
Step 5: Derive solutions in the form of traveling waves from the resulting ODE.
Let us examine a group of infinitesimal transformations characterized by a single parameter in Lie theory.
x * = x + ϵ Ψ 1 ( x , y , z , t , Υ ) + O ( ϵ 2 ) , y * = y + ϵ Ψ 2 ( x , y , z , t , Υ ) + O ( ϵ 2 ) , z * = z + ϵ Ψ 3 ( x , y , z , t , Υ ) + O ( ϵ 2 ) , t * = t + ϵ Ψ 4 ( x , y , z , t , Υ ) + O ( ϵ 2 ) , Υ * = Υ + ϵ Θ ( x , y , z , t , Υ ) + O ( ϵ 2 ) ,
where group parameter ϵ 1 and Ψ 1 , Ψ 2 , Ψ 3 , Ψ 4 , and Θ are the coefficient functions. An infinitesimal generator associated with the aforementioned Lie group can be expressed as follows:
F = Ψ 1 ( x , y , z , t , Υ ) x + Ψ 2 ( x , y , z , t , Υ ) y + Ψ 3 ( x , y , z , t , Υ ) z + Ψ 4 ( x , y , z , t , Υ ) t + Θ ( x , y , z , t , Υ ) Υ .
The invariance condition for Equation (4) with F becomes
P r ( 4 ) F Υ x t + a Υ x x + b ( Υ 2 ) x x + c Υ x x x x + d Υ y y + e Υ x y + h Υ x z + k Υ y z + m Υ z z Eq . ( 4 ) = 0 = 0 ,
where P r ( 4 ) F is the fourth prolongation of F , defined as
P r ( 4 ) F = F + Θ x Υ x + Θ x t Υ x t + Θ x x Υ x x + Θ y y Υ y y + Θ x y Υ x y + Θ x z Υ x z + Θ y z Υ y z + Θ z z Υ z z + Θ x x x x Υ x x x x ,
where
Θ x = D x ( Θ ) Υ x D x ( Ψ 1 ) Υ y D x ( Ψ 2 ) Υ z D x ( Ψ 3 ) Υ t D x ( Ψ 4 ) , Θ x t = D t ( Θ x ) Υ t x D t ( Ψ 1 ) Υ t y D t ( Ψ 2 ) Υ t z D t ( Ψ 3 ) Υ t t D t ( Ψ 4 ) , Θ x x = D x ( Θ x ) Υ x x D x ( Ψ 1 ) Υ x y D x ( Ψ 2 ) Υ x z D x ( Ψ 3 ) Υ x t D x ( Ψ 4 ) , Θ y y = D y ( Θ y ) Υ x y D y ( Ψ 1 ) Υ y y D y ( Ψ 2 ) Υ y z D y ( Ψ 3 ) Υ y t D y ( Ψ 4 ) , Θ x y = D y ( Θ x ) Υ x y D y ( Ψ 1 ) Υ y y D y ( Ψ 2 ) Υ y z D y ( Ψ 3 ) Υ y t D y ( Ψ 4 ) , Θ x z = D z ( Θ x ) Υ x z D z ( Ψ 1 ) Υ y z D z ( Ψ 2 ) Υ z z D z ( Ψ 3 ) Υ t z D z ( Ψ 4 ) , Θ y z = D z ( Θ y ) Υ x z D z ( Ψ 1 ) Υ y z D z ( Ψ 2 ) Υ z z D z ( Ψ 3 ) Υ t z D z ( Ψ 4 ) , Θ z z = D z ( Θ z ) Υ x z D z ( Ψ 1 ) Υ y z D z ( Ψ 2 ) Υ z z D z ( Ψ 3 ) Υ z t D z ( Ψ 4 ) , Θ x x x x = D x ( Θ x x x ) Υ x x x x D x ( Ψ 1 ) Υ x x x y D x ( Ψ 2 ) Υ x x x z D x ( Ψ 3 ) Υ x x x t D x ( Ψ 4 ) ,
Here D x , D y , D z , and D t are total derivatives with respect to x, y, z, and t; it can written as
D x = x + Υ x Υ + Υ x x Υ x + Υ x y Υ y + Υ x z Υ z + Υ x t Υ t , D y = y + Υ y Υ + Υ x y Υ x + Υ y y Υ y + Υ y z Υ z + Υ y t Υ t , D z = z + Υ z Υ + Υ x z Υ x + Υ y z Υ y + Υ z z Υ z + Υ z t Υ t , D t = t + Υ t Υ + Υ x t Υ x + Υ y t Υ y + Υ z t Υ z + Υ t t Υ t .
Combining system (7) and Equation (6) into Equation (5) produces a set of equations for determination. Solving this system results in the derivation of the following symmetries for Equation (4):
F 1 = D x , F 2 = D y , F 3 = D z , F 4 = D t , F 5 = G 1 ( t ) D Υ + 2 b G 1 ( t ) D x , F 6 = 2 d h t e k t 2 d z + k y 4 d m k 2 D x + t D z , F 7 = 2 e m t h k t + k z 2 m y 4 d m k 2 D x + t D y , F 8 = e z h y k D x + 2 d z k y k D y + k z 2 m y k D z , F 9 = G 2 ( t ) z D Υ 2 ( 4 m G 2 ( t ) d h 4 m G 2 ( t ) z d G 2 ( t ) h k 2 + 2 m G 2 ( t ) y k + G 2 ( t ) k 2 z 2 m G 2 ( t ) y k ) b 4 d m k 2 D x + ( 2 G 2 ( t ) b k ) d t D y 4 G 2 ( t ) b m D z , F 10 = Υ D Υ 3 2 t D t ( 8 a d k m t 2 a k 3 t 2 d h 2 k t 2 e 2 k m t + 2 e h k 2 t 8 d e m z + 2 d h k z + 8 d h m y + 4 d k m x + e k 2 z + 2 e k m y 3 h k 2 y k 3 x ) k ( 4 d m k 2 ) D x + 2 ( d z k y ) k D y 2 m y k D z , F 11 = G 3 ( t ) y D Υ + 1 4 d m k 2 ( 4 G 3 ( t ) b d k z 2 G 3 ( t ) b k 2 y 4 G 3 ( t ) b k d z + 8 G 3 ( t ) b d m y + 4 2 b ( 2 k G 3 ( t ) d h 2 k G 3 ( t ) d z k 2 G 3 ( t ) e + k 2 G 3 ( t ) y + 2 G 3 ( t ) d k z G 3 ( t ) k 2 y + 4 G 3 ( t ) d e m 2 G 3 d h k ) 4 d m k 2 d t ) m d ( 2 2 b ( 2 k G 3 ( t ) d h 2 k G 3 ( t ) d z k 2 G 3 ( t ) e + k 2 G 3 ( t ) y + 2 G 3 ( t ) d k z G 3 ( t ) k 2 y + 4 G 3 ( t ) d e m 2 G 3 d h k ) 4 d m k 2 ) d t ) k 2 ) D x + ( ( 4 G 3 ( t ) b d ) d t ) D y 2 G 3 ( t ) b k D z .
  • Commutator table:
    [ F i , F j ] F 1 F 2 F 3 F 4
    F 1 0000
    F 2 0000
    F 3 0000
    F 4 0000
    where [ F i , F j ] = F i ( F j ) F j ( F i ) , i = 1 , 2 , 3 , 4 .
  • Adjoint table:
    F 1 F 2 F 3 F 4
    F 1 F 1 F 2 F 3 F 4
    F 2 F 1 F 2 F 3 F 4
    F 3 F 1 F 2 F 3 F 4
    F 4 F 1 F 2 F 3 F 4
    where = A d ( ε F i , F j ) = F j ε [ F i , F j ] + ε 2 2 [ F i , [ F i , F j ] ] , .

Traveling Wave Structures

One can notice that the set of F 1 , F 2 , F 3 , F 4 forms an abelian subalgebra. Our objective is to identify the traveling wave structures in the given model, which are associated with translational symmetries of the form:
Υ = β 1 x + β 2 y + β 3 z + μ t .
Through the application of the above linear combination denoted as the Υ of symmetries, we establish the transformation in the following manner:
Υ ( x , y , z , t ) = Γ ( ζ ) , ζ = ( β 1 x + β 2 y + β 3 z ) + μ t .
Here, Γ ( ζ ) and μ represent the characteristics of the traveling wave, specifically referring to its shape and velocity. Additionally, β 1 , β 2 , β 3 serve as random parameters. Upon substituting the expression from Equation (8) into Equation (4), we obtain the resulting equation:
( c β 1 4 ) Γ ( 4 ) + ( β 1 μ + a β 1 2 + d β 2 2 + e β 1 β 2 + h β 1 β 3 + k β 2 β 3 + m β 3 2 ) Γ + 2 b β 1 2 Γ Γ + 2 b β 1 2 ( Γ ) 2 = 0 .
By integrating Equation (9) twice with respect to ζ , and setting the integration constants to zero, we obtain
( c β 1 4 ) Γ + ( β 1 μ + a β 1 2 + d β 2 2 + e β 1 β 2 + h β 1 β 3 + k β 2 β 3 + m β 3 2 ) Γ + b β 1 2 Γ 2 = 0 .

3. Analytical Solutions

In the following section, we examine the soliton solutions of Equation (4) by employing the NAE approach [12,18]. This approach encompasses an initial solution that can be expressed as follows:
Γ ( ζ ) = i = 0 j κ i Λ i Ω ( ζ ) ,
where κ i is a random real parameter such that κ j 0 . The positive integer j serves as the balancing number, determined through the principle of homogeneous balance. This principle entails balancing the highest order derivative with the nonlinear term in the ODEs. By setting the terms Γ and Γ 2 equal in Equation (10) as 2 j = j + 2 , it yields j = 2 . Utilizing this approach for j = 2 , the initial solution becomes
Γ ( ζ ) = κ 0 + κ 1 Λ Ω ( ζ ) + κ 2 Λ 2 Ω ( ζ ) .
Here, κ 0 , κ 1 , and κ 2 are arbitrary parameters to be determined. Additionally, an auxiliary function Ω = Ω ( ζ ) is introduced, which satisfies the following ODE:
Ω ( ζ ) = 1 ln ( Λ ) α Λ Ω ( ζ ) + ω + ρ Λ Ω ( ζ ) .
In this context, α , ρ , and ω represent arbitrary constants. Upon substituting the values from Equations (12) and (13) into Equation (10) and equating the coefficients of Λ Ω ( ξ ) to zero, a set of equations is obtained. Solving this system, with the assistance of software such as Mathematica, yields solutions for κ 0 , κ 1 , κ 2 , and a. The obtained solutions are as follows:
κ 0 = 6 c β 1 2 ρ α b , κ 1 = 6 c β 1 2 ρ ω b , κ 2 = 6 c β 1 2 ρ 2 b , a = d β 2 2 μ β 1 e β 1 β 2 h β 1 β 3 k β 2 β 3 c ω 2 β 1 4 + 4 c ρ α β 1 4 m β 3 2 β 1 2 .
By substituting the aforementioned values into Equation (12) and applying the transformation outlined in Equation (8), we derive the solutions for Equation (4).
  • Type 1: For ρ 0 and ω 2 4 ρ α < 0 , we have
    Υ 1 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 ρ + 4 α ρ ω 2 2 ρ tan 4 α ρ ω 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 ρ + 4 α ρ ω 2 2 ρ tan 4 α ρ ω 2 2 ζ 2 . Υ 1 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 ρ 4 α ρ ω 2 2 ρ cot 4 α ρ ω 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 ρ 4 α ρ ω 2 2 ρ cot 4 α ρ ω 2 2 ζ 2 .
  • Type 2: For ρ 0 and ω 2 4 ρ α > 0 , we have
    Υ 2 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 ρ ω 2 4 α ρ 2 ρ tanh ω 2 4 α ρ 2 ζ 6 c β 1 2 ρ 2 b ω 2 ρ ω 2 4 α ρ 2 ρ tanh ω 2 4 α ρ 2 ζ 2 . Υ 2 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 ρ ω 2 4 α ρ 2 ρ coth ω 2 4 α ρ 2 ζ 6 c β 1 2 ρ 2 b ω 2 ρ ω 2 4 α ρ 2 ρ coth ω 2 4 α ρ 2 ζ 2 .
  • Type 3: For ρ 0 , ρ = α , and ω 2 α 2 < 0 , we have
    Υ 3 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω ρ + ( ω 2 α 2 ) ρ tan ( ω 2 α 2 ) 2 ζ 6 c β 1 2 ρ 2 b ω ρ + ( ω 2 α 2 ) ρ tan ( ω 2 α 2 ) 2 ζ 2 . Υ 3 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω ρ + ( ω 2 α 2 ) ρ cot ( ω 2 α 2 ) 2 ζ 6 c β 1 2 ρ 2 b ω ρ + ( ω 2 α 2 ) ρ cot ( ω 2 α 2 ) 2 ζ 2 .
  • Type 4: For ρ 0 , ρ = α , and ω 2 α 2 > 0 , we have
    Υ 4 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω ρ + ( ω 2 α 2 ) ρ tanh ( ω 2 α 2 ) 2 ζ 6 c β 1 2 ρ 2 b ω ρ + ( ω 2 α 2 ) ρ tanh ( ω 2 α 2 ) 2 ζ 2 . Υ 4 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω ρ + ( ω 2 α 2 ) ρ coth ( ω 2 α 2 ) 2 ζ 6 c β 1 2 ρ 2 b ω ρ + ( ω 2 α 2 ) ρ coth ( ω 2 α 2 ) 2 ζ 2 .
  • Type 5: For ρ = α and ω 2 4 α 2 < 0 , we have
    Υ 5 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 α + 4 α 2 ω 2 2 α tan 4 α 2 ω 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 α + 4 α 2 ω 2 2 α tan 4 α 2 ω 2 2 ζ 2 . Υ 5 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 α 4 α 2 ω 2 2 α cot 4 α 2 ω 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 α 4 α 2 ω 2 2 α cot 4 α 2 ω 2 2 ζ 2 .
  • Type 6: For ρ = α and ω 2 4 α 2 > 0 , we have
    Υ 6 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 α ω 2 4 α 2 2 α tanh ω 2 4 α 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 α ω 2 4 α 2 2 α tanh ω 2 4 α 2 2 ζ 2 . Υ 6 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω 2 α ω 2 4 α 2 2 α coth ω 2 4 α 2 2 ζ 6 c β 1 2 ρ 2 b ω 2 α ω 2 4 α 2 2 α coth ω 2 4 α 2 2 ζ 2 .
  • Type 7: For ω 2 = 4 ρ α , we have
    Υ 7 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b 2 + ω ζ 2 ρ ζ 6 c β 1 2 ρ 2 b 2 + ω ζ 2 ρ ζ 2 .
  • Type 8: For ρ α < 0 , ω = 0 , and ρ 0 , we have
    Υ 8 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b α ρ tanh α ρ ζ 6 c β 1 2 ρ 2 b α ρ tanh α ρ ζ 2 . Υ 8 , 2 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b α ρ coth α ρ ζ 6 c β 1 2 ρ 2 b α ρ coth α ρ ζ 2 .
  • Type 9: For α = ρ , ω = 0 , we have
    Υ 9 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b 1 + e 2 ρ 1 e 2 ρ ζ 6 c β 1 2 ρ 2 b 1 + e 2 ρ 1 e 2 ρ ζ 2 .
  • Type 10: For ρ = ω = L , α = 0 , we have
    Υ 10 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b e L ζ 1 e L ζ 6 c β 1 2 ρ 2 b e L ζ 1 e L ζ 2 .
  • Type 11: For ω = α + ρ , we have
    Υ 11 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b 1 α e ( α ρ ) ζ 1 ρ e ( α ρ ) ζ 6 c β 1 2 ρ 2 b 1 α e ( α ρ ) ζ 1 ρ e ( α ρ ) ζ 2 .
  • Type 12: For ω = α ρ , we have
    Υ 12 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b e ( ρ + α ) ζ α e ( ρ + α ) ζ ρ 6 c β 1 2 ρ 2 b e ( ρ + α ) ζ α e ( ρ + α ) ζ ρ 2 .
  • Type 13: For α = ω = ρ 0 , we have
    Υ 13 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b 1 2 3 tan 3 2 α ζ 1 6 c β 1 2 ρ 2 b 1 2 3 tan 3 2 α ζ 1 2 .
  • Type 14: For α = 0 , we have
    Υ 14 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ω e ω ζ 1 ρ e ω ζ 6 c β 1 2 ρ 2 b ω e ω ζ 1 ρ e ω ζ 2 .
  • Type 15: For ω = α = 0 , we have
    Υ 15 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b 1 ρ ζ 6 c β 1 2 ρ 2 b 1 ρ ζ 2 .
  • Type 16: For ω = ρ = 0 , we have
    Υ 16 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b ( α ζ ) 6 c β 1 2 ρ 2 b ( α ζ ) 2 .
  • Type 17: For α = ρ , ω = 0 , we have
    Υ 17 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b tan ( α ζ ) 6 c β 1 2 ρ 2 b tan ( α ζ ) 2 .
  • Type 18: For α = ρ = 0 , we have
    Υ 18 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b sinh ( ω ζ ) + cosh ( ω ζ ) 6 c β 1 2 ρ 2 b sinh ( ω ζ ) + cosh ( ω ζ ) 2 .
  • Type 19: For α = ω = L , ρ = 0 , we have
    Υ 19 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b e L ζ 1 6 c β 1 2 ρ 2 b e L ζ 1 2 .
  • Type 20: For ρ = 0 , we have
    Υ 20 , 1 ( x , y , z , t ) = 6 c β 1 2 ρ α b 6 c β 1 2 ρ ω b e ω ζ p q 6 c β 1 2 ρ 2 b e ω ζ p q 2 .

Graphical Representation

In this part, we analyze a diverse set of acquired solutions, evaluating their unique features. By investigating the physical structures, we interpret the dynamic nature of these solutions concerning the EIWE (4). Employing the NAE method, this study generates a range of soliton structures. To deepen our comprehension of the dynamic wave behavior exhibited by the acquired solutions, we present graphical depictions in 3D, 2D, and density plots. Additionally, by employing suitable parameter values, we identify various soliton structures, including bright, kink, and dark solitons. In Figure 1a, the solution Υ 2 , 1 is visually presented in three dimensions in the x t -plane, accompanied by its two-dimensional and density plots. The parameter values for this display are as follows: β 1 = 2.4 , β 2 = 1.2 , β 3 = 2.5 , μ = 1.5 , b = 1 , c = 1.5 ω = 0.5 , α = 0.1 , ρ = 0.3 , y = 0.1 , and z = 2 . The 3D graph spans the ranges 5 x 5 and 5 t 5 . The 2D plot is between dependent variable Υ 2 , 1 and independent variable x within the interval 2 x 5 at various time values as t = 1 , 2 , 3 . In Figure 1b, 3D, density, and 2D visualizations of Υ 2 , 1 are presented in the x z -plane with the same parameters as in the previous case. The 3D graph is plotted within the ranges of 5 x 5 and 5 z 5 . The 2D plot illustrates the relationship between the dependent variable Υ 2 , 1 and the independent variable z within the interval 7 z 5 at x = 1 , 2 , 3 . In Figure 1c, a visualization of the same solution is presented in the y z -plane. The 3D graph is plotted within the ranges of 5 y 5 and 5 z 5 . The 2D plot illustrates the relationship between the dependent variable Υ 2 , 1 and the independent variable y within the interval 9 y 6 at z = 1 , 2 , 3 . In this context, the solutions manifest as bright solitons. In nonlinear optics, bright solitons are primary self-localized patterns found in optical fields within dispersive nonlinear mediums, such as optical fibers. They manifest as concentrated intensity peaks amid a uniform background.
Moving to Figure 2a, the solution Υ 4 , 1 is depicted in three dimensions in the x y -plane, along with its 2D and density graphs. The parameter values for this display are as follows: β 1 = 1.4 , β 2 = 1.2 , β 3 = 0.5 , μ = 1.5 , c = 1.5 , b = 1 , ω = 3 , α = 2 , ρ = 2 , t = 0.1 , and z = 2 . The 3D graph spans the ranges 8 x 8 and 6 y 6 . The 2D plot is between dependent variable Υ 4 , 1 and independent variable x within the interval 5 x 2 at y = 1 , 2 , 3 . In Figure 2b, the 3D, density, and 2D visualizations of Υ 4 , 1 are presented in the t z -plane with the same parameters as in the previous case. The 3D graph is plotted within the ranges of 15 t 15 and 5 z 5 . The 2D plot illustrates the relationship between the dependent variable Υ 4 , 1 and the independent variable t within the interval 3 t 8 at z = 1 , 2 , 3 . In Figure 2c, a visualization of same solution is presented in the x z -plane. The 3D graph is plotted within the ranges of 6 x 6 and 2 z 2 . The 2D plot illustrates the relationship between the dependent variable Υ 4 , 1 and the independent variable z within the interval 25 z 12 at x = 1 , 2 , 3 . In this case, the solutions are displayed as kink solitons. In nonlinear optics, kink waves are propagating waves that rise or incline from one asymptotic position to another, eventually stabilizing as they extend to infinity.
Furthermore, in Figure 3a, the solution Υ 6 , 1 is depicted in three dimensions in the x z -plane, along with its 2D and density graphs. The parameter values for this display are as follows: β 1 = 1 , β 2 = 1 , β 3 = 1 , μ = 1.5 , c = 1.5 , b = 1 , ω = 3 , α = 1 , ρ = 1 , t = 0.01 , and y = 1 . The 3D graph spans the ranges 5 x 5 and 5 z 5 . The 2D plot is between dependent variable Υ 6 , 1 and independent variable x within the interval 7 x 5 at z = 1 , 2 , 3 . In Figure 3b, the 3D, density, and 2D visualizations of Υ 6 , 1 are presented in the t z -plane with the same parameters as in the previous case. The 3D graph is plotted within the ranges of 3 t 3 and 3 z 3 . The 2D plot illustrates the relationship between the dependent variable Υ 6 , 1 and the independent variable z within the interval 5 z 7 at t = 1 , 2 , 3 . In Figure 3c, a visualization of the same solution is presented in the x y -plane. The 3D graph is plotted within the ranges of 2 x 2 and 1 y 2 . The 2D plot illustrates the relationship between the dependent variable Υ 6 , 1 and the independent variable x within the interval 4.8 x 1 at y = 1 , 2 , 3 . In this case, the solutions are displayed as dark solitons. Dark solitons attract considerable attention in optics due to their stable transmission. In fiber optics, a stable soliton typically emerges when there is a balance between Kerr nonlinearity and group velocity dispersion (GVD). Within the anomalous dispersion, solitons manifest as bright solitons, while in the normal dispersion, they are transmitted as dark solitons. All solutions exhibit practical applicability in optical communications. Furthermore, these solutions demonstrate the ability to maintain both their shape and velocity during prolonged travel over considerable distances.

4. Bifurcation Phenomena

In the current section, we analyze Equation (4) within the framework of bifurcation analysis [20,21,22]. Upon employing the Galilean transformation to Equation (10), we derive the planar system as follows:
d Γ d ζ = Z , d Z d ζ = σ 0 Γ σ 1 Γ 2 ,
where σ 0 = ( β 1 μ + a β 1 2 + d β 2 2 + e β 1 β 2 + h β 1 β 3 + k β 2 β 3 + m β 3 2 ) c β 1 4 , σ 1 = b c β 1 2 . The Hamiltonian of (34) is expressed as follows:
H ( Γ , Z ) = Z 2 2 + σ 0 Γ 2 2 + σ 1 Γ 3 3 ,
which satisfies
H Γ = Z , H Z = Γ .
In practice, Equation (34) defines a two-dimensional system, and the trajectories within its phase space are governed by the vector field associated with Equation (34). Hence, it is crucial to analyze different phase portraits of (34) by varying its parameters. Equilibria and periodic orbits typically experience changes in their local stability characteristics during a bifurcation. For this reason, we initially pinpoint the critical points P j ( Γ , Z ) (where j = 1, 2) for system (35). These points of equilibrium occur when d Γ d ζ = 0 and d Z d ζ = 0 , yielding two critical points determined by specific parameter values.
P 1 = ( 0 , 0 ) , P 2 = ( σ 0 σ 1 , 0 ) .
Furthermore, the Jacobian of the system will be
J ( Γ , Z ) = 0 1 σ 0 2 σ 1 Γ 0
The determinant and trace of (36) at the critical point P j are indicated by B and T , respectively, and these values are provided as follows:
T = t r a c e ( J ) | P j = 0 , B = d e t ( J ) | P j = σ 0 + 2 σ 1 Γ .
The fixed point is classified as a saddle when B < 0 , a central point when B > 0 and T = 0 , a cusp when B = 0 , and a node if B > 0 and T 2 4 B > 0 .
In the subsequent discussion, we depict symbols used for classifying distinct orbits within the phase portraits of a planner system (34) as follows:
  • A nonlinear periodic trajectory is denoted as N P T ( u * , v * ) .
  • A nonlinear homoclinic trajectory is denoted as N H T ( u * , v * ) .
Here, u * and v * refer to the total fixed points and the overall separation layers encompassed within the trajectory depicted in the phase portrait. This gives rise to various scenarios as follows:
  • Case 1: Let σ 0 > 0 and σ 1 > 0 .
    For β 1 = β 2 = β 3 = k = m = c = b = d = μ = 1 , a = e = h = 1 , σ 0 = 1 , and σ 1 = 1 , system (34) has two fixed points, P 1 = ( 0 , 0 ) and P 2 = ( 1 , 0 ) . In this case, P 1 is center and P 2 is saddle. These points are presented in Figure 4a.
  • Case 2: Let σ 0 < 0 and σ 1 < 0 .
    For β 1 = β 2 = β 3 = k = m = c = d = 1 , μ = a = b = e = h = 1 , σ 0 = 1 , and σ 1 = 1 , system (34) has two fixed points, P 1 = ( 0 , 0 ) and P 2 = ( 1 , 0 ) . In this case, P 1 is saddle and P 2 is center. These points are presented in Figure 4b.
  • Case 3: Let σ 0 < 0 and σ 1 > 0 .
    For β 1 = β 2 = β 3 = k = b = m = c = d = 1 , μ = a = e = h = 1 , σ 0 = 1 , and σ 1 = 1 , system (34) has fixed points as P 1 = ( 0 , 0 ) and P 2 = ( 1 , 0 ) . In this case, P 1 is saddle and P 2 is center. These points are presented in Figure 5a.
  • Case 4: Let σ 0 > 0 and σ 1 < 0 .
    For β 1 = β 2 = β 3 = k = m = c = d = μ = 1 , b = a = e = h = 1 , σ 0 = 1 , and σ 1 = 1 , system (34) has fixed points as P 1 = ( 0 , 0 ) and P 2 = ( 1 , 0 ) . In this case, P 1 is center and P 2 is saddle. These points are presented in Figure 5b.
  • Case 5: Let σ 0 = 0 and σ 1 > 0 .
    For β 1 = β 2 = β 3 = k = b = c = d = μ = h = 1 , m = 2 , a = e = 1 , σ 0 = 0 , and σ 1 = 1 , system (34) has only one fixed point as P 1 = ( 0 , 0 ) . In this case, P 1 is a cuspidal point, and it is presented in Figure 6a.
  • Case 6: Let σ 0 = 0 and σ 1 < 0 .
    For β 1 = β 2 = β 3 = k = c = d = μ = h = 1 , m = 2 , b = a = e = 1 , σ 0 = 0 , and σ 1 = 1 , system (34) has only one fixed point as P 1 = ( 0 , 0 ) . In this case, P 1 is a cuspidal point, and it is presented in Figure 6b.

5. Chaotic Phenomena

In this section, we introduce an outward periodic force into system (34) to examine the dynamics of both quasi-periodic and chaotic phenomena [20,21,22]. The modified system can be expressed in the following manner:
d Γ d ζ = Z , d Z d ζ = σ 0 Γ σ 1 Γ 2 + δ 0 cos ( G ) , d G d ζ = φ ,
The above system is autonomous with G = φ ζ . In the modified system (37), the perturbing force is characterized by two components referred to as δ 0 and φ . In this context, δ 0 represents the intensity, while φ signifies the frequency of the perturbed term introduced into the system (34). We have systematically examined and illustrated the chaotic dynamics of (37) through diverse tools, including phase plots, time plots, and Lyapunov exponents. To delve deeper into this issue, we examine the impact of both the intensity δ 0 and the frequency φ . We keep all other parameters constant at σ 0 = 3.5 and σ 1 = 1 during our analysis.
In Figure 7, a two-dimensional plot and a time plot are presented, deliberately opting for small values for both intensity and frequency, specifically, assigning δ 0 = 0.05 and φ = 0.02 . Under these chosen parameter values, periodic behavior is observed in system (37). Subsequently, a slight increment is made to these parameters, setting them to δ 0 = 0.8 and φ = π . During this adjustment, the system (37) experiences a subtle disturbance and transitions into a quasi-periodic state, as visually depicted in Figure 8. Continuing the analysis, the intensity is further elevated to δ 0 = 3.8 , and the frequency of the perturbing force is increased to φ = 2 π . In this scenario, system (37) displays increased irregularity, revealing a quasi-periodic and chaotic nature, as illustrated in Figure 9. Moving forward, the intensity and frequency are set to δ 0 = 8.5 and φ = 3 π . With these increased factors, system (37) demonstrates heightened randomness, unveiling a chaotic nature, as presented in Figure 10.
The Lyapunov exponent is employed to assess the rate at which nearby trajectories diverge within a dynamical system. It is a scalar metric that signifies the level of chaos within the system. A positive Lyapunov exponent denotes chaotic behavior, while a negative value indicates stability. A higher Lyapunov exponent value corresponds to a more chaotic system, showcasing a faster separation of neighboring trajectories. Subsequently, we visually represent the evolution of these exponents over time to gain insights into the dynamics of the system (37). In Figure 11, we illustrate the Lyapunov exponents obtained against the temporal evolution, aiming to identify the chaotic behavior of the modified dynamical system with the parameters σ 0 = 3.5 , σ 1 = 1 , δ 0 = 8.5 , φ = 3 π and initial condition (0.5, 0.5, 0.5).

6. Multistability

In the current part, our primary emphasis is on investigating the potential presence of multistability [20] within (37). When the system undergoes perturbation, it demonstrates a compelling ability to exhibit multistability. This suggests that it can display the simultaneous coexistence of various dynamic patterns, depending on specific parameter configurations and distinct initial conditions. These dynamic patterns encompass periodicity, quasi-periodicity, and chaos.
To visually depict this phenomenon of multistability, we refer to Figure 12, which illustrates the diverse responses of the system to varying initial values. In Figure 12a the red plot, corresponding to initial values (0.4, 0.3, 0), reveals a periodic nature. Meanwhile, the blue and green plots simultaneously demonstrate quasi-periodic behavior for initial values (−0.5, 0.1, 0) and (−0.03, 0.4, 0). In Figure 12b, the red plot indicates periodic behavior for (−0.5, −0.03, 0), while the blue and green plots exhibit chaotic phenomena at (−0.4, −0.3, 0) and (−0.03, 0.4, 0).

7. Sensitivity Analysis

In this section, we investigate how the suggested equation responds to changes in initial conditions [22,23]. To evaluate the sensitivity of the model, we analyze three distinct sets of initial conditions. The first set, represented by the green curve, corresponds to ( Γ , Z ) = ( 0.01 , 0.01 ) ; the second set, illustrated by the red curve, corresponds to ( Γ , Z ) = ( 0.05 , 0.03 ) ; and the third set, indicated by the pink curve, involves ( Γ , Z ) = ( 0.4 , 0.08 ) . In Figure 13a, two solutions are observed: one for the initial condition ( Γ , Z ) = ( 0.01 , 0.01 ) in green and the other for ( Γ , Z ) = ( 0.05 , 0.03 ) in red. Likewise, in Figure 13b, we present two solutions: ( Γ , Z ) = ( 0.01 , 0.01 ) in green and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink. Figure 13c also displays two solutions: ( Γ , Z ) = ( 0.05 , 0.03 ) in red and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink.
We conducted a comparative analysis using different initial values, specifically (0.01, 0.01), (0.05, 0.03), and (0.4, 0.08), as shown in Figure 13d. It is apparent that even slight modifications in the initial conditions can result in subtle shifts in the outcomes of a dynamic system. Therefore, we conclude that the proposed system demonstrates little sensitivity.

8. Conclusions

In summary, this study investigated the extended integrable wave equation, a widely utilized concept in physics, engineering, and diverse scientific disciplines. Our primary concern was to comprehensively analyze this equation, exploring various aspects including the extraction of Lie symmetries; solitons; an analysis of the phenomena of bifurcation; chaos; and an analysis of the sensitive nature of the suggested equation. Initially, the identification of the translational symmetries of the model was carried out through the utilization of the Lie approach. Then, the NAE method was applied to elucidate the soliton dynamics, and subsequent visualizations of the solutions were created utilizing Mathematica software, including three-dimensional, two-dimensional, and density plots.
Subsequently, an analysis of the dynamical nature of the model was conducted, encompassing various angles such as bifurcation, chaos, and sensitivity. The bifurcation phenomenon was analyzed at critical points within a dynamical system, accompanied by the application of an outward force, which unveiled the emergence of chaotic phenomena. Two-dimensional plots, time plots, multistability, and Lyapunov exponents were presented to illustrate these chaotic trajectories as depicted in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. Furthermore, the sensitivity of the investigated model was executed utilizing the RK4 method. This analysis confirms that the stability of the solution is minimally affected by small changes in initial conditions. The outcomes of our investigation constitute a significant addition to the domain, offering new perspectives on the relevance and applicability of the extended water equation. Going beyond the confines of soliton theory, our findings contribute to a more comprehensive exploration of dynamical systems.

Author Contributions

Conceptualization, M.B.R. and A.J.; methodology, M.B.R. and A.J.; software, J.M. and M.B.R.; validation, F.Z.D., M.B.R., and A.J.; formal analysis, M.B.R. and A.J.; investigation, M.B.R. and A.J.; data curation, F.Z.D. and J.M.; writing—original draft preparation, M.B.R. and A.J.; writing—review and editing, F.Z.D.; visualization, M.B.R. and A.J.; supervision, F.Z.D.; project administration, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This article has been produced with the financial support of the European Union under the REFRESH—Research Excellence For Region Sustainability and High-tech Industries project number CZ.10.03.01/00/22_003/0000048 via the Operational Programme Just Transition. Researchers Supporting Project number (RSPD2024R535), King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kumar, S.; Kumar, A.; Mohan, B. Evolutionary dynamics of solitary wave profiles and abundant analytical solutions to a (3 + 1)-dimensional burgers system in ocean physics and hydrodynamics. J. Ocean. Eng. Sci. 2021, 8, 1–14. [Google Scholar] [CrossRef]
  2. González-Gaxiola, O.; Biswas, A.; Moraru, L.; Alghamdi, A.A. Solitons in Neurosciences by the Laplace-Adomian Decomposition Scheme. Mathematics 2023, 11, 1080. [Google Scholar] [CrossRef]
  3. Kumar, S.; Rani, S.; Mann, N. Diverse analytical wave solutions and dynamical behaviors of the new (2 + 1)-dimensional Sakovich equation emerging in fluid dynamics. Eur. Phys. J. Plus 2022, 137, 1226. [Google Scholar] [CrossRef]
  4. Li, B.Q.; Ma, Y.L. Higher-order breathers and breather interactions for the AB system in fluids. Eur. Phys. J. Plus 2023, 138, 475. [Google Scholar] [CrossRef]
  5. Hosseini, K.; Akbulut, A.; Baleanu, D.; Salahshour, S.; Mirzazadeh, M.; Akinyemi, L. The geophysical KdV equation: Its solitons, complexiton, and conservation laws. GEM-Int. J. Geomath. 2022, 13, 12. [Google Scholar] [CrossRef]
  6. Wazwaz, A.M. Painlevé integrability and lump solutions for two extended (3 + 1)-and (2 + 1)-dimensional Kadomtsev-Petviashvili equations. Nonlinear Dyn. 2023, 111, 3623–3632. [Google Scholar] [CrossRef]
  7. Duran, S.; Karabulut, B. Nematicons in liquid crystals with Kerr Law by sub-equation method. Alex. Eng. J. 2022, 61, 1695–1700. [Google Scholar] [CrossRef]
  8. Vivas-Cortez, M.; Raza, N.; Kazmi, S.S.; Chahlaoui, Y.; Basendwah, G.A. A novel investigation of dynamical behavior to describe nonlinear wave motion in (3 + 1)-dimensions. Results Phys. 2023, 55, 107131. [Google Scholar] [CrossRef]
  9. Yuan, R.R.; Shi, Y.; Zhao, S.L.; Zhao, J.X. The combined KdV-mKdV equation: Bilinear approach and rational solutions with free multi-parameters. Results Phys. 2023, 55, 107188. [Google Scholar] [CrossRef]
  10. Rehman, H.; Said, G.S.; Amer, A.; Ashraf, H.; Tharwat, M.M.; Abdel-Aty, M.; Elazab, N.S.; Osman, M.S. Unraveling the (4 + 1)-dimensional Davey-Stewartson-Kadomtsev-Petviashvili equation: Exploring soliton solutions via multiple techniques. Alex. Eng. J. 2024, 90, 17–23. [Google Scholar] [CrossRef]
  11. Kumar, S.; Ma, W.X.; Kumar, A. Lie symmetries, optimal system and group-invariant solutions of the (3 + 1)-dimensional generalized KP equation. Chin. J. Phys. 2021, 69, 1–23. [Google Scholar] [CrossRef]
  12. Riaz, M.B.; Jhangeer, A.; Martinovic, J.; Kazmi, S.S. Dynamics and Soliton Propagation in a Modified Oskolkov Equation: Phase Plot Insights. Symmetry 2023, 15, 2171. [Google Scholar] [CrossRef]
  13. Wazwaz, A.M. A new (3 + 1)-dimensional Painlevé-integrable Sakovich equation: Multiple soliton solutions. Int. J. Numer. Methods Heat Fluid Flow 2021, 31, 3030–3035. [Google Scholar] [CrossRef]
  14. El-Nabulsi, R.A.; Anukool, W. A family of nonlinear Schrödinger equations and their solitons solutions. Chaos Solitons Fractals 2023, 166, 112907. [Google Scholar] [CrossRef]
  15. Hosseini, K.; Alizadeh, F.; Sadri, K.; Hinçal, E.; Akbulut, A.R.Z.U.; Alshehri, H.M.; Osman, M.S. Lie vector fields, conservation laws, bifurcation analysis, and Jacobi elliptic solutions to the Zakharov–Kuznetsov modified equal-width equation. Opt. Quantum Electron. 2024, 56, 506. [Google Scholar] [CrossRef]
  16. Akinyemi, L. Shallow ocean soliton and localized waves in extended (2 + 1)-dimensional nonlinear evolution equations. Phys. Lett. A 2023, 463, 128668. [Google Scholar] [CrossRef]
  17. Wazwaz, A.M.; Alhejaili, W.; El-Tantawy, S.A. Analytical study on two new (3 + 1)-dimensional Painlevé integrable equations: Kink, lump, and multiple soliton solutions in fluid mediums. Phys. Fluids 2023, 35, 093119. [Google Scholar] [CrossRef]
  18. Talafha, A.M.; Jhangeer, A.; Kazmi, S.S. Dynamical analysis of (4 + 1)-dimensional Davey Srewartson Kadomtsev Petviashvili equation by employing Lie symmetry approach. Ain Shams Eng. J. 2023, 14, 102537. [Google Scholar] [CrossRef]
  19. Lathrop, D. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  20. Kazmi, S.S.; Jhangeer, A.; Raza, N.; Alrebdi, H.I.; Abdel-Aty, A.H.; Eleuch, H. The Analysis of Bifurcation, Quasi-Periodic and Solitons Patterns to the New Form of the Generalized q-Deformed Sinh-Gordon Equation. Symmetry 2023, 15, 1324. [Google Scholar] [CrossRef]
  21. Rafiq, M.H.; Raza, N.; Jhangeer, A. Dynamic study of bifurcation, chaotic behavior and multi-soliton profiles for the system of shallow water wave equations with their stability. Chaos Solitons Fractals 2023, 171, 113436. [Google Scholar] [CrossRef]
  22. Hosseini, K.; Hinçal, E.; Ilie, M. Bifurcation analysis, chaotic behaviors, sensitivity analysis, and soliton solutions of a generalized Schrödinger equation. Nonlinear Dyn. 2023, 111, 17455–17462. [Google Scholar] [CrossRef]
  23. Arshed, S.; Raza, N.; Rahman, R.U.; Butt, A.R.; Huang, W.H. Sensitive behavior and optical solitons of complex fractional Ginzburg-Landau equation: A comparative paradigm. Results Phys. 2021, 28, 104533. [Google Scholar] [CrossRef]
  24. Kumar, S.; Kumar, D.; Kumar, A. Lie symmetry analysis for obtaining the abundant exact solutions, optimal system and dynamics of solitons for a higher-dimensional Fokas equation. Chaos Solitons Fractals 2021, 142, 110507. [Google Scholar] [CrossRef]
  25. Liu, H.; Li, J. Lie symmetry analysis and exact solutions for the short pulse equation. Nonlinear Anal. Theory Methods Appl. 2009, 71, 2126–2133. [Google Scholar] [CrossRef]
Figure 1. Depiction of the solution Υ 2 , 1 across various planes.
Figure 1. Depiction of the solution Υ 2 , 1 across various planes.
Symmetry 16 00608 g001
Figure 2. Depiction of the solution Υ 4 , 1 across various planes.
Figure 2. Depiction of the solution Υ 4 , 1 across various planes.
Symmetry 16 00608 g002
Figure 3. Depiction of the solution Υ 6 , 1 across various planes.
Figure 3. Depiction of the solution Υ 6 , 1 across various planes.
Symmetry 16 00608 g003
Figure 4. Phase plots for (34), when (a) σ 0 > 0 and σ 1 > 0 , and (b) σ 0 < 0 and σ 1 < 0 .
Figure 4. Phase plots for (34), when (a) σ 0 > 0 and σ 1 > 0 , and (b) σ 0 < 0 and σ 1 < 0 .
Symmetry 16 00608 g004
Figure 5. Phase plots for (34), when (a) σ 0 < 0 and σ 1 > 0 , and (b) σ 0 > 0 and σ 1 < 0 .
Figure 5. Phase plots for (34), when (a) σ 0 < 0 and σ 1 > 0 , and (b) σ 0 > 0 and σ 1 < 0 .
Symmetry 16 00608 g005
Figure 6. Phase plots for (34), when (a) σ 0 = 0 and σ 1 > 0 , and (b) σ 0 = 0 and σ 1 < 0 .
Figure 6. Phase plots for (34), when (a) σ 0 = 0 and σ 1 > 0 , and (b) σ 0 = 0 and σ 1 < 0 .
Symmetry 16 00608 g006
Figure 7. Detection of chaotic phenomena in the perturbed system (37).
Figure 7. Detection of chaotic phenomena in the perturbed system (37).
Symmetry 16 00608 g007
Figure 8. Detection of chaotic phenomena in the perturbed system (37).
Figure 8. Detection of chaotic phenomena in the perturbed system (37).
Symmetry 16 00608 g008
Figure 9. Detection of chaotic phenomena in the perturbed system (37).
Figure 9. Detection of chaotic phenomena in the perturbed system (37).
Symmetry 16 00608 g009
Figure 10. Detection of chaotic phenomena in the perturbed system (37).
Figure 10. Detection of chaotic phenomena in the perturbed system (37).
Symmetry 16 00608 g010
Figure 11. Detection of chaotic phenomena in the system (37) via Lyapunov exponent.
Figure 11. Detection of chaotic phenomena in the system (37) via Lyapunov exponent.
Symmetry 16 00608 g011
Figure 12. Detection of chaos for the system (37) through multistability analysis with σ 0 = 3.5 , σ 1 = 1 , δ 0 = 8.5 , and φ = 3 π .
Figure 12. Detection of chaos for the system (37) through multistability analysis with σ 0 = 3.5 , σ 1 = 1 , δ 0 = 8.5 , and φ = 3 π .
Symmetry 16 00608 g012
Figure 13. Sensitivity analysis across various initial values. (a) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green and ( Γ , Z ) = ( 0.05 , 0.03 ) in red. (b) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink. (c) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.05 , 0.03 ) in red and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink. (d) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green, ( Γ , Z ) = ( 0.05 , 0.03 ) in red, and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink.
Figure 13. Sensitivity analysis across various initial values. (a) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green and ( Γ , Z ) = ( 0.05 , 0.03 ) in red. (b) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink. (c) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.05 , 0.03 ) in red and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink. (d) Sensitivity analysis of the system (34) for initial values ( Γ , Z ) = ( 0.01 , 0.01 ) in green, ( Γ , Z ) = ( 0.05 , 0.03 ) in red, and ( Γ , Z ) = ( 0.4 , 0.08 ) in pink.
Symmetry 16 00608 g013
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

Riaz, M.B.; Jhangeer, A.; Duraihem, F.Z.; Martinovic, J. Analyzing Dynamics: Lie Symmetry Approach to Bifurcation, Chaos, Multistability, and Solitons in Extended (3 + 1)-Dimensional Wave Equation. Symmetry 2024, 16, 608. https://doi.org/10.3390/sym16050608

AMA Style

Riaz MB, Jhangeer A, Duraihem FZ, Martinovic J. Analyzing Dynamics: Lie Symmetry Approach to Bifurcation, Chaos, Multistability, and Solitons in Extended (3 + 1)-Dimensional Wave Equation. Symmetry. 2024; 16(5):608. https://doi.org/10.3390/sym16050608

Chicago/Turabian Style

Riaz, Muhammad Bilal, Adil Jhangeer, Faisal Z. Duraihem, and Jan Martinovic. 2024. "Analyzing Dynamics: Lie Symmetry Approach to Bifurcation, Chaos, Multistability, and Solitons in Extended (3 + 1)-Dimensional Wave Equation" Symmetry 16, no. 5: 608. https://doi.org/10.3390/sym16050608

APA Style

Riaz, M. B., Jhangeer, A., Duraihem, F. Z., & Martinovic, J. (2024). Analyzing Dynamics: Lie Symmetry Approach to Bifurcation, Chaos, Multistability, and Solitons in Extended (3 + 1)-Dimensional Wave Equation. Symmetry, 16(5), 608. https://doi.org/10.3390/sym16050608

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