Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 6629095 | https://doi.org/10.1155/2021/6629095

Mingyuan Li, Wei Zhang, Qiliang Wu, "Analytical and Numerical Results on Global Dynamics of the Generalized Boussinesq Equation with Cubic Nonlinearity and External Excitation", Mathematical Problems in Engineering, vol. 2021, Article ID 6629095, 11 pages, 2021. https://doi.org/10.1155/2021/6629095

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

Academic Editor: Diego Oliva
Received13 Nov 2020
Revised05 Mar 2021
Accepted22 May 2021
Published01 Jun 2021

Abstract

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.

1. 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 [46] 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) [1622]. 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 [2326] 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 [2729]. 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 aswith the periodic boundary conditionswhere , , 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:where and .

Thus, 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:::::

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 . Letting , equation (9) is transformed intowhere , , 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:

Equation (11) contains the invariant subspacewhich owns the resonant torus:

When , we havewith Hamiltonian

Assuming , it is observed that equation (14) is completely integrable and has a periodic solution [2326]as well as an analytical solution [2326]where

The symbol represents different parts of solutions being homoclinic to . Meanwhile, can also be seen as a heteroclinic orbit connecting different points in periodic orbit due to the phase shifts, as presented in Figure 1. Figure 2 depicts the heteroclinic structure in space , respectively. To analyze the dynamic characteristics of equation (11) near , the following transformations are introduced:and substituting equation (21) into equation (11), we havewith a Hamilton energy

3.1. First Measurement

Note that the first measurement, the Melnikov method [2225], shown in Figure 3(a), can be expressed as

To ensure the transversality of SM and USM holds, we have

Thus, the persistence of the heteroclinic structure in is confirmed.

3.2. Second Measurement

According to [2730], the second measurement, called as the geometric analysis (Figure 3(b)), is described as

Inserting equation (23) into equation (24) and using , the chaotic threshold can be derived aswhere . The chaotic threshold curves, i.e., black, red, and blue curves, denote the threshold curve (25) as , 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.

4. Numerical Simulation

To test the aforementioned chaotic threshold (25), the differential quadrature method [2729] 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 and , Figure 5 displays the bifurcation diagram, where the abscissa and ordinate represent the damping coefficient 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 chaotic wave in space is presented in Figure 6 as damping coefficient . Figure 7 illustrates the chaotic motion for the GBE when damping coefficient . Figure 7(a) shows the phase trajectory on the plane , Figure 7(b) presents the time response on the plane , Figure 7(c) is the phase trajectory on the plane , Figure 7(d) denotes the time response on the plane , and Figures 7(e) and 7(f) are the phase trajectories in 3D spaces and , respectively.

The chaotic behavior in 3D space as is presented in Figure 8. Figure 9 displays the chaotic phenomenon when the damping coefficient . 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 , 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].

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.

References

  1. V. Roeber and J. D. Bricker, “Destructive tsunami-like wave generated by surf beat over a coral reef during Typhoon Haiyan,” Nature Communications, vol. 6, p. 7854, 2015. View at: Publisher Site | Google Scholar
  2. T. Shimozono, Y. Tajima, A. B. Kennedy, H. Nobuoka, J. Sasaki, and S. Sato, “Combined infragravity wave and sea‐swell runup over fringing reefs by super typhoon Haiyan,” Journal of Geophysical Research: Oceans, vol. 120, no. 6, pp. 4463–4486, 2015. View at: Publisher Site | Google Scholar
  3. V. Roeber and K. F. Cheung, “Boussinesq-type model for energetic breaking waves in fringing reef environments,” Coastal Engineering, vol. 70, pp. 1–20, 2012. View at: Publisher Site | Google Scholar
  4. J. Weiss, M. Tabor, and G. Carnevale, “The Painlevé property for partial differential equations,” Journal of Mathematical Physics, vol. 24, no. 3, pp. 522–526, 1983. View at: Publisher Site | Google Scholar
  5. J. Weiss, “The Painlevé property for partial differential equations. II: Bäcklund transformation, Lax pairs, and the Schwarzian derivative,” Journal of Mathematical Physics, vol. 24, no. 6, pp. 1405–1413, 1983. View at: Publisher Site | Google Scholar
  6. J. Weiss, “The Painlevé property and Bäcklund transformations for the sequence of Boussinesq equations,” Journal of Mathematical Physics, vol. 26, no. 2, pp. 258–269, 1985. View at: Publisher Site | Google Scholar
  7. K. Mohammed, “Exact traveling wave solutions of the Boussinesq-Burgers equation,” Mathematical and Computer Modelling, vol. 49, pp. 66–671, 2009. View at: Google Scholar
  8. X. Yang, C. Zhao, and J. Cao, “Dynamics of the discrete coupled nonlinear Schrödinger-Boussinesq equations,” Applied Mathematics and Computation, vol. 219, no. 16, pp. 8508–8524, 2013. View at: Publisher Site | Google Scholar
  9. M. A. Helal, A. R. Seadawy, and M. H. Zekry, “Stability analysis of solitary wave solutions for the fourth-order nonlinear Boussinesq water wave equation,” Applied Mathematics and Computation, vol. 232, pp. 1094–1103, 2014. View at: Publisher Site | Google Scholar
  10. M. Awais and S. Ibrahim, “Nonlinear instability for the Boussinesq equations with diabatic forcing,” Nonlinear Analysis, vol. 166, pp. 1–18, 2018. View at: Publisher Site | Google Scholar
  11. Z. Yang, J. Tang, and Y. Shen, “Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model,” Applied Ocean Research, vol. 70, pp. 32–40, 2018. View at: Publisher Site | Google Scholar
  12. J. Cai, J. Chen, and B. Yang, “Efficient energy-preserving wavelet collocation schemes for the coupled nonlinear Schrödinger-Boussinesq system,” Applied Mathematics and Computation, vol. 357, pp. 1–11, 2019. View at: Publisher Site | Google Scholar
  13. S. Wang and X. Su, “The Cauchy problem for the dissipative Boussinesq equation,” Nonlinear Analysis: Real World Applications, vol. 45, pp. 116–141, 2019. View at: Publisher Site | Google Scholar
  14. K. Yildirim, “On the boundary control of a Boussinesq system,” Journal of Informatics and Mathematical Sciences, vol. 12, pp. 199–208, 2020. View at: Google Scholar
  15. K. Yildirim, “Active control of an improved Boussinesq system,” Mathematical Modelling of Natural Phenomena, vol. 15, p. 58, 2020. View at: Publisher Site | Google Scholar
  16. R. Tian, Q. Cao, and S. Yang, “The codimension-two bifurcation for the recent proposed SD oscillator,” Nonlinear Dynamics, vol. 59, no. 1-2, pp. 19–27, 2010. View at: Publisher Site | Google Scholar
  17. Q. L. Wu, W. Zhang, and E. H. Dowell, “Detecting multi-pulse chaotic dynamics of high-dimensional non-autonomous nonlinear system for circular mesh antenna,” International Journal of Non-linear Mechanics, vol. 102, pp. 25–40, 2018. View at: Publisher Site | Google Scholar
  18. Q. L. Wu and G. Y. Qi, “Viscoelastic string-beam coupled vibro-impact system: modeling and dynamic analysis,” European Journal of Mechanics A-Solid, vol. 82, Article ID 104012, 2020. View at: Publisher Site | Google Scholar
  19. R. L. Tian, Z. J. Zhao, and Y. Xu, “Variable scale-convex-peak method for weak signal detection,” Science China Technological Sciences, vol. 63, 2020. View at: Publisher Site | Google Scholar
  20. Q. Wu and G. Qi, “Quantum dynamics for Al-doped graphene composite sheet under hydrogen atom impact,” Applied Mathematical Modelling, vol. 90, pp. 1120–1129, 2021. View at: Publisher Site | Google Scholar
  21. Q. Wu, M. Yao, M. Li, D. Cao, and B. Bai, “Nonlinear coupling vibrations of graphene composite laminated sheets impacted by particles,” Applied Mathematical Modelling, vol. 93, pp. 75–88, 2021. View at: Publisher Site | Google Scholar
  22. Y. Niu, W. Zhang, and X. Y. Guo, “Free vibration of rotating pretwisted functionally graded composite cylindrical panel reinforced with graphene platelets,” European Journal of Mechanics- A Solids, vol. 77, Article ID 103798, 2019. View at: Publisher Site | Google Scholar
  23. Y. Li and D. W. McLaughlin, “Homoclinic orbits and chaos in discretized perturbed NLS systems: Part I. Homoclinic orbits,” Journal of Nonlinear Science, vol. 7, no. 3, pp. 211–269, 1997. View at: Publisher Site | Google Scholar
  24. Y. Li and S. Wiggins, “Homoclinic orbits and chaos in discretized perturbed NLS systems: Part II. Symbolic dynamics,” Journal of Nonlinear Science, vol. 7, no. 4, pp. 315–370, 1997. View at: Publisher Site | Google Scholar
  25. C. Zeng, “Homoclinic orbits for a perturbed nonlinear Schrödinger equation,” Communications on Pure and Applied Mathematics, vol. 53, no. 10, pp. 1222–1283, 2000. View at: Publisher Site | Google Scholar
  26. J. Shatah and C. Zeng, “Homoclinic orbits for the perturbed Sine-Gordon equation,” Communications on Pure and Applied Mathematics, vol. 53, no. 3, pp. 283–299, 2000. View at: Publisher Site | Google Scholar
  27. W. Zhang, Q. L. Wu, M. H. Yao, and E. H. Dowell, “Analysis on global and chaotic dynamics of nonlinear wave equations for truss core sandwich plate,” Nonlinear Dynamics, vol. 94, no. 1, pp. 21–37, 2018. View at: Publisher Site | Google Scholar
  28. Q. Wu and G. Qi, “Global dynamics of a pipe conveying pulsating fluid in primary parametrical resonance: analytical and numerical results from the nonlinear wave equation,” Physics Letters A, vol. 383, no. 14, pp. 1555–1562, 2019. View at: Publisher Site | Google Scholar
  29. Q. Wu and G. Qi, “Homoclinic bifurcations and chaotic dynamics of non-planar waves in axially moving beam subjected to thermal load,” Applied Mathematical Modelling, vol. 83, pp. 674–682, 2020. View at: Publisher Site | Google Scholar
  30. M. Li, Q. Wu, and B. Bai, “Size-dependent mechanics of viscoelastic carbon nanotubes: modeling, theoretical and numerical analysis,” Results in Physics, vol. 19, Article ID 103383, 2020. View at: Publisher Site | Google Scholar
  31. J. Boussinesq, “Théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal,” Journal de mathématiques pures et appliquées, vol. 2, pp. 55–108, 1872. View at: Google Scholar

Copyright © 2021 Mingyuan Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views120
Downloads331
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.