A Dissipation-Preserving Integrator for Damped Oscillatory Hamiltonian Systems

In this paper, based on discrete gradient, a dissipation-preserving integrator for weakly dissipative perturbations of oscillatory Hamiltonian system is established. The solution of this system is a damped nonlinear oscillator. Basically, lots of nonlinear oscillatory mechanical systems including frictional forces lend themselves to this approach. The new integrator gives a discrete analogue of the dissipation property of the original system. Meanwhile, since the integrator is based on the variation-of-constants formula for oscillatory systems, it preserves the oscillatory structure of the system. Some properties of the new integrator are derived. The convergence is analyzed for the implicit iterations based on the discrete gradient integrator, and it turns out that the convergence of the implicit iterations based on the new integrator is independent of $\|M\|$, where $M$ governs the main oscillation of the system and usually $\|M\|\gg1$. This significant property shows that a larger stepsize can be chosen for the new schemes than that for the traditional discrete gradient integrators when applied to the oscillatory Hamiltonian system. Numerical experiments are carried out to show the effectiveness and efficiency of the new integrator in comparison with the traditional discrete gradient methods in the scientific literature.     

Oscillation for a second-order neutral differential equation with impulses

We consider a certain type of second-order neutral delay differential systems and we establish two results concerning the oscillation of solutions after the system undergoes controlled abrupt perturbations (called impulses). As a matter of fact, some particular non-impulsive cases of the system are oscillatory already. Thus, we are interested in finding adequate impulse controls under which our system remains oscillatory.     

The self-induced rotational and oscillatory motions of an aerodynamic pendulum

The behaviour of an oscillatory mechanical system with alternating dissipation is considered taking an aerodynamic pendulum as an example. The phase portraits are investigated, their rearrangements are studied and the critical values of a parameter are determined. The equilibrium positions of the pendulum and the self induced rotational and oscillatory states of the motion are determined and their stability is investigated.     

Long-Time Oscillatory Energy Conservation of Total Energy-Preserving Methods for Highly Oscillatory Hamiltonian Systems

For an integrator when applied to a highly oscillatory system,the near conservation of the oscillatory energy over long times is an important aspect.In this paper,we study the long-time near conservation of oscillatory energy for the adapted average vector field(AAVF)method when applied to highly oscillatory Hamiltonian systems.This AAVF method is an extension of the average vector field method and preserves the total energy of highly oscillatory Hamiltonian systems exactly.This paper is devoted to analysing another important property of AAVF method,i.e.,the near conservation of its oscillatory energy in a long term.The long-time oscillatory energy conservation is obtained via constructing a modulated Fourier expansion of the AAVF method and deriving an almost invariant of the expansion.A similar result of the method in the multi-frequency case is also presented in this paper.     

A qualitative investigation of the oscillations of a pendulum with a periodically varying length and a mathematical model of a swing

The behaviour of the amplitude-frequency characteristics of families of periodic solutions, produced from the equilibrium position of a system, is established by a qualitative investigation of the equation of the oscillations of a pendulum, the length of which is an arbitrary periodic function of time. The non-local conditions for their stability and instability, expressed in terms of the amplitude and frequency of the oscillations, are obtained. The results are used when discussing the parametric and self-excited oscillatory model of a swing. In the parametric model the length of a swing is a specified periodic function of time, and in the self-excited oscillatory model it is a function of the phase coordinates of the system. For an appropriate choice of these functions, both systems have a common periodic solution. It is shown that the parametric model leads to an erroneous conclusion regarding the instability of the periodic mode, which is in fact realized in the oscillations of a swing, whereas the self-excited oscillatory model indicates its stability.     

血液流动与血管壁运动

本文讨论了哺乳动物循环系统的血液流动与血管壁运动之间的相互作用问题.在假定流动处于稳定的振荡流动情况下,导得了一组血液流动速度分布公式,压力分布公式以及约束应力公式,管壁位移公式.把Kuchar的公式从定常流动情况推广到非定常的振荡流动情况.文中还讨论了动脉血管壁的弹性效应问题.     

Dynamics of bidirectional associative memory networks with distributed delays and reaction–diffusion terms

In this paper, the periodic oscillatory solution and stability are investigated for a class of bidirectional associative memory neural networks with distributed delays and reaction–diffusion terms. By constructing a new Lyapunov functional, applying M-matrix theory and inequality technique, several novel sufficient conditions are derived to ensure the existence and uniqueness of periodic oscillatory solutions for bidirectional associative memory neural networks with distributed delays and reaction–diffusion terms, and all other solutions of this network converge exponentially to the unique periodic oscillatory solution. Moreover, the exponential convergence rate is estimated, which depends on the delay kernel functions and the system parameters. Two numerical examples are given to show the effectiveness of the obtained results. The results extend and improve the previously known results.     

Moving gap solitons in periodic potentials

We address the existence of moving gap solitons (traveling localized solutions) in the Gross–Pitaevskii equation with a small periodic potential. Moving gap solitons are approximated by the explicit solutions of the coupled‐mode system. We show, however, that exponentially decaying traveling solutions of the Gross–Pitaevskii equation do not generally exist in the presence of a periodic potential due to bounded oscillatory tails ahead and behind the moving solitary waves. The oscillatory tails are not accounted in the coupled‐mode formalism and are estimated by using techniques of spatial dynamics and local center‐stable manifold reductions. Existence of bounded traveling solutions of the Gross–Pitaevskii equation with a single bump surrounded by oscillatory tails on a large interval of the spatial scale is proven by using these techniques. Copyright © 2008 John Wiley & Sons, Ltd.     

Existence and stability of the relaxation cycle in a mathematical repressilator model

The three-dimensional nonlinear system of ordinary differential equations modeling the functioning of the simplest oscillatory genetic network, the so-called repressilator, is considered. The existence, asymptotics, and stability of the relaxation periodicmotion in this system are studied.     

Discrete nonlinear excitations on an inhomogeneous lattice in a spatially linear potential

Discrete nonlinear excitations on an inhomogeneous lattice in a spatially linear potential are investigated. Discrete breathers possessing multiple frequencies are always spatially localized and oscillatory in time, with oscillatory frequencies depending on an external potential. A number of phenomena associated with the translational symmetry of the system is addressed. The spatially localized, time periodic and quasi-periodic nonlinear coherent excitations comprising multiple-soliton components are presented with exact results for the corresponding energy, momentum and norm. The mean width of a pulse as a measure of localization is also studied. Project supported by the National Natural Science Foundation of China.     

Qualitative analysis and solutions of bounded traveling wave for B-BBM equation

In this paper, we apply the theory of planar dynamical systems to carry out qualitative analysis for the dynamical system corresponding to B-BBM equation, and obtain global phase portraits under various parameter conditions. Then, the relations between the behaviors of bounded traveling wave solutions and the dissipation coeffiicient μ are investigated. We find that a bounded traveling wave solution appears as a kink profile solitary wave solution when μ is more than the critical value obtained in this paper, while a bounded traveling wave solution appears as a damped oscillatory solution when μ is less than it. Furthermore, we explain the solitary wave solutions obtained in previous literature, and point out their positions in global phase portraits. In the meantime, approximate damped oscillatory solutions are given by means of undetermined coefficients method. Finally, based on integral equations that reflect the relations between the approximate damped oscillatory solutions and the implicit exact damped oscillatory solutions, error estimates for the approximate solutions are presented.     

On the Method of Neumann Series for Highly Oscillatory Equations

The main purpose of this paper is to describe and analyse techniques for the numerical solution of highily oscillatory ordinary differential equations by exployting a Neumann expansion. Once the variables in the differential system are changed with respect to a rapidly rotating frame of reference, the Neumann method becomes very effective indeed. However, this effectiveness rests upon suitable quadrature of highly oscillatory multivariate integrals, and we devote part of this paper to describe how to accomplish this to high accuracy with a modest computational effort.     

基于修正的Magnus方法的高振荡动力系统的数值积分方法

基于建立于一般线性动力系统上的Magnus数值积分方法,针对随时间而高频率振荡的二阶动力系统,给出了有效的修正Magnus数值积分算法．首先,将二阶动力系统重新表示为一阶系统的形式,通过引进新变量进行参考坐标变换,使动力系统的高振荡性质保留在新形式内;进而基于局部线性化技术用修正的Magnus方法求解新形式下的系统方程;最后,通过一系列数值实验说明了文中方法的有效性．     

Single-frequency oscillations of non-linear systems with distributed parameters

A weakly non-linear oscillatory system with distributed parameters is investigated An asymptotic method of constructing a solution, which describes the oscillatory motions of the single-mode (single-frequency) approximation, which is usually implemented in practical problems, is described and justified. Constructive sufficient conditions are formulated and the closeness of the approximate single-frequency solution to the exact solution in an asymptotically long time interval is proved. Possible extensions of the structure of the perturbing functions are considered and the case of the finite-mode approximation is investigated. Solutions of specific problems, which are of practical interest, are constructed to illustrate the effectiveness of the single-frequency approximation method.     

Generalized even and odd totally positive matrices

A generalization of the definition of an oscillatory matrix based on the theory of cones is given in this paper. The positivity and simplicity of all the eigenvalues of a generalized oscillatory matrix are proved. Classes of generalized even and odd oscillatory matrices are introduced. Spectral properties of the obtained matrices are studied. Criteria of generalized even and odd oscillation are given. Examples of generalized even and odd oscillatory matrices are presented.     

Fast,numerically stable computation of oscillatory integrals with stationary points

We present a numerically stable way to compute oscillatory integrals. For each additional frequency, only a small, well-conditioned linear system with a Hessenberg matrix must be solved, and the amount of work needed decreases as the frequency increases. Moreover, we can modify the method for computing oscillatory integrals with stationary points. This is the first stable algorithm for oscillatory integrals with stationary points which does not lose accuracy as the frequency increases and does not require deformation into the complex plane.     

Parameter estimation using dynamic optimization in limited nutrient-consumer hybrid dynamical systems

A limited nutrient consumer hybrid system which exhibits rich dynamics such as stability, instability and oscillatory characteristics is considered. The parameters that strongly influence the behavior of the system are identified and are estimated using the techniques of dynamic optimization. Two independent bio-control mechanisms are introduced into the system when the basic system exhibits instability characteristics. Values of the parameters representing these mechanisms are estimated at which the system stabilizes. Examples and simulations illustrate the main results.     

Bounded traveling waves of the Burgers-Huxley equation

In order to investigate bounded traveling waves of the Burgers-Huxley equation, bifurcations of codimension 1 and 2 are discussed for its traveling wave system. By reduction to center manifolds and normal forms we give conditions for the appearance of homoclinic solutions, heteroclinic solutions and periodic solutions, which correspondingly give conditions of existence for solitary waves, kink waves and periodic waves, three basic types of bounded traveling waves. Furthermore, their evolutions are discussed to investigate the existence of other types of bounded traveling waves, such as the oscillatory traveling waves corresponding to connections between an equilibrium and a periodic orbit and the oscillatory kink waves corresponding to connections of saddle-focus.     

Preconditioned defect corrected averaging for parabolic PDEs

Modelling physical systems with fast moving components leads to PDEs with highly oscillatory sources. Often, the time scale of the oscillation is much below the scale of the interesting variables. Time integrators must follow the scale of the fast motion, leading to long simulation times. For mechanical systems, methods like heterogeneous multiscaling and stroboscopic averaging are quite satisfactory. In case of semidiscretized PDEs, their advantage is limited. Here, we derive a smooth source term which generates a solution that coincides with the solution of the oscillatory system in stroboscopic points. The derivation involves the solution of a linear system with the solution operator of the PDE being the linear operator. Several preconditioners are developed and compared for those systems. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Onset of triply diffusive convection in a Maxwell fluid saturated porous layer

The linear stability of triply diffusive convection in a binary Maxwell fluid saturated porous layer is investigated. Applying the normal mode method theory, the criterion for the onset of stationary and oscillatory convection is obtained. The modified Darcy–Maxwell model is used as the analysis model, this allows us to make a thorough investigation of the processes of viscoelasticity and diffusions that causes the convection to set in through oscillatory rather than stationary. The effects of Vadasz number, normalized porosity parameter, relaxation parameter, Lewis number and solute Rayleigh number on the system are presented numerically and graphically.