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.

Various models of solar sub-surface 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.

Unveiling the evolution of toroidal field instability, known as Tayler instability, is essential to understand the strength and topology of the magnetic fields observed in early-type stars, in the core of the red giants, or in any stellar radiative zone. We want to study the non-linear evolution of the instability of a toroidal field stored in a stably stratified layer, in spherical symmetry and in the absence of rotation. In particular, we intend to quantify the suppression of the instability as a function of the Brunt-V\"ais\"ala ($\omega_{\rm BV}$) and the Alfv\'en ($\omega_{\rm A}$) frequencies. We use the MHD equations as implemented in the anelastic approximation in the EULAG-MHD code and perform a large series of numerical simulations of the instability exploring the parameter space for the $\omega_{\rm BV}$ and $\omega_{\rm A}$. We show that beyond a critical value gravity strongly suppress the instability, in agreement with the linear analysis. The intensity of the initial field also plays an important role: weaker fields show much slower growth rates. Moreover, in the case of very low gravity, the fastest growing modes have a large characteristic radial scale, at variance with the case of strong gravity, where the instability is characterized by horizontal displacements. Our results illustrate that the anelastic approximation can efficiently describe the evolution of toroidal field instability in stellar interiors. The suppression of the instability as a consequence of increasing values of $\omega_{\rm BV}$ might play a role to explain the magnetic desert in Ap/Bp stars since weak fields are only marginally unstable in the case of strong gravity.

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.

Recent observations of the magnetic field in pre-main sequence stars suggest that the magnetic field topology changes as a function of age. The presence of a tachocline could be an important factor in the development of magnetic field with higher multipolar modes. In this work we performed MHD simulations using the EULAG-MHD code to study the magnetic field generation and evolution in models that mimic stars at two evolutionary stages. The stratification for both stellar phases was computed by fitting stellar structure profiles obtained with the ATON stellar evolution code. The first stage is at $1.1 \mathrm{Myr}$, when the star is completely convective. The second stage is at $14 \mathrm{Myrs}$, when the star is partly convective, with a radiative core developed up to $30\%$ of the stellar radius. In this proceedings we present a preliminary analysis of the resulting mean-flows and magnetic field. The mean-flow analysis shown that the star rotate almost rigidly on the fully convective phase, whereas at the partially convective phase there is differential rotation with conical contours of iso-rotation. As for the mean magnetic field both simulations show similarities with respect to the field evolution. However, the topology of the magnetic field is different.

Parker's magnetostatic theorem extended to astrophysical magnetofluids with large magnetic Reynolds number supports ceaseless regeneration of current sheets and hence, spontaneous magnetic reconnections recurring in time. Consequently, a scenario is possible where the repeated reconnections provide an autonomous mechanism governing emergence of coherent structures in astrophysical magnetofluids. In this work, such a scenario is explored by performing numerical computations commensurate with the magnetostatic theorem. In particular, the computations explore the evolution of a flux-rope governed by repeated reconnections in a magnetic geometry resembling bipolar loops of solar corona. The revealed morphology of the evolution process, including onset and ascent of the rope, reconnection locations and the associated topology of the magnetic field lines, agrees with observations, and thus substantiates physical realisability of the advocated mechanism.

The phenomenon of solar "torsional oscillations" (TO) represents migratory zonal flows associated with the solar cycle. These flows are observed on the solar surface and, according to helioseismology, extend through the convection zone. We study the origin of the TO using results from a global MHD simulation of the solar interior that reproduces several of the observed characteristics of the mean-flows and magnetic fields. Our results indicate that the magnetic tension (MT) in the tachocline region is a key factor for the periodic changes in the angular momentum transport that causes the TO. The torque induced by the MT at the base of the convection zone is positive at the poles and negative at the equator. A rising MT torque at higher latitudes causes the poles to speed-up, whereas a declining negative MT torque at the lower latitudes causes the equator to slow-down. These changes in the zonal flows propagate through the convection zone up to the surface. Additionally, our results suggest that it is the magnetic field at the tachocline that modulates the amplitude of the surface meridional flow rather than the opposite as assumed by flux-transport dynamo models of the solar cycle.

The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were used to obtain these results. It appears today mandatory to assess their robustness with respect to the details of the numerics, and in particular to the treatment of turbulent sub-grid scales. We report on an ongoing comparison between two global models, the ASH and EULAG codes. In EULAG the sub-grid scales are treated implicitly by the numerical scheme, while in ASH their effect is generally modelled by using enhanced dissipation coefficients. We characterize the sub-grid scales effect in a turbulent convection simulation with EULAG. We assess their effect at each resolved scale with a detailed energy budget. We derive equivalent eddy-diffusion coefficients and use the derived diffusivities in twin ASH numerical simulations. We find a good agreement between the large-scale flows developing in the two codes in the hydrodynamic regime, which encourages further investigation in the magnetohydrodynamic regime for various dynamo solutions.

Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that the resulting mean-flows and dynamo properties like the growth rate, saturation energy and mode depend on the Rossby (Ro) number. For the first set of models either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field. Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo coefficients computed using FOSA, the field evolution in the oscillatory models without a tachocline seems to be consistent with dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachoclines the dynamics is more complex involving other transport mechanisms as well as tachocline instabilities.

To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal flows accelerated either in the direction of rotation (solar-like) or in the opposite direction (anti-solar). While the models with the solar-like differential rotation tend to produce multiple cells of meridional circulation, the models with anti-solar differential rotation result in only one or two meridional cells. Our simulations indicate that the rotation and large-scale flow patterns critically depend on the ratio between buoyancy and Coriolis forces. By including a subadiabatic layer at the bottom of the domain, corresponding to the stratification of a radiative zone, we reproduce a layer of strong radial shear similar to the solar tachocline. Similarly, enhanced superadiabaticity at the top results in a near-surface shear layer located mainly at lower latitudes. The models reveal a latitudinal entropy gradient localized at the base of the convection zone and in the stable region, which however does not propagate across the convection zone. In consequence, baroclinicity effects remain small and the rotation iso-contours align in cylinders along the rotation axis. Our results confirm the alignment of large convective cells along the rotation axis in the deep convection zone, and suggest that such "banana-cell" pattern can be hidden beneath the supergranulation layer.

Convective turbulent motions in the solar interior, as well as the mean flows resulting from them, determine the evolution of the solar magnetic field. With the aim to get a better understanding of these flows we study anelastic rotating convection in a spherical shell whose stratification resembles that of the solar interior. This study is done through numerical simulations performed with the EULAG code. Due to the numerical formulation, these simulations are known as implicit large eddy simulations (ILES), since they intrinsically capture the contribution of, non-resolved, small scales at the same time maximizing the effective Reynolds number. We reproduce some previous results and find a transition between buoyancy and rotation dominated regimes which results in anti-solar or solar like rotation patterns. Even thought the rotation profiles are dominated by Taylor-Proudman columnar rotation, we are able to reproduce the tachocline and a low latitude near-surface shear layer. We find that simulations results depend on the grid resolution as a consequence of a different sub-grid scale contribution.

We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett., 715, L133 (2010); Racine et al., Astrophys. J., 735, 46 (2011)) producing an axisymmetric large-scale magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like features exhibited by these oscillations, we carry out an analysis of the large-scale zonal dynamics. We demonstrate that simulated torsional oscillations are not driven primarily by the periodically-varying large-scale magnetic torque, as one might have expected, but rather via the magnetic modulation of angular-momentum transport by the large-scale meridional flow. This result is confirmed by a straightforward energy analysis. We also detect a fairly sharp transition in rotational dynamics taking place as one moves from the base of the convecting layers to the base of the thin tachocline-like shear layer formed in the stably stratified fluid layers immediately below. We conclude by discussing the implications of our analyses with regards to the mechanism of amplitude saturation in the global dynamo operating in the simulation, and speculate on the possible precursor value of torsional oscillations for the forecast of solar cycle characteristics.