Source author record

Dhrubaditya Mitra

Dhrubaditya Mitra 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

47works
18topics
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

47 published item(s)

preprint2022arXiv

Asymptotic interpretation of the Miles mechanism of wind-wave instability

When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air-water density ratio, $ρ_{\rm{a}}/ρ_{\rm{w}}\equivε<<1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $ε$, the eigenvalue has a finite imaginary part, which implies growth. Miles (1957) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically.

preprint2022arXiv

Chaos and irreversibility of a flexible filament in periodically-driven Stokes flow

The flow of Newtonian fluid at low Reynolds number is, in general, regular and time-reversible due to absence of nonlinear effects. For example, if the fluid is sheared by its boundary motion that is subsequently reversed, then all the fluid elements return to their initial positions. Consequently, mixing in microchannels happens solely due to molecular diffusion and is very slow. Here, we show, numerically, that the introduction of a single, freely-floating, flexible filament in a time-periodic linear shear flow can break reversibility and give rise to chaos due to elastic nonlinearities, if the bending rigidity of the filament is within a carefully chosen range. Within this range, not only the shape of the filament is spatiotemporally chaotic, but also the flow is an efficient mixer. Overall, we find five dynamical phases: the shape of a stiff filament is time-invariant -- either straight or buckled; it undergoes a period-two bifurcation as the filament is made softer; becomes spatiotemporally chaotic for even softer filaments but, surprisingly, the chaos is suppressed if bending rigidity is decreased further.

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.

preprint2022arXiv

Exact analytical solutions in inhomogeneous magnetic fields for linear asteroseismic waves

We solve for waves in an isothermal, stratified medium with a magnetic field that points along a direction perpendicular to that of gravity and varies exponentially in the direction of gravity. We find exact analytical solutions for two different cases: (a) waves propagating along the direction of the magnetic field and (b) waves propagating along the direction of the gravity. In each case, we find solutions in terms of either the hypergeometric functions or their confluent cousins. We solve the resultant transcendental dispersion relation numerically. The changes in the eigenfrequencies due to a spatially inhomogeneous Alfvén wave speed are significant -- there are more oscillations or leaking of nodes of the eigenfunctions in regions of lower Alfvén speed and wave reflection from regions of larger Alfvén speed. Such changes in the dispersion relation and the mode structures may allow detection of magnetic fields buried in the stellar interior.

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

A model-free, data-based forecast for sunspot cycle 25

The dynamic activity of the Sun, governed by its cycle of sunspots -- strongly magnetized regions that are observed on its surface -- modulate our solar system space environment creating space weather. Severe space weather leads to disruptions in satellite operations, telecommunications, electric power grids and air-traffic on polar routes. Forecasting the cycle of sunspots, however, has remained a challenging problem. We use reservoir computing -- a model-free, neural--network based machine-learning technique -- to forecast the upcoming solar cycle, sunspot cycle 25. The standard algorithm forecasts that solar cycle 25 is going to last about ten years, the maxima is going to appear in the year 2024 and the maximum number of sunspots is going to be 113 ($\pm15$). We also develop a novel variation of the standard algorithm whose forecasts for duration and peak timing matches that of the standard algorithm, but whose peak amplitude forecast is 124 ($\pm2$) -- within the upper bound of the standard reservoir computing algorithm. We conclude that sunspot cycle 25 is likely to be a weak, lower than average solar cycle, somewhat similar in strength to sunspot cycle 24.

preprint2020arXiv

Erosion of planetesimals by gas flow

The first stages of planet formation take place in protoplanetary disks that are largely made up of gas. Understanding how the gas affects planetesimals in the protoplanetary disk is therefore essential. In this paper, we discuss whether or not gas flow can erode planetesimals. We estimate how much shear stress is exerted onto the planetesimal surface by the gas as a function of disk and planetesimal properties. To determine whether erosion can take place, we compare this with previous measurements of the critical stress that a pebble-pile planetesimal can withstand before erosion begins. If erosion takes place, we estimate the erosion time of the affected planetesimals. We also illustrate our estimates with two-dimensional numerical simulations of flows around planetesimals using the lattice Boltzmann method. We find that the wall shear stress can overcome the critical stress of planetesimals in an eccentric orbit within the innermost regions of the disk. The high eccentricities needed to reach erosive stresses could be the result of shepherding by migrating planets. We also find that if a planetesimal erodes, it does so on short timescales. For planetesimals residing outside of $1 \ \rm{au}$, we find that they are mainly safe from erosion, even in the case of highly eccentric orbits.

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.

preprint2016arXiv

The effect of turbulent clustering on particle reactivity

The effect of turbulence on the heterogeneous (solid-fluid) reactions of solid particles is studied numerically with Direct Numerical Simulations (DNS). A simplified reaction system is used, where the solid-fluid reaction is represented by a single isothermal reaction step. It is found that, due to the clustering of particles by the isotropic turbulence, the overall reaction rate is entirely controlled by the turbulence for large Damköhler numbers. The particle clustering significantly slows down the reaction rate for increasing Damköhler numbers which reaches an asymptotic limit that can be analytically derived. This implies that the effect of turbulence on heterogeneously reacting particles should be included in models that are used in CFD simulations of e.g. char burnout in combustors or gasifiers. Such a model, based on the chemical and turbulent time scales, is here proposed for the heterogeneous reaction rate in the presence of turbulence.

preprint2016arXiv

Turbulent reconnection of magnetic bipoles in stratified turbulence

We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S>10^3, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings.

preprint2015arXiv

Bipolar magnetic spots from dynamos in stratified spherical shell turbulence

Recent work by Mitra et al. (2014) has shown that in strongly stratified forced two-layer turbulence with helicity and corresponding large-scale dynamo action in the lower layer, a magnetic field occurs in the upper layer in the form of sharply bounded bipolar magnetic spots. Here we extend this model to spherical wedge geometry covering the northern hemisphere up to 75° latitude and an azimuthal extent of 180°. The kinetic helicity and therefore also the large-scale magnetic field are strongest at low latitudes. For moderately strong stratification, several bipolar spots form that fill eventually the full longitudinal extent. At early times, the polarity of spots reflects the orientation of the underlying azimuthal field, as expected from Ω-shaped flux loops. At late times their tilt changes such that there is a radial field of opposite orientation at different latitudes separated by about 10°. Our model demonstrates the spontaneous formation of spots of sizes much larger than the pressure scale height. Their tendency to produce filling factors close to unity is argued to be reminiscent of highly active stars. We confirm that strong stratification and strong scale separation are essential ingredients behind magnetic spot formation, which appears to be associated with downflows at larger depths.

preprint2015arXiv

Rheology of confined non-Brownian suspensions

We study the rheology of confined suspensions of neutrally buoyant rigid monodisperse spheres in plane-Couette flow using Direct Numerical Simulations. We find that if the width of the channel is a (small) integer multiple of the sphere's diameter, the spheres self-organize into two-dimensional layers that slide on each other and the suspension's effective viscosity is significantly reduced. Each two-dimensional layer is found to be structurally liquid-like but their dynamics is frozen in time.

preprint2014arXiv

A microfluidic device to sort capsules by deformability: A numerical study

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device.

preprint2014arXiv

Ensemble Kalman filter data assimilation in a Babcock-Leighton solar dynamo model: an observation system simulation experiment for reconstructing meridional flow-speed

Accurate knowledge of time-variation in meridional flow-speed and profile is crucial for estimating a solar cycle's features, which are ultimately responsible for causing space climate variations. However, no consensus has been reached yet about the Sun's meridional circulation pattern observations and theories. By implementing an Ensemble Kalman Filter (EnKF) data assimilation in a Babcock-Leighton solar dynamo model using Data Assimilation Research Testbed (DART) framework, we find that the best reconstruction of time-variation in meridional flow-speed can be obtained when ten or more observations are used with an updating time of 15 days and a $\le 10\%$ observational error. Increasing ensemble-size from 16 to 160 improves reconstruction. Comparison of reconstructed flow-speed with "true-state" reveals that EnKF data assimilation is very powerful for reconstructing meridional flow-speeds and suggests that it can be implemented for reconstructing spatio-temporal patterns of meridional circulation.

preprint2014arXiv

Intense bipolar structures from stratified helical dynamos

We perform direct numerical simulations of the equations of magnetohydrodynamics with external random forcing and in the presence of gravity. The domain is divided into two parts: a lower layer where the forcing is helical and an upper layer where the helicity of the forcing is zero with a smooth transition in between. At early times, a large-scale helical dynamo develops in the bottom layer. At later times the dynamo saturates, but the vertical magnetic field continues to develop and rises to form dynamic bipolar structures at the top, which later disappear and reappear. Some of the structures look similar to $δ$ spots observed in the Sun. This is the first example of magnetic flux concentrations, owing to strong density stratification, from self-consistent dynamo simulations that generate bipolar, super-equipartition strength, magnetic structures whose energy density can exceeds the turbulent kinetic energy by even a factor of ten.

preprint2014arXiv

Particle energization through time-periodic helical magnetic fields

We solve for the motion of charged particles in a helical time-periodic ABC (Arnold-Beltrami-Childress) magnetic field. The magnetic field lines of a stationary ABC field with coefficients $A=B=C=1$ are chaotic, and we show that the motion of a charged particle in such a field is also chaotic at late times with positive Lyapunov exponent. We further show that in time-periodic ABC fields, the kinetic energy of a charged particle can increase indefinitely with time. At late times the mean kinetic energy grows as a power law in time with an exponent that approaches unity. For an initial distribution of particles, whose kinetic energy is uniformly distributed within some interval, the PDF of kinetic energy is, at late times, close to a Gaussian but with steeper tails.

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

Surface flux concentrations and spherical alpha-square dynamo

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $α^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $α$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $α$ effect dynamo.

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

A mean field dynamo from negative eddy diffusivity

Using direct numerical simulations, we verify that "flow IV" of Roberts (1972) exhibits dynamo action dominated by horizontally averaged large-scale magnetic field. With the test-field method we compute the turbulent magnetic diffusivity and find that it is negative and overcomes the molecular diffusivity, thus explaining quantitatively the large-scale dynamo for magnetic Reynolds numbers above $\approx8$. As expected for a dynamo of this type, but contrary to $α$-effect dynamos, the two horizontal field components grow independently of each other and have arbitrary amplitude ratios and phase differences. Small length scales of the mean magnetic field are shown to be stabilized by the turbulent magnetic diffusivity becoming positive at larger wavenumbers. Oscillatory decaying or growing solutions have also been found in certain wavenumber intervals and sufficiently large values of the magnetic Reynolds number. For magnetic Reynolds numbers below $\approx0.5$ the turbulent magnetic diffusivity is confirmed to be positive, as expected for all incompressible flows. Earlier claims of a dynamo driven by a modified Taylor-Green flow through negative eddy diffusivity could not be confirmed.

preprint2013arXiv

Can planetesimals form by collisional fusion?

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales.

preprint2013arXiv

Quantifying the effect of turbulent magnetic diffusion on the growth rate of the magneto-rotational instability

In astrophysics, turbulent diffusion is often used in place of microphysical diffusion to avoid resolving the small scales. However, we expect this approach to break down when time and length scales of the turbulence become comparable with other relevant time and length scales in the system. Turbulent diffusion has previously been applied to the magneto-rotational instability (MRI), but no quantitative comparison of growth rates at different turbulent intensities has been performed. We investigate to what extent turbulent diffusion can be used to model the effects of small-scale turbulence on the kinematic growth rates of the MRI, and how this depends on angular velocity and magnetic field strength. We use direct numerical simulations in three-dimensional shearing boxes with periodic boundary conditions in the spanwise direction and additional random plane-wave volume forcing to drive a turbulent flow at a given length scale. We estimate the turbulent diffusivity using a mixing length formula and compare with results obtained with the test-field method. It turns out that the concept of turbulent diffusion is remarkably accurate in describing the effect of turbulence on the growth rate of the MRI. No noticeable breakdown of turbulent diffusion has been found, even when time and length scales of the turbulence become comparable with those imposed by the MRI itself. On the other hand, quenching of turbulent magnetic diffusivity by the magnetic field is found to be absent. Turbulence reduces the growth rate of the MRI in a way that is the same as microphysical magnetic diffusion.

preprint2013arXiv

Scaling and intermittency in incoherent α-shear dynamo

We consider mean-field dynamo models with fluctuating αeffect, both with and without shear. The αeffect is chosen to be Gaussian white noise with zero mean and given covariance. We show analytically that the mean magnetic field does not grow, but, in an infinitely large domain, the mean-squared magnetic field shows exponential growth of the fastest growing mode at a rate proportional to the shear rate, which agrees with earlier numerical results of Yousef et al (2008) and recent analytical treatment by Heinemann et al (2011) who use a method different from ours. In the absence of shear, an incoherent α^2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third and fourth order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian.

preprint2013arXiv

Shear thickening in non-Brownian suspensions: an excluded volume effect

Shear thickening appears as an increase of the viscosity of a dense suspension with the shear rate, sometimes sudden and violent at high volume fraction. Its origin for noncolloidal suspension with non-negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate that fluid inertia causes a strong microstructure anisotropy that results in the formation of a shadow region with no relative flux of particles. We show that shear thickening at finite inertia can be explained as an increase of the effective volume fraction when considering the dynamically excluded volume due to these shadow regions

preprint2012arXiv

Active region formation through the negative effective magnetic pressure instability

The negative effective magnetic pressure instability operates on scales encompassing many turbulent eddies and is here discussed in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details in that their onset occurs at the same depth. This depth increases with increasing field strength, such that the maximum growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.

preprint2012arXiv

Breakdown of chiral symmetry during saturation of the Tayler instability

We study spontaneous breakdown of chiral symmetry during the nonlinear evolution of the Tayler instability. We start with an initial steady state of zero helicity. Within linearized perturbation calculations, helical perturbations of this initial state have the same growth rate for either sign of helicity. Direct numerical simulations (DNS) of the fully nonlinear equations, however, show that an infinitesimal excess of one sign of helicity in the initial perturbation gives rise to a saturated helical state. We further show that this symmetry breaking can be described by weakly nonlinear finite--amplitude equations with undetermined coefficients which can be deduced solely from symmetry consideration. By fitting solutions of the amplitude equations to data from DNS we further determine the coefficients of the amplitude equations.

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

Dynamics of Saturated Energy Condensation in Two-Dimensional Turbulence

In two-dimensional forced Navier-Stokes turbulence, energy cascades to the largest scales in the system to form a pair of coherent vortices known as the Bose condensate. We show, both numerically and analytically, that the energy condensation saturates and the system reaches a statistically stationary state. The time scale of saturation is inversely proportional to the viscosity and the saturation energy level is determined by both the viscosity and the force. We further show that, without sufficient resolution to resolve the small-scale enstrophy spectrum, numerical simulations can give a spurious result for the saturation energy level. We also find that the movement of the condensate is similar to the motion of an inertial particle with an effective drag force. Furthermore, we show that the profile of the saturated coherent vortices can be described by a Gaussian core with exponential wings.

preprint2012arXiv

Magnetic twist: a source and property of space weather

We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius in the northern hemisphere.

preprint2012arXiv

Rotational effects on the negative magnetic pressure instability

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

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.

preprint2011arXiv

Alleviating alpha quenching by solar wind and meridional flow

We study the ability of magnetic helicity expulsion to alleviate catastrophic $α$-quenching in mean field dynamos in two--dimensional spherical wedge domains. Motivated by the physical state of the outer regions of the Sun, we consider $α^2Ω$ mean field models with a dynamical $α$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow ("solar wind") and meridional circulation. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain parameter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field.

preprint2011arXiv

Detection of negative effective magnetic pressure instability in turbulence simulations

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

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

Magnetic helicity transport in the advective gauge family

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

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

preprint2011arXiv

Spontaneous chiral parity breaking by hydromagnetic buoyancy

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.

preprint2011arXiv

Spontaneous formation of flux concentrations in a stratified layer

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.

preprint2010arXiv

Alpha effect due to buoyancy instability of a magnetic layer

A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.

preprint2010arXiv

Oscillatory migrating magnetic fields in helical turbulence in spherical domains

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.

preprint2010arXiv

Plasmoid ejections driven by dynamo action underneath a spherical surface

We present a unified three-dimensional model of the convection zone and upper atmosphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity of the small--scale field is transported outwards and form large structures like magnetic clouds.

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.

preprint2009arXiv

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

preprint2009arXiv

Turbulent dynamos in spherical shell segments of varying geometrical extent

We use three-dimensional direct numerical simulations of the helically forced magnetohydrodynamic equations in spherical shell segments in order to study the effects of changes in the geometrical shape and size of the domain on the growth and saturation of large-scale magnetic fields. We inject kinetic energy along with kinetic helicity in spherical domains via helical forcing using Chandrasekhar-Kendall functions. We take perfect conductor boundary conditions for the magnetic field to ensure that no magnetic helicity escapes the domain boundaries. We find dynamo action giving rise to magnetic fields at scales larger than the characteristic scale of the forcing. The magnetic energy exceeds the kinetic energy over dissipative time scales, similar to that seen earlier in Cartesian simulations in periodic boxes. As we increase the size of the domain in the azimuthal direction we find that the nonlinearly saturated magnetic field organizes itself in long-lived cellular structures with aspect ratios close to unity. These structures tile the domain along the azimuthal direction, thus resulting in very small longitudinally averaged magnetic fields for large domain sizes. The scales of these structures are determined by the smallest scales of the domain, which in our simulations is usually the radial scale. We also find that increasing the meridional extent of the domains produces little qualitative change, except a marginal increase in the large-scale field. We obtain qualitatively similar results in Cartesian domains with similar aspect ratios.

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.