Source author record

A. Brandenburg

A. Brandenburg 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

65works
14topics
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

65 published item(s)

preprint2021arXiv

Spectral characterisation of inertial particle clustering in turbulence

Clustering of inertial particles is important for many types of astrophysical and geophysical turbulence, but it has been studied predominately for incompressible flows. Here we study compressible flows and compare clustering in both compressively (irrotationally) and vortically (solenoidally) forced turbulence. Vortically and compressively forced flows are driven stochastically either by solenoidal waves or by circular expansion waves, respectively. For compressively forced flows, the power spectrum of the density of inertial particles is a useful tool for displaying particle clustering relative to the fluid density enhancement. Power spectra are shown to be particularly sensitive for studying large-scale particle clustering, while conventional tools such as radial distribution functions are more suitable for studying small-scale clustering. Our primary finding is that particle clustering through shock interaction is particularly prominent in turbulence driven by spherical expansion waves. It manifests itself through a double-peaked distribution of spectral power as a function of Stokes number. The two peaks are associated with two distinct clustering mechanisms; shock interaction for smaller Stokes numbers and the centrifugal sling effect for larger values. The clustering of inertial particles is associated with the formation of caustics. Such caustics can only be captured in the Lagrangian description, which allows us to assess the relative importance of caustics in vortically and compressively forced turbulence. We show that the statistical noise resulting from the limited number of particles in the Lagrangian description can be removed from the particle power spectra, allowing us a more detailed comparison of the residual spectra. We focus on the Epstein drag law relevant for rarefied gases, but show that our findings apply also to the usual Stokes drag.

preprint2020arXiv

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.

preprint2016arXiv

Robustness of oscillatory $α^2$ dynamos in spherical wedges

We study the connection between spherical wedge and full spherical shell geometries using simple mean-field $α^2$ dynamos. We solve the equations for a one-dimensional time-dependent mean-field dynamo to examine the effects of varying the polar angle $θ_0$ between the latitudinal boundaries and the poles in spherical coordinates. We investigate the effects of turbulent magnetic diffusivity and $α$ effect profiles as well as different latitudinal boundary conditions to isolate parameter regimes where oscillatory solutions are found. Finally, we add shear along with a damping term mimicking radial gradients to study the resulting dynamo regimes. We find that the commonly used perfect conductor boundary condition leads to oscillatory $α^2$ dynamo solutions only if the wedge boundary is at least one degree away from the poles. Other boundary conditions always produce stationary solutions. By varying the profile of the turbulent magnetic diffusivity alone, oscillatory solutions are achieved with models extending to the poles, but the magnetic field is strongly concentrated near the poles and the oscillation period is very long. By introducing radial shear and a damping term mimicking radial gradients, we again see oscillatory dynamos, and the direction of drift follows the Parker--Yoshimura rule. Oscillatory solutions in the weak shear regime are found only in the wedge case with $θ_0 = 1^\circ$ and perfect conductor boundaries. A reduced $α$ effect near the poles with a turbulent diffusivity concentrated toward the equator yields oscillatory dynamos with equatorward migration and reproduces best the solutions in spherical wedges.

preprint2015arXiv

Magnetic flux concentrations from turbulent stratified convection

(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.

preprint2014arXiv

Confirmation of bistable stellar differential rotation profiles

(abridged) Context: Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can result in a slower equator and faster poles when the rotation is slow. Aims: We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one (slow equator, fast poles). We estimate the non-diffusive ($Λ$ effect) and diffusive (turbulent viscosity) contributions to the Reynolds stress. Methods: We present the results of three-dimensional numerical simulations of mildly turbulent convection in spherical wedge geometry. Here we apply a fully compressible setup which would suffer from a prohibitive time step constraint if the real solar luminosity was used. We regulate the convective velocities by varying the amount of heat transported by thermal conduction, turbulent diffusion, and resolved convection. Results: Increasing the efficiency of resolved convection leads to a reduction of the rotational influence on the flow and a sharp transition from solar-like to anti-solar differential rotation for Coriolis numbers around 1.3. We confirm the recent finding of a large-scale flow bistability: contrasted with running the models from an initial condition with unprescribed differential rotation, the initialization of the model with certain kind of rotation profile sustains the solution over a wider parameter range. Conclusions: Our results may have implications for real stars that start their lives as rapid rotators implying solar-like rotation in the early main-sequence evolution. As they slow down, they might be able to retain solar-like rotation for lower Coriolis numbers before switching to anti-solar rotation. This could partially explain the puzzling findings of anti-solar rotation profiles for models in the solar parameter regime.

preprint2014arXiv

Dynamical quenching with non-local alpha and downward pumping

In light of new results, the one-dimensional mean-field dynamo model of Brandenburg & Kapyla (2007) with dynamical quenching and a nonlocal Babcock-Leighton alpha effect is re-examined for the solar dynamo. We extend the one-dimensional model to include the effects of turbulent downward pumping (Kitchatinov & Olemskoy 2011), and to combine dynamical quenching with shear. We use both the conventional dynamical quenching model of Kleeorin & Ruzmaikin (1982) and the alternate one of Hubbard & Brandenburg (2011), and confirm that with varying levels of non-locality in the alpha effect, and possibly shear as well, the saturation field strength can be independent of the magnetic Reynolds number.

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

Magnetic flux concentrations in a polytropic atmosphere

Strongly stratified hydromagnetic turbulence has recently been identified as a candidate for explaining the spontaneous formation of magnetic flux concentrations by the negative effective magnetic pressure instability (NEMPI). Much of this work has been done for isothermal layers for which the density scale height is constant throughout. We now study the validity of earlier conclusions about the size and growth rate of magnetic structures in the case of polytropic layers, which scale height decreases sharply towards the surface. To allow for a continuous transition from isothermal to polytropic layers, we employ a generalization of the exponential function known as the q-exponential. Now, the top of the polytropic layer shifts with the polytropic index such that the scale height at some reference height is always the same. We use both mean-field and direct numerical simulations of forced stratified turbulence to determine the resulting flux concentrations. Magnetic structures begin to form at a depth where magnetic field strength is about 3-4% the local equipartition field strength with respect to the turbulent kinetic energy. Unlike the isothermal case where stronger fields can give rise to magnetic flux concentrations at larger depths, in the polytropic one the growth rates decreases for structures deeper down. For vertical fields, magnetic structures of super-equipartition strengths are formed because such fields survive downward advection, unlike NEMPI under horizontal magnetic fields. The horizontal cross-section of such structures is approximately circular. Results based on isothermal models can be applied locally to polytropic layers. For vertical fields, magnetic flux concentrations of super-equipartition strengths form, which supports suggestions that sunspot formation might be a shallow phenomenon.

preprint2014arXiv

Mean-field and direct numerical simulations of magnetic flux concentrations from vertical field

Strongly stratified hydromagnetic turbulence has previously been found to produce magnetic flux concentrations if the domain is large enough compared with the size of turbulent eddies. Mean-field simulations (MFS) using parameterizations of the Reynolds and Maxwell stresses show a negative effective magnetic pressure instability and have been able to reproduce many aspects of direct numerical simulations (DNS) regarding the growth rate of this large-scale instability, shape of the resulting magnetic structures, and their height as a function of magnetic field strength. Unlike the case of an imposed horizontal field, for a vertical one, magnetic flux concentrations of equipartition strength with the turbulence can be reached. This results in magnetic spots that are reminiscent of sunspots. Here we want to find out under what conditions magnetic flux concentrations with vertical field occur and what their internal structure is. We use a combination of MFS, DNS, and implicit large-eddy simulations to characterize the resulting magnetic flux concentrations in forced isothermal turbulence with an imposed vertical magnetic field. We confirm earlier results that in the kinematic stage of the large-scale instability the horizontal wavelength of structures is about 10 times the density scale height. At later times, even larger structures are being produced in a fashion similar to inverse spectral transfer in helically driven turbulence. Using turbulence simulations, we find that magnetic flux concentrations occur for different values of the Mach number between 0.1 and 0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are being confined by downflows along the tubes and corresponding inflow from the sides, which keep the field concentrated.

preprint2014arXiv

The alpha-effect in a turbulent liquid-metal plane Couette flow

We calculate the mean electromotive force in plane Couette flows of a nonrotating conducting fluid under the influence of a large-scale magnetic field for driven turbulence. A vertical stratification of the turbulence intensity results in an alpha effect owing to the presence of horizontal shear. Here we discuss the possibility of an experimental determination of the components of the alpha tensor using both quasilinear theory and nonlinear numerical simulations. For magnetic Prandtl numbers of the order of unity, we find that in the high-conductivity limit the alpha effect in the direction of the flow clearly exceeds the component in spanwise direction. In this limit, alpha runs linearly with the magnetic Reynolds number Rm while in the low-conductivity limit it runs with the product Rm*Re, where Re is the kinetic Reynolds number so that for given Rm the alpha effect grows with decreasing magnetic Prandtl number. For the small magnetic Prandtl numbers of liquid metals, a common value for the horizontal elements of the alpha tensor appears, which makes it unimportant whether the alpha effect is measured in the spanwise or streamwise directions. The resulting effect should lead to an observable voltage in both directions of about 0.5 mV for magnetic fields of 1 kgauss and velocity fluctuations of about 1 m/s in a channel of 50 cm height (independent of its width).

preprint2013arXiv

Astrophysical hydromagnetic turbulence

Recent progress in astrophysical hydromagnetic turbulence is being reviewed. The physical ideas behind the now widely accepted Goldreich--Sridhar model and its extension to compressible magnetohydrodynamic turbulence are introduced. Implications for cosmic ray diffusion and acceleration is being discussed. Dynamo-generated magnetic fields with and without helicity are contrasted against each other. Certain turbulent transport processes are being modified and often suppressed by anisotropy and inhomogeneities of the turbulence, while others are being produced by such properties, which can lead to new large-scale instabilities of the turbulent medium. Applications of various such processes to astrophysical systems are being considered.

preprint2013arXiv

Competition of rotation and stratification in flux concentrations

In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to spontaneously form local flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism, called the negative effective magnetic pressure instability (NEMPI), is of interest in connection with dynamo scenarios where most of the magnetic field resides in the bulk of the convection zone, and not at the bottom. Recent work using the mean-field hydromagnetic equations has shown that NEMPI becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1.}{Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the effective magnetic pressure. We also quantify the kinetic helicity from direct numerical simulations (DNS) and compare with earlier work.}{To calculate the rotational effect on the effective magnetic pressure we consider both DNS and analytical studies using the $τ$ approach. To study the effects of rotation on the development of NEMPI we use both DNS and mean-field calculations of the 3D hydromagnetic equations in a Cartesian domain.}{We find that the growth rates of NEMPI from earlier mean-field calculations are well reproduced with DNS, provided the Coriolis number is below about 0.06. In that case, kinetic and magnetic helicities are found to be weak. For faster rotation, dynamo action becomes possible. However, there is an intermediate range of rotation rates where dynamo action on its own is not yet possible, but the rotational suppression of NEMPI is being alleviated.}{Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours.}

preprint2013arXiv

Microphysics of cosmic ray driven plasma instabilities

Energetic nonthermal particles (cosmic rays, CRs) are accelerated in supernova remnants, relativistic jets and other astrophysical objects. The CR energy density is typically comparable with that of the thermal components and magnetic fields. In this review we discuss mechanisms of magnetic field amplification due to instabilities induced by CRs. We derive CR kinetic and magnetohydrodynamic equations that govern cosmic plasma systems comprising the thermal background plasma, comic rays and fluctuating magnetic fields to study CR-driven instabilities. Both resonant and non-resonant instabilities are reviewed, including the Bell short-wavelength instability, and the firehose instability. Special attention is paid to the longwavelength instabilities driven by the CR current and pressure gradient. The helicity production by the CR current-driven instabilities is discussed in connection with the dynamo mechanisms of cosmic magnetic field amplification.

preprint2013arXiv

Oscillatory large-scale dynamos from Cartesian convection simulations

We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo is expected to be of $α^2\varOmega$ type which is generally expected to be relevant for the Sun, whereas the latter case refers to $α^2$ dynamos which are more likely to occur in more rapidly rotating stars whose differential rotation is small. We perform a parameter study where the shear flow and the rotational influence are varied to probe the relative importance of both types of dynamos. Oscillatory solutions are preferred both in the kinematic and saturated regimes when the negative ratio of shear to rotation rates, $q\equiv -S/\varOmega$, is between 1.5 and 2, i.e., when shear and rotation are of comparable strengths. Other regions of oscillatory solutions are found with small values of $q$, i.e., when shear is weak in comparison to rotation, and in the regime of large negative $q$s, when shear is very strong in comparison to rotation. However, exceptions to these rules also appear so that for a given ratio of shear to rotation, solutions are non-oscillatory for small and large shear, but oscillatory in the intermediate range. Changing the boundary conditions from vertical field to perfect conductor ones changes the dynamo mode from oscillatory to quasi-steady. Furthermore, in many cases an oscillatory solution exists only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. However, the cases with rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic boundary conditions. The strengths of total and large-scale components of the magnetic field in the saturated state, however, are sensitive to the chosen boundary conditions.

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

Yoshizawa's cross-helicity effect and its quenching

A central quantity in mean-field magnetohydrodynamics is the mean electromotive force EMF, which in general depends on the mean magnetic field. It may however have a part independent of the mean magnetic field. Here we study an example of a rotating conducting body of turbulent fluid with non-zero cross-helicity, in which a contribution to the EMF proportional to the angular velocity occurs (Yoshizawa 1990). If the forcing is helical, it also leads to an alpha effect, and large-scale magnetic fields can be generated. For not too rapid rotation, the field configuration is such that Yoshizawa's contribution to the EMF is considerably reduced compared to the case without alpha effect. In that case, large-scale flows are also found to be generated.

preprint2012arXiv

Effects of stratification in spherical shell convection

We report on simulations of mildly turbulent convection in spherical wedge geometry with varying density stratification. We vary the density contrast within the convection zone by a factor of 20 and study the influence of rotation on the solutions. We demonstrate that the size of convective cells decreases and the anisotropy of turbulence increases as the stratification is increased. Differential rotation is found to change from anti-solar (slow equator) to solar-like (fast equator) at roughly the same Coriolis number for all stratifications. The largest stratification runs, however, are sensitive to changes of the Reynolds number. Evidence for a near-surface shear layer is found in runs with strong stratification and large Reynolds numbers.

preprint2012arXiv

Flux concentrations in turbulent convection

We present preliminary results from high resolution magneto-convection simulations where we find the formation of flux concentrations from an initially uniform magnetic field. We compute the effective magnetic pressure but find that the concentrations appear also in places where it is positive. The structures appear in roughly ten convective turnover times and live close to a turbulent diffusion time. The time scales are compatible with the negative effective magnetic pressure instability (NEMPI), although structure formation is not restricted to regions where the effective magnetic pressure is negative.

preprint2012arXiv

Mean-field transport in stratified and/or rotating turbulence

We investigate the mean electromotive force in the kinematic framework, that is, ignoring the back-reaction of the magnetic field on the fluid velocity, under the assumption of axisymmetric turbulence determined by the presence of either rotation, density stratification, or both. We use an analogous approach for the mean passive scalar flux. As an alternative to convection, we consider forced turbulence in an isothermal layer. When using standard ansatzes, the mean magnetic transport is then determined by nine, and the mean passive scalar transport by four coefficients. We give results for all these transport coefficients. We use the test-field method and the test-scalar method, where transport coefficients are determined by solving sets of equations with properly chosen mean magnetic fields or mean scalars. These methods are adapted to mean fields which may depend on all three space coordinates. We find the anisotropy of turbulent diffusion to be moderate in spite of rapid rotation or strong density stratification. Contributions to the mean electromotive force determined by the symmetric part of the gradient tensor of the mean magnetic field, which were ignored in several earlier investigations, turn out to be important. In stratified rotating turbulence, the $α$ effect is strongly anisotropic, suppressed along the rotation axis on large length scales, but strongly enhanced at intermediate length scales. Also the $\OO\times\meanJJ$ effect is enhanced at intermediate length scales. The turbulent passive scalar diffusivity is typically almost twice as large as the turbulent magnetic diffusivity. Both magnetic and passive scalar diffusion are slightly enhanced along the rotation axis, but decreased if there is gravity.

preprint2012arXiv

Negative effective magnetic pressure in turbulent convection

We investigate the effects of weakly and strongly stratified turbulent convection on the mean effective Lorentz force, and especially on the mean effective magnetic pressure. Earlier studies with isotropically forced non-stratified and stratified turbulence have shown that the contribution of the turbulence to the mean magnetic pressure is negative for mean horizontal magnetic fields that are smaller than the equipartition strength, so that the effective mean magnetic pressure that takes into account the turbulence effects, can be negative. Compared with earlier cases of forced turbulence with an isothermal equation of state, we find that the turbulence effect is similar to or even stronger in the present case of turbulent convection. This is argued to be due to the anisotropy of turbulence in the vertical direction. Another important difference compared with earlier studies is the presence of an evolution equation for the specific entropy. Mean-field modelling with entropy evolution indicates that the negative effective magnetic pressure can still lead to a large-scale instability which forms local flux concentrations, even though the specific entropy evolution tends to have a stabilizing effect when applied to a stably stratified (e.g., isothermal) layer. It is argued that this large-scale instability could be important for the formation of solar large-scale magnetic structures such as active regions.

preprint2012arXiv

New scaling for the alpha effect in slowly rotating turbulence

Using simulations of slowly rotating stratified turbulence, we show that the alpha effect responsible for the generation of astrophysical magnetic fields is proportional to the logarithmic gradient of kinetic energy density rather than that of momentum, as was previously thought. This result is in agreement with a new analytic theory developed in this paper for large Reynolds numbers. Thus, the contribution of density stratification is less important than that of turbulent velocity. The alpha effect and other turbulent transport coefficients are determined by means of the test-field method. In addition to forced turbulence, we also investigate supernova-driven turbulence and stellar convection. In some cases (intermediate rotation rate for forced turbulence, convection with intermediate temperature stratification, and supernova-driven turbulence) we find that the contribution of density stratification might be even less important than suggested by the analytic theory.

preprint2012arXiv

Properties of the negative effective magnetic pressure instability

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large-scale instability. In this study, hydromagnetic mean-field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two-dimensional. The growth rate is found to depend on a parameter beta_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter beta_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI.

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

Transport of angular momentum and chemical species by anisotropic mixing in stellar radiative interiors

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial-to-horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non-diffusive transport known as the Lambda-effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test-scalar method. As a consequence of the non-diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (< 100 Myr) between radiative core and convective envelope in young solar-type stars is much shorter compared to the time-scale of Lithium depletion (~1 Gyr).

preprint2012arXiv

Vorticity production and survival in viscous and magnetized cosmologies

We study the role of viscosity and the effects of a magnetic field on a rotating, self-gravitating fluid, using Newtonian theory and adopting the ideal magnetohydrodynamic approximation. Our results confirm that viscosity can generate vorticity in inhomogeneous environments, while the magnetic tension can produce vorticity even in the absence of fluid pressure and density gradients. Linearizing our equations around an Einstein-de Sitter cosmology, we find that viscosity adds to the diluting effect of the universal expansion. Typically, however, the dissipative viscous effects are confined to relatively small scales. We also identify the characteristic length bellow which the viscous dissipation is strong and beyond which viscosity is essentially negligible. In contrast, magnetism seems to favor cosmic rotation. The magnetic presence is found to slow down the standard decay-rate of linear vortices, thus leading to universes with more residual rotation than generally anticipated.

preprint2011arXiv

Detection of turbulent thermal diffusion of particles in numerical simulations

The phenomenon of turbulent thermal diffusion in temperature-stratified turbulence causing a non-diffusive turbulent flux of inertial and non-inertial particles in the direction of the turbulent heat flux is found using direct numerical simulations (DNS). In simulations with and without gravity, this phenomenon is found to cause a peak in the particle number density around the minimum of the mean fluid temperature for Stokes numbers less than 1, where the Stokes number is the ratio of particle Stokes time to turbulent Kolmogorov time at the viscous scale. Turbulent thermal diffusion causes the formation of large-scale inhomogeneities in the spatial distribution of inertial particles. The strength of this effect is maximum for Stokes numbers around unity, and decreases again for larger values. The dynamics of inertial particles is studied using Lagrangian modelling in forced temperature-stratified turbulence, whereas non-inertial particles and the fluid are described using DNS in an Eulerian framework.

preprint2011arXiv

Dissipation in dynamos at low and high magnetic Prandtl numbers

Using simulations of helically driven turbulence, it is shown that the ratio of kinetic to magnetic energy dissipation scales with the magnetic Prandtl number in power law fashion with an exponent of approximately 0.6. Over six orders of magnitude in the magnetic Prandtl number the magnetic field is found to be sustained by large-scale dynamo action of alpha-squared type. This work extends a similar finding for small magnetic Prandtl numbers to the regime of large magnetic Prandtl numbers. At large magnetic Prandtl numbers, most of the energy is dissipated viscously, lowering thus the amount of magnetic energy dissipation, which means that simulations can be performed at magnetic Reynolds numbers that are large compared to the usual limits imposed by a given resolution. This is analogous to an earlier finding that at small magnetic Prandtl numbers, most of the energy is dissipated resistively, lowering the amount of kinetic energy dissipation, so simulations can then be performed at much larger fluid Reynolds numbers than otherwise. The decrease in magnetic energy dissipation at large magnetic Prandtl numbers is discussed in the context of underluminous accretion found in some quasars.

preprint2011arXiv

From convective to stellar dynamos

Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar dynamos.

preprint2011arXiv

Mean-field closure parameters for passive scalar turbulence

Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution. We vary the Reynolds (Re) and Péclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.

preprint2011arXiv

Modeling spatio-temporal nonlocality in mean-field dynamos

When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean-field equations are solved for oscillatory alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.

preprint2011arXiv

Open and closed boundaries in large-scale convective dynamos

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect-conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime increases as Rm^(1/2). Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

preprint2011arXiv

Plasma flow vs. magnetic feature-tracking speeds in the Sun

We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.

preprint2011arXiv

Pumping velocity in homogeneous helical turbulence with shear

Using different analytical methods (the quasi-linear approach, the path-integral technique and tau-relaxation approximation) we develop a comprehensive mean-field theory for a pumping effect of the mean magnetic field in homogeneous non-rotating helical turbulence with imposed large-scale shear. The effective pumping velocity is proportional to the product of alpha effect and large-scale vorticity associated with the shear, and causes a separation of the toroidal and poloidal components of the mean magnetic field along the direction of the mean vorticity. We also perform direct numerical simulations of sheared turbulence in different ranges of hydrodynamic and magnetic Reynolds numbers and use a kinematic test-field method to determine the effective pumping velocity. The results of the numerical simulations are in agreement with the theoretical predictions.

preprint2011arXiv

Reynolds stress and heat flux in spherical shell convection

(abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong `banana cells'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor--Proudman balance.

preprint2011arXiv

Verification of Reynolds stress parameterizations from simulations

We determine the timescales associated with turbulent diffusion and isotropization in closure models using anisotropically forced and freely decaying turbulence simulations and to study the applicability of these models. We compare the results from anisotropically forced three-dimensional numerical simulations with the predictions of the closure models and obtain the turbulent timescales mentioned above as functions of the Reynolds number. In a second set of simulations, turning the forcing off enables us to study the validity of the closures in freely decaying turbulence. Both types of experiments suggest that the timescale of turbulent diffusion converges to a constant value at higher Reynolds numbers. Furthermore, the relative importance of isotropization is found to be about 2.5 times larger at higher Reynolds numbers than in the more viscous regime.

preprint2010arXiv

A high-order public domain code for direct numerical simulations of turbulent combustion

A high-order scheme for direct numerical simulations of turbulent combustion is discussed. Its implementation in the massively parallel and publicly available Pencil Code is validated with the focus on hydrogen combustion. Ignition delay times (0D) and laminar flame velocities (1D) are calculated and compared with results from the commercially available Chemkin code. The scheme is verified to be fifth order in space. Upon doubling the resolution, a 32-fold increase in the accuracy of the flame front is demonstrated. Finally, also turbulent and spherical flame front velocities are calculated and the implementation of the non-reflecting so-called Navier-Stokes Characteristic Boundary Condition is validated in all three directions.

preprint2010arXiv

Convective dynamos in spherical wedge geometry

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

preprint2010arXiv

Cross helicity and turbulent magnetic diffusivity in the solar convection zone

In a density-stratified turbulent medium the cross helicity <u'.B'> is considered as a result of the interaction of the velocity fluctuations and a large-scale magnetic field. By means of a quasilinear theory and by numerical simulations we find the cross helicity and the mean vertical magnetic field anti-correlated. In the high-conductivity limit the ratio of the helicity and the mean magnetic field equals the ratio of the magnetic eddy diffusivity and the (known) density scale height. The result can be used to predict that the cross helicity at the solar surface exceeds the value of 1 Gauss km/s. Its sign is anti-correlated with that of the radial mean magnetic field. Alternatively, we can use our result to determine the value of the turbulent magnetic diffusivity from observations of the cross helicity.

preprint2010arXiv

Magnetic helicity fluxes in an alpha-squared dynamo embedded in a halo

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical alpha quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that "softer" (halo) boundary conditions give a lower energy of the saturated mean magnetic field.

preprint2010arXiv

Surface appearance of dynamo-generated large-scale fields

Twisted magnetic fields are frequently seen to emerge above the visible surface of the Sun. This emergence is usually associated with the rise of buoyant magnetic flux structures. Here we ask how magnetic fields from a turbulent large-scale dynamo appear above the surface if there is no magnetic buoyancy. The computational domain is split into two parts. In the lower part, which we refer to as the turbulence zone, the flow is driven by an assumed helical forcing function leading to dynamo action. Above this region, which we refer to as the exterior, a nearly force-free magnetic field is computed at each time step using the stress-and-relax method. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Strong current sheets tend to form above the neutral line, where the vertical field component vanishes. Time series of the magnetic field structure show recurrent plasmoid ejections. The degree to which the exterior field is force free is estimated as the ratio of the dot product of current density and magnetic field strength to their respective rms values. This ratio reaches values of up to 95% in the exterior. A weak outward flow is driven by the residual Lorentz force.

preprint2010arXiv

Test-field method for mean-field coefficients with MHD background

Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using both the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both forces are chosen to have Roberts flow-like geometry. Results: Examples with an MHD background are studied where the previously used quasi-kinematic test-field method breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the same results as the imposed-field method, where the field-aligned component of the actual electromotive force from the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, $α$ is found to be quenched proportional to the fourth power of the field strength, regardless of the type of background studied.

preprint2010arXiv

The alpha effect in rotating convection with sinusoidal shear

Using three-dimensional convection simulations it is shown that a sinusoidal variation of horizontal shear leads to a kinematic αeffect with a similar sinusoidal variation. The effect exists even for weak stratification and arises owing to the inhomogeneity of turbulence and the presence of impenetrable vertical boundaries. This system produces large-scale magnetic fields that also show a sinusoidal variation in the streamwise direction. It is argued that earlier investigations overlooked these phenomena partly because of the use of horizontal averaging and also because measurements of αusing an imposed field combined with long time averages give erroneous results. It is demonstrated that in such cases the actual horizontally averaged mean field becomes non-uniform. The turbulent magnetic diffusion term resulting from such non-uniform fields can then no longer be neglected and begins to balance the αeffect.

preprint2010arXiv

The negative magnetic pressure effect in stratified turbulence

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux tubes at the bottom of the convection zone with field strengths two orders of magnitude above equipartition, remain numerically unverified at best. As such, in recent years the idea of a formation of sunspots near the top of the convection zone has generated some interest. The presence of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective transport coefficient formalism in the mean-field formulation. The question is what happens to these coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive behaviour. In the present work we investigate the formation of large-scale magnetic structures as a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in the presence of stratification. In direct numerical simulations of forced turbulence in a stratified box, we verify the existence of this effect. This phenomenon can cause formation of large-scale magnetic structures even from initially uniform large-scale magnetic field.

preprint2010arXiv

Turbulent front speed in the Fisher equation: dependence on Damkohler number

Direct numerical simulations and mean-field theory are used to model reactive front propagation in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken into account to estimate the front speed in cases when the Damkohler number is large. This effect is found to saturate the front speed to values comparable with the speed of the turbulent motions. By comparing with direct numerical simulations, it is found that the effective correlation time is much shorter than for non-reacting flows. The nonlinearity of the reaction term is found to make the front speed slightly faster.

preprint2010arXiv

Turbulent transport in hydromagnetic flows

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

preprint2009arXiv

Large-scale dynamos in rigidly rotating turbulent convection

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.

preprint2009arXiv

Mean electromotive force proportional to mean flow in mhd turbulence

In mean-field magnetohydrodynamics the mean electromotive force due to velocity and magnetic field fluctuations plays a crucial role. In general it consists of two parts, one independent of and another one proportional to the mean magnetic field. The first part may be nonzero only in the presence of mhd turbulence, maintained, e.g., by small-scale dynamo action. It corresponds to a battery, which lets a mean magnetic field grow from zero to a finite value. The second part, which covers, e.g., the alpha effect, is important for large-scale dynamos. Only a few examples of the aforementioned first part of mean electromotive force have been discussed so far. It is shown that a mean electromotive force proportional to the mean fluid velocity, but independent of the mean magnetic field, may occur in an originally homogeneous isotropic mhd turbulence if there are nonzero correlations of velocity and electric current fluctuations or, what is equivalent, of vorticity and magnetic field fluctuations. This goes beyond the Yoshizawa effect, which consists in the occurrence of mean electromotive forces proportional to the mean vorticity or to the angular velocity defining the Coriolis force in a rotating frame and depends on the cross-helicity defined by the velocity and magnetic field fluctuations. Contributions to the mean electromotive force due to inhomogeneity of the turbulence are also considered. Possible consequences of the above and related findings for the generation of magnetic fields in cosmic bodies are discussed.

preprint2009arXiv

Turbulent dynamos with shear and fractional helicity

Dynamo action owing to helically forced turbulence and large-scale shear is studied using direct numerical simulations. The resulting magnetic field displays propagating wave-like behavior. This behavior can be modelled in terms of an αΩdynamo. In most cases super-equipartition fields are generated. By varying the fraction of helicity of the turbulence the regeneration of poloidal fields via the helicity effect (corresponding to the α-effect) is regulated. The saturation level of the magnetic field in the numerical models is consistent with a linear dependence on the ratio of the fractional helicities of the small and large-scale fields, as predicted by a simple nonlinear mean-field model. As the magnetic Reynolds number (Rm) based on the wavenumber of the energy-carrying eddies is increased from 1 to 180, the cycle frequency of the large-scale field is found to decrease by a factor of about 6 in cases where the turbulence is fully helical. This is interpreted in terms of the turbulent magnetic diffusivity, which is found to be only weakly dependent on Rm.

preprint2008arXiv

Magnetic diffusivity tensor and dynamo effects in rotating and shearing turbulence

The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and opposite to the sign of the shear. This implies that dynamo action from the shear--current effect is impossible, except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the average, the components of the alpha tensor display Gaussian fluctuations around zero. These fluctuations are strong enough to drive an incoherent alpha--shear dynamo. The incoherent shear--current effect, on the other hand, is found to be subdominant.

preprint2007arXiv

The helicity constraint in spherical shell dynamos

The motivation for considering distributed large scale dynamos in the solar context is reviewed in connection with the magnetic helicity constraint. Preliminary accounts of 3-dimensional direct numerical simulations (in spherical shell segments) and simulations of 2-dimensional mean field models (in spherical shells) are presented. Interesting similarities as well as some differences are noted.

preprint2007arXiv

Turbulent viscosity and Lambda-effect from numerical turbulence models

Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress tensor and its parameterization in terms of turbulent viscosity and Lambda-effect. The turbulence is forced in an anisotropic fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here.

preprint2006arXiv

The nature of turbulence in OMC1 at the star forming scale: observations and simulations

Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions. Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence. Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.

preprint2004arXiv

Top quark pair production and decay at hadron colliders

In ongoing and upcoming hadron collider experiments, top quark physics will play an important role in testing the Standard Model and its possible extensions. In this work we present analytic results for the differential cross sections of top quark pair production in hadronic collisions at next-to-leading order in the QCD coupling, keeping the full dependence on the spins of the top quarks. These results are combined with the corresponding next-to-leading order results for the decay of polarized top quarks into dilepton, lepton plus jets, and all jets final states. As an application we predict double differential angular distributions which are due to the QCD-induced top quark spin correlations in the intermediate state. In addition to the analytic results, we give numerical results in terms of fit functions that can easily be used in an experimental analysis.

preprint2002arXiv

QCD-corrected spin analysing power of jets in decays of polarized top quarks

We present results for the differential distributions of jets from non-leptonic decays of polarized top quarks within the Standard Model, including QCD radiative corrections. Our work extends existing results which are only available for semileptonic top quark decays at the parton level. For $t(\uparrow)$ -> b-jet + 2 light jets we compute in particular the QCD-corrected top-spin analysing power of the b-quark jet and the least energetic light jet. The dependence of the results on the choice of the jet recombination scheme is found to be small. In addition we compute the spin analysing power of the thrust axis. Our results constitute a so far missing ingredient to analyse top quark production and subsequent non-leptonic decay at next-to-leading order in alpha_s, keeping the full information on the top quark polarization.

preprint1998arXiv

The spin density matrix of top quark pairs produced in electron-positron annihilation including QCD radiative corrections

We calculate the spin density matrix of top quark pairs for the reaction $e^+e^-\to t\bar{t}X$ to order $α_s$. As an application we show next-to-leading order results for a variety of spin observables for the $t\bar{t}$ system. These include the top quark and antiquark polarizations and $t\bar{t}$ spin-spin correlations as a function of the center-of-mass energy and of the top quark scattering angle for arbitrary longitudinal polarization of the electron/positron beam.

preprint1997arXiv

QCD Corrections to Decay Distributions of Neutral Higgs Bosons with (In)definite CP Parity

We compute the order $α_s$ QCD corrections to the density matrix for the decay of a neutral Higgs boson $ϕ$ with (in)definite CP parity into a quark antiquark pair, respectively the QED corrections for the decay into a pair of charged leptons. We classify and calculate single spin asymmetries and spin-spin correlations which are generated by the scalar and pseudoscalar Yukawa couplings. These spin effects can be traced in $ϕ\to τ^-τ^+$ and, for heavy Higgs bosons, in $ϕ\to t\bar{t}$. We also calculate resulting correlations among the final states and estimate, for the respective decay modes, the number of events needed to measure the Yukawa couplings with these correlations at the 3 sigma level.

preprint1995arXiv

Extraction of the pion distribution amplitude from polarized muon pair production

We consider the production of muon pairs from the scattering of pions on longitudinally polarized protons. We calculate the cross section and the single spin asymmetry for this process, taking into account pion bound state effects. We work in the kinematic region where the photon has a large longitudinal momentum fraction, which allows us to treat the bound state problem perturbatively. Our predictions are directly proportional to the pion distribution amplitude. A measurement of the polarized Drell-Yan cross section thus allows the determination of the shape of the pion distribution amplitude.

preprint1995arXiv

Probing Higgs Sector CP Violation with Top Quarks at a Photon Linear Collider

We investigate for two-Higgs doublet models with intrinsic CP violation in the scalar potential CP-nonconserving effects in unpolarized $γγ\rightarrow t \bar{t} $ for a range of neutral Higgs boson masses which includes resonant $φ$ production and the subsequent decay of $φ\rightarrow t \bar{t}$. The importance of taking into account, even in the resonant case, the interference with the nonresonant background is shown. Further, we propose and calculate three asymmetries which efficiently trace CP-violating effects in $γγ\rightarrow t \bar{t} $ using semileptonic $t \bar{t}$ decays.

preprint1993arXiv

Tracing CP violation in the production of top quark pairs by multiple TeV proton-proton collisions

We investigate the possibilities of searching for non-standard CP violation in $pp\to t\bar{t}X$ at multiple TeV collision energies. A general kinematic analysis of the underlying partonic production processes $gg\to t\bar{t}$ and $q\bar{q}\to t\bar{t}$ in terms of their density matrices is given. We evaluate the CP-violating parts of these matrices in two-Higgs doublet extensions of the standard model (SM) and give results for CP asymmetries at the parton level. We show that these asymmetries can be traced by measuring suitable observables constructed from energies and momenta of the decay products of $t$ and $\bar{t} $. We find CP-violating effects to be of the order of $10^{-3}$ and show that possible contaminations induced by SM interactions are savely below the expected signals.