This paper analytically and numerically presents global dynamics of the generalized Boussinesq equation (GBE) with cubic nonlinearity and harmonic excitation. The effect of the damping coefficient on the dynamical responses of the generalized Boussinesq equation is clearly revealed. Using the reductive perturbation method, an equivalent wave equation is then derived from the complex nonlinear equation of the GBE. The persistent homoclinic orbit for the perturbed equation is located through the first and second measurements, and the breaking of the homoclinic structure will generate chaos in a Smale horseshoe sense for the GBE. Numerical examples are used to test the validity of the theoretical prediction. Both theoretical prediction and numerical simulations demonstrate the homoclinic chaos for the GBE.
National Natural Science Foundation of China119022201129015211427801Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality1. Introduction
The field of global dynamics for infinite-dimensional systems is becoming increasingly attractive for interdisciplinary research. This 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–6] explored Lax pairs, Bäcklund transformation, and the Painlevé property of the BE. The 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–22]. Theoretically, 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. The 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]. The 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–26] to obtain the chaotic threshold of mechanical systems such as truss core sandwich plates, pipes conveying fluid, axial moving beams, and carbon nanotubes. The aforementioned literature focused on solid structure systems; however, analysis on global dynamics for fluid systems, e.g., GBE, remains lacking.
The 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. The reductive perturbation method (RPM) [27, 30] is employed to transform the complicated PDE of the GBE into a topologically equivalent wave equation. The 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–29]. The simulating results confirm the validity of the theoretical prediction and illustrate larger damping coefficient resulting in the more complicated dynamical behaviors of the GBE.
2. 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(1)qtt+μqt−qxx+q3+qxxxx=αcosωt,with the periodic boundary conditions(2)qt=qt+2π,qxt=qxt+2π,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:(3)qx,t,ε=εuX,Teiϕ+c.c.,ϕ=kx−ωt,where X=εx−ct and T=ε2t.
Thus, the derivatives used in the RPM are given as(4a)∂∂t⟶∂∂t−ε1/2c∂∂X+ε∂∂T,(4b)∂2∂t2⟶∂2∂t2−2ε1/2c∂2∂X∂t+εc2∂2∂X2+2ε∂2∂t∂T−2ε3/2c∂2∂X∂T+ε2∂2∂T2,(4c)∂∂x⟶∂∂x+ε1/2∂∂X,(4d)∂2∂x2⟶∂2∂x2+2ε1/2∂2∂x∂X+ε∂2∂X2,(4e)∂4∂x4⟶∂4∂x4+4ε1/2∂4∂x3∂X+6ε∂4∂x2∂X2+4ε3/2∂4∂x∂X3+ε2∂4∂X4.
Moreover, the following scale transformation is introduced:(5)μ⟶ε2μ.
By substituting equations (3)–(5) into equation (1) and removing the secular terms in the power of ε, the results are yielded as follows:
Oε:(6)ω2−k2−k4u=0.
Oε2:(7)cω−k+4k3∂u∂X=0.
Oε3:(8)−2iω∂u∂T+c2+1−k2∂2u∂X2−3k2u2u−iμωu=0.
Oε4:(9)2c∂2u∂X∂T+4ik∂3u∂X3−6ik∂u2u∂X−μc∂u∂X=0.
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 is −2icαe−iβω2T. Letting u⟶ue−iβω2T, equation (9) is transformed into(10)i∂u∂T=∂2u∂X2+βu2−ω2u+iμ2u−iα,where X∗=2k/cX, β=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). Therefore, by using equation (10), the nonlinear characteristics of the wave equation for the GBE can be explored.
3. Global Analysis
The dynamical characteristics of equation (10) are further studied in this section. A bookkeeping parameter ε is inserted and obtains the following equation:(11)i∂u∂T=∂2u∂X2+βu2−ω2u+iεμ2u−α.
Equation (11) contains the invariant subspace(12)∏=uX,T|∂Xu≡0,which owns the resonant torus:(13)S=u∈∏u≡ω.
When ε=0, we have(14)i∂u∂T=∂2u∂X2+βu2−ω2u,with Hamiltonian(15)Hu=∫01∂Xu2+βω2u2−β2u4dX.
Assuming β>0, it is observed that equation (14) is completely integrable and has a periodic solution [23–26](16)u0T=rexp−iβr2−ω2T+iT0,as well as an analytical solution [23–26](17)uh±X,T=cos2pcoshT¯−isin2psinhT¯±sinpcosXcoshT¯∓sinpcosXu0T,where(18)T¯=ϕT+T1,ϕ=2βr2−ω2,p=arctanσ.
The symbol ± represents different parts of solutions uh±X,T being homoclinic to u0T. Meanwhile, uh±X,T can also be seen as a heteroclinic orbit connecting different points in periodic orbit u0T due to the phase shifts, as presented in Figure 1. Figure 2 depicts the heteroclinic structure uh±X,T in space X,T,u, respectively. To analyze the dynamic characteristics of equation (11) near S, the following transformations are introduced:(19)u=ω2+εIexpiθ,τ=εT,and substituting equation (21) into equation (11), we have(20a)∂I∂τ=μ2ω2+εI−αω2+εIcosθ,(20b)∂θ∂τ=−βI+εαsinθω2+εI,with a Hamilton energy(21)EI,θ=β2I2+μ2ω2θ−ωαsinθ.
Schematic diagram of the homoclinic tubes.
The exact form of the heteroclinic orbit for equation (17).
3.1. First Measurement
Note that the first measurement, the Melnikov method [22–25], shown in Figure 3(a), can be expressed as(22)Δ=Re∫−∞+∞∫01iμ2∂ph∂Tph−α∂ph∂TdXdT=2ω2β−1μ2ϕω31−23ϕ2+αsin2ϕcos2θ.
(a) The first measurement. (b) The second measurement.
To ensure the transversality of SM WεsΠ and USM WuQε holds, we have(23)μ2ϕω31−23ϕ2+αsin2ϕcos2θ=0.
Thus, the persistence of the heteroclinic structure in WuQε∪WεsΠ is confirmed.
3.2. Second Measurement
According to [27–30], the second measurement, called as the geometric analysis (Figure 3(b)), is described as(24)μ2ϕω3+αsin2ϕcos2θ=0.
Inserting equation (23) into equation (24) and using cos2θ≤1, the chaotic threshold can be derived as(25)αμ≤12ω23sin2arctan2βω2−ω2ωarctan2βω2−ω23−2arctan22βω2−ω2,where β>1/2. The chaotic threshold curves, i.e., black, red, and blue curves, denote the threshold curve (25) as β=1.0, 2.0, and 3.0 in Figure 4, respectively. The chaotic motions of the GBE may be generated as the ratio α/μ locates below the curve.
The chaotic threshold curve is plotted with the parameter ratio α/μ via the parameter ω.
4. Numerical Simulation
To test the aforementioned chaotic threshold (25), the differential quadrature method [27–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. The bifurcation behaviors show that the increase of damping coefficient can enhance the complicated dynamics of the GBE in comparison to the results in [27].
The bifurcation diagram of the composite laminated plate is plotted as q via the damping coefficient μ.
The 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 q7,q˙7, Figure 7(b) presents the time response on the plane t,q7, Figure 7(c) is the phase trajectory on the plane q13,q˙13, Figure 7(d) denotes the time response on the plane t,q13, and Figures 7(e) and 7(f) are the phase trajectories in 3D spaces q7,q˙7,q13 and q13,q˙13,q7, respectively.
The chaotic wave in space x,t,q when damping coefficient μ=0.3.
The chaotic motion when damping coefficient μ=0.3. (a) The phase trajectory on plane q7,q˙7. (b) The phase trajectory on plane q13,q˙13. (c) The time response on plane t,q7. (d) The time response on plane t,q13. (e) The phase trajectory in 3D space q7,q˙7,q13. (f) The phase trajectory in 3D space q13,q˙13,q7.
The 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. The 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].
The chaotic wave in space x,t,q when damping coefficient μ=0.63.
The chaotic motion when damping coefficient μ=0.63.
The chaotic wave in space x,t,q when damping coefficient μ=0.71.
The chaotic motion when damping coefficient μ=0.71.
The chaotic wave in space x,t,q when damping coefficient μ=0.80.
The chaotic motion when damping coefficient μ=0.80.
5. Conclusions
The 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.
The 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 μ. Therefore, nonlinear wave responses for the GBE can be controlled by adjusting the damping coefficient.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank the National Natural Science Foundation of China (11902220, 11290152, and 11427801) and the Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality for supporting this study.
RoeberV.BrickerJ. D.Destructive tsunami-like wave generated by surf beat over a coral reef during Typhoon Haiyan20156785410.1038/ncomms88542-s2.0-84938765496ShimozonoT.TajimaY.KennedyA. B.NobuokaH.SasakiJ.SatoS.Combined infragravity wave and sea‐swell runup over fringing reefs by super typhoon Haiyan201512064463448610.1002/2015jc0107602-s2.0-84973560226RoeberV.CheungK. F.Boussinesq-type model for energetic breaking waves in fringing reef environments20127012010.1016/j.coastaleng.2012.06.0012-s2.0-84864075104WeissJ.TaborM.CarnevaleG.The Painlevé property for partial differential equations198324352252610.1063/1.525721WeissJ.The Painlevé property for partial differential equations. II: Bäcklund transformation, Lax pairs, and the Schwarzian derivative19832461405141310.1063/1.525875WeissJ.The Painlevé property and Bäcklund transformations for the sequence of Boussinesq equations198526225826910.1063/1.5266552-s2.0-0009231934MohammedK.Exact traveling wave solutions of the Boussinesq-Burgers equation20094966671YangX.ZhaoC.CaoJ.Dynamics of the discrete coupled nonlinear Schrödinger-Boussinesq equations2013219168508852410.1016/j.amc.2013.01.0532-s2.0-84876115469HelalM. A.SeadawyA. R.ZekryM. H.Stability analysis of solitary wave solutions for the fourth-order nonlinear Boussinesq water wave equation20142321094110310.1016/j.amc.2014.01.0662-s2.0-84894451304AwaisM.IbrahimS.Nonlinear instability for the Boussinesq equations with diabatic forcing201816611810.1016/j.na.2017.08.0012-s2.0-85032796730YangZ.TangJ.ShenY.Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model201870324010.1016/j.apor.2017.09.0012-s2.0-85037359731CaiJ.ChenJ.YangB.Efficient energy-preserving wavelet collocation schemes for the coupled nonlinear Schrödinger-Boussinesq system201935711110.1016/j.amc.2019.03.0582-s2.0-85063758468WangS.SuX.The Cauchy problem for the dissipative Boussinesq equation20194511614110.1016/j.nonrwa.2018.06.0122-s2.0-85049648847YildirimK.On the boundary control of a Boussinesq system202012199208YildirimK.Active control of an improved Boussinesq system2020155810.1051/mmnp/2020024TianR.CaoQ.YangS.The codimension-two bifurcation for the recent proposed SD oscillator2010591-2192710.1007/s11071-009-9517-92-s2.0-77949424206WuQ. L.ZhangW.DowellE. H.Detecting multi-pulse chaotic dynamics of high-dimensional non-autonomous nonlinear system for circular mesh antenna2018102254010.1016/j.ijnonlinmec.2018.03.0062-s2.0-85044468039WuQ. L.QiG. Y.Viscoelastic string-beam coupled vibro-impact system: modeling and dynamic analysis20208210401210.1016/j.euromechsol.2020.104012TianR. L.ZhaoZ. J.XuY.Variable scale-convex-peak method for weak signal detection20206310.1007/s11431-019-1530-4WuQ.QiG.Quantum dynamics for Al-doped graphene composite sheet under hydrogen atom impact2021901120112910.1016/j.apm.2020.10.025WuQ.YaoM.LiM.CaoD.BaiB.Nonlinear coupling vibrations of graphene composite laminated sheets impacted by particles202193758810.1016/j.apm.2020.12.008NiuY.ZhangW.GuoX. Y.Free vibration of rotating pretwisted functionally graded composite cylindrical panel reinforced with graphene platelets20197710379810.1016/j.euromechsol.2019.1037982-s2.0-85067612900LiY.McLaughlinD. W.Homoclinic orbits and chaos in discretized perturbed NLS systems: Part I. Homoclinic orbits19977321126910.1007/bf026780882-s2.0-0000038726LiY.WigginsS.Homoclinic orbits and chaos in discretized perturbed NLS systems: Part II. Symbolic dynamics19977431537010.1007/bf026781412-s2.0-0041716979ZengC.Homoclinic orbits for a perturbed nonlinear Schrödinger equation200053101222128310.1002/1097-0312(200010)53:10<1222::aid-cpa2>3.0.co;2-fShatahJ.ZengC.Homoclinic orbits for the perturbed Sine-Gordon equation200053328329910.1002/(sici)1097-0312(200003)53:3<283::aid-cpa1>3.0.co;2-2ZhangW.WuQ. L.YaoM. H.DowellE. H.Analysis on global and chaotic dynamics of nonlinear wave equations for truss core sandwich plate2018941213710.1007/s11071-018-4343-62-s2.0-85053534268WuQ.QiG.Global dynamics of a pipe conveying pulsating fluid in primary parametrical resonance: analytical and numerical results from the nonlinear wave equation2019383141555156210.1016/j.physleta.2019.02.0192-s2.0-85061915272WuQ.QiG.Homoclinic bifurcations and chaotic dynamics of non-planar waves in axially moving beam subjected to thermal load20208367468210.1016/j.apm.2020.03.013LiM.WuQ.BaiB.Size-dependent mechanics of viscoelastic carbon nanotubes: modeling, theoretical and numerical analysis20201910338310.1016/j.rinp.2020.103383BoussinesqJ.Théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal1872255108