What sets the magnetic field strength and cycle period in solar-type stars?

Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two regimes of activity, the active and inactive branches, also have been identified. For both of them, the longer the rotation period, the longer the activity cycle. Using global dynamo simulations of solar like stars with Rossby numbers between ~0.4 and ~2, this paper explores the relevance of rotational shear layers in determining these observational properties. Our results, consistent with non-linear alpha^2-Omega dynamos, show that the total magnetic field strength is independent of the rotation period. Yet at surface levels, the origin of the magnetic field is determined by Ro. While for Ro<1 it is generated in the convection zone, for Ro>1 strong toroidal fields are generated at the tachocline and rapidly emerge towards the surface. In agreement with the observations, the magnetic cycle period increases with the rotational period. However, a bifurcation is observed for Ro~1, separating a regime where oscillatory dynamos operate mainly in the convection zone, from the regime where the tachocline has a predominant role. In the latter the cycles are believed to result from the periodic energy exchange between the dynamo and the magneto-shear instabilities developing in the tachocline and the radiative interior.Comment: 43 pages, 14 figures, accepted for publication in The Astrophysical Journa

Multi-scale waves in sound-proof global simulations with EULAG

EULAG is a computational model for simulating flows across a wide range of scales and physical scenarios. A standard option employs an anelastic approximation to capture nonhydrostatic effects and simultaneously filter sound waves from the solution. In this study, we examine a localized gravity wave packet generated by instabilities in Held-Suarez climates. Although still simplified versus the Earth’s atmosphere, a rich set of planetary wave instabilities and ensuing radiated gravity waves can arise. Wave packets are observed that have lifetimes ≤ 2 days, are negligibly impacted by Coriolis force, and do not show the rotational effects of differential jet advection typical of inertia-gravity waves. Linear modal analysis shows that wavelength, period, and phase speed fit the dispersion equation to within a mean difference of ∼ 4%, suggesting an excellent fit. However, the group velocities match poorly even though a propagation of uncertainty analysis indicates that they should be predicted as well as the phase velocities. Theoretical arguments suggest the discrepancy is due to nonlinearity — a strong southerly flow leads to a critical surface forming to the southwest of the wave packet that prevents the expected propagation

3D MHD Modeling of the Impact of Subsurface Stratification on the Solar Dynamo

Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram - from a North-South symmetric pattern to a staggered anti-symmetric emergence

A Comparison of Two Shallow Water Models with Non-Conforming Adaptive Grids: classical tests

In an effort to study the applicability of adaptive mesh refinement (AMR) techniques to atmospheric models an interpolation-based spectral element shallow water model on a cubed-sphere grid is compared to a block-structured finite volume method in latitude-longitude geometry. Both models utilize a non-conforming adaptation approach which doubles the resolution at fine-coarse mesh interfaces. The underlying AMR libraries are quad-tree based and ensure that neighboring regions can only differ by one refinement level. The models are compared via selected test cases from a standard test suite for the shallow water equations. They include the advection of a cosine bell, a steady-state geostrophic flow, a flow over an idealized mountain and a Rossby-Haurwitz wave. Both static and dynamics adaptations are evaluated which reveal the strengths and weaknesses of the AMR techniques. Overall, the AMR simulations show that both models successfully place static and dynamic adaptations in local regions without requiring a fine grid in the global domain. The adaptive grids reliably track features of interests without visible distortions or noise at mesh interfaces. Simple threshold adaptation criteria for the geopotential height and the relative vorticity are assessed.Comment: 25 pages, 11 figures, preprin

What sets the magnetic field strength and cycle period in solar-type stars?

Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two regimes of activity, the active and inactive branches, also have been identified. For both of them, the longer the rotation period, the longer the activity cycle. Using global dynamo simulations of solar like stars with Rossby numbers between ~0.4 and ~2, this paper explores the relevance of rotational shear layers in determining these observational properties. Our results, consistent with non-linear alpha^2-Omega dynamos, show that the total magnetic field strength is independent of the rotation period. Yet at surface levels, the origin of the magnetic field is determined by Ro. While for Ro<1 it is generated in the convection zone, for Ro>1 strong toroidal fields are generated at the tachocline and rapidly emerge towards the surface. In agreement with the observations, the magnetic cycle period increases with the rotational period. However, a bifurcation is observed for Ro~1, separating a regime where oscillatory dynamos operate mainly in the convection zone, from the regime where the tachocline has a predominant role. In the latter the cycles are believed to result from the periodic energy exchange between the dynamo and the magneto-shear instabilities developing in the tachocline and the radiative interior.Comment: 43 pages, 14 figures, accepted for publication in The Astrophysical Journa

Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows

We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable

First analysis of anisotropic flow with Lee--Yang zeroes

We report on the first analysis of directed and elliptic flow with the new method of Lee--Yang zeroes. Experimental data are presented for Ru+Ru reactions at 1.69 AGeV measured with the FOPI detector at SIS/GSI. The results obtained with several methods, based on the event-plane reconstruction, on Lee--Yang zeroes, and on multi-particle cumulants (up to 5th order) applied for the first time at SIS energies, are compared. They show conclusive evidence that azimuthal correlations between nucleons and composite particles at this energy are largely dominated by anisotropic flow.Comment: 5 pages, 3 figures, submitted to Phys. Rev. C Rapid Co

Meridional circulation dynamics in a cyclic convective dynamo

Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to slowing down the axisymmetric poleward horizontal component. The behavior of this large-scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large-scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid-latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.Centra-ISTGRPS-UdeMNatural Sciences and Engineering Research Council of CanadaNational Science FoundationUniversity of the Algarveinfo:eu-repo/semantics/publishedVersio