BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

P. J. Käpylä

P. J. Käpylä appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-phphysics.flu-dynphysics.plasm-phastro-ph.COastro-ph.GAastro-ph.IM

Catalog footprint

What is connected

30works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Magnetism, rotation, and nonthermal emission in cool stars -- Average magnetic field measurements in 292 M dwarfs

Stellar dynamos generate magnetic fields that are of fundamental importance to the variability and evolution of Sun-like and low-mass stars, and for the development of their planetary systems. We report measurements of surface-average magnetic fields in 292 M dwarfs from a comparison with radiative transfer calculations; for 260 of them, this is the first measurement of this kind. Our data were obtained from more than 15,000 high-resolution spectra taken during the CARMENES project. They reveal a relation between average field strength, <B>, and Rossby number, $Ro$, resembling the well-studied rotation-activity relation. Among the slowly rotating stars, we find that magnetic flux, $Φ_\textrm{B}$, is proportional to rotation period, $P$, and among the rapidly rotating stars that average surface fields do not grow significantly beyond the level set by the available kinetic energy. Furthermore, we find close relations between nonthermal coronal X-ray emission, chromospheric H$α$ and Ca H&K emission, and magnetic flux. Taken together, these relations demonstrate empirically that the rotation-activity relation can be traced back to a dependence of the magnetic dynamo on rotation. We advocate the picture that the magnetic dynamo generates magnetic flux on the stellar surface proportional to rotation rate with a saturation limit set by the available kinetic energy, and we provide relations for average field strengths and nonthermal emission that are independent of the choice of the convective turnover time. We also find that Ca H&K emission saturates at average field strengths of $\langle B \rangle \approx 800$ G while H$α$ and X-ray emission grow further with stronger fields in the more rapidly rotating stars. This is in conflict with the coronal stripping scenario predicting that in the most rapidly rotating stars coronal plasma would be cooled to chromospheric temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electrodynamics of turbulent fluids with fluctuating electric conductivity

The influence of fluctuating conductivity on the coefficients known from the mean-field electrodynamics is considered. If the conductivity fluctuations are assumed as uncorrelated with the turbulent velocity field then only the effective magnetic diffusivity of the fluid is reduced and the decay time of a large-scale magnetic field is increased. If the fluctuations of conductivity and flow are correlated in a certain direction then an additional diamagnetic pumping effect results transporting magnetic field in opposite direction to the resistivity flux vector $\langle η'\vec{u}'\rangle$. Even for homogeneous turbulence fields in the presence of rotation an alpha effect appears. With the characteristic values of the outer core of the Earth or the solar convection zone, however, the dynamo number of the alpha effect never reaches supercritical values to operate as an $α^2$-dynamo.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing turbulent closure models with convection simulations

We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors, corresponding to locating the computational domain at different latitudes ranging from the pole to the equator of the star. We perform three-dimensional numerical simulations in the same parameter regimes for comparison. The free parameters appearing in the closure models are calibrated by two fitting methods using simulation data. Unique determination of the closure parameters is possible only in the non-rotating case or when the system is placed at the pole. In the other cases the fit procedures yield somewhat differing results. The quality of the closure is tested by substituting the resulting coefficients back into the closure model and comparing with the simulation results. To eliminate the possibilities that the results obtained depend on the aspect ratio of the simulation domain or suffer from too small Rayleigh numbers we performed runs varying these parameters. The simulation data for the Reynolds stress and heat fluxes broadly agree with previous compressible simulations. The closure works fairly well with slow and fast rotation but its quality degrades for intermediate rotation rates. We find that the closure parameters depend not only on rotation rate but also on latitude. The weak dependence on Rayleigh number and the aspect ratio of the domain indicates that our results are generally valid

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Starspots due to large-scale vortices in rotating turbulent convection

We study the generation of large-scale vortices in rotating turbulent convection by means of Cartesian direct numerical simulations. We find that for sufficiently rapid rotation, cyclonic structures on a scale large in comparison to that of the convective eddies, emerge, provided that the fluid Reynolds number exceeds a critical value. For slower rotation, cool cyclonic vortices are preferred, whereas for rapid rotation, warm anti-cyclonic vortices are favored. In some runs in the intermediate regime both types of cyclones co-exist for thousands of convective turnover times. The temperature contrast between the vortices and the surrounding atmosphere is of the order of five per cent. We relate the simulation results to observations of rapidly rotating late-type stars that are known to exhibit large high-latitude spots from Doppler imaging. In many cases, cool spots are accompanied with spotted regions with temperatures higher than the average. In this paper, we investigate a scenario according to which the spots observed in the temperature maps could have a non-magnetic origin due to large-scale vortices in the convection zones of the stars.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Effects of Rotation and Input Energy Flux on Convective Overshooting

We study convective overshooting by means of local 3D convection calculations. Using a mixing length model of the solar convection zone (CZ) as a guide, we determine the Coriolis number (Co), which is the inverse of the Rossby number, to be of the order of ten or larger at the base of the solar CZ. Therefore we perform convection calculations in the range Co = 0...10 and interpret the value of Co realised in the calculation to represent a depth in the solar CZ. In order to study the dependence on rotation, we compute the mixing length parameters alpha_T and alpha_u relating the temperature and velocity fluctuations, respectively, to the mean thermal stratification. We find that the mixing length parameters for the rapid rotation case, corresponding to the base of the solar CZ, are 3-5 times smaller than in the nonrotating case. Introducing such depth-dependent alpha into a solar structure model employing a non-local mixing length formalism results in overshooting which is approximately proportional to alpha at the base of the CZ. Although overshooting is reduced due to the reduced alpha, a discrepancy with helioseismology remains due to the steep transition to the radiative temperature gradient. In comparison to the mixing length models the transition at the base of the CZ is much gentler in the 3D models. It was suggested recently (Rempel 2004) that this discrepancy is due to the significantly larger (up to seven orders of magnitude) input energy flux in the 3D models in comparison to the Sun and solar models, and that the 3D calculations should be able to approach the mixing length regime if the input energy flux is decreased by a moderate amount. We present results from local convection calculations which support this conjecture.

0 citations0 reviewsNo code linkedOpen access