Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control chemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

Aerodynamic noise from wave-turbulence interaction

Abstract. In order to estimate the acoustic energy scattered when a unit volume of free turbulence, such as in free jets, interacts with a plane steady sound wave, theoretical expressions are derived for two simple models of turbulence: eddy model and isotropic model. The effect of convection by mean motion of the energy-bearing eddies on the incident sound wave and on the sound generated from wave-turbulence interaction is taken into account. Finally, by means of a representative calculation,the directionality pattern and Mach number dependence of the noise so generated is discussed.

Low-dimensional chaos and wave turbulence in plasmas

We investigated drift-wave turbulence in the plasma edge of a small tokamak by considering solutions of the Hasegawa-Mima equation involving three interacting modes in Fourier space. The resulting low-dimensional dynamics presented periodic as well as chaotic evolution of the Fourier-mode amplitudes, and we performed the control of chaotic behaviour through the application of a fourth resonant wave of small amplitude.

SYNCHRONIZATION OF CHAOS AND THE TRANSITION TO WAVE TURBULENCE

We investigated the transition to wave turbulence in a spatially extended three-wave interacting model, where a spatially homogeneous state undergoing chaotic dynamics undergoes spatial mode excitation. The transition to this weakly turbulent state can be regarded as the loss of synchronization of chaos of mode oscillators describing the spatial dynamics.

Spectral bifurcations in dispersive wave turbulence

Dispersive wave turbulence is studied numerically for a class of one-dimensional nonlinear wave equations. Both deterministic and random (white noise in time) forcings are studied. Four distinct stable spectra are observed—the direct and inverse cascades of weak turbulence (WT) theory, thermal equilibrium, and a fourth spectrum (MMT; Majda, McLaughlin, Tabak). Each spectrum can describe long-time behavior, and each can be only metastable (with quite diverse lifetimes)—depending on details of nonlinearity, forcing, and dissipation. Cases of a long-live MMT transient state dcaying to a state with WT spectra, and vice-versa, are displayed. In the case of freely decaying turbulence, without forcing, both cascades of weak turbulence are observed. These WT states constitute the clearest and most striking numerical observations of WT spectra to date—over four decades of energy, and three decades of spatial, scales. Numerical experiments that study details of the composition, coexistence, and transition between spectra are then discussed, including: (i) for deterministic forcing, sharp distinctions between focusing and defocusing nonlinearities, including the role of long wavelength instabilities, localized coherent structures, and chaotic behavior; (ii) the role of energy growth in time to monitor the selection of MMT or WT spectra; (iii) a second manifestation of the MMT spectrum as it describes a self-similar evolution of the wave, without temporal averaging; (iv) coherent structures and the evolution of the direct and inverse cascades; and (v) nonlocality (in k-space) in the transferral process.

Phase transitions in wave turbulence

We consider turbulence within the Gross-Pitaevsky model and look into the creation of a coherent condensate via an inverse cascade originating at small scales. The growth of the condensate leads to a spontaneous breakdown of statistical symmetries of overcondensate fluctuations: First, isotropy is broken, then a series of phase transitions marks the changing symmetry from twofold to threefold to fourfold. We describe respective anisotropic flux flows in the k space. At the highest level reached, we observe a short-range positional and long-range orientational order (as in a hexatic phase). In other words, the more one pumps the system, the more ordered the system becomes. The phase transitions happen when the system is pumped by an instability term and does not occur when pumped by a random force. We thus demonstrate nonuniversality of an inverse-cascade turbulence with respect to the nature of small-scale forcing.

Optical wave turbulence

We review recent progress in optical wave turbulence with a specific focus on the fast growing field of fibre lasers. Weak irregular nonlinear interactions between a large number of resonator modes are responsible for practically important characteristics of fibre lasers such as spectral broadening of radiation. Wave turbulence is a fundamental nonlinear phenomenon which occurs in a variety of nonlinear wave-bearing physical systems. The experimental impediments and the computationally intensive nature of simulating of hydrodynamic or plasma wave turbulence often make it rather challenging to collect a significant number of statistical data The study of turbulent wave behaviour in optical devices offers quite a unique opportunity to collect an enormous amount of data on statistical properties of wave turbulence using high-speed, high precision optical measurements during a relatively short period of time. We present recent theoretical, numerical and experimental results on optical wave turbulence in fibre lasers ranging from weak to strong developed turbulence for different signs of fibre dispersion. Furthermore, we report on our studies of spectral wave condensate in fibre lasers that make interdisciplinary links with a number of other research fields.

Optical wave turbulence

We review recent progress in optical wave turbulence with a specific focus on the fast growing field of fibre lasers. Weak irregular nonlinear interactions between a large number of resonator modes are responsible for practically important characteristics of fibre lasers such as spectral broadening of radiation. Wave turbulence is a fundamental nonlinear phenomenon which occurs in a variety of nonlinear wave-bearing physical systems. The experimental impediments and the computationally intensive nature of simulating of hydrodynamic or plasma wave turbulence often make it rather challenging to collect a significant number of statistical data The study of turbulent wave behaviour in optical devices offers quite a unique opportunity to collect an enormous amount of data on statistical properties of wave turbulence using high-speed, high precision optical measurements during a relatively short period of time. We present recent theoretical, numerical and experimental results on optical wave turbulence in fibre lasers ranging from weak to strong developed turbulence for different signs of fibre dispersion. Furthermore, we report on our studies of spectral wave condensate in fibre lasers that make interdisciplinary links with a number of other research fields.

Exact solution for the energy spectrum of Kelvin-wave turbulence in superfluids

We study the statistical and dynamical behavior of turbulent Kelvin waves propagating on quantized vortices in superfluids and address the controversy concerning the energy spectrum that is associated with these excitations. Finding the correct energy spectrum is important because Kelvin waves play a major role in the dissipation of energy in superfluid turbulence at near-zero temperatures. In this paper, we show analytically that the solution proposed by [L’vov and Nazarenko, JETP Lett. 91, 428 (2010)] enjoys existence, uniqueness, and regularity of the prefactor. Furthermore, we present numerical results of the dynamical equation that describes to leading order the nonlocal regime of the Kelvin-wave dynamics. We compare our findings with the analytical results from the proposed local and nonlocal theories for Kelvin-wave dynamics and show an agreement with the nonlocal predictions. Accordingly, the spectrum proposed by L’vov and Nazarenko should be used in future theories of quantum turbulence. Finally, for weaker wave forcing we observe an intermittent behavior of the wave spectrum with a fluctuating dissipative scale, which we interpreted as a finite-size effect characteristic of mesoscopic wave turbulence.

One-dimensional optical wave turbulence:experiment and theory

We present a review of the latest developments in one-dimensional (1D) optical wave turbulence (OWT). Based on an original experimental setup that allows for the implementation of 1D OWT, we are able to show that an inverse cascade occurs through the spontaneous evolution of the nonlinear field up to the point when modulational instability leads to soliton formation. After solitons are formed, further interaction of the solitons among themselves and with incoherent waves leads to a final condensate state dominated by a single strong soliton. Motivated by the observations, we develop a theoretical description, showing that the inverse cascade develops through six-wave interaction, and that this is the basic mechanism of nonlinear wave coupling for 1D OWT. We describe theory, numerics and experimental observations while trying to incorporate all the different aspects into a consistent context. The experimental system is described by two coupled nonlinear equations, which we explore within two wave limits allowing for the expression of the evolution of the complex amplitude in a single dynamical equation. The long-wave limit corresponds to waves with wave numbers smaller than the electrical coherence length of the liquid crystal, and the opposite limit, when wave numbers are larger. We show that both of these systems are of a dual cascade type, analogous to two-dimensional (2D) turbulence, which can be described by wave turbulence (WT) theory, and conclude that the cascades are induced by a six-wave resonant interaction process. WT theory predicts several stationary solutions (non-equilibrium and thermodynamic) to both the long- and short-wave systems, and we investigate the necessary conditions required for their realization. Interestingly, the long-wave system is close to the integrable 1D nonlinear Schrödinger equation (NLSE) (which contains exact nonlinear soliton solutions), and as a result during the inverse cascade, nonlinearity of the system at low wave numbers becomes strong. Subsequently, due to the focusing nature of the nonlinearity, this leads to modulational instability (MI) of the condensate and the formation of solitons. Finally, with the aid of the probability density function (PDF) description of WT theory, we explain the coexistence and mutual interactions between solitons and the weakly nonlinear random wave background in the form of a wave turbulence life cycle (WTLC).

Optical wave turbulence and wave condensation in a nonlinear optical experiment

We present theory, numerical simulations and experimental observations of a 1D optical wave system. We show that this system is of a dual cascade type, namely, the energy cascading directly to small scales, and the photons or wave action cascading to large scales. In the optical context the inverse cascade is particularly interesting because it means the condensation of photons. We show that the cascades are induced by a six-wave resonant interaction process described by weak turbulence theory. We show that by starting with weakly nonlinear randomized waves as an initial condition, there exists an inverse cascade of photons towards the lowest wavenumbers. During the cascade nonlinearity becomes strong at low wavenumbers and, due to the focusing nature of the nonlinearity, it leads to modulational instability resulting in the formation of solitons. Further interaction of the solitons among themselves and with incoherent waves leads to the final condensate state dominated by a single strong soliton. In addition, we show the existence of the direct energy cascade numerically and that it agrees with the wave turbulence prediction.

Wave and turbulence measurements at a tidal energy site

This work presents the analysis of wave and turbulence measurements collected at a tidal energy site. A new method is introduced to produce more consistent and rigorous estimations of the velocity fluctuations power spectral densities. An analytical function is further proposed to fit the observed spectra and could be input to the numerical models predicting power production and structural loading on tidal turbines. Another new approach is developed to correct for the effect of the Doppler noise on the high frequencies power spectral densities. The analysis of velocity time series combining wave and turbulent contributions demonstrates that the turbulent motions are coherent throughout the water column, rendering the wave coherence-based methods not applicable to our dataset. To avoid this problem, an alternative approach relying on the pressure data collected by the ADCP is introduced and shows appreciable improvement in the wave-turbulence separation.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Coronal loops heated by turbulence-driven Alfven waves

Li, Xing, Habbal, S. R., 'Coronal loops heated by turbulence-driven Alfven waves', The Astrophysical Journal, (2003) 598(2) pp.L125-L128 RAE2008

Evolution of slow electrostatic shock into a plasma shock mediated by electrostatic turbulence

The collision of two plasma clouds at a speed that exceeds the ion acoustic speed can result in the formation of shocks. This phenomenon is observed not only in astrophysical scenarios, such as the propagation of supernova remnant (SNR) blast shells into the interstellar medium, but also in laboratory-based laser-plasma experiments. These experiments and supporting simulations are thus seen as an attractive platform for small-scale reproduction and study of astrophysical shocks in the laboratory. We model two plasma clouds, which consist of electrons and ions, with a 2D particle-in-cell simulation. The ion temperatures of both clouds differ by a factor of ten. Both clouds collide at a speed that is realistic for laboratory studies and for SNR shocks in their late evolution phase, like that of RCW86. A magnetic field, which is orthogonal to the simulation plane, has a strength that is comparable to that of SNR shocks. A forward shock forms between the overlap layer of both plasma clouds and the cloud with cooler ions. A large-amplitude ion acoustic wave is observed between the overlap layer and the cloud with hotter ions. It does not steepen into a reverse shock because its speed is below the ion acoustic speed. A gradient of the magnetic field amplitude builds up close to the forward shock as it compresses the magnetic field. This gradient gives rise to an electron drift that is fast enough to trigger an instability. Electrostatic ion acoustic wave turbulence develops ahead of the shock, widens its transition layer, and thermalizes the ions, but the forward shock remains intact. © 2014 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.