© ESO, 2014

1. Introduction

High resolution observations of the solar surface continue to provide new insights into the structure and evolution of small-scale magnetic fields in regions of quiet Sun. In any highly conducting fluid, overturning eddies tend to expel magnetic flux (Weiss 1966), so it is unsurprising that the magnetic field distribution at the solar surface is strongly influenced by the local convective motions. On the scale of the convectively-driven granulation (a typical granule has a width of approximately 1 Mm), the near-surface flow pattern naturally leads to the expulsion of vertical magnetic flux from the granular interiors, which then becomes concentrated into localised structures within the intergranular lanes (see, for example, Lin & Rimmele 1999). Larger scales of convective motion also play an important role in determining the spatial distribution of quiet Sun magnetic fields. The boundaries of supergranules, which have typical diameters of approximately 30 Mm (see e.g. Rieutord & Rincon 2010), are associated with a network of strong vertical magnetic flux concentrations (Simon & Leighton 1964). Furthermore, long-lived magnetic features in the internetwork regions tend to be advected preferentially towards the supergranular boundaries (de Wijn et al. 2008; Orozco Suárez et al. 2012). Between the granular and the supergranular scales, many studies have also suggested the presence of mesogranular convective cells (see, e.g., November et al. 1981; Muller et al. 1992; Shine et al. 2000; Domínguez Cerdeña 2003; Yelles Chaouche et al. 2011). Whether or not mesogranulation is really an intrinsic scale of convection in the near-surface layers of the Sun is still a matter of some debate (see, e.g., Rieutord et al. 2010; Rieutord & Rincon 2010; Yelles Chaouche et al. 2011; Katsukawa & Orozco Suárez 2012). Nevertheless, observations of magnetic fields do seem to indicate that magnetic flux concentrations accumulate preferentially at mesogranular boundaries (Domínguez Cerdeña 2003; de Wijn et al. 2005; Lites et al. 2008; Ishikawa & Tsuneta 2011; Yelles Chaouche et al. 2011).

In the quiet Sun, a near-surface magnetic structure with a characteristic field strength of approximately 400 G would have a magnetic energy density that is comparable in magnitude to the mean kinetic energy density of the surrounding non-magnetic convection (Galloway et al. 1977). It is perhaps surprising, therefore, that the peak vertical magnetic field strengths in regions of quiet Sun can often exceed a kG (see, for example, Stenflo 1973; Sánchez Almeida & Lites 2000; Domínguez Cerdeña et al. 2006a,b; Orozco Suárez et al. 2007). The existence of such strong fields is usually attributed to a phenomenon known as convective collapse. This magnetic field intensification process is driven by vertical convective motions along the field lines, which naturally lead to the formation of partially-evacuated magnetic flux concentrations (within which the magnetic pressure is high enough to be comparable to the surrounding gas pressure). This phenomenon has been studied extensively within the context of simplified models of thin flux tubes (see, for example, Webb & Roberts 1978; Spruit 1979) as well as in more realistic magnetoconvection simulations (Grossmann-Doerth et al. 1998; Bushby et al. 2008). Furthermore, recent Hinode observations seem to confirm that this is indeed the mechanism that is responsible for the production of kG-strength magnetic fields in regions of quiet Sun (Nagata et al. 2008). In addition to these locally intense vertical magnetic flux concentrations, regions of quiet Sun also contain significant quantities of predominantly horizontal magnetic flux (Orozco Suárez et al. 2007; Lites et al. 2008). Indeed, within the field of view that was investigated by Lites et al. (2008), the spatially averaged horizontal flux was found to be approximately five times larger than the corresponding (unsigned) mean vertical flux. Where the magnetic field is strongly inclined to the vertical (usually near the edges of granules), measurements indicate a typical field strength of the order of 100 G (Orozco Suárez et al. 2007). It is interesting to note that this value is comparable in magnitude to the estimated intrinsic field strength of the unresolved turbulent magnetic fields in these regions (Trujillo Bueno et al. 2004).

Several studies have shown that magnetoconvection models can produce a magnetic field distribution that is similar to that observed in the quiet Sun (Bushby & Houghton 2005; Khomenko et al. 2005; Stein & Nordlund 2006). However, it has been suggested that a near-surface dynamo, driven by the local convective motions, is responsible for producing a significant part of the observed small-scale magnetic activity, and recent observations seem to support this idea (see, e.g., Buehler et al. 2013). Although solar-like parameter regimes cannot yet be studied (due to computational constraints), numerical simulations of convectively-driven dynamo action in a Cartesian domain can be regarded as an idealised representation of the dynamo process in a small region of quiet Sun. Meneguzzi & Pouquet (1989) and Cattaneo (1999) established that non-rotating Boussinesq convection in an electrically-conducting fluid can act as an efficient dynamo provided that the magnetic Reynolds number is large. The resulting magnetic field distribution is highly disordered, although small-scale vertical magnetic flux concentrations do indeed accumulate within the convective downflows, as observed in the quiet Sun. A number of recent studies have investigated dynamo action in fully compressible convection. One approach to this problem is to simulate the dynamo in a model of quiet Sun convection which includes as many of the relevant physical processes as possible (such as radiative transfer and a realistic equation of state), and models of this type have produced results which compare favourably to observations (Abbett 2007; Vögler & Schüssler 2007; Schüssler & Vögler 2008; Pietarila Graham et al. 2009; Danilovic et al. 2010; Schüssler 2013). Other studies have considered dynamo action in much simpler models of compressible convection in a polytropic layer (or, in some cases, systems of more than one polytropic layers), including no non-essential physical effects (Käpylä et al. 2008; Brummell et al. 2010; Bushby et al. 2010, 2011, 2012). Minimal models of this form are well suited to parametric surveys. Furthermore, the findings of Moll et al. (2011) suggest that additional physical processes such as radiative transfer probably do not play a major role in the operation of the dynamo, so it is instructive to study the properties of these comparatively simple models.

Most previous simulations of dynamo action in compressible convection were carried out in relatively small computational domains, so are not able to investigate the effects of convective motions on scales that are larger than the basic granulation. In large domains, mesogranular cells have been observed in hydrodynamic Boussinesq convection (Cattaneo et al. 2001). Furthermore, when this flow was allowed to drive a dynamo, regions of vertical magnetic flux were found to accumulate preferentially at the corners of mesogranules (as observed in the quiet Sun). Simulations of fully compressible convection have produced rather mixed results in terms of mesogranulation. Some hydrodynamic calculations have produced mesogranulation with a clearly-defined horizontal spatial scale (see, for example, Rincon et al. 2005), whilst other studies have produced mesogranular cells with no intrinsic scale (Matloch et al. 2010). In a very recent study, Schüssler (2013) found a mesogranular pattern whose scale was determined (at least partially) by the growth rate of the associated dynamo. So there are at least some suggestions in the literature that mesogranulation might be strongly model-dependent. This is clearly an area in which further study is needed.

In a recent paper, Bushby et al. (2012) investigated the dynamo properties of compressible convection in a large aspect ratio domain. Like Rincon et al. (2005) they found clear evidence for mesogranulation in their hydrodynamic flow, with a well-defined peak in the kinetic energy spectrum at the mesogranular scale. Furthermore, the calculations of Bushby et al. (2012) suggest that the presence of mesogranulation is beneficial for dynamo action, at least in the kinematic regime (in comparable regions of parameter space, the dynamo is less efficient in domains that are too small to accommodate mesogranules). However, only a single flow was investigated in that study, so it is unclear how dependent these conclusions are upon the system parameters. In this paper, we carry out a systematic survey (subject to computational constraints) of the Reynolds number dependence of this system, focusing particularly upon the consequences for mesogranulation, before going on to consider the dynamo properties of the resultant flows. The model setup is described in the next section. We then discuss the hydrodynamic properties of the convective flows in Sect. 3, before moving on to describe the results from a series of dynamo calculations. Finally, our conclusions are presented in Sect. 4, where we discuss the relevance of our results to quiet Sun magnetism.

2. Model setup and governing equations

The model is an idealised representation of convection in the near-surface layers of a region of quiet Sun. We consider the evolution of a layer of electrically-conducting, compressible fluid, which is heated from below. The fluid properties are characterised by various parameters, including the thermal conductivity, K, the shear viscosity, μ, the magnetic diffusivity, η, the magnetic permeability, μ₀, and the specific heat capacities at constant density and pressure, c_v and c_p (respectively). All of these quantities are assumed to be constant throughout the fluid, which occupies a Cartesian domain of dimensions 0 ≤ x,y ≤ λd and 0 ≤ z ≤ d. The axes are defined in such a way that the z-axis points vertically downwards, parallel to the constant gravitational acceleration, g = gẑ, which implies that z = 0 corresponds to the upper surface of the domain. Periodic boundary conditions are enforced in each of the horizontal directions, whilst (in this idealised model) the upper and lower bounding surfaces are assumed to be impermeable and stress-free. Both the upper and lower boundaries are held at fixed, uniform temperatures (T₀ and T₀ + ΔT respectively, with ΔT > 0). We also assume that the horizontal components of the magnetic field vector vanish at z = 0 and z = d, which implies that the field is constrained to be vertical at the upper and lower boundaries.

Under these assumptions, the governing equations for the density ρ(x,t), the fluid velocity u(x,t), the temperature T(x,t) and magnetic field B(x,t) can be written in the following dimensionless form (see, for example, Bushby et al. 2008, for more details): where the pressure P(x,t) satisfies the equation of state for a perfect gas, P = ρT, the components of the rate of strain tensor, τ, are defined by (5)whilst the magnetic field satisfies ∇·B = 0. Note that all lengths have been scaled by d, which implies that the domain has unit depth in this dimensionless system, whilst d/(R_∗T₀)^1/2 has been used to scale time. The temperature has been scaled by T₀, whilst the density of the fluid has been scaled by ρ₀ (the density at the upper surface in the absence of any motion). Finally, the velocity field has been scaled by , whilst the magnetic field has been scaled by (μ₀R_∗ρ₀T₀)^1/2. In addition to the aspect ratio, λ, the evolution of the system is then determined by various dimensionless parameters. The ratio of specific heat capacities is defined by γ = c_p/c_v, whilst m = (gd/R_∗ΔT) − 1 is the polytropic index. The temperature difference across the domain is determined by the stratification parameter, θ = ΔT/T₀. Defining to be the dimensionless thermal diffusivity, the remaining parameters are the Prandtl number, σ = μc_p/K, and ζ₀ = ηc_pρ₀/K, which represents the ratio of the magnetic to the thermal diffusivities at the upper surface of the domain.

Given the complexity of the governing equations, we solve these numerically using an updated version of a well-tested code (Matthews et al. 1995). Horizontal derivatives are evaluated in Fourier space, whilst a fourth-order finite difference scheme is used to approximate the vertical derivatives. The equations are evolved in time using an explicit third-order Adams-Bashforth scheme with a variable timestep. A poloidal-toroidal decomposition is used to ensure that the magnetic field remains solenoidal. The code is parallelised using MPI, which enables us to investigate convectively-driven dynamos in large computational domains. In the absence of a magnetic field, the governing equations have an equilibrium solution, given by (6)which corresponds to a hydrostatic polytropic layer. All of the hydrodynamic simulations in this paper are initialised by perturbing this equilibrium state via the introduction of a small, random perturbation to the temperature distribution.

3. Numerical results

	Fig. 1 Temperature distribution in the horizontal plane defined by z = 0.1, for the λ = 10 cases: B1 (top left), B2 (top right), B3 (bottom left) and B4 (bottom right).
Open with DEXTER

3.1. Specifying the model parameters

Throughout this study, we fix the ratio of specific heat capacities to be γ = 5/3, which is the appropriate value for a monatomic gas. Fixing the polytropic index to be m = 1 implies that the layer is superadiabatically stratified, whilst a thermal stratification of θ = 3 implies that the temperature varies by a factor of 4 across the domain. We choose a fixed value of κ = 0.00548, whilst σ varies in the range 0.375 ≤ σ ≤ 1.5. These particular parameter choices imply that the mid-layer Rayleigh number (7)takes values in the range 2 × 10⁵ ≤ Ra ≤ 8 × 10⁵. In a large aspect ratio domain, the critical Rayleigh number for the onset of convection in the absence of a magnetic field is Ra_crit ≈ 900, so these parameter choices ensure that the resulting convection will be vigorous and time-dependent. The Reynolds number of the system is a measurable parameter which provides some indication of the level of turbulence within the flow. Using the depth of the domain (which, as we have already stated, equals unity in these dimensionless units) as the characteristic length-scale, we define the mid-layer Reynolds number to be (8)where U_rms is the total rms velocity, and ρ_mid is the mean density at the mid-layer of the domain. When considering the dynamo properties of this system, we shall often refer to the magnetic Reynolds number, which is defined by (9)whilst the (mid-layer) magnetic Prandtl number is defined to be (10)For a given flow, both Rm and Pm can be varied simultaneously by changing the value of ζ₀. Dynamo action is expected only in the high Rm (low ζ₀) regime. Finally, following Cattaneo (1999), we define the convective turnover time to be t_conv = 1/U_rms.

The hydrodynamic parameters for the simulation set that will be considered in this paper are defined in Table 1. In order to investigate the dependence of the flows upon the size of the computational domain, calculations have been carried out for λ = 4, λ = 10 and λ = 20. Although the flows in cases C1 and C2 are adequately resolved for these values of σ, it was not possible to carry out λ = 20 calculations at lower σ (and consequently higher Re) due to computational constraints. For a given value of σ, it is clear that Re and t_conv are only weakly dependent upon λ. Indeed, in this respect, the λ = 10 and λ = 20 cases are remarkably similar, whilst the small differences in the λ = 4 case suggest that the corresponding flow is weakly-constrained by the size of the domain in the λ = 4 calculations. Of course, this is a very crude comparison, so further analysis is needed in order to analyse the extent to which these flows are dependent upon the choice of aspect ratio.

3.2. Granulation and mesogranulation in hydrodynamic convection

The convective flows for the λ = 10 cases (B1 – B4) are illustrated in Fig. 1. All simulations produce highly time-dependent flows, with a near-surface pattern of convection that is comparable to the granular pattern that is observed at the surface of the Sun (with broad, warm upflows surrounded by a network of cooler, narrow downflows). In qualitative terms, it is clear that there tends to be more small-scale structure in the granular pattern at higher Reynolds numbers. Figure 2 shows the temperature distribution in two different horizontal planes for one of the λ = 20 cases (C2). In the z = 0.1 plot in the upper part of Fig. 2, the near-surface granulation pattern is qualitatively similar to that observed in the corresponding λ = 10 case (B2). However there is some evidence for the existence of structure on larger scales in the temperature distribution at the mid-plane. Only the strongest convective downflows from the surface layers penetrate to this depth – if we were to argue that these strong downflows correspond to the boundaries of mesogranules, then this mid-plane temperature distribution would provide some visual indication of the corresponding mesogranular scale. However, such an assertion clearly needs to be tested by carrying out a more quantitative analysis of these hydrodynamic flows.

	Fig. 2 Temperature distribution in two different horizontal planes for the λ = 20, σ = 1.0 case (C2). The upper plot shows the z = 0.1 plane, whilst the lower plot shows the z = 0.5 plane.
Open with DEXTER

To determine conclusively whether or not mesogranulation is present in these simulations, we calculated the time- and depth-averaged kinetic energy spectra, based upon the horizontal wavenumber, k_H. These kinetic energy spectra are shown in Fig. 3, where the wavenumber in the lower aspect ratio cases has been scaled for consistency with the corresponding λ = 20 spectra. In a previous paper (Bushby et al. 2012) we investigated the depth-dependence of the kinetic energy spectra for cases A2 and B2, concluding that signatures of mesogranulation (if it is present) can be seen at all depths. We therefore decided to focus upon depth-averaged spectra here to simplify the analysis. Various conclusions can be drawn from these kinetic energy spectra. Firstly, there is a clear “mesogranular” peak at low wavenumber (k_H ≈ 5–6) in all of the λ = 10 and λ = 20 cases. This corresponds to the horizontal scale that can be identified visually in the mid-plane temperature distribution that is illustrated in the lower part of Fig. 2. Furthermore, although it is present for all Reynolds numbers, the peak at this mesogranular scale seems to be more pronounced in the higher Re cases, which (unsurprisingly) also exhibit more structure at higher k_H. Regardless of the value of Re, it is clear that the flow is constrained by the size of the box in the λ = 4 cases, where the width of the domain is comparable to the mesogranular scale. There are no obvious indications of geometrical constraints in any of the other cases. Indeed, given the similarities between the comparable spectra in the λ = 10 and λ = 20 cases, we would argue that a λ = 10 domain should be large enough for the purposes of this particular study. We shall therefore not consider the dynamo properties of the λ = 20 simulations in this paper, deferring all further considerations of this case to a future study.

	Fig. 3 Kinetic energy spectra as a function of the horizontal wavenumber for the hydrodynamic convective flows that are defined in Table 1. Top: λ = 4; middle: λ = 10; bottom: λ = 20. In the upper two plots, the wavenumber has been normalised so that the scales are comparable in all cases.
Open with DEXTER

Mesogranulation is clearly a robust feature of this system, which raises questions about its origin and stability. Even though the peak in the kinetic energy energy spectrum is more pronounced in the higher Reynolds number cases, it is also intriguing to note that the horizontal scale of the mesogranular cells seems to be only rather weakly-dependent upon Re. Whilst we are not yet able to speculate about how this horizontal scale is determined, differences in the behaviour of the vertical vorticity across the various simulations may provide some clues regarding the origin of mesogranulation (and it has already been suggested that the vorticity distribution may be playing an important role in this regard, see e.g. Cattaneo et al. 2001). We show in Fig. 4 the time-averaged probability density function of the vertical component of the vorticity in the horizontal plane at z = 0.1, for simulations B2 and B4, whilst Fig. 5 shows a snapshot of the distribution of the vertical vorticity in the same horizontal plane for the same two cases. In the higher Re case (B4) rare events with very large absolute values of the vertical vorticity can be seen in the probability density function. The strongest localised vortices that can be seen in the lower part of Fig. 5 are coincident with the corners of the mesogranules (something that can easily be verified by using a passive tracer method, see e.g. Bushby et al. 2012), so we would argue that it is probable that these features play a role in the formation and maintenance of the mesogranular pattern. If this causal connection does indeed exist, the fact that these concentrated vortices are located at the edges of the mesogranules, where large-scale straining motions can amplify them, could explain why mesogranules persist for many turnover times. The weaker vortices that can be seen in the B2 case would also then be consistent with the fact that the mesogranulation is less pronounced in this lower Re case. It is worth noting here that the horizontal components of the vorticity seem to exhibit a much weaker dependence upon the Reynolds number, so the observed Re-dependence in the mesogranular pattern is more likely to be attributable to changes in the vertical vorticity profile. The dependence of the vorticity dynamics upon the aspect ratio also tends to support this hypothesis: presumably as a consequence of geometric constraints, wide tails in the probability density function for the vertical vorticity are only seen in the large λ cases in which mesogranulation is observed.

	Fig. 4 Time-averaged probability density function of the vertical vorticity close to the upper boundary (z = 0.1) for simulations B2 (dotted line) and B4 (solid line).
Open with DEXTER

	Fig. 5 Vertical vorticity distribution in the z = 0.1 horizontal plane. The upper plot corresponds to simulation B2, whilst the lower plot shows the results from simulation B4.
Open with DEXTER

3.3. Dynamo action in the kinematic regime

Having investigated the hydrodynamic properties of these flows, we now turn our attention to their dynamo properties. We start off each dynamo calculation by introducing a seed magnetic field into fully-developed hydrodynamic convection. This seed field is initially of low amplitude and is purely vertical, with a high wavenumber sinusoidal dependence upon both x and y. Whether the total magnetic energy grows or decays is then determined uniquely by the value of ζ₀ (or, equivalently, Rm), which can be specified as the field is introduced. For low Rm (high ζ₀), the evolution of the field is dominated by effects of ohmic dissipation, which implies that there is no dynamo. Dynamo action only becomes possible when the magnetic Reynolds number exceeds some critical value, Rm_crit. In the dynamo regime, the magnetic energy can exhibit strong fluctuations although, whilst the field is still weak, the overall trend is for exponential growth. In numerical simulations, this phase of kinematic growth can be prolonged indefinitely by artificially removing the Lorentz force terms from the system of equations. Here, we adopt this approach to determine how Rm_crit (and how the Rm-dependence of the dynamo growth rate) depends upon the other parameters in the system. Particularly close to Rm_crit, where there can be significant long-term fluctuations in the magnetic energy, it is necessary to evolve these calculations over many convective turnover times in order to obtain reliable estimates for the kinematic growth rates.

	Fig. 6 Growth rate of the magnetic energy in the kinematic dynamo regime. Top: Cases A1, A3 and A4. Bottom: Cases A2 and B2. In each plot, the growth rate has been scaled by the convective turnover time, tconv.
Open with DEXTER

Figure 6 shows the kinematic growth rates as a function of Rm for cases A1–A4 and B2. Note that the growth rates have been been scaled by the convective turnover time, t_conv, so as to ensure some degree of comparability between the different flows. Focusing initially upon the upper plot in Fig. 6, which shows the growth rates for three of the λ = 4 cases, we see that the flow with the lowest Reynolds number (A1) is the most efficient dynamo. At comparable values of Rm, the growth rates in the higher Re calculations are rather similar, although (in the absence of further calculations at higher Reynolds numbers) it would certainly be premature to suggest that we are approaching an asymptotic regime in which the behaviour is independent of Re. The lower plot in Fig. 6 shows the Rm-dependence of the growth rate for cases A2 and B2. These calculations are comparable in all respects except for the aspect ratio of the domain. At similar values of Rm the growth rates are systematically higher in the larger aspect ratio case, although the most dramatic differences are observed at low Rm, where (at comparable parameter values) the magnetic energy decays much more rapidly in the lower aspect ratio case. It is interesting to note that the growth rate appears to have a logarithmic dependence upon Rm in cases A1 and A2 (see also Bushby et al. 2011), whilst the growth rate curves in the other cases are better approximated by straight lines (at least for the range of values for Rm that has been considered). Logarithmic growth rate curves are not unheard of in the context of turbulent small-scale dynamos (Rogachevskii & Kleeorin 1997). However, direct comparison with previous numerical studies of dynamo action in compressible flows (e.g. Haugen et al. 2004) is complicated by the fact that, for computational convenience, we hold Re constant whilst varying Rm, which implies that the magnetic Prandtl number is not constant. Having said that, it is worth noting that our scaled growth rates at “high” Rm are comparable to those found by Haugen et al. (2004) in their study of dynamo action in forced (non-helical) turbulence.

	Fig. 7 Critical magnetic Reynolds number as a function of the Reynolds number (top) and the magnetic Prandtl number (bottom). The circles correspond to the λ = 4 cases, whilst the stars represent the λ = 10 calculations.
Open with DEXTER

Our estimates for Rm_crit as a function of the Reynolds number and the magnetic Prandtl number are shown in Fig. 7. The error bars in the upper plot give some indication of the level of precision of these estimates, which are bounded above and below by dynamo calculations which have positive and negative growth rates respectively. The main conclusion that can be drawn from the upper plot is that (at comparable values of Re) the critical magnetic Reynolds number is always smaller in the larger aspect ratio calculations. This indicates that it is always easier to drive a dynamo in the larger domain, where mesogranulation is present. Bushby et al. (2012) came to a similar conclusion although they only considered a single flow – the key point here is that this conclusion seems to be independent of the Reynolds number, which suggests that it is a robust feature of this system. The lower plot in Fig. 7 indicates that dynamo action is always easier to excite at higher values of the magnetic Prandtl number, with systematically higher values of Rm_crit in the lower Pm cases. The relative inefficiency of turbulent dynamos at low Pm is well known (Boldyrev & Cattaneo 2004; Schekochihin et al. 2005), so this result is unsurprising. It is worth noting here that estimates suggest that Pm ≈ 10^-6 in the solar photosphere (see, e.g., Ossendrijver 2003). Due to computational constraints, the smallest value of the magnetic Prandtl number for which we have found a dynamo is Pm ≈ 0.5, so we are clearly not in a position to say anything definitive about Rm_crit for the convection in the near-surface layers of the Sun. However, the local magnetic Reynolds number in this part of the Sun is very large (Rm ≈ 10⁶, see Ossendrijver 2003) so it would be surprising if a local dynamo was not operating in this region.

3.4. Filtered velocity field

	Fig. 8 Horizontal velocity in the x direction in the horizontal plane z = 0.1 for the λ = 10, σ = 1.0 case (B2). The upper plot shows the complete velocity field, the middle plot shows the large-scale velocity field (where short wavelengths are removed in Fourier space) and the lower plot shows the small-scale velocity field (where long wavelengths are removed from the flow).
Open with DEXTER

Fig. 9 Magnetic energy versus time for a kinematic dynamo simulation for the B2 case (ζ0 = 0.15). The solid line corresponds to the full dynamo problem in which the induction equation is solved using the complete velocity field. The dotted line corresponds to the dynamo that is driven by the small-scale flow (only modes such that kH > 10 are retained in the velocity field that is used to solve the induction equation), whilst the dash-dotted line corresponds to the case in which only the large-scale flow (where only modes such that kH < 10 are retained) has been used.

Open with DEXTER

These kinematic dynamo results show that the dynamo is more efficient in a larger domain. This suggests that the mesogranular component of the flow may be playing a crucial role in promoting dynamo action in the λ = 10 cases. After the mesogranular peak, the kinetic energy spectra that are shown in Fig. 3 decrease monotonically with increasing k_H. The fact that there is no secondary peak corresponding to a granular scale means that there is no clear separation in scales between the motions at large and small scales. Nevertheless, it is of interest to try to determine the ways in which components of the flow at different spatial scales influence the dynamo. In order to investigate this issue, we filter the velocity field in Fourier space. At each time-step within the numerical scheme, the velocity field is Fourier-transformed in the two horizontal directions. The filtered velocity field is obtained by multiplying the depth-dependent Fourier coefficients by the following mask function (11)where k_c is an arbitrary cut-off wavenumber. By suppressing the energy contained in wavenumbers larger than k_c, we isolate the large-scale component of the velocity field. Similarly, one can extract the small-scale component of the velocity field by ignoring wavenumbers smaller than k_c. In this approach, the full velocity field is then used in all of the governing equations except in the induction equation, Eq. (4), where the filtered velocity field is used to determine the evolution of the magnetic field. A similar filtering method has already been used by Hughes & Proctor (2013) to study the effect of shear on convectively-driven dynamos and by Tobias & Cattaneo (2008) to study fast dynamos. This simple filtering process does not preserve coherent structures in physical space, since the mask function only depends on the wavenumber and not on the phase of the Fourier coefficients. However, despite its limitations, this approach should give us some insights into the role of mesogranulation. In the following, we use the same parameters as for simulation B2 in Table 1. We choose a cut-off wavenumber of k_c = 10 which lies to the right of the mesogranular peak in the kinetic energy spectrum that is shown in Fig. 3.

Fig. 10 Snapshot of a statistically-steady nonlinear dynamo calculation (corresponding to case B4, with ζ0 = 0.15 – see Table 2). Top: the vertical magnetic field distribution in the horizontal planes defined by z = 0.1 (left) and z = 0.5 (right), where the colour table has been saturated at Bz =  ± Beq. Middle: as above, but here the contour levels have been saturated at ± 0.25Beq to highlight some of the small-scale structure in the field distribution. Bottom: the corresponding temperature distributions for z = 0.1 (left) and z = 0.5 (right).

Open with DEXTER

The effects of this filtering process are illustrated in Fig. 8, which shows snapshots of the full and filtered velocity fields for this simulation, whilst Fig. 9 shows the time-evolution of the magnetic energy for the corresponding kinematic dynamo calculations (for ζ₀ = 0.15). Although the magnetic energy grows at a slightly slower rate for the large-scale flow than it does in the unfiltered case, the dynamo efficiencies are fairly comparable. However, the small-scale flow is only a marginal kinematic dynamo for these parameter values. This disparity between the dynamo growth rates for the large-scale and the small-scale flows cannot be attributed to there being more kinetic energy in the large-scale motions. In fact the rms velocity associated with the large-scale flow is approximately the same as for the small-scale flow (U_rms ≈ 0.25), which implies that the effective magnetic Reynolds number is approximately the same for both filtered flows. Therefore these results do seem to confirm the notion that the flows at scales that are comparable to the mesogranular scale are more effective than the small-scale flows when it comes to driving a kinematic dynamo in this particular system.

Fig. 11 Time-averaged probability density function for the vertical component of the magnetic field at 4 different depths for the statistically-steady, nonlinear dynamo calculation that is illustrated in Fig. 10. The different curves correspond to z = 0 (solid line), z = 0.1 (dash-dotted line), z = 0.5 (dashed line) and z = 0.9 (dotted line). At each depth, the magnetic field has been normalised by the local equipartition field strength, Beq.

Open with DEXTER

3.5. Nonlinear results and comparison to observations

Having established the conditions under which this system can excite a kinematic dynamo, we now return to the full nonlinear problem in which the magnetic field is allowed to exert a dynamical influence upon the flow. One particular nonlinear calculation is illustrated in Fig. 10, which shows a snapshot of the temperature and magnetic field distribution in two different horizontal planes after the dynamo has evolved into a statistically-steady state. This dynamo is based upon case B4, which is the highest Reynolds number λ = 10 simulation that has been considered in this paper, whilst a value of ζ₀ = 0.15 was chosen when the initial seed magnetic field was introduced into the flow. In this nonlinear state, the rms velocity (U_rms = 0.390) is approximately 5% smaller than that of the original hydrodynamic flow. This implies that the corresponding Reynolds numbers are Re = 458 and Rm = 475 (with Pm = 1.04). The near-surface temperature distribution that is shown in Fig. 10 is qualitatively similar to the one illustrated in Fig. 1, whilst the mid-layer temperature distribution is again indicative of the presence of mesogranulation. The fact that the temperature distribution is rather similar to that of the corresponding hydrodynamic flow suggests that the magnetic field only has a local (rather than global) effect upon the convective dynamics. The contour levels in the corresponding magnetic field plots are related to the time-averaged equipartition field strength, , at that particular depth (where the overbar corresponds to a horizontal spatial average). In the upper plots in Fig. 10, the colour table has been saturated at ± B_eq, whilst the colour table in the middle plots has been saturated at ± 0.25B_eq. Comparing the near-surface magnetic field (in particular) with the mid-layer temperature distribution, it is evident that the stronger vertical magnetic fields tend to be associated with the mesogranular boundaries, which correspond to the locations of the strongest convective downflows. This is, therefore, at least qualitatively consistent with observations of quiet Sun magnetic fields. Near the upper surface, the intergranular lanes tend to be associated with weaker magnetic fields, whilst very weak fields can also be seen in some of the granular interiors. These presumably correspond to loops of field that are advected into the surface layers by the convective upflows. As illustrated in Fig. 10, the weak magnetic field is almost “space filling” at the mid-plane of the domain, where a complex, mixed-polarity field distribution can be observed.

Still focusing upon the calculation that is illustrated in Fig. 10, Fig. 11 shows the corresponding (time-averaged) probability density functions (PDFs) for the vertical component of the magnetic field at four different depths within the domain. At each depth, the magnetic field strength has again been normalised by the local equipartition magnetic field strength. Near the surface, the PDFs take the form of stretched exponential distributions, with the tail extending well into the “super-equipartition” range (with a peak field that is roughly three times the equipartition value). At lower values of Re, similar PDFs can be found in nonlinear dynamo calculations at comparable values of Rm (Bushby et al. 2012) so this tendency to produce super-equipartition fields does not appear to be strongly dependent upon either Re or Pm (at least for computationally accessible values of these parameters). To relate these peak field strengths to solar observations, we need to bear in mind that B_eq ≈ 400 G near the solar surface (Galloway et al. 1977). The strongest vertical magnetic fields in these simulations would therefore correspond to the kG-strength fields that are observed in the quiet Sun. As discussed in the Introduction, the mechanism for the formation of these strong fields is related to the well-known convective collapse instability, and the resultant concentrations of vertical magnetic flux are partially evacuated. For this calculation, the minimum surface density within one of these magnetic regions is (typically) of order of 30% of the mean surface density, although it can be substantially smaller than this. Looking again at Fig. 11, it is clear that the PDFs at z = 0 and z = 0.1 are rather similar, although slightly stronger peak fields (relative to B_eq) are found at z = 0.1. Even at z = 0.5 the peak field strength can significantly exceed the equipartition value. In the deeper layers, the PDF is well-approximated by a standard exponential distribution in which the peak field strength is comparable to B_eq. Note that all of these PDFs are almost symmetric about B_z/B_eq = 0. This implies a mixed polarity field distribution at all depths with no preference for either polarity. This result reflects the absence of any imposed magnetic fields in these calculations.

When comparing these calculations to observations of magnetic fields in the quiet Sun, it is also of interest to measure the mean (unsigned) vertical flux across the surface of the domain, . In the study of Lites et al. (2008), they measured the mean value of the modulus of the line of sight component of B to be approximately 11G, which is slightly less than 3% of the local equipartition field strength. Defining to be the surface value of the mean of | B_z | and B_eq,0 to be the corresponding equipartition field strength, this simulation gives a time-averaged value of , which is remarkably similar to the value quoted by Lites et al. (2008). Given the idealised nature of this model, and the fact that we are comparing numerical results to lower resolution observations, we should not read too much into the fact that we have such close agreement in this case. However, it is clearly encouraging that the mean value of | B_z | is of the correct order of magnitude to be consistent with observations. For completeness, it is worth noting that the quantity is actually an increasing function of depth until just below the mid-layer (reaching a peak value of approximately 0.075) before decreasing again with depth in the lower part of the domain. The time-average of the ratio of the total magnetic energy (ME) to the total kinetic energy (KE) is of the same order of magnitude as (here ME/KE = 0.030).

One of the main aims of this study was to assess the Re dependence of these convectively-driven dynamos. We therefore carried out a range of nonlinear calculations (which are summarised in Table 2). The presence of partially-evacuated regions (which, for various reasons, reduces the critical time-step for the stability of the numerical scheme) greatly increases the computational demands relative to the corresponding kinematic calculations. We therefore carried out a more restricted parameter survey for the full nonlinear problem. However, it is still possible to identify clear trends in the parameter-dependence of the statistics that were measured. In the λ = 4 calculations, we chose values of ζ₀ that produced dynamos with similar values of Rm. In that case, both and ME/KE increase with increasing Pm (or, equivalently, decreasing Re). So, at comparable values of Rm, the nonlinear efficiency of the dynamo is greatest in the higher Pm/lower Re cases. This can be attributed to the fact that Rm_crit is lower in the higher Pm cases, which implies that the dynamos in these cases are further above onset than those in the lower Pm regime. For computational convenience, we adopted a slightly different approach in the more computationally expensive λ = 10 cases in that we fixed the input parameter ζ₀ = 0.15 (as opposed to adjusting it to attain a dynamo at a certain value of Rm). Although in all of these cases, there is again a tendency for larger ME/KE ratios in the higher Pm cases. The two dynamo simulations corresponding to A4 and B4 are comparable in every respect except for the differing aspect ratios. The higher level of efficiency in the larger aspect ratio B4 case is again simply a consequence of the fact that it is further from onset than the corresponding A4 calculation.

4. Summary and discussion

Motivated by observations of convection and small-scale magnetic fields in regions of quiet Sun, we have investigated the dynamo properties of compressible convection in a large aspect ratio domain. From a purely hydrodynamic perspective, the first key result from this study is that mesogranulation seems to be a robust feature of models of this type, although a wide domain is needed in order to accommodate these large-scale features. Furthermore, over the range of computationally accessible values of the Reynolds number, we have found that the peak at the mesogranular scale in the kinetic energy spectrum is more pronounced in the higher Re cases. All of the convective flows that were considered in this paper are capable of driving a dynamo, provided that the magnetic Reynolds number exceeds some critical value. This critical value tends to increase with increasing values of Re. However Rm_crit is always smaller in the λ = 10 cases than it is in comparable calculations in the λ = 4 domain, and for a given value of Rm > Rm_crit the kinematic growth rates are always higher in the λ = 10 case. In addition to this aspect ratio dependence, our analysis of the dynamo properties of velocity fields that have been filtered in Fourier space also seems to confirm the suggestion of Bushby et al. (2012) that mesogranulation is beneficial for dynamo action. In the nonlinear regime, near-surface vertical magnetic flux concentrations tend to accumulate preferentially at the boundaries of mesogranules, where peak field strengths can significantly exceed the equipartition field strength (an observation that is consistent with the kG-strength fields that are observed in the quiet Sun). The horizontally-averaged values of | B_z | at the upper surface of the domain are also of the correct order of magnitude to be comparable to observed values. In general, the nonlinear efficiency of the dynamo (the ratio of the total magnetic energy to the total kinetic energy) increases with increasing values of the magnetic Prandtl number, although the magnetic energy never exceeds 10% of the kinetic energy in the parameter regimes that are investigated in this paper.

It should be stressed again that it is not possible to carry out simulations of this type in realistic solar-like parameter regimes, so some care is needed when relating results to observations. However, our findings do tend to favour the existence of a genuine mesogranular scale of motion in the near-surface convective flows. As we have noted, the Re-dependence of the vertical vorticity distribution provides some support for the hypothesis that there is a relationship between these vortices and the mesogranulation, although this causal connection is still rather speculative. Nevertheless, the fact that mesogranulation appears to be absent in rotating turbulent convection (in which vortices tend to align themselves with the rotation axis, see, e.g., Brummell et al. 1998), provides some support for the idea that the vortices may be playing an important role in this context. Another issue that needs to addressed here is the question of why some simulations naturally produce mesogranulation whilst others do not. Some numerical schemes rely upon artificial viscosities, which can have a strong dissipative effect at small scales. If this strongly influences the vorticity distribution, which is structured on small scales, this may explain why calculations of this type do not tend to produce mesogranulation with a well-defined spatial scale. However, this is a hypothesis that needs to be tested, so this is an obvious area for future work.

Having established that the convective flows in this idealised model can produce dynamo-generated magnetic fields that compare favourably to observations, the next step will be to relax some of the simplifying assumptions that have been made. Whilst the assumption of periodic boundary conditions in the horizontal direction is perfectly reasonably for a local model of this type, we should not rule out the possibility that the upper and lower boundary conditions may be playing an important role in determining the behaviour of the system. Some calculations assume open upper and lower boundaries, and we have carried out some test calculations for an analogous system. Although results are still at a rather preliminary stage, they do seem to indicate the presence of mesogranulation, which suggests that this phenomenon is not strongly-dependent upon the choice of hydrodynamic boundary conditions. This system also allows us to relax some of the constraints due to the magnetic boundary conditions. An initial dynamo calculation suggests that the near-surface vertical magnetic field distribution is qualitatively similar to that observed in the more idealised calculations that are described in this paper (despite the presence of additional horizontal magnetic fields at the same level). This preliminary calculation therefore indicates that the operation of the dynamo may not be strongly-dependent upon the choice of magnetic field boundary conditions, although further work needs to be done in order to test this hypothesis. Finally, it is worth noting that the convective domain that is described in this paper is relatively shallow in terms of pressure scale heights. Increasing the number of scale heights in the system may have an effect upon the efficiency of the dynamo, but again further study is needed in order to investigate this issue.

Acknowledgments

All calculations were carried out on the UK MHD Consortium computing facilities at St Andrews and Warwick. B. Favier thanks the Cambridge Newton Trust for financial support.

References

All Tables

All Figures

	Fig. 1 Temperature distribution in the horizontal plane defined by z = 0.1, for the λ = 10 cases: B1 (top left), B2 (top right), B3 (bottom left) and B4 (bottom right).
Open with DEXTER
In the text

	Fig. 2 Temperature distribution in two different horizontal planes for the λ = 20, σ = 1.0 case (C2). The upper plot shows the z = 0.1 plane, whilst the lower plot shows the z = 0.5 plane.
Open with DEXTER
In the text

	Fig. 3 Kinetic energy spectra as a function of the horizontal wavenumber for the hydrodynamic convective flows that are defined in Table 1. Top: λ = 4; middle: λ = 10; bottom: λ = 20. In the upper two plots, the wavenumber has been normalised so that the scales are comparable in all cases.
Open with DEXTER
In the text

	Fig. 4 Time-averaged probability density function of the vertical vorticity close to the upper boundary (z = 0.1) for simulations B2 (dotted line) and B4 (solid line).
Open with DEXTER
In the text

	Fig. 5 Vertical vorticity distribution in the z = 0.1 horizontal plane. The upper plot corresponds to simulation B2, whilst the lower plot shows the results from simulation B4.
Open with DEXTER
In the text

	Fig. 6 Growth rate of the magnetic energy in the kinematic dynamo regime. Top: Cases A1, A3 and A4. Bottom: Cases A2 and B2. In each plot, the growth rate has been scaled by the convective turnover time, tconv.
Open with DEXTER
In the text

	Fig. 7 Critical magnetic Reynolds number as a function of the Reynolds number (top) and the magnetic Prandtl number (bottom). The circles correspond to the λ = 4 cases, whilst the stars represent the λ = 10 calculations.
Open with DEXTER
In the text

	Fig. 8 Horizontal velocity in the x direction in the horizontal plane z = 0.1 for the λ = 10, σ = 1.0 case (B2). The upper plot shows the complete velocity field, the middle plot shows the large-scale velocity field (where short wavelengths are removed in Fourier space) and the lower plot shows the small-scale velocity field (where long wavelengths are removed from the flow).
Open with DEXTER
In the text

	Fig. 9 Magnetic energy versus time for a kinematic dynamo simulation for the B2 case (ζ0 = 0.15). The solid line corresponds to the full dynamo problem in which the induction equation is solved using the complete velocity field. The dotted line corresponds to the dynamo that is driven by the small-scale flow (only modes such that kH > 10 are retained in the velocity field that is used to solve the induction equation), whilst the dash-dotted line corresponds to the case in which only the large-scale flow (where only modes such that kH < 10 are retained) has been used.
Open with DEXTER
In the text

	Fig. 10 Snapshot of a statistically-steady nonlinear dynamo calculation (corresponding to case B4, with ζ0 = 0.15 – see Table 2). Top: the vertical magnetic field distribution in the horizontal planes defined by z = 0.1 (left) and z = 0.5 (right), where the colour table has been saturated at Bz =  ± Beq. Middle: as above, but here the contour levels have been saturated at ± 0.25Beq to highlight some of the small-scale structure in the field distribution. Bottom: the corresponding temperature distributions for z = 0.1 (left) and z = 0.5 (right).
Open with DEXTER
In the text

	Fig. 11 Time-averaged probability density function for the vertical component of the magnetic field at 4 different depths for the statistically-steady, nonlinear dynamo calculation that is illustrated in Fig. 10. The different curves correspond to z = 0 (solid line), z = 0.1 (dash-dotted line), z = 0.5 (dashed line) and z = 0.9 (dotted line). At each depth, the magnetic field has been normalised by the local equipartition field strength, Beq.
Open with DEXTER
In the text