Dhrubaditya Mitra

Convective dynamos in spherical wedge geometry

Astronomy Unit, School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London E14NS, United KingdomReceived 2009 Sep 9, accepted 2009 Nov 16Published online 2009 Dec 30Key words Sun: magnetic fields – magnetohydrodynamics (MHD)Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation isgenerated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained onlyfor sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields arefound to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesiandomains is made.

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.

Dynamo-driven plasmoid ejections above a spherical surface

Aims: We extend earlier models of turbulent dynamos with an upper, nearly force-free exterior to spherical geometry, and study how flux emerges from lower layers to the upper ones without being dri ...Magnetic buoyancy is believed to drive the transport of magnetic flux tubes from the convection zone to the surface of the Sun. The magnetic fields form twisted loop-like structures in the solar atmosphere. In this paper we use helical forcing to produce a large-scale dynamo-generated magnetic field, which rises even without magnetic buoyancy. A two layer system is used as computational domain where the upper part represents the solar atmosphere. Here, the evolution of the magnetic field is solved with the stress--and--relax method. Below this region a magnetic field is produced by a helical forcing function in the momentum equation, which leads to dynamo action. We find twisted magnetic fields emerging frequently to the outer layer, forming arch-like structures. In addition, recurrent plasmoid ejections can be found by looking at space--time diagrams of the magnetic field. Recent simulations in spherical coordinates show similar results.

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

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.

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 timescales, 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.

Is Multiscaling an Artifact in the Stochastically Forced Burgers Equation

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

Alpha effect and diffusivity in helical turbulence with shear

Aims. We study the dependence of turbulent transport coefficients, such as the components of the α tensor (αij) and the turbulent magnetic diffusivity tensor (ηij), on shear and magnetic Reynolds number in the presence of helical forcing. Methods. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. Results. We find that with increasing shear the diagonal components of αij quench, whereas those of ηij increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or γ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of αij, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

Active Region Formation through the Negative Effective Magnetic Pressure Instability

The negative effective magnetic-pressure instability operates on scales encompassing many turbulent eddies, which correspond to convection cells in the Sun. This instability is discussed here 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, for example the onset of the instability occurs at the same depth. This depth increases with increasing field strength, such that the 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.

Spontaneous Formation of Magnetic Flux Concentrations in Stratified Turbulence

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNSs) 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 the formation of large-scale inhomogeneous magnetic structures from an initially uniform magnetic field. Earlier investigations of this instability in DNSs of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma β 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, L2/ηt, where ηt is the turbulent diffusivity and L is the scale of the domain. DNSs 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.