Local Bifurcations in a Nonlinear Model of a Bioreactor

This paper is partially supported by the Bulgarian Science Fund under grant Nr. DO 02– 359/2008.

Analyse des bifurcations dans un modèle du flutter auriculaire

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Analyse des bifurcations dans un modèle du flutter auriculaire

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Big bang bifurcations in von Bertalanffy’s dynamics with strong and weak Allee effects

The main purpose of this work was to study population dynamic discrete models in which the growth of the population is described by generalized von Bertalanffy's functions, with an adjustment or correction factor of polynomial type. The consideration of this correction factor is made with the aim to introduce the Allee effect. To the class of generalized von Bertalanffy's functions is identified and characterized subclasses of strong and weak Allee's functions and functions with no Allee effect. This classification is founded on the concepts of strong and weak Allee's effects to population growth rates associated. A complete description of the dynamic behavior is given, where we provide necessary conditions for the occurrence of unconditional and essential extinction types. The bifurcation structures of the parameter plane are analyzed regarding the evolution of the Allee limit with the aim to understand how the transition from strong Allee effect to no Allee effect, passing through the weak Allee effect, is realized. To generalized von Bertalanffy's functions with strong and weak Allee effects is identified an Allee's effect region, to which is associated the concepts of chaotic semistability curve and Allee's bifurcation point. We verified that under some sufficient conditions, generalized von Bertalanffy's functions have a particular bifurcation structure: the big bang bifurcations of the so-called box-within-a-box type. To this family of maps, the Allee bifurcation points and the big bang bifurcation points are characterized by the symmetric of Allee's limit and by a null intrinsic growth rate. The present paper is also a significant contribution in the framework of the big bang bifurcation analysis for continuous 1D maps and unveil their relationship with the explosion birth and the extinction phenomena.

Stability analysis based on birfurcation theory of the DSTATCOM operating in current control mode

This paper presents the stability analysis for a distribution static compensator (DSTATCOM) that operates in current control mode based on bifurcation theory. Bifurcations delimit the operating zones of nonlinear circuits and, hence, the capability to compute these bifurcations is of important interest for practical design. A control design for the DSTATCOM is proposed. Along with this control, a suitable mathematical representation of the DSTATCOM is proposed to carry out the bifurcation analysis efficiently. The stability regions in the Thevenin equivalent plane are computed for different power factors at the point of common coupling. In addition, the stability regions in the control gain space, as well as the contour lines for different Floquet multipliers are computed. It is demonstrated through bifurcation analysis that the loss of stability in the DSTATCOM is due to the emergence of a Neimark bifurcation. The observations are verified through simulation studies.

Planar radial spots in a three-component FitzHugh-Nagumo system

Localized planar patterns arise in many reaction-diffusion models. Most of the paradigm equations that have been studied so far are two-component models. While stationary localized structures are often found to be stable in such systems, travelling patterns either do not exist or are found to be unstable. In contrast, numerical simulations indicate that localized travelling structures can be stable in three-component systems. As a first step towards explaining this phenomenon, a planar singularly perturbed three-component reaction-diffusion system that arises in the context of gas-discharge systems is analysed in this paper. Using geometric singular perturbation theory, the existence and stability regions of radially symmetric stationary spot solutions are delineated and, in particular, stable spots are shown to exist in appropriate parameter regimes. This result opens up the possibility of identifying and analysing drift and Hopf bifurcations, and their criticality, from the stationary spots described here.

Pulse dynamics in a three-component system : existence analysis

In this article, we analyze the three-component reaction-diffusion system originally developed by Schenk et al. (PRL 78:3781–3784, 1997). The system consists of bistable activator-inhibitor equations with an additional inhibitor that diffuses more rapidly than the standard inhibitor (or recovery variable). It has been used by several authors as a prototype three-component system that generates rich pulse dynamics and interactions, and this richness is the main motivation for the analysis we present. We demonstrate the existence of stationary one-pulse and two-pulse solutions, and travelling one-pulse solutions, on the real line, and we determine the parameter regimes in which they exist. Also, for one-pulse solutions, we analyze various bifurcations, including the saddle-node bifurcation in which they are created, as well as the bifurcation from a stationary to a travelling pulse, which we show can be either subcritical or supercritical. For two-pulse solutions, we show that the third component is essential, since the reduced bistable two-component system does not support them. We also analyze the saddle-node bifurcation in which two-pulse solutions are created. The analytical method used to construct all of these pulse solutions is geometric singular perturbation theory, which allows us to show that these solutions lie in the transverse intersections of invariant manifolds in the phase space of the associated six-dimensional travelling wave system. Finally, as we illustrate with numerical simulations, these solutions form the backbone of the rich pulse dynamics this system exhibits, including pulse replication, pulse annihilation, breathing pulses, and pulse scattering, among others.

Discrete phase-locked loop systems and spreadsheets

This paper demonstrates the use of a spreadsheet in exploring non-linear difference equations that describe digital control systems used in radio engineering, communication and computer architecture. These systems, being the focus of intensive studies of mathematicians and engineers over the last 40 years, may exhibit extremely complicated behaviour interpreted in contemporary terms as transition from global asymptotic stability to chaos through period-doubling bifurcations. The authors argue that embedding advanced mathematical ideas in the technological tool enables one to introduce fundamentals of discrete control systems in tertiary curricula without learners having to deal with complex machinery that rigorous mathematical methods of investigation require. In particular, in the appropriately designed spreadsheet environment, one can effectively visualize a qualitative difference in the behviour of systems with different types of non-linear characteristic.

Decision making in social neurobiological systems modeled as transitions in dynamic pattern formation

Extant models of decision making in social neurobiological systems have typically explained task dynamics as characterized by transitions between two attractors. In this paper, we model a three-attractor task exemplified in a team sport context. The model showed that an attacker–defender dyadic system can be described by the angle x between a vector connecting the participants and the try line. This variable was proposed as an order parameter of the system and could be dynamically expressed by integrating a potential function. Empirical evidence has revealed that this kind of system has three stable attractors, with a potential function of the form V(x)=−k1x+k2ax2/2−bx4/4+x6/6, where k1 and k2 are two control parameters. Random fluctuations were also observed in system behavior, modeled as white noise εt, leading to the motion equation dx/dt = −dV/dx+Q0.5εt, where Q is the noise variance. The model successfully mirrored the behavioral dynamics of agents in a social neurobiological system, exemplified by interactions of players in a team sport.

Infinite Dimensional Slow Modulations In A Delayed Model For Orthogonal Cutting Vibrations

We apply the method of multiple scales (MMS) to a well known model of regenerative cutting vibrations in the large delay regime. By ``large'' we mean the delay is much larger than the time scale of typical cutting tool oscillations. The MMS upto second order for such systems has been developed recently, and is applied here to study tool dynamics in the large delay regime. The second order analysis is found to be much more accurate than first order analysis. Numerical integration of the MMS slow flow is much faster than for the original equation, yet shows excellent accuracy. The main advantage of the present analysis is that infinite dimensional dynamics is retained in the slow flow, while the more usual center manifold reduction gives a planar phase space. Lower-dimensional dynamical features, such as Hopf bifurcations and families of periodic solutions, are also captured by the MMS. Finally, the strong sensitivity of the dynamics to small changes in parameter values is seen clearly.

Study of reproducibility of human arterial plaque reconstruction and its effects on stress analysis based on multispectral in vivo magnetic resonance imaging

Purpose: To quantify the uncertainties of carotid plaque morphology reconstruction based on patient-specific multispectral in vivo magnetic resonance imaging (MRI) and their impacts on the plaque stress analysis. Materials and Methods: In this study, three independent investigators were invited to reconstruct the carotid bifurcation with plaque based on MR images from two subjects to study the geometry reconstruction reproducibility. Finite element stress analyses were performed on the carotid bifurcations, as well as the models with artificially modified plaque geometries to mimic the image segmentation uncertainties, to study the impacts of the uncertainties to the stress prediction. Results: Plaque reconstruction reproducibility was generally high in the study. The uncertainties among interobservers are around one or the subpixel level. It also shows that the predicted stress is relatively less sensitive to the arterial wall segmentation uncertainties, and more affected by the accuracy of lipid region definition. For a model with lipid core region artificially increased by adding one pixel on the lipid region boundary, it will significantly increase the maximum Von Mises Stress in fibrous cap (>100%) compared with the baseline model for all subjects. Conclusion: The current in vivo MRI in the carotid plaque could provide useful and reliable information for plaque morphology. The accuracy of stress analysis based on plaque geometry is subject to MRI quality. The improved resolution/quality in plaque imaging with newly developed MRI protocols would generate more realistic stress predictions.

Nonlinear secondary whirl of an overhung rotor

Gravity critical speeds of rotors have hitherto been studied using linear analysis, and ascribed to rotor stiffness asymmetry. Here, we study an idealized asymmetric nonlinear overhung rotor model of Crandall and Brosens, spinning close to its gravity critical speed.Nonlinearities arise from finite displacements, and the rotor's staticlateral deflection under gravity is taken as small. Assuming small asymmetry and damping, slow modulations of whirl amplitudes are studied using the method of multiple scales. Inertia asymmetry appears only at second order. More interestingly, even without stiffness asymmetry, the gravity-induced resonance survives through geometric nonlinearities. The gravity resonant forcing does not influence the resonant mode at leading order, unlike the typical resonant oscillations. Nevertheless,the usual phenomena of resonances, namely saddle-node bifurcations, jump phenomena and hysteresis, are all observed. An unanticipated periodic solution branch is found. In the three-dimensional space oftwo modal coefficients and a detuning parameter, the full set of periodic solutions is found to be an imperfect version of three mutually intersecting curves: a straight line,a parabola and an ellipse.

Infinite dimensional slow modulations in a well known delayed model for cutting tool vibrations

We apply the method of multiple scales (MMS) to a well-known model of regenerative cutting vibrations in the large delay regime. By ``large'' we mean the delay is much larger than the timescale of typical cutting tool oscillations. The MMS up to second order, recently developed for such systems, is applied here to study tool dynamics in the large delay regime. The second order analysis is found to be much more accurate than the first order analysis. Numerical integration of the MMS slow flow is much faster than for the original equation, yet shows excellent accuracy in that plotted solutions of moderate amplitudes are visually near-indistinguishable. The advantages of the present analysis are that infinite dimensional dynamics is retained in the slow flow, while the more usual center manifold reduction gives a planar phase space; lower-dimensional dynamical features, such as Hopf bifurcations and families of periodic solutions, are also captured by the MMS; the strong sensitivity of the slow modulation dynamics to small changes in parameter values, peculiar to such systems with large delays, is seen clearly; and though certain parameters are treated as small (or, reciprocally, large), the analysis is not restricted to infinitesimal distances from the Hopf bifurcation.

Dynamics of a driven magneto-martensitic ribbon

We investigate the dynamics of a sinusoidally driven ferromagnetic martensitic ribbon by adopting a recently introduced model that involves strain and magnetization as order parameters. Retaining only the dominant mode of excitation we reduce the coupled set of partial differential equations for strain and magnetization to a set of coupled ordinary nonlinear equations for the strain and magnetization amplitudes. The equation for the strain amplitude takes the form of parametrically driven oscillator. Finite strain amplitude can only be induced beyond a critical value of the strength of the magnetic field. Chaotic response is seen for a range of values of all the physically interesting parameters. The nature of the bifurcations depends on the choice of temperature relative to the ordering of the Curie and the martensite transformation temperatures. We have studied the nature of response as a function of the strength and frequency of the magnetic field, and magneto-elastic coupling. In general, the bifurcation diagrams with respect to these parameters do not follow any standard route. The rich dynamics exhibited by the model is further illustrated by the presence of mixed mode oscillations seen for low frequencies. The geometric structure of the mixed mode oscillations in the phase space has an unusual deep crater structure with an outer and inner cone on which the orbits circulate. We suggest that these features should be seen in experiments on driven magneto-martensitic ribbons. (C) 2014 Elsevier B. V. All rights reserved.

Time scales of the stick-slip dynamics of the peeling of an adhesive tape

The stick-slip dynamics of the peeling of an adhesive tape is characterized by bifurcations that have been experimentally well studied. In this work, we investigate the time scale in which the the stick-slips happen leading to the bifurcations. This is fundamental to understanding the triboluminescence and acoustic emissions associated with the bifurcations. We establish a relationship between the time scale of the bifurcations and the inherent mathematical structure of the peeling dynamics by studying a characteristic time quantity associated with the dynamics.