From Canards of Folded Singularities to Torus Canards in a Forced van der Pol Equation

In this article, we study canard solutions of the forced van der Pol equation in the relaxation limit for low-, intermediate-, and high-frequency periodic forcing. A central numerical observation made herein is that there are two branches of canards in parameter space which extend across all positive forcing frequencies. In the low-frequency forcing regime, we demonstrate the existence of primary maximal canards induced by folded saddle nodes of type I and establish explicit formulas for the parameter values at which the primary maximal canards and their folds exist. Then, we turn to the intermediate- and high-frequency forcing regimes and show that the forced van der Pol possesses torus canards instead. These torus canards consist of long segments near families of attracting and repelling limit cycles of the fast system, in alternation. We also derive explicit formulas for the parameter values at which the maximal torus canards and their folds exist. Primary maximal canards and maximal torus canards correspond geometrically to the situation in which the persistent manifolds near the family of attracting limit cycles coincide to all orders with the persistent manifolds that lie near the family of repelling limit cycles. The formulas derived for the folds of maximal canards in all three frequency regimes turn out to be representations of a single formula in the appropriate parameter regimes, and this unification confirms the central numerical observation that the folds of the maximal canards created in the low-frequency regime continue directly into the folds of the maximal torus canards that exist in the intermediate- and high-frequency regimes. In addition, we study the secondary canards induced by the folded singularities in the low-frequency regime and find that the fold curves of the secondary canards turn around in the intermediate-frequency regime, instead of continuing into the high-frequency regime. Also, we identify the mechanism responsible for this turning. Finally, we show that the forced van der Pol equation is a normal form-type equation for a class of single-frequency periodically driven slow/fast systems with two fast variables and one slow variable which possess a non-degenerate fold of limit cycles. The analytic techniques used herein rely on geometric desingularisation, invariant manifold theory, Melnikov theory, and normal form methods. The numerical methods used herein were developed in Desroches et al. (SIAM J Appl Dyn Syst 7:1131–1162, 2008, Nonlinearity 23:739–765 2010).     

Comparative dynamics of chains of coupled van der Pol equations and coupled systems of van der Pol equations

Theoretical and Mathematical Physics - We consider chains of van der Pol equations closed into a ring and chains of systems of two first-order van der Pol equations. We assume that the couplings...     

On the motion of a generalized van der Pol oscillator

A generalized van der Pol oscillator is considered, with positive real power nonlinearities in the restoring and damping force, including fractional powers. An analytical approach based on the Krylov–Bogoliubov method is adjusted to derive analytical expressions for the amplitude of a limit cycle for small values of the damping coefficient. These expressions are also derived for some integer power nonlinearities in the equation of motion and the results obtained compared with the existing results from the literature. Relaxation oscillations are studied for larger values of the damping coefficient. Matched asymptotic expansions are used and the influence of the powers of the restoring and damping force on the period of these oscillations is investigated. It is shown that not only can the period increase with the damping power, but it can also have a decreasing trend for some cases and the condition for this to hold is obtained.     

Stable Subharmonics of the Forced Van der Pol Equation

The Van der Pol equation with periodic forcing term is investigatednumerically for a range of values of the parameters. The equationis integrated using the fourth order Runge-Kutta method. Twodifferent types of subharmonics have been found, namely "normal"subharmonics and "abnormal" subharmonics. Finally, the criticalvalue b = is investigated.     

参数激励Duffing-VanderPol振子的动力学响应及反馈控制

研究了Duffing-Van der Pol振子的主参数共振响应及其时滞反馈控制问题．依平均法和对时滞反馈控制项Taylor展开的截断得到的平均方程表明,除参数激励的幅值和频率外,零解的稳定性只与原方程中线性项的系数和线性反馈有关,但周期解的稳定性还与原方程中非线性项的系数和非线性反馈有关．通过调整反馈增益和时滞,可以使不稳定的零解变得稳定．非零周期解可能通过鞍结分岔和Hopf分岔失去稳定性,但选择合适的反馈增益和时滞,可以避免鞍结分岔和Hopf分岔的发生．数值仿真的结果验证了理论分析的正确性．     

Time Optimal Control and the Van der Pol Oscillator

The problem considered is the control of the Van der Pol oscillator,where the control appears as a forcing term on the right-handside of the differential equation, the "goal" is the equilibriumstate and the "cost" is the time of transfer from an initialstate to the equilibrium state. Qualitative methods and some simple analytic procedures includinga Liapunov type stability theorem are used to determine domainsof controllability and switching curves.     

Construction of approximate periodic solutions to a modified van der Pol oscillator

A new iteration scheme is proposed and applied for the modified van der Pol oscillator. A simple and effective iteration procedure to search for the periodic solutions of the equation is given. This procedure is a powerful tool for the determination of the approximate frequencies and periodic solutions of the nonlinear differential equations. The solutions obtained using the present iteration procedure are in good agreement with the numerical integration obtained by a fourth order Runge–Kutta method, which shows the applicability of the procedure.     

Horizontal fast excitation in delayed van der Pol oscillator

This paper investigates the interaction effect of horizontal fast harmonic parametric excitation and time delay on self-excited vibration in van der Pol oscillator. We apply the method of direct partition of motion to derive the main autonomous equation governing the slow dynamic of the oscillator. The method of averaging is then performed on the slow dynamic to obtain a slow flow which is analyzed for equilibria and periodic motion. This analysis provides analytical approximations of regions in parameter space where periodic self-excited vibrations can be eliminated. A numerical study is performed on the original oscillator and compared to analytical approximations. It was shown that in the delayed case, horizontal fast harmonic excitation can eliminate undesirable self-excited vibrations for moderate values of the excitation frequency. In contrast, the case without delay requires large excitation frequency to eliminate such motions. This work has application to regenerative behavior in high-speed milling.     

Optimal feedback control of the forced van der Pol system

A simple feedback control strategy for chaotic systems is investigated using the forced van der Pol system as an example. The strategy regards chaos control as an optimization problem, where the maximum magnitude Floquet multiplier of a target unstable periodic orbit (UPO) is used as a cost function that needs to be minimized. Thus, the method obtains the optimal control gain in terms of the stability of the target UPO. This strategy was recently proposed for the proportional feedback control (PFC) method. Here, it is extended to the highly popular delayed feedback control (DFC) method. Since the DFC method treats the system as a delay-differential equation whose phase space is infinite-dimensional, the characteristic multipliers are found through a truncation in the number of delayed states. Control of a target UPO is achieved for several values of the forcing amplitude. We compare the DFC and PFC methods in terms of stability of the controlled orbit, steady state error and control effort.     

Existence of limit cycles for coupled van der Pol equations

In this paper, we consider the existence of limit cycles of coupled van der Pol equations by using S¹-degree theory due to Dylawerski et al. (see Ann. Polon. Math. 62 (1991) 243).     

On limit cycle approximations in the van der Pol oscillator

This paper applies bifurcation analysis to the well-known van der Pol oscillator to obtain approximations of its periodic solutions in the nearly sinusoidal regime. A frequency domain method based on harmonic balance approximations is used for small values of the bifurcation parameter. Moreover, a comparison with some other frequency domain approaches is also given. Finally, a total harmonic distortion is computed using the information provided by the frequency domain approach.     

Models and numerical schemes for generalized van der Pol equations

This paper presents three generalizations of the van der Pol equation (VDPE) using newly proposed three new generalized K-, A- and B-operators. These operators allow kernel to be arbitrary. As a result, these operators provide a greater generalization of the VDPE than the fractional integral and differential operators do. Like the original VDPE, the generalized van der Pol equations (GVDPEs) are also nonlinear equations, and in most cases, they can not be solved analytically. Numerical algorithms are presented and used to solve the GVDPEs. Results for several examples are presented to demonstrate the effectiveness of the numerical algorithms, and to examine the behavior of the GVDPEs and the limit cycles associated with them. Although the numerical algorithms have been used to solve the GVDPEs only, they can also be used to solve many other generalized oscillators and generalized differential equations. This will be considered in the future.     

Numerical Scaling Invariance Applied to the van der Pol Model

In this study we use the van der Pol model to explain a novel numerical application of scaling invariance. The model in point is not invariant to a scaling group of transformations, but by introducing an embedding parameter we are able to recover it from an extended model which is invariant to an extended scaling group. As well known, within a similarity analysis we can define a family of solutions from a computed one, so that the solution of a target problem can be obtained by rescaling the solution of a reference problem. The main idea is to use scaling invariance and numerical analysis to find a reference problem easier to solve, from a numerical viewpoint, than the target problem. This allows us to save human efforts and computational resources every time we have to solve a challenging problem. We test our approach using three stiff solvers available within the most recent releases of MATLAB. Independently from the solver used, by employing the described scaling invariance we are able to significantly reduce the computational cost of the numerical solution of the van der Pol model.      

Synchronization of van der Pol oscillator and Chen chaotic dynamical system

This paper addresses the synchronization problem of two different electronic circuits by using nonlinear control function. This technique is applied to achieve synchronization for the stable van der Pol oscillator and Chen chaotic dynamical system. Numerical simulations results are given to demonstrate the effectiveness of the proposed control method.     

Invariant manifolds of the Bonhoeffer–van der Pol oscillator

The stable and unstable manifolds of a saddle fixed point (SFP) of the Bonhoeffer–van der Pol oscillator are numerically studied. A correspondence between the existence of homoclinic tangencies (which are related to the creation or destruction of Smale horseshoes) and the chaos observed in the bifurcation diagram is described. It is observed that in the non-chaotic zones of the bifurcation diagram, there may or may not be Smale horseshoes, but there are no homoclinic tangencies.     

Hypergeometric first integrals of the Duffing and van der Pol oscillators

The autonomous Duffing oscillator, and its van der Pol modification, are known to admit time-dependent first integrals for specific values of parameters. This corresponds to the existence of Darboux polynomials, and in fact more can be shown: that there exist Liouvillian first integrals which do not depend on time. They can be expressed in terms of the Gauss and Kummer hypergeometric functions, and are neither analytic, algebraic nor meromorphic. A criterion for this to happen in a general dynamical system is formulated as well.     

Stability and bifurcation analysis in the delay-coupled van der Pol oscillators

In this paper, the dynamics of a system of two van der Pol equations with a finite delay are investigated. We show that there exist the stability switches and a sequence of Hopf bifurcations occur at the zero equilibrium when the delay varies. Using the theory of normal form and the center manifold theorem, the explicit expression for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived.