Wave propagation in ventricular muscle is rendered highly anisotropic by the intramural rotation of the fiber This rotational anisotropy is especially important because it can produce a twist of electrical vortices, which measures the rate of rotation (in degree/mm) of activation wavefronts in successive planes perpendicular to a line of phase singularity, or filament This twist can then significantly alter the dynamics of the filament This paper explores this dynamics via numerical simulation After a review of the literature, we present modeling tools that include: (i) a simplified ionic model with three membrane currents that approximates well the restitution properties and spiral wave behavior of more complex ionic models of cardiac action potential (Beeler-Reuter and others), and (ii) a semi-implicit algorithm for the fast solution of monodomain cable equations with rotational anisotropy We then discuss selected results of a simulation study of vortex dynamics in a parallelepipedal slab of ventricular muscle of varying wall thickness (S) and fiber rotation rate (theta(z)) The main finding is that rotational anisotropy generates a sufficiently large twist to destabilize a single transmural filament and cause a transition to a wave turbulent state characterized by a high density of chaotically moving filaments This instability is manifested by the propagation of localized disturbances along the filament and has no previously known analog in isotropic excitable media These disturbances correspond to highly twisted and distorted regions of filament, or "twistons," that create vortex rings when colliding with the natural boundaries of the ventricle Moreover, when sufficiently twisted, these rings expand and create additional filaments by further colliding with boundaries This instability mechanism is distinct from the commonly invoked patchy failure or wave breakup that is not observed here during the initial instability For modified Beeler-Reuter-like kinetics with stable reentry in two dimensions, decay into turbulence occurs in the left ventricle in about one second above a critical wall thickness in the range of 4-6 mm that matches experiment However this decay is suppressed by uniformly decreasing excitability Specific experiments to test these results, and a method to characterize the filament density during fibrillation are discussed Results are contrasted with other mechanisms of fibrillation and future prospects are summarized (c)1998 American Institute of Physics

Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation.

A variety of observations of intense, long-lived oceanic vortices are interpreted as examples of a distinct phenomenon, which is given the name Submesoscale, Coherent Vortices (SCV's). The distinguishing characteristics of SCV's are defined and illustrated by example, and a survey is made of the different SCV types presently known. On the basis of extant theoretical and modeling solutions, interpretations are made of the dynamics associated with SCV existence, movement, endurance, interactions with other currents, generation, and contributions to the transport of chemical properties in the ocean.

Submesoscale, coherent vortices in the ocean

Forced turbulence in a rotating frame is studied using numerical simulations in a triply periodic box. The random forcing is three dimensional and localized about an intermediate wavenumber kf. The results show that energy is transferred to scales larger than the forcing scale when the rotation rate is large enough. The scaling of the energy spectrum approaches E(k)∝k−3 for k<kf. Almost all of the energy for k<kf lies in the two-dimensional (2D) plane perpendicular to the rotation z-axis, and thus the large-scale motions are quasi-2D with E(k)≈E(kh,kz=0), where kh and kz are, respectively, the horizontal and vertical components of the wavevector. The large scales consist of cyclonic vortices. Possible mechanisms responsible for the two-dimensionalization are discussed. The development of the 2D spectrum E(kh,kz=0)∝kh−3 is analogous to the dynamics of β-plane turbulence leading to the Rhines spectrum E(ky,kx=0)∝ky−5.

Transfer of energy to two-dimensional large scales in forced, rotating three-dimensional turbulence

A quantitative experimental study of the two-dimensional inverse energy cascade is presented. The flow is electrically driven in a horizontal layer of mercury and three-dimensional perturbations are suppressed by means of a uniform magnetic field, so that the flow can be well approximated by a two-dimensional Navier–Stokes equation with a steady forcing term and a linear friction due to the Hartmann layer. Turbulence is produced by the instability of a periodic square network of 36 electrically driven alternating vortices. The inverse cascade is limited at large scales, either by the linear friction or by the finite size of the domain, depending on the experimental parameters. In the first case, spectra are measured and the corresponding two-dimensional Kolmogorov constant is in the range 3–7. In the second case, a condensation of the turbulent energy in the lowest mode, corresponding to a spontaneous mean global rotation, is observed. Such a condensation was predicted by Kraichnan (1967) from statistical thermodynamics arguments, but without the symmetry breaking. Random reversals of the rotation sense, owing to turbulent fluctuations, are more and more sparse as friction is decreased. The lowest mode fluctuations and the small scales are statistically independent.

Experimental study of the two-dimensional inverse energy cascade in a square box

The structures of the vorticity fields in several homogeneous irrotational straining flows and a homogeneous turbulent shear flow were examined using a database generated by direct numerical simulation of the unsteady Navier-Stokes equations. In all cases, strong evidence was found for the presence of coherent vortical structures. The initially isotropic vorticity fields were rapidly affected by imposed mean strain and the rotational component of mean shear and developed accordingly. In the homogeneous turbulent shear-flow cases, the roll-up of mean vorticity into characteristic hairpin vortices was clearly observed, supporting the view that hairpin vortices are an important vortical structure in all turbulent shear flows; the absence of mean shear in the homogeneous irrotational straining flows precludes the presence of hairpin-like vortices.

The structure of the vorticity field in homogeneous turbulent flows

A turbulent field is produced with an oscillating grid in a deep, rotating tank. Near the grid, the Rossby number is kept large, 0(3-33), and the turbulence is locally unaffected by rotation. Away from the grid, the scale of the turbulence increases, the r.m.s. turbulent velocity decreases, and rotation becomes increasingly important. The flow field changes dramatically at a local Rossby number of about 0.20, and thereafter remains independent of depth. The flow consists of concentrated vortices having axes approximately parallel to the rotation axis, and extending throughout the depth of the fluid above the turbulent Ekman layer. The number of vortices per unit area is a function of the grid Rossby number. The local vorticity within cores can be a factor of 50 larger than the tank vorticity 2Ω. The total relative circulation contained in the vortices remains, however, a small fraction of the tank circulation.The concentrated vortex cores support waves consisting of helical distortions, which travel along the axes of individual vortices. Isolated, travelling waves seem well-described by the vortex-soliton theory of Hasimoto (1972). The nonlinear waves transport mass, momentum and energy from the vigorously turbulent region near the grid to the rotation-dominated flow above. Interactions between waves, which are frequent occurrences, almost always result in a local breakdown of the vortex core, and small-scale turbulence production. Usually the portions of broken core reform within ½−1 rotation periods, but occasionally cores are destroyed and reformed on a much longer timescale.

https://www.cambridge.org/core/services/aop-cambridge-core/content/view/S0022112082003462

Y. Gagne

Papers

Turbulence and waves in a rotating tank