Bifurcation analysis of steady states and limit cycles in a thermal pulse combustor model

Oscillatory behaviour of state variables is desirable in pulse combustors, as properly designed pulsations lead to improved performances, such as higher thermal efficiency and lower emissions compared to steady combustors. In the present work, we perform a systematic investigation of the stability of steady states and limit cycles of a pulse combustor model through numerical continuation. Different bifurcation parameters such as tailpipe friction factor, wall temperature, convective heat transfer coefficient, inlet temperature and inlet fuel mass fraction are varied to identify the complete ranges of those parameters at which limit cycles can be expected. This analysis identifies alternative stable periodic regimes in parameter space (e.g. friction factor). In addition, a few performance indicators such as amplitude of oscillations, cycle-averaged heat transfer and cycle-averaged specific thrust are compared between different ranges of friction factor yielding limit cycle oscillations. The comparison clearly shows that, depending upon the application, friction factor can be chosen from different regimes. The time-integration of the model is also performed to support the bifurcation results obtained from numerical continuation, wherever appropriate. The complete stability margin of limit cycle oscillations for those bifurcation parameters can be useful for improved design of the combustor and for determining the optimal operating conditions of pulse combustors.     

From simple to complex oscillatory behavior in metabolic and genetic control networks

We present an overview of mechanisms responsible for simple or complex oscillatory behavior in metabolic and genetic control networks. Besides simple periodic behavior corresponding to the evolution toward a limit cycle we consider complex modes of oscillatory behavior such as complex periodic oscillations of the bursting type and chaos. Multiple attractors are also discussed, e.g., the coexistence between a stable steady state and a stable limit cycle (hard excitation), or the coexistence between two simultaneously stable limit cycles (birhythmicity). We discuss mechanisms responsible for the transition from simple to complex oscillatory behavior by means of a number of models serving as selected examples. The models were originally proposed to account for simple periodic oscillations observed experimentally at the cellular level in a variety of biological systems. In a second stage, these models were modified to allow for complex oscillatory phenomena such as bursting, birhythmicity, or chaos. We consider successively (1) models based on enzyme regulation, proposed for glycolytic oscillations and for the control of successive phases of the cell cycle, respectively; (2) a model for intracellular Ca(2+) oscillations based on transport regulation; (3) a model for oscillations of cyclic AMP based on receptor desensitization in Dictyostelium cells; and (4) a model based on genetic regulation for circadian rhythms in Drosophila. Two main classes of mechanism leading from simple to complex oscillatory behavior are identified, namely (i) the interplay between two endogenous oscillatory mechanisms, which can take multiple forms, overt or more subtle, depending on whether the two oscillators each involve their own regulatory feedback loop or share a common feedback loop while differing by some related process, and (ii) self-modulation of the oscillator through feedback from the system's output on one of the parameters controlling oscillatory behavior. However, the latter mechanism may also be viewed as involving the interplay between two feedback processes, each of which might be capable of producing oscillations. Although our discussion primarily focuses on the case of autonomous oscillatory behavior, we also consider the case of nonautonomous complex oscillations in a model for circadian oscillations subjected to periodic forcing by a light-dark cycle and show that the occurrence of entrainment versus chaos in these conditions markedly depends on the wave form of periodic forcing. (c) 2001 American Institute of Physics.     

Bursting Ca2+ Oscillations and Synchronization in Coupled Cells

A mathematical model proposed by Grubelnk et al. [Biophys. Chem. 94 （2001） 59] is employed to study the physiological role of mitochondria and the cytosolic proteins in generating complex Ca^2＋ oscillations, lntracellulax bursting calcium oscillations of point-point, point-cycle and two-folded limit cycle types are observed and explanations are given based on the fast/slow dynamical analysis, especially for point-cycle and two-folded limit cycle types, which have not been reported before. Furthermore, synchronization of coupled bursters of Ca^2＋ oscillations via gap junctions and the effect of bursting types on synchronization of coupled cells are studied. It is argued that bursting oscillations of point-point type may be superior to achieve synchronization than that of point cycle type.     

Phase synchronization between collective rhythms of globally coupled oscillator groups: noisy identical case

We theoretically investigate the collective phase synchronization between interacting groups of globally coupled noisy identical phase oscillators exhibiting macroscopic rhythms. Using the phase reduction method, we derive coupled collective phase equations describing the macroscopic rhythms of the groups from microscopic Langevin phase equations of the individual oscillators via nonlinear Fokker-Planck equations. For sinusoidal microscopic coupling, we determine the type of the collective phase coupling function, i.e., whether the groups exhibit in-phase or antiphase synchronization. We show that the macroscopic rhythms can exhibit effective antiphase synchronization even if the microscopic phase coupling between the groups is in-phase, and vice versa. Moreover, near the onset of collective oscillations, we analytically obtain the collective phase coupling function using center-manifold and phase reductions of the nonlinear Fokker-Planck equations.     

Coherent periodic activity in excitatory Erdös-Renyi neural networks: the role of network connectivity

In this article, we investigate the role of connectivity in promoting coherent activity in excitatory neural networks. In particular, we would like to understand if the onset of collective oscillations can be related to a minimal average connectivity and how this critical connectivity depends on the number of neurons in the networks. For these purposes, we consider an excitatory random network of leaky integrate-and-fire pulse coupled neurons. The neurons are connected as in a directed Erdo?s-Renyi graph with average connectivity scaling as a power law with the number of neurons in the network. The scaling is controlled by a parameter γ, which allows to pass from massively connected to sparse networks and therefore to modify the topology of the system. At a macroscopic level, we observe two distinct dynamical phases: an asynchronous state corresponding to a desynchronized dynamics of the neurons and a regime of partial synchronization (PS) associated with a coherent periodic activity of the network. At low connectivity, the system is in an asynchronous state, while PS emerges above a certain critical average connectivity (c). For sufficiently large networks, (c) saturates to a constant value suggesting that a minimal average connectivity is sufficient to observe coherent activity in systems of any size irrespectively of the kind of considered network: sparse or massively connected. However, this value depends on the nature of the synapses: reliable or unreliable. For unreliable synapses, the critical value required to observe the onset of macroscopic behaviors is noticeably smaller than for reliable synaptic transmission. Due to the disorder present in the system, for finite number of neurons we have inhomogeneities in the neuronal behaviors, inducing a weak form of chaos, which vanishes in the thermodynamic limit. In such a limit, the disordered systems exhibit regular (non chaotic) dynamics and their properties correspond to that of a homogeneous fully connected network for any γ-value. Apart for the peculiar exception of sparse networks, which remain intrinsically inhomogeneous at any system size.     

Bifurcation analysis of the travelling waveform of FitzHugh-Nagumo nerve conduction model equation

The FitzHugh-Nagumo model for travelling wave type neuron excitation is studied in detail. Carrying out a linear stability analysis near the equilibrium point, we bring out various interesting bifurcations which the system admits when a specific Z(2) symmetry is present and when it is not. Based on a center manifold reduction and normal form analysis, the Hopf normal form is deduced. The condition for the onset of limit cycle oscillations is found to agree well with the numerical results. We further demonstrate numerically that the system admits a period doubling route to chaos both in the presence as well as in the absence of constant external stimuli. (c) 1997 American Institute of Physics.     

Contagion effects in a chartist-fundamentalist model with time delays

In this paper two models of speculative markets are developed to study the effects of feedback mechanisms in financial markets. In the first model, a crash market model couples a linear chartist-fundamentalist model with time delays with a log-periodic market index I(t) through direct coupling. Numerical solutions to the model show that asset prices exhibit significant persistence as a result of the coupling to the log-periodic market index. An extension to include endogenous wealth dynamics shows that the chartists benefit from the persistent dynamics induced by the coupling. The second model is a two-asset model represented by a 2-dimensional delay-differential equation. Asset one price exhibits limit cycle dynamics while in the second market asset prices follow stable damped oscillations. The markets are coupled through a diffusive coupling term. Solutions to the coupled model show that the dynamics of asset two changes fundamentally with the price now exhibiting a limit cycle. The stable converging dynamics is replaced with limit cycle oscillations around the fundamental.     

Collective phase sensitivity

The collective phase response to a macroscopic external perturbation of a population of interacting nonlinear elements exhibiting collective oscillations is formulated for the case of globally coupled oscillators. The macroscopic phase sensitivity is derived from the microscopic phase sensitivity of the constituent oscillators by a two-step phase reduction. We apply this result to quantify the stability of the macroscopic common-noise-induced synchronization of two uncoupled populations of oscillators undergoing coherent collective oscillations.     

Mode-to-mode energy transfers in convective patterns

We investigate the energy transfer between various Fourier modes in a low-dimensional model for thermal convection. We have used the formalism of mode-to-mode energy transfer rate in our calculation. The evolution equations derived using this scheme is the same as those derived using the hydrodynamical equations for thermal convection in Boussinesq fluids. Numerical and analytical studies of this model show that convective rolls appear as the Rayleigh number R is raised above its critical value R _c. Further increase of Rayleigh number generates rolls in the perpendicular directions as well, and we obtain a dynamic asymmetric square pattern. This pattern is due to Hopf bifurcation. There are two sets of limit cycles corresponding to the two competing asymmetric square patterns. When the Rayleigh number is increased further, the limit cycles become unstable simultaneously, and chaotic motion sets in. The onset of chaos is via intermittent route. The trajectories wander for quite a long time almost periodically before jumping irregularly to one of the two ghost limit cycles.     

AEROELASTIC OSCILLATIONS OF A SEESAW-TYPE OSCILLATOR UNDER STRONG WIND CONDITIONS

This paper is concerned with one-degree-of-freedom aeroelastic oscillations of a seesaw- type structure in a steady wind flow. Here it is assumed that strong wind conditions induce nonlinear aeroelastic stiffness forces that are of the same order of magnitude as the structural stiffness forces. As a model equation for the aeroelastic behaviour of the seesaw-type structure, a strongly nonlinear self-excited oscillator is obtained. The bifurcation and the stability of limit cycles for this equation are studied using a special perturbation method. Both the case with linear structural stiffness and the case with nonlinear structural stiffness are studied. For both cases is assumed a general cubic approximation to describe the aerodynamic coefficient. Conditions for the existence, the stability, and the bifurcation of limit cycles are given.     

Coupling-induced oscillations in nonhomogeneous, overdamped, bistable systems

Coupling-induced oscillations in a homogeneous network of overdamped bistable systems have been previously studied both theoretically and experimentally for a system of N (odd) elements, unidirectionally coupled in a ring topology. In this work, we extend the analysis of this system to include a network of nonhomogeneous elements with respect to the parameter that controls the topology of the potential function and the bistability of each element. In particular, we quantify the effects of the nonhomogeneity on the onset of oscillations and the response of the network to external (assumed to be constant and very small) perturbations, using our (recently developed) coupled core fluxgate magnetometer as a representative system. The potential applications of this work include signal detection and characterization for a large class of sensor systems.     

Nonlinear Dynamics of Exciton-Polariton Condensates in a One-Dimensional Lattice

We consider the problem of exciton-polariton-condensate formation in a semiconductor microcavity in the strong coupling regime. The condensate is confined in a one-dimensional periodic potential and coupled to an exciton reservoir that is formed by the external cw pumping. The condensate dynamics is studied in the center and in the edges of the Brillouin zone (BZ). Modeling the formation of the condensate from weak initial noise shows that besides steady states the macroscopic oscillations of polariton density can also occur. Within the framework of the mean field approach, we obtain important analytic relations for the condensate eigenstates using the developed simplified model for three coupled spatial harmonics. A numerical analysis verifies the correctness of the analytic results.     

On intrinsic oscillation in radiation-affected diffusion flames

A linear stability analysis is conducted to study the onset of near-limit flame oscillation with radiative heat loss in 1-D chambered planar flames using multi-scale activation-energy asymptotics. The oscillatory instability near the radiation-induced extinction limit at large Damköhler numbers is identified, in additional to the one near the kinetic limit at small Damköhler numbers. It is shown that radiative loss assumes a similar role as varying the thermal diffusivity of the reactants. Thus, flame oscillation near the radiative limit is still thermal-diffusive in nature although it may develop under unity Lewis numbers. The unstable range of Damköhler numbers near the radiative limit shows quite similar parametric dependence on the Lewis numbers of reactants, Le_F and Le_O, the stoichiometry, ?, and the radiative loss as that near the kinetic limit. They both increase monotonically with Le_O and ? and increase then decrease with Le_F. Increasing radiative loss extends the parameter range under which flame oscillations may develop. However, they show different dependence on the temperature difference between the supplying reactants. Unless radiative loss approaches its maximum value the system can sustain, flame oscillation near the radiative limit is only possible within a limited range of ΔT, whereas it is promoted monotonically with decreasing ΔT near the kinetic limit. Furthermore, while radiative loss shows small effect on the nondimensional oscillation frequency, the dimensional frequency of flame oscillations near the radiative limit can be substantially smaller than that near the kinetic limit.     

On the dimension of deterministic and random Cantor-like sets,symbolic dynamics,and the Eckmann-Ruelle Conjecture

In this paper we unify and extend many of the known results on the dimension of deterministic and random Cantor-like sets in ? ⁿ , and apply these results to study some problems in dynamical systems. In particular, we verify the Eckmann-Ruelle Conjecture for equilibrium measures for Hölder continuous conformal expanding maps and conformal Axiom A^# (topologically hyperbolic) homeomorphims. We also construct a Hölder continuous Axiom A^# homeomorphism of positive topological entropy for which the unique measure of maximal entropy is ergodic and has different upper and lower pointwise dimensions almost everywhere. this example shows that the non-conformal Hölder continuous version of the Eckmann-Ruelle Conjecture is false. The Cantor-like sets we consider are defined by geometric constructions of different types. The vast majority of geometric constructions studied in the literature are generated by a finite collection ofp maps which are either contractions or similarities and are modeled by the full shift onp symbols (or at most a subshift of finite type). In this paper we consider much more general classes of geometric constructions: the placement of the basic sets at each step of the construction can be arbitrary, and they need not be disjoint. Moreover, our constructions are modeled by arbitrary symbolic dynamical systems. The importance of this is to reveal the close and nontrivial relations between the statistical mechanics (and especially the absence of phase transitions) of the symbolic dynamical system underlying the geometric construction and the dimension of its limit set. This has not been previously observed since no phase transitions can occur for subshifts of finite type. We also consider nonstationary constructions, random constructions (determined by an arbitrary ergodic stationary distribution), and combinations of the above.     

Self-sustained oscillations in a low-current discharge with a semiconductor serving as a cathode and ballast resistor: II. Theory

Current self-sustained oscillations in a Townsend discharge are studied analytically. The proposed method for solving simple equations is applicable even when the Townsend coefficient of electron reproduction g (the main factor governing the oscillations) cannot be determined theoretically. The coefficient g is related to the discharge current-voltage characteristic, which can be obtained experimentally or from numerical simulations. Self-oscillating solutions (limit cycles) are found under various conditions. The mechanisms governing the excitation and stabilization of the solutions obtained are interpreted. It is shown that the waveform of the undamped oscillations may change significantly (the current peaks are smoothed, and the oscillation period decreases) when a weak constant cathode current, which is presumably related to the emission caused by slowly diffusing metastable molecules, is included in the equations.     

Energy law preserving C finite element schemes for phase field models in two-phase flow computations

We use the idea in [33] to develop the energy law preserving method and compute the diffusive interface (phase-field) models of Allen–Cahn and Cahn–Hilliard type, respectively, governing the motion of two-phase incompressible flows. We discretize these two models using a C⁰ finite element in space and a modified midpoint scheme in time. To increase the stability in the pressure variable we treat the divergence free condition by a penalty formulation, under which the discrete energy law can still be derived for these diffusive interface models. Through an example we demonstrate that the energy law preserving method is beneficial for computing these multi-phase flow models. We also demonstrate that when applying the energy law preserving method to the model of Cahn–Hilliard type, un-physical interfacial oscillations may occur. We examine the source of such oscillations and a remedy is presented to eliminate the oscillations. A few two-phase incompressible flow examples are computed to show the good performance of our method.     

Critical clusters in interdependent economic sectors

In this paper we develop a data-driven hierarchical clustering methodology to group the economic sectors of a country in order to highlight strongly coupled groups that are weakly coupled with other groups. Specifically, we consider an input-output representation of the coupling among the sectors and we interpret the relation among sectors as a directed graph; then we recursively apply the spectral clustering methodology over the graph, without a priori information on the number of groups that have to be obtained. In order to do this, we resort to the eigengap criterion, where a suitable number of groups is selected automatically based on the intensity and structure of the coupling among the sectors. We validate the proposed methodology considering a case study for Italy, inspecting how the coupling among clusters and sectors changes from the year 1995 to 2011, showing that in the years the Italian structure underwent deep changes, becoming more and more interdependent, i.e., a large part of the economy has become tightly coupled.