These lecture notes are based on a tutorial given in 2017 at a plasma physics winter school in Les Houches. Their aim is to provide a self-contained graduate-student level introduction to the theory and modelling of the dynamo effect in turbulent fluids and plasmas, blended with a review of current research in the field. The primary focus is on the physical and mathematical concepts underlying different (turbulent) branches of dynamo theory, with some astrophysical, geophysical and experimental context disseminated throughout the document. The text begins with an introduction to the rationale, observational and historical roots of the subject, and to the basic concepts of magnetohydrodynamics relevant to dynamo theory. The next two sections discuss the fundamental phenomenological and mathematical aspects of (linear and nonlinear) small- and large-scale MHD dynamos. These sections are complemented by an overview of a selection of current active research topics in the field, including the numerical modelling of the geo- and solar dynamos, shear dynamos driven by turbulence with zero net helicity, and MHD-instability-driven dynamos such as the magnetorotational dynamo. The difficult problem of a unified, self-consistent statistical treatment of small and large-scale dynamos at large magnetic Reynolds numbers is also discussed throughout the text. Finally, an excursion is made into the relatively new but increasingly popular realm of magnetic-field generation in weakly-collisional plasmas. A short discussion of the outlook and challenges for the future of the field concludes the presentation.

Supergranulation is a fluid-dynamical phenomenon taking place in the solar photosphere, primarily detected in the form of a vigorous cellular flow pattern with a typical horizontal scale of approximately 30--35~megameters, a dynamical evolution time of 24--48~h, a strong 300--400~m/s (rms) horizontal flow component and a much weaker 20--30~m/s vertical component. Supergranulation was discovered more than sixty years ago, however, explaining its physical origin and most important observational characteristics has proven extremely challenging ever since, as a result of the intrinsic multiscale, nonlinear dynamical complexity of the problem concurring with strong observational and computational limitations. Key progress on this problem is now taking place with the advent of 21st-century supercomputing resources and the availability of global observations of the dynamics of the solar surface with high spatial and temporal resolutions. This article provides an exhaustive review of observational, numerical and theoretical research on supergranulation, and discusses the current status of our understanding of its origin and dynamics, most importantly in terms of large-scale nonlinear thermal convection, in the light of a selection of recent findings.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

The Sun provides us with the only spatially well-resolved astrophysical example of turbulent thermal convection. While various aspects of solar photospheric turbulence, such as granulation (one-Megameter horizontal scale), are well understood, the questions of the physical origin and dynamical organization of larger-scale flows, such as the 30-Megameters supergranulation and flows deep in the solar convection zone, remain largely open in spite of their importance for solar dynamics and magnetism. Here, we present a new critical global observational characterization of multiscale photospheric flows and subsequently formulate an anisotropic extension of the Bolgiano-Obukhov theory of hydrodynamic stratified turbulence that may explain several of their distinctive dynamical properties. Our combined analysis suggests that photospheric flows in the horizontal range of scales between supergranulation and granulation have a typical vertical correlation scale of 2.5 to 4 Megameters and operate in a strongly anisotropic, self-similar, nonlinear, buoyant dynamical regime. While the theory remains speculative at this stage, it lends itself to quantitative comparisons with future high-resolution acoustic tomography of subsurface layers and advanced numerical models. Such a validation exercise may also lead to new insights into the asymptotic dynamical regimes in which other, unresolved turbulent anisotropic astrophysical fluid systems supporting waves or instabilities operate.

In Keplerian accretion disks, turbulence and magnetic fields may be jointly excited through a subcritical dynamo process involving the magnetorotational instability (MRI). High-resolution simulations exhibit a tendency towards statistical self-organization of MRI dynamo turbulence into large-scale cyclic dynamics. Understanding the physical origin of these structures, and whether they can be sustained and transport angular momentum efficiently in astrophysical conditions, represents a significant theoretical challenge. The discovery of simple periodic nonlinear MRI dynamo solutions has recently proven useful in this respect, and has notably served to highlight the role of turbulent magnetic diffusion in the seeming decay of the dynamics at low magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a common regime in accretion disks. The connection between these simple structures and the statistical organization reported in turbulent simulations remained elusive, though. Here, we report the numerical discovery in moderate aspect ratio Keplerian shearing boxes of new periodic, incompressible, three-dimensional nonlinear MRI dynamo solutions with a larger dynamical complexity reminiscent of such simulations. These "chimera" cycles are characterized by multiple MRI-unstable dynamical stages, but their basic physical principles of self-sustainment are nevertheless identical to those of simpler cycles found in azimuthally elongated boxes. In particular, we find that they are not sustained at low Pm either due to subcritical turbulent magnetic diffusion. These solutions offer a new perspective into the transition from laminar to turbulent instability-driven dynamos, and may prove useful to devise improved statistical models of turbulent accretion disk dynamos.

A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfv\'en wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at $\beta\gtrsim1$ and by magnetosonic/whistler fluctuations at lower $\beta$. The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma $\beta$, and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions.

Magnetic fields pervade the entire Universe and affect the formation and evolution of astrophysical systems from cosmological to planetary scales. The generation and dynamical amplification of extragalactic magnetic fields through cosmic times, up to $\mu$Gauss levels reported in nearby galaxy clusters, near equipartition with kinetic energy of plasma motions and on scales of at least tens of kiloparsecs, is a major puzzle largely unconstrained by observations. A dynamo effect converting kinetic flow energy into magnetic energy is often invoked in that context, however extragalactic plasmas are weakly collisional (as opposed to magnetohydrodynamic fluids), and whether magnetic-field growth and sustainment through an efficient turbulent dynamo instability is possible in such plasmas is not established. Fully kinetic numerical simulations of the Vlasov equation in a six-dimensional phase space necessary to answer this question have until recently remained beyond computational capabilities. Here, we show by means of such simulations that magnetic-field amplification via a dynamo instability does occur in a stochastically-driven, non-relativistic subsonic flow of initially unmagnetized collisionless plasma. We also find that the dynamo self-accelerates and becomes entangled with kinetic instabilities as magnetization increases. The results suggest that such a plasma dynamo may be realizable in laboratory experiments, support the idea that intracluster medium (ICM) turbulence may have significantly contributed to the amplification of cluster magnetic fields up to near-equipartition levels on a timescale shorter than the Hubble time, and emphasize the crucial role of multiscale kinetic physics in high-energy astrophysical plasmas.

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

Slow dynamical changes in magnetic-field strength and invariance of the particles' magnetic moments generate ubiquitous pressure anisotropies in weakly collisional, magnetized astrophysical plasmas. This renders them unstable to fast, small-scale mirror and firehose instabilities, which are capable of exerting feedback on the macroscale dynamics of the system. By way of a new asymptotic theory of the early nonlinear evolution of the mirror instability in a plasma subject to slow shearing or compression, we show that the instability does not saturate quasilinearly at a steady, low-amplitude level. Instead, the trapping of particles in small-scale mirrors leads to nonlinear secular growth of magnetic perturbations, $\delta B/B \propto t^{2/3}$. Our theory explains recent collisionless simulation results, provides a prediction of the mirror evolution in weakly collisional plasmas and establishes a foundation for a theory of nonlinear mirror dynamics with trapping, valid up to $\delta B/B =O(1)$.

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, nonlinear magnetohydrodynamic process whose study is relevant to the understanding of accretion and magnetic field generation in astrophysics. Transition to this form of dynamo is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. This suggests that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in various astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows.

We briefly review the observations of the solar photosphere and pinpoint some open questions related to the magnetohydrodynamics of this layer of the Sun. We then discuss the current modelling efforts, addressing among other problems, that of the origin of supergranulation.

The nature of dynamo action in shear flows prone to magnetohydrodynamic instabilities is investigated using the magnetorotational dynamo in Keplerian shear flow as a prototype problem. Using direct numerical simulations and Newton's method, we compute an exact time-periodic magnetorotational dynamo solution to the three-dimensional dissipative incompressible magnetohydrodynamic equations with rotation and shear. We discuss the physical mechanism behind the cycle and show that it results from a combination of linear and nonlinear interactions between a large-scale axisymmetric toroidal magnetic field and non-axisymmetric perturbations amplified by the magnetorotational instability. We demonstrate that this large scale dynamo mechanism is overall intrinsically nonlinear and not reducible to the standard mean-field dynamo formalism. Our results therefore provide clear evidence for a generic nonlinear generation mechanism of time-dependent coherent large-scale magnetic fields in shear flows and call for new theoretical dynamo models. These findings may offer important clues to understand the transitional and statistical properties of subcritical magnetorotational turbulence.

Weakly collisional plasmas dynamically develop pressure anisotropies with respect to the magnetic field. These anisotropies trigger plasma instabilities at scales just above the ion Larmor radius \rho_i and much below the mean free path \lambda_{mfp}. They have growth rates of a fraction of the ion cyclotron frequency - much faster than either the global dynamics or local turbulence. The instabilities dramatically modify the transport properties and, therefore, the macroscopic dynamics of the plasma. Their nonlinear evolution drives pressure anisotropies towards marginal stability, controlled by the plasma beta \beta_i. Here this nonlinear evolution is worked out for the simplest analytically tractable example - the parallel firehose instability. In the nonlinear regime, both analytical theory and the numerical solution predict secular growth of magnetic fluctuations. They develop a k^{-3} spectrum, extending from scales somewhat larger than \rho_i to the maximum scale that grows secularly with time (~t^{1/2}); the relative pressure anisotropy (\pperp-\ppar)/\ppar tends to the marginal value -2/\beta_i. The marginal state is achieved via changes in the magnetic field, not particle scattering. When a parallel ion heat flux is present, the firehose mutates into the new gyrothermal instability (GTI), which continues to exist up to firehose-stable values of pressure anisotropy, which can be positive and are limited by the heat flux. The nonlinear evolution of the GTI also features secular growth of magnetic fluctuations, but the spectrum is eventually dominated by modes around the scale ~\rho_i l_T/\lambda_{mfp}, where l_T is the scale of the parallel temperature variation. Implications for momentum and heat transport are speculated about. This study is motivated by the dynamics of galaxy cluster plasmas.

Both global dynamics and turbulence in magnetized weakly collisional cosmic plasmas are described by general magnetofluid equations that contain pressure anisotropies and heat fluxes that must be calculated from microscopic plasma kinetic theory. It is shown that even without a detailed calculation of the pressure anisotropy or the heat fluxes, one finds the macroscale dynamics to be generically unstable to microscale Alfvenically polarized fluctuations. Two instabilities are considered in detail: the parallel firehose instability (including the finite-Larmor-radius effects that determine the fastest growing mode) and the gyrothermal instability (GTI). The latter is a new result - it is shown that a parallel ion heat flux destabilizes Alfvenically polarized fluctuations even in the absence of the negative pressure anisotropy required for the firehose. The main conclusion is that both pressure anisotropies and heat fluxes trigger plasma microinstabilities and, therefore, their values will likely be set by the nonlinear evolution of these instabilities. Ideas for understanding this nonlinear evolution are discussed. It is argued that cosmic plasmas will generically be "three-scale systems," comprising global dynamics, mesoscale turbulence and microscale plasma fluctuations. The astrophysical example of cool cores of galaxy clusters is considered and it is noted that observations point to turbulence in clusters being in a marginal state with respect to plasma microinstabilities and so it is the plasma microphysics that is likely to set the heating and conduction properties of the intracluster medium. In particular, a lower bound on the scale of temperature fluctuations implied by the GTI is derived.

The aim of this work is to give new observational constraints on solar surface flows by determining the horizontal scale dependence of the velocity and intensity fields, as represented by their power spectra, and to offer some theoretical guidelines to interpret these spectra. We use long time series of images taken by SOT/Hinode and reconstruct both horizontal (by granule tracking) and vertical (by Doppler effect) velocity fields in a field of view 75x75Mm^2. At small sub-granulation scales, the kinetic energy spectral density associated with vertical motions exhibits a k^{-10/3}-like spectrum, while the intensity fluctuation spectrum follows a k^{-3} or k^{-17/3}-like spectrum at the two continuum levels investigated (525 and 450 nm respectively). We discuss the physical origin of these scalings and argue that they provide a direct observational signature of buoyancy-driven turbulent dynamics in a strongly thermally diffusive regime. In the mesogranulation range and up to a scale of 25Mm, we find that the vertical velocity field amplitude decreases like L^{-3/2} with the horizontal scale L. This behaviour corresponds to a k^2 spectral power law. Still in the mesoscale range, [2.5, 10]Mm, we find that intensity fluctuations in the blue continuum follow a power law in k^2. We show that granule tracking cannot sample scales below 2.5Mm. We locate the supergranulation energy peak around 30Mm and show that the emergence of a pore erases this spectral peak. Thanks to a scale height estimate, we find that supergranular flows are shallow. (abridged abstract).

Aims: The interactions of velocity scales on the Sun's surface, from granulation to supergranulation are still not understood, nor are their interaction with magnetic fields. We thus aim at giving a better description of dynamics in the mesoscale range which lies between the two scales mentioned above. Method: We analyse a 48h high-resolution time sequence of the quiet Sun photosphere at the disk center obtained with the Solar Optical Telescope onboard Hinode. The observations, which have a field of view of 100 \arcsec$\times$ 100 \arcsec, typically contain four supergranules. We monitor in detail the motion and evolution of granules as well as those of the radial magnetic field. Results: This analysis allows us to better characterize Trees of Fragmenting Granules issued from repeated fragmentation of granules, especially their lifetime statistics. Using floating corks advected by measured velocity fields, we show their crucial role in the advection of the magnetic field and in the build up of the network. Finally, thanks to the long duration of the time series, we estimate that the turbulent diffusion coefficient induced by horizontal motion is approximately $430 \mathrm{km}^2 \mathrm{s}^{-1}$. Conclusions: These results demonstrate that the long living families contribute to the formation of the magnetic network and suggest that supergranulation could be an emergent length scale building up as small magnetic elements are advected and concentrated by TFG flows. Our estimate for the magnetic diffusion associated with this horizontal motion might provide a useful input for mean-field dynamo models.

Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed shear, gamma ~ S, while the characteristic spatial scale of the field is l_b ~ S^(-1/2). The effect is quite general: earlier results for the nonrotating case by Yousef et al. 2008 (PRL 100, 184501) are extended to shearing boxes with Keplerian rotation; it is also shown that the shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of magnetic fields in astrophysical systems.

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non-rotating wall-bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non-axisymmetric instabilities of shear-induced axisymmetric toroidal velocity or magnetic fields, such as Kelvin-Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high-Reynolds number shear flows could act as a large-scale driving mechanism of turbulent flows that would in turn generate an independent small-scale dynamo.

A small fraction of intermediate-mass main sequence (A and B type) stars have strong, organised magnetic fields. The large majority of such stars, however, show no evidence for magnetic fields, even when observed with very high precision. In this paper we describe a simple model, motivated by qualitatively new observational results, that provides a natural physical explanation for the small fraction of observed magnetic stars.

We have investigated a sample of 28 well-known spectroscopically-identified magnetic Ap/Bp stars, with weak, poorly-determined or previously undetected magnetic fields, with the aim of exploring the weak part of the magnetic field distribution of Ap/Bp stars. Using the MuSiCoS and NARVAL spectropolarimeters we have obtained 282 LSD Stokes V signatures of our 28 sample stars. All stars were detected, showing clearly that when observed with sufficient precision, all firmly classified Ap/Bp stars show detectable surface magnetic fields. To better characterise the surface magnetic field intensities and geometries of the sample, we have inferred the dipolar field intensity and the magnetic obliquity. The distribution of derived dipole strengths for these stars exhibits a plateau at about 1 kG, falling off to larger and smaller field strengths. Remarkably, in this sample of stars selected for their presumably weak magnetic fields, we find only 2 stars for which the derived dipole strength is weaker than 300 G. We interpret this "magnetic threshold" as a critical value necessary for the stability of large-scale magnetic fields, and develop a simple quantitative model that is able to approximately reproduce the observed threshold characteristics. This scenario leads to a natural explanation of the small fraction of intermediate-mass magnetic stars. It may also explain the near-absence of magnetic fields in more massive B and O-type stars.

The problem of two-dimensional steady nonlinear dynamics in plane Couette flow is revisited using homotopy from either plane Poiseuille flow or from plane Couette flow perturbed by a small symmetry-preserving identity operator. Our results show that it is not possible to obtain the nonlinear plane Couette flow solutions reported by Cherhabili and Ehrenstein [Eur. J. Mech. B/Fluids, 14, 667 (1995)] using their Poiseuille-Couette homotopy. We also demonstrate that the steady solutions obtained by Mehta and Healey [Phys. Fluids, 17, 4108 (2005)] for small symmetry-preserving perturbations are influenced by an artefact of the modified system of equations used in their paper. However, using a modified version of their model does not help to find plane Couette flow solution in the limit of vanishing symmetry-preserving perturbations either. The issue of the existence of two-dimensional nonlinear steady states in plane Couette flow remains unsettled.

A three-dimensional nonlinear dynamo process is identified in rotating plane Couette flow in the Keplerian regime. It is analogous to the hydrodynamic self-sustaining process in non-rotating shear flows and relies on the magneto-rotational instability of a toroidal magnetic field. Steady nonlinear solutions are computed numerically for a wide range of magnetic Reynolds numbers but are restricted to low Reynolds numbers. This process may be important to explain the sustenance of coherent fields and turbulent motions in Keplerian accretion disks, where all its basic ingredients are present.

Subcritical transition to turbulence in Keplerian accretion disks is still a controversial issue and some theoretical progress is required in order to determine whether or not this scenario provides a plausible explanation for the origin of angular momentum transport in non-magnetized accretion disks. Motivated by the recent discoveries of exact nonlinear steady self-sustaining solutions in linearly stable non-rotating shear flows, we attempt to compute similar solutions in Rayleigh-stable rotating plane Couette flows and to identify transition mechanisms in such flows by combining nonlinear continuation methods and asymptotic theory. We obtain exact nonlinear solutions for Rayleigh-stable cyclonic regimes but show that it is not possible to compute solutions for Rayleigh-stable anticyclonic regimes, including Keplerian flow, using similar techniques. We also present asymptotic descriptions of these various problems at large Reynolds numbers that provide some insight into the differences between the non-rotating and Rayleigh-stable anticyclonic regimes and derive some necessary conditions for mechanisms analogous to the non-rotating self-sustaining process to be present in flows on the Rayleigh line. Our results demonstrate that subcritical transition mechanisms cannot be identified in wall-bounded Rayleigh-stable anticyclonic shear flows by transposing directly the phenomenology of subcritical transition in cyclonic and non-rotating wall-bounded shear flows. Asymptotic developments, however, leave open the possibility that nonlinear self-sustaining solutions may exist in unbounded or periodic flows on the Rayleigh line. These could serve as a starting point to discover solutions in Rayleigh-stable flows, but the nonlinear stability of Keplerian accretion disks remains to be determined.

We present the results of a very large aspect ratio (42.6) numerical simulation of fully compressible turbulent convection in a polytropic atmosphere, and focus on the properties of large-scale flows. Mesoscale patterns dominate the turbulent energy spectrum. We show that these structures, which had already been observed in Boussinesq simulations by Cattaneo et al. (2001), have a genuine convective origin and do not result directly from collective interactions of the smaller scales of the flow, even though their growth is strongly affected by nonlinear transfers. If this result is relevant to the solar photosphere, it suggests that the dominant convective mode below the Sun's surface may be at mesoscales.

Classical theories of turbulence do not describe accurately inertial range scaling laws in turbulent convection and notably fail to model the shape of the turbulent spectrum of solar photospheric convection. To understand these discrepancies, a detailed study of scale-by-scale budgets in turbulent Rayleigh-B\'enard convection is presented, with particular emphasis placed on anisotropy and inhomogeneity. A generalized Kolmogorov equation applying to convection is derived and its various terms are computed using numerical simulations of turbulent Boussinesq convection. The analysis of the isotropic part of the equation shows that the third-order velocity structure function is significantly affected by buoyancy forcing and large-scale inhomogeneities. Anisotropic contributions to this equation are also shown to be comparable to their isotropic counterpart at moderate to large scales. Implications of these results for convection in the solar photosphere, mesogranulation and supergranulation are discussed.

This paper examines the consistency of the exact scaling laws for isotropic MHD turbulence in numerical simulations with large magnetic Prandtl numbers Pm and with Pm=1. The exact laws are used to elucidate the structure of the magnetic and velocity fields. Despite the linear scaling of certain third-order correlation functions, the situation is not analogous to the case of Kolmogorov turbulence. The magnetic field is adequately described by a model of stripy (folded) field with direction reversals at the resistive scale. At currently available resolutions, the cascade of kinetic energy is short-circuited by the direct exchange of energy between the forcing-scale motions and the stripy magnetic fields. This nonlocal interaction is the defining feature of isotropic MHD turbulence.

This paper provides a detailed study of scale-by-scale budgets in turbulent Rayleigh-B\'enard convection and aims at testing the applicability of Kolmogorov (1941) and Bolgiano (1959) theories for this flow. Particular emphasis is laid on anisotropic and inhomogeneous effects: the SO(3) decomposition of structure functions (Arad et al 1999) and a method of description of inhomogeneities proposed by Danaila et al (2001) are used to derive inhomogeneous and anisotropic generalizations of Kolmogorov and Yaglom equations applying to RB convection. The various terms in these equations are computed using data from a DNS of turbulent Boussinesq convection at $\rayleigh=10^6$ and $\prandtl=1$ with aspect ratio A=5. The analysis of the isotropic component demonstrates that the shape of the third-order velocity structure function is significantly influenced by buoyancy forcing and large-scale inhomogeneities, while the mixed third-order structure function appearing in Yaglom equation exhibits a clear scaling exponent 1 in a small range of scales. The magnitudes of the various low $\ell$ degree anisotropic components of the equations are also estimated and are shown to be comparable to their isotropic counterparts at moderate to large scales. Finally, a qualitative analysis shows that the influence of buoyancy forcing at scales smaller than the Bolgiano scale is likely to remain important up to $\rayleigh=10^9$, thus preventing Kolmogorov scalings from showing up in convective flows at lower Rayleigh numbers.

We carry out numerical and mathematical investigations of shear Alfven waves inside of a spherical shell filled with an incompressible conducting fluid, and bathed in a strong dipolar magnetic field. We focus on axisymmetric toroidal and non-axisymmetric modes, in continuation of a previous work by Rincon & Rieutord (2003, A&A, 398, 663). Analytical expressions are obtained for toroidal eigenmodes and their corresponding frequencies at low diffusivities. These oscillations behave like magnetic shear layers, in which the magnetic poles play a key role, and hence become singular when diffusivities vanish. It is also demonstrated that non-axisymmetric modes are split into two categories, namely poloidal or toroidal types, following similar asymptotic behaviours as their axisymmetric counterparts when the diffusivities become arbitrarily small.

We present arguments which show that in all likelihood mesogranulation is not a true scale of solar convection but the combination of the effects of both highly energetic granules, which give birth to strong positive divergences (SPDs) among which we find exploders, and averaging effects of data processing. The important role played by SPDs in horizontal velocity fields appears in the spectra of these fields where the scale $\sim$4 Mm is most energetic; we illustrate the effect of averaging with a one-dimensional toy model which shows how two independent non-moving (but evolving) structures can be transformed into a single moving structure when time and space resolution are degraded. The role of SPDs in the formation of the photospheric network is shown by computing the advection of floating corks by the granular flow. The coincidence of the network bright points distribution and that of the corks is remarkable. We conclude with the possibility that supergranulation is not a proper scale of convection but the result of a large-scale instability of the granular flow, which manifests itself through a correlation of the flows generated by SPDs.