Search

Active-fluid turbulence has been found in bacterial suspensions, but not so far in their algal counterparts. We present the first experimental evidence for turbulence in dense algal suspensions of Chlamydomonas reinhardtii. We carry out a detailed analysis of the statistical properties of the flow present in these cell suspensions and show that they are quantitatively distinct from their counterparts in two-dimensional fluid and bacterial turbulence. Both kinetic-energy and density spectra of the fluid flow in algal turbulence show power-law regimes with unique scaling exponents. The fluid velocity probability distribution function (PDF) is strongly non-Gaussian and the length dependence of the PDF of fluid-velocity increments indicates small-scale intermittency. We compare and contrast our results with recent theoretical predictions for active-scalar turbulence and active glasses. Overall, our results highlight that active turbulence can arise, even in absence of orientational instabilities, so it is not limited to bacterial suspensions but it can also be found in many biological systems with free-swimming micro-organisms., Comment: 9 pages (main article), 5 figures (main article), 3 pages (supplementary information), 3 figures (supplementary information)

Recent results suggest that highly active, chaotic, non-equilibrium states of living fluids might share much in common with high Reynolds number, inertial turbulence. We now show, by using a hydrodynamical model, the onset of intermittency and the consequent multiscaling of Eulerian and Lagrangian structure functions as a function of the bacterial activity. Our results bridge the worlds of low and high Reynolds number flows as well as open up intriguing possibilities of what makes flows intermittent.

We demonstrate the existence of interface-induced turbulence, an emergent nonequilibrium statistically steady state (NESS) with spatiotemporal chaos, which is induced by interfacial fluctuations in low-Reynolds-number binary-fluid mixtures. We uncover the properties of this NESS via direct numerical simulations (DNSs) of cellular flows in the Cahn-Hilliard-Navier-Stokes (CHNS) equations for binary fluids. We show that, in this NESS, the shell-averaged energy spectrum $E(k)$ is spread over several decades in the wavenumber $k$ and it exhibits a power-law region, indicative of turbulence but without a conventional inertial cascade. To characterize the statistical properties of this turbulence, we compute, in addition to $E(k)$, the time series $e(t)$ of the kinetic energy and its power spectrum, scale-by-scale energy transfer as a function of $k$, and the energy dissipation resulting from interfacial stresses. Furthermore, we analyze the mixing properties of this low-Reynolds-number turbulence via the mean-square displacement (MSD) of Lagrangian tracer particles, for which we demonstrate diffusive behavior at long times, a hallmark of strong mixing in turbulent flows., Comment: 9 pages, 6 figures

Energy cascades lie at the heart of the dynamics of turbulent flows. In a recent study of turbulence in fluids with odd-viscosity [de Wit \textit{et al.}, Nature \textbf{627}, 515 (2024)], the two-dimensionalization of the flow at small scales leads to the arrest of the energy cascade and selection of an intermediate scale, between the forcing and the viscous scales. To investigate the generality of this phenomenon, we study a shell model that is carefully constructed to have three-dimensional turbulent dynamics at small wavenumbers and two-dimensional turbulent dynamics at large wavenumbers. The large scale separation that we can achieve in our shell model allows us to examine clearly the interplay between these dynamics, which leads to an arrest of the energy cascade at a transitional wavenumber and an associated accumulation of energy at the same scale. Such pile-up of energy around the transitional wavenumber is reminiscent of the formation of condensates in two-dimensional turbulence, \textit{but, in contrast, it occurs at intermediate wavenumbers instead of the smallest wavenumber, Comment: Typos corrected

We uncover a new type of turbulence -- activity-induced homogeneous and isotropic turbulence in a model that has been employed to investigate motility-induced phase separation (MIPS) in a system of microswimmers. The active Cahn-Hilliard-Navier-Stokes equations (CHNS), also called active model H, provides a natural theoretical framework for our study. In this CHNS model, a single scalar order parameter $\phi$, positive (negative) in regions of high (low) microswimmer density, is coupled with the velocity field $\bm u$. The activity of the microswimmers is governed by an activity parameter $\zeta$ that is positive for \textit{extensile} swimmers and negative for \textit{contractile} swimmers. With extensile swimmers, this system undergoes complete phase separation, which is similar to that in binary-fluid mixtures. By carrying out pseudospectral direct numerical simulations (DNSs), we show, for the first time, that this model (a) develops an emergent nonequilibrium, but statistically steady, state (NESS) of active turbulence, for the case of contractile swimmers, if $\zeta$ is sufficiently large and negative and (b) this turbulence arrests the phase separation. We quantify this suppression by showing how the coarsening-arrest length scale does not grow indefinitely, with time $t$, but saturates at a finite value at large times. We characterise the statistical properties of this active-scalar turbulence by employing energy spectra and fluxes and the spectrum of $\phi$. For sufficiently high Reynolds numbers, the energy spectrum $\mathcal E(k)$ displays an inertial range, with a power-law dependence on the wavenumber $k$. We demonstrate that, in this range, the flux $\Pi(k)$ assumes a nearly constant, negative value, which indicates that the system shows an inverse cascade of energy, even though energy injection occurs over a wide range of wavenumbers in our active-CHNS model., Comment: 8 pages, 5 figures

We carry out a detailed study of dynamic multiscaling in the turbulent nonequilibrium, but statistically steady, state of the stochastically forced one-dimensional Burgers equation. We introduce the concept of $\textit{interval collapse times}$ $\tau_{\rm col}$, the time taken for an interval of length $\ell$, demarcated by a pair of Lagrangian tracers, to collapse at a shock. By calculating the dynamic scaling exponent of the order-$p$ moment of $\tau_{\rm col}$, we show that (a) there is $\textit{not one but an infinity of characteristic time scales}$ and (b) the probability distribution function of $\tau_{\rm col}$ is non-Gaussian and has a power-law tail. Our study is based on (a) a theoretical framework that allows us to obtain dynamic-multiscaling exponents analytically, (b) extensive direct numerical simulations, and (c) a careful comparison of the results of (a) and (b). We discuss possible generalizations of our work to dimensions $d >1 $, for the stochastically forced Burgers equation, and to other compressible flows that exhibit turbulence with shocks.

We use the mean-bacterial-velocity model to investigate the \textit{irreversibility} of two-dimensional (2D) \textit{bacterial turbulence} and to compare it with its 2D fluid-turbulence counterpart. We carry out extensive direct numerical simulations of Lagrangian tracer particles that are advected by the velocity field in this model. Our work uncovers an important, qualitative way in which irreversibility in bacterial turbulence is different from its fluid-turbulence counterpart: For large positive (or large but negative) values of the \textit{friction} (or \textit{activity}) parameter, the probability distribution functions of energy increments, along tracer trajectories, or the power are \textit{positively} skewed; so irreversibility in bacterial turbulence can lead, on average, to \textit{particles gaining energy faster than they lose it}, which is the exact opposite of what is observed for tracers in 2D fluid turbulence.

We use pseudospectral direct numerical simulations (DNSs) to solve the three-dimensional (3D) Hall-Vinen-Bekharevich-Khalatnikov (HVBK) model of superfluid Helium. We then explore the statistical properties of inertial particles, in both coflow and counterflow superfluid turbulence (ST) in the 3D HVBK system; particle motion is governed by a generalization of the Maxey-Riley-Gatignol equations. We first characterize the anisotropy of counterflow ST by showing that there exist large vortical columns. The light particles show confined motion as they are attracted towards these columns and they form large clusters; by contrast, heavy particles are expelled from these vortical regions. We characterise the statistics of such inertial particles in 3D HVBK ST: (1) The mean angle $\Theta(\tau)$, between particle positions, separated by the time lag $\tau$, exhibits two different scaling regions in (a) dissipation and (b) inertial ranges, for different values of the parameters in our model; in particular, the value of $\Theta(\tau)$, at large $\tau$, depends on the magnitude of ${\bf U}_{ns}$. (2) The irreversibility of 3D HVBK turbulence is quantified by computing the statistics of energy increments for inertial particles. (3) The probability distribution function (PDF) of energy increments is of direct relevance to recent experimental studies of irreversibility in superfluid turbulence; we find, in agreement with these experiments, that, for counterflow ST, the skewness of this PDF is less pronounced than its counterparts for coflow ST or for classical-fluid turbulence., Comment: 15 pages, 13 figures

The three-dimensional (3D) Hall magnetohydrodynamics (HMHD) equations are often used to study turbulence in the solar wind. Some earlier studies have investigated the statistical properties of 3D HMHD turbulence by using simple shell models or pseudospectral direct numerical simulations (DNSs) of the 3D HMHD equations; these DNSs have been restricted to modest spatial resolutions and have covered a limited parameter range. To explore the dependence of 3D HMHD turbulence on the Reynolds number $Re$ and the ion-inertial scale $d_{i}$, we have carried out detailed pseudospectral DNSs of the 3D HMHD equations and their counterparts for 3D MHD ($d_{i} = 0$). We present several statistical properties of 3D HMHD turbulence, which we compare with 3D MHD turbulence by calculating (a) the temporal evolution of the energy-dissipation rates and the energy, (b) the wave-number dependence of fluid and magnetic spectra, (c) the probability distribution functions (PDFs) of the cosines of the angles between various pairs of vectors, such as the velocity and the magnetic field, and (d) various measures of the intermittency in 3D HMHD and 3D MHD turbulence.

Infarction- or ischaemia-induced cardiac fibrosis can be arrythmogenic. We use mathematcal models for diffuse fibrosis ($\mathcal{DF}$), interstitial fibrosis ($\mathcal{IF}$), patchy fibrosis ($\mathcal{PF}$), and compact fibrosis ($\mathcal{CF}$) to study patterns of fibrotic cardiac tissue that have been generated by new mathematical algorithms. We show that the fractal dimension $\mathbb{D}$, the lacunarity $\mathcal{L}$, and the Betti numbers $\beta_0$ and $\beta_1$ of such patterns are \textit{fibrotic-tissue markers} that can be used to characterise the arrhythmogenicity of different types of cardiac fibrosis. We hypothesize, and then demonstrate by extensive \textit{in silico} studies of detailed mathematical models for cardiac tissue, that the arrhytmogenicity of fibrotic tissue is high when $\beta_0$ is large and the lacunarity parameter $b$ is small., Comment: Manuscript: 6 pages 7 figures Supplemental Material: 4 pages 6 figures

Ventricular tachycardia (VT) and ventricular fibrillation (VF) are lethal rhythm disorders, which are associated with the occurrence of abnormal electrical scroll waves in the heart. Given the technical limitations of imaging and probing, the in situ visualization of these waves inside cardiac tissue remains a challenge. Therefore, we must, perforce, rely on in-silico simulations of scroll waves in mathematical models for cardiac tissue to develop an understanding of the dynamics of these waves in mammalian hearts. We use direct numerical simulations of the Hund-Rudy-Dynamic (HRD) model, for canine ventricular tissue, to examine the interplay between electrical scroll-waves and conduction and ionic inhomogeneities, in anatomically realistic canine ventricular geometries with muscle-fiber architecture. We find that millimeter-sized, distributed, conduction inhomogeneities cause a substantial decrease in the scroll wavelength, thereby increasing the probability for wave breaks; by contrast, single, localized, medium-sized ($\simeq $ cm) conduction inhomogeneities, exhibit the potential to suppress wave breaks or enable the self-organization of wave fragments into stable, intact scrolls. We show that ionic inhomogeneities, both distributed or localised, suppress scroll-wave break up. The dynamics of a stable rotating wave is not affected significantly by such inhomogeneities, except at high concentrations of distributed inhomogeneities, which can cause a partial break up of scroll waves. Our results indicate that inhomogeneities in the canine ventricular tissue are less arrhythmogenic than inhomogeneities in porcine ventricular tissue, for which an earlier in silico study has shown that the inhomogeneity-induced suppression of scroll waves is a rare occurrence.

We study the spreading of viruses, such as SARS-CoV-2, by airborne aerosols, via a new first-passage-time problem for Lagrangian tracers that are advected by a turbulent flow: By direct numerical simulations of the three-dimensional (3D) incompressible, Navier-Stokes equation, we obtain the time $t_R$ at which a tracer, initially at the origin of a sphere of radius $R$, crosses the surface of the sphere \textit{for the first time}. We obtain the probability distribution function $\mathcal{P}(R,t_R)$ and show that it displays two qualitatively different behaviors: (a) for $R \ll L_{\rm I}$, $\mathcal{P}(R,t_R)$ has a power-law tail $\sim t_R^{-\alpha}$, with the exponent $\alpha = 4$ and $L_{\rm I}$ the integral scale of the turbulent flow; (b) for $l_{\rm I} \lesssim R $, the tail of $\mathcal{P}(R,t_R)$ decays exponentially. We develop models that allow us to obtain these asymptotic behaviors analytically. We show how to use $\mathcal{P}(R,t_R)$ to develop social-distancing guidelines for the mitigation of the spreading of airborne aerosols with viruses such as SARS-CoV-2.

We obtain the von K\'arm\'an-Howarth relation for the stochastically forced three-dimensional Hall-Vinen-Bekharvich-Khalatnikov (3D HVBK) model of superfluid turbulence in Helium ($^4$He) by using the generating-functional approach. We combine direct numerical simulations (DNSs) and analyitcal studies to show that, in the statistically steady state of homogeneous and isotropic superfluid turbulence, in the 3D HVBK model, the probability distribution function (PDF) $P(\gamma)$, of the ratio $\gamma$ of the magnitude of the normal fluid velocity and superfluid velocity, has power-law tails that scale as $P(\gamma) \sim \gamma^3$, for $\gamma \ll 1$, and $P(\gamma) \sim \gamma^{-3}$, for $\gamma \gg 1$. Furthermore, we show that the PDF $P(\theta)$, of the angle $\theta$ between the normal-fluid velocity and superfluid velocity exhibits the following power-law behaviors: $P(\theta)\sim \theta$ for $\theta \ll \theta_*$ and $P(\theta)\sim \theta^{-4}$ for $\theta_* \ll \theta \ll 1$, where $\theta_*$ is a crossover angle that we estimate. From our DNSs we obtain energy, energy-flux, and mutual-friction-transfer spectra, and the longitudinal-structure-function exponents for the normal fluid and the superfluid, as a function of the temperature $T$, by using the experimentally determined mutual-friction coefficients for superfluid Helium $^4$He, so our results are of direct relevance to superfluid turbulence in this system., Comment: 13 pages, 3 figures

Two-dimensional magnetohydrodynamics (2D MHD), forced at (a) large length scales or (b) small length scales, displays turbulent, but statistically steady, states with widely different statistical properties. We present a systematic, comparative study of these two cases (a) and (b) by using direct numerical simulations (DNSs). We find that, in case (a), there is energy equipartition between the magnetic and velocity fields, whereas, in case (b), such equipartition does not exist. By computing various probability distribution functions (PDFs), we show that case (a) displays extreme events that are much less common in case (b).

We carry out extensive direct numerical simulations (DNSs) to investigate the interaction of active particles and fields in the two-dimensional (2D) Gross-Pitaevskii (GP) superfluid, in both simple and turbulent flows. The particles are active in the sense that they affect the superfluid even as they are affected by it. We tune the mass of the particles, which is an important control parameter. At the one-particle level, we show how light, neutral, and heavy particles move in the superfluid, when a constant external force acts on them; in particular, beyond a critical velocity, at which a vortex-antivortex pair is emitted, particle motion can be periodic or chaotic. We demonstrate that the interaction of a particle with vortices leads to dynamics that depends sensitively on the particle characteristics. We also demonstrate that assemblies of particles and vortices can have rich, and often turbulent spatiotemporal evolution. In particular, we consider the dynamics of the following illustrative initial configurations: (a) one particle placed in front of a translating vortex-antivortex pair; (b) two particles placed in front of a translating vortex-antivortex pair; (c) a single particle moving in the presence of counter-rotating vortex clusters; and (d) four particles in the presence of counter-rotating vortex clusters. We compare our work with earlier studies and examine its implications for recent experimental studies in superfluid Helium and Bose-Einstein condensates., Comment: 24 figures

The Beale-Kato-Majda theorem contains a single criterion that controls the behaviour of solutions of the $3D$ incompressible Euler equations. Versions of this theorem are discussed in terms of the regularity issues surrounding the $3D$ incompressible Euler and Navier-Stokes equations together with a phase-field model for the statistical mechanics of binary mixtures called the $3D$ Cahn-Hilliard-Navier-Stokes (CHNS) equations. A theorem of BKM-type is established for the CHNS equations for the full parameter range. Moreover, for this latter set, it is shown that there exists a Reynolds number and a bound on the energy-dissipation rate that, remarkably, reproduces the $Re^{3/4}$ upper bound on the inverse Kolmogorov length normally associated with the Navier-Stokes equations alone. An alternative length-scale is introduced and discussed, together with a set of pseudo-spectral computations on a $128^{3}$ grid., Comment: 3 figures and 3 tables

Replicator equation -- a paradigm equation in evolutionary game dynamics -- mathematizes the frequency dependent selection of competing strategies vying to enhance their fitness (quantified by the average payoffs) with respect to the average fitnesses of the evolving population under consideration. In this paper, we deal with two discrete versions of the replicator equation employed to study evolution in a population where any two players, interaction is modeled by a two-strategy symmetric normal-form game. There are twelve distinct classes of such games, each typified by a particular ordinal relationship among the elements of the corresponding payoff matrix. Here we find the sufficient conditions for the existence of asymptotic solutions of the replicator equations such that the solutions -- fixed points, periodic orbits, and chaotic trajectories -- are all strictly physical, meaning that the frequency of any strategy lies inside the closed interval zero to one at all times. Thus, we elaborate which of the twelve types of games are capable of showing meaningful physical solutions and for which of the two types of replicator equation. Subsequently, we introduce the concept of the weight of fitness deviation that is the scaling factor in a positive affine transformation connecting two payoff matrices such that the corresponding one-shot games have exactly same Nash equilibria and evolutionary stable states. The weight also quantifies how much the excess of fitness of a strategy over the average fitness of the population affects the per capita change in the frequency of the strategy. Intriguingly, the weight's variation is capable of making the Nash equilibria and the evolutionary stable states useless by introducing strict physical chaos in the replicator dynamics based on the normal-form game., Comment: Published in Chaos. This version has some typos corrected

We present the first study of the multiscaling of time-dependent velocity and magnetic-field structure functions in homogeneous, isotropic magnetohydrodynamic (MHD) turbulence in three dimensions. We generalize the formalism that has been developed for analogous studies of time-dependent structure functions in fluid turbulence to MHD. By carrying out detailed numerical studies of such time-dependent structure functions in a shell model for three-dimensional MHD turbulence, we obtain both equal-time and dynamic scaling exponents., Comment: 9 pages; 5 figures; Published in Phys. Rev. E

We build on recent developments in the study of fluid turbulence [Gibbon \textit{et al.} Nonlinearity 27, 2605 (2014)] to define suitably scaled, order-$m$ moments, $D_m^{\pm}$, of $\omega^\pm= \omega \pm j$, where $\omega$ and $j$ are, respectively, the vorticity and current density in three-dimensional magnetohydrodynamics (MHD). We show by mathematical analysis, for unit magnetic Prandtl number $P_M$, how these moments can be used to identify three possible regimes for solutions of the MHD equations; these regimes are specified by inequalities for $D_m^{\pm}$ and $D_1^{\pm}$. We then compare our mathematical results with those from our direct numerical simulations (DNSs) and thus demonstrate that 3D MHD turbulence is like its fluid-turbulence counterpart insofar as all solutions, which we have investigated, remain in \textit{only one of these regimes}; this regime has depleted nonlinearity. We examine the implications of our results for the exponents $q^{\pm}$ that characterize the power-law dependences of the energy spectra $\mathcal{E}^{\pm}(k)$ on the wave number $k$, in the inertial range of scales. We also comment on (a) the generalization of our results to the case $P_M \neq 1$ and (b) the relation between $D_m^{\pm}$ and the order-$m$ moments of gradients of hydrodynamic fields, which are used in characterizing intermittency in turbulent flows., Comment: 14 pages, 3 figures

We examine the multiscaling behavior of the normal- and superfluid-velocity structure functions in three-dimensional superfluid turbulence by using a shell model for the three-dimensional (3D) Hall-Vinen-Bekharevich-Khalatnikov (HVBK) equations. Our 3D-HVBK shell model is based on the Gledzer-Okhitani-Yamada (GOY) shell model. We examine the dependence of the multiscaling exponents on the normal-fluid fraction and the mutual-friction coefficients. Our extensive study of the 3D-HVBK shell model shows that the multiscaling behavior of the velocity structure functions in superfluid turbulence is more complicated than it is in fluid turbulence., Comment: 12 pages, 6 figures

We uncover universal statistical properties of the trajectories of heavy inertial particles in three-dimensional, statistically steady, homogeneous, and isotropic turbulent flows by extensive direct numerical simulations. We show that the probability distribution functions (PDFs) $P(\phi)$, of the angle $\phi$ between the Eulerian velocity ${\bf u}$ and the particle velocity ${\bf v}$, at this point and time, shows a power-law region in which $P(\phi) \sim \phi^{-\gamma}$, with a new universal exponent $\gamma \simeq 4$. Furthermore, the PDFs of the trajectory curvature $\kappa$ and modulus $\theta$ of the torsion $\vartheta$ have power-law tails that scale, respectively, as $P(\kappa) \sim \kappa^{-h_\kappa}$, as $\kappa \to \infty$, and $P(\theta) \sim \theta^{-h_\theta}$, as $\theta \to \infty$, with exponents $h_\kappa \simeq 2.5$ and $h_\theta \simeq 3$ that are universal to the extent that they do not depend on the Stokes number ${\rm St}$ (given our error bars). We also show that $\gamma$, $h_\kappa$ and $h_\theta$ can be obtained by using simple stochastic models. We characterize the complexity of heavy-particle trajectories by the number $N_{\rm I}(t,{\rm St})$ of points (up until time $t$) at which $\vartheta$ changes sign. We show that $n_{\rm I}({\rm St}) \equiv \lim_{t\to\infty} \frac{N_{\rm I}(t,{\rm St})}{t} \sim {\rm St}^{-\Delta}$, with $\Delta \simeq 0.4$ a universal exponent.

We present algorithms for the ab-initio determination of the temperature ($T$) dependence of the mutual-friction coefficients $\alpha$ and $\alpha'$ and the normal-fluid density $\rho_{\rm n}$ in the two-dimensional (2D) Galerkin-truncated Gross-Pitaevskii system. Our algorithms enable us to determine $\alpha(T)$, even though fluctuations in 2D are considerably larger than they are in 3D. We also examine the implications of our measurements of $\alpha'(T)$ for the Iordanskii force, whose existence is often questioned., Comment: For associated videos, see http://youtu.be/yUxRhLDeGcI and http://youtu.be/gMp_Rj_aMns

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normal-fluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli. Our study leads to specific suggestions for experiments.

We present a detailed direct numerical simulation of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic (2D MHD) turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work sig- nificantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distin- guishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream-function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies.

The periodic $3D$ Navier-Stokes equations are analyzed in terms of dimensionless, scaled, $L^{2m}$-norms of vorticity $D_{m}$ ($1 \leq m < \infty$). The first in this hierarchy, $D_{1}$, is the global enstrophy. Three regimes naturally occur in the $D_{1}-D_{m}$ plane. Solutions in the first regime, which lie between two concave curves, are shown to be regular, owing to strong nonlinear depletion. Moreover, numerical experiments have suggested, so far, that all dynamics lie in this heavily depleted regime \cite{DGGKPV13}\,; new numerical evidence for this is presented. Estimates for the dimension of a global attractor and a corresponding inertial range are given for this regime. However, two more regimes can theoretically exist. In the second, which lies between the upper concave curve and a line, the depletion is insufficient to regularize solutions, so no more than Leray's weak solutions exist. In the third, which lies above this line, solutions are regular, but correspond to extreme initial conditions. The paper ends with a discussion on the possibility of transition between these regimes., Comment: 9 figures and 3 tables

We study the statistical properties of orientation and rotation dynamics of elliptical tracer particles in two-dimensional, homogeneous and isotropic turbulence by direct numerical simulations. We consider both the cases in which the turbulent flow is generated by forcing at large and intermediate length scales. We show that the two cases are qualitatively different. For large-scale forcing, the spatial distribution of particle orientations forms large-scale structures, which are absent for intermediate-scale forcing. The alignment with the local directions of the flow is much weaker in the latter case than in the former. For intermediate-scale forcing, the statistics of rotation rates depends weakly on the Reynolds number and on the aspect ratio of particles. In contrast with what is observed in three-dimensional turbulence, in two dimensions the mean-square rotation rate decreases as the aspect ratio increases., Comment: 5 pages, 6 figures

We show that a shell-model version of the three-dimensional Hall-magnetohydrodynamic (3D Hall-MHD) equations provides a natural theoretical model for investigating the multiscaling behaviors of velocity and magnetic structure functions. We carry out extensive numerical studies of this shell model, obtain the scaling exponents for its structure functions, in both the low-$k$ and high-$k$ power-law ranges of 3D Hall-MHD, and find that the extended-self-similarity (ESS) procedure is helpful in extracting the multiscaling nature of structure functions in the high-$k$ regime, which otherwise appears to display simple scaling. Our results shed light on intriguing solar-wind measurements., Comment: 7 pages, 6 figures

The issue of intermittency in numerical solutions of the 3D Navier-Stokes equations on a periodic box $[0,\,L]^{3}$ is addressed through four sets of numerical simulations that calculate a new set of variables defined by $D_{m}(t) = \left(\varpi_{0}^{-1}\Omega_{m}\right)^{\alpha_{m}}$ for $1 \leq m \leq \infty$ where $\alpha_{m}= \frac{2m}{4m-3}$ and $\left[\Omega_{m}(t)\right]^{2m} = L^{-3}\I |\bom|^{2m}dV$ with $\varpi_{0} = \nu L^{-2}$. All four simulations unexpectedly show that the $D_{m}$ are ordered for $m = 1\,,...,\,9$ such that $D_{m+1} < D_{m}$. Moreover, the $D_{m}$ squeeze together such that $D_{m+1}/D_{m}\nearrow 1$ as $m$ increases. The first simulation is of very anisotropic decaying turbulence\,; the second and third are of decaying isotropic turbulence from random initial conditions and forced isotropic turbulence at constant Grashof number respectively\,; the fourth is of very high Reynolds number forced, stationary, isotropic turbulence at up to resolutions of $4096^{3}$.

We undertake a systematic, direct numerical simulation (DNS) of the two-dimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. First, there are transients that depend on the initial conditions. In the second regime, power-law scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power-laws and the extents of the scaling regions change with time and depended on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers $k > k_c$ and partial thermalization takes place for modes with $k < k_c$; the self-truncation wave number $k_c(t)$ depends on the initial conditions and it grows either as a power of $t$ or as $\log t$. Finally, in the fourth regime, complete-thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms., Comment: 30 pages, 12 figures

An energy-spectrum bottleneck, a bump in the turbulence spectrum between the inertial and dissipation ranges, is shown to occur in the non-turbulent, one-dimensional, hyperviscous Burgers equation and found to be the Fourier-space signature of oscillations in the real-space velocity, which are explained by boundary-layer-expansion techniques. Pseudospectral simulations are used to show that such oscillations occur in velocity correlation functions in one- and three-dimensional hyperviscous hydrodynamical equations that display genuine turbulence., Comment: 5 pages, 2 figures

We obtain, by extensive direct numerical simulations, time-dependent and equal-time structure functions for the vorticity, in both quasi-Lagrangian and Eulerian frames, for the direct-cascade regime in two-dimensional fluid turbulence with air-drag-induced friction. We show that differ- ent ways of extracting time scales from these time-dependent structure functions lead to different dynamic-multiscaling exponents, which are related to equal-time multiscaling exponents by different classes of bridge relations; for a representative value of the friction we verify that, given our error bars, these bridge relations hold., Comment: Published in Phys. Rev. Lett, Supplemental materials added as appendix

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter $\Lambda$ to distinguish between vortical and extensional regions. We then use a direct numerical simulation (DNS) of the two-dimensional, incompressible Navier--Stokes equation with Ekman friction to study probability distribution functions (PDFs) of the persistence times of vortical and extensional regions by employing both Eulerian and Lagrangian measurements. We find that, in the Eulerian case, the persistence-time PDFs have exponential tails; by contrast, this PDF for Lagrangian particles, in vortical regions, has a power-law tail with an exponent $\theta=2.9\pm0.2$., Comment: consistent with the published version

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, isotropic fluid turbulence. We find clear manifestations of dissipation-reduction phenomena: On the addition of polymers to the turbulent fluid, we obtain a reduction in the energy dissipation rate, a significant modification of the fluid energy spectrum, especially in the deep-dissipation range, a suppression of small-scale intermittency, and a decrease in small-scale vorticity filaments. We also compare our results with recent experiments and earlier DNS studies of decaying fluid turbulence with polymer additives., Comment: consistent with the published version

We present the results of our detailed pseudospectral direct numerical simulation (DNS) studies, with up to $1024^3$ collocation points, of incompressible, magnetohydrodynamic (MHD) turbulence in three dimensions, without a mean magnetic field. Our study concentrates on the dependence of various statistical properties of both decaying and statistically steady MHD turbulence on the magnetic Prandtl number ${\rm Pr_M}$ over a large range, namely, $0.01 \leq {\rm Pr_M} \leq 10$. We obtain data for a wide variety of statistical measures such as probability distribution functions (PDFs) of moduli of the vorticity and current density, the energy dissipation rates, and velocity and magnetic-field increments, energy and other spectra, velocity and magnetic-field structure functions, which we use to characterise intermittency, isosurfaces of quantities such as the moduli of the vorticity and current, and joint PDFs such as those of fluid and magnetic dissipation rates. Our systematic study uncovers interesting results that have not been noted hitherto. In particular, we find a crossover from larger intermittency in the magnetic field than in the velocity field, at large ${\rm Pr_M}$, to smaller intermittency in the magnetic field than in the velocity field, at low ${\rm Pr_M}$. Furthermore, a comparison of our results for decaying MHD turbulence and its forced, statistically steady analogue suggests that we have strong universality in the sense that, for a fixed value of ${\rm Pr_M}$, multiscaling exponent ratios agree, at least within our errorbars, for both decaying and statistically steady homogeneous, isotropic MHD turbulence., Comment: 49 pages,33 figures

We carry out systematic and high-resolution studies of dynamo action in a shell model for magnetohydrodynamic (MHD) turbulence over wide ranges of the magnetic Prandtl number $Pr_{\rm M}$ and the magnetic Reynolds number $Re_{\rm M}$. Our study suggests that it is natural to think of dynamo onset as a nonequilibrium, first-order phase transition between two different turbulent, but statistically steady, states. The ratio of the magnetic and kinetic energies is a convenient order parameter for this transition. By using this order parameter, we obtain the stability diagram (or nonequilibrium phase diagram) for dynamo formation in our MHD shell model in the $(Pr^{-1}_{\rm M}, Re_{\rm M})$ plane. The dynamo boundary, which separates dynamo and no-dynamo regions, appears to have a fractal character. We obtain hysteretic behavior of the order parameter across this boundary and suggestions of nucleation-type phenomena., Comment: 9 pages, 9 figures

To develop an understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that, as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially the temporal oscillations are periodic; this periodic behaviour becomes more and more complicated, with increasing Reynolds number, until the film enters a spatially disordered nonequilibrium statistical steady that is turbulent. We study this sequence of transitions by using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g., \Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and which lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterise short- and long-range order in the turbulent and crystalline phases, respectively., Comment: 31 pages, 56 figures, movie files not included

We present an introductory overview of several challenging problems in the statistical characterisation of turbulence. We provide examples from fluid turbulence in three and two dimensions, from the turbulent advection of passive scalars, turbulence in the one-dimensional Burgers equation, and fluid turbulence in the presence of polymer additives., Comment: 34 pages, 31 figures

We present a detailed direct numerical simulation (DNS) designed to investigate the combined effects of walls and Ekman friction on turbulence in forced soap films. We concentrate on the forward-cascade regime and show how to extract the isotropic parts of velocity and vorticity structure functions and thence the ratios of multiscaling exponents. We find that velocity structure functions display simple scaling whereas their vorticity counterparts show multiscaling; and the probability distribution function of the Weiss parameter $\Lambda$, which distinguishes between regions with centers and saddles, is in quantitative agreement with experiments., Comment: 4 pages, 6 figures

It is shown that the use of a high power $\alpha$ of the Laplacian in the dissipative term of hydrodynamical equations leads asymptotically to truncated inviscid \textit{conservative} dynamics with a finite range of spatial Fourier modes. Those at large wavenumbers thermalize, whereas modes at small wavenumbers obey ordinary viscous dynamics [C. Cichowlas et al. Phys. Rev. Lett. 95, 264502 (2005)]. The energy bottleneck observed for finite $\alpha$ may be interpreted as incomplete thermalization. Artifacts arising from models with $\alpha > 1$ are discussed., Comment: 4 pages, 2 figures, Phys. Rev. Lett. in press

We systematise the study of dynamic multiscaling of time-dependent structure functions in different models of passive-scalar and fluid turbulence. We show that, by suitably normalising these structure functions, we can eliminate their dependence on the origin of time at which we start our measurements and that these normalised structure functions yield the same linear bridge relations that relate the dynamic-multiscaling and equal-time exponents for statistically steady turbulence. We show analytically, for both the Kraichnan Model of passive-scalar turbulence and its shell model analogue, and numerically, for the GOY shell model of fluid turbulence and a shell model for passive-scalar turbulence, that these exponents and bridge relations are the same for statistically steady and decaying turbulence. Thus we provide strong evidence for dynamic universality, i.e., dynamic-multiscaling exponents do not depend on whether the turbulence decays or is statistically steady., Comment: 29 pages, 18 figures. For figures with better resolution, please download the preprint available at the website of the International Collaboration of Turbulence Research

We give an overview of the progress that has been made in recent years in understanding the dynamic multiscaling of homogeneous, isotropic turbulence and related problems. We emphasise the similarity of this problem with the dynamic scaling of time-dependent correlation functions in the vicinity of a critical point in, e.g., a spin system. The universality of dynamic-multiscaling exponents in fluid turbulence is explored by detailed simulations of the GOY shell model for fluid turbulence., Comment: 13 pages, 10 figures, Submitted for the Proceedings of the Statphys 23 Conference

Cardiac arrhythmias such as ventricular tachycardia (VT) or ventricular fibrillation (VF) are the leading cause of death in the industrialised world. There is a growing consensus that these arrhythmias arise because of the formation of spiral waves of electrical activation in cardiac tissue; unbroken spiral waves are associated with VT and broken ones with VF. Several experimental studies have been carried out to determine the effects of inhomogeneities in cardiac tissue on such arrhythmias. We give a brief overview of such experiments, and then an introduction to partial-differential-equation models for ventricular tissue. We show how different types of inhomogeneities can be included in such models, and then discuss various numerical studies, including our own, of the effects of these inhomogeneities on spiral-wave dynamics. The most remarkable qualitative conclusion of our studies is that the spiral-wave dynamics in such systems depends very sensitively on the positions of these inhomogeneities., Comment: 19 pages, 11 figures. To be published in the Proceedings of the Conference on Complex Dynamics of Physiological Systems: From Heart to Brain held at Presidency College, Kolkata, India, 12-14 February 2007; to appear in the Life Science Series of Springer, The Netherlands

The existence of drag reduction by polymer additives, well established for wall-bounded turbulent flows, is controversial in homogeneous, isotropic turbulence. To settle this controversy we carry out a high-resolution direct numerical simulation (DNS) of decaying, homogeneous, isotropic turbulence with polymer additives. Our study reveals clear manifestations of drag-reduction-type phenomena: On the addition of polymers to the turbulent fluid we obtain a reduction in the energy dissipation rate, a significant modification of the fluid energy spectrum especially in the deep-dissipation range, a suppression of small-scale intermittency, and a decrease in small-scale vorticity filaments., Comment: 4 pages, 3 figures

We present the first study of the dynamic scaling or multiscaling of passive-scalar and passive-vector turbulence. For the Kraichnan version of passive-scalar and passive-vector turbulence we show analytically, in both Eulerian and quasi-Lagrangian frameworks, that simple dynamic scaling is obtained but with different dynamic exponents. By developing the multifractal model we show that dynamic multiscaling occurs in passive-scalar turbulence only if the advecting velocity field is itself multifractal. We substantiate our results by detailed numerical simulations in shell models of passive-scalar advection., Comment: published version

We present results from a systematic numerical study of decaying turbulence in a dilute polymer solution by using a shell-model version of the FENE-P equations. Our study leads to an appealing definition of drag reduction for the case of decaying turbulence. We exhibit several new results, such as the potential-energy spectrum of the polymer, hitherto unobserved features in the temporal evolution of the kinetic-energy spectrum, and characterize intermittency in such systems. We compare our results with the GOY shell model for fluid turbulence., Comment: 4 pages RevTex4, 7 figures, accepted for publication in Physical Review E

We study turbulence in the one-dimensional Burgers equation with a white-in-time, Gaussian random force that has a Fourier-space spectrum $\sim 1/k$, where $k$ is the wave number. From very-high-resolution numerical simulations, in the limit of vanishing viscosity, we find evidence for multiscaling of velocity structure functions which cannot be falsified by standard tests. We find a new artifact in which logarithmic corrections can appear disguised as anomalous scaling and conclude that bifractal scaling is likely., Comment: previous version replaced, with major changes

We show that different ways of extracting time scales from time-dependent velocity structure functions lead to different dynamic-multiscaling exponents in fluid turbulence. These exponents are related to equal-time multiscaling exponents by different classes of bridge relations which we derive. We check this explicitly by detailed numerical simulations of the GOY shell model for fluid turbulence. Our results can be generalized to any system in which both equal-time and time-dependent structure functions show multiscaling., Comment: 4 pages 4 figures

We derive relations for the decay of the kinetic and magnetic energies and the growth of the Taylor and integral scales in unforced, incompressible, homogeneous and isotropic three-dimensional magnetohydrodynamic (3DMHD) turbulence with power-law initial energy spectra. We also derive bounds for the decay of the cross- and magnetic helicities. We then present results from systematic numerical studies of such decay both within the context of an MHD shell model and direct numerical simulations (DNS) of 3DMHD. We show explicitly that our results about the power-law decay of the energies hold for times $t

We consider the effects of a shear on velocity fluctuations in a flow. The shear gives rise to a transient amplification that not only influences the amplitude of perturbations but also their time correlations. We show that, in the presence of white noise, time correlations of transversal velocity components are exponential and that correlations of the longitudinal components are exponential with an algebraic prefactor. Cross correlations between transversal and downstream components are strongly asymmetric and provide a clear indication of non-normal amplification. We suggest experimental tests of our predictions., Comment: 6 pages, 2 figures, gzipped tar file

Ventricular fibrillation, the major reason behind sudden cardiac death, is turbulent cardiac electrical activity in which rapid, irregular disturbances in the spatiotemporal electrical activation of the heart makes it incapable of any concerted pumping action. Methods of controlling ventricular fibrillation include electrical defibrillation as well as injected medication. Electrical defibrillation, though widely used, involves subjecting the whole heart to massive, and often counterproductive, electrical shocks. We propose a defibrillation method that uses a very-low-amplitude shock (of order mV) applied for a brief duration (of order 100 ms) and over a coarse mesh of lines on our model ventricle., Comment: 11 pages revtex, 2 postscript figures, 1 jpeg figure