bursts of activity , be they regular or irregular , are a common occurrence in physical and biological systems . in recent years several models of bursting behavior in hydrodynamical systems have been described using ideas from dynamical systems theory . in this article
we review these and then describe a new mechanism ( moehlis & knobloch  @xcite , knobloch & moehlis  @xcite ) which provides an explanation for the bursting behavior observed in experiments on convection in @xmath0he/@xmath1he mixtures by sullivan & ahlers  @xcite .
this mechanism ope - rates naturally in systems with broken d@xmath2 symmetry undergoing a hopf bifurcation from a trivial state .
this symmetry , the symmetry group of a square , may be present because of the geometry of the system under consideration ( for example , the shape of the container ) but also appears in large aspect ratio systems with reflection symmetry ( landsberg & knobloch  @xcite ) . in either case
bursting arises as a result of the nonlinear interaction between two nearly degenerate modes with different symmetries .
this article is an expanded version of knobloch & moehlis @xcite .
as detailed further below , bursts come in many different forms , distinguished by their dynamic range ( i.e. , the range of amplitudes during each burst ) , duration , and recurrence properties . particularly important for the purposes of the present article
is the question of whether the observed bursts occur close to the threshold of a primary instability or whether they are found far from threshold . in the former case
a dynamical systems approach is likely to be successful : in this regime the spatial structure usually resembles the eigenfunctions of the linear problem and it is likely that only a small number of degrees of freedom participate in the burst .
such bursts take place fundamentally in the time domain with their spatial manifestation of secondary importance , in contrast to _
pulses _ which are structures localized in both time _ and _ space ; the latter are not considered here . since the equations governing the evolution of primary instabilities are often highly symmetric ( see crawford & knobloch @xcite ) global bifurcations are likely to occur and these serve as potential candidates for bursting mechanisms .
in contrast , bursts found far from threshold usually involve many degrees of freedom but even here some progress is sometimes possible .      the presence of coherent structures in a turbulent boundary layer is well established ( see , e.g. , robinson  @xcite and the collection of articles edited by panton  @xcite ) .
the space - time evolution of these structures is often characterized by intermittent bursting events involving low speed streamwise `` streaks '' of fluid .
specifically , let @xmath3 , and @xmath4 be the streamwise , wall normal , and spanwise directions with associated velocity components @xmath5 , and @xmath6 , respectively ; here @xmath7 is the mean flow . in a `` burst ''
the streak breaks up and low speed fluid moves upward away from the wall ( @xmath8 ) ; this is followed by a `` sweep '' in which fast fluid moves downward towards the wall ( @xmath9 ) .
after the burst / sweep cycle the streak reforms , often with a lateral spanwise shift . + a low - dimensional model of the burst / sweep cycle was developed by aubry et al .
@xcite ; further details and later references may be found in holmes et al .
@xcite . to construct such a model
the authors used a karhunen - love decomposition of the data to identify an energetically dominant empirical set of eigenfunctions , hereafter `` modes '' .
the original study of aubry et al .
@xcite used experimental data for pipe flow with @xmath10 , while later studies used data for channel flow from direct numerical simulation with @xmath11 and large eddy simulations with @xmath12 ( holmes et al .
the model was constructed by projecting the navier - stokes equation onto this basis and consists of a set of coupled odes for the amplitudes of these modes .
the fixed points of these equations are to be associated with the presence of coherent structures .
there are two types , related by half - wavelength spanwise translation .
+ numerical integration of the model reveals that these fixed points are typically unstable and that they are connected by a heteroclinic cycle .
in such a cycle the trajectory alternately visits the vicinities of the two unstable fixed points . in the model of aubry
@xcite this heteroclinic cycle is found to be structurally stable , i.e. , it persists over a _ range _ of parameter values .
this is a consequence of the o(2 ) symmetry of the equations inherited from periodic boundary conditions in the spanwise direction .
moreover , for the parameter values of interest this cycle is _ attracting _
, i.e. , it attracts all nearby trajectories .
since the transition from one fixed point to the other corresponds to a spanwise translation by half a wavelength , the recurrent excursions along such a heteroclinic cycle can be identified with the burst / sweep cycle described above . however , since this cycle is attracting , the time between successive bursts will increase as time progresses .
this is not observed and aubry et al .
@xcite appeal to the presence of a random pressure term modeling the effect of the outer fluid layer to kick the trajectory from the heteroclinic cycle . in the language of busse @xcite such a pressure term results in a _
statistical _ @xmath13 limit cycle , with the bursting events occurring randomly in time but with a well - defined mean rate .
the resulting temporal distribution of the burst events is characterized by a strong exponential tail , matching experimental observations .
attracting structurally stable heteroclinic cycles occur in a number of problems of this type , i.e. , mode interaction problems with o(2 ) symmetry ( armbruster et al .
@xcite , proctor & jones  @xcite , melbourne et al .
@xcite , steindl & troger  @xcite , krupa @xcite , hirschberg & knobloch @xcite ) .
experimental studies of plane couette flow and poiseuille ( pipe ) flow have shown that at high enough values of the reynolds number @xmath14 the basic laminar flow becomes turbulent . however , these laminar flows are linearly stable at all @xmath14 ( see , e.g. ,  drazin & reid  @xcite ) .
consequently , the transition to turbulence in these systems must arise from finite ( i.e. , not infinitesimal ) perturbations to the basic flow ; such transitions are called subcritical since the turbulent state exists for values of @xmath14 for which the laminar state is stable .
much recent work ( reviewed in baggett & trefethen  @xcite ) has emphasized the importance of the fact that the linear operators @xmath15 governing the evolution of perturbations of these flows are non - normal , i.e. , @xmath16 , where @xmath17 is the adjoint of @xmath15 .
linear systems with a non - normal @xmath15 can exhibit transient growth even though the laminar state is linearly stable ; if the growth is large enough , nonlinearities in the system may then trigger a transition to turbulence .
analysis of low - dimensional models supported by numerical simulations suggests that the minimum perturbation amplitude @xmath18 that results in turbulence scales as @xmath19 for some @xmath20 ( baggett & trefethen  @xcite ) .
thus , @xmath18 decreases rapidly with increasing @xmath14 , and the experimentally determined @xmath14 for transition should be that value at which @xmath18 is roughly equal to perturbations due to noise or imperfections in the system .
+ the turbulence excited by finite amplitude perturbations in shear flows often takes the form of turbulent spots which can move , grow , split , and merge ( see , e.g. , daviaud et al .
turbulent spots can also burst intermittently .
for example , bottin et al .
@xcite observed intermittent turbulent bursts in plane couette flow for @xmath21 ; these bursts were triggered by a spanwise wire through the center of the channel , and could be localized in space by introducing a bead into the central plane .
we focus here on this burst - like behavior , also found in low - dimensional models of the subcritical transition to turbulence .
the important issue in such models is the attractor to which the flow settles after an appropriate finite perturbation to the basic laminar flow .
this will depend crucially on nonlinear terms in the equations , a point emphasized , for example , by waleffe  @xcite and dauchot & manneville  @xcite .
if this attractor is a fixed point , the system settles into a steady state in which bursts do not occur . on the other hand , if this attractor is a limit cycle it may be appropriate to interpret the resulting behavior as burst - like .
this is the case in the model studied in waleffe  @xcite in which the amplitudes of streamwise streaks , streamwise rolls , and the streak instability can periodically undergo short - lived explosive growths at the expense of the mean shear .
waleffe  @xcite refers to this as a self - sustaining process .
since the laminar state is linearly stable , such limit cycles can not bifurcate off the laminar state . instead ,
as pointed out in waleffe  @xcite , they may be born in a hopf bifurcation ( from a fixed point other than that corresponding to the laminar state ) or a homoclinic bifurcation in which a trajectory homoclinic to a fixed point forms at some value of @xmath14 such that for smaller ( larger ) values of @xmath14 a limit cycle exists ( does not exist ) , or vice versa .
near such a homoclinic bifurcation , bursts are expected because the trajectory spends a long time near the fixed point , bursting away and returning in a periodic or chaotic fashion , depending on the eigenvalues of the fixed point .
other generic codimension one bifurcations leading to the appearance or disappearance of limit cycles are saddle - node bifurcations of limit cycles and saddle - nodes in which a pair of fixed points appears on the limit cycle ( guckenheimer et al .
limit cycles could also come in from or go off to infinity .
in other models chaotic attractors have been found ( see , e.g. ,  gebhardt & grossman @xcite and baggett et al .
@xcite ) and these may also give rise to burst - like behavior .
the relation of the bursts described in this section ( which occur for moderate values of @xmath14 , e.g. , @xmath22 ) to those which occur in the turbulent boundary layer ( which occur for high values of @xmath14 ) is not completely clear .
however , the similarities in the phenomenology have led waleffe  @xcite to propose that the self - sustaining process will continue to have importance in the near - wall region of high @xmath14 flows .
another mechanism involving heteroclinic connections , distinct from that discussed in section [ section : turb_bl ] , has been investigated by newell et al .
@xcite as a possible model for spatio - temporal intermittency in turbulent flow .
the authors suggest that such systems may be viewed as nearly hamiltonian except during periods of localized intense dissipation .
a related `` punctuated hamiltonian '' approach to the evolution of two - dimensional turbulence has met with considerable success ( carnevale et al .  @xcite and weiss & mcwilliams  @xcite ) .
for their description newell et al .
divide the instantaneous states of the flow into two categories , a turbulent soup ( ts ) characterized by weak coherence , and a singular ( s ) state characterized by strong coherence , and suppose that the ts and s states are generalized saddles in an appropriate phase space .
furthermore , they suppose that in the hamiltonian limit the unstable manifold of ts ( s ) intersects transversally the stable manifold of s ( ts ) .
if the constant energy surfaces are noncompact ( i.e. ,  unbounded ) , the evolution of the hamiltonian system may take the system into regions of phase space with very high ( `` infinite '' ) velocities and small scales .
these regions are identified with the s states and high dissipation .
in such a scenario the strong dissipation events are therefore identified with excursions along heteroclinic connections to infinity .
perturbations to the system ( such as the addition of dissipative processes ) may prevent the trajectory from actually reaching infinity , but this underlying unperturbed structure implies that large excursions are still possible .
newell et al .
apply these ideas to the two - dimensional nonlinear schrdinger equation ( nlse ) with perturbations in the form of special driving and dissipative terms which act at large and small scales , respectively .
here s consists of `` filament '' solutions to the unperturbed nlse which become singular in finite time and represent coherent structures which may occur at any position in the flow field .
when the solution is near s a large portion of the energy is in small scales ; for the perturbed equations the dissipative term then becomes important so that the filament solution is approached but collapses before it is reached .
this leads to a spatially and temporally random occurrence of localized burst - like events for the perturbed equation .
the rate of attraction at s is determined by the faster than exponential rate at which the filament becomes singular , while the rate of repulsion at s is governed by the dissipative process and hence is unrelated to the rate of attraction .
this bursting mechanism shares characteristics with that described in kaplan et al .
@xcite in which solutions of a single complex ginzburg - landau equation with periodic boundary conditions undergo faster than exponential bursting due to a destabilizing nonlinearity and collapse due to strong nonlinear dispersion ( see also bretherton & spiegel  @xcite ) .
a study of a generalization of burger s equation modeling nonlocality effects suggests the presence of burst - like events through a similar scenario ( makarenko et al .
@xcite ) .      the kolmogorov flow @xmath23 is an exact solution of the two - dimensional incompressible navier - stokes equation with unidirectional forcing @xmath24 at wavenumber @xmath25 : @xmath26 . with increasing reynolds number @xmath27
this flow becomes unstable , and direct numerical simulation with @xmath28-periodic boundary conditions shows that for moderately high reynolds numbers and @xmath29 the resulting flow is characterized by intermittent bursting ( she @xcite , nicolaenko & she @xcite , armbruster et al .
a burst occurs when the system evolves from a coherent vortex - like modulated traveling wave ( mw ) to a spatially disordered state following transfer of energy from large to small scales .
the system then relaxes to the vicinity of another symmetry - related mw state , and the process continues with bursts occurring irregularly but with a well - defined mean period .
+ the details of what actually happens appear to depend on the value of @xmath25 because @xmath25 determines the symmetry of the nonlinear equation describing the evolution of the perturbation streamfunction @xmath30 about the kolmogorov flow .
although there is no compelling reason for it , all simulations of this equation have been performed with @xmath28-periodic boundary conditions in both directions . with these boun - dary conditions
this equation has a symmetry which is the semi - direct product of the dihedral group d@xmath31 ( generated by the actions @xmath32 and @xmath33 ) and the group so(2 ) ( representing the symmetry under translations @xmath34const ) . in the simplest case , @xmath35 , this symmetry group is isomorphic to the direct product of the group o(2 ) of rotations and reflections of a circle and the group z@xmath36 representing reflections in the @xmath37-direction . unfortunately , when @xmath35 the kolmogorov flow with @xmath28-periodic boundary conditions is not unstable for any value of @xmath14 ( meshalkin & sinai  @xcite , green  @xcite , marchioro  @xcite ) , and one is forced to consider @xmath29 .
armbruster et al .  @xcite analyzed carefully the @xmath38 case and showed that while a heteroclinic cycle between the mw states does form , it is not structurally stable ; the bursts are therefore not produced by a mechanism of the type described in section [ section : turb_bl ] .
it is possible , however , that the onset of bursting is associated with a symmetry - increasing bifurcation at @xmath39 ( see , e.g. , rucklidge & matthews  @xcite ) .
this would explain why the system stays in a single mw state for @xmath14 just below @xmath40 but visits the vicinity of different but symmetry - related mw states for @xmath14 just above @xmath40 .
however , despite much work a detailed understanding of the bursts in this system remains elusive .
+ an alternative approach is to consider periodic domains with different periodicities in the two directions . in particular ,
if we consider the domain @xmath41 with @xmath29 the symmetry group remains o(2)@xmath42z@xmath36 but sufficiently long perturbations now grow .
the unstable modes are either even or odd under the reflection @xmath43 with respect to a suitable origin .
mode interaction between these two steady modes can result in a sequence of transitions summarized in figure [ phil ] ( hirschberg & knobloch  @xcite ) : the kolmogorov flow loses stability to an even mode ( z ) , followed by a steady state bifurcation to a mixed parity state ( mm@xmath44 ) .
since each of these states is defined to within a translation in @xmath45 modulo @xmath28 , we say that it forms a _ circle _ of states .
the mm@xmath44 state then loses stability in a further steady state bifurcation to a traveling wave ( tw ) which in turn loses stability at a hopf bifurcation to a mw .
the mw two - torus terminates in a collision with the two circles of pure parity states , forming an attracting structurally stable heteroclinic cycle connecting them and their quarter - wavelength translates ( hirschberg & knobloch  @xcite ) .    in this regime
the behavior would resemble that found in the numerical simulations , with higher modes kicking the system away from this cycle .
indeed this sequence of transitions echoes the results obtained by she  @xcite and nicolaenko & she  @xcite for @xmath46 . while it is likely that the @xmath35 scenario is relevant to these calculations because of the tendency towards an inverse cascade in these two - dimensional systems , we emphasize that , as mentioned above , a careful analysis of the @xmath38 case by armbruster et al .
@xcite shows that while a heteroclinic cycle of the required type does indeed form , it is not structurally stable . the case @xmath47 has also been studied ( platt et al .
@xcite ) and a similar sequence of transitions found .
undoubtedly simulations on the domain @xmath41 would shed new light on the problem , cf .  posch & hoover  @xcite .
+ the group o(2)@xmath42z@xmath36 also arises in convection in rotating straight channels ( knobloch  @xcite ) and in natural convection in a vertical slot ( xin et al .
@xcite ) . in both of these cases
the linear eigenfunctions are either even or odd with respect to a rotation by @xmath48 and the systems exhibit similar transitions .
related behavior has recently been observed in three - dimensional magnetoconvection with periodic boundary conditions on a square lattice . unlike the kolmogorov
flow the evolution equation describing a steady state secondary instability of a square pattern is isotropic in the horizontal
. numerical simulations ( rucklidge et al .
@xcite ) show intermittent breakdown of a square pattern followed by its restoration modulo translation .
the taylor - couette system consists of concentric cylinders enclosing a fluid - filled annulus ( see , e.g. , tagg @xcite ) .
the cylinders can be rotated independently . in the counterrotating regime the first state consists of spiral vortices of either odd or even parity with respect to midheight . slightly above onset
the flow resembles interpenetrating spiral ( ips ) flow which may be intermittently interrupted by bursts of turbulence which fill the entire flow field ( hamill  @xcite , coughlin et al .
@xcite ) . with periodic boundary conditions in the axial direction numerical simulations by coughlin et al .
@xcite show that the ips flow is temporally chaotic and consists of coexisting modes with different axial and azimuthal wavenumbers .
this flow is confined primarily to the vicinity of the inner cylinder where the axisymmetric base flow is subject to an inviscid rayleigh instability .
coughlin et al .
@xcite conclude that the onset of turbulence is correlated with a secondary instability of one of the coexisting modes of the ips flow , namely the basic spiral vortex flow with azimuthal wavenumber @xmath49 . indeed , when this mode is taken as the initial flow for parameters chosen such that the full ips flow undergoes bursts , a secondary hopf bifurcation from this state with the same azimuthal wavenumber but four times the axial wavelength is identified .
the secondary instability grows in amplitude and ultimately provides a finite amplitude perturbation to the inviscidly stable flow near the outer cylinder , triggering a turbulent burst throughout the whole apparatus . during a burst small scales
are generated throughout the apparatus leading to a rapid collapse of the turbulence and restoration of the ips flow ; the bursting process then repeats roughly periodically in time ( see figure  [ coughlin ] ) .
+ as discussed in section [ section : our_bursts.1 ] , in a finite taylor - couette apparatus there is a natural mechanism for generating bursts close to onset .
this mechanism operates in the regime @xmath50 , where @xmath18 measures the fractional distance above threshold for the primary instability and @xmath51 is the aspect ratio of the annulus ( landsberg & knobloch  @xcite , renardy  @xcite ) . for larger @xmath18
the influence of the boundaries no longer extends throughout the apparatus and is confined to boundary layers near the top and bottom . in this regime the dynamics in the bulk may be described by imposing periodic boundary conditions with period that is a multiple of the wavelength of the primary instability .
the success of the simulation of the observed turbulent bursts using such periodic boundary conditions ( coughlin et al .
@xcite ) suggests that these bursts occur too far above threshold to be explained by the mechanism described in section [ section : our_bursts.1 ] below .
this suggestion is supported by figure [ tc_bursts_onset ] which compares the location of the regime where this mechanism may be expected to operate with @xmath52 , the experimental value of @xmath18 for the onset of bursts , as a function of @xmath51 .
the latter is obtained using the approximation @xmath53 , where @xmath54 and @xmath55 are the inner cylinder reynolds numbers for the onset of ips flow and bursts , and @xmath56 ( hamill  @xcite ) .
this is a good approximation for strongly counterrotating cylinders because the ips flow sets in for inner cylinder reynolds number only about 1% above the primary onset to spiral vortices ( hamill  @xcite ) .
the figure suggests that the observed bursts may fall within the range of validity of this theory for @xmath51 only slightly smaller than those used in the experiments ( @xmath57 ,  @xmath58 ) ; for such bursts the presence of endwalls should become significant .
the term intermittency refers to occasional , non - periodic switching between different types of behavior ; such transitions may be viewed as bursts . in this section
we describe different types of intermittency , dividing them into three classes .
it is often the case that fluid dynamical systems may be modelled by low - dimensional poincar maps .
fixed points of such maps correspond to limit cycles of an appropriate dynamical system ; when these lose stability intermittency may result . without loss of generality
we assume that for @xmath59 the system is in a `` laminar '' state corresponding to a stable fixed point of the map , and the onset of intermittency occurs at @xmath60 . for @xmath61 the state of the system will resemble the stable laminar state which existed for @xmath59 for some time but
will intermittently undergo `` bursts '' away from this state , followed by `` reinjections '' to the vicinity of the laminar state .
the three types of intermittency identified by pomeau & manneville @xcite correspond to three different ways in which the laminar state ceases to exist or loses stability at @xmath62 : type i intermittency results when a stable and unstable limit cycle annihilate in a saddle - node bifurcation , type ii results when the limit cycle loses stability in a subcritical hopf bifurcation , while in type iii the limit cycle loses stability in an inverse period - doubling bifurcation .
type i and type iii intermittency have been observed , for example , in rayleigh - bnard convection ( see berg et al .
@xcite and dubois et al .
@xcite , respectively ) , while type ii intermittency has been observed in a hot wire experiment ( ringuet et al .
a different type of intermittency , called type x , was observed by price & mullin @xcite in a variant of the taylor - couette system ; this is similar to type i intermittency but involves a hysteretic transition due to the nature of the reinjection .
another type of intermittency , called type v , was introduced in bauer et al .
@xcite and he et al .
@xcite ; this is similar in spirit to type i intermittency but involves one - dimensional maps which are nondifferentiable or discontinuous .
these different types of intermittency have distinct properties , such as the scaling behavior of the average time between bursts with @xmath63 .
a crisis is a sudden change in a strange attractor as a parameter is varied ( grebogi et al .
there are three types of crises , and without loss of generality we assume that the crisis occurs as @xmath64 is increased through @xmath65 . in a boundary crisis , at @xmath66 the strange attractor collides with a coexisting unstable periodic orbit which lies on the boundary of the basin of attraction of the strange attractor .
this leads to the destruction of the strange attractor , but chaotic transients will still exist .
intermittency is not associated with this crisis . in an interior crisis
the strange attractor collides with a coexisting unstable periodic orbit at @xmath66 , but here this leads to a widening rather than the destruction of the strange attractor . for @xmath64
slightly larger than @xmath65 the trajectory stays near the region of phase space occupied by the strange attractor before the crisis for long times , but intermittently bursts into a new region . in an attractor merging crisis two strange attractors with basins of attraction separated by a basin boundary
are present when @xmath67 . at @xmath65
the two strange attractors simultaneously touch the basin boundary and `` merge '' to form a larger strange attractor for @xmath68 .
such a crisis is often associated with a symmetry - increasing bifurcation ( chossat & golubitsky @xcite ) .
for @xmath68 there is a single strange attractor on which the trajectory intermittently switches between states resembling the distinct strange attractors which existed for @xmath69 .
a manifold in phase space is invariant if every initial condition on the manifold generates an orbit that remains on the manifold .
invariant manifolds often exist due to symmetries , but this is not necessary .
there are several mechanisms for intermittency involving strange attractors on an invariant manifold .    first
, suppose that as a bifurcation parameter @xmath64 is increased through @xmath70 , the strange attractor on the invariant manifold loses stability transverse to the manifold ; this is called a blowout bifurcation ( ott & sommerer  @xcite ) .
suppose that the dynamics _ within _
the invariant subspace do not depend on @xmath64 .
such a system thus has a skew - product structure ( platt et al .
@xcite ) , and @xmath64 is called a normal parameter ( ashwin et al .
if the blowout bifurcation is `` supercritical '' ( ott & sommerer  @xcite ) then for @xmath64 just above @xmath70 a trajectory will spend a long time near the invariant manifold , intermittently bursting away from it , only to return due to the presence of a reinjection mechanism .
this scenario is known as on - off intermittency , where the `` off '' ( `` on '' ) state corresponds to the system being near ( away from ) the invariant manifold ( see , e.g. , platt et al .
@xcite , venkataramani et al .
recently it has been shown that in appropriate circumstances a blowout bifurcation can lead to a structurally stable ( possibly attracting ) heteroclinic cycle between chaotic invariant sets ( ashwin & rucklidge @xcite ) .    for systems which lack a skew - product structure and
thus are governed by non - normal parameters , in - out intermittency is possible ( ashwin et al .
@xcite , covas et al .
@xcite ) . here
the attraction and repulsion to the invariant subspace are controlled by _
different _ dynamics .
this occurs when the attractor for dynamics restricted to the invariant subspace is smaller than the intersection of the attractor for the full dynamics with the invariant subspace , and may be viewed as a generalization of on - off intermittency in which the attraction to and repulsion from the invariant subspace are controlled by the _
same _ dynamics .
another mechanism for bursting occurs when the strange attractor on the invariant manifold attracts typical orbits near the surface but is unstable in the sense that there are unstable periodic orbits embedded within the chaotic set which are transversely repelling .
if the trajectory comes near such an unstable periodic orbit , there will be a burst away from the invariant surface .
such bursts may occur intermittently if noise is present or small changes ( called `` mismatch '' ) are made to the dynamical system that destroy the invariant surface .
the instability of the first of these orbits is known as a bubbling transition ( see , e.g. , ashwin et al .
@xcite , venkataramani et al .
@xcite ) .      in neural systems
, bursting refers to the switching of an observable such as a voltage or chemical concentration between an active state characterized by rapid ( spike ) oscillations and a rest state .
models of such bursting typically involve singularly perturbed vector fields in which system variables are classified as being `` fast '' or `` slow '' depending on whether or not they change significantly over the duration of a single spike .
the slow variables may then be thought of as slowly varying parameters for the equations describing the fast variables ( rinzel  @xcite , bertram et al .
@xcite , wang & rinzel  @xcite , guckenheimer et al .
as the slow variables evolve , the state of the system in the fast variables may change from a stable periodic orbit ( corresponding to the active state ) to a stable fixed point ( corresponding to the rest state ) and vice versa ; such transitions are often associated with a region of bistability for the periodic orbit and the fixed point , but need not be . mechanisms by which such transitions can occur repeatedly have been classified ( rinzel  @xcite , bertram et al .
@xcite , wang & rinzel  @xcite ) .
behavior of the time interval between successive spikes near a transition from the active to the rest state is discussed by guckenheimer et al .
@xcite ; in this paper the presence of a subcritical hopf - homoclinic bifurcation is also identified as a possible mechanism for the transition from an active to a rest state .
we now describe a bursting mechanism which involves the interaction between oscillatory modes in systems with approximate d@xmath2 symmetry , where d@xmath2 is the symmetry group of the square .
this mechanism can lead to bursts of large dynamic range very close to the instability onset ( moehlis & knobloch  @xcite , knobloch & moehlis  @xcite ) and is expected to be relevant in many different systems with approximate d@xmath2 symmetry . this symmetry may be present for obvious or subtle reasons , as the following discussion demonstrates .
consider binary fluid convection in a system of large but finite aspect ratio .
if the separation ratio is sufficiently negative the system is overstable , i.e. , the primary instability is via a hopf bifurcation .
this is the case for the @xmath0he/@xmath1he mixture used by sullivan & ahlers  @xcite in their experiment carried out in a rectangular container @xmath71 with @xmath72 . in this experiment sullivan & ahlers
observed that immediately above threshold ( @xmath73 ) convective heat transport may take place in a sequence of irregular bursts of large dynamic range despite constant heat input ( see figure  [ sullivan ] ) .    in this system the presence of sidewalls
destroys translation symmetry in the @xmath45 direction which would be present if the system were unbounded , but with identical boundary conditions at the sidewalls the system retains a reflection symmetry about @xmath74 ; the primary modes are thus either even or odd with respect to this reflection ( dangelmayr & knobloch @xcite ) . numerical simulations of the appropriate partial differential equations suggest that the bursts observed in the experiments involve the interaction between the first odd and even modes of the system ( jacqmin & heminger @xcite ; see also batiste et al .
@xcite as described in section [ section : linear ] ) .
thus , to describe the dynamical behavior near threshold we suppose that the perturbation from the conduction state takes the form @xmath75 where @xmath76 , @xmath77 .
following landsberg & knobloch @xcite we now derive amplitude equations describing the evolution of @xmath78 and @xmath79 using symmetry arguments .
to do this , we first briefly review the topic of bifurcations in systems with symmetry ( see , e.g. , golubitsky et al .
@xcite and crawford & knobloch @xcite ) .
suppose that @xmath80 where @xmath81 and @xmath82 represent dependent variables and system parameters , respectively .
let @xmath83 describe a linear group action on the dependent variables .
we say that if @xmath84 for all @xmath83 then ( [ equivariant ] ) is _ equivariant _ with respect to the group @xmath85 .
this is equivalent to the statement that if @xmath86 is a solution to ( [ equivariant ] ) then so is @xmath87 .
for example , if a system is equivariant under left - right reflections and a right - moving wave exists as a solution , then a reflection - related left - moving wave will also exist as a solution . for the amplitudes @xmath78 and @xmath79 the requirement that a reflected state ( obtained by letting @xmath88 in ( [ perturb ] ) )
also be a state of the system gives the requirement that the amplitude equations be equivariant with respect to the group action @xmath89 moreover , as argued by landsberg & knobloch @xcite , the equations for the formally infinite system can not distinguish between the two modes , i.e. , in this limit the amplitude equations must also be equivariant with respect to the group action @xmath90 which we call an interchange symmetry .
these two operations generate together the symmetry group d@xmath2 . for a container with large but finite length ,
this symmetry will be weakly broken ; in particular , the even and odd modes typically become unstable at slightly different rayleigh numbers and with slightly different frequencies ( see section [ section : linear ] ) .
the resulting equations are thus close to those for a 1:1 resonance , but with a special structure dictated by the proximity to d@xmath2 symmetry .
finally , we may put the equations for @xmath78 and @xmath79 into normal form by performing a series of near - identity nonlinear transformations of the dependent variables so as to simplify the equations as much as possible ( see , e.g. , guckenheimer & holmes @xcite ) .
the normal form equations have the additional symmetry ( elphick et al .
@xcite ) @xmath91 which may be interpreted as a phase shift symmetry . truncating the resulting equations at third order
we obtain @xmath92 z_+ + a ( |z_+|^2 + |z_-|^2 ) z_+ \nonumber \\   & & + b |z_+|^2 z_+ + c
\bar{z}_+ z_- ^2   \label{z+}\end{aligned}\ ] ] @xmath93 z_- + a ( |z_+|^2 + |z_-|^2 ) z_- \nonumber \\   & & + b |z_-|^2 z_- + c
\bar{z}_- z_+ ^2 . \label{z-}\end{aligned}\ ] ] here @xmath94 measures the difference in frequency between the two modes at onset , and @xmath95 measures the difference in their linear growth rates . under appropriate nondegeneracy conditions ( which we assume here ) we may neglect all interchange symmetry - breaking contributions to the nonlinear terms . in
the following we consider the regime in which @xmath64 , @xmath95 , and @xmath94 are all of the same order ; in the large aspect ratio binary fluid convection context this will occur when these quantities are all @xmath96 ( see section [ section : linear ] ) .
we will see that when @xmath95 and/or @xmath94 are nonzero , ( [ z+],[z- ] ) have bursting solutions .
thus , the bursting mechanism may be viewed as an interaction between spontaneous symmetry breaking ( in which the trivial conduction state loses stability to a convecting state with less symmetry ) and forced symmetry breaking ( in which the presence of sidewalls make @xmath95 and/or @xmath94 nonzero , thereby breaking the d@xmath2 symmetry ) .
the introduction of small symmetry - breaking terms is also responsible for the possibility of complex dynamics in other systems that would otherwise behave in a regular manner ( dangelmayr & knobloch @xcite , lauterbach & roberts @xcite , knobloch @xcite , hirschberg & knobloch @xcite ) .    to identify the bursts we introduce the change of variables @xmath97 and a new time - like variable @xmath98 defined by @xmath99 . in terms of these variables
( [ z+],[z- ] ) become @xmath100 \label{rho } \\ & & - 2 ( \lambda + \delta \lambda \cos \theta ) \rho^2 \nonumber \\
\frac{d \theta}{d \tau } & = & \sin \theta [ \cos \theta ( -b_r + c_r \cos 2 \phi ) - c_i \sin 2 \phi ] \nonumber \\ & & - 2 \delta \lambda \;\rho \sin \theta \label{theta } \\ \frac{d \phi}{d \tau } & = & \cos \theta ( b_i - c_i \cos 2 \phi )   - c_r \sin 2 \phi + 2 \delta \omega \;\rho , \label{phi}\end{aligned}\ ] ] where @xmath101 , etc . there is also a decoupled equation for @xmath102 so that fixed points and periodic solutions of equations ( [ rho]-[phi ] ) correspond , respectively , to periodic solutions and two - tori in equations ( [ z+],[z- ] ) .    in the following we measure the amplitude of the disturbance by @xmath103 ; thus @xmath104 corresponds to _ infinite _ amplitude states . equations
( [ rho]-[phi ] ) show that the restriction to the invariant subspace @xmath105 is equivalent to taking @xmath106 in ( [ theta],[phi ] ) .
the resulting d@xmath2-symmetric problem has three generic types of fixed points ( swift @xcite ) :    @xmath107 @xmath108 solutions with @xmath109    @xmath107 @xmath110 solutions with @xmath111    @xmath107 @xmath112 solutions with @xmath113 .    in the binary fluid context
the @xmath108 , @xmath110 and @xmath112 solutions represent mixed parity traveling wave states localized near one of the container walls , mixed parity chevron ( or counterpropagating ) states , and pure even ( @xmath114 ) or odd ( @xmath115 ) parity chevron states , respectively ( landsberg & knobloch @xcite ) . such states are shown in figure [ uvw ] using the approximate eigenfunctions @xmath116 where @xmath117 , @xmath118 and @xmath119 .    depending on the coefficients @xmath120 ,  @xmath121 and
@xmath122 the subspace @xmath123 may contain additional fixed points and/or limit cycles ( swift @xcite ) . in our scenario ,
a burst occurs for @xmath124 when a trajectory follows the stable manifold of a fixed point ( or a limit cycle ) @xmath125 that is _ unstable _ within @xmath123 .
the instability within @xmath123 then kicks the trajectory towards another fixed point ( or limit cycle ) @xmath126 .
if this point has an unstable @xmath127 eigenvalue the trajectory escapes from @xmath123 towards a finite amplitude ( @xmath128 ) state , forming a burst .
if @xmath129 and/or @xmath130 this state may itself be unstable to perturbations of type @xmath131 and the process then repeats .
this bursting behavior is thus associated with a codimension one heteroclinic cycle between the _ infinite amplitude _ solutions
@xmath131 and @xmath132 ( moehlis & knobloch  @xcite , knobloch & moehlis  @xcite ) .
examples of such cycles are shown in figure [ hetero_fig ] .    since in such cycles the trajectory reaches infinity in finite time
the heteroclinic cycle actually describes bursts of _ finite _ duration ( moehlis & knobloch @xcite ) .
+ for such a heteroclinic cycle to form it is required that at least one of the branches in the d@xmath2-symmetric system be subcritical ( @xmath131 ) and one supercritical ( @xmath132 ) .
based on the @xmath0he/@xmath1he experiments , we focus on parameter values for which the @xmath108 solutions are subcritical and the @xmath110 , @xmath112 solutions supercritical when @xmath133 ( moehlis & knobloch @xcite ) . when @xmath129 and/or @xmath130 , two types of oscillations in @xmath134 are possible :    @xmath107 rotations ( see figure [ rotation ] )    @xmath107 librations ( see figure [ libration ] ) .
-0.2 in    -0.2 in    for @xmath124 these give rise , under appropriate conditions , to sequences of large amplitude bursts arising from repeated excursions towards the infinite amplitude ( @xmath104 ) @xmath108 solutions .
irregular bursts are also readily generated : figure [ chaotic_rotation ] shows bursts arising from chaotic rotations .
figure [ bds ] provides a partial summary of the different solutions of ( [ rho]-[phi ] ) and their stability properties ; much of the complexity revealed in these figures is due to the shilnikov - like properties of the heteroclinic cycle ( moehlis & knobloch  @xcite , knobloch & moehlis  @xcite ) .
we now focus on the physical manifestation of the bursts . in figure
[ space_time ] we show the solutions of figures [ rotation ] and [ libration ] in the form of space - time plots using the approximate eigenfunctions ( [ efunctions ] ) .
the bursts in figure [ space_time](a ) are generated as a result of successive visits to _
different _ ( but symmetry - related ) infinite amplitude @xmath108 solutions , cf .
figure  [ rotation ] ; in figure  [ space_time](b ) the generating trajectory makes repeated visits to the _ same _ infinite amplitude @xmath108 solution , cf .
figure  [ libration ] .
the former state is typical of the blinking state identified in binary fluid and doubly diffusive convection in rectangular containers ( kolodner et al .
@xcite , steinberg et al .
@xcite , predtechensky et al .
it is likely that the irregular bursts reported in sullivan & ahlers  @xcite are due to such a state .
the latter is a new state which we call a _ winking _ state ; winking states may be stable but often coexist with stable chevron - like states which are more likely to be observed in experiments in which the rayleigh number is ramped upwards ( cf .
figure [ bds ] ) . for other values of @xmath95 and @xmath94
it is also possible to find stable chaotic winking states and states which are neither purely blinking nor purely winking ( see figure [ new_solutions ] ) .
the bursts described above are the result of oscillations in amplitude between two modes of opposite parity and `` frozen '' spatial structure .
consequently the above burst mechanism applies in systems in which bursts occur very close to threshold .
this occurs not only in the convection experiments already mentioned but also in the mathematically identical ( counterrotating ) taylor - couette system where counterpropagating spiral vortices play the same role as traveling waves in convection ( andereck et al .
@xcite , pierce & knobloch @xcite ) . in slender systems , such as the convection system described above or a long taylor - couette apparatus ,
a large aspect ratio @xmath51 is required for the presence of the approximate d@xmath2 symmetry .
if the size of the d@xmath2 symmetry - breaking terms @xmath129 , @xmath135 is increased too much the bursts fade away and are replaced by smaller amplitude , higher frequency states ( see figure [ fade ] ) .    indeed , if @xmath136 , averaging eliminates the @xmath122 terms responsible for the bursts ( landsberg & knobloch @xcite ) . from these considerations , we conclude that bursts will not be present if @xmath51 is too small or @xmath18 too large .
however , the mechanism is quite robust and even for @xmath136 it may still be possible to choose @xmath64 values so that bursts of large dynamic range occur ( moehlis & knobloch  @xcite ) .
+ it is possible that the burst amplitude can become large enough for secondary instabilities not captured by the ansatz ( [ perturb ] ) to be triggered .
such instabilities could occur on very different scales and result in _ turbulent _ rather than just large amplitude bursts .
however , it should be emphasized that the physical amplitude of the bursts is @xmath137 and so approaches zero as @xmath138 , cf .
eq .  ( [ perturb ] ) .
thus despite their large dynamic range ( cf .
figure [ very_large ] ) , the bursts are fully and correctly described by the asymptotic expansion that leads to ( [ z+],[z- ] ) .
in particular , as shown in moehlis & knobloch  @xcite , the mechanism is robust with respect to the addition of small fifth order terms .
however , the effects of including d@xmath2-symmetry - breaking terms in the cubic terms in ( [ z+],[z- ] ) have not been analyzed ; these terms can dominate the symmetry - breaking terms retained in the linear terms when @xmath127 is small .      in view of the motivation for studying systems with approximate d@xmath2 symmetry described above , it is of interest to examine carefully the properties of the linear stability problem for binary fluid convection in finite containers . in the boussinesq approximation
this system is described by the nondimensionalized equations ( clune & knobloch  @xcite ) @xmath139 \hat{\bf z } + \sigma \nabla^2 { \bf u},\\ \partial_t \theta + ( { \bf u } \cdot
\nabla ) \theta & = & w + \nabla^2 \theta , \\ \partial_t \eta + ( { \bf u}\cdot \nabla)\eta & = & \tau \nabla^2 \eta + \nabla^2 \theta,\end{aligned}\ ] ] together with the incompressibility condition @xmath140 here @xmath141 is the velocity field in @xmath142 coordinates , @xmath143 , @xmath144 and @xmath122 denote the departures of the pressure , temperature and concentration fields from their conduction profiles , and @xmath145 .
these equations are to be solved in the rectangular domain @xmath146 .
the system is specified by four dimensionless parameters in addition to the aspect ratio @xmath51 : the separation ratio @xmath147 , the prandtl and lewis numbers @xmath148 , @xmath98 , and the rayleigh number @xmath149 .
the boundary conditions appropriate to the experiments are no - slip everywhere , with the temperature fixed at the top and bottom and no sideways heat flux .
the final set of boundary conditions is provided by the requirement that there is no mass flux through any of the boundaries .
the boundary conditions are thus @xmath150 @xmath151 here @xmath152 denotes the boundary of @xmath153 .
figure [ linear ] shows the results of solving the _ linear _ problem describing the stability properties of the conduction state @xmath154 for parameter values used by sullivan & ahlers  @xcite in their @xmath0he/@xmath1he experiment : @xmath155 , @xmath156 , @xmath157 .
the figure shows the neutral stability curves and corresponding frequencies for the first _ four _ modes in the range @xmath158 for @xmath159 and @xmath160 .
observe that when @xmath161 is sufficiently small the first two families of neutral curves are separated by a gap that is much larger than the amplitude of the `` braids '' within each family ( figure [ linear](a ) ) .
this is typical of what happens in rayleigh - bnard convection with non - neumann boundary conditions ( hirschberg & knobloch @xcite ) and makes it easy to justify projecting the fluid equations onto the first two modes that become unstable .
however , the situation is not so simple .
this is because in the case of overstability this behavior does not persist for all @xmath51 or all values of @xmath161 . for larger values of these parameters
the results take instead the form shown in figures [ linear](c , d ) which show the linear stability results for @xmath160 and the same range of values of @xmath51 as figures  [ linear](a , b ) .
the modes from the different families now cross and the first unstable mode belongs to successively higher and higher families when extrapolated to small @xmath161 ( batiste et al .
@xcite ) . figure  [ linear](c ) shows the crossing of two even modes involving a _
double hopf bifurcation ( figure  [ linear](d ) ) . as discussed in detail in batiste et al .
@xcite the transition between these two types of behavior is mediated by a _ resonant _
1:1 mode crossing at a somewhat smaller value of @xmath161 .
the experimental value of the separation ratio from sullivan & ahlers  @xcite , @xmath160 , therefore corresponds to the `` crossing '' case and the projection of the equations onto two modes can not be rigorously justified except in the neighborhood of mode crossing points .
+ we denote the growth rates and frequencies of the modes @xmath162 by @xmath163 and @xmath164 . for large aspect ratios ,
the mode frequencies must go like @xmath165 .
the fact that the frequency curves in figure  [ linear](d ) are essentially parallel `` straight lines '' implies that @xmath166 .
therefore , @xmath167 for large @xmath51 ( batiste et al .
moreover , as argued in landsberg & knobloch  @xcite , the parabolic minimum of the neutral stability curve leads to the expectation that @xmath168 .
thus , in the range @xmath169 , @xmath64 , @xmath95 , and @xmath135 are all of the _ same _ order as @xmath170 , as required for the applicability of equations ( [ z+ ] , [ z- ] ) .
of course , close enough to the mode crossing point @xmath171 , and in this region averaging methods can be used to eliminate the @xmath172 terms from the mode interaction equations ( landsberg & knobloch @xcite ) .
however , for typical values of @xmath95 it appears likely that the system is correctly described by equations ( [ z+ ] , [ z- ] ) , as hypothesized in landsberg & knobloch @xcite and moehlis & knobloch @xcite .
+ the first odd and even temperature eigenfunctions for @xmath160 are shown in figure  [ efunction_plot ] in the form of a space - time diagram , with time increasing upward .    as in the approximate expression ( [ efunctions ] )
the eigenfunction consists of waves propagating outwards from the center of the container .
the eigenfunction amplitude has a local minimum at the center and increases outwards , peaking near the sidewalls .
this type of eigenfunction was anticipated by cross @xcite and is characteristic of eigenfunctions in systems with _ positive _ group velocity ( although , strictly speaking , in a finite system one can not define a group velocity since the allowed wave number is quantized by the sidewalls as well as being nonuniform ) .
however , for the present purposes the most important observation is that for aspect ratios as large as this , the odd and even eigenfunctions are essentially indistinguishable , as hypothesized by landsberg & knobloch  @xcite .
there are a number of other systems of interest where an approximate d@xmath2 symmetry arises in a natural way and the bursting mechanism described in section [ section : our_bursts.1 ] may be relevant .
these include overstable convection in small aspect ratio containers with nearly square cross - section ( armbruster @xcite ) and more generally any partial differential equation on a nearly square domain describing the evolution of an oscillatory instability , cf .
ashwin & mei  @xcite .
other systems in which our bursting mechanism might be detected are electrohydrodynamic convection in liquid crystals ( silber et al .
@xcite ; t. peacock , private communication ) , lasers ( feng et al .
@xcite ) , spring - supported fluid - conveying tubes ( steindl & troger  @xcite ) , and dynamo theories of magnetic field generation in the sun ( knobloch & landsberg  @xcite , knobloch et al .
@xcite ) .
perhaps more interesting is the possibility that large scale spatial modulation due to distant walls may produce bursting in a fully nonlinear state with d@xmath2 symmetry undergoing a symmetry - breaking hopf bifurcation . as an example we envisage a steady pattern of fully nonlinear two - dimensional rolls . with periodic boundary conditions with period four times the basic roll period , the roll pattern has d@xmath2 symmetry since the pattern is preserved under spatial translations by 1/4 period and a reflection .
if such a pattern undergoes a secondary hopf bifurcation with a _ spatial _ floquet multiplier @xmath173 , the hopf bifurcation breaks d@xmath2 symmetry .
if the invariance of the basic pattern under translations by 1/4 period is only approximate ( this would be the case if the roll amplitude varied on a slow spatial scale ) , the d@xmath2 symmetry itself would be weakly broken and the new mechanism described above could operate . also of interest is the faraday system in a nearly square container . in this system gravity - capillary waves
are excited on the surface of a viscous fluid by vertical vibration of the container , usually as a result of a subharmonic resonance .
simonelli & gollub @xcite studied the effect of changing the shape of the container from a square to a slightly rectangular container , focusing on the @xmath174 , @xmath175 interaction in this system .
these modes are degenerate in a square container and only pure and mixed modes were found in this case . in a slightly rectangular container the degeneracy between these modes
is broken , however , and in this case a region of quasiperiodic and chaotic behavior was present near onset .
when these oscillations first appear they take the form of relaxation oscillations in which the surface of the fluid remains flat for a long time before a `` large wave grows , reaches a maximum , and decays , all in a time short compared with the period '' .
the duration of the spikes is practically independent of the forcing amplitude , while the interspike period appears to diverge as the forcing amplitude _
decreases_. the spikes themselves possess the characteristic asymmetry seen in figures  [ rotation ] and [ libration ] .
this behavior occurs when the forcing frequency lies _ below _ the resonance frequency of the square container , i.e. , precisely when d@xmath2-symmetric problem has a _ subcritical _ branch .
irregular bursts are also found , depending on parameters , but these are distinct from the chaotic states found by nagata @xcite far from threshold and present even in a square container .
in this article we have seen that there are many different mechanisms responsible for bursting in hydrodynamical systems .
the table below summarizes the different mechanisms described in terms of properties that are most relevant to hydrodynamics .
thus no single mechanism can be expected to provide a universal explanation for all observations .
the bursts found experimentally in taylor - couette flow ( cf . section 2.5 ) and large aspect ratio binary fluid convection ( cf . section 3.1 ) occur very close to the threshold of a primary instability and thus have the greatest potential for a successful dynamical systems interpretation of the type emphasized here .
we have seen , however , that even for fully developed turbulent boundary layers at very large reynolds numbers a dynamical systems approach can be profitable ( cf .
section 2.1 ) .
+ although nearly all of the mechanisms we have described rely on the presence of global bifurcations , there are important differences among them .
for example , the bursts in the wall region of a turbulent boundary layer described in section [ section : turb_bl ] are due to a ( structurally stable ) heteroclinic cycle connecting fixed points with _ finite _ amplitude ; such a cycle leads to bursts with a limited dynamic range .
in contrast , in the mechanism of section [ section : our_bursts.1 ] the dynamic range is unlimited .
moreover , the role of the fixed points is different : in the former the bursts are associated with the _ excursions _ between the fixed points while in the latter the bursts are associated _ with _ the fixed points . because of the asymptotic stability of the cycle the time between successive bursts in the turbulent boundary layer will increase without bound unless the stochastic pressure term is included ; such a stochastic term is not required in the mechanism of section [ section : our_bursts.1 ] .
in particular , in this mechanism the duration of the bursts remains finite despite the fact that they are associated with a heteroclinic connection .
this is because of the faster than exponential escape to `` infinity '' that is typical of this mechanism .
this is so also for the mechanism described in section [ section : hci ] .
both of these mechanisms involve global connections to infinity and hence are capable of describing bursts of arbitrarily large dynamic range .
the models of the subcritical transition to turbulence and various types of intermittency also produce bursts of finite duration but rely on global reinjection which produces bursts of bounded amplitude .
knobloch , e. & moehlis , j. , bursting mechanisms for hydrodynamical systems , _ pattern formation in continuous and coupled systems : a survey volume _ , eds . m. golubitsky , d. luss & s.h .
strogatz , springer , new york , pp . 157174 , 1999 .
waleffe , f. , how streamwise rolls and streaks self - sustain in a shear flow , _ self - sustaining mechanisms of wall turbulence _ , ed .
panton , computational mechanics publications , southampton , pp . 309332 , 1997 .
she , z.s . ,
large - scale dynamics and transition to turbulence in the two - dimensional kolmogorov flow , _ current trends in turbulence research _ , eds . h. branover , m. mond & y. unger , aiaa , washington , pp . 374396 , 1988 .
nicolaenko , b. & she , z.s . , symmetry - breaking homoclinic chaos in kolmogorov flows , _ nonlinear world , int .
workshop on nonlinear and turbulent processes in physics , kiev 1989 _ eds .
baryakhtar , v.m .
chernousenko , n.s .
erokhin , a.g .
sitenko & v.e .
zakharov , world scientific , singapore , pp . 602618 , 1990 .
nicolaenko , b. & she , z.s . , coherent structures , homoclinic cycles and vorticity explosions in navier - stokes flows , _ topological fluid mechanics _ , ed .
moffatt , cambridge university press , cambridge , pp . 265277 , 1990 .
nicolaenko , b. & she , z.s . , turbulent bursts , inertial sets and symmetry - breaking homoclinic cycles in periodic navier - stokes flows , _ turbulence in fluid flows : a dynamical systems approach _ , eds .
g.r . sell , c. foias & r. temam , springer - verlag , new york , pp . 123136 , 1993 .
armbruster , d. , nicolaenko , b. , smaoui , n. & chossat , p. , analyzing bifurcations in the kolmogorov flow equations , _ dynamics , bifurcation and symmetry _
p. chossat , kluwer , dordrecht , pp . 1134 , 1994 .
meshalkin , l.d .
& sinai , ia.g .
, investigation of the stability of a stationary solution of a system of equations for the plane movement of an incompressible viscous fluid , _ j. appl .
_ , * 25 * , pp .
17001705 , 1961 .
rinzel , j. , a formal classification of bursting mechanisms in excitable systems , _ proceedings of the international congress of mathematicians _
a.m. gleason , american mathematical society , providence , pp . 15781594 , 1987 .
rinzel , j. , a formal classification of bursting mechanisms in excitable systems , _ mathematical topics in population biology , morphogenesis and neurosciences _ ,
e. teramoto & m. yamaguti , lecture notes in biomathematics vol .
71 , springer - verlag , berlin , pp . 267281 , 1987 .
knobloch , e. , system symmetry breaking and shilnikov dynamics , _ pattern formation : symmetry methods and applications _ , eds . j. chadam , m. golubitsky , w. langford & b. wetton , american mathematical society , providence , pp . 271279 , 1996 .
pierce , r.d .
& knobloch , e. , spiral vortices in finite cylinders , _ ordered and turbulent patterns in taylor - couette flow
andereck & f. hayot , nato asi series b 297 , plenum , new york , pp . 8390 , 1992 .
armbruster , d. , codimension 2 bifurcation in binary convection with square symmetry , _ nonlinear evolution of spatio - temporal structures in dissipative continuous systems _ ,
busse & l. kramer , plenum press , new york , pp . 385398 , 1990 .
steindl a. & troger , h. , one and two - parameter bifurcations to divergence and flutter in the three - dimensional motions of a fluid conveying viscoelastic tube with d@xmath2 symmetry , _ nonlinear dynam .
_ , * 8 * , pp .
161178 , 1995 .