
Physica D: Nonlinear Phenomena [SJR: 1.049] [HI: 102] [3 followers] Follow Hybrid journal (It can contain Open Access articles) ISSN (Print) 01672789 Published by Elsevier [3043 journals] 
 A nonperturbative analytic expression of signal amplification factor in
stochastic resonance Authors: Asish Kumar Dhara
Pages: 1  44
Abstract: Publication date: 1 April 2017
Source:Physica D: Nonlinear Phenomena, Volume 344
Author(s): Asish Kumar Dhara
We put forward a nonperturbative scheme to calculate the response of an overdamped bistable system driven by a Gaussian white noise and perturbed by a weak monochromatic force (signal) analytically. The formalism takes into account infinite number of perturbation terms of a perturbation series with amplitude of the signal as an expansion parameter. The contributions of infinite number of relaxation modes of the stochastic dynamics to the response are also taken into account in this formalism. A closed form analytic expression of the response is obtained. Only the knowledge of the first nontrivial eigenvalue and the lowest eigenfunction of the unperturbed Fokker–Planck operator are needed to evaluate the response. The response calculated from the derived analytic expression matches fairly well with the numerical results.
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2016.11.002
Issue No: Vol. 344 (2017)
 Authors: Asish Kumar Dhara
 Elastic and inelastic collisions of swarms
 Authors: Dieter Armbruster; Stephan Martin; Andrea Thatcher
Pages: 45  57
Abstract: Publication date: 1 April 2017
Source:Physica D: Nonlinear Phenomena, Volume 344
Author(s): Dieter Armbruster, Stephan Martin, Andrea Thatcher
Scattering interactions of swarms in potentials that are generated by an attraction–repulsion model are studied. In free space, swarms in this model form a welldefined steady state describing the translation of a stable formation of the particles whose shape depends on the interaction potential. Thus, the collision between a swarm and a boundary or between two swarms can be treated as (quasi)particle scattering. Such scattering experiments result in internal excitations of the swarm or in bound states, respectively. In addition, varying a parameter linked to the relative importance of damping and potential forces drives transitions between elastic and inelastic scattering of the particles. By tracking the swarm’s center of mass, a refraction rule is derived via simulations relating the incoming and outgoing directions of a swarm hitting the wall. Iterating the map derived from the refraction law allows us to predict and understand the dynamics and bifurcations of swarms in square boxes and in channels.
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2016.11.008
Issue No: Vol. 344 (2017)
 Authors: Dieter Armbruster; Stephan Martin; Andrea Thatcher
 Swarming in bounded domains
 Authors: Dieter Armbruster; Sébastien Motsch; Andrea Thatcher
Pages: 58  67
Abstract: Publication date: 1 April 2017
Source:Physica D: Nonlinear Phenomena, Volume 344
Author(s): Dieter Armbruster, Sébastien Motsch, Andrea Thatcher
The Vicsek model is a prototype for the emergence of collective motion. In free space, it is characterized by a swarm of particles all moving in the same direction. Since this dynamic does not include attraction among particles, the swarm, while aligning in velocity space, has no spatial coherence. Adding specular reflection at the boundaries generates global spatial coherence of the swarms while maintaining its velocity alignment. We investigate numerically how the geometry of the domain influences the Vicsek model using three type of geometry: a channel, a disk and a rectangle. Varying the parameters of the Vicsek model (e.g. noise levels and influence horizons), we discuss the mechanisms that generate spatial coherence and show how they create new dynamical solutions of the swarming motions in these geometries. Several observables are introduced to characterize the simulated patterns (e.g. mass profile, center of mass, connectivity of the swarm).
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2016.11.009
Issue No: Vol. 344 (2017)
 Authors: Dieter Armbruster; Sébastien Motsch; Andrea Thatcher
 Selfassembly of a filament by curvatureinducing proteins
 Authors: James Kwiecinski; S. Jonathan Chapman; Alain Goriely
Pages: 68  80
Abstract: Publication date: 1 April 2017
Source:Physica D: Nonlinear Phenomena, Volume 344
Author(s): James Kwiecinski, S. Jonathan Chapman, Alain Goriely
We explore a simplified macroscopic model of membrane shaping by means of curvaturesensing BAR proteins. Equations describing the interplay between the shape of a freely floating filament in a fluid and the adhesion kinetics of proteins are derived from mechanical principles. The constant curvature solutions that arise from this system are studied using weakly nonlinear analysis. We show that the stability of the filament’s shape is completely characterized by the parameters associated with protein recruitment and establish that in the bistable regime, proteins aggregate on the filament forming regions of high and low curvatures. This pattern formation is then followed by phasecoarsening that resolves on a timescale dependent on protein diffusion and drift across the filament, which contend to smooth and maintain the pattern respectively. The model is generalized for multiple species of BAR proteins and we show that the stability of the assembled shape is determined by a competition between proteins attaching on opposing sides.
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2016.12.001
Issue No: Vol. 344 (2017)
 Authors: James Kwiecinski; S. Jonathan Chapman; Alain Goriely
 Corrigendum to “Fluctuations induced transition of localization of
granular objects caused by degrees of crowding” [Physica D 336 (2016)
3946] Authors: Soutaro Oda; Yoshitsugu Kubo; ChwenYang Shew; Kenichi Yoshikawa
First page: 81
Abstract: Publication date: 1 April 2017
Source:Physica D: Nonlinear Phenomena, Volume 344
Author(s): Soutaro Oda, Yoshitsugu Kubo, ChwenYang Shew, Kenichi Yoshikawa
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2017.02.006
Issue No: Vol. 344 (2017)
 Authors: Soutaro Oda; Yoshitsugu Kubo; ChwenYang Shew; Kenichi Yoshikawa
 Lowdimensional reducedorder models for statistical response and
uncertainty quantification: Barotropic turbulence with topography Authors: Di Qi; Andrew J. Majda
Pages: 7  27
Abstract: Publication date: 15 March 2017
Source:Physica D: Nonlinear Phenomena, Volume 343
Author(s): Di Qi, Andrew J. Majda
A lowdimensional reducedorder statistical closure model is developed for quantifying the uncertainty to changes in forcing in a barotropic turbulent system with topography involving interactions between smallscale motions and a largescale mean flow. Imperfect model sensitivity is improved through a recent mathematical strategy for calibrating model errors in a training phase, where information theory and linear statistical response theory are combined in a systematic fashion to achieve the optimal model parameters. Statistical theories about a Gaussian invariant measure and the exact statistical energy equations are also developed for the truncated barotropic equations that can be used to improve the imperfect model prediction skill. A stringent paradigm model of 57 degrees of freedom is used to display the feasibility of the reducedorder methods. This simple model creates largescale zonal mean flow shifting directions from westward to eastward jets with an abrupt change in amplitude when perturbations are applied, and prototype blocked and unblocked patterns can be generated in this simple model similar to the real natural system. Principal statistical responses in mean and variance can be captured by the reducedorder models with desirable accuracy and efficiency with only 3 resolved modes. An even more challenging regime with nonGaussian equilibrium statistics using the fluctuation equations is also tested in the reducedorder models with accurate prediction using the first 5 resolved modes. These reducedorder models also show potential for uncertainty quantification and prediction in more complex realistic geophysical turbulent dynamical systems.
PubDate: 20170209T06:01:16Z
DOI: 10.1016/j.physd.2016.11.006
Issue No: Vol. 343 (2017)
 Authors: Di Qi; Andrew J. Majda
 Isolating blocks as computational tools in the circular restricted
threebody problem Authors: Rodney L. Anderson; Robert W. Easton; Martin W. Lo
Pages: 38  50
Abstract: Publication date: 15 March 2017
Source:Physica D: Nonlinear Phenomena, Volume 343
Author(s): Rodney L. Anderson, Robert W. Easton, Martin W. Lo
Isolating blocks may be used as computational tools to search for the invariant manifolds of orbits and hyperbolic invariant sets associated with libration points while also giving additional insight into the dynamics of the flow in these regions. We use isolating blocks to investigate the dynamics of objects entering the Earth–Moon system in the circular restricted threebody problem with energies close to the energy of the L 2 libration point. Specifically, the stable and unstable manifolds of Lyapunov orbits and the hyperbolic invariant set around the libration points are obtained by numerically computing the way orbits exit from an isolating block in combination with a bisection method. Invariant spheres of solutions in the spatial problem may then be located using the resulting manifolds.
PubDate: 20170209T06:01:16Z
DOI: 10.1016/j.physd.2016.10.004
Issue No: Vol. 343 (2017)
 Authors: Rodney L. Anderson; Robert W. Easton; Martin W. Lo
 Finitetime thin film rupture driven by modified evaporative loss
 Authors: Hangjie Ji; Thomas P. Witelski
Pages: 1  15
Abstract: Publication date: 1 March 2017
Source:Physica D: Nonlinear Phenomena, Volume 342
Author(s): Hangjie Ji, Thomas P. Witelski
Rupture is a nonlinear instability resulting in a finitetime singularity as a film layer approaches zero thickness at a point. We study the dynamics of rupture in a generalized mathematical model of thin films of viscous fluids with modified evaporative effects. The governing lubrication model is a fourthorder nonlinear parabolic partial differential equation with a nonconservative loss term. Several different types of finitetime singularities are observed due to balances between conservative and nonconservative terms. Nonselfsimilar behavior and two classes of selfsimilar rupture solutions are analyzed and validated against high resolution PDE simulations.
PubDate: 20170114T04:47:15Z
DOI: 10.1016/j.physd.2016.10.002
Issue No: Vol. 342 (2017)
 Authors: Hangjie Ji; Thomas P. Witelski
 Stability on timedependent domains: convective and dilution effects
 Authors: R. Krechetnikov; E. Knobloch
Pages: 16  23
Abstract: Publication date: 1 March 2017
Source:Physica D: Nonlinear Phenomena, Volume 342
Author(s): R. Krechetnikov, E. Knobloch
We explore nearcritical behavior of spatially extended systems on timedependent spatial domains with convective and dilution effects due to domain flow. As a paradigm, we use the Swift–Hohenberg equation, which is the simplest nonlinear model with a nonzero critical wavenumber, to study dynamic pattern formation on timedependent domains. A universal amplitude equation governing weakly nonlinear evolution of patterns on timedependent domains is derived and proves to be a generalization of the standard Ginzburg–Landau equation. Its key solutions identified here demonstrate a substantial variety–spatially periodic states with a timedependent wavenumber, steady spatially nonperiodic states, and pulsetrain solutions–in contrast to extended systems on timefixed domains. The effects of domain flow, such as bifurcation delay due to domain growth and destabilization due to oscillatory domain flow, on the Eckhaus instability responsible for phase slips in spatially periodic states are analyzed with the help of both local and global stability analyses. A nonlinear phase equation describing the approach to a phaseslip event is derived. Detailed analysis of a phase slip using multiple time scale methods demonstrates different mechanisms governing the wavelength changing process at different stages.
PubDate: 20170114T04:47:15Z
DOI: 10.1016/j.physd.2016.10.003
Issue No: Vol. 342 (2017)
 Authors: R. Krechetnikov; E. Knobloch
 Stability of the phase motion in racetrack microtons
 Authors: Yu.A. Kubyshin; O. Larreal; R. RamírezRos; T.M. Seara
Abstract: Publication date: Available online 18 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Yu.A. Kubyshin, O. Larreal, R. RamírezRos, T.M. Seara
We model the phase oscillations of electrons in racetrack microtrons by means of an area preserving map with a fixed point at the origin, which represents the synchronous trajectory of a reference particle in the beam. We study the nonlinear stability of the origin in terms of the synchronous phase —the phase of the synchronous particle at the injection. We estimate the size and shape of the stability domain around the origin, whose main connected component is enclosed by an invariant curve. We describe the evolution of the stability domain as the synchronous phase varies. We also clarify the role of the stable and unstable invariant curves of some hyperbolic (fixed or periodic) points.
PubDate: 20170320T19:31:07Z
DOI: 10.1016/j.physd.2017.03.001
 Authors: Yu.A. Kubyshin; O. Larreal; R. RamírezRos; T.M. Seara
 Markovian properties of velocity increments in boundary layer turbulence
 Authors: Murat Tutkun
Abstract: Publication date: Available online 16 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Murat Tutkun
Markovian properties of the turbulent velocity increments in a flat plate boundary layer at Re θ = 19100 are investigated using hotwire anemometry measurements of the streamwise velocity component in a wind tunnel. Increments of the longitudinal velocities at different wallnormal positions show that the flow exhibits Markovian properties when the separation between different scales, or the MarkovEinstein coherence length, is on the order of the Taylor microscale, λ . The results indicate that Markovian nature of turbulence evolves across the boundary layer showing certain characteristics pertaining to the distance to the wall. The connection between the Markovian properties of turbulent boundary layer and existence of the spectral gap is explored. Markovianity of the process is also discussed in relation to the nonlocal nonlinear versus local nonlinear transfer of energy, triadic interactions and dissipation.
PubDate: 20170320T19:31:07Z
DOI: 10.1016/j.physd.2017.03.002
 Authors: Murat Tutkun
 Targeted energy transfer in laminar vortexinduced vibration of a sprung
cylinder with a nonlinear dissipative rotator Authors: Antoine Blanchard; Lawrence A. Bergman; Alexander F. Vakakis
Abstract: Publication date: Available online 16 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Antoine Blanchard, Lawrence A. Bergman, Alexander F. Vakakis
We computationally investigate the dynamics of a linearlysprung circular cylinder immersed in an incompressible flow and undergoing transverse vortexinduced vibration (VIV), to which is attached a rotational nonlinear energy sink (NES) consisting of a mass that freely rotates at constant radius about the cylinder axis, and whose motion is restrained by a rotational linear viscous damper. The inertial coupling between the rotational motion of the attached mass and the rectilinear motion of the cylinder is “essentially nonlinear”, which, in conjunction with dissipation, allows for oneway, nearly irreversible targeted energy transfer (TET) from the oscillating cylinder to the nonlinear dissipative attachment. At the intermediate Reynolds number R e = 100 , the NESequipped sprung cylinder undergoes repetitive cycles of slowly decaying oscillations punctuated by intervals of chaotic instabilities. During the slowly decaying portion of each cycle, the dynamics of the cylinder is regular and, for large enough values of the ratio ε of the NES mass to the total mass (i.e., NES mass plus cylinder mass), can lead to significant vortex street elongation with partial stabilization of the wake. As ε approaches zero, no such vortex elongation is observed and the wake patterns appear similar to that for a sprung cylinder with no NES. We apply proper orthogonal decomposition (POD) to the velocity flow field during a slowly decaying portion of the solution and show that, in situations where vortex elongation occurs, the NES, though not in direct contact with the surrounding fluid, has a drastic effect on the underlying flow structures, imparting significant and continuous passive redistribution of energy among POD modes. We construct a PODbased reducedorder model for the lift coefficient to characterize energy transactions between the fluid and the cylinder throughout the slowly decaying cycle. We introduce a quantitative signed measure of the work done by the fluid on the cylinder over one quasiperiod of the slowly decaying response and find that vortex elongation is associated with a sign change of that measure, indicating that a reversal of the direction of energy transfer, with the cylinder “leaking energy back” to the flow, is responsible for partial stabilization and elongation of the wake. We interpret these findings in terms of the spatial structure and energy distribution of the POD modes, and relate them to the mechanism of transient resonance capture into a slow invariant manifold of the fluid–structure interaction dynamics.
PubDate: 20170320T19:31:07Z
DOI: 10.1016/j.physd.2017.03.003
 Authors: Antoine Blanchard; Lawrence A. Bergman; Alexander F. Vakakis
 Fluctuations, response, and resonances in a simple atmospheric model
 Authors: Andrey Gritsun; Valerio Lucarini
Abstract: Publication date: Available online 14 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Andrey Gritsun, Valerio Lucarini
We study the response of a simple quasigeostrophic barotropic model of the atmosphere to various classes of perturbations affecting its forcing and its dissipation using the formalism of the Ruelle response theory. We investigate the geometry of such perturbations by constructing the covariant Lyapunov vectors of the unperturbed system and discover in one specific case–orographic forcing–a substantial projection of the forcing onto the stable directions of the flow. This results into a resonant response shaped as a Rossbylike wave that has no resemblance to the unforced variability in the same range of spatial and temporal scales. Such a climatic surprise corresponds to a violation of the fluctuationdissipation theorem, in agreement with the basic tenets of nonequilibrium statistical mechanics. The resonance can be attributed to a specific group of rarely visited unstable periodic orbits of the unperturbed system. Our results reinforce the idea of using basic methods of nonequilibrium statistical mechanics and highdimensional chaotic dynamical systems to approach the problem of understanding climate dynamics.
PubDate: 20170320T19:31:07Z
DOI: 10.1016/j.physd.2017.02.015
 Authors: Andrey Gritsun; Valerio Lucarini
 Standing and travelling waves in a spherical brain model: The Nunez model
revisited Authors: S. Visser; R. Nicks; O. Faugeras; S. Coombes
Abstract: Publication date: Available online 8 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): S. Visser, R. Nicks, O. Faugeras, S. Coombes
The Nunez model for the generation of electroencephalogram (EEG) signals is naturally described as a neural field model on a sphere with spacedependent delays. For simplicity, dynamical realisations of this model either as a damped wave equation or an integrodifferential equation, have typically been studied in idealised one dimensional or planar settings. Here we revisit the original Nunez model to specifically address the role of spherical topology on spatiotemporal pattern generation. We do this using a mixture of Turing instability analysis, symmetric bifurcation theory, center manifold reduction and direct simulations with a bespoke numerical scheme. In particular we examine standing and travelling wave solutions using normal form computation of primary and secondary bifurcations from a steady state. Interestingly, we observe spatiotemporal patterns which have counterparts seen in the EEG patterns of both epileptic and schizophrenic brain conditions.
PubDate: 20170312T21:31:41Z
DOI: 10.1016/j.physd.2017.02.017
 Authors: S. Visser; R. Nicks; O. Faugeras; S. Coombes
 Controlling coexisting attractors of an impacting system via linear
augmentation Authors: Yang Liu; Joseph Páez Chávez
Abstract: Publication date: Available online 8 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Yang Liu, Joseph Páez Chávez
This paper studies the control of coexisting attractors in an impacting system via a recently developed control law based on linear augmentation. Special attention is given to two control issues in the framework of multistable engineering systems, namely, the switching between coexisting attractors without altering the system’s main parameters and the avoidance of grazinginduced chaotic responses. The effectiveness of the proposed control scheme is confirmed numerically for the case of a periodically excited, soft impact oscillator. Our analysis shows how pathfollowing techniques for nonsmooth systems can be used in order to determine the optimal control parameters in terms of energy expenditure due to the control signal and transient behavior of the control error, which can be applied to a broad range of engineering problems.
PubDate: 20170312T21:31:41Z
DOI: 10.1016/j.physd.2017.02.018
 Authors: Yang Liu; Joseph Páez Chávez
 Construction of Darboux coordinates and PoincaréBirkhoff normal forms in
noncanonical Hamiltonian systems Authors: Andrej Junginger; Jörg Main; Günter Wunner
Abstract: Publication date: Available online 7 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Andrej Junginger, Jörg Main, Günter Wunner
We demonstrate a general method to construct Darboux coordinates via normal form expansions in noncanonical Hamiltonian system obtained from e.g. a variational approach to quantum systems. The procedure serves as a tool to naturally extract canonical coordinates out of the variational parameters and at the same time to transform the energy functional into its PoincaréBirkhoff normal form. The method is general in the sense that it is applicable for arbitrary degrees of freedom, in arbitrary orders of the local expansion, and it is independent of the precise form of the Hamilton operator. The method presented allows for the general and systematic investigation of quantum systems in the vicinity of fixed points, which e.g. correspond to ground, excited or transition states. Moreover, it directly allows to calculate classical and quantum reaction rates by applying transition state theory.
PubDate: 20170312T21:31:41Z
DOI: 10.1016/j.physd.2017.02.014
 Authors: Andrej Junginger; Jörg Main; Günter Wunner
 Assessment of the effects of azimuthal mode number perturbations upon the
implosion processes of fluids in cylinders Authors: Michael Lindstrom
Abstract: Publication date: Available online 1 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Michael Lindstrom
Fluid instabilities arise in a variety of contexts and are often unwanted results of engineering imperfections. In one particular model for a magnetized target fusion reactor, a pressure wave is propagated in a cylindrical annulus comprised of a dense fluid before impinging upon a plasma and imploding it. Part of the success of the apparatus is a function of how axiallysymmetric the final pressure pulse is upon impacting the plasma. We study a simple model for the implosion of the system to study how imperfections in the pressure imparted on the outer circumference grow due to geometric focusing. Our methodology entails linearizing the compressible Euler equations for mass and momentum conservation about a cylindrically symmetric problem and analyzing the perturbed profiles at different mode numbers. The linearized system gives rise to singular shocks and through analyzing the perturbation profiles at various times, we infer that high mode numbers are dampened through the propagation. We also study the Linear KleinGordon equation in the context of stability of linear cylindrical wave formation whereby highly oscillatory, bounded behaviour is observed in a far field solution.
PubDate: 20170307T21:14:06Z
DOI: 10.1016/j.physd.2017.02.012
 Authors: Michael Lindstrom
 KdV cnoidal waves in a traffic flow model with periodic boundaries
 Authors: L.L. Hattam
Abstract: Publication date: Available online 1 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): L.L. Hattam
An optimalvelocity (OV) model describes car motion on a single lane road. In particular, near to the boundary signifying the onset of traffic jams, this model reduces to a perturbed Kortewegde Vries (KdV) equation using asymptotic analysis. Previously, the KdV soliton solution has then been found and compared to numerical results (see Muramatsu and Nagatani (1999)). Here, we instead apply modulation theory to this perturbed KdV equation to obtain at leading order, the modulated cnoidal wave solution. At the next order, the Whitham equations are derived, which have been modified due to the equation perturbation terms. Next, from this modulation system, a family of spatially periodic cnoidal waves are identified that characterise vehicle headway distance. Then, for this set of solutions, we establish the relationship between the wave speed, the modulation term and the driver sensitivity. This analysis is confirmed with comparisons to numerical solutions of the OV model. As well, the longtime behaviour of these solutions is investigated.
PubDate: 20170307T21:14:06Z
DOI: 10.1016/j.physd.2017.02.010
 Authors: L.L. Hattam
 Controlling roughening processes in the stochastic KuramotoSivashinsky
equation Authors: S.N. Gomes; S. Kalliadasis; D.T. Papageorgiou; G.A. Pavliotis; M. Pradas
Abstract: Publication date: Available online 1 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): S.N. Gomes, S. Kalliadasis, D.T. Papageorgiou, G.A. Pavliotis, M. Pradas
We present a novel control methodology to control the roughening processes of semilinear parabolic stochastic partial differential equations in one dimension, which we exemplify with the stochastic KuramotoSivashinsky equation. The original equation is split into a linear stochastic and a nonlinear deterministic equation so that we can apply linear feedback control methods. Our control strategy is then based on two steps: first, stabilize the zero solution of the deterministic part and, second, control the roughness of the stochastic linear equation. We consider both periodic controls and point actuated ones, observing in all cases that the second moment of the solution evolves in time according to a powerlaw until it saturates at the desired controlled value.
PubDate: 20170307T21:14:06Z
DOI: 10.1016/j.physd.2017.02.011
 Authors: S.N. Gomes; S. Kalliadasis; D.T. Papageorgiou; G.A. Pavliotis; M. Pradas
 A novel route to a Hopf bifurcation scenario in switched systems with
deadzone Authors: P. Kowalczyk
Abstract: Publication date: Available online 1 March 2017
Source:Physica D: Nonlinear Phenomena
Author(s): P. Kowalczyk
Planar switched systems with deadzone are analyzed. In particular, we consider the effects of a perturbation which is applied to a linear control law and, due to the perturbation, the control changes from purely positional to positionvelocity control. This type of a perturbation leads to a novel Hopflike discontinuity induced bifurcation. We show that this bifurcation leads to the creation of a small scale limit cycle attractor, which scales as the square root of the bifurcation parameter. We then investigate numerically a planar switched system with a positional feedback law, deadzone and time delay in the switching function. Using the same parameter values as for the switched system without time delay in the switching function, we show a Hopflike bifurcation scenario which exhibits a qualitative and a quantitative agreement with the scenario analyzed for the nondelayed system.
PubDate: 20170307T21:14:06Z
DOI: 10.1016/j.physd.2017.02.007
 Authors: P. Kowalczyk
 Hopf and Homoclinic Bifurcations on the sliding vector field of switching
systems in R3: A case study in power electronics Authors: Rony Cristiano; Tiago Carvalho; Durval J. Tonon; Daniel J. Pagano
Abstract: Publication date: Available online 24 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Rony Cristiano, Tiago Carvalho, Durval J. Tonon, Daniel J. Pagano
In this paper, Hopf and homoclinic bifurcations that occur in the sliding vector field of switching systems in R 3 are studied. In particular, a dcdc boost converter with sliding mode control and washout filter is analyzed. This device is modelled as a threedimensional Filippov system, characterized by the existence of sliding movement and restricted to the switching manifold. The operating point of the converter is a stable pseudoequilibrium and it undergoes a subcritical Hopf bifurcation. Such a bifurcation occurs in the sliding vector field and creates, in this field, an unstable limit cycle. The limit cycle is connected to the switching manifold and disappears when it touches the visibleinvisible twofold point, resulting in an homoclinic loop which itself closes in this twofold point. The study of these dynamic phenomena that can be found in different power electronic circuits controlled by sliding mode control strategies are relevant from the viewpoint of the global stability and robustness of the control design.
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2017.02.005
 Authors: Rony Cristiano; Tiago Carvalho; Durval J. Tonon; Daniel J. Pagano
 Axisymmetric pulse train solutions in narrowgap spherical Couette flow
 Authors: Adam Child; Rainer Hollerbach; Evy Kersalé
Abstract: Publication date: Available online 24 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Adam Child, Rainer Hollerbach, Evy Kersalé
We numerically compute the flow induced in a spherical shell by fixing the outer sphere and rotating the inner one. The aspect ratio ϵ = ( r o − r i ) / r i is set at 0.04 and 0.02, and in each case the Reynolds number measuring the inner sphere’s rotation rate is increased to ∼ 10% beyond the first bifurcation from the basic state flow. For ϵ = 0.04 the initial bifurcations are the same as in previous numerical work at ϵ = 0.154 , and result in steady one and twovortex states. Further bifurcations yield travelling wave solutions similar to previous analytic results valid in the ϵ → 0 limit. For ϵ = 0.02 the steady onevortex state no longer exists, and the first bifurcation is directly to these travelling wave solutions, consisting of pulse trains of Taylor vortices travelling toward the equator from both hemispheres, and annihilating there in distinct phaseslip events. We explore these timedependent solutions in detail, and find that they can be both equatorially symmetric and asymmetric, as well as periodic or quasiperiodic in time.
PubDate: 20170301T05:09:55Z
DOI: 10.1016/j.physd.2017.02.009
 Authors: Adam Child; Rainer Hollerbach; Evy Kersalé
 Assimilating Eulerian and Lagrangian data in trafficflow models
 Authors: Chao Xia; Courtney Cochrane; Joseph DeGuire; Gaoyang Fan; Emma Holmes; Melissa McGuirl; Patrick Murphy; Jenna Palmer; Paul Carter; Laura Slivinski; Björn Sandstede
Abstract: Publication date: Available online 21 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Chao Xia, Courtney Cochrane, Joseph DeGuire, Gaoyang Fan, Emma Holmes, Melissa McGuirl, Patrick Murphy, Jenna Palmer, Paul Carter, Laura Slivinski, Björn Sandstede
Data assimilation of traffic flow remains a challenging problem. One difficulty is that data come from different sources ranging from stationary sensors and camera data to GPS and cell phone data from moving cars. Sensors and cameras give information about traffic density, while GPS data provide information about the positions and velocities of individual cars. Previous methods for assimilating Lagrangian data collected from individual cars relied on specific properties of the underlying computational model or its reformulation in Lagrangian coordinates. These approaches make it hard to assimilate both Eulerian density and Lagrangian positional data simultaneously. In this paper, we propose an alternative approach that allows us to assimilate both Eulerian and Lagrangian data. We show that the proposed algorithm is accurate and works well in different traffic scenarios and regardless of whether ensemble Kalman or particle filters are used. We also show that the algorithm is capable of estimating parameters and assimilating real traffic observations and synthetic observations obtained from microscopic models.
PubDate: 20170221T09:49:54Z
DOI: 10.1016/j.physd.2017.02.004
 Authors: Chao Xia; Courtney Cochrane; Joseph DeGuire; Gaoyang Fan; Emma Holmes; Melissa McGuirl; Patrick Murphy; Jenna Palmer; Paul Carter; Laura Slivinski; Björn Sandstede
 Complex predator invasion waves in a HollingTanner model with nonlocal
prey interaction Authors: A. Bayliss; V.A. Volpert
Abstract: Publication date: Available online 20 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): A. Bayliss, V.A. Volpert
We consider predator invasions for the nonlocal HollingTanner model. Predators are introduced in a small region adjacent to an extensive predatorfree region. In its simplest form an invasion front propagates into the predatorfree region with a predatorprey coexistence state displacing the predatorfree state. However, patterns may form in the wake of the invasion front due to instability of the coexistence state. The coexistence state can be subject to either oscillatory or cellular instability, depending on parameters. Furthermore, the oscillatory instability can be either at zero wave number or finite wave number. In addition, the (unstable) predatorfree state can be subject to additional cellular instabilities when the extent of the nonlocality is sufficiently large. We perform numerical simulations that demonstrate that the invasion wave may have a complex structure in which different spatial regions exhibit qualitatively different behaviors. These regions are separated by relatively narrow transition regions that we refer to as fronts. We also derive analytic approximations for the speeds of the fronts and find qualitative and quantitative agreement with the results of computations.
PubDate: 20170221T09:49:54Z
DOI: 10.1016/j.physd.2017.02.003
 Authors: A. Bayliss; V.A. Volpert
 Elementary solutions for a model Boltzmann equation in one dimension and
the connection to grossly determined solutions Authors: Thomas E. Carty
Abstract: Publication date: Available online 20 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Thomas E. Carty
The Fouriertransformed version of the time dependent slipflow model Boltzmann equation associated with the linearized BGK model is solved in order to determine the solution’s asymptotics. The ultimate goal of this paper is to demonstrate that there exists a robust set of solutions to this model Boltzmann equation that possess a special property that was conjectured by Truesdell and Muncaster: that solutions decay to a subclass of the solution set uniquely determined by the initial mass density of the gas called the grossly determined solutions. First we determine the spectrum and eigendistributions of the associated homogeneous equation. Then, using Case’s method of elementary solutions, we find analytic timedependent solutions to the model Boltzmann equation for initial data with a specialized compact support condition under the Fourier transform. In doing so, we show that the spectrum separates the solutions into two distinct parts: one that behaves as a set of transient solutions and the other limiting to a stable subclass of solutions. Thus, we demonstrate that for gas flows with this specialized initial density condition, in time all gas flows for the one dimensional model Boltzmann equation act as grossly determined solutions.
PubDate: 20170221T09:49:54Z
DOI: 10.1016/j.physd.2017.02.008
 Authors: Thomas E. Carty
 Bright and dark solitons in the unidirectional long wave limit for the
energy transfer on anharmonic crystal lattices Authors: Luis A. CisnerosAke; José F. Solano Peláez
Abstract: Publication date: Available online 9 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Luis A. CisnerosAke, José F. Solano Peláez
The problem of energy transportation along a cubic anharmonic crystal lattice, in the unidirectional long wave limit, is considered. A detailed process, in the discrete lattice equations, shows that unidirectional stable propagating waves for the continuum limit produces a coupled system between a nonlinear Schrödinger (NLS) equation and the KortewegdeVries (KdV) equation. The traveling wave formalism provides a diversity of exact solutions ranging from the classical Davydov’s soliton (subsonic and supersonic) of the first and second kind to a class consisting in the coupling between the KdV soliton and dark solitons containing the typical ones (similar to the darkgray soliton in the standard defocusing NLS) and a new kind in the form of a twohump dark soliton. This family of exact solutions are numerically tested, by means of the pseudo spectral method, in our NLSKdV system.
PubDate: 20170215T09:34:42Z
DOI: 10.1016/j.physd.2017.02.001
 Authors: Luis A. CisnerosAke; José F. Solano Peláez
 The LyapunovKrasovskii theorem and a sufficient criterion for local
stability of isochronal synchronization in networks of delaycoupled
oscillators Authors: J.M.V. Grzybowski; E.E.N. Macau; T. Yoneyama
Abstract: Publication date: Available online 9 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): J.M.V. Grzybowski, E.E.N. Macau, T. Yoneyama
This paper presents a selfcontained framework for the stability assessment of isochronal synchronization in networks of chaotic and limitcycle oscillators. The results were based on the LyapunovKrasovskii theorem and they establish a sufficient condition for local synchronization stability of as a function of the system and network parameters. With this in mind, a network of mutually delaycoupled oscillators subject to direct selfcoupling is considered and then the resulting error equations are blockdiagonalized for the purpose of studying their stability. These error equations are evaluated by means of analytical stability results derived from the LyapunovKrasovskii theorem. The proposed approach is shown to be a feasible option for the investigation of local stability of isochronal synchronization for a variety of oscillators coupled through linear functions of the state variables under a given undirected graph structure. This ultimately permits the systematic identification of stability regions within the highdimensionality of the network parameter space. Examples of applications of the results to a number of networks of delaycoupled chaotic and limitcycle oscillators are provided, such as Lorenz, Rössler, Cubic Chua’s circuit, Van der Pol oscillator and the HindmarshRose neuron.
PubDate: 20170215T09:34:42Z
DOI: 10.1016/j.physd.2017.01.005
 Authors: J.M.V. Grzybowski; E.E.N. Macau; T. Yoneyama
 Wave fronts and cascades of soliton interactions in the periodic two
dimensional Volterra system Authors: Rhys Bury; Alexander V. Mikhailov; Jing Ping Wang
Abstract: Publication date: Available online 6 February 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Rhys Bury, Alexander V. Mikhailov, Jing Ping Wang
In the paper we develop the dressing method for the solution of the twodimensional periodic Volterra system with a period N . We derive soliton solutions of arbitrary rank k and give a full classification of rank 1 solutions. We have found a new class of exact solutions corresponding to wave fronts which represent smooth interfaces between two nonlinear periodic waves or a periodic wave and a trivial (zero) solution. The wave fronts are nonstationary and they propagate with a constant average velocity. The system also has soliton solutions similar to breathers, which resembles soliton webs in the KP theory. We associate the classification of soliton solutions with the Schubert decomposition of the Grassmannians Gr R ( k , N ) and Gr C ( k , N ) .
PubDate: 20170209T06:01:16Z
DOI: 10.1016/j.physd.2017.01.003
 Authors: Rhys Bury; Alexander V. Mikhailov; Jing Ping Wang
 The stability spectrum for elliptic solutions to the focusing NLS equation
 Authors: Bernard Deconinck; Benjamin L. Segal
Abstract: Publication date: Available online 23 January 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Bernard Deconinck, Benjamin L. Segal
We present an analysis of the stability spectrum of all stationary elliptictype solutions to the focusing Nonlinear Schrödinger equation (NLS). An analytical expression for the spectrum is given. From this expression, various quantitative and qualitative results about the spectrum are derived. Specifically, the solution parameter space is shown to be split into four regions of distinct qualitative behavior of the spectrum. Additional results on the stability of solutions with respect to perturbations of an integer multiple of the period are given, as well as a procedure for approximating the greatest real part of the spectrum.
PubDate: 20170128T05:15:28Z
DOI: 10.1016/j.physd.2017.01.004
 Authors: Bernard Deconinck; Benjamin L. Segal
 Reducedspace Gaussian Process Regression for datadriven probabilistic
forecast of chaotic dynamical systems Authors: Zhong Yi Wan; Themistoklis P. Sapsis
Abstract: Publication date: Available online 5 January 2017
Source:Physica D: Nonlinear Phenomena
Author(s): Zhong Yi Wan, Themistoklis P. Sapsis
We formulate a reducedorder strategy for efficiently forecasting complex highdimensional dynamical systems entirely based on data streams. The first step of our method involves reconstructing the dynamics in a reducedorder subspace of choice using Gaussian Process Regression (GPR). GPR simultaneously allows for reconstruction of the vector field and more importantly, estimation of local uncertainty. The latter is due to (i) local interpolation error and (ii) truncation of the highdimensional phase space. This uncertainty component can be analytically quantified in terms of the GPR hyperparameters. In the second step we formulate stochastic models that explicitly take into account the reconstructed dynamics and their uncertainty. For regions of the attractor which are not sufficiently sampled for our GPR framework to be effective, an adaptive blended scheme is formulated to enforce correct statistical steady state properties, matching those of the real data. We examine the effectiveness of the proposed method to complex systems including the Lorenz 96, the KuramotoSivashinsky, as well as a prototype climate model. We also study the performance of the proposed approach as the intrinsic dimensionality of the system attractor increases in highly turbulent regimes.
PubDate: 20170114T04:47:15Z
DOI: 10.1016/j.physd.2016.12.005
 Authors: Zhong Yi Wan; Themistoklis P. Sapsis
 Characterizing complex networks through statistics of Möbius
transformations Authors: Vladimir Jaćimović; Aladin Crnkić
Abstract: Publication date: Available online 31 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Vladimir Jaćimović, Aladin Crnkić
It is wellknown now that dynamics of large populations of globally (alltoall) coupled oscillators can be reduced to lowdimensional submanifolds (WS transformation and OA ansatz). Marvel et al. (2009) described an intriguing algebraic structure standing behind this reduction: oscillators evolve by the action of the group of Möbius transformations. Of course, dynamics in complex networks of coupled oscillators is highly complex and not reducible. Still, closer look unveils that even in complex networks some (possibly overlapping) groups of oscillators evolve by Möbius transformations. In this paper we study properties of the network by identifying Möbius transformations in the dynamics of oscillators. This enables us to introduce some new (statistical) concepts that characterize the network. In particular, the notion of coherence of the network (or subnetwork) is proposed. This conceptual approach is meaningful for the broad class of networks, including those with timedelayed, noisy or mixed interactions. In this paper several simple (random) graphs are studied illustrating the meaning of the concepts introduced in the paper.
PubDate: 20170105T04:02:55Z
DOI: 10.1016/j.physd.2016.12.007
 Authors: Vladimir Jaćimović; Aladin Crnkić
 Classical quasisteady state reduction—A mathematical
characterization Authors: Alexandra Goeke; Sebastian Walcher; Eva Zerz
Abstract: Publication date: Available online 30 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Alexandra Goeke, Sebastian Walcher, Eva Zerz
We discuss parameter dependent polynomial ordinary differential equations that model chemical reaction networks. By classical quasisteady state (QSS) reduction we understand the following familiar (heuristically motivated) mathematical procedure: Set the rate of change for certain (a priori chosen) variables equal to zero and use the resulting algebraic equations to obtain a system of smaller dimension for the remaining variables. This procedure will generally be valid only for certain parameter ranges. We start by showing that the reduction is accurate if and only if the corresponding parameter is what we call a QSS parameter value, and that the reduction is approximately accurate if and only if the corresponding parameter is close to a QSS parameter value. The QSS parameter values can be characterized by polynomial equations and inequations, hence parameter ranges for which QSS reduction is valid are accessible in an algorithmic manner. A defining characteristic of a QSS parameter value is that the algebraic variety defined by the QSS relations is invariant for the differential equation. A closer investigation of the associated systems shows the existence of further invariant sets; here singular perturbations enter the picture in a natural manner. We compare QSS reduction and singular perturbation reduction, and show that, while they do not agree in general, they do, up to lowest order in a small parameter, for a quite large and relevant class of examples. This observation, in turn, allows the computation of QSS reductions even in cases where an explicit resolution of the polynomial equations is not possible.
PubDate: 20170105T04:02:55Z
DOI: 10.1016/j.physd.2016.12.002
 Authors: Alexandra Goeke; Sebastian Walcher; Eva Zerz
 Cascades of alternating pitchfork and flip bifurcations in Hbridge
inverters Authors: Viktor Avrutin; Zhanybai T. Zhusubaliyev; Erik Mosekilde
Abstract: Publication date: Available online 28 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Viktor Avrutin, Zhanybai T. Zhusubaliyev, Erik Mosekilde
Power electronic DC/AC converters (inverters) play an important role in modern power engineering. These systems are also of considerable theoretical interest because their dynamics is influenced by the presence of two vastly different forcing frequencies. As a consequence, inverter systems may be modeled in terms of piecewise smooth maps with an extremely high number of switching manifolds. We have recently shown that models of this type can demonstrate a complicated bifurcation structure associated with the occurrence of border collisions. Considering the example of a PWM Hbridge singlephase inverter, the present paper discusses a number of unusual phenomena that can occur in piecewise smooth maps with a very large number of switching manifolds. We show in particular how smooth (pitchfork and flip) bifurcations may form a macroscopic pattern that stretches across the overall bifurcation structure. We explain the observed bifurcation phenomena, show under which conditions they occur, and describe them quantitatively by means of an analytic approximation.
PubDate: 20170105T04:02:55Z
DOI: 10.1016/j.physd.2016.12.008
 Authors: Viktor Avrutin; Zhanybai T. Zhusubaliyev; Erik Mosekilde
 Spiral disk packings
 Authors: Yoshikazu Yamagishi; Takamichi Sushida
Abstract: Publication date: Available online 26 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Yoshikazu Yamagishi, Takamichi Sushida
It is shown that van Iterson’s metric for disk packings, proposed in 1907 in the study of a centric model of spiral phyllotaxis, defines a bounded distance function in the plane. This metric is also related to the bifurcation of Voronoi tilings for logarithmic spiral lattices, through the continued fraction expansion of the divergence angle. The phase diagrams of disk packings and Voronoi tilings for logarithmic spirals are dual graphs to each other. This gives a rigorous proof that van Iterson’s diagram in the centric model is connected and simply connected. It is a nonlinear analog of the duality between the phase diagrams for disk packings and Voronoi tilings on the linear lattices, having the modular group symmetry.
PubDate: 20170105T04:02:55Z
DOI: 10.1016/j.physd.2016.12.003
 Authors: Yoshikazu Yamagishi; Takamichi Sushida
 Complicated quasiperiodic oscillations and chaos from driven
piecewiseconstant circuit: Chenciner bubbles do not necessarily occur via
simple phaselocking Authors: Tri Quoc Truong; Tadashi Tsubone; Munehisa Sekikawa; Naohiko Inaba
Pages: 1  9
Abstract: Publication date: 15 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 341
Author(s): Tri Quoc Truong, Tadashi Tsubone, Munehisa Sekikawa, Naohiko Inaba
We analyze the complex quasiperiodic oscillations and chaos generated by two coupled piecewiseconstant hysteresis oscillators driven by a rectangular wave force. Oscillations generate Arnol’d resonance webs wherein lower dimensional resonance tongues extend such as that of a web in numerous directions. We provide the fundamental tools for bifurcation analysis of nonautonomous piecewiseconstant oscillators. To optimize the outstanding simplicity of piecewiseconstant circuits, we formulate a generalized procedure for calculating the Lyapunov exponents in nonautonomous piecewiseconstant dynamics. The Lyapunov exponents in these dynamics can be calculated with a precision approximately similar to that of maps. We observe twoparameter Lyapunov diagrams near the fundamental resonance region called Chenciner bubbles, wherein the oscillation frequencies of the two oscillators and the force are synchronized with a ratio of 1:1:1. Inevitably, the hysteresis considerably distorts the Chenciner bubbles. This result suggests that the Chenciner bubbles do not necessarily occur due to simple phaselocking of twodimensional tori that can be explained by homeomorphism on the circle. Furthermore, we observe the Farey sequence in the experimental measurements.
PubDate: 20161226T03:35:57Z
DOI: 10.1016/j.physd.2016.09.008
Issue No: Vol. 341 (2016)
 Authors: Tri Quoc Truong; Tadashi Tsubone; Munehisa Sekikawa; Naohiko Inaba
 Gaussian noise and the twonetwork frustrated Kuramoto model
 Authors: Andrew B. Holder; Mathew L. Zuparic; Alexander C. Kalloniatis
Pages: 10  32
Abstract: Publication date: 15 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 341
Author(s): Andrew B. Holder, Mathew L. Zuparic, Alexander C. Kalloniatis
We examine analytically and numerically a variant of the stochastic Kuramoto model for phase oscillators coupled on a general network. Two populations of phased oscillators are considered, labelled ‘Blue’ and ‘Red’, each with their respective networks, internal and external couplings, natural frequencies, and frustration parameters in the dynamical interactions of the phases. We disentangle the different ways that additive Gaussian noise may influence the dynamics by applying it separately on zero modes or normal modes corresponding to a Laplacian decomposition for the subgraphs for Blue and Red. Under the linearisation ansatz that the oscillators of each respective network remain relatively phasesynchronised centroids or clusters, we are able to obtain simple closedform expressions using the Fokker–Planck approach for the dynamics of the average angle of the two centroids. In some cases, this leads to subtle effects of metastability that we may analytically describe using the theory of ratchet potentials. These considerations are extended to a regime where one of the populations has fragmented in two. The analytic expressions we derive largely predict the dynamics of the nonlinear system seen in numerical simulation. In particular, we find that noise acting on a more tightly coupled population allows for improved synchronisation of the other population where deterministically it is fragmented.
PubDate: 20161226T03:35:57Z
DOI: 10.1016/j.physd.2016.09.009
Issue No: Vol. 341 (2016)
 Authors: Andrew B. Holder; Mathew L. Zuparic; Alexander C. Kalloniatis
 Microorganism billiards
 Authors: Saverio E. Spagnolie; Colin Wahl; Joseph Lukasik; JeanLuc Thiffeault
Pages: 33  44
Abstract: Publication date: 15 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 341
Author(s): Saverio E. Spagnolie, Colin Wahl, Joseph Lukasik, JeanLuc Thiffeault
Recent experiments and numerical simulations have shown that certain types of microorganisms “reflect” off of a flat surface at a critical angle of departure, independent of the angle of incidence. The nature of the reflection may be active (cell and flagellar contact with the surface) or passive (hydrodynamic) interactions. We explore the billiardlike motion of a body with this empirical reflection law inside a regular polygon and show that the dynamics can settle on a stable periodic orbit or can be chaotic, depending on the swimmer’s departure angle and the domain geometry. The dynamics are often found to be robust to the introduction of weak random fluctuations. The Lyapunov exponent of swimmer trajectories can be positive or negative, can have extremal values, and can have discontinuities depending on the degree of the polygon. A passive sorting device is proposed that traps swimmers of different departure angles into separate bins. We also study the external problem of a microorganism swimming in a patterned environment of square obstacles, where the departure angle dictates the possibility of trapping or diffusive trajectories.
PubDate: 20161226T03:35:57Z
DOI: 10.1016/j.physd.2016.09.010
Issue No: Vol. 341 (2016)
 Authors: Saverio E. Spagnolie; Colin Wahl; Joseph Lukasik; JeanLuc Thiffeault
 A theory of synchrony for active compartments with delays coupled through
bulk diffusion Authors: Bin Xu; Paul C. Bressloff
Pages: 45  59
Abstract: Publication date: 15 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 341
Author(s): Bin Xu, Paul C. Bressloff
We extend recent work on the analysis of synchronization in a pair of biochemical oscillators coupled by linear bulk diffusion, in order to explore the effects of discrete delays. More specifically, we consider two wellmixed, identical compartments located at either end of a bounded, onedimensional domain. The compartments can exchange signaling molecules with the bulk domain, within which the signaling molecules undergo diffusion. The concentration of signaling molecules in each compartment is modeled by a delay differential equation (DDE), while the concentration in the bulk medium is modeled by a partial differential equation (PDE) for diffusion. Coupling in the resulting PDE–DDE system is via flux terms at the boundaries. Using linear stability analysis, numerical simulations and bifurcation analysis, we investigate the effect of diffusion on the onset of a supercritical Hopf bifurcation. The direction of the Hopf bifurcation is determined by numerical simulations and a winding number argument. Near a Hopf bifurcation point, we find that there are oscillations with two possible modes: inphase and antiphase. Moreover, the critical delay for oscillations to occur increases with the diffusion coefficient. Our numerical results suggest that the selection of the inphase or antiphase oscillation is sensitive to the diffusion coefficient, time delay and coupling strength. For slow diffusion and weak coupling both modes can coexist, while for fast diffusion and strong coupling, only one of the modes is dominant, depending on the explicit choice of DDE.
PubDate: 20161226T03:35:57Z
DOI: 10.1016/j.physd.2016.10.001
Issue No: Vol. 341 (2016)
 Authors: Bin Xu; Paul C. Bressloff
 On Wright’s generalized Bessel kernel
 Authors: Lun Zhang
Pages: 92  119
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): Lun Zhang
In this paper, we consider the Wright’s generalized Bessel kernel K ( α , θ ) ( x , y ) defined by θ x α ∫ 0 1 J α + 1 θ , 1 θ ( u x ) J α + 1 , θ ( ( u y ) θ ) u α d u , α > − 1 , θ > 0 , where J a , b ( x ) = ∑ j = 0 ∞ ( − x ) j j ! Γ ( a + b j ) , a ∈ C , b > − 1 , is Wright’s generalization of the Bessel function. This nonsymmetric kernel, which generalizes the classical Bessel kernel (corresponding to θ = 1 ) in random matrix theory, is the hard edge scaling limit of the correlation kernel for certain Muttalib–Borodin ensembles. We show that, if θ is rational, i.e., θ = m n with m , n ∈ N , g c d ( m , n ) = 1 , and α > m − 1 − m n , the Wright’s generalized Bessel kernel is integrable in the sense of Its–Izergin–Korepin–Slavnov. We then come to the Fredholm determinant of this kernel over the union of several scaled intervals, which can also be interpreted as the gap probability (the probability of finding no particles) on these intervals. The integrable structure allows us to obtain a system of coupled partial differential equations associated with the corresponding Fredholm determinant as well as a Hamiltonian interpretation. As a consequence, we are able to represent the gap probability over a single interval ( 0 , s ) in terms of a solution of a system of nonlinear ordinary differential equations.
PubDate: 20161218T03:09:08Z
DOI: 10.1016/j.jat.2016.09.002
Issue No: Vol. 213 (2016)
 Authors: Lun Zhang
 Windows of opportunity for synchronization in stochastically coupled maps
 Authors: Olga Golovneva; Russell Jeter Igor Belykh Maurizio Porfiri
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): Olga Golovneva, Russell Jeter, Igor Belykh, Maurizio Porfiri
Several complex systems across science and engineering display on–off intermittent coupling among their units. Most of the current understanding of synchronization in switching networks relies on the fast switching hypothesis, where the network dynamics evolves at a much faster time scale than the individual units. Recent numerical evidence has demonstrated the existence of windows of opportunity, where synchronization may be induced through nonfast switching. Here, we study synchronization of coupled maps whose coupling gains stochastically switch with an arbitrary switching period. We determine the role of the switching period on synchronization through a detailed analytical treatment of the Lyapunov exponent of the stochastic dynamics. Through closedform expressions and numerical findings, we demonstrate the emergence of windows of opportunity and elucidate their nontrivial relationship with the stability of synchronization under static coupling. Our results are expected to provide a rigorous basis for understanding the dynamic mechanisms underlying the emergence of windows of opportunity and leverage nonfast switching in the design of evolving networks.
Graphical abstract
PubDate: 20161218T03:09:08Z
 Authors: Olga Golovneva; Russell Jeter Igor Belykh Maurizio Porfiri
 Arnold tongues in a billiard problem in nonlinear and nonequilibrium
systems Authors: Tomoyuki Miyaji
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): Tomoyuki Miyaji
We study a billiard problem in nonlinear and nonequilibrium systems. This is motivated by the motions of a traveling spot in a reaction–diffusion system (RDS) in a rectangular domain. We consider a fourdimensional dynamical system, defined by ordinary differential equations. This was first derived by S.I. Ei et al. (2006), based on a reduced system on the center manifold in a neighborhood of a pitchfork bifurcation of a stationary spot for the RDS. In contrast to the classical billiard problem, this defines a dynamical system that is dissipative rather than conservative, and has an attractor. According to previous numerical studies, the attractor of the system changes depending on parameters such as the aspect ratio of the domain. It may be periodic, quasiperiodic, or chaotic. In this paper, we elucidate that it results from parameters crossing Arnold tongues and that the organizing center is a Hopf–Hopf bifurcation of the trivial equilibrium.
PubDate: 20161218T03:09:08Z
 Authors: Tomoyuki Miyaji
 Corrigendum to “Dynamics and stability of a discrete breather in a
harmonically excited chain with vibroimpact onsite potential” [Physica
D 292–293 (2015) 8–28] Authors: Nathan Perchikov; O.V. Gendelman
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): Nathan Perchikov, O.V. Gendelman
PubDate: 20161218T03:09:08Z
 Authors: Nathan Perchikov; O.V. Gendelman
 Timing variation in an analytically solvable chaotic system
 Authors: J.N. Blakely; M.S. Milosavljevic N.J. Corron
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): J.N. Blakely, M.S. Milosavljevic, N.J. Corron
We present analytic solutions for a chaotic dynamical system that do not have the regular timing characteristic of recently reported solvable chaotic systems. The dynamical system can be viewed as a first order filter with binary feedback. The feedback state may be switched only at instants defined by an external clock signal. Generalizing from a period one clock, we show analytic solutions for period two and higher period clocks. We show that even when the clock ‘ticks’ randomly the chaotic system has an analytic solution. These solutions can be visualized in a stroboscopic map whose complexity increases with the complexity of the clock. We provide both analytic results as well as experimental data from an electronic circuit implementation of the system. Our findings bridge the gap between the irregular timing of well known chaotic systems such as Lorenz and Rossler and the well regulated oscillations of recently reported solvable chaotic systems.
Graphical abstract
PubDate: 20161218T03:09:08Z
 Authors: J.N. Blakely; M.S. Milosavljevic N.J. Corron
 Databased stochastic model reduction for the Kuramoto–Sivashinsky
equation Authors: Fei Kevin; Lin Alexandre Chorin
Abstract: Publication date: 1 February 2017
Source:Physica D: Nonlinear Phenomena, Volume 340
Author(s): Fei Lu, Kevin K. Lin, Alexandre J. Chorin
The problem of constructing databased, predictive, reduced models for the Kuramoto–Sivashinsky equation is considered, under circumstances where one has observation data only for a small subset of the dynamical variables. Accurate prediction is achieved by developing a discretetime stochastic reduced system, based on a NARMAX (Nonlinear Autoregressive Moving Average with eXogenous input) representation. The practical issue, with the NARMAX representation as with any other, is to identify an efficient structure, i.e., one with a small number of terms and coefficients. This is accomplished here by estimating coefficients for an approximate inertial form. The broader significance of the results is discussed.
PubDate: 20161218T03:09:08Z
 Authors: Fei Kevin; Lin Alexandre Chorin
 Numerical bifurcation for the capillary Whitham equation
 Authors: Filippo Remonato; Henrik Kalisch
Abstract: Publication date: Available online 6 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Filippo Remonato, Henrik Kalisch
The socalled Whitham equation arises in the modeling of free surface water waves, and combines a generic nonlinear quadratic term with the exact linear dispersion relation for gravity waves on the free surface of a fluid with finite depth. In this work, the effect of incorporating capillarity into the Whitham equation is in focus. The capillary Whitham equation is a nonlocal equation similar to the usual Whitham equation, but containing an additional term with a coefficient depending on the Bond number which measures the relative strength of capillary and gravity effects on the wave motion. A spectral collocation scheme for computing approximations to periodic traveling waves for the capillary Whitham equation is put forward. Numerical approximations of periodic traveling waves are computed using a bifurcation approach, and a number of bifurcation curves are found. Our analysis uncovers a rich structure of bifurcation patterns, including subharmonic bifurcations, as well as connecting and crossing branches. Indeed, for some values of the Bond number, the bifurcation diagram features distinct branches of solutions which intersect at a secondary bifurcation point. The same branches may also cross without connecting, and some bifurcation curves feature selfcrossings without selfconnections.
PubDate: 20161211T02:34:16Z
DOI: 10.1016/j.physd.2016.11.003
 Authors: Filippo Remonato; Henrik Kalisch
 One and twodimensional bright solitons in inhomogeneous defocusing
nonlinearities with an antisymmetric periodic gain and loss Authors: Dengchu Guo; Jing Xiao; Linlin Gu; Hongzhen Jin; Liangwei Dong
Abstract: Publication date: Available online 2 December 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Dengchu Guo, Jing Xiao, Linlin Gu, Hongzhen Jin, Liangwei Dong
We address that various branches of bright solitons exist in a spatially inhomogeneous defocusing nonlinearity with an imprinted antisymmetric periodic gainloss profile. The spectra of such systems with a purely imaginary potential never become complex and thus the paritytime symmetry is unbreakable. The mergence between pairs of soliton branches occurs at a critical gainloss strength, above which no soliton solutions can be found. Intriguingly, which pair of soliton branches will merge together can be changed by varying the modulation frequency of gain and loss. Most branches of onedimensional solitons are stable in wide parameter regions. We also provide the first example of twodimensional bright solitons with unbreakable paritytime symmetry.
PubDate: 20161204T02:03:09Z
DOI: 10.1016/j.physd.2016.11.005
 Authors: Dengchu Guo; Jing Xiao; Linlin Gu; Hongzhen Jin; Liangwei Dong
 Excitability, mixedmode oscillations and transition to chaos in a
stochastic ice ages model Authors: D.V. Alexandrov; I.A. Bashkirtseva; L.B. Ryashko
Abstract: Publication date: Available online 30 November 2016
Source:Physica D: Nonlinear Phenomena
Author(s): D.V. Alexandrov, I.A. Bashkirtseva, L.B. Ryashko
Motivated by an important geophysical significance, we consider the influence of stochastic forcing on a simple threedimensional climate model previously derived by Saltzman and Sutera. A nonlinear dynamical system governing three physical variables, the bulk ocean temperature, continental and marine ice masses, is analyzed in deterministic and stochastic cases. It is shown that the attractor of deterministic model is either a stable equilibrium or a limit cycle. We demonstrate that the process of continental ice melting occurs with a noisedependent time delay as compared with marine ice melting. The paleoclimate cyclicity which is near 100 ky in a wide range of model parameters abruptly increases in the vicinity of a bifurcation point and depends on the noise intensity. In a zone of stable equilibria, the 3D climate model under consideration is extremely excitable. Even for a weak random noise, the stochastic trajectories demonstrate a transition from small to largeamplitude stochastic oscillations (SLASO). In a zone of stable cycles, SLASO transitions are analyzed too. We show that such stochastic transitions play an important role in the formation of a mixedmode paleoclimate scenario. This mixedmode dynamics with the intermittency of large and smallamplitude stochastic oscillations and coherence resonance are investigated via analysis of interspike intervals. A tendency of dynamic paleoclimate to abrupt and rapid glaciations and deglaciations as well as its transition from order to chaos with increasing noise are shown.
PubDate: 20161204T02:03:09Z
DOI: 10.1016/j.physd.2016.11.007
 Authors: D.V. Alexandrov; I.A. Bashkirtseva; L.B. Ryashko
 Optical dispersive shock waves in defocusing colloidal media
 Authors: X. An; T.R. Marchant; N.F. Smyth
Abstract: Publication date: Available online 24 November 2016
Source:Physica D: Nonlinear Phenomena
Author(s): X. An, T.R. Marchant, N.F. Smyth
The propagation of an optical dispersive shock wave, generated from a jump discontinuity in light intensity, in a defocussing colloidal medium is analysed. The equations governing nonlinear light propagation in a colloidal medium consist of a nonlinear Schrödinger equation for the beam and an algebraic equation for the medium response. In the limit of low light intensity, these equations reduce to a perturbed higher order nonlinear Schrödinger equation. Solutions for the leading and trailing edges of the colloidal dispersive shock wave are found using modulation theory. This is done for both the perturbed nonlinear Schrödinger equation and the full colloid equations for arbitrary light intensity. These results are compared with numerical solutions of the colloid equations.
PubDate: 20161127T16:42:30Z
DOI: 10.1016/j.physd.2016.11.004
 Authors: X. An; T.R. Marchant; N.F. Smyth
 Spatiotemporal control to eliminate cardiac alternans using isostable
reduction Authors: Dan Wilson; Jeff Moehlis
Abstract: Publication date: Available online 16 November 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Dan Wilson, Jeff Moehlis
Cardiac alternans, an arrhythmia characterized by a beattobeat alternation of cardiac action potential durations, is widely believed to facilitate the transition from normal cardiac function to ventricular fibrillation and sudden cardiac death. Alternans arises due to an instability of a healthy period1 rhythm, and most dynamical control strategies either require extensive knowledge of the cardiac system, making experimental validation difficult, or are model independent and sacrifice important information about the specific system under study. Isostable reduction provides an alternative approach, in which the response of a system to external perturbations can be used to reduce the complexity of a cardiac system, making it easier to work with from an analytical perspective while retaining many of its important features. Here, we use isostable reduction strategies to reduce the complexity of partial differential equation models of cardiac systems in order to develop energy optimal strategies for the elimination of alternans. Resulting control strategies require significantly less energy to terminate alternans than comparable strategies and do not require continuous state feedback.
PubDate: 20161120T14:57:33Z
DOI: 10.1016/j.physd.2016.11.001
 Authors: Dan Wilson; Jeff Moehlis
 Macroscopic heat transport equations and heat waves in nonequilibrium
states Authors: Yangyu Guo; David Jou; Moran Wang
Abstract: Publication date: Available online 16 November 2016
Source:Physica D: Nonlinear Phenomena
Author(s): Yangyu Guo, David Jou, Moran Wang
Heat transport may behave as wave propagation when the time scale of processes decreases to be comparable to or smaller than the relaxation time of heat carriers. In this work, a generalized heat transport equation including nonlinear, nonlocal and relaxation terms is proposed, which sums up the CattaneoVernotte, dualphaselag and phonon hydrodynamic models as special cases. In the frame of this equation, the heat wave propagations are investigated systematically in nonequilibrium steady states, which were usually studied around equilibrium states. The phase (or front) speed of heat waves is obtained through a perturbation solution to the heat differential equation, and found to be intimately related to the nonlinear and nonlocal terms. Thus, potential heat wave experiments in nonequilibrium states are devised to measure the coefficients in the generalized equation, which may throw light on understanding the physical mechanisms and macroscopic modeling of nanoscale heat transport.
PubDate: 20161120T14:57:33Z
DOI: 10.1016/j.physd.2016.10.005
 Authors: Yangyu Guo; David Jou; Moran Wang