Source author record

Prasad Perlekar

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

29 published item(s)

preprint2026arXiv

Bridging Filtering and Point-Splitting Approaches for Variable-Density Flows

Energy transfer in turbulent flows is typically described either through correlation functions, via the Kármán-Howarth-Monin relation, or through a scale-by-scale budget of filtered energy (Frisch 1995). For constant-density turbulence, the equivalence between these two descriptions is well understood. In compressible turbulence, however, several definitions of filtered energy exist, and for most of them the associated formulation in terms of correlation functions is unclear. We develop a general framework to determine the multi-point correlation functions corresponding to any specified filtered energy. As a demonstration, we show that the Favre filtered energy--defined as the ratio of the squared filtered momentum to the filtered density--and the terms in its budget can be written as an infinite series of multi-point correlation functions. We validate the framework numerically using three-dimensional buoyancy-driven bubbly flows.

preprint2021arXiv

Phase ordering, topological defects, and turbulence in the 3D incompressible Toner-Tu equation

We investigate phase ordering dynamics of the incompressible Toner-Tu equation in three dimensions. We show that the phase ordering proceeds via defect merger events and the dynamics is controlled by the Reynolds number Re. At low Re, the dynamics is similar to that of the Ginzburg-Landau equation. At high Re, turbulence controls phase ordering. In particular, we observe a forward energy cascade from the coarsening length scale to the dissipation scale.

preprint2021arXiv

Rate of formation of caustics in heavy particles advected by turbulence

The rate of collision and the relative velocities of the colliding particles in turbulent flows is a crucial part of several natural phenomena, e.g., rain formation in warm clouds and planetesimal formation in a protoplanetary disks. The particles are often modeled as passive, but heavy and inertial. Within this model, large relative velocities emerge due to formation of singularities (caustics) of in the gradient matrix of the velocities of the particles. Using extensive direct numerical simulations of heavy particles in both two (direct and inverse cascade) and three dimensional turbulent flows we calculate the rate of formation of caustics, $J$ as a function of the Stokes number (${\rm St}$).The best approximation to our data is $J \sim \exp(-C/{\rm St})$, in the limit ${\rm St} \to 0 $ where $C$ is a non-universal constant.

preprint2020arXiv

Fixation in Competing Populations: Diffusion and Strategies for Survival

How should dispersal strategies be chosen to increase the likelihood of survival of a species? We obtain the answer for the spatially extended versions of three well-known models of two competing species with unequal diffusivities. Though identical at the mean-field level, the three models exhibit drastically different behaviour leading to different optimal strategies for survival, with or without a selective advantage for one species. With conserved total particle number, dispersal has no effect on survival probability. With a fluctuating number, faster dispersal is advantageous if intra-species competition is present, while moving slower is the optimal strategy for the disadvantaged species if there is no intra-species competition: it is imperative to include fluctuations to properly formulate survival strategies.

preprint2019arXiv

Liquid velocity fluctuations and energy spectra in three-dimensional buoyancy driven bubbly flows

We present a direct numerical simulation (DNS) study of pseudo-turbulence in buoyancy driven bubbly flows for a range of Reynolds ($\Rey$) and Atwood ($\At$) numbers. We study the probability distribution function of the horizontal and vertical liquid velocity fluctuations and find them to be in quantitative agreement with the experiments. The energy spectrum shows the $k^{-3}$ scaling at high $\Rey$ and becomes steeper on reducing the $\Rey$. To investigate the spectral transfers in the flow, we derive the scale-by-scale energy budget equation. Our analysis shows that, for scales smaller than the bubble diameter, the net production because of the surface tension and the kinetic energy flux balances viscous dissipation to give the $k^{-3}$ scaling of the energy spectrum for both low and high $\At$.

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.

preprint2016arXiv

Coherent structures and extreme events in rotating multiphase turbulent flows

By using direct numerical simulations (DNS) at unprecedented resolution we study turbulence under rotation in the presence of simultaneous direct and inverse cascades. The accumulation of energy at large scale leads to the formation of vertical coherent regions with high vorticity oriented along the rotation axis. By seeding the flow with millions of inertial particles, we quantify -for the first time- the effects of those coherent vertical structures on the preferential concentration of light and heavy particles. Furthermore, we quantitatively show that extreme fluctuations, leading to deviations from a normal-distributed statistics, result from the entangled interaction of the vertical structures with the turbulent background. Finally, we present the first-ever measurement of the relative importance between Stokes drag, Coriolis force and centripetal forces along the trajectories of inertial particles. We discover that vortical coherent structures lead to unexpected diffusion properties for heavy and light particles in the directions parallel and perpendicular to the rotation axis.

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.

preprint2015arXiv

Clustering of vertically constrained passive particles in homogeneous, isotropic turbulence

We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale, that maximizes the clustering of the particles.

preprint2015arXiv

On the Peterlin approximation for turbulent flows of polymer solutions

We study the impact of the Peterlin approximation on the statistics of the end-to-end separation of poly- mers in a turbulent flow. The FENE and FENE-P models are numerically integrated along a large number of Lagrangian trajectories resulting from a direct numerical simulation of three-dimensional homogeneous isotropic turbulence. Although the FENE-P model yields results in qualitative agreement with those of the FENE model, quantitative differences emerge. The steady-state probability of large extensions is overesti- mated by the FENE-P model. The alignment of polymers with the eigenvectors of the rate-of-strain tensor and with the direction of vorticity is weaker when the Peterlin approximation is used. At large Weissenberg numbers, both the correlation times of the extension and of the orientation of polymers are underestimated by the FENE-P model.

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.

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.

preprint2013arXiv

Attracting fixed points for heavy particles in the vicinity of a vortex pair

We study the behaviour of heavy inertial particles in the flow field of two like-signed vortices. In a frame co-rotating with the two vortices, we find that stable fixed points exist for these heavy inertial particles; these stable frame-fixed points exist only for particle Stokes number $St<St_{cr}$. We estimate $St_{cr}$ and compare this with direct numerical simulations, and find that the addition of viscosity increases the $St_{cr}$ slightly. We also find that the fixed points become more stable with increasing $St$ until they abruptly disappear at $St=St_{cr}$. These frame-fixed points are between fixed points and limit cycles in character.

preprint2013arXiv

Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [ Phys. Rev. E 47 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed.

preprint2013arXiv

Spinodal decomposition in homogeneous and isotropic turbulence

We study the competition between domain coarsening in a symmetric binary mixtures below the critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in presence of turbulence. The physics of the process shares remarkable similarities with the behaviour of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although in the absence of flow the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

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.

preprint2012arXiv

Population dynamics in compressible flows

Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another

preprint2012arXiv

Statistics of polymer extensions in turbulent channel flow

We present direct numerical simulations of turbulent channel flow with passive Lagrangian polymers. To understand the polymer behavior we investigate the behavior of infinitesimal line elements and calculate the probability distribution function (PDF) of finite-time Lyapunov exponents and from them the corresponding Cramer's function for the channel flow. We study the statistics of polymer elongation for both the Oldroyd-B model (for Weissenberg number $\Wi <1$) and the FENE model. We use the location of the minima of the Cramer's function to define the Weissenberg number precisely such that we observe coil-stretch transition at $\Wi\approx1$. We find agreement with earlier analytical predictions for PDF of polymer extensions made by Balkovsky, Fouxon and Lebedev [Phys. Rev. Lett., 84, 4765 (2000).] for linear polymers (Oldroyd-B model) with $\Wi<1$ and by Chertkov [Phys. Rev. Lett., 84, 4761 (2000).] for nonlinear FENE-P model of polymers. For $\Wi>1$ (FENE model) the polymer are significantly more stretched near the wall than at the center of the channel where the flow is closer to homogenous isotropic turbulence. Furthermore near the wall the polymers show a strong tendency to orient along the stream-wise direction of the flow but near the centerline the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows.

preprint2012arXiv

Unified framework for a side-by-side comparison of different multicomponent algorithms: lattice Boltzmann vs. phase field model

Lattice Boltzmann Models (LBM) and Phase Field Models (PFM) are two of the most widespread approaches for the numerical study of multicomponent fluid systems. Both methods have been successfully employed by several authors but, despite their popularity, still remains unclear how to properly compare them and how they perform on the same problem. Here we present a unified framework for the direct (one-to-one) comparison of the multicomponent LBM against the PFM. We provide analytical guidelines on how to compare the Shan-Chen (SC) lattice Boltzmann model for non-ideal multicomponent fluids with a corresponding free energy (FE) lattice Boltzmann model. Then, in order to properly compare the LBM vs. the PFM, we propose a new formulation for the free energy of the Cahn-Hilliard/Navier-Stokes equations. Finally, the LBM model is numerically compared with the corresponding phase field model solved by means of a pseudo-spectral algorithm. This work constitute a first attempt to set the basis for a quantitative comparison between different algorithms for multicomponent fluids. We limit our scope to the few of the most common variants of the two most widespread methodologies, namely the lattice Boltzmann model (SC and FE variants) and the Phase Field Model.

preprint2011arXiv

Cumulative compressibility effects on population dynamics in turbulent flows

Bacteria and plankton populations living in oceans and lakes reproduce and die under the in- fluence of turbulent currents. Turbulent transport interacts in a complex way with the dynamics of populations because the typical reproduction time of microorganism is within the inertial range of turbulent time scales. In the present manuscript we quantitatively investigate the effect of flow compressibility on the dynamics of populations. While a small compressibility can be induced by several physical mechanisms, like density mismatch or the finite size of microorganisms with respect to the fluid turbulence, its effect on the the carrying capacity of the ecosystem can be dramatic. We report, for the first time, how a small compressibility can produce a sizeable reduction in the carrying capacity, due to an integrated effect made possible by the long replication times of the organisms with respect to turbulent time scales. A full statistical quantification of the fluctuations of population concentration field leads to data collapse over a broad range in parameter space.

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.

preprint2011arXiv

Droplet size distribution in homogeneous isotropic turbulence

We study the physics of droplet breakup in a statistically stationary homogeneous and isotropic turbulent flow by means of high resolution numerical investigations based on the multicomponent lattice Boltzmann method. We verified the validity of the criterion proposed by Hinze (1955) for droplet breakup and we measured the full probability distribution function (pdf) of droplets radii at different Reynolds numbers and for different volume fraction. By means of a Lagrangian tracking we could follow individual droplets along their trajectories, define a local Weber number based on the velocity gradients and study its cross-correlation with droplet deformation.

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$.

preprint2010arXiv

Droplet breakup in homogeneous and isotropic turbulence

This fluid dynamics video shows the breakup of a droplet in a stationary homogeneous and isotropic turbulent flow. We consider droplets with the same density of the transporting fluid. The droplets and the fluid are numerically modelled by means of a multicompo- nent Lattice-Boltzmann method. The turbulent fluid is maintained through a large scale stirring force and the radius of stable droplets, for the parameters in our simulation, is larger than the Kolmogorov scale. Events of droplet deformation, break-up and aggregation are clearly visible from the movie. With the present database droplet evo- lution can be studied from both an Eulerian and Lagrangian point of view. The Kolmogorov-Hinze criteria for droplets break-up can be tested also by means of simulations with different viscosity contrast between the two components.

preprint2010arXiv

Population dynamics at high Reynolds number

We study the statistical properties of population dynamics evolving in a realistic two-dimensional compressible turbulent velocity field. We show that the interplay between turbulent dynamics and population growth and saturation leads to quasi-localization and a remarkable reduction in the carrying capacity. The statistical properties of the population density are investigated and quantified via multifractal scaling analysis. We also investigate numerically the singular limit of negligibly small growth rates and delocalization of population ridges triggered by uniform advection.

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.

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.

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.