BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Rahul Pandit

Rahul Pandit appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynnlin.CDcond-mat.stat-mechcond-mat.quant-gasBiological Physicscond-mat.softastro-ph.SRTissues and Organsastro-ph.GAastro-ph.HEchao-dyncond-matmath-phmath.MPnlin.PSphysics.comp-phquant-ph

Catalog footprint

What is connected

43works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2023arXiv

Inertial Particles in Superfluid Turbulence: Coflow and Counterflow

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 $Θ(τ)$, between particle positions, separated by the time lag $τ$, 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 $Θ(τ)$, at large $τ$, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Activity-induced droplet propulsion and multifractality

We develop a minimal hydrodynamic model, without an orientational order parameter, for assemblies of contractile swimmers encapsulated in a droplet of a binary-fluid emulsion. Our model uses two coupled scalar order parameters, $ϕ$ and $ψ$, which capture, respectively, the droplet interface and the activity of the contractile swimmers inside this droplet. These order parameters are also coupled to the velocity field $\bm u$. At low activity, our model yields a self-propelling droplet whose center of mass $(CM)$ displays rectilinear motion, powered by the spatiotemporal evolution of the field $ψ$, which leads to a time-dependent vortex dipole at one end of the droplet. As we increase the activity, this $CM$ shows chaotic super-diffusive motion, which we characterize by its mean-square displacement; and the droplet interface exhibits multifractal fluctuations, whose spectrum of exponents we calculate. We explore the implications of our results for experiments on active droplets of contractile swimmers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamic multiscaling in stochastically forced Burgers turbulence

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}$ $τ_{\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 $τ_{\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 $τ_{\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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Insights from a pseudospectral study of a potentially singular solution of the three-dimensional axisymmetric incompressible Euler equation

We develop a Fourier-Chebyshev pseudospectral direct numerical simulation (DNS) to examine a potentially singular solution of the radially bounded, three-dimensional (3D), axisymmetric Euler equations [G. Luo and T.Y. Hou, Proc. Natl. Acad. Sci. USA, 111.36 (2014)]. We demonstrate that: (a) the time of singularity is preceded, in any spectrally truncated DNS, by the formation of oscillatory structures called tygers, first investigated in the one-dimensional (1D) Burgers and two-dimensional (2D) Euler equations; (b) the analyticity-strip method can be generalized to obtain an estimate for the (potential) singularity time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Irreversiblity in Bacterial Turbulence: Insights from the Mean-Bacterial-Velocity Model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rotating self-gravitating Bose-Einstein condensates with a crust: a minimal model for pulsar glitches

We develop a minimal model for \textit{pulsar glitches} by introducing a solid-crust potential in the three-dimensional (3D) Gross-Pitaevskii-Poisson equation (GPPE), which we have used earlier to study gravitationally bound Bose-Einstein Condensates (BECs), i.e., bosonic stars. In the absence of the crust potential, we show that, if we rotate such a bosonic star, it is threaded by vortices. We then show, via extensive direct numerical simulations (DNSs), that the interaction of these vortices with the crust potential yields (a) stick-slip dynamics and (b) dynamical glitches. We demonstrate that, if enough momentum is transferred to the crust from the bosonic star, then the vortices are expelled from the star and the crust's angular momentum $J_c$ exhibits features that can be interpreted naturally as glitches. From the time series of $J_c$, we compute the cumulative probability distribution functions (CPDFs) of event sizes, event durations, and waiting times. We show that these CPDFs have signatures of self-organized criticality (SOC), which have been seen in observations on pulsar glitches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The spreading of viruses by airborne aerosols: lessons from a first-passage-time problem for tracers in turbulent flows

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^{-α}$, with the exponent $α= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Deep-learning-assisted detection and termination of spiral- and broken-spiral waves in mathematical models for cardiac tissue

Unbroken and broken spiral waves, in partial-differential-equation (PDE) models for cardiac tissue, are the mathematical analogs of life-threatening cardiac arrhythmias, namely, ventricular tachycardia (VT) and ventricular-fibrillation (VF). We develop a (a) deep-learning method for the detection of unbroken and broken spiral waves and (b) the elimination of such waves, e.g., by the application of low-amplitude control currents in the cardiac-tissue context. Our method is based on a convolutional neural network (CNN) that we train to distinguish between patterns with spiral waves S and without spiral waves NS. We obtain these patterns by carrying out extensive direct numerical simulations (DNS) of PDE models for cardiac tissue in which the transmembrane potential V, when portrayed via pseudocolor plots, displays patterns of electrical activation of types S and NS. We then utilize our trained CNN to obtain, for a given pseudocolor image of V, a heat map that has high intensity in the regions where this image shows the cores of spiral waves. Given this heat map, we show how to apply low-amplitude Gaussian current pulses to eliminate spiral waves efficiently. Our in silico results are of direct relevance to the detection and elimination of these arrhythmias because our elimination of unbroken or broken spiral waves is the mathematical analog of low-amplitude defibrillation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The one-dimensional Kardar-Parisi-Zhang and Kuramoto-Sivashinsky universality class: limit distributions

Tracy-Widom and Baik-Rains distributions appear as universal limit distributions for height fluctuations in the one-dimensional Kardar-Parisi-Zhang (KPZ) \textit{stochastic} partial differential equation (PDE). We obtain the same universal distributions in the spatiotemporally chaotic, nonequilibrium, but statistically steady state (NESS) of the one-dimensional Kuramoto-Sivashinsky (KS) \textit{deterministic} PDE, by carrying out extensive pseudospectral direct numerical simulations to obtain the spatiotemporal evolution of the KS height profile $h(x,t)$ for different initial conditions. We establish, therefore, that the statistical properties of the 1D KS PDE in this state are in the 1D KPZ universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Comparisons of wave dynamics in Hodgkin-Huxley and Markov-state formalisms for the Sodium (Na) channel in some mathematical models for human cardiac tissue

We compare and contrast spiral- and scroll-wave dynamics in five different mathematical models for cardiac tissue. The first is the TP06 model, due to ten Tusscher and Panfilov, which is based on the Hodgkin-Huxley formalism; the remaining four are Markov-state models, MM1 WT and MM2 WT, for the wild-type (WT) Na channel, and MM1 MUT and MM2 MUT, for the mutant Na channel. Our results are based on extensive direct numerical simulations of waves of electrical activation in these models, in two- and three-dimensional (2D and 3D) homogeneous simulation domains and also in domains with localised heterogeneities, either obstacles with randomly distributed inexcitable regions or mutant cells in a wild-type background. Our study brings out the sensitive dependence of spiral- and scroll-wave dynamics on these five models and the parameters that define them. We also explore the control of spiral-wave turbulence in these models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A regularity criterion for solutions of the three-dimensional Cahn-Hilliard-Navier-Stokes equations and associated computations

We consider the 3D Cahn-Hilliard equations coupled to, and driven by, the forced, incompressible 3D Navier-Stokes equations. The combination, known as the Cahn-Hilliard-Navier-Stokes (CHNS) equations, is used in statistical mechanics to model the motion of a binary fluid. The potential development of singularities (blow-up) in the contours of the order parameter $ϕ$ is an open problem. To address this we have proved a theorem that closely mimics the Beale-Kato-Majda theorem for the $3D$ incompressible Euler equations [Beale et al. Commun. Math. Phys., Commun. Math. Phys., ${\rm 94}$, $ 61-66 ({\rm 1984})$]. By taking an $L^{\infty}$ norm of the energy of the full binary system, designated as $E_{\infty}$, we have shown that $\int_{0}^{t}E_{\infty}(τ)\,dτ$ governs the regularity of solutions of the full 3D system. Our direct numerical simulations (DNSs), of the 3D CHNS equations, for (a) a gravity-driven Rayleigh Taylor instability and (b) a constant-energy-injection forcing, with $128^3$ to $512^3$ collocation points and over the duration of our DNSs, confirm that $E_{\infty}$ remains bounded as far as our computations allow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Binary-Fluid Turbulence: Signatures of Multifractal Droplet Dynamics and Dissipation Reduction

We present an extensive direct numerical simulation of statistically steady, homogeneous, isotropic turbulence in two-dimensional, binary-fluid mixtures with air-drag-induced friction by using the Cahn-Hilliard-Navier-Stokes equations. We choose parameters, e.g., the surface tension, such that we have a droplet of the minority phase moving inside a turbulent background of the majority phase. We characterize the deformation of the droplet and show that it displays multifractal dynamics. The probability distribution functions of the components of the acceleration of the center of mass of the droplet exhibit wide, non-Gaussian tails. Our study reveals that the droplet enhances the energy spectrum $E(k)$ when the wavenumber $k$ is large; this enhancement leads to dissipation reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamic multiscaling in magnetohydrodynamic turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Instability of spiral and scroll waves in the presence of a gradient in the fibroblast density: the effects of fibroblast-myocyte coupling

Fibroblast-myocyte coupling can modulate electrical-wave dynamics in cardiac tissue. In diseased hearts, the distribution of fibroblasts is heterogeneous, so there can be gradients in the fibroblast density (henceforth we call this GFD) especially from highly injured regions, like infarcted or ischemic zones, to less-wounded regions of the tissue. Fibrotic hearts are known to be prone to arrhythmias, so it is important to understand the effects of GFD in the formation and sustenance of arrhythmic re- entrant waves, like spiral or scroll waves. Therefore, we investigate the effects of GFD on the stability of spiral and scroll waves of electrical activation in a state-of-the- art mathematical model for cardiac tissue in which we also include fibroblasts. By introducing GFD in controlled ways, we show that spiral and scroll waves can be unstable in the presence of GFDs because of regions with varying spiral or scroll-wave frequency ω, induced by the GFD. We examine the effects of the resting membrane potential of the fibroblast and the number of fibroblasts attached to the myocytes on the stability of these waves. Finally, we show that the presence of GFDs can lead to the formation of spiral waves at high-frequency pacing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sticking Transition in a Minimal Model for the Collisions of Active Particles in Quantum Fluids

Particles of low velocity, travelling without dissipation in a superfluid, can interact and emit sound when they collide. We propose a minimal model in which the equations of motion of the particles, including a short-range repulsive force, are self-consistently coupled with the Gross-Pitaevskii equation. We use this model to demonstrate the existence of an effective superfluid-mediated attractive interaction between the particles; and we study numerically the collisional dynamics of particles as a function of their incident kinetic energy and the length-scale of the repulsive force. We find a transition from almost elastic to completely inelastic (sticking) collisions as the parameters are tuned. We find that aggregation and clustering result from this sticking transition in multi-particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of fibroblasts on wave dynamics in a mathematical model for human ventricular tissue

We present systematic numerical studies of electrical-wave propagation in two-dimensional (2D) and three-dimensional (3D) mathematical models, for human, ventricular tissue with myocyte cells that are attached (a) regularly and (b) randomly to distributed fibroblasts. In both these cases we show that there is a parameter regime in which single rotating spiral- and scroll-wave states (RS) retain their integrity and do not evolve to a state ST that displays spatiotemporal chaos and turbulence. However, in another range of parameters, we observe a transition from ST to RS states in both 2D or 3D domains and for both cases (a) and (b). Our studies show that the ST-RS transition and rotation period of a spiral or scroll wave in the RS state depends on (i) the coupling strength between myocytes and fibroblasts and (ii) the number of fibroblasts attached to myocytes. We conclude that myocyte-fibroblast coupling strength and the number of fibroblasts are more important for the ST-RS transition than the precise way in which fibroblasts are distributed over myocyte tissue.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two-dimensional, homogeneous, isotropic fluid turbulence with polymer additives

We present the most extensive direct numerical simulations, attempted so far, of statistically steady, homogeneous, isotropic turbulence in two-dimensional fluid films with air-drag-induced friction and with polymer additives. Our study reveals that the polymers (a) reduce the total fluid energy, enstrophy, and palinstrophy, (b) modify the fluid energy spectrum both in inverse- and forward-cascade regimes, (c) reduce small-scale intermittency, (d) suppress regions of large vorticity and strain rate, and (e) stretch in strain-dominated regions. We compare our results with earlier experimental studies; and we propose new experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscaling in superfluid turbulence: A shell-model study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Elliptical Tracers in Two-dimensional, Homogeneous, Isotropic Fluid Turbulence: the Statistics of Alignment, Rotation, and Nematic Order

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Homogeneous Isotropic Superfluid Turbulence in Two Dimensions: Inverse and Forward Cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phases, transitions, and patterns in the one-dimensional Extended Bose-Hubbard model

We carry out an extensive study of the phase diagram of the extended Bose Hubbard model, with a mean filling of one boson per site, in one dimension by using the density matrix renormalization group and show that it can have Superfluid (SF), Mott-insulator (MI), density-wave (DW) and Haldane-insulator (HI) phases depending on the precise value of filling and how edge states are handled. We show that the critical exponents and central charge for the HI-DW, MI-HI and SF-MI transitions are consistent with those for models in the two-dimensional Ising, Gaussian, and Berezinskii-Kosterlitz-Thouless (BKT) universality classes, respectively; and we suggest that the SF-HI transition may be more exotic than a simple BKT transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Random-Phase-Approximation Excitation Spectra for Bose-Hubbard Models

We obtain the excitation spectra of the following three generalized Bose-Hubbard (BH) models: (1) a two-species generalization of the spinless BH model, (2) a single-species, spin-1 BH model, and (3) the extended Bose-Hubbard model (EBH) for spinless interacting bosons of one species. In all the phases of these models we provide a unified treatment of random-phase-approximation (RPA) excitation spectra. These spectra have gaps in all the MI phases and gaps in the DW phases in the EBH model; they are gapless in all the SF phases in these models and in the SS phases in the EBH model. We obtain the dependence of (a) gaps $Δ$ and (b) the sound velocity $u_s$ on the parameters of these models and examine $Δ$ and $u_s$ as these systems go through phase transitions. At the SF-MI transitions in the spin-1 BH model, $u_s$ goes to zero continuously (discontinuously) for MI phases with an odd (even) number of bosons per site; this is consistent with the natures of these transition in mean-field theory. In the SF phases of these models, our excitation spectra agree qualitatively, at weak couplings, with those that can be obtained from Gross-Pitaevskii-type models. We compare the results of our work with earlier studies of related models and discuss implications for experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Regimes of nonlinear depletion and regularity in the 3D Navier-Stokes equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical Properties of the Intrinsic Geometry of Heavy-particle Trajectories in Two-dimensional, Homogeneous, Isotropic Turbulence

We obtain, by extensive direct numerical simulations, trajectories of heavy inertial particles in two-dimensional, statistically steady, homogeneous, and isotropic turbulent flows, with friction. We show that the probability distribution function $\mathcal{P}(κ)$, of the trajectory curvature $κ$, is such that, as $κ\to \infty$, $\mathcal{P}(κ) \sim κ^{-h_{\rm r}}$, with $h_{\rm r} = 2.07 \pm 0.09$. The exponent $h_{\rm r}$ is universal, insofar as it is independent of the Stokes number ($\rm{St}$) and the energy-injection wave number. We show that this exponent lies within error bars of their counterparts for trajectories of Lagrangian tracers. We demonstrate that the complexity of heavy-particle trajectories can be characterized by the number $N_{\rm I}(t,{\rm St})$ of inflection points (up until time $t$) in the trajectory and $n_{\rm I} ({\rm St}) \equiv \lim_{t\to\infty} \frac{N_{\rm I} (t,{\rm St})}{t} \sim {\rm St}^{-Δ}$, where the exponent $Δ= 0.33 \pm0.02$ is also universal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistics of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superfluid Mutual-friction Coefficients from Vortex Dynamics in the Two-dimensional Galerkin-truncated Gross-Pitaevskii Equation

We present algorithms for the ab-initio determination of the temperature ($T$) dependence of the mutual-friction coefficients $α$ and $α'$ and the normal-fluid density $ρ_{\rm n}$ in the two-dimensional (2D) Galerkin-truncated Gross-Pitaevskii system. Our algorithms enable us to determine $α(T)$, even though fluctuations in 2D are considerably larger than they are in 3D. We also examine the implications of our measurements of $α'(T)$ for the Iordanskii force, whose existence is often questioned.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universal Statistical Properties of Inertial-particle Trajectories in Three-dimensional, Homogeneous, Isotropic, Fluid Turbulence

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(ϕ)$, of the angle $ϕ$ 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(ϕ) \sim ϕ^{-γ}$, with a new universal exponent $γ\simeq 4$. Furthermore, the PDFs of the trajectory curvature $κ$ and modulus $θ$ of the torsion $\vartheta$ have power-law tails that scale, respectively, as $P(κ) \sim κ^{-h_κ}$, as $κ\to \infty$, and $P(θ) \sim θ^{-h_θ}$, as $θ\to \infty$, with exponents $h_κ\simeq 2.5$ and $h_θ\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 $γ$, $h_κ$ and $h_θ$ 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}^{-Δ}$, with $Δ\simeq 0.4$ a universal exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multiscaling in Hall-Magnetohydrodynamic Turbulence: Insights from a Shell Model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bose-Hubbard Models in Confining Potentials: An Inhomogeneous Mean-Field Theory

We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we develop an inhomogeneous mean-field theory. Our results for the BH model with one type of spinless bosons agrees quantitatively with quantum Monte Carlo (QMC) simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculation on experimentally realistic, large 3D systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also generalize our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons we obtain rich phase diagrams with a variety of SF and MI phases and associated shells, when we include a quadratic confining potential. For the spin-1 BH model we show, in a representative case, that the system can display alternating shells of polar SF and MI phases; and we make interesting predictions for experiments in such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamic Multiscaling in Two-dimensional Fluid Turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phases, transitions, and boundary conditions in a model of interacting bosons

We carry out an extensive study of the phase diagrams of the extended Bose Hubbard model, with a mean filling of one boson per site, in one dimension by using the density matrix renormalization group and show that it contains Superfluid (SF), Mott-insulator (MI), density-wave (DW) and Haldane-insulator (HI) phases. We show that the critical exponents and central charges for the HI-DW, MI-HI and SF-MI transitions are consistent with those for models in the two-dimensional Ising, Gaussian, and Berezinskii-Kosterlitz-Thouless (BKT) universality classes, respectively; and we suggest that the SF-HI transition may be more exotic than a simple BKT transition. We show explicitly that different boundary conditions lead to different phase diagrams.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Real-space Manifestations of Bottlenecks in Turbulence Spectra

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Inhomogeneous Extended Bose-Hubbard Model: A Mean-Field Theory

We develop an inhomogeneous mean-field theory for the extended Bose-Hubbard model with a quadratic, confining potential. In the absence of this potential, our mean-field theory yields the phase diagram of the homogeneous extended Bose-Hubbard model. This phase diagram shows a superfluid (SF) phase and lobes of Mott-insulator (MI), density-wave (DW), and supersolid (SS) phases in the plane of the chemical potential (μ) and on-site repulsion (U); we present phase diagrams for representative values of V, the repulsive energy for bosons on nearest-neighbor sites. We demonstrate that, when the confining potential is present, superfluid and density-wave order parameters are nonuniform; in particular, we obtain, for a few representative values of parameters, spherical shells of SF, MI, DW, and SS phases. We explore the implications of our study for experiments on cold-atom dipolar condensates in optical lattices in a confining potential

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Direct numerical simulations of statistically steady, homogeneous, isotropic fluid turbulence with polymer additives

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Persistence Problem in Two-Dimensional Fluid Turbulence

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter $Λ$ 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 $θ=2.9\pm0.2$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Systematics of the magnetic-Prandtl-number dependence of homogeneous, isotropic magnetohydrodynamic turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamo Onset as a First-Order Transition: Lessons from a Shell Model for Magnetohydrodynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Manifestations of Drag Reduction by Polymer Additives in Decaying, Homogeneous, Isotropic Turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Turbulence and Multiscaling in the Randomly Forced Navier Stokes Equation

We present an extensive pseudospectral study of the randomly forced Navier-Stokes equation (RFNSE) stirred by a stochastic force with zero mean and a variance $\sim k^{4-d-y}$, where $k$ is the wavevector and the dimension $d = 3$. We present the first evidence for multiscaling of velocity structure functions in this model for $y \geq 4$. We extract the multiscaling exponent ratios $ζ_p/ζ_2$ by using extended self similarity (ESS), examine their dependence on $y$, and show that, if $y = 4$, they are in agreement with those obtained for the deterministically forced Navier-Stokes equation ($3d$NSE). We also show that well-defined vortex filaments, which appear clearly in studies of the $3d$NSE, are absent in the RFNSE.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The Exotic Barium Bismuthates

We review the remarkable properties, including superconductivity, charge-density-wave ordering, and metal-insulator transitions, of lead- and potassium-doped barium bismuthate. We discuss some of the early theoretical studies of these systems. Our recent theoretical work, on the negative-$U\/$, extended-Hubbard model for these systems, is also described. Both the large- and intermediate-$U\/$ regimes of this model are examined, using mean-field and random-phase approximations, particularly with a view to fitting various experimental properties of these bismuthates. On the basis of our studies, we point out possibilities for exotic physics in these systems. We also emphasize the different consequences of electronic and phonon-mediated mechanisms for the negative $U.\/$ We show that, for an electronic mechanism, the \secin \,\,phases of these bismuthates must be unique, with their transport properties {\it dominated by charge $\pm 2e$ Cooperon bound states}. This can explain the observed difference between the optical and transport gaps. We propose other experimental tests for this novel mechanism of charge transport and comment on the effects of disorder.

0 citations0 reviewsNo code linkedOpen access