The present study aims at studying the flow and field produced by a stellar radiative zone which is initially made to rotate differentially in the presence of a large-scale poloidal magnetic field threading the whole domain. We focus both on the axisymmetric configurations produced by the initial winding-up of the magnetic field lines and on the possible instabilities of those configurations. We aim in particular at assessing the role of the stratification at stabilising the system. We perform 2D and 3D global Boussinesq numerical simulations started from an initial radial or cylindrical differential rotation and a large-scale poloidal magnetic field. Under the conditions of a large rotation frequency compared to the Alfv\'en frequency, a magnetic configuration strongly dominated by its toroidal component is built. The parameters of the simulations are chosen to respect the ordering of time scales of a typical stellar radiative zone. In this framework, the axisymmetric evolution is studied by varying the relative effects of the thermal diffusion, the Br\"unt-V\"ais\"al\"a frequency, the rotation and the initial poloidal field strength. We find that the axisymmetric state only depends on $t_{es}/t_{Ap}$, the ratio between the Eddington-Sweet circulation time scale and the Alfv\'en time scale. A scale analysis of the Boussinesq equations allows us to recover this result. In the cylindrical case, a magneto-rotational instability develops when the thermal diffusivity is sufficiently high to enable the favored wavenumbers to be insensitive to the effects of the stable stratification. In the radial case, the magneto-rotational instability is driven by the latitudinal shear created by the back-reaction of the Lorentz force on the flow. Increasing the level of stratification then leaves the growth rate of the instability mainly unaffected while its horizontal length scale grows.

Magnetic fields vary in complexity for different stars. The stability of dipolar magnetic fields is known to depend on different quantities, e.g., the stellar rotation, the stratification, and the intensity of convective motions. Here, we study the dipolar stability in a system with an inner stable zone. We present preliminary results of dynamo simulations using the Rayleigh number as a control parameter. The stiffness of the stable zone is accordingly varied to keep a constant ratio of the Brunt-Vaisala frequency to the angular velocity. Similarly to the completely convective spherical shell, we find that a transition exists between a regime where the magnetic field is dipolar to a multipolar regime when the Rossby number is increased. The value of the Rossby number at the transition is very close to the one of the fully convective case.

We report the analysis, conducted as part of the MaTYSSE programme, of a spectropolarimetric monitoring of the ~0.8 Myr, ~1.4 MSun disc-less weak-line T Tauri star V410 Tau with the ESPaDOnS instrument at the Canada-France-Hawaii Telescope and NARVAL at the T\'elescope Bernard Lyot, between 2008 and 2016. With Zeeman-Doppler Imaging, we reconstruct the surface brightness and magnetic field of V410 Tau, and show that the star is heavily spotted and possesses a ~550 G relatively toroidal magnetic field. We find that V410 Tau features a weak level of surface differential rotation between the equator and pole ~5 times weaker than the solar differential rotation. The spectropolarimetric data exhibit intrinsic variability, beyond differential rotation, which points towards a dynamo-generated field rather than a fossil field. Long-term variations in the photometric data suggest that spots appear at increasing latitudes over the span of our dataset, implying that, if V410 Tau has a magnetic cycle, it would have a period of more than 8 years. Having derived raw radial velocities (RVs) from our spectra, we filter out the stellar activity jitter, modeled either from our Doppler maps or using Gaussian Process Regression. Thus filtered, our RVs exclude the presence of a hot Jupiter-mass companion below ~0.1 au, which is suggestive that hot Jupiter formation may be inhibited by the early depletion of the circumstellar disc, which for V410 Tau may have been caused by the close (few tens of au) M dwarf stellar companion.

Magnetohydrodynamic interactions between plasma flows and magnetic fields is fundamental to the origin and sustenance of the 11-year sunspot cycle. These processes are intrinsically three-dimensional (3D) in nature. Our goal is to construct a 3D solar dynamo model that on the one hand captures the buoyant emergence of tilted bipolar sunspot pairs, and on the other hand produces cyclic large-scale field reversals mediated via surface flux-transport processes -- that is, the Babcock-Leighton mechanism. We perform kinematic dynamo simulations where the prescribed velocity field is a combination of solar-like differential rotation and meridional circulation, along with a parametrized turbulent diffusivity. We use a novel methodology for modeling magnetic buoyancy through field-strength-dependent 3D helical up-flows that results in the formation of tilted bipolar sunspots. The bipolar spots produced in our simulations participate in the process of poloidal-field generation through the Babcock-Leighton mechanism, resulting in self-sustained and periodic large-scale magnetic field reversal. Our parameter space study varying the amplitude of the meridional flow, the convection zone diffusivity, and parameters governing the efficiency of the magnetic buoyancy mechanism reveal their relative roles in determining properties of the sunspot cycle such as amplitude, period, and dynamical memory relevant to solar cycle prediction. We also derive a new dynamo number for the Babcock-Leighton solar dynamo mechanism which reasonably captures our model dynamics. This study elucidates the relative roles of different flux-transport processes in the Sun's convection zone in determining the properties and physics of the sunspot cycle and could potentially lead to realistic, data-driven 3D dynamo models for solar-activity predictions and exploration of stellar magnetism and starspot formation in other stars.

This study aims at investigating the ability of strong interactions between magnetic field concentrations during their rise through the convection zone to produce complex active regions at the solar surface. To do so, we perform numerical simulations of buoyant magnetic structures evolving and interacting in a model solar convection zone. We first produce a 3D model of rotating convection and then introduce idealized magnetic structures close to the bottom of the computational domain. These structures possess a certain degree of field line twist and they are made buoyant on a particular extension in longitude. The resulting twisted $\Omega$-loops will thus evolve inside a spherical convective shell possessing large-scale mean flows. We present results on the interaction between two such loops with various initial parameters (mainly buoyancy and twist) and on the complexity of the emerging magnetic field. In agreement with analytical predictions, we find that if the loops are introduced with opposite handedness and same axial field direction or same handedness but opposite axial field, they bounce against each other. The emerging region is then constituted of two separated bipolar structures. On the contrary, if the loops are introduced with the same direction of axial and peripheral magnetic fields and if sufficiently close, they merge while rising. This more interesting case produces complex magnetic structures, with a high degree of non-neutralized currents, especially when the convective motions act significantly on the magnetic field. This indicates that those interactions could be good candidates to produce eruptive events like flares or CMEs.

We present a three-dimensional numerical model for the generation and evolution of the magnetic field in the solar convection zone, in which sunspots are produced and contribute to the cyclic reversal of the large-scale magnetic field. We then assess the impact of this dynamo-generated field on the structure of the solar corona and solar wind. This model solves the induction equation in which the velocity field is prescribed. This velocity field is a combination of a solar-like differential rotation and meridional circulation. We develop an algorithm that enables the magnetic flux produced in the interior to be buoyantly transported towards the surface to produce bipolar spots. We find that those tilted bipolar magnetic regions contain a sufficient amount of flux to periodically reverse the polar magnetic field and sustain dynamo action. We then track the evolution of these magnetic features at the surface during a few consecutive magnetic cycles and analyze their effects on the topology of the corona and on properties of the solar wind (distribution of streamers and coronal holes, and of slow and fast wind streams) in connection with current observations of the Sun.

We present in this work the development of a solar data assimilation method based on an axisymmetric mean field dynamo model and magnetic surface data, our mid-term goal is to predict the solar quasi cyclic activity. Here we focus on the ability of our algorithm to constrain the deep meridional circulation of the Sun based on solar magnetic observations. To that end, we develop a variational data assimilation technique. Within a given assimilation window, the assimilation procedure minimizes the differences between data and the forecast from the model, by finding an optimal meridional circulation in the convection zone, and an optimal initial magnetic field, via a quasi-Newton algorithm. We demonstrate the capability of the technique to estimate the meridional flow by a closed-loop experiment involving 40 years of synthetic, solar-like data. By assimilating the synthetic magnetic proxies annually, we are able to reconstruct a (stochastic) time-varying meridional circulation which is also slightly equatorially asymmetric. We show that the method is robust in estimating a flow whose level of fluctuation can reach 30% about the average, and that the horizon of predictive capability of the method is of the order of 1 cycle length.

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow which best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires the knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique, for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations, by varying their spatial density, temporal density, noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step towards using this technique to i) better constrain the physical processes occurring inside the Sun and ii) better predict solar activity on decadal time scales.

We review recent advances and results in enhancing and developing helioseismic analysis methods and in solar data assimilation. In the first part of this paper we will focus on selected developments in time-distance and global helioseismology. In the second part, we review the application of data assimilation methods on solar data. Relating solar surface observations as well as helioseismic proxies with solar dynamo models by means of the techniques from data assimilation is a promising new approach to explore and to predict the magnetic activity cycle of the Sun.

An extremely weak circularly polarized signature was recently discovered in spectral lines of the chemically peculiar Am star Sirius A. A weak surface magnetic field was proposed to account for the observed polarized signal, but the shape of the phase-averaged signature, dominated by a prominent positive lobe, is not expected in the standard theory of the Zeeman effect. We aim at verifying the presence of weak circularly polarized signatures in two other bright Am stars, beta UMa and theta Leo, and investigating the physical origin of Sirius-like polarized signals further. We present here a set of deep spectropolarimetric observations of beta UMa and theta Leo, observed with the NARVAL spectropolarimeter. We analyzed all spectra with the Least Squares Deconvolution multiline procedure. To improve the signal-to-noise ratio and detect extremely weak signatures in Stokes V profiles, we co-added all available spectra of each star (around 150 observations each time). Finally, we ran several tests to evaluate whether the detected signatures are consistent with the behavior expected from the Zeeman effect. The line profiles of the two stars display circularly polarized signatures similar in shape and amplitude to the observations previously gathered for Sirius A. Our series of tests brings further evidence of a magnetic origin of the recorded signal. These new detections suggest that very weak magnetic fields may well be present in the photospheres of a significant fraction of intermediate-mass stars. The strongly asymmetric Zeeman signatures measured so far in Am stars (featuring a dominant single-sign lobe) are not expected in the standard theory of the Zeeman effect and may be linked to sharp vertical gradients in photospheric velocities and magnetic field strengths.

We study the magnetic flux carried by pores located outside active regions with sunspots and investigate their possible contribution to the reversal of the global magnetic field of the Sun. We find that they contain a total flux of comparable amplitude to the total magnetic flux contained in polar caps. The pores located at distances of 40--100~Mm from the closest active region have systematically the correct sign to contribute to the polar cap reversal. These pores can predominantly be found in bipolar magnetic regions. We propose that during grand minima of solar activity, such a systematic polarity trend, akin to a weak magnetic (Babcock-Leighton-like) source term could still be operating but was missed by the contemporary observers due to the limited resolving power of their telescopes.

Recent spectropolarimetric surveys of main-sequence intermediate-mass stars have exhibited a dichotomy in the distribution of the observed magnetic field between the kG dipoles of Ap/Bp stars and the sub-Gauss magnetism of Vega and Sirius. We would like to test whether this dichotomy is linked to the stability versus instability of large-scale magnetic configurations in differentially rotating radiative zones. We computed the axisymmetric magnetic field obtained from the evolution of a dipolar field threading a differentially rotating shell. A full parameter study including various density profiles and initial and boundary conditions was performed with a 2D numerical code. We then focused on the ratio between the toroidal and poloidal components of the magnetic field and discuss the stability of the configurations dominated by the toroidal component using local stability criteria and insights from recent 3D numerical simulations. The numerical results and a simple model show that the ratio between the toroidal and the poloidal magnetic fields is highest after an Alfv\'en crossing time of the initial poloidal field. For high density contrasts, this ratio converges towards an asymptotic value that can thus be extrapolated to realistic stellar cases. We then consider the stability of the magnetic configurations to non-axisymmetric perturbations and find that configurations dominated by the toroidal component are likely to be unstable if the shear strength is significantly higher than the poloidal Alfv\'en frequency. An expression for the critical poloidal field below which magnetic fields are likely to be unstable is found and is compared to the lower bound of Ap/Bp magnetic fields.

The question of the origin and evolution of magnetic fields in stars possessing a radiative envelope, like the A-type stars, is still regarded as a challenge for stellar physics. Those zones are likely to be differentially rotating, which suggests that strong interactions between differential rotation and magnetic fields could be at play. We numerically compute the joint evolution of the magnetic and velocity fields in a 3D spherical shell starting from an initial profile for the poloidal magnetic field and differential rotation. The poloidal magnetic field is initially wound-up by the differential rotation to produce a toroidal field which becomes unstable. In the particular setup studied here where the differential rotation is dominant, the magneto-rotational instability is triggered. The growth rate of the instability depends mainly on the initial rotation rate, while the background state typically oscillates over a poloidal Alfv\'en time. We thus find that the axisymmetric magnetic configuration is strongly modified by the instability only if the ratio between the poloidal Alfv\'en frequency and the rotation rate is sufficiently small. An enhanced transport of angular momentum is found in the most unstable cases: the typical time to flatten the rotation profile is then much faster than the diffusion time scale. We conclude that the magneto-rotational instability is always favored (over the Tayler instability) in unstratified spherical shells when an initial poloidal field is sheared by a sufficiently strong cylindrical differential rotation. A possible application to the magnetic desert observed among A stars is given. We argue that the dichotomy between stars exhibiting strong axisymmetric fields (Ap stars) and those harboring a sub-Gauss magnetism could be linked to the threshold for the instability.

In a uniformly rotating fluid, inertial waves propagate along rays that are inclined to the rotation axis by an angle that depends on the wave frequency. In closed domains, multiple reflections from the boundaries may cause inertial waves to focus on to particular structures known as wave attractors. Such structures have previously been studied from a theoretical point of view, in laboratory experiments, in linear numerical calculations and in some recent numerical simulations. In the present paper, two-dimensional direct numerical simulations of an inertial wave attractor are presented. In the linear regime, we first recover the results of the linear calculations and asymptotic theory of Ogilvie (2005) who considered a prototypical problem involving the focusing of linear internal waves into a narrow beam centred on a wave attractor in a steady state. The velocity profile of the beam and its scalings with the Ekman number, as well as the asymptotic value of the dissipation rate, are found to be in agreement with the linear theory. We also find that, as the beam builds up around the wave attractor, the power input by the applied force reaches its limiting value more rapidly than the dissipation rate, which saturates only when the beam has reached its final thickness. In the nonlinear regime, the beam is strongly affected and becomes unstable to a subharmonic instability. This instability transfers energy to secondary waves possessing shorter wavelengths and lower frequencies. The onset of the instability of a narrow inertial wave beam is investigated by means of a separate linear analysis and the results, such as the onset of the instability, are found to be consistent with the global simulations of the wave attractor. The excitation of such secondary waves described theoretically in this work has also been seen in recent laboratory experiments on internal gravity waves.

We present three-dimensional numerical simulations of a magnetic loop evolving in either a convectively stable or unstable rotating shell. The magnetic loop is introduced in the shell in such a way that it is buoyant only in a certain portion in longitude, thus creating an \Omega-loop. Due to the action of magnetic buoyancy, the loop rises and develops asymmetries between its leading and following legs, creating emerging bipolar regions whose characteristics are similar to the ones of observed spots at the solar surface. In particular, we self-consistently reproduce the creation of tongues around the spot polarities, which can be strongly affected by convection. We moreover emphasize the presence of ring-shaped magnetic structures around our simulated emerging regions, which we call "magnetic necklace" and which were seen in a number of observations without being reported as of today. We show that those necklaces are markers of vorticity generation at the periphery and below the rising magnetic loop. We also find that the asymmetry between the two legs of the loop is crucially dependent on the initial magnetic field strength. The tilt angle of the emerging regions is also studied in the stable and unstable cases and seems to be affected both by the convective motions and the presence of a differential rotation in the convective cases.

We here report calculations of magnetic buoyancy instabilities of a sheared magnetic layer where two separate zones are unstable. The idea is to study the possible generation of large-scale helical structures which could then rise through a stellar convection zone and emerge at the surface to create active regions. The calculations shown here are a follow-up of the work of Favier et al. (2012) where the instability developed in a weakly magnetized atmosphere, consisting of a uniform field oriented in a different direction from the unstable layer below. Here, the top layer representing the atmosphere is itself unstable to buoyancy instabilities and thus quickly creates a more complex magnetic configuration with which the layer below will interact. We also find in this case that the accumulation of magnetic tension between the two unstable layers favors the creation of large-scale helical structures.

We consider the three-dimensional instability of a layer of horizontal magnetic field in a polytropic atmosphere where, contrary to previous studies, the field lines in the initial state are not unidirectional. We show that if the twist is initially concentrated inside the unstable layer, the modifications of the instability reported by several authors (see e.g. Cattaneo et al. (1990)) are only observed when the calculation is restricted to two dimensions. In three dimensions, the usual interchange instability occurs, in the direction fixed by the field lines at the interface between the layer and the field-free region. We therefore introduce a new configuration: the instability now develops in a weakly magnetised atmosphere where the direction of the field can vary with respect to the direction of the strong unstable field below, the twist being now concentrated at the upper interface. Both linear stability analysis and non-linear direct numerical simulations are used to study this configuration. We show that from the small-scale interchange instability, large-scale twisted coherent magnetic structures are spontaneously formed, with possible implications to the formation of active regions from a deep-seated solar magnetic field.

We aim at a highly sensitive search for weak magnetic fields in main sequence stars of intermediate mass, by scanning classes of stars with no previously reported magnetic members. After detecting a weak magnetic field on the normal, rapidly rotating A-type star Vega, we concentrate here on the bright star Sirius A, taken as a prototypical, chemically peculiar, moderately rotating Am star. We employed the NARVAL and ESPaDOnS high-resolution spectropolarimeters to collect 442 circularly polarized spectra, complemented by 60 linearly polarized spectra. Using a list of about 1,100 photospheric spectral lines, we computed a cross correlation line profile from every spectrum, leading to a signal-to-noise ratio of up to 30,000 in the polarized profile. We report the repeated detection of circularly polarized, highly asymmetric signatures in the line profiles, interpreted as Zeeman signatures of a large-scale photospheric magnetic field, with a line-of-sight component equal to $0.2 \pm 0.1$ G. This is the first polarimetric detection of a surface magnetic field on an Am star. Using rough estimates of the physical properties of the upper layers of Sirius A, we suggest that a dynamo operating in the shallow convective envelope cannot account for the field strength reported here. Together with the magnetic field of Vega, this result confirms that a new class of magnetic objects exists among non Ap/Bp intermediate-mass stars, and it may indicate that a significant fraction of tepid stars are magnetic.

We study the connections between the sun's convection zone and the evolution of the solar wind and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code (STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in good agreement with the available observations. The magnetic polarity reversal happens at different paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal distribution of the sources of open flux and solar wind and the geometry of the Alfv\'en surface. Wind sources and braking torque application zones also vary accordingly.

We have developed a variational data assimilation technique for the Sun using a toy {\alpha}{\Omega} dynamo model. The purpose of this work is to apply modern data assimilation techniques to solar data using a physically based model. This work represents the first step toward a complete variational model of solar magnetism. We derive the adjoint {\alpha}{\Omega} dynamo code and use a minimization procedure to invert the spatial dependence of key physical ingredients of the model. We find that the variational technique is very powerful and leads to encouraging results that will be applied to a more realistic model of the solar dynamo. We also note that the continuous adjoint mean-field induction equation is derived in the appendix.

The Sun is a magnetic star whose cyclic activity is thought to be linked to internal dynamo mechanisms. A combination of numerical modelling with various levels of complexity is an efficient and accurate tool to investigate such intricate dynamical processes. We investigate the role of the magnetic buoyancy process in 2D Babcock-Leighton dynamo models, by modelling more accurately the surface source term for poloidal field. Methods. To do so, we reintroduce in mean-field models the results of full 3D MHD calculations of the non-linear evolution of a rising flux tube in a convective shell. More specifically, the Babcock-Leighton source term is modified to take into account the delay introduced by the rise time of the toroidal structures from the base of the convection zone to the solar surface. We find that the time delays introduced in the equations produce large temporal modulation of the cycle amplitude even when strong and thus rapidly rising flux tubes are considered. Aperiodic modulations of the solar cycle appear after a sequence of period doubling bifurcations typical of non-linear systems. The strong effects introduced even by small delays is found to be due to the dependence of the delays on the magnetic field strength at the base of the convection zone, the modulation being much less when time delays remain constant. We do not find any significant influence on the cycle period except when the delays are made artificially strong. A possible new origin of the solar cycle variability is here revealed. This modulated activity and the resulting butterfly diagram are then more compatible with observations than what the standard Babcock-Leighton model produces.

Aims: To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field dynamo models with stellar activity observations Methods: 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational data. The novelty is that these 2-D mean field dynamo models incorporate scaling laws deduced from 3-D hydrodynamical simulations for the influence of rotation rate on the amplitude and profile of the meridional circulation. These models make also use of observational scaling laws for the variation of differential rotation with rotation rate. Results: We find that Babcock-Leighton flux transport dynamo models are able to reproduce the change in topology of the magnetic field (i.e. toward being more toroidal with increasing rotation rate) but seem to have difficulty reproducing the cycle period vs activity period correlation observed in solar-like stars if a monolithic single cell meridional flow is assumed. It may however be possible to recover the $P_{cyc}$ vs $P_{rot}$ relation with more complex meridional flows, if the profile changes in a particular assumed manner with rotation rate. Conclusions: The Babcock-Leighton flux transport dynamo model based on single cell meridional circulation does not reproduce the $P_{cyc}$ vs $P_{rot}$ relation unless the amplitude of the meridional circulation is assumed to increase with rotation rate which seems to be in contradiction with recent results obtained with 3-D global simulations.

We present the first 3D MHD study in spherical geometry of the non-linear dynamical evolution of magnetic flux tubes in a turbulent rotating convection zone. We study numerically the rise of magnetic toroidal flux ropes from the base of a modelled convection zone up to the top of our computational domain where bipolar patches are formed. We compare the dynamical behaviour of flux tubes in a fully convective shell possessing self-consistently generated mean flows such as meridional circulation and differential rotation, with reference calculations done in a quiet isentropic zone. We find that two parameters influence the tubes during their rise through the convection zone: the initial field strength and amount of twist, thus confirming previous findings in Cartesian geometry. Further, when the tube is sufficiently strong with respect to the equipartition field, it rises almost radially independently of the initial latitude (either low or high). By contrast, weaker field cases indicate that downflows and upflows control the rising velocity of particular regions of the rope and could in principle favour the emergence of flux through Omega-loop structures. For these latter cases, we focus on the orientation of bipolar patches and find that sufficiently arched structures are able to create bipolar regions with a predominantly East-West orientation. Meridional flow seems to determine the trajectory of the magnetic rope when the field strength has been significantly reduced near the top of the domain. Finally differential rotation makes it more difficult for tubes introduced at low latitudes to reach the top of the domain.

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code. We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes.

The Sun is a magnetic star whose magnetism and cyclic activity is linked to the existence of an internal dynamo. We aim to understand the establishment of the solar magnetic 22-yr cycle, its associated butterfly diagram and field parity selection through numerical simulations of the solar global dynamo. Inspired by recent observations and 3D simulations that both exhibit multicellular flows in the solar convection zone, we seek to characterise the influence of various profiles of circulation on the behaviour of solar mean-field dynamo models. We are using 2-D mean field flux transport Babcock-Leighton numerical models in which we test several types of meridional flows: 1 large single cell, 2 cells in radius and 4 cells per hemisphere. We confirm that adding cells in latitude tends to speed up the dynamo cycle whereas adding cells in radius more than triples the period. We find that the cycle period in the four cells model is less sensitive to the flow speed than in the other simpler meridional circulation profiles studied. Moreover, our studies show that adding cells in radius or in latitude seems to favour the parity switching to a quadrupolar solution. According to our numerical models, the observed 22-yr cycle and dipolar parity is easily reproduced by models including multicellular meridional flows. On the contrary, the resulting butterfly diagram and phase relationship between the toroidal and poloidal fields are affected to a point where it is unlikely that such multicellular meridional flows persist for a long period of time inside the Sun, without having to reconsider the model itself.