Chaotic Phenomena, Sensitivity Analysis, Bifurcation Analysis, and New Abundant Solitary Wave Structures of The Two Nonlinear Dynamical Models in Industrial Optimization

: In this research, we discussed the different chaotic phenomena, sensitivity analysis, and bifurcation analysis of the planer dynamical system by considering the Galilean transformation to the Lonngren wave equation (LWE) and the (2 + 1)-dimensional stochastic Nizhnik–Novikov–Veselov System (SNNVS). These two important equations have huge applications in the fields of modern physics, especially in the electric signal in data communication for LWE and the mechanical signal in a tunnel diode for SNNVS. A different chaotic nature with an additional perturbed term was also highlighted. Concerning the theory of the planer dynamical system, the bifurcation analysis incorporating phase portraits of the dynamical systems of the declared equations was performed. Additionally, a sensitivity analysis was used to monitor the sensitivity of the mentioned equations. Also, we extracted new, abundant solitary wave structures with the graphical phenomena of the mentioned nonlinear mathematical models. By conducting an expansion method on the abovementioned equations, we generated three types of soliton structures, which are rational function, trigonometric function, and hyperbolic function. By simulating the 3D, contour, and 2D graphs of these obtained solitons, we scrutinized the behavior of the waves affecting the nonlinear terms. The figures show that the solitary waves obtained from LWE are efficient in analyzing electromagnetic wave signals in the cable lines, and the solitary waves from SNNVS are essential in any stochastic system like a sound wave. Moreover, by taking some values of the parameters, we found some interesting soliton shapes, such as compaction soliton, singular periodic solution, bell-shaped soliton, anti-kink-shaped soliton, one-sided kink-shaped soliton, and some flat kink-shaped solitons, etc. This article will have a great impact on nonlinear science due to the new solitary wave structures with different complex phenomena, sensitivity analysis, and bifurcation analysis.


Introduction
Nonlinear partial differential equations (NPDEs) make an indispensable impact on giving the model for characterizing the complex natural occurrence of many scientific and engineering phenomena, namely, optical fibers, mechanics, plasma physics, reaction models, chemical kinematics, fluid dynamics, relativistic physics, biological science, and so on .It is highly important to investigate the progressive wave-like solution for the best perception of NLPDEs and their application in real life.Lately, different kinds of techniques have been exhibited for generating numerical and analytical demonstrations by many experts, such as the Riccati equation method [1], the F-expansion technique [2,3], the auxiliary equation method [4,5], the Jacobi elliptic function method [6,7], the direct algebraic function technique [8,9], the Cole-Hopf conversion technique [10,11], the tanhfunction method [12,13], the Backlund transform technique [14,15], the Hirota's bilinear technique [16][17][18], the exp(−ϕ(ξ))-expansion method [19,20], the generalized Kudryashov method [21,22], the homotopy exploration technique [23], the homogeneous balance technique [24][25][26], the variational iteration method [27], the sine cosine algorithm [28], and the G ′ G ′ +G+A -expansion technique [29-31].In addition, the (G ′ /G)-expansion technique was introduced by Wang et al. for describing the outcomes of NLPDEs [32].Many investigators have applied this technique in their investigative works [33][34][35][36][37][38][39][40].In addition to these, many investigators aspired to look for other alternatives that are more diligent and better than the former technique.Recently, Iqbal et al. [41] and Miah et al. [42] have applied the (G ′ /G, 1/G)-expansion technique to attain new soliton solutions for the mZK equation and the Gerdjikov-Ivanov equation and have applied the abundant closed-form wave solution to some NLEEs in mathematical physics, respectively.Moreover, many scientists [43][44][45][46][47][48] have chosen it as a mechanism to obtain the solution for NLPDEs.In this program, we detail and engage the (G ′ /G, 1/G)-expansion method to analyze the LWE and the (2 + 1)-dimensional SNNVS of the nonlinear electrical transmission line and the propagation of the stochastic systems, which have not been yet been tried with our method.
Lately, a great deal of research has been completed on nonlinear systems.Many disciplines, especially engineering, telecommunications, and ecology alongside economics, have explored bifurcation as well as chaotic nature, which is an appealing nonlinear phenomenon.We can better understand the process in which systems move from a stable state to an unstable state by having an in-depth knowledge of bifurcations.When the behavior of a dynamical system depends extensively on the initial condition, it is categorized as chaotic.Chaos theory, which has implications in ecology, weather forecasting, and financial markets, emphasizes the innate difficulty and variability of dynamic systems.Moreover, the Runge-Kutta technique is exercised for the sensitivity analysis of a dynamical system.The entire analysis provides insightful information on the suggested dynamical systems, which is amazing and useful for engineering and mathematical physics, as well as for developing a better comprehension of real-world occurrences and supporting control alongside anticipation.
The novelty of our research work is that our obtained results from the LWE, also known as the wave equation in general, can give the fundamental concept in physics and engineering that describes the behavior of waves.Also, our obtained solutions can be used in numerous applications across various fields, including acoustics; the attained solutions of the LWE can be used to study the propagation of sound waves in different mediums.This is crucial in designing acoustic systems in concert halls, auditoriums, and soundproofing materials.In seismology, seismologists can use these solutions to analyze seismic waves generated by earthquakes.This helps in understanding the Earth's interior structure and predicting the impact of seismic events.In electromagnetics, our obtained results can be applied to describe the behavior of electromagnetic waves, including light, radio waves, and microwaves.Our results can also be used in designing antennas, optical fibers, and other communication systems.In engineering, engineers may use our results from the LWE to analyze and design structures subjected to dynamic loads, such as vibrations and oscillations.This is crucial in fields like civil engineering, mechanical engineering, and aerospace engineering.In medical imaging, techniques like ultrasound and MRI rely on the principles of wave propagation.Understanding the wave equation is essential in developing and interpreting the images produced by these medical imaging modalities.In oceanography, oceanographers may use the results to study ocean waves and currents, which are essential for understanding climate patterns, coastal erosion, and navigation safety.In quantum mechanics, the results can describe the behavior of quantum particles such as electrons and photons.This is fundamental for understanding the behavior of matter at the atomic and subatomic levels.In optics, these wave solutions are crucial for understanding phenomena like diffraction, interference, and the polarization of light.This understanding can be applied in designing optical instruments, lenses, and imaging systems.In material science, understanding how waves propagate through materials is crucial for studying properties like elasticity, conductivity, and thermal conductivity.This knowledge is applied in various industries, including construction, electronics, and aerospace.In geophysics, geophysicists can use our results to study various phenomena such as seismic exploration for oil and gas, groundwater mapping, and environmental monitoring.
These are just a few examples of how the obtained results can be applied in real-life scenarios across different disciplines.The equation's versatility and fundamental nature make it a cornerstone of modern physics and engineering.The LWE is of the following form [49,50], with two arbitrary constants α and β.The (2 + 1)-dimensional SNNVS is a mathematical model that describes certain physical phenomena in a specific mathematical framework.Our obtained results from the SNNVS have potential applications in various fields of science and technology.
In nonlinear optics, the obtained solutions can be used to describe the behavior of optical solitons in certain nonlinear optical media.Understanding and controlling solitons are important for developing technologies such as optical communication systems and fiber optics.In plasma physics, the obtained solutions may have applications in modeling certain behaviors in plasma physics, such as wave propagation and turbulence in plasmas.This is relevant in fields like nuclear fusion research and astrophysics.In fluid dynamics, the equation and its results governing fluid flow often exhibit nonlinear behavior similar to that described by the SNNVS.Understanding such behavior is crucial in various engineering applications, including weather prediction, aerodynamics, and oceanography.The obtained solutions can be related to certain aspects of quantum field theory and integrable systems.While the direct applications in everyday life may not be obvious, advancements in theoretical physics can lead to technological breakthroughs in the future.In mathematical physics research, the study of our soliton solutions to the SNNVS contributes to theoretical developments in mathematical physics.While this may not have immediate practical applications, it advances our understanding of fundamental physical principles, which can eventually lead to technological innovations.In materials science, understanding the mathematical properties of our results can provide insights into the behavior of complex materials, such as ferroelectric or ferromagnetic materials, which have applications in electronic devices, sensors, and data storage.In biophysics, nonlinear dynamics play a role in modeling certain biological processes, such as the propagation of nerve impulses or the dynamics of biochemical reactions.Insights gained from studying nonlinear systems like the obtained solutions of SNNVS can contribute to understanding biological phenomena.
Overall, while the obtained solutions of (2 + 1)-dimensional SNNVS may not have direct applications in everyday life, they are a valuable mathematical model that contributes to our understanding of complex physical phenomena and has potential applications in various scientific and technological fields.Now, the (2 + 1)-dimensional SNNVS with the Ito sense of augmentative noise [51,52] is as follows: This paper has been arranged in this fashion: In Section 1, we give a brief introduction.We narrate the expansion method in Section 2. In Section 3, the applications of the LWE for chaotic, sensitivity, bifurcation, and soliton solutions are introduced.In Section 4, the applications of the SNNVS for chaotic, sensitivity, bifurcation, and soliton solutions are introduced.We give the graphical declaration in Section 5.In Section 6, we give the conclusion of this article.
Step 2. We write a solution in Equation (15) in the function of Q(τ) and R(τ) as the following, and the parameters a m (m = 0, 1, 2, . . ., T), b m (m = 1, 2, . . ., T), c, λ, and l will be gained through step 3.The balance number is obtained by applying the homogenous balance rules between the heights power nonlinear term and the heights differential term.In Equation ( 16), we put the obtained balance number, and after inserting it in Equation ( 15) and then using Equations ( 6), ( 8), (10), and ( 12), we change the left part of Equation ( 15) in Q and R, where the degree of Q is not greater than one and take it from zero up to any integral value.
Step 3. By equating the same index of the expression from the two parts, we obtain a pattern of algebraic solutions in a m, b m , c, λ, l, c 1 , and c 2 .Now, using computational software (we use Mathematica 13.1), we gain the results of obtaining an algebraic system.By putting the obtained results in Equation ( 16), we attain the outcomes of Equation (15).Now, by setting τ = (mx + ny − pt) in the outcomes of Equation ( 15), we attain the results of Equation (13).For more details of the mentioned method, visit [41][42][43][44][45][46][47][48].

Application of the LWE
In this section, the (G ′ /G, 1/G)-expansion method is engaged to extract the wave soliton solutions and to consider the following transformation, where 'c ′ implies the wave velocity.Applying this transformation, we convert the LWE presented in Equation (1) to the following ODE, the prime implicates the derivative regarding τ.

Chaotic Analysis of LWE
Here, we adapt Equation ( 18) within the planer dynamical system by considering the Galilean transformation [53][54][55][56], The perturbation term ζ 1 cos(ζ 2 t), where ζ 1 stands for amplitude and ζ 2 stands for the frequency, is introduced in the dynamical system of Equation ( 19) so that we can understand the chaotic character of the announced equation, In this instance, the chaotic character of the above dynamical system is graphed for many specific values of the parameters.In Figure 1a 1a,b, we find out complex dynamics, and in Figure 1c, we have the limit cycle.In Figure 2a, we observe ringlet dynamics alongside Figure 2b as well as Figure 2c, which also depict strange dynamics.In addition, the limit cycle is displayed in Figure 3a; likewise, ringlet dynamics are also displayed in Figure 3b,c.Some strong and clear chaotic phenomena are designed that may be advantageous in mathematical physics and modern engineering.

𝑅 𝑊 𝑢 𝑡
. ( In this instance, the chaotic character of the above dynamical system is graphed for many specific values of the parameters.In Figure 1a 2.In Figure 1a,b, we find out complex dynamics, and in Figure 1c, we have the limit cycle.In Figure 2a, we observe ringlet dynamics alongside Figure 2b as well as Figure 2c, which also depict strange dynamics.In addition, the limit cycle is displayed in Figure 3a; likewise, ringlet dynamics are also displayed in Figure 3b,c.Some strong and clear chaotic phenomena are designed that may be advantageous in mathematical physics and modern engineering.
In this instance, the chaotic character of the above dynamical system is graphed for many specific values of the parameters.In Figure 1a 2.In Figure 1a,b, we find out complex dynamics, and in Figure 1c, we have the limit cycle.In Figure 2a, we observe ringlet dynamics alongside Figure 2b as well as Figure 2c, which also depict strange dynamics.In addition, the limit cycle is displayed in Figure 3a; likewise, ringlet dynamics are also displayed in Figure 3b,c.Some strong and clear chaotic phenomena are designed that may be advantageous in mathematical physics and modern engineering.

Sensitivity Analysis of LWE
In this subsection, our aim is to observe the sensitivity nature of the dynamic system stated in Equation ( 19).To complete this activity, we resolve the ensuing dynamical system in Equation ( 19) with the aid of the Runge-Kutta technique.The particular values of the parameters  = 5,  = 5, and  = 0.3 are picked up for this instance.Furthermore, we assume the initial condition of the dynamical system as follows [53][54][55][56]: (i) (0) = 0.1 and  = 0;

Sensitivity Analysis of LWE
In this subsection, our aim is to observe the sensitivity nature of the dynamic system stated in Equation (19).To complete this activity, we resolve the ensuing dynamical system in Equation (19) with the aid of the Runge-Kutta technique.The particular values of the parameters α = 5, β = 5, and c = 0.3 are picked up for this instance.Furthermore, we assume the initial condition of the dynamical system as follows [53][54][55][56]: The outcome is illustrated in Figure 4 for the abovementioned initial conditions, in which the blue curve represents the dynamics of class u, and the orange curve signifies the dynamics of class R. It is demonstrated in Figure 4 that small alterations to the initial states have a major impact on the dynamical system.

Sensitivity Analysis of LWE
In this subsection, our aim is to observe the sensitivity nature of the dynamic syste stated in Equation (19).To complete this activity, we resolve the ensuing dynamical sy tem in Equation ( 19) with the aid of the Runge-Kutta technique.The particular values the parameters  5,  5, and  0.3 are picked up for this instance.Furthermore, w assume the initial condition of the dynamical system as follows [53][54][55][56] The outcome is illustrated in Figure 4 for the abovementioned initial conditions, which the blue curve represents the dynamics of class , and the orange curve signifi the dynamics of class .It is demonstrated in Figure 4 that small alterations to the init states have a major impact on the dynamical system.

Bifurcation and Phase Portrait Analysis of LWE
The planer dynamical system in Equation ( 19) yields the subsequent Hamiltonian function as follows [53][54][55][56]: Here, the Hamiltonian constant is represented by  1 .We have the equilibrium positions at (0, 0) and − W 2 W 1 , 0 , and we subsequently solve the next system, One can evaluate the Jacobian of the dynamical system in Equation ( 19): From the Planar dynamical systems theory, we have the following circumstances [53][54][55][56]: 1.
In the succeeding part, we obtain many situations that are discussed in depth for various conceivable parameter selections.
One can evaluate the Jacobian of the dynamical system in Equation ( 19): From the Planar dynamical systems theory, we have the following circumstances [53][54][55][56]: 1.The equilibrium position ,  denotes a saddle point, whereas  1 ,  0; 2. The equilibrium position ,  stands for a center point, whereas  1 ,  0; 3. The equilibrium position ,  characterizes a cuspid point, whereas  1 ,  0. In the succeeding part, we obtain many situations that are discussed in depth for various conceivable parameter selections.
Instance 1:  0 and  0 In Figure 5a, the equilibrium position 0, 0 implies a saddle point that, together with the equilibrium position 0.9375, 0 , implies a center point by selecting the parameters  1,  1, and  4. Instance 2:  0 and  0 In Figure 5b, the equilibrium position 0, 0 signifies a saddle point alongside the equilibrium position 1.7188, 0 , which signifies a center point upon choosing the value of the parameters  5,  2, and  0.

Soliton Solutions of LWE
In Equation ( 18), we apply the balance rule and obtain the balance number  2. After placing this in Equation ( 16), we obtain Instance 2: W 1 > 0 and W 2 > 0 In Figure 5b, the equilibrium position (0, 0) signifies a saddle point alongside the equilibrium position (−1.7188, 0),which signifies a center point upon choosing the value of the parameters α = 5, β = −2, and c = 0.8.

Soliton Solutions of LWE
In Equation (18), we apply the balance rule and obtain the balance number T = 2.After placing this in Equation ( 16), we obtain where the functions Q(τ) and R(τ) are stated in Equations ( 5) and (6).According to the signs of λ, the three kinds of basic outcomes of Equation ( 18) are as follows: Case-1.For λ > 0.

Application of the SNNVS
Here, we assign the following transformation: Here, τ = mx + ny − pt and A, B, C are deductive functions with the non-zero constants m, n, p, and 'σ' is the noise stability parameter.Placing Equation (39) into Equations ( 2) and (3), we obtain, , where σ 2 2 t is the it ô term.Now, we have the characteristic relations, which are given below: By integrating (41) with regard to τ and taking the constant of integration as equal to zero, we obtain mA = nB and nA = mC.

Chaotic Analysis of SNNVS
Here, we utilize the Galilean transformation for Equation (45), which supplies us with the ensuing planer dynamical system [53][54][55][56], wherein To accomplish the chaotic nature of the declared equation, the perturbation term ζ 3 cos(ζ 4 t), having the magnitude ζ 3 and the frequency ζ 4 , is added into the dynamical system declared in Equation (46): We illustrate the chaotic character of the advocated dynamical system along with the abundant values of the parameters.In Figure 6a, s wherein   ,   .To accomplish the chaotic nature of the declared equation, the perturbation term     , having the magnitude  and the frequency  , is added into the dynamical system declared in Equation ( 46):

𝑅 𝑊 𝐴 𝑡
. (47) We illustrate the chaotic character of the advocated dynamical system along with the abundant values of the parameters.In Figure 6a 2.Now, it is evident from the demonstrated figure that complex dynamics are presented in Figure 6a, and also that the limit cycles are presented in Figure 6b,c.Moreover, strange dynamics are offered in Figure 7a,c; also, in Figure 7c, one can see the periodic dynamics.In addition, Figure 8a,b supplies us the surprising nature along with Figure 8c, which gives us the periodic dynamics.From the above observation, we are confident that this analysis detects insights into the dynamical system and that it might be valuable as well as advantageous in modern physics, advanced mathematics, and innovative engineering.

Sensitivity Analysis of SNNVS
Here, the sensitivity of the dynamical system in Equation ( 46) is reviewed.To finalize its sensitivity, we choose the definite value of the parameters  2,  2,  1,  1,  1,  1,  1, and  1.We also set the initial conditions for the selected dynamical system as follows [53][54][55][56]: Now, it is evident from the demonstrated figure that complex dynamics are presented in Figure 6a, and also that the limit cycles are presented in Figure 6b,c.Moreover, strange dynamics are offered in Figure 7a,c; also, in Figure 7c, one can see the periodic dynamics.In addition, Figure 8a,b supplies us the surprising nature along with Figure 8c, which gives us the periodic dynamics.From the above observation, we are confident that this analysis detects insights into the dynamical system and that it might be valuable as well as advantageous in modern physics, advanced mathematics, and innovative engineering.

Sensitivity Analysis of SNNVS
Here, the sensitivity of the dynamical system in Equation ( 46) is reviewed.To finalize its sensitivity, we choose the definite value of the parameters s = 2, q = 2, p = 1, r = 1, k = 1, m = 1, n = 1, and β = 1.We also set the initial conditions for the selected dynamical system as follows [53][54][55][56]: (i) A(0) = 0.2 and R = 0; (ii) A(0) = 0 and R = 0.2; (iii) A(0) = 0.3 and R = 0; (iv) A(0) = 0 and R = 0.3 In Figure 9, the sensitivity nature of the proposed system for numerous situations of initial positions are depicted, where the dynamics of class A and the dynamics of class R are symbolized by blue curves and orange curves, respectively.From a close inspection of Figure 9, we realize a major change in the dynamical system's nature, where one could employ a tiny alteration in the initial position.
Herein, the Hamiltonian constant is identified by  .Now, the equilibrium stan points 0, 0 , , 0 can be determined after the resolution of the dynamical system
The equilibrium standpoint (A, R) refer to a saddle point, while D 2 (A, R) < 0; 2.
Instance 3:  0 and  0 In Figure 10c, the equilibrium standpoint 0, 0 exhibits a center point, and the equilibrium standpoint 1.3333, 0 exhibits a saddle point by picking the parameters as

Soliton Solutions of SNNVS
Balancing the principle in Equation ( 45) is as follows: where the functions   and  are explained in Equations ( 5) and (6).Subject to the signs of λ, we obtain three cases: Case-1.For positive values, i.e.,  0: By doing two times differentiation in Equation (48) and engaging Equations ( 6) and (8), we convert the back side of Equation (48) in  and  variables.Now, balancing the coefficients in the acquired polynomial and taking zero, we gain a system of equations in the constants  ,  ,  ,  ,  , , and .After finishing this action and employing the program Mathematica, we have two sets of results: Set I:  Instance 2: W 3 < 0 and W 4 < 0 In Figure 10b, the equilibrium standpoint (0, 0) displays a center point together with the equilibrium standpoint (−0.5556, 0), which displays a saddle point by opting the value of the parameters as m = 1, n = −1, k = 2, r = −1, s = 1, p = −5, and q = 1.

Soliton Solutions of SNNVS
Balancing the principle in Equation ( 45) is as follows: where the functions Q(τ) and R(τ) are explained in Equations ( 5) and (6).Subject to the signs of λ, we obtain three cases: Case-1.For positive values, i.e., λ > 0: By doing two times differentiation in Equation (48) and engaging Equations ( 6) and ( 8), we convert the back side of Equation (48) in Q and R variables.Now, balancing the coefficients in the acquired polynomial and taking zero, we gain a system of equations in the constants a 0 , a 1 , a 2 , b 1 , b 2 , λ, and c.After finishing this action and employing the program Mathematica, we have two sets of results: Set I: where λ = M 1 ; M 1 < 0. Taking c 2 = 0 and l = 0 but c 1 ̸ = 0, where M 1 < 0.
Remark 1. Forλ = 0, we cannot find any convenient solution using this method.

Graphical Representation of Soliton Solutions
In this section, we give the picturesque description of the gained solutions of the LWE and the (2 + 1)-dimensional SNNVS.The attained results are trigonometric, hyperbolic, and rational function types.The acquired results of Equations ( 25), ( 26), ( 28), ( 29), ( 52), (53), and ( 56) are of trigonometric function; the results in Equations ( 33), ( 35), (59), and (61) are of hyperbolic function; and the results in Equations ( 37) and ( 38) are of algebraic function types.For the values taken as constant, we obtain the graph of the above explicit solutions in 3D, 2D, and contour forms as given in Figures 11-17.
In Figure 11, we plot Equation ( 25) by taking the values α = −0.01,β = 2, and c = 1 in 3D and 2D (t = 0, 1, 5, 10 with the range −5 ≤ x ≤ 5) shapes, along with their contour shapes, by giving a compaction soliton within the range x ∈ [0, 2], t ∈ [0, 2].From the contour plot, we observe that the energy density is distributed uniformly in the region x ∈ [1, 2], t ∈ [0, 5], and from 2D shape we see that as time increases, figures with equal amplitudes and equal widths are displaced from left to right.
with their contour shapes, by giving a compaction soliton within the range  ∈ 5, 5 ,  ∈ 0, 3 .From the contour plot, we observe that the energy density is at maximum in the region  ∈ 4, 2 ,  ∈ 1, 3 , and from 2D shape, we see that the amplitude and width of the figures decrease as time increases.(i) (ii) (iii)  with their contour shapes, by giving a compaction soliton within the range  ∈ 5, 5 ,  ∈ 0, 3 .From the contour plot, we observe that the energy density is at maximum in the region  ∈ 4, 2 ,  ∈ 1, 3 , and from 2D shape, we see that the amplitude and width of the figures decrease as time increases.

Conclusions
In the recent study, we successfully investigated some important analyses, chaotic phenomena, sensitivity, and bifurcation to the LWE and the SNNVS.In chaotic phenomena, we plotted different physically applicable patterns such as limit cycle, complex dynamics, ringlet dynamics, strange dynamics, periodic dynamics, etc.For both mathematical models, we applied sensitivity analyses, and we discovered that these systems' outcomes are very sensitive to their initial states.Bifurcation analysis was exploited to figure out the equilibrium points of the relevant systems.Also, we extracted new, abundant wave soliton solutions of the mentioned equations by applying the double-variable expansion method.The graphical demonstration of the soliton solutions implies that the formations of the attained solutions are of compaction soliton, parabolic-shaped soliton, singular periodic solution, bell-shaped soliton, anti-kink-shaped soliton, one-sided kinkshaped soliton, flat kink-shaped solitons, etc.Our obtained results may be used to explore the naturalistic phenomena subject to the nonlinear partial differential equation.Finally, we have checked our results by back substitute through Mathematica software.The obtained solutions might be important for exposing the naturalistic occurrence arising in science and engineering, which can amplify the impact on posterior research work.This In Figure 12, we plot Equation ( 29) by taking the values α = 2, β = 1, and c = 1 in 3D and 2D (t = 0, 5, 10 with the range −5 ≤ x ≤ 5) shapes, along with their contour shapes, by giving a singular periodic solution within the range x ∈ [−5, 5], t ∈ [0, 5].From the contour plot, we observe that the energy density is distributed uniformly in the vertical region with x ∈ [2.5, 3], and from 2D shape, we see that as time increases, figures with equal amplitudes and equal widths are displaced from left to right.
In Figure 13, we plot Equation ( 28) by taking the values α = 2, β = 1, and c = 1 in 3D and 2D (t = 0, 5, 10 with the range −10 ≤ x ≤ 10) shapes, along with their contour shapes, by giving a bell-shaped soliton within the range x ∈ [−5, 10], t ∈ [0, 5].From the contour plot, we observe that the energy density is distributed discontinuously in the inclined region x ∈ [−2, 2], and from 2D shape, we see that as time increases, figures with equal amplitudes and equal widths are displaced from left to right.
In Figure 14, we plot Equation (35) by taking the values α = 1, β = 1, and c = 0.1 in 3D and 2D (t = 0, 5, 10 with the range −5 ≤ x ≤ 5) shapes, along with their contour shapes, by giving an anti-kink-shaped soliton within the range x ∈ [0, 5], t ∈ [0, 5].From the contour plot, we observe that the energy density is distributed uniformly except the vertical strip of the region x ∈ [0, 0.1], t ∈ [0, 5], and from 2D shape, we see that the amplitude remains the same, but width of the figures increases as time increases.

Conclusions
In the recent study, we successfully investigated some important analyses, chaotic phenomena, sensitivity, and bifurcation to the LWE and the SNNVS.In chaotic phenomena, we plotted different physically applicable patterns such as limit cycle, complex dynamics, ringlet dynamics, strange dynamics, periodic dynamics, etc.For both mathematical models, we applied sensitivity analyses, and we discovered that these systems' outcomes are very sensitive to their initial states.Bifurcation analysis was exploited to figure out the equilibrium points of the relevant systems.Also, we extracted new, abundant wave soliton solutions of the mentioned equations by applying the double-variable expansion method.The graphical demonstration of the soliton solutions implies that the formations of the attained solutions are of compaction soliton, parabolic-shaped soliton, singular periodic solution, bell-shaped soliton, anti-kink-shaped soliton, one-sided kink-shaped soliton, flat kink-shaped solitons, etc.Our obtained results may be used to explore the naturalistic phenomena subject to the nonlinear partial differential equation.Finally, we have checked our results by back substitute through Mathematica software.The obtained solutions might be important for exposing the naturalistic occurrence arising in science and engineering, which can amplify the impact on posterior research work.This article will have great impact on nonlinear science due to the new solitary wave structures with different complex phenomena, sensitivity analysis, and bifurcation analysis.

Figure 1 .Figure 2 .
Figure 1.Chaotic patterns of the offered dynamical system with specific values of the parameters together with   2 ⁄ .

Figure 1 .
Figure 1.Chaotic patterns of the offered dynamical system with specific values of the parameters together with ζ 2 = π/2.

Figure 1 .Figure 2 .
Figure 1.Chaotic patterns of the offered dynamical system with specific values of the parameters together with   2 ⁄ .

Figure 2 .
Figure 2. Chaotic patterns of the offered dynamical system with several values of the parameters together with ζ 2 = π.

Figure 1 .Figure 2 .Figure 3 .
Figure 1.Chaotic patterns of the offered dynamical system with specific values of the parameters together with  =  2 ⁄ .

Figure 3 .
Figure 3. Chaotic patterns of the offered dynamical system with numerous definite values of the parameters together with ζ 2 = 2π.

Figure 3 .
Figure 3. Chaotic patterns of the offered dynamical system with numerous definite values of t parameters together with  2.

Figure 4 .
Figure 4. Sketch of the sensitivity analysis of the proposed dynamical system upon considering se eral values of the parameters alongside some initial conditions: (a)  0 0.1 and  0 , (  0 0 and  0.1, (c)  0 0.15 and  0 and (d)  0 0 and  0.15.

8 .Figure 5 .
Figure 5. Bifurcation and phase portrait analyses of the proposed system using different values of the parameters together with miscellaneous states for  and  .

Figure 5 .
Figure 5. Bifurcation and phase portrait analyses of the proposed system using different values of the parameters together with miscellaneous states for W 1 and W 2 .

Figure 6 .
Figure 6.Chaotic patterns of the mentioned dynamical system along with many values of the parameters, including   2 ⁄ .

Figure 7 .Figure 7 .Figure 8 .
Figure 7. Chaotic patterns of the mentioned dynamical system together with diverse values of the parameters, including ζ 4 = π.

Figure 8 .
Figure 8. Chaotic patterns of the mentioned dynamical system with several values of the parameters, including ζ 4 = 2π.

Figure 9 .
Figure 9. Sketch of the sensitivity analysis of the proposed dynamical system made by choosi different values of the parameters with specific initial conditions: (a)  0 0.2 and  0, (  0 0 and  0.2, (c)  0 0.3 and  0, and (d)  0 0 and  0.3.

Figure 10 .
Figure 10.Bifurcation and phase portrait analyses of the proposed system using different values of the parameters together with miscellaneous states for  and  .

Figure 10 .
Figure 10.Bifurcation and phase portrait analyses of the proposed system using different values of the parameters together with miscellaneous states for W 3 and W 4 .

Figure 11 .
Figure 11.The (i) 3D, (ii) contour, and (iii) 2D shapes that have compaction solitons for Equation (25) with the given set of parameters.

Figure 12 .
Figure 12.The (i) 3D, (ii) contour, and (iii) 2D shapes that have singular periodic solitons for Equation (29) with the given set of parameters.

Figure 11 .
Figure 11.The (i) 3D, (ii) contour, and (iii) 2D shapes that have compaction solitons for Equation (25) with the given set of parameters.

Figure 11 .
Figure 11.The (i) 3D, (ii) contour, and (iii) 2D shapes that have compaction solitons for Equation (25) with the given set of parameters.

Figure 12 .
Figure 12.The (i) 3D, (ii) contour, and (iii) 2D shapes that have singular periodic solitons for Equation (29) with the given set of parameters.

Figure 12 .
Figure 12.The (i) 3D, (ii) contour, and (iii) 2D shapes that have singular periodic solitons for Equation (29) with the given set of parameters.Mathematics 2024, 12, x FOR PEER REVIEW 19 of 23

Figure 13 .
Figure 13.The (i) 3D, (ii) contour, and (iii) 2D shapes that have bell-shaped solitons for Equation (33) with the given set of parameters.

Figure 13 .
Figure 13.The (i) 3D, (ii) contour, and (iii) 2D shapes that have bell-shaped solitons for Equation (33) with the given set of parameters.

Figure 16 .
Figure 16.The (i) 3D, (ii) contour, and (iii) 2D shapes that have flat kink-shaped solitons for Equation (53) with the given set of parameters.

Figure 16 .
Figure 16.The (i) 3D, (ii) contour, and (iii) 2D shapes that have flat kink-shaped solitons for Equation (53) with the given set of parameters.

Figure 16 .
Figure 16.The (i) 3D, (ii) contour, and (iii) 2D shapes that have flat kink-shaped solitons for Equation (53) with the given set of parameters.

Figure 17 .
Figure 17.The (i) 3D, (ii) contour, and (iii) 2D shapes that have compaction solitons for Equation (56) with the given set of parameters.

Figure 17 .
Figure 17.The (i) 3D, (ii) contour, and (iii) 2D shapes that have compaction solitons for Equation (56) with the given set of parameters.