Effects of constraint curvature on structural instability: tensile buckling and multiple bifurcations

Bifurcation of an elastic structure crucially depends on the curvature of the constraints against which the ends of the structure are prescribed to move, an effect which deserves more attention than it has received so far. In fact, we show theoretically and we provide definitive experimental verification that an appropriate curvature of the constraint over which the end of a structure has to slide strongly affects buckling loads and can induce: (i.) tensile buckling; (ii.) decreasing- (softening), increasing- (hardening), or constant-load (null stiffness) postcritical behaviour; (iii.) multiple bifurcations, determining for instance two bifurcation loads (one tensile and one compressive) in a single-degree-of-freedom elastic system. We show how to design a constraint profile to obtain a desired postcritical behaviour and we provide the solution for the elastica constrained to slide along a circle on one end, representing the first example of an inflexional elastica developed from a buckling in tension. These results have important practical implications in the design of compliant mechanisms and may find applications in devices operating in quasi-static or dynamic conditions.

Measurements of triggering and transient growth in a model lean-premixed gas turbine combustor

Instability triggering and transient growth of thermoacoustic oscillations were experimentally investigated in combination with linear/nonlinear flame transfer function (FTF) methodology in a model lean-premixed gas turbine combustor operated with CH 4 and air at atmospheric pressure. A fully premixed flame with 10kW thermal power and an equivalence ratio of 0.60 was chosen for detailed characterization of the nonlinear transient behaviors. Flame transfer functions were experimentally determined by simultaneous measurements of inlet velocity fluctuations and heat release rate oscillations using a constant temperature anemometer and OH */CH * chemiluminescence emissions, respectively. The phase-resolved variation of the local flame structure at a limit cycle was measured by planar laser-induced fluorescence of OH. Simultaneous measurements of inlet velocity, OH */CH * emission, and acoustic pressure were performed to investigate the temporal evolution of the system from a stable to a limit cycle operation. This measurement allows us to describe an unsteady instability triggering event in terms of several distinct stages: (i) initiation of a small perturbation, (ii) exponential amplification, (iii) saturation, (iv) nonlinear evolution of the perturbations towards a new unstable periodic state, (v) quasi-steady low-amplitude periodic oscillation, and (vi) fully-developed high-amplitude limit cycle oscillation. Phase-plane portraits of instantaneous inlet velocity and heat release rate clearly show the presence of two different attractors. Depending on its initial position in phase space at infinitesimally small amplitude, the system evolves towards either a high-amplitude oscillatory state or a low-amplitude oscillatory state. This transient phenomenon was analyzed using frequency- and amplitude-dependent damping mechanisms, and compared to subcritical and supercritical bifurcation theories. The results presented in this paper experimentally demonstrate the hypothesis proposed by Preetham et al. based on analytical and computational solutions of the nonlinear G-equation [J. Propul. Power 24 (2008) 1390-1402]. Good quantitative agreement was obtained between measurements and predictions in terms of the conditions for the onset of triggering and the amplitude of triggered combustion instabilities. © 2011 The Combustion Institute.

Progress towards the design and numerical analysis of a 3D microchannel biochip separator

This paper reports the design and numerical analysis of a three-dimensional biochip plasma blood separator using computational fluid dynamics techniques. Based on the initial configuration of a two-dimensional (2D) separator, five three-dimensional (3D) microchannel biochip designs are categorically developed through axial and plenary symmetrical expansions. These include the geometric variations of three types of the branch side channels (circular, rectangular, disc) and two types of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming that steady-state flow has been established, the behaviour of the blood fluid in the devices is algebraically analysed and numerically modelled. The roles of the relevant microchannel mechanisms, i.e. bifurcation, constriction and bending channel, on promoting the separation process are analysed based on modelling results. The differences among the different 3D implementations are compared and discussed. The advantages of 3D over 2D separator in increasing separation volume and effectively depleting cell-free layer fluid from the whole cross section circumference are addressed and illustrated. © 2011 John Wiley & Sons, Ltd.

Biofluid behaviour in 3D microchannel systems: Numerical analysis and design development of 3D microchannel biochip separators

This paper describes the design and development cycle of a 3D biochip separator and the modelling analysis of flow behaviour in the biochip microchannel features. The focus is on identifying the difference between 2D and 3D implementations as well as developing basic forms of 3D microfluidic separators. Five variants, based around the device are proposed and analysed. These include three variations of the branch channels (circular, rectangular, disc) and two variations of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming steady state flow has been established, the efficiencies of the flow between the main and side channels for the different designs are analysed and compared with regard to relevant biomicrofluidic laws or effects (bifurcation law, Fahraeus effect, cell-free phenomenon, bending channel effect and laminar flow behaviour). The modelling results identify flow features in microchannels, a constriction and bifurcations and show detailed differences in flow fields between the various designs. The manufacturing process using injection moulding for the initial base case design is also presented and discussed. The work reported here is supported as part of the UK funded 3D-MINTEGRATION project. © 2010 IEEE.

A Balance Equation Determines a Switch in Neuronal Excitability

We use the qualitative insight of a planar neuronal phase portrait to detect an excitability switch in arbitrary conductance-based models from a simple mathematical condition. The condition expresses a balance between ion channels that provide a negative feedback at resting potential (restorative channels) and those that provide a positive feedback at resting potential (regenerative channels). Geometrically, the condition imposes a transcritical bifurcation that rules the switch of excitability through the variation of a single physiological parameter. Our analysis of six different published conductance based models always finds the transcritical bifurcation and the associated switch in excitability, which suggests that the mathematical predictions have a physiological relevance and that a same regulatory mechanism is potentially involved in the excitability and signaling of many neurons. © 2013 Franci et al.

An organizing center in a planar model of neuronal excitability

This paper studies the excitability properties of a generalized FitzHugh-Nagumo model. The model differs from the classical FitzHugh-Nagumo model in that it accounts for the effect of cooperative gating variables such as activation of calcium currents. Excitability is explored by unfolding a pitchfork bifurcation that is shown to organize five different types of excitability. In addition to the three classical types of neuronal excitability, two novel types are described and distinctly associated to the presence of cooperative variables. © 2012 Society for Industrial and Applied Mathematics.

Delayed control of a Moore-Greitzer axial compressor model

Several feedback control laws have appeared in the literature concerning the stabilization of the nonlinear Moore-Greitzer axial compression model. Motivated by magnitude and rate limitations imposed by the physical implementation of the control law, Larsen et al. studied a dynamic implementation of the S-controller suggested by Sepulchre and Kokotović. They showed the potential benefit of implementing the S-controller through a first-order lag: while the location of the closed-loop equilibrium achieved with the static control law was sensitive to poorly known parameters, the dynamic implementation resulted in a small limit cycle at a very desirable location, insensitive to parameter variations. In this paper, we investigate the more general case when the control is applied with a time delay. This can be seen as an extension of the model with a first-order lag. The delay can either be a result of system constraints or be deliberately implemented to achieve better system behavior. The resulting closed-loop system is a set of parameter-dependent delay differential equations. Numerical bifurcation analysis is used to study this model and investigate whether the positive results obtained for the first-order model persist, even for larger values of the delay.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

Lock-in and quasiperiodicity in a forced hydrodynamically self-excited jet

The ability of hydrodynamically self-excited jets to lock into strong external forcing is well known. Their dynamics before lock-in and the specific bifurcations through which they lock in, however, are less well known. In this experimental study, we acoustically force a low-density jet around its natural global frequency. We examine its response leading up to lock-in and compare this to that of a forced van der Pol oscillator. We find that, when forced at increasing amplitudes, the jet undergoes a sequence of two nonlinear transitions: (i) from periodicity to T{double-struck}2 quasiperiodicity via a torus-birth bifurcation; and then (ii) from T{double-struck}2 quasiperiodicity to 1:1 lock-in via either a saddle-node bifurcation with frequency pulling, if the forcing and natural frequencies are close together, or a torus-death bifurcation without frequency pulling, but with a gradual suppression of the natural mode, if the two frequencies are far apart. We also find that the jet locks in most readily when forced close to its natural frequency, but that the details contain two asymmetries: the jet (i) locks in more readily and (ii) oscillates more strongly when it is forced below its natural frequency than when it is forced above it. Except for the second asymmetry, all of these transitions, bifurcations and dynamics are accurately reproduced by the forced van der Pol oscillator. This shows that this complex (infinite-dimensional) forced self-excited jet can be modelled reasonably well as a simple (three-dimensional) forced self-excited oscillator. This result adds to the growing evidence that open self-excited flows behave essentially like low-dimensional nonlinear dynamical systems. It also strengthens the universality of such flows, raising the possibility that more of them, including some industrially relevant flames, can be similarly modelled. © 2013 Cambridge University Press.

A new species of Allocreadium (Trematoda : Allocreadiidae) from freshwater fishes in the Danjiangkou Reservoir in China

A new species of Allocreadium, Allocreadium danjiangensis n. sp., is described from the intestine of several species of freshwater fish, including Abbottina rivularis (Basilewsky, 1855), Sarcocheilichthys nigripinnis nigripinns (Gunther, 1873), Gnathopogon argentatus (Sauvage et Dabry 1874), Opsariichthys uncirostris bidens (Gunther, 1873), and Erythroculter mongolicus mongolicus (Basilewsky, 1855) (Cyprinidae) from the Danjiangkou Reservoir in central China. The main morphological characters of the new species are as follows: vitelline follicles numerous, extending from the level of acetabulum to posterior extremity, distributed over both sides around the ceca; cirrus sac relatively large, developed, lying obliquely anterior to the acetabulum, extending from the level of the intestinal bifurcation to the central level of acetabulum, and overlapping left or right cecal; and ovary much smaller than testes, generally close to or even overlapping the anterior border of anterior testis. Observation by scanning electron microscopy shows only 2 kinds of tegumental formations, i.e., papillae and tubercles, instead of 3 types of tegumental formations, i.e., papillae, bosses, and minute sensor receptors observed on other species of the Allocreadiidae. The tegumental striations of the present species vary on the different parts of the body. In addition, a new structure, identified as the "groove" with a tonguelike tubercle, was observed on the inner wall of acetabulum.

Influence of exotic Mn on magnetic properties of YCo5.0-xMnxGa7.0 (0.05 <= x <= 3.0)

The YCo5.0-xMnxGa7.0 compounds crystallize with the ScFe6Ga6-type structure. The lattice of YCo5.0-xMnxGa7.0 expands with the increase of the Mn content for 0.05 <= x <= 2.5, but the lattice of YCo2.0Mn3.0Ga7.0 shrinks compared with YCo2.5Mn2.5Ga7.0. The shrinkage of the lattice is attributed to the magnetostriction of YCo2.0Mn3.0Ga7.0. The substitution of Mn for Co forms magnetic clusters in the antiferromagnetic matrix. The magnetic frustration results in the spin-glass-like behavior for 0.8 <= x <= 1.5 and the difference between zero-field-cooling (ZFC) and field-cooling (FC) magnetizations for x = 2.0, 2.5, and 3.0. A stable long-range magnetic ordering appears among the Mn-centered magnetic clusters with the ordering temperature 110 K for x = 2.0. The hump in the thermomagnetization of YCo3.0Mn2.0Ga7.0 can be attributed to the competitive effects between the thermal fluctuation and the enhanced magnetic interaction. Both the hump and the bifurcation between the ZFC and the FC magnetizations of YCo3.0Mn2.0Ga7.0 occur at lower temperatures as the applied field increases. On the two-step magnetization curve of YCo3.0Mn2.0Ga7.0, the inflection point at 4000 Oe is due to the coercive field, and the magnetic moments in the clusters are tilted to the applied field above 4000 Oe. The magnetic ordering temperature is further increased to 210 K for x = 2.5 and to 282 K for x = 3.0. The spontaneous magnetization of YCo2.0Mn3.0Ga7.0 is 0.575 mu B/f.u. at 5 K with a canted magnetic structure.

Structures of energy spectrum and persistent currents in mesoscopic rings with radial potential barrier in magnetic field

The energy spectrum and the persistent currents are calculated for finite-width mesoscopic annular structures with radial potential barrier in the presence of a magnetic field. The introduction of the tunneling barrier leads to the creation of extra edge states around the barrier and the occurrence of oscillatory structures superimposed on the bulk Landau level plateaus in the energy spectrum. We found that the Fermi energy E-F increases with the number of electrons N emerging many kinks. The single eigenstate persistent current exhibits complicated structures with vortex-like texture, ''bifurcation'', and multiple ''furcation'' patterns as N is increased. The total currents versus N display wild fluctuations.

Evolution of Localized Damage Zone in Heterogeneous Media

Evolution of localized damage zone is a key to catastrophic rupture in heterogeneous materials. In the present article, the evolutions of strain fields of rock specimens are investigated experimentally. The observed evolution of fluctuations and autocorrelations of strain fields under uniaxial compression demonstrates that the localization of deformation always appears ahead of catastrophic rupture. In particular, the localization evolves pronouncedly with increasing deformation in the rock experiments. By means of the definition of the zone with high strain rate and likely damage localization, it is found that the size of the localized zone decreases from the sample size at peak load to an eventual value. Actually, the deformation field beyond peak load is bound to suffer bifurcation, namely an elastic unloading part and a continuing but localized damage part will co-exist in series in a specimen. To describe this continuous bifurcation and localization process observed in experiments, a model on continuum mechanics is developed. The model can explain why the decreasing width of localized zone can lead stable deformation to unstable, but it still has not provided the complete equations governing the evolution of the localized zone.

Development and transition of spiral wave in the coupled Hindmarsh-Rose neurons in two-dimensional space

The dynamics and the transition of spiral waves in the coupled Hindmarsh-Rose (H-R) neurons in two-dimensional space are investigated in the paper. It is found that the spiral wave can be induced and developed in the coupled HR neurons in two-dimensional space, with appropriate initial values and a parameter region given. However, the spiral wave could encounter instability when the intensity of the external current reaches a threshold value of 1.945. The transition of spiral wave is found to be affected by coupling intensity D and bifurcation parameter r. The spiral wave becomes sparse as the coupling intensity increases, while the spiral wave is eliminated and the whole neuronal system becomes homogeneous as the bifurcation parameter increases to a certain threshold value. Then the coupling action of the four sub-adjacent neurons, which is described by coupling coefficient D', is also considered, and it is found that the spiral wave begins to breakup due to the introduced coupling action from the sub-adjacent neurons (or sites) and together with the coupling action of the nearest-neighbour neurons, which is described by the coupling intensity D.

Topological solitons in the CPN model

We investigate the solitons in the CPN supercript stop model in terms of the decomposition of gauge potential. Based on the phi-mapping topological current theory, the charge and position of solitons is determined by the properties of the typical component. Furthermore, the motion and the bifurcation of multi-soliton is discussed. And the knotted solitons in high dimension is explored also.