Periodic wave solutions of coupled integrable dispersionless equations by residue harmonic balance

We introduce the residue harmonic balance method to generate periodic solutions for nonlinear evolution equations. A PDE is firstly transformed into an associated ODE by a wave transformation. The higher-order approximations to the angular frequency and periodic solution of the ODE are obtained analytically. To improve the accuracy of approximate solutions, the unbalanced residues appearing in harmonic balance procedure are iteratively considered by introducing an order parameter to keep track of the various orders of approximations and by solving linear equations. Finally, the periodic solutions of PDEs result. The proposed method has the advantage that the periodic solutions are represented by Fourier functions rather than the sophisticated implicit functions as appearing in most methods.     

Oscillatory region and asymptotic solution of fractional van der Pol oscillator via residue harmonic balance technique

In this paper, a powerfully analytical technique is proposed for predicting and generating the steady state solution of the fractional differential system based on the method of harmonic balance. The zeroth-order approximation using just one Fourier term is applied to predict the parametric function for the boundary between oscillatory and non-oscillatory regions of the fractional van der Pol oscillator. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear algebraic equations to improve the accuracy of the solutions successively. The highly accurate solutions to the angular frequency and limit cycle of fractional van der Pol oscillator are obtained and compared. The results reveal that the technique described in this paper is very effective and simple for obtaining asymptotic solution of nonlinear system having fractional order derivative.     

Newton–harmonic balancing approach for accurate solutions to nonlinear cubic–quintic Duffing oscillators

This paper presents a new approach for solving accurate approximate analytical higher-order solutions for strong nonlinear Duffing oscillators with cubic–quintic nonlinear restoring force. The system is conservative and with odd nonlinearity. The new approach couples Newton’s method with harmonic balancing. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton’s method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of nonlinear algebraic equations without analytical solution. Using the approach, accurate higher-order approximate analytical expressions for period and periodic solution are established. These approximate solutions are valid for small as well as large amplitudes of oscillation. In addition, it is not restricted to the presence of a small parameter such as in the classical perturbation method. Illustrative examples are presented to verify accuracy and explicitness of the approximate solutions. The effect of strong quintic nonlinearity on accuracy as compared to cubic nonlinearity is also discussed.     

The effect of time delay on the dynamics of an SEIR model with nonlinear incidence

In this paper, a SEIR epidemic model with nonlinear incidence rate and time delay is investigated in three cases. The local stability of an endemic equilibrium and a disease-free equilibrium are discussed using stability theory of delay differential equations. The conditions that guarantee the asymptotic stability of corresponding steady-states are investigated. The results show that the introduction of a time delay in the transmission term can destabilize the system and periodic solutions can arise through Hopf bifurcation when using the time delay as a bifurcation parameter. Applying the normal form theory and center manifold argument, the explicit formulas determining the properties of the bifurcating periodic solution are derived. In addition, the effect of the inhibitory effect on the properties of the bifurcating periodic solutions is studied. Numerical simulations are provided in order to illustrate the theoretical results and to gain further insight into the behaviors of delayed systems.     

Forward residue harmonic balance for autonomous and non-autonomous systems with fractional derivative damping

Both the autonomous and non-autonomous systems with fractional derivative damping are investigated by the harmonic balance method in which the residue resulting from the truncated Fourier series is reduced iteratively. The first approximation using a few Fourier terms is obtained by solving a set of nonlinear algebraic equations. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear algebraic equations to improve the accuracy and increase the number of Fourier terms of the solutions successively. Multiple solutions, representing the occurrences of jump phenomena, supercritical pitchfork bifurcation and symmetry breaking phenomena are predicted analytically. The interactions of the excitation frequency, the fractional order, amplitude, phase angle and the frequency amplitude response are examined. The forward residue harmonic balance method is presented to obtain the analytical approximations to the angular frequency and limit cycle for fractional order van der Pol oscillator. Numerical results reveal that the method is very effective for obtaining approximate solutions of nonlinear systems having fractional order derivatives.     

On Bifurcating Periodic Solutions at Low Frequency

A study is made of differential equations which have low-frequency periodic bifurcating solutions. It is shown that the classical Poincaré-Hopf technique for constructing periodic solutions encounters difficulties in the low-frequency limit, and that the branches are determined by nonlinear, rather than linear, balances. Two types of models are investigated: one is autonomous and the other non-autonomous, with a forcing term of small, fixed frequency. The latter model is relevant to several fluid-dynamical situations.     

Fractional derivative and time delay damper characteristics in Duffing–van der Pol oscillators

In this paper, we investigate the damping characteristics of two Duffing–van der Pol oscillators having damping terms described by fractional derivative and time delay respectively. The residue harmonic balance method is presented to find periodic solutions. No small parameter is assumed. Highly accurate limited cycle frequency and amplitude are captured. The results agree well with the numerical solutions for a wide range of parameters. Based on the obtained solutions, the damping effects of these two oscillators are investigated. When the system parameters are identical, the steady state responses and their stability are qualitatively different. The initial approximations are obtained by solving a few harmonic balance equations. They are improved iteratively by solving linear equations of increasing dimension. The second-order solutions accurately exhibit the dynamical phenomena when taking the fractional derivative and time delay as bifurcation parameters respectively. When damping is described by time delay, the stable steady state response is more complex because time delay takes past history into account implicitly. Numerical examples taking time delay and fractional derivative are respectively given for feature extraction and convergence study.     

Global residue harmonic balance method to periodic solutions of a class of strongly nonlinear oscillators

A new approach, namely the global residue harmonic balance method, was advanced to determine the accurate analytical approximate periodic solution of a class of strongly nonlinear oscillators. A class of nonlinear jerk equation containing velocity-cubed and velocity times displacements-squared was taken as a typical example. Unlike other harmonic balance methods, all the former residual errors are introduced in the present approximation to improve the accuracy. Comparison of the result obtained using this approach with the exact one and simplicity and efficiency of the proposed procedure. The method can be easily extended to other strongly nonlinear oscillators.     

Local Hopf bifurcation and global existence of periodic solutions in a kind of physiological system

The dynamics of a physiological control systems described by a first-order nonlinear delay differential equations are investigated. we proved that a sequence of Hopf bifurcations occur at the positive equilibrium as the delay increases. Explicit algorithm for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived, using the theory of normal form and center manifold. Global existence of periodic solutions are established using a global Hopf bifurcation result due to Wu [Symmetric functional differential equations and neural networks with memory, Trans. Amer. Math. Soc. 350 (1998) 4799–4838].     

Nonlinear transversely vibrating beams by the improved energy balance method and the global residue harmonic balance method

Mathematical modeling of many engineering systems such as beam structures often leads to nonlinear ordinary or partial differential equations. Nonlinear vibration analysis of the beam structures is very important in mechanical and industrial applications. This paper presents the high order frequency-amplitude relationship for nonlinear transversely vibrating beams with odd and even nonlinearities using the improved energy balance method and the global residue harmonic balance method. The accuracy of the energy balance method is improved based on combining features of collocation method and Galerkin–Petrov method, and an improved harmonic balance method is proposed which is called the global residue harmonic balance method. Unlike other harmonic balance methods, all the former global residual errors are introduced in the present approximation to improve the accuracy. Finally, preciseness of the present analytic procedures is evaluated in contrast with numerical calculations methods, giving excellent results.     

A Volterra integral formulation for determining the periodic solutions of some autonomous, nonlinear, third-order ordinary differential equations

A Volterra integral formulation based on the introduction of a term proportional to the velocity times the square of the (unknown) frequency of oscillation, a new independent variable equal to the original one times the (unknown) frequency of oscillation, the method of variation of parameters and series expansions of both the solution and the frequency of oscillation, is used to determine the periodic solutions to three nonlinear, autonomous, third-order, ordinary differential equations. It is shown that the first term of the series expansion of the frequency of oscillation coincides with that determined from a first-order harmonic balance procedure, whereas the two-term approximation to the frequency of oscillation is shown to be more accurate than that of a parameter perturbation procedure for the second equation considered in this paper. For the third equation, it is shown that the two-term approximation presented in this paper is more accurate than the corresponding one for one of the parameter perturbation methods, and for initial velocities less than one, for the other parameter perturbation approach.     

Periodic Solution for Strongly Nonlinear Vibration Systems by He’s Energy Balance Method

This paper applies He’s Energy balance method (EBM) to study periodic solutions of strongly nonlinear systems such as nonlinear vibrations and oscillations. The method is applied to two nonlinear differential equations. Some examples are given to illustrate the effectiveness and convenience of the method. The results are compared with the exact solution and the comparison showed a proper accuracy of this method. The method can be easily extended to other nonlinear systems and can therefore be found widely applicable in engineering and other science.     

On the approximation of periodic solutions of equations with bounded nonlinearities

We present a method for obtaining a closed-form approximation to harmonic periodic solutions of a class of ordinary differential equations containing a bounded, Lipschitzian nonlinearity and a sinusoidal forcing term. Our technique replaces the continuous nonlinearity with a suitable step-function nonlinearity and uses the Phase-Shift-Averaging Method to write the solution of the piecewise-linear problem in closed form.     

Analytical and approximate solutions to autonomous,nonlinear, third-order ordinary differential equations

Analytical solutions to autonomous, nonlinear, third-order nonlinear ordinary differential equations invariant under time and space reversals are first provided and illustrated graphically as functions of the coefficients that multiply the term linearly proportional to the velocity and nonlinear terms. These solutions are obtained by means of transformations and include periodic as well as non-periodic behavior. Then, five approximation methods are employed to determine approximate solutions to a nonlinear jerk equation which has an analytical periodic solution. Three of these approximate methods introduce a linear term proportional to the velocity and a book-keeping parameter and employ a Linstedt–Poincaré technique; one of these techniques provides accurate frequencies of oscillation for all the values of the initial velocity, another one only for large initial velocities, and the last one only for initial velocities close to unity. The fourth and fifth techniques are based on the Galerkin procedure and the well-known two-level Picard’s iterative procedure applied in a global manner, respectively, and provide iterative/sequential approximations to both the solution and the frequency of oscillation.     

Steady state bifurcation of a periodically excited system under delayed feedback controls

This paper investigates the steady state bifurcation of a periodically excited system subject to time-delayed feedback controls by the combined method of residue harmonic balance and polynomial homotopy continuation. Three kinds of delayed feedback controls are considered to examine the effects of different delayed feedback controls and delay time on the steady state response. By means of polynomial homotopy continuation, all the possible steady state solutions corresponding the third-order superharmonic and second-subharmonic responses are derived analytically, i.e. without numerical integration. It is found that the delayed feedback changes the bifurcating curves qualitatively and possibly eliminates the saddle-node bifurcation during resonant. The delayed position-velocity coupling and the delayed velocity feedback controls can destabilize the steady state responses. Coexisting periodic solutions, period-doubling bifurcation and even chaos are found in these control systems. The neighborhood of the periodic solutions is verified numerically in the phase portraits. The various effects of time delay on the steady state response are investigated. Many new phenomena are observed.     

Almost periodic solutions for undamped nonhomogeneous delay-differential equations

We first establish a result giving conditions that certain undamped delay differential equations with almost periodic time dependence have unique almost periodic solutions. Using this result we obtain conditions that a second order scalar nonlinear delay differential equation with almost periodic forcing will have a unique almost periodic solution having saddle-type stability properties. These results use the method of averaging.

     

Approximate methods based on order reduction for the periodic solutions of nonlinear third-order ordinary differential equations

Four approximate methods based on order reduction, the introduction of a book-keeping parameter and power series expansions for the solution and the frequency of oscillation are used to analyze three autonomous, nonlinear, third-order ordinary differential equations which have analytical periodic solutions. The first method introduces the velocity in both sides of the equation if this (linear) term is not present. The second one is based on the first one but employs a new independent variable, whereas, in the third and fourth techniques, a term equal to the velocity times the square of the unknown frequency of oscillation is introduced in both sides of the equation. The third method uses the original independent variable, whereas the fourth one employs a new independent variable which depends linearly on the unknown frequency of oscillation. It is shown that the second method provides accurate solutions only for initial velocities close to unity, whereas the third one is found to yield very accurate results for the first and second equations considered here and only for large initial velocities for the third one. The fourth technique provides as accurate results as or more accurate results than parameter-perturbation techniques which deal with the third-order equations directly and are based on the expansion of certain constants that appear in the differential equations in terms of a book-keeping parameter.     

Stability Switches, Hopf Bifurcations, and Spatio-temporal Patterns in a Delayed Neural Model with Bidirectional Coupling

The dynamical behavior of a delayed neural network with bi-directional coupling is investigated by taking the delay as the bifurcating parameter. Some parameter regions are given for conditional/absolute stability and Hopf bifurcations by using the theory of functional differential equations. As the propagation time delay in the coupling varies, stability switches for the trivial solution are found. Conditions ensuring the stability and direction of the Hopf bifurcation are determined by applying the normal form theory and the center manifold theorem. We also discuss the spatio-temporal patterns of bifurcating periodic oscillations by using the symmetric bifurcation theory of delay differential equations combined with representation theory of Lie groups. In particular, we obtain that the spatio-temporal patterns of bifurcating periodic oscillations will alternate according to the change of the propagation time delay in the coupling, i.e., different ranges of delays correspond to different patterns of neural activities. Numerical simulations are given to illustrate the obtained results and show the existence of bursts in some interval of the time for large enough delay.     

Asymptotic analysis and accurate approximate solutions for strongly nonlinear conservative symmetric oscillators

A new accurate iterative and asymptotic method is introduced to construct analytical approximate solutions to strongly nonlinear conservative symmetric oscillators. The method is based on applying a second-order expansion with the harmonic balance method and it excludes the requirement of solving a set of coupled nonlinear algebraic equations. Newton's iteration or the linearized model may be readily deduced by considering only the first-order terms in the model. According to this iterative approach, only Fourier series expansions of restoring force function, its first- and second-order derivatives for each iteration are required. It is concluded here that by only one single iteration, very brief and yet accurate analytical approximate solutions can be attained. Three physical examples are solved and accurate solutions are presented to illustrate the physics of the system and the effectiveness of the proposed asymptotic method.     

Sine‐cosine wavelet method for fractional oscillator equations

Purpose In this article, a novel computational method is introduced for solving the fractional nonlinear oscillator differential equations on the semi‐infinite domain. The purpose of the proposed method is to get better and more accurate results. Design/methodology/approach The proposed method is the combination of the sine‐cosine wavelets and Picard technique. The operational matrices of fractional‐order integration for sine‐cosine wavelets are derived and constructed. Picard technique is used to convert the fractional nonlinear oscillator equations into a sequence of discrete fractional linear differential equations. Operational matrices of sine‐cosine wavelets are utilized to transformed the obtained sequence of discrete equations into the systems of algebraic equations and the solutions of algebraic systems lead to the solution of fractional nonlinear oscillator equations. Findings The convergence and supporting analysis of the method are investigated. The operational matrices contains many zero entries, which lead to the high efficiency of the method, and reasonable accuracy is achieved even with less number of collocation points. Our results are in good agreement with exact solutions and more accurate as compared with homotopy perturbation method, variational iteration method, and Adomian decomposition method. Originality/value Many engineers can utilize the presented method for solving their nonlinear fractional models.