BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Petri J. Käpylä

Petri J. Käpylä contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.HEphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2022arXiv

Origin of eclipsing time variations: Contributions of different modes of the dynamo-generated magnetic field

The possibility to detect circumbinary planets and to study stellar magnetic fields through eclipsing time variations (ETVs) in binary stars has sparked an increase of interest in this area of research. We revisit the connection between stellar magnetic fields and the gravitational quadrupole moment $Q_{xx}$ and compare different dynamo-generated ETV models with our simulations. We present magnetohydrodynamical simulations of solar mass stars with rotation periods of 8.3, 1.2, and 0.8 days and perform a detailed analysis of the magnetic and quadrupole moment using spherical harmonic decomposition. The extrema of $Q_{xx}$ are associated with changes in the magnetic field structure. This is evident in the simulation with a rotation period of 1.2 days. Its magnetic field has a more complex behavior than in the other models, as the large-scale nonaxisymmetric field dominates throughout the simulation and the axisymmetric component is predominantly hemispheric. This triggers variations in the density field that follow the magnetic field asymmetry with respect to the equator, affecting the $zz$ component of the inertia tensor, and thus modulating $Q_{xx}$. The magnetic field of the two other runs are less variable in time and more symmetric with respect to the equator, such that the variations in the density are weaker, and therefore only small variations in $Q_{xx}$ are seen. If interpreted via the classical Applegate mechanism (tidal locking), the quadrupole moment variations obtained in the current simulations are about two orders of magnitude below those deduced from observations of post-common-envelope binaries. However, if no tidal locking is assumed, our results are compatible with the observed ETVs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulations of fully convective M dwarfs: dynamo action with varying magnetic Prandtl numbers

M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an important role in the dynamo action; dynamos for low-$\rm Pr_M$ and large-$\rm Pr_M$ have very different properties. We performed three-dimensional magnetohydrodynamical (MHD) numerical simulations with the ``star-in-a-box'' model using stellar parameters for an M5 dwarf with $0.21 M_{\odot}$. We found that the dynamo solutions are sensitive to ${\rm Pr_M}$. The simulations at this rotation period present periodic cycles of the large-scale magnetic field up to ${\rm Pr_M} \leq 2$; for higher values the cycles disappear and irregular solutions to arise. Our results are consistent with previous studies and suggest that the dynamos operating in fully convective stars behave similarly as those in partially convective stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solar-like dynamos and rotational scaling of cycles from star-in-a-box simulations

Magnetohydrodynamic star-in-a-box simulations of convection and dynamos in a solar-like star with different rotation rates are presented. These simulations produce solar-like differential rotation with a fast equator and slow poles, and magnetic activity that resembles that of the Sun with equatorward migrating activity at the surface. Furthermore, the ratio of rotation to cycle period is almost constant as the rotation period decreases in the limited sample considered here. This is reminiscent of the suggested inactive branch of stars from observations and differs from most earlier simulation results from spherical shell models. While the exact excitation mechanism of the dynamos in the current simulations is not yet clear, it is plausible that the greater freedom that the magnetic field has due to the inclusion of the radiative core and regions exterior to the star are important in shaping the dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transition from anti-solar to solar-like differential rotation: Dependence on Prandtl number

(abridged) Context: Late-type stars such as the Sun rotate differentially due to the interaction of turbulent convection and rotation. Aims: The aim of the study is to investigate the effects of the thermal Prandtl number on the transition from anti-solar (slow equator, fast poles) to solar-like (fast equator, slow poles) differential rotation. Methods: Three-dimensional hydrodynamic and magnetohydrodynamic simulations in semi-global spherical wedge geometry are used to model convection zones of solar-like stars. Results: The overall convective velocity amplitude increases as the Prandtl number decreases in accordance with earlier studies. The transition from anti-solar to solar-like differential rotation is insensitive to the Prandtl number for Prandtl numbers below unity but for Prandtl numbers greater than unity, solar-like differential rotation becomes significantly harder to excite. Magnetic fields and more turbulent regimes with higher fluid and magnetic Reynolds numbers help in achieving solar-like differential rotation in near-transition cases. Solar-like differential rotation occurs only in cases with radially outward angular momentum transport at the equator. The dominant contribution to such outward transport near the equator is due to prograde propagating thermal Rossby waves. Conclusions: The differential rotation is sensitive to the Prandtl number only for large Prandtl numbers in the parameter regime explored in the current study. Magnetic fields have a greater effect on the differential rotation, although the inferred presence of a small-scale dynamo does not lead to drastically different results in the present study. The dominance of the thermal Rossby waves in the simulations is puzzling given the non-detection in the Sun. The current simulations are shown to be incompatible with the currently prevailing mean-field theory of differential rotation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Common dynamo scaling in slowly rotating young and evolved stars

One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence to stellar dynamos by including evolved, post main sequence stars in the analysis of the rotation-activity relation. These stars rotate very slowly compared with main sequence stars, but exhibit similar activity levels. We show that the two evolutionary stages fall together in the rotation-activity diagram and form a single sequence in the unsaturated regime in relation only to Rossby numbers derived from stellar models, confirming earlier preliminary results that relied on a more simplistic parametrisation of the convective turnover time. This mirrors recent results of fully convective M dwarfs, which likewise fall on the same rotation-activity sequence as partially convective solar-type stars. Our results demonstrate that turbulence plays a crucial role in driving stellar dynamos and suggest that there is a common turbulence-related dynamo mechanism explaining the magnetic activity of all late-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of small-scale dynamo and compressibility on the $Λ$ effect

The $Λ$ effect describes a rotation-induced non-diffusive contribution to the Reynolds stress. It is commonly held responsible for maintaining the observed differential rotation of the Sun and other late-type stars. Here the sensitivity of the $Λ$ effect to small-scale magnetic fields and compressibility is studied by means of forced turbulence simulations either with anisotropic forcing in fully periodic cubes or in density-stratified domains with isotropic forcing. Effects of small-scale magnetic fields are studied in cases where the magnetic fields are self-consistently generated by a small-scale dynamo. The results show that small-scale magnetic fields lead to a quenching of the $Λ$ effect which is milder than in cases where also a large-scale field is present. The effect of compressibility on the $Λ$ effect is negligible in the range of Mach numbers from 0.015 to 0.8. Density stratification induces a marked anisotropy in the turbulence and a vertical $Λ$ effect if the forcing scale is roughly two times larger than the density scale height.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Overshooting in simulations of compressible convection

(abridged) Context: Convective motions overshooting to regions that are formally convectively stable cause extended mixing. Aims: To determine the scaling of overshooting depth ($d_{\rm os}$) at the base of the convection zone as a function of imposed energy flux ($\mathscr{F}_{\rm n}$) and to estimate the extent of overshooting at the base of the solar convection zone. Methods: Three-dimensional Cartesian simulations of compressible non-rotating convection with unstable and stable layers are used. The simulations use either a fixed heat conduction profile or a temperature and density dependent formulation based on Kramers opacity law. The simulations cover a range of almost four orders of magnitude in the imposed flux. Results: A smooth heat conduction profile (either fixed or through Kramers opacity law) leads to a relatively shallow power law with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.08}$ for low $\mathscr{F}_{\rm n}$. A fixed step-profile of the heat conductivity at the bottom of the convection zone leads to a somewhat steeper dependency with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.12}$. Experiments with and without subgrid-scale entropy diffusion revealed a strong dependence on the effective Prandtl number which is likely to explain the steep power laws as a function of $\mathscr{F}_{\rm n}$ reported in the literature. Furthermore, changing the heat conductivity artificially below the convection zone is shown to lead to substantial underestimation of overshooting depth. Conclusions: Extrapolating from the results obtained with smooth heat conductivity profiles suggest that the overshooting depth for the solar flux is of the order of $0.2$ pressure scale heights at the base of the convection zone which is two to four times higher than estimates from helioseismology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence

(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magneto-hydrodynamical origin of eclipsing time variations in post-common-envelope binaries for solar mass secondaries

Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cyclic magnetic activity due to turbulent convection in spherical wedge geometry

We report on simulations of turbulent, rotating, stratified, magnetohydrodynamic convection in spherical wedge geometry. An initially small-scale, random, weak-amplitude magnetic field is amplified by several orders of magnitude in the course of the simulation to form oscillatory large-scale fields in the saturated state of the dynamo. The differential rotation is solar-like (fast equator), but neither coherent meridional poleward circulation nor near-surface shear layer develop in these runs. In addition to a poleward branch of magnetic activity beyond 50 degrees latitude, we find for the first time a pronounced equatorward branch at around 20 degrees latitude, reminiscent of the solar cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of enhanced stratification on equatorward dynamo wave propagation

We present results from simulations of rotating magnetized turbulent convection in spherical wedge geometry representing parts of the latitudinal and longitudinal extents of a star. Here we consider a set of runs for which the density stratification is varied, keeping the Reynolds and Coriolis numbers at similar values. In the case of weak stratification, we find quasi-steady dynamo solutions for moderate rotation and oscillatory ones with poleward migration of activity belts for more rapid rotation. For stronger stratification, the growth rate tends to become smaller. Furthermore, a transition from quasi-steady to oscillatory dynamos is found as the Coriolis number is increased, but now there is an equatorward migrating branch near the equator. The breakpoint where this happens corresponds to a rotation rate that is about 3-7 times the solar value. The phase relation of the magnetic field is such that the toroidal field lags behind the radial field by about $π/2$, which can be explained by an oscillatory $α^2$ dynamo caused by the sign change of the $α$-effect about the equator. We test the domain size dependence of our results for a rapidly rotating run with equatorward migration by varying the longitudinal extent of our wedge. The energy of the axisymmetric mean magnetic field decreases as the domain size increases and we find that an $m=1$ mode is excited for a full $2π$ azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spoke-like differential rotation in a convective dynamo with a coronal envelope

We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ times the value at the bottom of the spherical shell. We associate this region with the solar corona. We find solar-like differential rotation with radial contours of constant rotation rate, together with a near-surface shear layer. This non-cylindrical rotation profile is caused by a non-zero latitudinal entropy gradient that offsets the Taylor--Proudman balance through the baroclinic term. The meridional circulation is multi-cellular with a solar-like poleward flow near the surface at low latitudes. In most of the cases, the mean magnetic field is oscillatory with equatorward migration in two cases. In other cases, the equatorward migration is overlaid by stationary or even poleward migrating mean fields.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ejections of magnetic structures above a spherical wedge driven by a convective dynamo with differential rotation

We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists of a convection zone and an extended coronal region to $\approx1.5$ times the radius of the sphere. The wedge contains a quarter of the azimuthal extent of the sphere and $150\degr$ in latitude. The magnetic field is self-consistently generated by the turbulent motions due to convection beneath the surface. Magnetic fields are found to emerge at the surface and are ejected to the coronal part of the domain. These ejections occur at irregular intervals and are weaker than in earlier work. We tentatively associate these events with coronal mass ejections on the Sun, even though our model of the solar atmosphere is rather simplistic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Angular momentum transport in convectively unstable shear flows

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetorotational instability driven dynamos at low magnetic Prandtl numbers

Numerical simulations of the magnetorotational instability (MRI) with zero initial net flux in a non-stratified isothermal cubic domain are used to demonstrate the importance of magnetic boundary conditions.In fully periodic systems the level of turbulence generated by the MRI strongly decreases as the magnetic Prandtl number (Pm), which is the ratio of kinematic viscosity and magnetic diffusion, is decreased. No MRI or dynamo action below Pm=1 is found, agreeing with earlier investigations. Using vertical field conditions, which allow the generation of a net toroidal flux and magnetic helicity fluxes out of the system, the MRI is found to be excited in the range 0.1 < Pm < 10, and that the saturation level is independent of Pm. In the vertical field runs strong mean-field dynamo develops and helps to sustain the MRI.

0 citations0 reviewsNo code linkedOpen access