Localised edge states nucleate turbulence in extended plane Couette cells

We study the turbulence transition of plane Couette flow in large domains where localised perturbations are observed to generate growing turbulent spots. Extending previous studies on the boundary between laminar and turbulent dynamics we determine invariant structures intermediate between laminar and turbulent flow. In wide but short domains we find states that are localised in spanwise direction, and in wide and long domains the states are also localised in downstream direction. These localised states act as critical nuclei for the transition to turbulence in spatially extended domains.Comment: 15 pages, 5 figure

Role of Fiber Orientation in Atrial Arrythmogenesis

Electrical wave-front propagation in the atria is determined largely by local fiber orientation. Recent study suggests that atrial fibrillation (AF) progresses with enhanced anisotropy. In this work, a 3D rabbit atrial anatomical model at 20 × 20 × 20 μm3 resolution with realistic fiber orientation was constructed based on the novel contrast-enhanced micro-CT imaging. The Fenton-Karma cellular activation model was adapted to reproduce rabbit atrial action potential period of 80 ms. Diffusivities were estimated for longitudinal and transverse directions of the fiber orientation respectively. Pacing was conducted in the 3D anisotropic atrial model with a reducing S2 interval to facilitate initiation of atrial arrhythmia. Multiple simulations were conducted with varying values of diffusion anisotropy and stimulus locations to evaluate the role of anisotropy in initiating AF. Under physiological anisotropy conditions, a rapid right atrial activation was followed by the left atrial activation. Excitation waves reached the atrio-ventricular border where they terminated. Upon reduction of conduction heterogeneity, re-entry was initiated by the rapid pacing and the activation of both atrial chambers was almost simultaneous. Myofiber orientation is an effective mechanism for regulating atrial activation. Modification of myoarchitecture is proarrhythmic

Turbulence transition and the edge of chaos in pipe flow

The linear stability of pipe flow implies that only perturbations of sufficient strength will trigger the transition to turbulence. In order to determine this threshold in perturbation amplitude we study the \emph{edge of chaos} which separates perturbations that decay towards the laminar profile and perturbations that trigger turbulence. Using the lifetime as an indicator and methods developed in (Skufca et al, Phys. Rev. Lett. {\bf 96}, 174101 (2006)) we show that superimposed on an overall $1/\Re$-scaling predicted and studied previously there are small, non-monotonic variations reflecting folds in the edge of chaos. By tracing the motion in the edge we find that it is formed by the stable manifold of a unique flow field that is dominated by a pair of downstream vortices, asymmetrically placed towards the wall. The flow field that generates the edge of chaos shows intrinsic chaotic dynamics.Comment: 4 pages, 5 figure

Turbulence transition in the asymptotic suction boundary layer

We study the transition to turbulence in the asymptotic suction boundary layer (ASBL) by direct numerical simulation. Tracking the motion of trajectories intermediate between laminar and turbulent states we can identify the invariant object inside the laminar-turbulent boundary, the edge state. In small domains, the flow behaves like a travelling wave over short time intervals. On longer times one notes that the energy shows strong bursts at regular time intervals. During the bursts the streak structure is lost, but it reforms, translated in the spanwise direction by half the domain size. Varying the suction velocity allows to embed the flow into a family of flows that interpolate between plane Couette flow and the ASBL. Near the plane Couette limit, the edge state is a travelling wave. Increasing the suction, the travelling wave and a symmetry-related copy of it undergo a saddle-node infinite-period (SNIPER) bifurcation that leads to bursting and discrete-symmetry shifts. In wider domains, the structures localize in the spanwise direction, and the flow in the active region is similar to the one in small domains. There are still periodic bursts at which the flow structures are shifted, but the shift-distance is no longer connected to a discrete symmetry of the flow geometry. Two different states are found by edge tracking techniques, one where structures are shifted to the same side at every burst and one where they are alternatingly shifted to the left and to the right.Comment: Conference TSFP8, Poitiers 2013. TSFP-8 conference proceedings 2013, http://www.tsfp-conference.org/proceedings

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric- field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency - amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnol'd tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.Comment: 9 pages, LaTex, RevTex. 12 uuencoded figures (1.8M) should be requested by e-mail from the autho

Dynamics of Oscillators Coupled by a Medium with Adaptive Impact

In this article we study the dynamics of coupled oscillators. We use mechanical metronomes that are placed over a rigid base. The base moves by a motor in a one-dimensional direction and the movements of the base follow some functions of the phases of the metronomes (in other words, it is controlled to move according to a provided function). Because of the motor and the feedback, the phases of the metronomes affect the movements of the base while on the other hand, when the base moves, it affects the phases of the metronomes in return. For a simple function for the base movement (such as $y = \gamma_{x} [r \theta_1 + (1 - r) \theta_2]$ in which $y$ is the velocity of the base, $\gamma_{x}$ is a multiplier, $r$ is a proportion and $\theta_1$ and $\theta_2$ are phases of the metronomes), we show the effects on the dynamics of the oscillators. Then we study how this function changes in time when its parameters adapt by a feedback. By numerical simulations and experimental tests, we show that the dynamic of the set of oscillators and the base tends to evolve towards a certain region. This region is close to a transition in dynamics of the oscillators; where more frequencies start to appear in the frequency spectra of the phases of the metronomes

On the self-sustained nature of large-scale motions in turbulent Couette flow

Large-scale motions in wall-bounded turbulent flows are frequently interpreted as resulting from an aggregation process of smaller-scale structures. Here, we explore the alternative possibility that such large-scale motions are themselves self-sustained and do not draw their energy from smaller-scale turbulent motions activated in buffer layers. To this end, it is first shown that large-scale motions in turbulent Couette flow at Re=2150 self-sustain even when active processes at smaller scales are artificially quenched by increasing the Smagorinsky constant Cs in large eddy simulations. These results are in agreement with earlier results on pressure driven turbulent channels. We further investigate the nature of the large-scale coherent motions by computing upper and lower-branch nonlinear steady solutions of the filtered (LES) equations with a Newton-Krylov solver,and find that they are connected by a saddle-node bifurcation at large values of Cs. Upper branch solutions for the filtered large scale motions are computed for Reynolds numbers up to Re=2187 using specific paths in the Re-Cs parameter plane and compared to large-scale coherent motions. Continuation to Cs = 0 reveals that these large-scale steady solutions of the filtered equations are connected to the Nagata-Clever-Busse-Waleffe branch of steady solutions of the Navier-Stokes equations. In contrast, we find it impossible to connect the latter to buffer layer motions through a continuation to higher Reynolds numbers in minimal flow units

Influence of cardiac tissue anisotropy on re-entrant activation in computational models of ventricular fibrillation

The aim of this study was to establish the role played by anisotropic diffusion in (i) the number of filaments and epicardial phase singularities that sustain ventricular fibrillation in the heart, (ii) the lifetimes of filaments and phase singularities, and (iii) the creation and annihilation dynamics of filaments and phase singularities. A simplified monodomain model of cardiac tissue was used, with membrane excitation described by a simplified 3-variable model. The model was configured so that a single re-entrant wave was unstable, and fragmented into multiple re-entrant waves. Re-entry was then initiated in tissue slabs with varying anisotropy ratio. The main findings of this computational study are: (i) anisotropy ratio influenced the number of filaments Sustaining simulated ventricular fibrillation, with more filaments present in simulations with smaller values of transverse diffusion coefficient, (ii) each re-entrant filament was associated with around 0.9 phase singularities on the surface of the slab geometry, (iii) phase singularities were longer lived than filaments, and (iv) the creation and annihilation of filaments and phase singularities were linear functions of the number of filaments and phase singularities, and these relationships were independent of the anisotropy ratio. This study underscores the important role played by tissue anisotropy in cardiac ventricular fibrillation