Analytical and Numerical Results on Global Dynamics of the Generalized Boussinesq Equation with Cubic Nonlinearity and External Excitation

Beijing Key Laboratory of Nonlinear Vibrations and Strength of Mechanical Structures, Beijing University of Technology, Beijing 100124, China College of Mechanical Engineering, Beijing University of Technology, Beijing 100124, China Big Data Applied Research and Collaborative Innovation Center, Inner Mongolia University of Finance and Economics, Hohhot 010070, Inner Mongolia, China School of Statistics and Mathematics, Inner Mongolia University of Finance and Economics, Hohhot 010070, Inner Mongolia, China School of Artificial Intelligence, Tiangong University, Tianjin 300387, China Tianjin Key Laboratory of Advanced Technology of Electrical Engineering and Energy, Tiangong University, Tianjin 300387, China


Introduction
e field of global dynamics for infinite-dimensional systems is becoming increasingly attractive for interdisciplinary research. is paper aims to develop a theoretical analysis method of global dynamics for infinite-dimensional dynamical systems. An analytical method is conducted to detect the chaotic waves and nonlinear dynamical behaviors of the generalized Boussinesq equation (GBE) with cubic nonlinearity and external excitation. As a typical model in physics and mechanics, the Boussinesq equation (BE) possesses both dispersive and nonlinear attributes [1]. BE is a powerful tool to predict coastal hydrodynamics such as wave propagation, wave-current interaction, wave breaking, and nearshore circulation [2], even in extreme weather conditions such as typhoons and windstorms [3]. Many scholars have theoretically analyzed the characteristics of the BE. Weiss and his coworkers [4][5][6] explored Lax pairs, Bäcklund transformation, and the Painlevé property of the BE. e homogeneous balance method was extended to study the solution of traveling waves for the Boussinesq-Burgers equation [7]. Yang et al. [8] analyzed the global attractors and asymptotic behaviors for Schrödinger-Boussinesq equations. Based on the Hamiltonian structure, second-and fourth-order energy-preserving wavelet collocation schemes were proposed by Helal et al. [9] for the nonlinear coupled Schrödinger-Boussinesq model. Awais and Ibrahim [10] carried out studies on water waves in a fourth-order Boussinesq nonlinear model and analyzed its stability for soliton solutions. Considering the impact of diabatic force, the instability of Boussinesq equations was investigated [11]. In [12], a numerical analysis was conducted to investigate the sensitivities of coastal wave attenuation resulting from incident wave height, period, water depth, and vegetation configurations in a nonlinear model of Boussinesq equations. Wang and Su [13] studied the finite time blow-up and solvability of solutions in all space dimensions for the Cauchy problem of a Boussinesq equation with dissipation. Yildirim analyzed, theoretically and numerically, the boundary control problem [14] and active control [15] for the nonlinear Boussinesq system.
Traditionally, infinite-dimensional systems are described as PDEs (partial differential equations), and then the PDEs are truncated into ODEs (ordinary differential equations) [16][17][18][19][20][21][22]. eoretically, ODEs are topologically equivalent to PDEs, but the dynamical phenomena cannot be perfectly exhibited. Hence, how to improve an analytical method of global dynamics for infinite-dimensional systems must be resolved still by scientists and scholars. In recent decades, a level of improvement in this area has been achieved. e first studies of the global perturbation technique for nonlinear PDEs were conducted by Li et al. [23,24], Zeng [25], and Shatah and Zeng [26]. e persistent homoclinic orbits of nonlinear Schrödinger (NLS) and sine-Gordon (SG) equations were constructed by using invariant manifolds, foliations, and Melnikov analysis. Furthermore, Zhang et al. [27], Wu and Qi [28,29], and Li et al. [30] improved the analytic theory for infinite-dimensional systems [23][24][25][26] to obtain the chaotic threshold of mechanical systems such as truss core sandwich plates, pipes conveying fluid, axial moving beams, and carbon nanotubes. e aforementioned literature focused on solid structure systems; however, analysis on global dynamics for fluid systems, e.g., GBE, remains lacking. e objective of this study is to derive the analytical chaotic threshold of the GBE and indicate how the damping coefficient can clearly affect the dynamical responses of the GBE. e reductive perturbation method (RPM) [27,30] is employed to transform the complicated PDE of the GBE into a topologically equivalent wave equation. e first and second measurements are then applied to explore the homoclinic chaos and the chaotic threshold for the wave equation. Afterward, numerical examples are also carried out using the differential quadrature method [27][28][29]. e simulating results confirm the validity of the theoretical prediction and illustrate larger damping coefficient resulting in the more complicated dynamical behaviors of the GBE.

Motion Equation and Perturbation Analysis
Boussinesq [31] assumed that the horizontal velocity is constant along the direction of water depth, and the vertical velocity is linearly distributed, obtaining the one-dimensional nonlinear water wave governing equation. Based on literature [31], the nonlinear governing equation of the cubic nonlinearity and harmonic excited GBE can be proposed as with the periodic boundary conditions where μ, α, and ω stand for the damping coefficient, amplitude, and frequency of the external force, respectively.
To carry out the analysis on the complicated dynamical responses of equation (1), the RPM is applied to simplify the formula, and the following formal solution of equation (1) yields: us, the derivatives used in the RPM are given as Moreover, the following scale transformation is introduced: By substituting equations (3)-(5) into equation (1) and removing the secular terms in the power of ε, the results are yielded as follows: 2c In order to derive the perturbed equation possessing the same form in [23,24], equation (9) is integrated with respect to X, and the integral constant where X * � ����� � (2k/c)X, β � 3k/c, and the asterisk is dropped for convenience.
Hence, wave equation (10) has the same form in [23,24] and is thus topologically equivalent to equation (1). erefore, by using equation (10), the nonlinear characteristics of the wave equation for the GBE can be explored.

Global Analysis
e dynamical characteristics of equation (10) are further studied in this section. A bookkeeping parameter ε is inserted and obtains the following equation: Equation (11) contains the invariant subspace which owns the resonant torus: When ε � 0, we have with Hamiltonian Assuming β > 0, it is observed that equation (14) is completely integrable and has a periodic solution [23][24][25][26] as well as an analytical solution [23][24][25][26] where e symbol ± represents different parts of solutions u ± h (X, T) being homoclinic to u 0 (T). Meanwhile, u ± h (X, T) can also be seen as a heteroclinic orbit connecting different points in periodic orbit u 0 (T) due to the phase shifts, as presented in Figure 1. Figure 2 depicts the heteroclinic structure u ± h (X, T) in space (X, T, |u|), respectively. To analyze the dynamic characteristics of equation (11) near S, the following transformations are introduced: and substituting equation (21) into equation (11), we have with a Hamilton energy Mathematical Problems in Engineering

First Measurement.
Note that the first measurement, the Melnikov method [22][23][24][25], shown in Figure 3(a), can be expressed as To ensure the transversality of SM W s ε (Π) and USM W u (Q ε ) holds, we have us, the persistence of the heteroclinic structure in W u (Q ε ) ∪ W s ε (Π) is confirmed.

Numerical Simulation
To test the aforementioned chaotic threshold (25), the differential quadrature method [27][28][29] in MATLAB is employed using equation (1) to study the influence of the damping coefficient μ on the nonlinear dynamical behaviors of the GBE with cubic nonlinearity. Letting α � 1.5 and ω � 1.0, Figure 5 displays the bifurcation diagram, where the abscissa and ordinate represent the damping coefficient μ ∈ [0, 1] and displacement of the GBE, respectively. e bifurcation behaviors show that the increase of damping coefficient can enhance the complicated dynamics of the GBE in comparison to the results in [27]. e chaotic wave in space (x, t, q) is presented in Figure 6 as damping coefficient μ � 0.3. Figure 7 illustrates the chaotic motion for the GBE when damping coefficient μ � 0.3. Figure 7(a) shows the phase trajectory on the plane (q 7 , _ q 7 ), Figure 7 13 ) and (q 13 , _ q 13 , q 7 ), respectively. e chaotic behavior in 3D space (x, t, q) as μ � 0.63 is presented in Figure 8. Figure 9 displays the chaotic phenomenon when the damping coefficient μ � 0.63. When μ is increased to 0.71, Figures 10 and 11 illustrate the chaotic motion and chaotic wave for the GBE, respectively. When the damping coefficient μ � 0.80, the chaotic motion and chaotic wave are depicted, respectively, in Figures 12 and 13. e observed results from Figures 6 to 13 show that the amplitude of the GBE is larger for smaller μ. However, the higher peak value of the truss core sandwich plate system is exhibited when the damping coefficient is sufficiently small or large [27].      Mathematical Problems in Engineering 9

Conclusions
e global perturbation method for the infinite-dimensional system was used to discuss global dynamics for the GBE in this work. Using the RPM, an analytical wave equation was obtained by simplifying the complicated nonlinear equation of the GBE. Afterward, the persistent homoclinic structure for wave equations was constructed by using the first and second measurements. Finally, quantitative analysis was conducted to examine the theoretical results using the differential quadrature method. e numerical studies illustrate that the damping coefficient μ is essential to determine the dynamical responses of the GBE. Different bifurcation phenomena are shown with changing damping coefficient μ. erefore, nonlinear wave responses for the GBE can be controlled by adjusting the damping coefficient.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.