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

Mahendra K. Verma

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

ResearcherUnclaimed source record
physics.flu-dynphysics.plasm-phnlin.CDastro-ph.SRphysics.comp-phastro-ph.EPcond-mat.stat-mechnlin.PSastro-phcond-mathep-thnlin.SIphysics.geo-phphysics.hist-phphysics.space-ph

Catalog footprint

What is connected

41works
15topics
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

41 published item(s)

preprint2022arXiv

Equilibrium states of Burgers and KdV equations

We simulate KdV and dissipation-less Burgers equations using delta-correlated random noise as initial condition. We observe that the energy fluxes of the two equations remain zero throughout, thus indicating their equilibrium nature. We characterize the equilibrium states using Gaussian probability distribution for the real space field, and using Boltzmann distribution for the modal energy. We show that the single soliton of the KdV equation too exhibits zero energy flux, hence it is in equilibrium. We argue that the energy flux is a good measure for ascertaining whether a system is in equilibrium or not.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scalable Multi-node Fast Fourier Transform on GPUs

In this paper, we present the details of our multi-node GPU-FFT library, as well its scaling on Selene HPC system. Our library employs slab decomposition for data division and MPI for communication among GPUs. We performed GPU-FFT on $1024^3$, $2048^3$, and $4096^3$ grids using a maximum of 512 A100 GPUs. We observed good scaling for $4096^3$ grid with 64 to 512 GPUs. We report that the timings of multicore FFT of $1536^3$ grid with 196608 cores of Cray XC40 is comparable to that of GPU-FFT of $2048^3$ grid with 128 GPUs. The efficiency of GPU-FFT is due to the fast computation capabilities of A100 card and efficient communication via NVlink.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent Drag Reduction in Magnetohydrodynamic Turbulence and Dynamo from Energy Flux Perspectives

In this review, we describe turbulent drag reduction in a variety of flows using a universal framework of energy flux. In a turbulent flow with dilute polymers and magnetic field, the kinetic energy injected at large scales cascades to the velocity field at intermediate scales, as well as to the polymers and magnetic field at all scales. Consequently, the kinetic energy flux, $ Π_u(k) $, is suppressed in comparison to the pure hydrodynamic turbulence. We argue that the suppression of $Π_u(k)$ is an important factor in the reduction of the inertial force $\langle {\bf u \cdot \nabla u} \rangle$ and \textit{turbulent drag}. This feature of turbulent drag reduction is observed in polymeric, magnetohydrodynamic, quasi-static magnetohydrodynamic, and stably-stratified turbulence, and in dynamos. In addition, it is shown that turbulent drag reduction in thermal convection is due to the smooth thermal plates, similar to the turbulent drag reduction over bluff bodies. In all these flows, turbulent drag reduction often leads to a strong large-scale velocity in the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Euler Turbulence and thermodynamic equilibrium

We perform a unique direct numerical simulation of Euler turbulence using delta-correlated velocity field as an initial condition, and report a full range of $k^2$ and $k$ energy spectra for 3D and 2D flows respectively, zero energy flux, and Maxwell-Boltzmann distribution for the velocity field. These are direct verification of the predictions of the absolute equilibrium theory of turbulence. For a coherent vortex as an initial condition, Euler turbulence transitions from a mixture of nonequilibrium-equilibrium state to a equilibrium state through a process called thermalization. In this letter, we present a model for thermalization in Euler turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kolmogorov flow: Linear Stability and Energy Transfers in a minimal low-dimensional model

In this paper, we derive a four-mode model for the Kolmogorov flow by employing Galerkin truncation and Craya-Herring basis for the decomposition of velocity field. After this, we perform a bifurcation analysis of the model. Though our low-dimensional model has fewer modes than the past models, it captures the essential features of the primary bifurcation of the Kolmogorov flow. For example, it reproduces the critical Reynolds number for the supercritical pitchfork bifurcation and the flow structures of the past works. We also demonstrate energy transfers from intermediate scales to large scales. We perform direct numerical simulations of the Kolmogorov flow and show that our model predictions match with the numerical simulations very well.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revisiting Reynolds and Nusselt numbers in turbulent thermal convection

In this paper, we extend Grossmann and Lohse's (GL) model [Phys. Rev. Lett. {\bf 86}, 3316 (2001)] for the predictions of Reynolds number (Re) and Nusselt number (Nu) in turbulent Rayleigh-Bénard convection (RBC). Towards this objective, we use functional forms for the prefactors of the dissipation rates in the bulk and the boundary layers. The functional forms arise due to inhibition of nonlinear interactions in the presence of walls and buoyancy compared to free turbulence, along with a deviation of viscous boundary layer profile from Prandtl-Blasius theory. We perform 60 numerical runs on a three-dimensional unit box for a range of Rayleigh numbers (Ra) and Prandtl numbers (Pr) and determine the aforementioned functional forms using machine learning. The revised predictions are in better agreement with the past numerical and experimental results than those of the GL model, especially for extreme Prandtl numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent drag reduction in magnetohydrodynamic and quasi-static magnetohydrodynamic turbulence

In hydrodynamic turbulence, the kinetic energy injected at large scales cascades to the inertial range, leading to a constant kinetic energy flux. In contrast, in magnetohydrodynamic (MHD) turbulence, a fraction of kinetic energy is transferred to the magnetic energy. Consequently, for the same kinetic energy injection rate, the kinetic energy flux in MHD turbulence is reduced compared to its hydrodynamic counterpart. This leads to relative weakening of the nonlinear term ($\la | ({\bf u}\cdot \nabla) {\bf u} | \ra$, where ${\bf u}$ is the velocity field) and turbulent drag, but strengthening of the velocity field in MHD turbulence. We verify the above using shell model simulations of hydrodynamic and MHD turbulence. Quasi-static MHD turbulence too exhibits turbulent drag reduction similar to MHD turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Transient flows and reorientations of large-scale convection in a cubic cell

The transient processes of a turbulent large-scale convective circulation (LSC) in a cubic cell are investigated using large-eddy simulations for Rayleigh number $\Ray=10^8$ and Prandtl number $\Pran=0.7$. For the first time, we have explicitly shown that LSC is accompanied by large-scale azimuthal flows with non-zero total angular momentum. It is also shown that solid-body rotation of the entire fluid is not realized. It is found that correlation between rotation of LSC plane and the mean azimuthal motion is high during quasiperiodic oscillations of LSC near the diagonal plane and relatively weak during LSC reorientations. We propose a new plausible scenario for the reorientations of the LSC in a cube that does not involve a mean azimuthal flow. Instead of a single-roll, we introduce the superposition of a pair of large-scale orthogonal quasi-two-dimensional (Q2D) rolls and the reorientation of the LSC occurs as a result of the cessation of one of the Q2D rolls. This scenario is consistent with all known experimental and numerical data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Turbulence in a stably stratified fluid: Onset of global anisotropy as a function of the Richardson number

It is necessary to introduce an external forcing to induce turbulence in a stably stratified fluid. The Heisenberg eddy viscosity technique should in this case suffice to calculate a space-time averaged quantity like the global anisotropy parameter as a function of the Richardson number. We find analytically that the anisotropy increases linearly with the Richardson number, with a small quadratic correction. A numerical simulation of the complete equations shows the linear behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of large-scale quantities in Rayleigh-Bénard convection

In this paper we estimate the relative strengths of various terms of the Rayleigh-Bénard equations. Based on these estimates and scaling analysis, we derive a general formula for the large-scale velocity, $U$, or the Péclet number that is applicable for arbitrary Rayleigh number $\mathrm{Ra}$ and Prandtl number $\mathrm{Pr}$. Our formula fits reasonably well with the earlier simulation and experimental results. Our analysis also shows that the wall-bounded convection has enhanced viscous force compared to free turbulence. We also demonstrate how correlations deviate the Nusselt number scaling from the theoretical prediction of $\mathrm{Ra}^{1/2}$ to the experimentally observed scaling of nearly $\mathrm{Ra}^{0.3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamos at extreme magnetic Prandtl numbers: Insights from shell models

We present an MHD shell model suitable for computation of various energy fluxes of magnetohydrodynamic turbulence for very small and very large magnetic Prandtl numbers $\mathrm{Pm}$; such computations are inaccessible to direct numerical simulations. For small $\mathrm{Pm}$, we observe that both kinetic and magnetic energy spectra scale as $k^{-5/3}$ in the inertial range, but the dissipative magnetic energy scales as $k^{-11/3}\exp(-k/k_η)$. Here, the kinetic energy at large length scale feeds the large-scale magnetic field that cascades to small-scale magnetic field, which gets dissipated by Joule heating. The large-$\mathrm{Pm}$ dynamo has a similar behaviour except that the dissipative kinetic energy scales as $k^{-13/3}$. For this case, the large-scale velocity field transfers energy to the large-scale magnetic field, which gets transferred to small-scale velocity and magnetic fields; the energy of the small-scale magnetic field also gets transferred to the small-scale velocity field, and the energy thus accumulated is dissipated by the viscous force.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scaling of large-scale quantities in Rayleigh-Bénard convection

We derive a formula for the Péclet number ($\mathrm{Pe}$) by estimating the relative strengths of various terms of the momentum equation. Using direct numerical simulations in three dimensions we show that in the turbulent regime, the fluid acceleration is dominated by the pressure gradient, with relatively small contributions arising from the buoyancy and the viscous term, in the viscous regime, acceleration is very small due to a balance between the buoyancy and the viscous term. Our formula for $\mathrm{Pe}$ describes the past experiments and numerical data quite well. We also show that the ratio of the nonlinear term and the viscous term is $\mathrm{Re} \mathrm{Ra}^{-0.14}$, where $\mathrm{Re}$ and $\mathrm{Ra}$ are Reynolds and Rayleigh numbers respectively, and that the viscous dissipation rate $ε_u = (U^3/d) \mathrm{Ra}^{-0.21}$, where $U$ is the root mean square velocity and $d$ is the distance between the two horizontal plates. The aforementioned decrease in nonlinearity compared to free turbulence arises due to the wall effects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Shell Model for Buoyancy-Driven Turbulence

In this paper we present a unified shell model for stably stratified and convective turbulence. Numerical simulation of this model for stably stratified flow shows Bolgiano-Obukhbov scaling in which the kinetic energy spectrum varies as $k^{-11/5}$. We also observe a dual scaling ($k^{-11/5}$ and $k^{-5/3}$) for a limited range of parameters. The shell model of convective turbulence yields Kolmogorov's spectrum. These results are consistent with the energy flux and energy feed due to buoyancy, and are in good agreement with direct numerical simulations of Kumar {\em et al.} [Phys. Rev. E {\bf 90}, 023016 (2014)].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy transfers in dynamos with small magnetic Prandtl numbers

We perform numerical simulation of dynamo with magnetic Prandtl number $\mathrm{Pm} =0.2$ on $1024^3$ grid, and compute the energy fluxes and the shell-to-shell energy transfers. These computations indicate that the magnetic energy growth takes place mainly due to the energy transfers from large-scale velocity field to large-scale magnetic field and that the magnetic energy flux is forward. The steady-state magnetic energy is much smaller than the kinetic energy, rather than equipartition; this is because the magnetic Reynolds number is near the dynamo transition regime. We also contrast our results with those for dynamo with $\mathrm{Pm} =20$ and decaying dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Flow reversals in turbulent convection with free-slip walls

We perform numerical simulations of turbulent convection for infinite Prandtl number with free-slip walls, and study the dynamics of flow reversals. We show interesting correlations between the flow reversals and the nonlinear interactions among the large-scale flow structures represented by the modes $(1,1), (2,1), (3,1)$ and some others. After a flow reversal, the odd modes, e.g. $(1,1), (3,1)$, switch sign, but the even modes, e.g. $(2,2)$, retain their sign. The mixed modes $(1,2)$ and $(2,1)$ fluctuate around zero. Using the properties of the modes and their interactions, we show that they form a Klein four-group $Z_2 \times Z_2$. We also show that for the free-slip boundary condition, the corner rolls and vortex reconnection are absent during a flow reversal, in contrast to active role played by them in flow reversals for the no-slip boundary condition. We argue that the flow reversals with the no-slip and free-slip boundary conditions are different because they are induced by nonlinearities $ ({\bf u} \cdot \nabla) {\bf u}$ and $({\bf u} \cdot \nabla) θ$ respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anisotropic energy transfers in quasi-static magnetohydrodynamic turbulence

We perform direct numerical simulations of quasi-static magnetohydrodynamic turbulence, and compute various energy transfers including the ring-to-ring and conical energy transfers, and the energy fluxes of the perpendicular and parallel components of the velocity field. We show that the rings with higher polar angles transfer energy to ones with lower polar angles. For large interaction parameters, the dominant energy transfer takes place near the equator (polar angle $θ\approx \fracπ{2}$). The energy transfers are local both in wavenumbers and angles. The energy flux of the perpendicular component is predominantly from higher to lower wavenumbers (inverse cascade of energy), while that of the parallel component is from lower to higher wavenumbers (forward cascade of energy). Our results are consistent with earlier results, which indicate quasi two-dimensionalization of quasi-static magnetohydrodynamic (MHD) flows at high interaction parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of reversals and condensates in 2D Kolmogorov flows

We present direct numerical simulations of the different two-dimensional flow regimes generated by a constant spatially periodic forcing balanced by viscous dissipation and large scale drag with a dimensionless damping rate $1/Rh$. The linear response to the forcing is a $6\times6$ square array of counter-rotating vortices, which is stable when the Reynolds number $Re$ or $Rh$ are small. After identifying the sequence of bifurcations that lead to a spatially and temporally chaotic regime of the flow when $Re$ and $Rh$ are increased, we study the transitions between the different turbulent regimes observed for large $Re$ by varying $Rh$. A large scale circulation at the box size (the condensate state) is the dominant mode in the limit of vanishing large scale drag ($Rh$ large). When $Rh$ is decreased, the condensate becomes unstable and a regime with random reversals between two large scale circulations of opposite signs is generated. It involves a bimodal probability density function of the large scale velocity that continuously bifurcates to a Gaussian distribution when $Rh$ is decreased further.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy spectrum of Buoyancy-driven Flows

Using high-resolution direct numerical simulation and arguments based on the kinetic energy flux $Π_u$, we demonstrate that for stably stratified flows, the kinetic energy spectrum $E_u(k) \sim k^{-11/5}$, the entropy spectrum $E_θ(k) \sim k^{-7/5}$, and $Π_u(k) \sim k^{-4/5}$ (Bolgiano-Obukhov scaling). This scaling is due to the depletion of kinetic energy because of buoyancy. For weaker buoyancy in stratified flows, $E_u(k)$ follows Kolmgorov's spectrum with a constant energy flux. We also argue that for Rayleigh Bénard convection, the Bolgiano-Obukhov scaling will not hold for the bulk flow due to the positive energy supply by buoyancy and non-decreasing $Π_u(k)$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy Spectrum of Buoyancy-Driven Turbulence

Using high-resolution direct numerical simulation and arguments based on the kinetic energy flux $Π_u$, we demonstrate that for stably stratified flows, the kinetic energy spectrum $E_u(k) \sim k^{-11/5}$, the entropy spectrum $E_θ(k) \sim k^{-7/5}$, and $Π_u(k) \sim k^{-4/5}$, consistent with the Bolgiano-Obukhov scaling. This scaling arises due to the conversion of kinetic energy to the potential energy by buoyancy. For weaker buoyancy, this conversion is weak, hence $E_u(k)$ follows Kolmogorov's spectrum with a constant energy flux. For Rayleigh Bénard convection, we show that the energy supply rate by buoyancy is positive, which leads to an increasing $Π_u(k)$ with $k$, thus ruling out Bolgiano-Obukhov scaling for the convective turbulence. Our numerical results show that convective turbulence for unit Prandt number exhibits a constant $Π_u(k)$ and $E_u(k) \sim k^{-5/3}$ for a narrow band of wavenumbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy transfers and magnetic energy growth in small-scale dynamo

In this letter we investigate the dynamics of magnetic energy growth in small-scale dynamo by studying energy transfers, mainly energy fluxes and shell-to-shell energy transfers. We perform dynamo simulations for magnetic Prandtl number $\mathrm{Pm}=20$ on $1024^3$ grid using pseudospectral method. We demonstrate that the magnetic energy growth is caused by nonlocal energy transfers from the large-scale or forcing-scale velocity field to small-scale magnetic field. The peak of these energy transfers move towards lower wavenumbers as dynamo evolves, which is the reason why the integral scale of the magnetic field increases with time. The energy transfers $U2U$ (velocity to velocity) and $B2B$ (magnetic to magnetic) are forward and local.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling quasi-static magnetohydrodynamic turbulence with variable energy flux

In quasi-static MHD, experiments and numerical simulations reveal that the energy spectrum is steeper than Kolmogorov's $k^{-5/3}$ spectrum. To explain this observation, we construct turbulence models based on variable energy flux, which is caused by the Joule dissipation. In the first model, which is applicable to small interaction parameters, the energy spectrum is a power law, but with a spectral exponent steeper than -5/3. In the other limit of large interaction parameters, the second model predicts an exponential energy spectrum and flux. The model predictions are in good agreement with the numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling of heat flux and energy spectrum for "very large" Prandtl number convection

Under the limit of infinite Prandtl number, we derive analytical expressions for the large-scale quantities, e.g., Péclet number Pe, Nusselt number Nu, and rms value of the temperature fluctuations $θ_\mathrm{rms}$. We complement the analytical work with direct numerical simulations, and show that $\mathrm{Nu} \sim \mathrm{Ra}^γ$ with $γ\approx (0.30-0.32)$, $\mathrm{Pe} \sim \mathrm{Ra}^η$ with $η\approx (0.57-0.61)$, and $θ_\mathrm{rms} \sim \mathrm{const}$. The Nusselt number is observed to be an intricate function of $\mathrm{Pe}$, $θ_\mathrm{rms}$, and a correlation function between the vertical velocity and temperature. Using the scaling of large-scale fields, we show that the energy spectrum $E_u(k)\sim k^{-13/3}$, which is in a very good agreement with our numerical results. The entropy spectrum $E_θ(k)$ however exhibits dual branches consisting of $k^{-2}$ and $k^0$ spectra; the $k^{-2}$ branch corresponds to the Fourier modes $\hatθ(0,0,2n)$, which are approximately $-1/(2n π)$. The scaling relations for Prandtl number beyond $10^2$ match with those for infinite Prandtl number.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Similarities between 2D and 3D convection for large Prandtl number

Using direct numerical simulations of Rayleigh-Bénard convection (RBC), we perform a comparative study of the spectra and fluxes of energy and entropy, and the scaling of large-scale quantities for large and infinite Prandtl numbers in two (2D) and three (3D) dimensions. We observe close similarities between the 2D and 3D RBC, in particular the kinetic energy spectrum $E_u(k) \sim k^{-13/3}$, and the entropy spectrum exhibits a dual branch with a dominant $k^{-2}$ spectrum. We showed that the dominant Fourier modes in the 2D and 3D flows are very close. Consequently, the 3D RBC is quasi two-dimensional, which is the reason for the similarities between the 2D and 3D RBC for large- and infinite Prandtl numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong anisotropy in quasi-static magnetohydrodynamic turbulence for high interaction parameters

We simulate forced quasi-static magnetohydrodynamic turbulence and investigate the anisotropy, energy spectrum, and energy flux of the flow, specially for large interaction parameters ($N$). We show that the angular dependence of the energy spectrum is well quantified using Legendre polynomials. For large $N$, the energy spectrum is exponential. Our direct computation of energy flux reveals an inverse cascade of energy at low wavenumbers, similar to that in two-dimensional turbulence. We observe the flow be two-dimensional (2D) for moderate $N$ ($N \sim 20$), and two-dimensional three-component (2D-3C) type for $N \ge 27$. In our forced simulation, the transition from 2D to 2D-3C occurs at higher value of $N$ than Favier et al., ["On the two-dimensionalization of quasistatic magnetohydrodynamic turbulence," Phys. Fluids 22, 075104 (2010)] who employ decaying simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamo in Protostar

In this paper, we estimate the magnetic Reynolds number of a typical protostar before and after deuterium burning, and claim for the existence of dynamo process in both the phases, because the magnetic Reynolds number of the protostar far exceeds the critical magnetic Reynolds number for dynamo action. Using the equipartition of kinetic and magnetic energies, we estimate the steady-state magnetic field of the protostar to be of the order of kilo-gauss, which is in good agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Supercriticality to subcriticality in dynamo transitions

Evidence from numerical simulations suggest that the nature of dynamo transition changes from supercritical to subcritical as the magnetic Prandtl number is decreased. To explore this interesting crossover we first use direct numerical simulations to investigate the hysteresis zone of a subcritical Taylor-Green dynamo. We establish that a well defined boundary exists in this hysteresis region which separates dynamo states from the purely hydrodynamic solution. We then propose simple dynamo models which show similar crossover from supercritical to subcritical dynamo transition as a function of the magnetic Prandtl number. Our models show that the change in the nature of dynamo transition is connected to the stabilizing or de-stabilizing influence of governing non-linearities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variable energy flux in quasi-static magnetohydrodynamic turbulence

Experiments and numerical simulations show that the energy spectrum of the magnetohydrodynamic turbulence in the quasi-static limit deviates from Kolmogorov's $k^{-5/3}$ spectrum as the external magnetic field, or equivalently the interaction parameter, is increased. To explain this phenomena, we construct an analytical turbulence model with variable energy flux that arises due to the Lorentz-force induced dissipation. The energy spectra computed using our model for various interaction parameters are in qualitative agreement with earlier experimental and numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Benchmarking and scaling studies of a pseudospectral code Tarang for turbulence simulations

Tarang is a general-purpose pseudospectral parallel code for simulating flows involving fluids, magnetohydrodynamics, and Rayleigh-Bénard convection in turbulence and instability regimes. In this paper we present code validation and benchmarking results of Tarang. We performed our simulations on $1024^3$, $2048^3$, and $4096^3$ grids using the {\em HPC system} of IIT Kanpur and {\em Shaheen} of KAUST. We observe good "weak" and "strong" scaling for Tarang on these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Field correlations and the ultimate regime of turbulent convection

Using direct numerical simulations of Rayleigh-Bénard convection (RBC) under free-slip boundary condition, we show that the normalized correlation function between the vertical velocity field and the temperature field, as well as the normalized viscous dissipation rate, scales as $Ra^{-0.22}$ for moderately large Rayleigh number $Ra$. This scaling accounts for the Nusselt number ($Nu$) exponent to be around 0.3 observed in experiments. Numerical simulations also reveal that the above normalized correlation functions are constants for the convection simulation under periodic boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flow reversals in turbulent convection via vortex reconnections

We employ detailed numerical simulations to probe the mechanism of flow reversals in two-dimensional turbulent convection. We show that the reversals occur via vortex reconnection of two attracting corner rolls having same sign of vorticity, thus leading to major restructuring of the flow. Large fluctuations in heat transport are observed during the reversal due to this flow reconfiguration. The flow configurations during the reversals have been analyzed quantitatively using large-scale modes. Using these tools, we also show why flow reversals occur for a restricted range of Rayleigh and Prandt numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Variable enstrophy flux and energy spectrum in two-dimensional turbulence with Ekman friction

Experiments and numerical simulations reveal that in the forward cascade regime, the energy spectrum of two-dimensional turbulence with Ekman friction deviates from Kraichnan's prediction of $k^{-3}$ power spectrum. In this letter we explain this observation using an analytic model based on variable enstrophy flux arising due to Ekman friction. We derive an expression for the enstrophy flux which exhibits a logarithmic dependence in the inertial range for the Ekman-friction dominated flows. The energy spectrum obtained using this enstrophy flux shows a power law scaling for large Reynolds number and small Ekman friction, but has an exponential behaviour for large Ekman friction and relatively small Reynolds number.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bistability and chaos in Taylor-Green dynamo

Using direct numerical simulations we study dynamo action under the Taylor-Green forcing with Prandtl number less than one. We observe bistability with a weak magnetic field branch and a strong magnetic field branch. Both the dynamo branches undergo subcritical dynamo transition. We also observe host of dynamo states including constant, periodic, quasiperiodic, and chaotic magnetic fields. One of the chaotic state originates through a quasiperiodic route with phase locking, while another chaotic attractor appears to follow Newhouse-Ruelle-Takens route to chaos. We also observe intermittent transitions among quasiperiodic and chaotic states for a given Taylor-Green forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics and symmetries of flow reversals in turbulent convection

Based on direct numerical simulations and symmetry arguments, we show that the large-scale Fourier modes are useful tools to describe the flow structures and dynamics of flow reversals in Rayleigh-Bénard convection (RBC). We observe that during the reversals, the amplitude of one of the large-scale modes vanishes, while another mode rises sharply, very similar to the "cessation-led" reversals observed earlier in experiments and numerical simulations. We find anomalous fluctuations in the Nusselt number during the reversals. Using the structures of the RBC equations in the Fourier space, we deduce two symmetry transformations that leave the equations invariant. These symmetry transformations help us in identifying the reversing and non-reversing Fourier modes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamics of zero-Prandtl number convection near the onset

In this paper we present various convective states of zero-Prandtl number Rayleigh-Bénard convection using direct numerical simulations (DNS) and a 27-mode low-dimensional model containing the energetic modes of DNS. The origin of these convective states have been explained using bifurcation analysis. The system is chaotic at the onset itself with three coexisting chaotic attractors that are born at two codimension-2 bifurcation points. One of the bifurcation points with a single zero eigenvalue and a complex pair $(0,\pm i ω)$ generates chaotic attractors and associated periodic, quasiperiodic, and phase-locked states that are related to the wavy rolls observed in experiments and simulations. The frequency of the wavy rolls are in general agreement with $ω$ of the above eigenvalue of the stability matrix. The other bifurcation point with a double zero eigenvalue produces the other set of chaotic attractors and ordered states such as squares, asymmetric squares, oscillating asymmetric squares, relaxation oscillations with intermediate squares, some of which are common to the 13-mode model of Pal et al.(2009).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy spectra and fluxes for Rayleigh-Benard convection

We compute the spectra and fluxes of the velocity and temperature fields in Rayleigh-Bénard convection in turbulent regime for a wide range of Prandtl numbers using pseudo-spectral simulations on $512^3$ grids. Our spectral and flux results support the Kolmogorov-Obukhov (KO) scaling for zero Prandtl number and low Prandtl number ($P=0.02$) convection. The KO scaling for the velocity field in zero-Prandtl number and low-Prandtl number convection is because of the weak buoyancy in the inertial range (buoyancy is active only at the very low wavenumbers). We also observe that for intermediate Prandtl numbers ($P=0.2$) the KO scaling fits better with the numerical results than the Bolgiano-Obukhov (BO) scaling. For large Prandtl number ($P=6.8$), the spectra and flux results are somewhat inconclusive on the validity of the KO or BO scaling, yet the BO scaling is preferred over the KO scaling for these cases. The numerical results for P=1 is rather inconclusive.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Bifurcation and chaos in zero Prandtl number convection

We present the detailed bifurcation structure and associated flow patterns near the onset of zero Prandtl number Rayleigh Bénard convection. We employ both direct numerical simulation and a low-dimensional model ensuring qualitative agreement between the two. Various flow patterns originate from a stationary square observed at a higher Rayleigh number through a series of bifurcations starting from a pitchfork followed by a Hopf and finally a homoclinic bifurcation as the Rayleigh number is reduced to the critical value. Global chaos, intermittency, and crises are observed near the onset.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Patterns and bifurcations in low-Prandtl number Rayleigh-Benard convection

We present a detailed bifurcation structure and associated flow patterns for low-Prandtl number ($P=0.0002, 0.002, 0.005, 0.02$) Rayleigh-Bénard convection near its onset. We use both direct numerical simulations and a 30-mode low-dimensional model for this study. We observe that low-Prandtl number (low-P) convection exhibits similar patterns and chaos as zero-P convection \cite{pal:2009}, namely squares, asymmetric squares, oscillating asymmetric squares, relaxation oscillations, and chaos. At the onset of convection, low-P convective flows have stationary 2D rolls and associated stationary and oscillatory asymmetric squares in contrast to zero-P convection where chaos appears at the onset itself. The range of Rayleigh number for which stationary 2D rolls exist decreases rapidly with decreasing Prandtl number. Our results are in qualitative agreement with results reported earlier.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Energy transfer and locality in magnetohydrodynamic turbulence

The shell-to-shell energy transfer rates for magnetohydrodynamic (MHD) turbulence are computed analytically, which shows local energy transfer rates from velocity to velocity, velocity to magnetic, magnetic to velocity, and magnetic to magnetic fields for nonhelical MHD in the inertial range. It is also found that for kinetic-energy dominated MHD fluid, there is a preferential shell-to-shell energy transfer from kinetic to magnetic energy; the transfer is reversed for magnetic-energy dominated MHD fluid. This property is the reason for the asymptotic value of Alfven ratio to be close to 0.5. The analytical results are in close agreement with recent numerical results. When magnetic and kinetic helicities are turned on, the helical contributions are opposite to the corresponding nonhelical contributions. The helical energy transfers have significant nonlocal components.

0 citations0 reviewsNo code linkedOpen access