Van Der Pol Oscillator Research Articles

Influence of the magnetic flux on the dynamics of a self-sustaining system: analytical, numerical and analogical investigations

This paper investigates the nonlinear dynamics of a ferroelectric enzyme-substrate reaction modeled by the birhythmic van der Pol oscillator coupled to the magnetic flux. We derive the equilibrium points and study their stability. We analyze some bifurcation structures and the variation of the Lyapunov exponents. The phenomena of symmetric attractors and the anti-monotonicity are observed. By increasing the magnetic flux, we find that the equilibrium points are stable, tends to control chaotic regimes, and affects regular and quasi-regular ones. As the magnetic flux increases, the amplitude of the oscillations around the equilibrium points decreases and the amplitude of the limit cycles at the Hopf bifurcation tends to disappear. Further increasing the magnetic flux gives rise to chaotic dynamics. The electrical circuit and analogical simulations are derived using the PSpice software. The agreement between analogical and numerical results is acceptable.

Characterization of cracking-healing cycle of asphalt mixture based on air-coupled second harmonic measurement using Duffing-van der pol oscillator

Second harmonic generation (SHG) technique using Rayleigh surface wave has been successfully attempted for the non-destruction evaluation of asphalt mixture during the self-healing process. The air-coupled transducer placement is desirable for the convenient implementation of SHG for the large-scale detection. Nevertheless, the pronounced attenuation of ultrasonic wave by air-coupled propagation will cause the second harmonic induced by microcracks to be barely detectable, particularly under the influence of inevitable noise. In this study, a novel approach by combining Duffing oscillator and van der pol equation is developed to identify and quantify the second harmonic in noise-contained signals based on the phase trajectory and periodic trajectory area exponent. The multi-physics finite element model is first established to simulate the air-asphalt propagation of Rayleigh waves, and proof-of-concept numerical simulations are carried out to estimate the second harmonic amplitude of noise-contaminated signals. The air-coupled SHG experiments are conducted, and the approach of the Duffing-van der pol oscillator is applied to extract the second harmonic in experimental data, and the self-healing behavior of asphalt mixture is analyzed through the variation of the nonlinear parameter of SHG measurements. The results indicate that the contactless SHG method combined with the Duffing-van der pol oscillator can effectively realize the convenient and reliable estimation of the nonlinear ultrasonic signal during the cycle of asphalt deterioration and healing, which is important for the property characterization of asphalt mixture considering its high attenuation to the ultrasonic waves.

Dynamics and energy harvesting from parametrically coupled self-excited electromechanical oscillator

The investigated parametrically coupled electromechanical structure is composed of a mechanical Duffing oscillator whose mass sits on a moving belt surface. The driving electrical network is a van der Pol oscillator whose aim is to actuate the attached DC motor to provide some rotatry unbalances and parametric coupling in the vibrating structure. The coupled oscillator is applied to energy harvesting and overcomes the limitation of low energy generation associated with a single oscillator of this kind. The system was solved analytically and validated by numerical methods. The global dynamics of the structure were investigated, and nonlinear phenomena such as Neimark–Sacker bifurcation, discontinuity-induced bifurcation, grazing–sliding, and bifurcation to multiple tori were identified. These nonlinear behaviors affect the harvested energy at bifurcation points, resulting in jumps from one energy level to another. In addition to harnessing the highest energy under hard parametric coupling, the coupling ensures that higher and more useful energy is harvested over a wider range of belt speeds. Finally, the qualitative validation of the numerical concept by experimental setup verifies the workings of the model.

The masking technique for forced nonlinear oscillator stability behavior analysis using the non-perturbative approach

The study utilized the masking technique to explore the stability behavior of a forced nonlinear oscillator through the non-perturbative approach, with a particular focus on a Van der Pol oscillator subjected to external force, characterized by both cubic and quadratic nonlinearities. The application of the non-perturbative method (NPM) in conjunction with the masking technique was a pivotal aspect of this research, transforming the inherently non-homogeneous, nonlinear system into a homogeneous linear system. This transformation was crucial as it simplified the complex dynamics of the system, rendering it more amenable to analysis. Through this method, the research successfully established the system’s overall frequency, meticulously accounting for the impact of the periodic external force. The study also identified a distinct type of resonance response, where the system’s frequency incorporates the excited frequency in a nonlinear relationship. The masking technique proved to be an invaluable tool for examining the stability behavior of forced vibrations in oscillators via the NPM, providing profound insights into stability under external forces and enhancing the understanding and control of oscillatory behaviors in nonlinear dynamical systems. A critical confirmation of the current methodology is provided by the remarkable agreement found between the numerical solution and the provided analytical solution. This agreement shows that the analytical method produces trustworthy predictions and appropriately describes the system’s behavior. The plotted stability diagrams, which demonstrate that the model’s simulation of stability behavior is consistent with observed events, particularly resonance phenomena, offer further validity for the findings. In the resonance case, the effects of the damping coefficient and the external force’s magnitude are significant. The results of the analysis show that an increase in the damping coefficient has a destabilizing effect that causes unstable zones to expand. In contrast, in the resonance state, the quadratic and cubic nonlinearity factors both contribute to stabilization. Understanding how various system factors impact stability dynamics particularly in relation to resonance phenomena is made easier with the help of this insight.

Model-Free Optimal Vibration Control of a Nonlinear System Based on Deep Reinforcement Learning

Optimal control of nonlinear vibration requires precise knowledge of the system and the solution to Hamilton–Jacobi–Bellman (HJB) equation. However, in practical engineering applications, acquiring precise system parameters poses challenges, and the analytical solutions for the HJB equation are difficult to obtain. In this paper, two reinforcement learning algorithms, Twin Delayed Deep Deterministic Policy Gradient (TD3) algorithm and Soft Actor-Critic (SAC) algorithm, are employed to train neural network-based optimal controllers for the van der Pol vibration system in the presence of unknown system parameters. To validate their performance, the controllers undergo testing in a series of experiments, including assessments of free vibration, frequency sweep excitation, and Gaussian noise excitation. The results indicate that both the TD3-trained and SAC-trained neural network-based controller are capable of proficiently suppress the vibration of the van der Pol oscillator. Additionally, these two model-free controllers can approximate the optimal control law which solved based on the dynamic model of the nonlinear system.

Two-dimensional Coupled Asymmetric van der Pol Oscillator

A reliable algorithm based on an adaptation of the standard differential transform method (DTM) is presented, which is the multi-step differential transform method (Ms-DTM) since it may be difficult to directly apply differential transform method (DTM) to obtain the series solutions for the present two- dimensional coupled asymmetric van der Pol oscillator. The solutions of a two dimensional coupled asymmetric van der Pol oscillator were obtained by Ms-DTM. Figurative comparisons between the Ms-DTM and the classical fourth order Runge-Kutta method (RK4) are given. The obtained results reveal that the proposed technique is a promising tool to solve the considered van der Pol oscillator and yield same information on the phase portrait confirming the stability of the system, effectively. It can be said that the considered approach can be easily extended to other nonlinear van der Pol oscillator systems and therefore is widely applicable in engineering and other sciences.

Frequency lock-in control and mitigation of nonlinear vortex-induced vibrations of an airfoil structure using a conserved-mass linear vibration absorber

We investigate the effectiveness of a vibration absorber on the vortex-induced vibration response of turbine blades during the frequency lock-in phase. A reduced-order model of a turbine blade and van der Pol oscillator is used to represent the fluid–blade interaction caused by the vortex shedding. A spring-mass-damper system is considered to model the vibration absorber. The advantage of the vibration absorber is demonstrated by simulating the nonlinear coupled four-degree-of-freedom aeroelastic system for the different sets of system parameters. We observe the dominance of a nonlinear vibration absorber over the linear vibration absorber only for the higher coupling parameter values. The analytical solution of the nonlinear coupled system is obtained through the method of multiple scales for the case of 1:1 internal resonance to identify the critical design parameters of the vibration absorber. We observe the high sensitivity of the system's frequency response to the distance of the vibration absorber from the elastic axis, along with the absorber's damping, stiffness, and mass. Finally, we perform a parametric analysis on the lock-in of the stability region to better understand the effect of the vibration absorber on the instability region.

A multi-modal nonlinear dynamic model to investigate time-domain responses of a micro-cantilever in fluids

In this current work, a new nonlinear dynamic model based on the forced Van der Pol oscillator is introduced to demonstrate the time-domain sensitivities of the micro-cantilever to the varying properties of the surrounding fluids. Effects of diverse multi-frequency excitations on the hydrodynamically forced displacements are investigated for the Glycerol-water solutions with different concentrations. Driving forces at the eigenmode frequencies are applied simultaneously to actuate the micro-cantilever in multi-modal operations. The hydrodynamic force induces notable variations in the observables of high-frequency steady-state vibrations. To illustrate, the frequency of the displacements decreases with increasing dynamic viscosity and density of the fluids (among 55% and 85% Glycerol-water solutions) in bimodal- and trimodal-frequency excitations. Essentially, the observable responses are often used to distinguish the surrounding fluids in which the micro-cantilever operates. In addition, steady-state observables are achieved at only particular eigenmodes in single- and multi-frequency operations. It is highlighted that the periodic oscillations are obtained for the first and second eigenmodes with the highest value of forced Van der Pol parameter (μ = 1030). Clearly, higher eigenmodes require different values of the nonlinearity parameter to acquire periodic vibrations in multi-modal operations. In general, achieving steady-state observables is substantially critical in quantifying sensitivity to varying fluid properties. For instance, the vibration frequency of around 7.33 MHz and amplitude of around 0.03 pm are obtained at the first eigenmode for 75% Glycerol-water solution in tetra-modal operations. Note that femtometer amplitudes of deflections can be measured using quantum-enhanced AFM techniques. The frequency responses obtained in this work are compared with the measured ones in the literature and the results show satisfactory agreements. Therefore, a novel multi-modal nonlinear dynamic model enables to quantify observable sensitivity to micro-rheological properties at higher eigenmodes of the micro-cantilever.

Computer-algebraic approach to first differential approximation: Van der Pol oscillator

First differential approximation has been used to analyze various numerical methods for solving systems of ordinary differential equations. This has made it possible to estimate the stiffness of the ODE system that models the oscillations of the Van der Pol oscillator and the error of the method as well as to propose simple criteria for choosing a calculation step. The presented methods allow one to perform efficient calculations using computer algebra systems.

Impact of coupling on the road to synchronization of two coupled Van der Pol oscillators

Synchronized dynamics of two coupled van der Pol oscillators x and y balances between the attracting forces and the prevention of synchrony represented by the coupling μ and the half difference of their natural frequencies Δω. The current paper investigates two regimes of such balance: (i) μ slightly exceeds Δω, and (ii) μ is large. In case of regime (i), the oscillators are shown to quickly attain the neighborhood of the limit cycle of a complex geometry; their pre-limit behavior can be characterized by the vertical mismatch x−y. If μ→∞, case (ii), the shape of the solutions of the underlying differential equations promptly follows that of the limit cycle, which is the solution of a single van der Pol equation, but the tendency to the limit cycle is slow, being proportional to 1/μ. Further, the mismatch between the oscillators disappears as μ increases. During the transition between the regimes occurred with the growth of μ, the vertical mismatch and the time required to solutions to reach the neighborhood of the limit cycle exhibit the Λ- and V-shapes. These shapes attain the maximum and, respectively, minimum at moderate values of μ that were earlier suggested by studies of synchronization between solar activity components.

Theoretical Study of Forced Van Der Pol Oscillator Equation Using Multiple Two-timing Regular Parameter Perturbation and Asymptotic Expansion Techniques

This paper presents the theoretical study of forced Van der Pol oscillator equation. Oscillatory systems are studied to know measures that can reduce the amplitude of oscillation of the oscillatory system. Here, multiple two-timing regular parameter perturbation is applied since it is a kind of perturbation among other perturbation techniques that enables the study of the behaviour of a system under certain conditions. Asymptotic expansion technique was also applied. Excel Microsoft was used to analyse the uniformly valid asymptotic solution of the Van der Pol oscillator equation obtained. The uniformly valid asymptotic solution in the independent variable obtained, showed that damping alters the amplitude of the oscillatory system thereby affecting its motion. Increase in damping decreases the amplitude of oscillation of the system. With damping incorporated in the system though very small damping, the amplitude of oscillation reduces with time.

Optimal Parameter Estimation Techniques for Complex Nonlinear Systems

AbstractAccurate parameter estimation and state identification within nonlinear systems are fundamental challenges addressed by optimization techniques. This paper fills a critical gap in previous research by investigating tailored optimization methods for parameter estimation in nonlinear system modeling, with a particular emphasis on chaotic dynamical systems. We introduce and compare three optimization methods: a gradient-based iterative algorithm, the Levenberg-Marquardt algorithm, and the Nelder-Mead simplex method. These methods are strategically employed to simplify complex nonlinear optimization problems, rendering them more manageable. Through a comprehensive exploration of the performance of these methods in determining parameters across diverse systems, including the van der Pol oscillator, the Rössler system, and pharmacokinetic modeling, our study revealed that the accuracy and reliability of the Nelder-Mead simplex method were consistent. The Nelder-Mead simplex algorithm emerged as a powerful tool, that consistently outperforms alternative methods in terms of root mean squared error (RMSE) and convergence reliability. Visualizations of trajectory comparisons and parameter convergence under various noise levels further emphasize the algorithm’s robustness. These studies suggest that the Nelder-Mead simplex method has potential as a valuable tool for parameter estimation in chaotic dynamical systems. Our study’s implications extend beyond theoretical considerations, offering promising insights for parameter estimation techniques in diverse scientific fields reliant on nonlinear system modeling.

Λ -Fractional Forced Van der Pol Oscillator with Horizon

Extending the investigation of the modified Van der Pol equation with its chaotic behaviour into the context of Λ -fractional analysis, the solutions of the Λ-fractional Van der Pol oscillator are studied under the influence of a horizon restricting the region of non-locality. The various solutions to the equation depend upon the horizon, however, the chaotic behaviour of the oscillator mainly depends upon the fractional order. A comparison of the chaotic behaviour of the various Van der Pol oscillator systems is presented.

Computer-Algebraic Approach to First Differential Approximations: Van der Pol Oscillator

Computer-Algebraic Approach to First Differential Approximations: Van der Pol Oscillator

Local Lipschitz Continuity in the Initial Value and Strong Completeness for Nonlinear Stochastic Differential Equations

Recently, Hairer et al. (2015) showed that there exist stochastic differential equations (SDEs) with infinitely often differentiable and globally bounded coefficient functions whose solutions fail to be locally Lipschitz continuous in the strong L p L^p -sense with respect to the initial value for every p ∈ ( 0 , ∞ ] p\in (0,\infty ] . In this article we provide conditions on the coefficient functions of the SDE and on p ∈ ( 0 , ∞ ] p \in (0,\infty ] that are sufficient for local Lipschitz continuity in the strong L p L^p -sense with respect to the initial value and we establish explicit estimates for the local Lipschitz continuity constants. In particular, we prove local Lipschitz continuity in the initial value for several nonlinear stochastic ordinary and stochastic partial differential equations in the literature such as the stochastic van der Pol oscillator, Brownian dynamics, the Cox-Ingersoll-Ross processes and the Cahn-Hilliard-Cook equation. As an application of our estimates, we obtain strong completeness for several nonlinear SDEs.

Novel approach for solving higher-order differential equations with applications to the Van der Pol and Van der Pol–Duffing equations

BackgroundNumerical methods are used to solve differential equations, but few are effective for nonlinear ordinary differential equations (ODEs) of order higher than one. This paper proposes a new method for such ODEs, based on Taylor series expansion. The new method is a second-order method for second-order ODEs, and it is equivalent to the central difference method, a well-known method for solving differential equations. The new method is also simple to implement for higher-order differential equations. The proposed technique was applied to solve the Van der Pol and Van der Pol–Duffing equations. It is stable over a wide range of nonlinearity and produces accurate and reliable results. For the self-excitation Van der Pol equation, the proposed technique was applied with different values of nonlinear damping.ResultsThe results were compared with those obtained using the ODE15s solver in MATLAB. The two sets of results showed excellent agreement. For the forced Van der Pol–Duffing equation, the proposed technique was applied with different values of exciting force amplitude and frequency. It was found that for certain conditions, the solution obtained using the proposed technique differed from that obtained using ODE15s.ConclusionsThe solution obtained using the proposed technique showed good agreement with the solutions obtained using ODE45 and Runge–Kutta fourth order. The results show that the proposed approach is very simple to apply and produces acceptable error. It is a powerful and versatile tool for solving of high-order nonlinear differential equations accurately.

Combined vibration analysis of a top tensioned riser with the geometrical nonlinearity

A coupling model representing the interactions between the riser and vortices is proposed, of which the time-varying top tension and the internal nonlinear axial force due to the bending of the riser are considered. In addition, the van der Pol oscillator is introduced to simulate the time-varying characteristic of the vortices. The motion equations are derived and the first-order mode approximation is obtained with the Galerkin approach. The multiple scale method is applied to study the steady-state solutions of the system. Effects of system parameters (the structural damping, nonlinearity, the amplitude ratio of the varying-tension, the shedding frequency) on the responses are investigated in detail. The synchronous and non-synchronous motions between the structure and wake oscillator are studied. The bifurcation diagrams for the responses in terms of the amplitude ratio and the nonlinear parameter are obtained in the large parameter ranges. The results show that responses with different topological properties including quasi-periodic, periodic and chaotic solutions can occur under different parameter conditions. A pair of chaotic attractors can be found for the large parameter range of the amplitude ratio. The joint effects of the amplitude ratio and the nonlinear parameter on the responses in the large parameter range are studied as well, which are further verified by the phase projections and Poincaré sections. The results show that the topological properties of responses can be controlled by the amplitude ratio. These results can be helpful to the dynamic design of the riser in practice.

Oscillator Simulation with Deep Neural Networks

The motivation behind this study is to overcome the complex mathematical formulation and time-consuming nature of traditional numerical methods used in solving differential equations. It seeks an alternative approach for more efficient and simplified solutions. A Deep Neural Network (DNN) is utilized to understand the intricate correlations between the oscillator’s variables and to precisely capture their dynamics by being trained on a dataset of known oscillator behaviors. In this work, we discuss the main challenge of predicting the behavior of oscillators without depending on complex strategies or time-consuming simulations. The present work proposes a favorable modified form of neural structure to improve the strategy for simulating linear and nonlinear harmonic oscillators from mechanical systems by formulating an ANN as a DNN via an appropriate oscillating activation function. The proposed methodology provides the solutions of linear and nonlinear differential equations (DEs) in differentiable form and is a more accurate approximation as compared to the traditional numerical method. The Van der Pol equation with parametric damping and the Mathieu equation are adopted as illustrations. Experimental analysis shows that our proposed scheme outperforms other numerical methods in terms of accuracy and computational cost. We provide a comparative analysis of the outcomes obtained through our proposed approach and those derived from the LSODA algorithm, utilizing numerical techniques, Adams–Bashforth, and the Backward Differentiation Formula (BDF). The results of this research provide insightful information for engineering applications, facilitating improvements in energy efficiency, and scientific innovation.

Explosive and semi-explosive transitions in parametrically perturbed systems

The phenomenon of transition between steady states, namely amplitude death and oscillation death, is identified in certain models of coupled parametrically driven nonlinear oscillators. Intriguingly, the nature of the transition is explosive. It is significant to note that this observation of explosive transition has been made in systems of identical coupled oscillators with repulsive mean-field coupling. In particular, we find that when the coupling strength increases, there are semi-explosive transitions in parametrically perturbed duffing oscillators and explosive transitions in parametrically perturbed van der Pol oscillators. We obtain analytical expressions for the critical bifurcation points through linear stability analysis. Detailed numerical studies have been performed to confirm the phenomenon.

Analysis of Successive Doubly Nested Mixed-Mode Oscillations

In previous works [Inaba & Kousaka, 2020; Inaba & Tsubone, 2020; Inaba et al., 2023], significant bifurcation structures referred to as nested Mixed-Mode Oscillations (MMOs) were found to be present in forced Bonhoeffer–van der Pol (BVP) oscillators. It is well known that unnested Mixed-Mode Oscillation-Incrementing Bifurcations (MMOIBs) can generate [Formula: see text] oscillations (i.e. [Formula: see text] followed by [Formula: see text] repeated [Formula: see text] times) for successive values of [Formula: see text], where [Formula: see text] and [Formula: see text] are adjacent fundamental simple MMOs, e.g. [Formula: see text] and [Formula: see text], where [Formula: see text] is an integer. Furthermore, it has been confirmed that MMOIBs can generate nested MMOs. Let two adjacent unnested MMOIB-generated MMOs be denoted [Formula: see text] and [Formula: see text]. Then, singly nested MMOIBs can generate [Formula: see text] for successive values of [Formula: see text], i.e. [Formula: see text] followed by [Formula: see text] repeated [Formula: see text] times, between the [Formula: see text]- and [Formula: see text]-generating regions. The sequential generation of singly nested MMOs has been investigated in detail in previous work [Ito et al., 2021]. Nested MMOs can, however, be at least doubly nested. In this study, we investigate doubly nested MMOs considering a constrained nonautonomous BVP oscillator containing an idealized diode. Based on the observed dynamics of this system, Poincaré return maps are rigorously constructed in one dimension. Therefore, we can solve the successive saddle-node bifurcations using a nested (i.e. double-loop) bisection method. We track 60 successive doubly nested MMOIBs and we do not rule out the possibility that the 58 scaling constants corresponding to the MMOIB intervals converge to unity. We note that because we solve the bifurcation equation avoiding the use of a method that requires the careful selection of the initial conditions (e.g. the Newton–Raphson), we can accurately track the saddle-node bifurcations without missing any doubly nested MMO sequences.