Principal resonance responses of SDOF systems with small fractional derivative damping under narrow-band random parametric excitation

The principal resonance responses of nonlinear single-degree-of-freedom (SDOF) systems with lightly fractional derivative damping of order α (0 < α < 1) subject to the narrow-band random parametric excitation are investigated. The method of multiple scales is developed to derive two first order stochastic differential equation of amplitude and phase, and then to examine the influences of fractional order and intensity of random excitation on the first-order and second-order moment. As an example, the stochastic Duffing oscillator with fractional derivative damping is considered. The effects of detuning frequency parameter, the intensity of random excitation and the fractional order derivative damping on stability are studied through the largest Lyapunov exponent. The corresponding theoretical results are well verified through direct numerical simulations. In addition, the phenomenon of stochastic jump is analyzed for parametric principal resonance responses via finite differential method. The stochastic jump phenomena indicates that the most probable motion is around the larger non-trivial branch of the amplitude response when the intensity of excitation is very small, and the probable motion of amplitude responses will move from the larger non-trivial branch to trivial branch with the increasing of the intensity of excitation. Such stochastic jump can be considered as bifurcation.     

Approximate solutions of a nonlinear oscillator typified as a mass attached to a stretched elastic wire by the homotopy perturbation method

The homotopy perturbation method is used to solve the nonlinear differential equation that governs the nonlinear oscillations of a system typified as a mass attached to a stretched elastic wire. The restoring force for this oscillator has an irrational term with a parameter λ that characterizes the system (0 ? λ ? 1). For λ = 1 and small values of x, the restoring force does not have a dominant term proportional to x. We find this perturbation method works very well for the whole range of parameters involved, and excellent agreement of the approximate frequencies and periodic solutions with the exact ones has been demonstrated and discussed. Only one iteration leads to high accuracy of the solutions and the maximal relative error for the approximate frequency is less than 2.2% for small and large values of oscillation amplitude. This error corresponds to λ = 1, while for λ < 1 the relative error is much lower. For example, its value is as low as 0.062% for λ = 0.5.     

Stick balancing with reflex delay in case of parametric forcing

The effect of parametric forcing on a PD control of an inverted pendulum is analyzed in the presence of feedback delay. The stability of the time-periodic and time-delayed system is determined numerically using the first-order semi-discretization method in the 5-dimensional parameter space of the pendulum’s length, the forcing frequency, the forcing amplitude, the proportional and the differential gains. It is shown that the critical length of the pendulum (that can just be balanced against the time-delay) can significantly be decreased by parametric forcing even if the maximum forcing acceleration is limited. The numerical analysis showed that the critical stick length about 30 cm corresponding to the unforced system with reflex delay 0.1 s can be decreased to 18 cm with keeping maximum acceleration below the gravitational acceleration.     

Eigenvalue problems for fractional ordinary differential equations

The eigenvalue problems are considered for the fractional ordinary differential equations with different classes of boundary conditions including the Dirichlet, Neumann, Robin boundary conditions and the periodic boundary condition. The eigenvalues and eigenfunctions are characterized in terms of the Mittag–Leffler functions. The eigenvalues of several specified boundary value problems are calculated by using MATLAB subroutine for the Mittag–Leffler functions. When the order is taken as the value 2, our results degenerate to the classical ones of the second-ordered differential equations. When the order α satisfies 1 < α < 2 the eigenvalues can be finitely many.     

Homotopy analysis of unsteady boundary layer flow adjacent to permeable stretching surface in a porous medium

This communication deals with the unsteady boundary layer flow of a viscous fluid in porous medium started due to the impulsively stretching of the plane wall. The wall is assumed to be porous so that suction or injection is possible. Complete analytic solution which is uniformly valid for all the dimensionless times 0 ⩽ τ < 0 in the whole spatial region 0 ⩽ η < ∞ is obtained by a purely analytic technique, namely the homotopy analysis method. Results are discussed through graphs.     

Dynamics of landslide model with time delay and periodic parameter perturbations

In present paper, we analyze the dynamics of a single-block model on an inclined slope with Dieterich–Ruina friction law under the variation of two new introduced parameters: time delay T_d and initial shear stress μ. It is assumed that this phenomenological model qualitatively simulates the motion along the infinite creeping slope. The introduction of time delay is proposed to mimic the memory effect of the sliding surface and it is generally considered as a function of history of sliding. On the other hand, periodic perturbation of initial shear stress emulates external triggering effect of long-distant earthquakes or some non-natural vibration source. The effects of variation of a single observed parameter, T_d or μ, as well as their co-action, are estimated for three different sliding regimes: β < 1, β = 1 and β > 1, where β stands for the ratio of long-term to short-term stress changes. The results of standard local bifurcation analysis indicate the onset of complex dynamics for very low values of time delay. On the other side, numerical approach confirms an additional complexity that was not observed by local analysis, due to the possible effect of global bifurcations. The most complex dynamics is detected for β < 1, with a complete Ruelle–Takens–Newhouse route to chaos under the variation of T_d, or the co-action of both parameters T_d and μ. These results correspond well with the previous experimental observations on clay and siltstone with low clay fraction. In the same regime, the perturbation of only a single parameter, μ, renders the oscillatory motion of the block. Within the velocity-independent regime, β = 1, the inclusion and variation of T_d generates a transition to equilibrium state, whereas the small oscillations of μ induce oscillatory motion with decreasing amplitude. The co-action of both parameters, in the same regime, causes the decrease of block’s velocity. As for β > 1, highly-frequent, limit-amplitude oscillations of initial stress give rise to oscillatory motion. Also for β > 1, in case of perturbing only the initial shear stress, with smaller amplitude, velocity of the block changes exponentially fast. If the time delay is introduced, besides the stress perturbation, within the same regime, the co-action of T_d (T_d < 0.1) and small oscillations of μ induce the onset of deterministic chaos.     

GNGA for general regions: Semilinear elliptic PDE and crossing eigenvalues

We consider the semilinear elliptic PDE Δu + f(λ, u) = 0 with the zero-Dirichlet boundary condition on a family of regions, namely stadions. Linear problems on such regions have been widely studied in the past. We seek to observe the corresponding phenomena in our nonlinear setting. Using the Gradient Newton Galerkin Algorithm (GNGA) of Neuberger and Swift, we document bifurcation, nodal structure, and symmetry of solutions. This paper provides the first published instance where the GNGA is applied to general regions. Our investigation involves both the dimension of the stadions and the value λ as parameters. We find that the so-called crossings and avoided crossings of eigenvalues as the dimension of the stadions vary influences the symmetry and variational structure of nonlinear solutions in a natural way.     

The stability of some classes of systems in standard Bogolyubov form

The stability of the equilibrium of quasilinear systems in standard Bogolyubov form is investigated. Classes of systems are distinguished out for which it is possible to determine the threshold value ?₀ of the small parameter ? which ensures qualitative agreement between the solutions of the initial systems of equations and the solutions of the averaged system corresponding it in an infinite time interval when ? < ?₀.     

Wave propagation in nonlocal elastic continua modelled by a fractional calculus approach

In the present paper, the wave propagation in one-dimensional elastic continua, characterized by nonlocal interactions modeled by fractional calculus, is investigated. Spatial derivatives of non-integer order 1 < α < 2 are involved in the governing equation, which is solved by fractional finite differences. The influence of long-range interactions is then analyzed as α varies: the resonant frequencies and the standing waves of a nonlocal bar are evaluated and the deviations from the classical (local) ones, recovered by imposing α = 2, are discussed.     

Coupled nonlinear waves in two-dimensional lattice

For one-dimensional nonlinear lattices, such as Toda lattice, it has been extensively studied. By considering the nonlinear effects of two-dimensional lattice, we set up the equation of motion for each particles (atoms, molecules or ions). For small amplitude and long wavelength nonlinear waves in this system, both the linear dispersion relation and the coupled Korteweg de Vries (KdV) equation are obtained. The simple soliton solution is obtained. If the nonlinear lattice is symmetric in the x and y directions, It is noted that there are two kinds of solitons. one is that propagates in either x or y directions, (1, 0) or (0, 1), the other is that propagates in the direction of (1, 1). It is in agreements with that of one-dimensional lattice. The different properties are investigated for different nonlinear interacting potentials, such as Toda potential, Morse potential and LJ potential.     

Numerical simulation of magnetohydrodynamic buoyancy-induced flow in a non-isothermally heated square enclosure

Numerical simulation of magnetohydrodynamic (MHD) buoyancy-induced heat transfer and fluid flow has been analyzed in a non-isothermally heated square enclosure using finite volume method. The bottom wall of enclosure were heated and cooled with a sinusoidal function and top wall was cooled isothermally. Vertical walls of the enclosure were adiabatic. Effects of Rayleigh number (Ra = 10⁴, 10⁵ and 10⁶), Hartman number (Ha = 0, 50 and 100) and amplitude of sinusoidal function (n = 0.25, 0.5 and 1) on temperature and flow fields were analyzed. It was observed that heat transfer was decreased with increasing Hartmann number and decreasing value of amplitude of sinusoidal function.     

Bifurcations of travelling wave solutions in the (N + 1)-dimensional sine–cosine-Gordon equations

In this paper, the (N + 1)-dimensional sine–cosine-Gordon equations are studied. The existence of solitary wave, kink and anti-kink wave, and periodic wave solutions are proved, by using the method of bifurcation theory of dynamical systems. All possible bounded exact explicit parametric representations of the above travelling solutions are obtained.     

Effect of hot ion temperature on obliquely propagating ion-acoustic solitons and double layers in an auroral plasma

Properties of obliquely propagating ion-acoustic solitons and double layers in a magnetized auroral plasma composed of hot adiabatic ions and two types of, cool and hot Maxwellian electrons are studied using Sagdeev pseudo-potential technique and assuming the quasi-neutrality condition. The new and surprising result which emerges from the model is that in contrast to the case of cold ions where ion-acoustic solitons and double layers are found for subsonic Mach numbers only, the hot ions case allows these nonlinear structures to exist for both subsonic and supersonic Mach number regimes. The double layers exist at lower angle of propagation as hot ion temperature is increased. The soliton electric field amplitudes are increased but their width and pulse duration are decreased with the increase in hot ion temperature. For the auroral zone parameters, the maximum electric field amplitude, width, pulse duration and speed for the solitons come out to be in the range ∼ (0.3–15) mV/m, ∼ (195–455) m, (7–20) ms and (22–26) km/s, respectively. The results seem to be in agreement with the Viking satellite observations in the auroral zone.     

Mathematical analysis of a two-patch model of tuberculosis disease with staged progression

The spread of tuberculosis is studied through a two-patch epidemiological system SE₁ ? E_nI which incorporates migrations from one patch to another just by susceptible individuals. Our model is consider with bilinear incidence and migration between two patches, where infected and infectious individuals cannot migrate from one patch to another, due to medical reasons. The existence and uniqueness of the associated endemic equilibria are discussed. Quadratic forms and Lyapunov functions are used to show that when the basic reproduction ratio is less than one, the disease-free equilibrium (DFE) is globally asymptotically stable, and when it is greater than one there exists in each case a unique endemic equilibrium (boundary equilibria and endemic equilibrium) which is globally asymptotically stable. Numerical simulation results are provided to illustrate the theoretical results.     

On generating (2 + 1)-dimensional hierarchies of evolution equations

Two isospectral problems are constructed with the help of a 6-dimensional Lie algebra. By using the Tu scheme, a (1 + 1)-dimensional expanding integrable couplings of the KdV hierarchy is obtained and the corresponding Hamiltonian structure is established. In addition, the 2-order matrix operators proposed by Tuguizhang are extended to the case where some 4-order matrices are given. Based on the extension, a new hierarchy of 2 + 1 dimensions is obtained by the Hamiltonian operator of the above (1 + 1)-dimensional case and the TAH scheme. The new hierarchy of 2 + 1 dimensions can be reduced to a coupled (2 + 1)-dimensional nonlinear equation and furthermore it can be reduced to the (2 + 1)-dimensional KdV equation which has important physics applications. The Hamiltonian structure for the (2 + 1)-dimensional hierarchy is derived with the aid of an extended trace identity. To the best of our knowledge, generating the (2 + 1)-dimensional equation hierarchies by virtue of the TAH scheme has not been studied in detail except to previous little work by Tu et al.     

Mathematical model for coupled roll and yaw motions of a floating body in regular waves under resonant and non-resonant conditions

The prediction of resonance is very important with respect to the vessels stability in the early stages of design. In this paper, an efficient modeling approach is presented to determine coupled roll and yaw motions of a symmetric and slender floating body when the influences of small amplitude regular waves are dominant. The angular motions described in time domain by considering all internal and external forces are transformed into frequency domain to obtain motion characteristics. We adopt a semi-analytical treatment to obtain roll and yaw motions and derive system instability due to roll resonance. To compute hydrodynamic forces, we employ strip theory method where frequency dependent sectional added-mass, damping and restoring coefficients are derived from the Frank’s close-fit curve. Numerical experiments carried out for a vessel of mass 19,190 ton under the action of wave of frequencies 0.56 and 0.76 rad/s with zero and non-zero initial conditions are reported and the effect of various parameters on system stability is investigated. Model results indicate that damping factor (ς) plays a pivotal role when wave encountering frequency (ω) and undamped natural frequency (β) are nearly equal. The essence of this study lies in the efficient modeling technique to evaluate damping factor and critical encountering frequency regime for a given ship particulars when experimentally derived resonance zone is absent.     

The new airfoil model NACA0015, modal analysis and flutter properties

The experimental airfoil model NACA0015 was used to study aeroelastic phenomena during self-excited profile vibration. It provides data for control of aeroelastic calculation programs at subsonic speeds of the stream. The model movability is two-dimensional with two-degree of freedom dynamic system, one in pitch and the second in plunge, and is proposed to be a dynamic system having two near corresponding eigenfrequencies. To quantitatively evaluate flow field using interferometry, a special test section design and profile was constructed. It utilized a large visual field for the optical system together with the option of changing support stiffness for both degrees of freedom. In this paper experimental results from the range of Reynolds numbers Re = (2.63–2.83) 10⁵ are published. The identified eigenvalues and eigenmodes for zero flow velocity are compared with measured flutter properties (frequency, modes and time evolutions) of the airfoil.     

Neural networks analysis of free laminar convection heat transfer in a partitioned enclosure

The feasibility of using neural networks (NNs) to predict the complete thermal and flow variables throughout a complicated domain, due to free convection, is demonstrated. Attention is focused on steady, laminar, two-dimensional, natural convective flow within a partitioned cavity. The objective is to use NN (trained on a database generated by a CFD analysis of the problem of a partitioned enclosure) to predict new cases; thus saving effort and computation time. Three types of NN are evaluated, namely General Regression NNs, Polynomial NNs, and a versatile design of Backpropagation neural networks. An important aspect of the study was optimizing network architecture in order to achieve best performance. For each of the three different NN architectures evaluated, parametric studies were performed to determine network parameters that best predict the flow variables.A CFD simulation software was used to generate a database that covered the range of Rayleigh number Ra = 10⁴–5 × 10⁶. The software was used to calculate the temperature, the pressure, and the horizontal and vertical components of flow speed. The results of the CFD were used for training and testing the neural networks (NN). The robustness of the trained NNs was tested by applying them to a “production” data set (1500 patterns for Ra = 8 × 10⁴ and 1500 patterns for Ra = 3 × 10⁶), which the networks have never been “seen” before. The results of applying the technique on the “production” data set show excellent prediction when the NNs are properly designed. The success of the NN in accurately predicting free convection in partitioned enclosures should help reduce analysis-time and effort. Neural networks could potentially help solve some cases in which CFD fails to solve because of numerical instability.     

Delayed transitions in non-linear replicator networks: About ghosts and hypercycles

In this paper we analyze delayed transition phenomena associated to extinction thresholds in a mean field model for hypercycles composed of three and four units, respectively. Hence, we extend a previous analysis carried out with the two-membered hypercycle [see Sardanyés J, Solé RV. Ghosts in the origins of life? Int J Bifurcation Chaos 2006;16(9), in press]. The models we analyze show that, after the tangent bifurcation, these hypercycles also leave a ghost in phase space. These ghosts, which actually conserve the dynamical properties of the coalesced coexistence fixed point, delay the flows before hypercycle extinction. In contrast with the two-component hypercycle, both ghosts show a plateau in the delay as ϕ → 0, thus displacing the power-law dependence to higher values of ϕ, in which the scaling law is now given by τ ∼ ϕ^β, with β = −1/3 (where τ is the delay and ϕ = ϵ − ϵ_c, the parametric distance above the extinction bifurcation point). These results suggest that the presence of the ghost is a general property of hypercycles. Such ghosts actually cause a memory effect which might increase hypercycle survival chances in fluctuating environments.     

Time-delayed control to suppress the nonlinear vibrations of a horizontally suspended Jeffcott-rotor system

Time-delay is an unavoidable phenomenon in active control systems. Measuring of the system states, processing of the measured signals, executing the control laws, conditioning and enforcing the control actions are the main reasons of time-delayed systems. This paper studies the vibration control of a horizontally suspended Jeffcott-rotor system having cubic and quadratic nonlinearities via time-delayed position-velocity controller. The intervals of the time-delays (τ₁ and τ₂) at which the system response is stable has been studied. The τ₁ − τ₂ plane is constructed to illustrate the area at which the system solutions are stable. The influences of the controller gains on the stable-solutions area in τ₁ − τ₂ plane are explored. The analysis revealed that the time-delay increases the vibration amplitudes and can destabilize the system solution in the case of negative position feedback control, while at positive position feedback control it improves the vibration suppression performance. The time-delays mechanism in stabilizing and destabilizing the dynamical systems is explained. Then, we proposed a simple and concrete method to determine the optimal value for time-delays that can improve the vibrations suppression efficiency. The acquired analytical results are confirmed numerically and the optimal working conditions of the system are concluded. Finally, a comparison with the papers that published previously is included.