Unidirectionally coupled map lattice as a model for open flow systems

Topological early warning signals: Quantifying varying routes to extinction in a spatially distributed population model

Understanding and predicting critical transitions in spatially explicit ecological systems is particularly challenging due to their complex spatial and temporal dynamics and high dimensionality. Here, we explore changes in population distribution patterns during a critical transition (an extinction event) using computational topology. Computational topology allows us to quantify certain features of a population distribution pattern, such as the level of fragmentation. We create population distribution patterns via a simple coupled patch model with Ricker map growth and nearest neighbors dispersal on a two dimensional lattice. We observe two dominant paths to extinction within the explored parameter space that depend critically on the dispersal rate d and the rate of parameter drift, Δϵ. These paths to extinction are easily topologically distinguishable, so categorization can be automated. We use this population model as a theoretical proof-of-concept for the methodology, and argue that computational topology is a powerful tool for analyzing dynamical changes in systems with noisy data that are coarsely resolved in space and/or time. In addition, computational topology can provide early warning signals for chaotic dynamical systems where traditional statistical early warning signals would fail. For these reasons, we envision this work as a helpful addition to the critical transitions prediction toolbox.

Characteristic distribution of finite-time Lyapunov exponents for chimera states

Our fascination with chimera states stems partially from the somewhat paradoxical, yet fundamental trait of identical, and identically coupled, oscillators to split into spatially separated, coherently and incoherently oscillating groups. While the list of systems for which various types of chimeras have already been detected continues to grow, there is a corresponding increase in the number of mathematical analyses aimed at elucidating the fundamental reasons for this surprising behaviour. Based on the model systems, there are strong indications that chimera states may generally be ubiquitous in naturally occurring systems containing large numbers of coupled oscillators - certain biological systems and high-Tc superconducting materials, for example. In this work we suggest a new way of detecting and characterising chimera states. Specifically, it is shown that the probability densities of finite-time Lyapunov exponents, corresponding to chimera states, have a definite characteristic shape. Such distributions could be used as signatures of chimera states, particularly in systems for which the phases of all the oscillators cannot be measured directly. For such cases, we suggest that chimera states could perhaps be detected by reconstructing the characteristic distribution via standard embedding techniques, thus making it possible to detect chimera states in systems where they could otherwise exist unnoticed.

Hierarchy and polysynchrony in an adaptive network

We describe a simple adaptive network of coupled chaotic maps. The network reaches a stationary state (frozen topology) for all values of the coupling parameter, although the dynamics of the maps at the nodes of the network can be nontrivial. The structure of the network shows interesting hierarchical properties and in certain parameter regions the dynamics is polysynchronous: Nodes can be divided in differently synchronized classes but, contrary to cluster synchronization, nodes in the same class need not be connected to each other. These complicated synchrony patterns have been conjectured to play roles in systems biology and circuits. The adaptive system we study describes ways whereby this behavior can evolve from undifferentiated nodes.

Pulses of chaos synchronization in coupled map chains with delayed transmission

Pulses of synchronization in chaotic coupled map lattices are discussed in the context of transmission of information. Synchronization and desynchronization propagate along the chain with different velocities which are calculated analytically from the spectrum of convective Lyapunov exponents. Since the front of synchronization travels slower than the front of desynchronization, the maximal possible chain length for which information can be transmitted by modulating the first unit of the chain is bounded.

Spatial chaos of traveling waves has a given velocity

We study the complexity of stable waves in unidirectional bistable coupled map lattices as a test tube to spatial chaos of traveling patterns in open flows. Numerical calculations reveal that, grouping patterns into sets according to their velocity, at most one set of waves has positive topological entropy for fixed parameters. By using symbolic dynamics and shadowing, we analytically determine velocity-dependent parameter domains of existence of pattern families with positive entropy. These arguments provide a method to exhibit chaotic sets of stable waves with arbitrary velocity in extended systems.

Destabilization patterns in chains of coupled oscillators

We describe the mechanism of destabilization in a chain of identical coupled oscillators. Along with the transition from stationary to oscillatory behavior of the single oscillator, the network undergoes a complicated bifurcation scenario including the coexistence of multiple periodic orbits with different frequencies, spatial patterns, and modulation instabilities. This scenario, which is similar to the well-known Eckhaus scenario in spatially extended systems, occurs here also in the case of purely convective unidirectional coupling, and hence it cannot be explained as a simple discretization of its spatially continuous counterpart. Although the number of coexisting periodic orbits grows with the number of oscillators, we are able to treat this problem independently of the actual size of the network by investigating the limiting equations for the related spectral problems.

Chaotic spatial bifurcation by complex coupling

A spatial bifurcation (a transition from stationary to oscillatory regime) in a chain of unidirectionally coupled phase systems is studied. It is shown that complication of coupling terms can make this bifurcation spatially chaotic in contrast to the previously observed "regular" and "predictable" type. It is demonstrated that the found type of spatial bifurcation corresponds to a smooth (predictable) manifold in the parameter space, while its spatial location gets actually unpredictable being governed by regularities of chaotic behavior. We infer that complex collective dynamics may arise in networks with plain architecture and simple dynamics of individual elements if nontrivial coupling is realized.

Hydrodynamic Lyapunov modes in coupled map lattices

In this paper, numerical and analytical results are presented which indicate that hydrodynamic Lyapunov modes (HLMs) also exist for coupled map lattices (CMLs). The dispersion relations for the HLMs of CMLs are found to fall into two different universality classes. It is characterized by lambda approximately k for coupled standard maps and lambda approximately k2 for coupled circle maps. The conditions under which HLMs can be observed are discussed. The role of the Hamiltonian structure, conservation laws, translational invariance, and damping is elaborated. Our results are as follows: (1) The Hamiltonian structure is not a necessary condition for the existence of HLMs. (2) Conservation laws or the translational invariance alone cannot guarantee the existence of HLMs. (3) Including a damping term in the system of coupled Hamiltonian maps does not destroy the HLMs. The lambda-k dispersion relation of HLMs, however, changes to the universality class with lambda-k2 under damping. In contrast, no HLMs survives in the system of coupled circle maps under damping. (4) An on-site potential destroys the HLMs. (5) The study of zero-value Lyapunov exponents (LEs) and associated Lyapunov vectors (LVs) shows that translational invariance and conservation laws play different roles in the tangent space dynamics. (6) The dynamics of the coordinate and momentum parts of LVs in Hamiltonian systems are related but different. Furthermore, numerical results for a two-dimensional system show that the appearance of HLMs in CMLs is not restricted to the one-dimensional case.

Dynamical behavior of hydrodynamic Lyapunov modes in coupled map lattices

In our previous study of hydrodynamic Lyapunov modes (HLMs) in coupled map lattices, we found that there are two classes of systems with different lambda-k dispersion relations. For coupled circle maps we found the quadratic dispersion relations lambda approximately k2 and lambda approximately k for coupled standard maps. Here, we carry out further numerical experiments to investigate the dynamic Lyapunov vector (LV) structure factor which can provide additional information on the Lyapunov vector dynamics. The dynamic LV structure factor of coupled circle maps is found to have a single peak at omega=0 and can be well approximated by a single Lorentzian curve. This implies that the hydrodynamic Lyapunov modes in coupled circle maps are nonpropagating and show only diffusive motion. In contrast, the dynamic LV structure factor of coupled standard maps possesses two visible sharp peaks located symmetrically at +/- omega. The spectrum can be well approximated by the superposition of three Lorentzian curves centered at omega=0 and +/-omegau, respectively. In addition, the omega-k dispersion relation takes the form omegau=cuk for k --> 2pi/L. These facts suggest that the hydrodynamic Lyapunov modes in coupled standard maps are propagating. The HLMs in the two classes of systems are shown to have different dynamical behavior besides their difference in spatial structure. Moreover, our simulations demonstrate that adding damping to coupled standard maps turns the propagating modes into diffusive ones alongside a change of the lambda-k dispersion relation from lambda approximately k to lambda approximately k2. In cases of weak damping, there is a crossover in the dynamic LV structure factors; i.e., the spectra with smaller k are akin to those of coupled circle maps while the spectra with larger k are similar to those of coupled standard maps.

A self-learning coupled map lattice for vortex shedding in cable and cylinder wakes

A coupled map lattice (CML) with self-learning features is developed to model flow over freely vibrating cables and stationary cylinders at low Reynolds numbers. Coupled map lattices that combine a series of low-dimensional circle maps with a diffusion model have been used previously to predict qualitative features of these flows. However, the simple nature of these CML models implies that there will be unmodeled wake features if a detailed, quantitative comparison is made with laboratory or simulated wake flows. Motivated by a desire to develop an improved CML model, we incorporate self-learning features into a new CML that is first trained to precisely estimate wake patterns from a target numerical simulation. A new convective-diffusive map that includes additional wake dynamics is developed. The new self-learning CML uses an adaptive estimation scheme (multivariable least-squares algorithm). Studies of this approach are conducted using wake patterns from a Navier-Stokes solution (spectral element-based NEKTAR simulation) of freely vibrating cable wakes at Reynolds numbers Re=100. It is shown that the self-learning model accurately and efficiently estimates the simulated wake patterns. The self-learning scheme is then successfully applied to vortex shedding patterns obtained from experiments on stationary cylinders. This constitutes a first step toward the use of the self-learning CML as a wake model in flow control studies of laboratory wake flows.

Amplification of noise in a cascade chemical reaction

Networks of chemical reactions have been given much attention recently. However, dynamical aspects of such networks remain to be elucidated. In this paper, we study a cascade chemical reaction, consisting of a series of downstream-coupled Brusselators. Along the cascade of reaction, small fluctuations naturally existing in the concentration of chemical species are amplified. Such amplification of small noise leads to the formation of chemical oscillations in the downstream chemical species. The amplification rate of small noise in the concentration along the cascade is studied and the method to calculate the amplification rate analytically is developed. It is also shown that the nonlinear evolution of the chemical oscillation in the downstream reaction strongly depends on the frequency of the initial inlet chemical concentration.

Spiral waves in a coupled network of sine-circle maps

A coupled two-dimensional lattice of sine-circle maps is investigated numerically as a simple model for coupled network of nonlinear oscillators under a spatially uniform, temporally periodic, external forcing. Various patterns, including quasiperiodic spiral waves, periodic, banded spiral waves in several different polygonal shapes, and domain patterns, are observed. The banded spiral waves and domain patterns match well with the results of earlier experimental studies. Several transitions are analyzed. Among others, the source-sink transition of a quasiperiodic spiral wave and the cascade of "side-doubling" bifurcations of polygonal spiral waves are of great interest.

Bicritical scaling behavior in unidirectionally coupled oscillators

We study the scaling behavior of period doublings in a system of two unidirectionally coupled parametrically forced pendulums near a bicritical point where two critical lines of period-doubling transition to chaos in both subsystems meet. When crossing a bicritical point, a hyperchaotic attractor with two positive Lyapunov exponents appears, i.e., a transition to hyperchaos occurs. Varying the control parameters of the two subsystems, the unidirectionally coupled parametrically forced pendulums exhibit multiple period-doubling transitions to hyperchaos. For each transition to hyperchaos, using both a "residue-matching" renormalization group method and a direct numerical method, we make an analysis of the bicritical scaling behavior. It is thus found that the second response subsystem exhibits a new type of non-Feigenbaum scaling behavior, while the first drive subsystem is in the usual Feigenbaum critical state. The universality of the bicriticality is also examined for several different types of unidirectional couplings.

Bicritical behavior of period doublings in unidirectionally coupled maps

We study the scaling behavior of period doublings in two unidirectionally coupled one-dimensional maps near a bicritical point where two critical lines of period-doubling transition to chaos in both subsystems meet. Note that the bicritical point corresponds to a border of chaos in both subsystems. For this bicritical case, the second response subsystem exhibits a type of non-Feigenbaum critical behavior, while the first drive subsystem is in the Feigenbaum critical state. Using two different methods, we make the renormalization-group analysis of the bicritical behavior and find the corresponding fixed point of the renormalization transformation with two relevant eigenvalues. The scaling factors obtained by the renormalization-group analysis agree well with those obtained by a direct numerical method.