Some Novel Analytical Approximations to the (Un)damped Duffing–Mathieu Oscillators

Alyousef, Haifa A.; Salas, Alvaro H.; Alkhateeb, Sadah A.; El-Tantawy, S. A.

doi:https://doi.org/10.1155/2022/2715767

Journal of Mathematics

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Methods in Dynamical Systems and Applications in Engineering

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2715767 | https://doi.org/10.1155/2022/2715767

Some Novel Analytical Approximations to the (Un)damped Duffing–Mathieu Oscillators

Haifa A. Alyousef,¹Alvaro H. Salas,²Sadah A. Alkhateeb,³and S. A. El-Tantawy^4,5

Academic Editor: Fairouz Tchier

Received24 Dec 2021

Accepted27 Mar 2022

Published05 May 2022

Abstract

Some novel exact solutions and approximations to the damped Duffing–Mathieu-type oscillator with cubic nonlinearity are obtained. This work is divided into two parts: in the first part, some exact solutions to both damped and undamped Mathieu oscillators are obtained. These solutions are expressed in terms of the Mathieu functions of the first kind. In the second part, the equation of motion to the damped Duffing–Mathieu equation (dDME) is solved using some effective and highly accurate approaches. In the first approach, the nonintegrable dDME with cubic nonlinearity is reduced to the integrable dDME with linear term having undermined optimal parameter (maybe called reduced method). Using a suitable technique, we can determine the value of the optimal parameter and then an analytical approximation is obtained in terms of the Mathieu functions. In the second approach, the ansatz method is employed for deriving an analytical approximation in terms of trigonometric functions. In the third approach, the homotopy perturbation technique with the extended Krylov–Bogoliubov–Mitropolskii (HKBM) method is applied to find an analytical approximation to the dDME. Furthermore, the dDME is solved numerically using the Runge–Kutta (RK) numerical method. The comparison between the analytical and numerical approximations is carried out. All obtained approximations can help a large number of researchers interested in studying the nonlinear oscillations and waves in plasma physics and many other fields because many evolution equations related to the nonlinear waves and oscillations in a plasma can be reduced to the family of Mathieu-type equation, Duffing-type equation, etc.

1. Introduction

Several physical and natural phenomena related to biology, chemistry, physics, engineering problems, and so on can be modelled by both ordinary differential equations (ODEs) and partial differential equations (PDEs) for studying the nonlinear self-excited oscillators [1–12]. Also, in many real-life problems, some internal and external forces that can affect the system under consideration cannot be neglected. For example, the friction and collisional force and many others that affect on the motion of particles, whether in solid, liquid, gas or plasma physics, cannot be neglected. Therefore, these forces must be included in the mathematical models that will be used for studying the natural and physical problem, such as investigating the nonlinear oscillations in various plasma models [13–23]. Interest in the study of nonlinear oscillations in a plasma is due to its many potential applications. Nowadays, plasma processing is seen as an important and effective technology which has been able to enter into all modern industries. In addition, plasma had a great credit for the modern technology in electronics, medicine, agriculture, biomedicine, automobiles, optics, aerospace, telecommunications, solar energy, polymers, papers, textiles, etc. [16, 24]. Therefore, we focus our attention on the applications of the family of Duffing-type equation for modelling the nonlinear oscillations in a plasma. The Duffing-type oscillator is one of the most popular differential equations that has spread widely due to its ability to explain many nonlinear phenomena in various fields of science and in mechanical systems and engineering. This equation is a mathematical model described by a second order of ODEs with a nonlinear spring force. It is used for describing the motion of (un)damped oscillator rather than simple harmonic motion. Motivated by potential engineering and plasma physics applications in addition to many others applications in electronic circuits and microcontroller, the family of Duffing equation such as (un)damped Duffing equation with lower and higher-order nonlinearities, (un)damped Duffing–Helmholtz equation, and (un)damped Mathieu–Duffing equation have received wide attention due to their ability for investigating the mechanism of a rigid pendulum oscillator, oscillations in different plasma models, and so on. This family of second-order differential equations has provided useful and successful models for investigating the nonlinear oscillations and chaotic nature. The biggest challenge in the study of the dynamics of nonlinear mechanical systems is to find some real solutions (including the analytical and numerical solutions) to the evolution equations that are used for describing the characteristics of nonlinear phenomena under consideration. Accordingly, studying the solutions of many equations of motion to various oscillators is one of the most difficult tasks facing many researchers.

The solutions associated to mentioned evolution equations and many other related equations have been studied extensively due to the fact that such equations arise in a variety of realistic problems. For instance, the fluid equations of electron-ion unmagnetized cold plasma were reduced to Mathieu equation in order to investigate the electron waves [13]. During this study, the authors assumed that the density perturbations of the plasma species are only time-dependent functions and do not depend on space. Also, the basic set of fluid equations to a multicomponent complex plasma consisting of inertial two types of dust grains (both positive and negative charges) as well as inertialess Maxwellian species including electrons and ions were reduced to a Mathieu-type equation for studying the excitation of dust-acoustic oscillations [14]. Moreover, a nonlinear Van der Pol–Mathieu-type equation was derived for the dust grain density in order to investigate the dynamics of dust-acoustic oscillations in a dusty plasma consisting of inertial Boltzmann distributed species (electrons and ions) and inertialess dust grains [15]. A modified Van der Pol–Duffing oscillator with forced term was used in the study and was used in modelling the dynamics of nonlinear oscillations in different plasma models [16]. More recently, the multistage method was used for solving the damping Duffing equation with forced term in order to model the oscillations in a complex unmagnetized plasma [17].

Due to the importance of the family of Duffing-type equation and motivated by the mentioned studies, we focus our attention on the analysis of the so-called (un)damped Mathieu–Duffing oscillator with a twin-well potential [25].andfor studying the vibrating/oscillating behavior of systems described by (1) and (2), where are, respectively, the stiffness coefficients of the linear and nonlinear terms and represents the coefficients of damping term. The longitudinal loading is periodic, and are the frequency and excitation strength of the periodic loading, respectively. The total energy of the undamped Mathieu–Duffing oscillator (DMO) according to (1) is defined by , where and are, respectively, the unperturbed and perturbation in the Hamiltonian of (1). There is another form to the equation of periodic motion which is called (un)damped Mathieu–Helmholtz oscillator.

The analytical solution to the damped Mathieu–Helmholtz oscillator (3) was obtained using the finite Fourier series expansion [26]. Note that in evolution equation (3), the nonlinear term is different from the nonlinear term in equations (1) and (2).

The objectives of our study are to find some novel solutions to the (un)damped Duffing–Mathieu-type oscillator, under the initial conditions and . Two cases for Duffing–Mathieu-type oscillator will be discussed. In the first case, we will get some exact solutions to (un)damped Mathieu equation in terms of the Mathieu functions of the first kind. In the second case, the damped Duffing–Mathieu oscillator (dDMO) will be solved analytically and numerically using some different approaches. In the first approach, the cubic nonlinear term in equation (1) is replaced by the linear term , where the constant represents an optimal parameter. Then, the nonintegrable dDMO reduces to an integrable one which has an exact solution but with undermined parameter . Using a suitable technique, we can determine the value of the optimal parameter . Thus, we can get an analytical approximation to the dDMO (2) in terms of the Mathieu functions. For the second approach, the ansatz method with the help of the solution to the undamped Duffing oscillator is employed to derive an analytical approximation to the dDMO (2) in the form of trigonometric functions. Furthermore, the homotopy perturbation technique with the extended Krylov–Bogoliubov–Mitropolskii (KBM) which is called HKBM method is also devoted for solving the dDME (2) for arbitrary physical parameters [27, 28].

2. Mathematical Analysis

Here, we proceed to find some approximate solutions to both undamped and damped Duffing–Mathieu oscillators (1) and (2), respectively. Below we discuss the different approaches for solving the mentioned equations.

2.1. An Exact Solution to (Un)damped Mathieu Equation

Both undamped Mathieu equation, i.e., (1) for , and the damped Mathieu equation, i.e., (2) for , have exact solutions in the form of Mathieu functions. First, let us find the solution of the following damped Mathieu equation:with subjected to the initial conditions (ICs): and .

Using the following MATHEMATICA commandwe can get the exact solution to (4) as follows:withwhere MathieuS and MathieuC are the Mathieu functions of the first kind or sometimes called sine-elliptic and cosine-elliptic, respectively. For , the damped Mathieu equation (4) reduces to the following undamped Mathieu equation:and solution (6) reduces to the following one:

The exact solution (6) is compared with RK numerical solution as shown in Figures 1(a) and 1(b) for different values of .

(a)

(b)

2.2. Some Analytical Approximations to the Damped Duffing–Mathieu Oscillator

Here, we proceed to discuss two techniques (the hybrid expansion method and the ansatz method) for finding some analytical approximations to dDMO (2). For studying dDMO (2), first we rewrite equation (2) in the following new i.v.p.

2.2.1. First Approach: Reduced Method

For and , we may obtain simple approximation to the i.v.p. (10) by replacing the cubic term by the linear term , where the constant which is used as an optimal parameter to reduce the residual error. Accordingly, dDMO (2) of cubic nonlinearity reduces to the following dDMO with linear term. Thus, we can replace the i.v.p. (10) by the following new i.v.p.

Thus, the exact solution to the i.v.p. (11) is expressed bywith

The residual is defined as

A suitable value of can be obtained by solving the equation for some , say . Making use of the Padé approximate technique, for and , we getwith

For a given , the residual error of the approximation (12) to the i.v.p. (10) is defined as

The optimal value of for the parameter on the interval is defined as

Let us apply the obtained approximation (12) for investigating the properties of the damping oscillations to the dDMO (2) at different values to the physical parameters . The profile of the approximation (12) using the values of (given in equation (15)) and (given in (18)) is displayed in Figures 2 and 3 at and , respectively. The obtained results showed that this approximation gives results with good and acceptable accuracy.

(a)

(b)

(a)

(b)

2.2.2. The Ansatz Method for Solving dDMO

Now, we can summarize the main points to get some approximations to the i.v.p. (10) in the following steps.

Step 1. Let us assume the following ansatz:where the time-dependent function can be determined later and the function represents the solution of the following i.v.p.

Step 2. Inserting ansatz (19) into i.v.p. (10), we get

Step 3. For vanishing the coefficient of in equation (21), we get

Step 4. Integrating equation (22) and with the help of , we havewhere stands for the EllipticE function.
Now, the value of has been determined but the solution of the i.v.p. (20) needs to be determined. Thus, we are faced with two things: either we use one of the solutions found in the literature [29] or try to find another solution in the form of trigonometric functions.

Step 5. In this step, we proceed to find a solution to i.v.p. (20) in the form of trigonometric functions. Without loss of generality, i.v.p. (20) can be redefined aswhere i.v.p. (24) is the same as i.v.p. (20) with , , and .

Step 6. Our objective here is to derive another solution that does not involve elliptic functions but an elementary solution. To do that, we assume , and for , we getwhere and is undermined function.

Step 7. By substituting solution (25) into i.v.p. (24), we have

Step 8. For vanishing the coefficient of in equation (26): , we haveand by solving (27) with , we getwithwhere .

Step 9. The number is obtained from the condition and it is a solution to the quarticwhere the number is a free/optimal parameter that is chosen in order to minimize the residual error. Its default value is .

Step 10. Finally, the trigonometric approximation to i.v.p. (10) is obtained:

Step 11. Also, we can solve i.v.p. (20) using RK numerical method and then replacing (given in equation (23)). The following MATHEMATICA command is introduced for this purpose:Both analytical and numerical approximations (31) and (33) to i.v.p. (10) are, respectively, plotted against the RK numerical solution as illustrated in Figures 4 and 5. Also, at , the maximum global distance of both approximations (31) and (33) is estimated for different values to asThe obtained results show the high accuracy and efficiency of the obtained approximations (31) and (33). Moreover, these approximations are stable against long time and for arbitrary values of the physical parameters.

(a)

(b)

(a)

(b)

2.3. The Homotopy Extended Krylov–Bogoliubov–Mitropolskii Method

The homotopy extended Krylov–Bogoliubov–Mitropolskii (HKBM) method may be used for solving both conservative and nonconservative oscillators. Based on this method (more details can be found in [27, 28]), i.v.p. (10) can be redefined aswhere indicates the solution of i.v.p. (35) while the solution of the dDMO (2) is obtained for . For , and , i.v.p. (35) can be written in the following reduced form:

According to the HKBM method, the following ansatz solution is introduced:where each is a periodic function in , and both amplitude and phase are assumed to vary with time and subject to the conditionswhere and .

Inserting ansatz solution (37) and using (38) and after several tedious calculations, we can determine the unknown time-dependent functions . To avoid the so-called secularity, we choose only the solution that does not contain nor . For , we getand

By solving system (39), we have

The first-order approximate solution is obtained for :where the values of are defined in (41) while the constants and can be obtained from the initial conditions.

The comparison between the HKBM first-order approximate solution (37) and the RK numerical simulation is reported as shown in Figures 6(a) and 6(b) for and , respectively. Also, both HKBM first-order approximate solution (42) and RK numerical simulation are, respectively, compared with each other for weak and strong damping as illustrated in Figures 7(a) and 7(b). Furthermore, at and for different values to , the maximum global distance error to the HKBM first-order approximate solution (42) is estimated as

(a)

(b)

(a)

(b)

It is clear that the HKBM first-order approximate solution (42) is characterized by high accuracy and more stability at long time.

3. Conclusion

Given the importance of nonlinear oscillations in plasma physics and engineering and their strong connection to the family of the Duffing-type oscillator, in this work, some exact solutions to the damped and undamped Mathieu equations as well as some analytical approximations to the damped Duffing–Mathieu oscillator (dDMO) using different approaches have been obtained. The exact solutions to both damped and undamped Mathieu equation have been obtained in the terms of Mathieu functions of the first kind. These solutions are numerically compared with the Runge–Kutta (RK) numerical simulation. It was observed that both exact and numerical solutions are completely matched with each other in the whole time interval. On the other hand, the dDMO has been solved using some different approaches. In the first one, the nonintegrable dDMO with cubic nonlinear term has been reduced to an integrable dDMO with linear term in which is undermined optimal parameter. The kappa optimal parameter has been determined using a suitable technique as we discussed in the text above. After determining the kappa optimal parameter, a highly accurate analytical approximation has been obtained in terms of the Mathieu functions. In the second approach, a highly accurate analytical approximation has been derived in detail in terms of trigonometric functions using the ansatz method. In the third technique, the homotopy extended Krylov–Bogoliubov–Mitropolskii (HKBM) method was used for getting an effective analytical approximation. Furthermore, the dDMO has been analyzed numerically using the RK numerical method. The comparison between all obtained approximations and the RK numerical solutions has been carried out. Moreover, the maximum global distance error in the whole time interval to all obtained approximations has been estimated. All obtained approximations are characterized by the high accuracy and efficiency in addition to being more stable for a long time.

3.1. Future Work

We may solve the following oscillators using of the methods described in this paper:

3.1.1. Future Idea I

Cubic-quintic Duffing–Mathieu equation:

3.1.2. Future Idea II

Forced damped Duffing–Mathieu equation:

3.1.3. Future Idea III

The forced Van der Pol–Duffing oscillator:and many others oscillators.

Data Availability

The data generated or analyzed during this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed equally to this study and approved the final version of the manuscript.

Acknowledgments

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project (grant no. PNURSP2022R17), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

References

A.-M. Wazwaz, Partial Differential Equations and Solitary Waves Theory, Springer, Berlin, Germany, 2009.
A.-M. Wazwaz, Partial Differential Equations: Methods and Applications, A.A.Balkema Publishers, Avereest, Netherlands, 2002.
S. A. Almutlak and S. A. El-Tantawy, “On the approximate solutions of a damped nonplanar modified Korteweg-de Vries equation for studying dissipative cylindrical and spherical solitons in plasmas,” Results in Physics, vol. 23, Article ID 104034, 2021.
View at: Publisher Site | Google Scholar
W. Albalawi, A. H. Salas, S. A. El-Tantawy, and A. A. A.-R. Youssef, “Approximate analytical and numerical solutions to the damped pendulum oscillator: Newton-Raphson and moving boundary methods,” Journal of Taibah University for Science, vol. 15, no. 1, pp. 479–485, 2021.
View at: Publisher Site | Google Scholar
N. H. Aljahdaly, S. A. El-Tantawy, A.-M. Wazwaz, and H. A. Ashi, “Adomian decomposition method for modelling the dissipative higher-order rogue waves in a superthermal collisional plasma,” Journal of Taibah University for Science, vol. 15, no. 1, pp. 971–983, 2021.
View at: Publisher Site | Google Scholar
S. Rajasekar, S. Parthasarathy, and M. Lakshmanan, “Prediction of horseshoe chaos in BVP and DVP oscillators,” Chaos, Solitons & Fractals, vol. 2, no. 3, pp. 271–280, 1992.
View at: Publisher Site | Google Scholar
K. Johannessen, “An analytical solution to the equation of motion for the damped nonlinear pendulum,” European Journal of Physics, vol. 35, no. 3, Article ID 035014, 2014.
View at: Publisher Site | Google Scholar
E. Esmailzadeh and G. Nakhaie-jazar, “Periodic solution of a Mathieu-Duffing type equation,” International Journal of Non-linear Mechanics, vol. 32, no. 5, pp. 905–912, 1997.
View at: Publisher Site | Google Scholar
M. Zamora, “A note on the periodic solutions of a Mathieu–Duffing type equations,” Mathematische Nachrichten, vol. 290, 2016.
View at: Publisher Site | Google Scholar
V. Belovodskiy, “About the new periodic solutions of the mathieu-duffing equation in the principal zone of instability,” Mechanics, vol. 33, pp. 38–42, 2010.
View at: Google Scholar
E. Babilio, “The duffing-mathieu equation arising from dynamics of post-buckled beams,” Nonlinear Dynamics of Structures, Systems and Devices, Springer, Berlin, Germany, pp. 267–275, 2020.
View at: Publisher Site | Google Scholar
A. Y. T. Leung, Z. Guo, and H. X. Yang, “Transition curves and bifurcations of a class of fractional mathieu-type equations,” International Journal of Bifurcation and Chaos, vol. 22, no. 11, Article ID 1250275, 2012.
View at: Publisher Site | Google Scholar
G. Brodin and L. Stenflo, “A simple electron plasma wave,” Physics Letters A, vol. 381, no. 11, pp. 1033–1035, 2017.
View at: Publisher Site | Google Scholar
I. Kourakis, M. Momeni, and P. K. Shukla, “A Mathieu equation for dust charge dynamics in multi-component dusty plasmas,” in Proceedings of the 33rd EPS Conference, vol. 30I, Rome, Italy, June 2006.
View at: Google Scholar
M. Momeni, I. Kourakis, M. Moslehi-Fard, and P. K. Shukla, “A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas,” Journal of Physics A: Mathematical and Theoretical, vol. 40, no. 24, pp. F473–F481, 2007.
View at: Publisher Site | Google Scholar
C. H. Miwadinou, L. A. Hinviy, A. V. Monwanouz, and J. B. Chabi Orou, “Nonlinear dynamics of plasma oscillations modeled by a forced modified Van der Pol-Duffing oscillator,” 2013, http://arxiv.org/abs/1308.6132v1.
View at: Google Scholar
N. H. Aljahdaly and S. A. El-Tantawy, “On the multistage differential transformation method for analyzing damping duffing oscillator and its applications to plasma physics,” Mathematics, vol. 9, no. 4, p. 432, 2021.
View at: Publisher Site | Google Scholar
A. H. Salas, S. A. El-Tantawy, and E. Jairo, “On the approximate and analytical solutions to the fifth-order Duffing oscillator and its physical applications,” Waves in Random and Complex Media, pp. 1–21, 2021.
View at: Publisher Site | Google Scholar
A. H. Salas S, S. A. El-Tantawy, and M. R. Alharthi, “Novel solutions to the (un)damped Helmholtz-Duffing oscillator and its application to plasma physics: moving boundary method,” Physica Scripta, vol. 96, no. 10, Article ID 104003, 2021.
View at: Publisher Site | Google Scholar
S. A. El-Tantawy, A. H. Salas, and M. R. Alharthi, “On the analytical and numerical solutions of the damped nonplanar Shamel Korteweg-de Vries Burgers equation for modeling nonlinear structures in strongly coupled dusty plasmas: multistage homotopy perturbation method,” Physics of Fluids, vol. 33, no. 4, Article ID 043106, 2021.
View at: Publisher Site | Google Scholar
S. A. El-Tantawy, A. H. Salas, and M. R. Alharthi, “On the analytical solutions of the forced damping duffing equation in the form of weierstrass elliptic function and its applications,” Mathematical Problems in Engineering, vol. 2021, Article ID 6678102, 9 pages, 2021.
View at: Publisher Site | Google Scholar
S. A. El-Tantawy, A. H. Salas, and M. R. Alharthi, “A new approach for modelling the damped Helmholtz oscillator: applications to plasma physics and electronic circuits,” Communications in Theoretical Physics, vol. 73, no. 3, Article ID 035501, 2021.
View at: Publisher Site | Google Scholar
A. H. Salas, S. A. El-Tantawy, and N. H. Aljahdaly, “An exact solution to the quadratic damping strong nonlinearity duffing oscillator,” Mathematical Problems in Engineering, vol. 2021, Article ID 8875589, 8 pages, 2021.
View at: Publisher Site | Google Scholar
J. Proud, Plasma Processing of Materials: Scientific Opportunities and Technologies Challenges, National Academy Press, Washington, DC, USA, 1991.
A. C. J. Luo, “Chaotic motion in the generic separatrix band of a Mathieu-Duffing oscillator with a twin-well potential,” Journal of Sound and Vibration, vol. 248, no. 3, pp. 521–532, 2001.
View at: Publisher Site | Google Scholar
A. C. J. Luo and B. Yu, “Analytical periodic motions of a parametric oscillator with quadratic nonlinearity,” in Proceedings of the 26th Conference on Mechanical Vibration and Noise, Buffalo, NY, USA, August 2014.
View at: Publisher Site | Google Scholar
M. A. Uddin, M. W. Ullah, and R. S. Bipasha, “An approximate analytical technique for solving second order strongly nonlinear generalized duffing equation with small damping,” Journal of Bangladesh Academy of Sciences, vol. 39, no. 1, pp. 103–114, 2015.
View at: Publisher Site | Google Scholar
A. H. Salas, “Analytical approximant to a quadratically damped duffing oscillator,” The Scientific World Journal, vol. 2022, Article ID 3131253, 10 pages, 2022.
View at: Publisher Site | Google Scholar
A. H. Salas and S. A. El-Tantawy, “On the approximate solutions to a damped harmonic oscillator with higher-order nonlinearities and its application to plasma physics: semi-analytical solution and moving boundary method,” The European Physical Journal Plus, vol. 135, no. 10, p. 833, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Haifa A. Alyousef 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.

PDF Download Citation

Download other formats

Order printed copies