Mixed-mode oscillations in a multiple time scale phantom bursting system

In this work we study mixed mode oscillations in a model of secretion of GnRH (Gonadotropin Releasing Hormone). The model is a phantom burster consisting of two feedforward coupled FitzHugh-Nagumo systems, with three time scales. The forcing system (Regulator) evolves on the slowest scale and acts by moving the slow nullcline of the forced system (Secretor). There are three modes of dynamics: pulsatility (transient relaxation oscillation), surge (quasi steady state) and small oscillations related to the passage of the slow nullcline through a fold point of the fast nullcline. We derive a variety of reductions, taking advantage of the mentioned features of the system. We obtain two results; one on the local dynamics near the fold in the parameter regime corresponding to the presence of small oscillations and the other on the global dynamics, more specifically on the existence of an attracting limit cycle. Our local result is a rigorous characterization of small canards and sectors of rotation in the case of folded node with an additional time scale, a feature allowing for a clear geometric argument. The global result gives the existence of an attracting unique limit cycle, which, in some parameter regimes, remains attracting and unique even during passages through a canard explosion.Comment: 38 pages, 16 figure

Dynamic modeling and stability analysis of a nonlinear system with primary resonance

In recent years, there has been growing interest in the study of nonlinear phenomena. This is due to the modernization of structures related to the need of using lighter, more resistant and flexible materials. Thus, this work aims to study the behavior of a mechanical system with two degrees of freedom with nonlinear characteristics in primary resonance. The structure consists of the main system connected to a secondary system to act as a Nonlinear Dynamic Vibration Absorber, which partially or fully absorbs the vibrational energy of the system. The numerical solutions of the problem are obtained using the Runge-Kutta methods of the 4th order and approximate analytical solutions are obtained using the Multiple Scales Method. Then, the approximation error between the two solutions is analyzed. Using the aforementioned perturbation method, the responses for the ordinary differential equations of the first order can be determined, which describe the modulation amplitudes and phases. Thus, the solution in steady state and the stability are studied using the frequency response. Furthermore, the behavior of the main system and the absorber are investigated through numerical simulations, such as responses in the time domain, phase planes and Poincaré map; which shows that the system displays periodic, quasi-periodic and chaotic movements. The dynamic behavior of the system is analyzed using the Lyapunov exponent and the bifurcation diagram is presented to better summarize all the possible behaviors as the force amplitude varies. In general, the main characteristics of a dynamic system that experiences the chaotic response will be identified

Recovering ‘lost’ information in the presence of noise: application to rodent–predator dynamics.

A Hamiltonian approach is introduced for the reconstruction of trajectories and models of complex stochastic dynamics from noisy measurements. The method converges even when entire trajectory components are unobservable and the parameters are unknown. It is applied to reconstruct nonlinear models of rodent–predator oscillations in Finnish Lapland and high-Arctic tundra. The projected character of noisy incomplete measurements is revealed and shown to result in a degeneracy of the likelihood function within certain null-spaces. The performance of the method is compared with that of the conventional Markov chain Monte Carlo (MCMC) technique

Global bifurcations to subcritical magnetorotational dynamo action in Keplerian shear flow

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, non-linear magnetohydrodynamic process whose study is relevant to the understanding of accretion processes and magnetic field generation in astrophysics. Transition to this form of dynamo action is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we build on recent work on the two problems to investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are found to be supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. The results suggest that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the physical conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in a variety of astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows

Frequency domain and time domain analysis of thermoacoustic oscillations with wave-based acoustics

Many thermoacoustic systems exhibit rich nonlinear behaviour. Recent studies show that this nonlinear dynamics can be well captured by low-order time domain models that couple a level set kinematic model for a laminar flame, the $G$-equation, with a state-space realization of the linearized acoustic equations. However, so far the $G$-equation has been coupled only with straight ducts with uniform mean acoustic properties, which is a simplistic configuration. In this study, we incorporate a wave-based model of the acoustic network, containing area and temperature variations and frequency-dependent boundary conditions. We cast the linear acoustics into state-space form using a different approach from that in the existing literature. We then use this state-space form to investigate the stability of the thermoacoustic system, both in the frequency and time domains, using the flame position as a control parameter. We observe frequency-locked, quasiperiodic and chaotic oscillations. We identify the location of Neimark–Sacker bifurcations with Floquet theory. We also find the Ruelle–Takens–Newhouse route to chaos with nonlinear time series analysis techniques. We highlight important differences between the nonlinear response predicted by the frequency domain and the time domain methods. This reveals deficiencies with the frequency domain technique, which is commonly used in academic and industrial studies of thermoacoustic systems. We then demonstrate a more accurate approach based on continuation analysis applied to time domain techniques.This project was funded by the European Research Council through Project ALORS 2590620.This is the author accepted manuscript. The final version is available from Cambridge University Press via http://dx.doi.org/10.1017/jfm.2015.13