This is a repository copy of Jets and large-scale vortices in rotating Rayleigh-Bénard convection.
|
|
- Annice Glenn
- 5 years ago
- Views:
Transcription
1 This is a repository copy of Jets and large-scale vortices in rotating Rayleigh-Bénard convection. White Rose Research Online URL for this paper: Version: Accepted Version Article: Guervilly, C and Hughes, DW orcid.org/ (2017) Jets and large-scale vortices in rotating Rayleigh-Bénard convection. Physical Review Fluids, 2 (11) ISSN X American Physical Society. This is an author produced version of a paper published in Physical Review Fluids. Uploaded in accordance with the publisher's self-archiving policy. Reuse Unless indicated otherwise, fulltext items are protected by copyright with all rights reserved. The copyright exception in section 29 of the Copyright, Designs and Patents Act 1988 allows the making of a single copy solely for the purpose of non-commercial research or private study within the limits of fair dealing. The publisher or other rights-holder may allow further reproduction and re-use of this version - refer to the White Rose Research Online record for this item. Where records identify the publisher as the copyright holder, users can verify any specific terms of use on the publisher s website. Takedown If you consider content in White Rose Research Online to be in breach of UK law, please notify us by ing eprints@whiterose.ac.uk including the URL of the record and the reason for the withdrawal request. eprints@whiterose.ac.uk
2 Jets and large-scale vortices in rotating Rayleigh-Bénard convection Céline Guervilly 1 & David W. Hughes 2 1 School of Mathematics, Statistics and Physics, Newcastle University, Newcastle Upon Tyne, NE1 7RU, UK 2 Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK November 10, 2017 Abstract One of the most prominent dynamical features of turbulent, rapidly-rotating convection is the formation of large-scale coherent structures, driven by Reynolds stresses resulting from the small-scale convective flows. In spherical geometry, such structures consist of intense zonal flows that are invariant along the rotation axis. In planar geometry, long-lived, depthinvariant structures also form at large scales, but, in the absence of horizontal anisotropy, they consist of vortices that grow to the domain size. In this work, through the introduction of horizontal anisotropy into a numerical model of planar rotating convection by the adoption of unequal horizontal box sizes (i.e. L x L y, where the xy-plane is horizontal), we investigate whether unidirectional flows and large-scale vortices can coexist. We find that only a small degree of anisotropy is required to bring about a transition from dynamics dominated by persistent large-scale vortices to dynamics dominated by persistent unidirectional flows parallel to the shortest horizontal direction. When the anisotropy is sufficiently large, the unidirectional flow consists of multiple jets, generated on a timescale smaller than a global viscous timescale, thus signifying that the upscale energy transfer does not spontaneously feed the largest available mode in the system. That said, the multiple jets merge on much longer timescales. Large-scale vortices of size comparable with L x systematically form in the flanks of the jets and can be persistent or intermittent. This indicates that large-scale vortices, either coexisting with jets or not, are a robust dynamical feature of planar rotating convection. 1 Introduction In the liquid core of planets, heat is transported primarily by convection. Much effort has been dedicated to the study of the properties of convection, in particular because of its crucial implications for the thermal evolution of planets and the generation of magnetic fields (e.g. Jones, 2007; Aurnou et al., 2015). Convection is strongly affected by the rapid rotation of the planet via the action of the Coriolis force, which tends to inhibit small-scale motions along the direction of the rotation axis. Owing to the very low fluid viscosity, flows in planetary cores are turbulent, although the nonlinear inertial effects are relatively weak compared with the Coriolis force. Rotationally-constrained convection in planetary cores is thus characterised by small Ekman numbers (Ek = ν/2ωd 2, where ν is the fluid kinematic viscosity, Ω the rotation rate and d is the size of the system), large Reynolds numbers (Re = Ud/ν, where U is a typical flow velocity at the lengthscale d) and small Rossby numbers (Ro = EkRe = U/2Ωd). Under these conditions, the convective flows take the form of tall narrow columns aligned with the rotation axis. The concerted action of the Reynolds stresses resulting from the small-scale convective columns drives large-scale coherent structures, which can become prominent dynamical features. The occurence of these large-scale structures in experimental and numerical models of turbulent flows (e.g. Heimpel et al., 2005; Read et al., 2015) attracts much interest because they are 1
3 observed in many geophysical and astrophysical bodies (e.g. Porco et al., 2003; Dyudina et al., 2008). Although the large-scale flows are usually directed perpendicularly to the direction of gravity, they can nonetheless affect the convective heat transport by shearing the convective cells or by modifying the local properties of the convection (e.g. Guervilly et al., 2014; Yadav et al., 2016). Understanding the formation and long term evolution of the large-scale flows is therefore an essential part of the description of a convective system. Furthermore, largescale flows can promote the generation of coherent magnetic fields at the system size while the underlying small-scale turbulent flows may, of themselves, be able to generate only fluctuating small-scale magnetic fields (Guervilly et al., 2015, 2017). In spherical rotating convection, the large-scale flows consist of intense azimuthal and axisymmetric (i.e. zonal) flows that are invariant along the rotation axis (e.g. Gilman, 1977). The presence of the curved spherical boundaries is crucial for the production of such zonal flows (Busse & Hood, 1982). In planar geometry in a horizontally periodic domain, long-lived, depthinvariant structures also form at large scales when the horizontal box sizes are equal. However, the absence of horizontal anisotropy in this case ensures that no persistent unidirectional flow emerges spontaneously; the large-scale structures therefore consist of vortices that grow to the box size (e.g. Chan, 2007). These large-scale vortices (LSVs) form whether the rotation axis is aligned or inclined with respect to the direction of gravity (Chan & Mayr, 2013). In this work, we are interested in the possible relation between unidirectional flows and LSVs and whether these two types of large-scale coherent structures can coexist. In simulations of rotating spherical convection, the formation of tall vortices at scales larger than the typical convective size has not been observed. In a Cartesian domain, the size of the LSVs is limited only by the horizontal extent of the periodic domain. If the upscale energy transfer were allowed to continue to larger scales, the LSV would eventually feel either the latitudinal variation of the Coriolis parameter (the β-effect) or the slope of the vertical boundaries. The absence of LSVs in spherical models suggests that the large-scale dynamics prefers unidirectional flows rather than vortices when possible. The convective flows in a turbulent rotating system are three-dimensional (3D), but strongly anisotropic owing to the rapid rotation. Results from two-dimensional (2D) turbulence might thus be helpful in explaining or predicting some aspects of the dynamics of convective systems on timescales longer than the rotation period when Ek 1, Ro 1. In particular, the LSVs obtained in 3D convective systems bear a strong resemblance to the condensates that form in 2D turbulence from an upscale kinetic energy transfer (Rubio et al., 2014). However, the analogy with 2D turbulence remains only partial in this particular case because the convective flows that feed energy into the LSVs are anisotropic but certainly not depth-invariant (Guervilly et al., 2014). The idea that unidirectional flows are preferred to vortices on the large scales once the horizontal symmetry is broken is supported by studies of forced 2D turbulence in which the gradual increase of the β-effect leads first to the weakening of the coherent vortices, and eventually to their disappearance, while unidirectional flows emerge (Maltrud & Vallis, 1991). However, a counter-example to the weakening of vortices in the presence of symmetry breaking is demonstrated by the recent study of Frishman et al. (2017), who considered a forced 2D turbulence model with anisotropy introduced by using a rectangular (non-square) periodic domain; with this geometry, they found that both unidirectional flows and vortices emerge at large scales and coexist on long timescales. Here we study whether this remarkable result can carry over to three-dimensional flows. The two main objectives of the present paper are to determine whether the LSVs that form in rotating planar convection can persist when the symmetry between the horizontal directions is broken and also to determine whether unidirectional flows and LSVs can coexist on long timescales. There are a number of ways of breaking the horizontal symmetry. Considering convection with a tilted rotation axis, thereby modelling mid latitudes, is an obvious physically motivated approach. Here though, we choose to introduce a distinction between the horizontal 2
4 directions in a three-dimensional Cartesian model of rotating Boussinesq convection by considering rectangular, as opposed to square, horizontal domains. Horizontal anisotropy is thus imposed on the system by allowing larger lengthscales to exist in one horizontal direction than the other. This formulation is computationally straightforward to implement, thus allowing a thorough exploration of parameter space, and also allows us to make comparisons with the forced 2D turbulence study of Frishman et al. (2017). The layout of the paper is as follows. The mathematical formulation of the problem is briefly described in 2. The formation of unidirectional flows and LSVs is described in 3.1. The box aspect ratio required for the transition from LSVs to unidirectional flows is quantified in 3.2 and the selection of the lengthscale of the large-scale flow is described in 3.3. The occurence of bistable states is presented in 3.4. In 3.5, we establish the domain of existence of the large-scale flows. Finally, in 3.6, we discuss the effect of the large-scale flows on the convective heat transfer. A conclusion is given in 4. 2 Mathematical formulation We use a local planar model of rotating Boussinesq convection, as in(guervilly et al., 2014). The computational domain is three-dimensional and periodic in the horizontal directions. A vertical temperature difference, T, is imposed across the layer of depth d. In the horizontal plane, the size of the domain is L x d in the x direction and L y d in the y direction. The gravitational field is uniform, g = ge z. The rotation vector is Ωe z. The fluid has kinematic viscosity ν, thermal diffusivity κ, density ρ and thermal expansion coefficient α, all of which are constant. Lengths are scaled with d, times with 1/(2Ω) and temperature with T. The system of dimensionless governing equations thus becomes u t +u u+e z u = p+ek 2 u+ RaEk2 θe z, Pr (1) u = 0, (2) θ t +u θ u z = Ek Pr 2 θ, (3) where u = (u x,u y,u z ) is the velocity field, p the pressure, and θ the temperature perturbation relative to a linear background profile. The dimensionless parameters are the Rayleigh number, Ra = αg Td 3 /κν, the Ekman number, Ek = ν/2ωd 2, and the Prandtl number, Pr = ν/κ. The upper and lower boundaries are taken to be perfect thermal conductors, impermeable and stress free. Stress-free boundary conditions, where boundary friction is absent, favour the emergence of large-scale flows of large amplitude in both planar (Stellmach et al., 2014; Kunnen et al., 2016) and spherical geometry (e.g. Yadav et al., 2016). Equations (1) (3) are solved using the pseudospectral code described in detail in (Cattaneo et al., 2003). 3 Results 3.1 Jets and large-scale vortices For most of the results discussed in the paper, the parameters are fixed as Ek = 10 5, Ra = , Pr = 1, L x = 1, with variations only in L y. The Ekman and Rayleigh numbers are varied only in 3.3 and 3.4, in order to integrate the system over many viscous timescales, and in 3.5, in order to determine the domain of existence of the large-scale flows; furthermore, we perform one run with L x = 2 and L y = 4, which we discuss in the conclusion. The value of the Rayleigh number Ra = corresponds to 7.4Ra c, where Ra c is the critical Rayleigh number at the linear onset of convection (Chandrasekhar, 1961). 3
5 L y Re x Re y Re z Re Nu Γ α n N x N y N z Table 1: Simulations run at varying L y for Ek = 10 5, Ra = and L x = 1. The value n denotesthenumberofpositive(or, equivalently, negative)jetsofthemeanvelocityu x. Thefinal column lists the numerical resolutions employed, with N i indicating the number of collocation points in the direction i ˆRe x ˆRey ˆRez time Figure 1: Time series of the r.m.s. values of the three velocity components for L y = 8. The grey dashed lines indicate the times of the snapshots in figure 2. Table 1 contains some global quantities obtained from simulations with L y in the range 1 L y 8, together with the numerical resolutions used. The root mean square (r.m.s.) amplitudeofthevelocityismeasuredbythereynoldsnumberre = u 1/2 Ek 1, andther.m.s. amplitude of the i-component of the velocity by Re i = u 2 i 1/2 Ek 1, where the angle brackets denote volume and time averages. The inclusion of the Ekman number in these definitions is due to our choice of dimensional units. In order to confirm that the resolution was adequate, and also to verify that the results were not critically dependent on the number of spectral modes adopted, provided this was sufficiently large, we simulated the case of L y = 3 with both 512 and 768 modes in the y direction. As can be seen from table 1, the differences between the two simulations are small Initial growth of the convective and large-scale flows Unless stated otherwise, the simulations described were started from a small perturbation to the basic state of rest. Figure 1 shows the time series of ˆRe x, ˆRe y, and ˆRe z, which are defined as 4
6 Figure 2: Snapshots of the depth-averaged axial vorticity taken at the times indicated by the grey dashed lines in figure 1. the Reynolds number Re i except with only a spatial average, for L y = 8. The figure shows the evolution only for dimensionless time units, i.e of the global viscous timescale; the full simulation was actually run for one viscous timescale. We first describe the initial growing phase of the kinetic energy; the saturated phase is detailed in the following subsections. Initially, the kinetic energy grows exponentially, with the velocity characterised by small horizontal scales. This short period of rapid growth in the kinetic energy is followed by an abrupt fall once the nonlinear terms become important. At this stage the convection remains small scale with a horizontal scale of approximately 0.1 for the chosen parameters. Subsequently the kinetic energy grows again, this time on a much longer timescale, as the flow moves to large horizontal scales. Figure 2 shows snapshots of the depth-averaged axial vorticity (denoted by ω z z, where i indicates an average along the direction i) at the different times indicated in figure 1. This representation focuses on the formation of large-scale depth-invariant flows; the convective flows, whose axial vorticity is mainly anti-symmetric with respect to the mid-plane (Chandrasekhar, 1961), are not visible on these plots. During the period of slow evolution to large-scale flows, vortices of large horizontal scale compared with the convective scale first become established across the domain (figure 2a), with movement of the vortices in both horizontal directions. The growth of these vortices corresponds to an increase in both ˆRe x and ˆRe y, with ˆRe z saturating at a lower value; these large-scale vortices consist essentially of depth-invariant horizontal flows. Vortex mergers occur and the LSVs grow to a scale comparable with the shortest horizontal box size (figure 2b). Subsequently, a number of jets become established, parallel to the shortest horizontal direction (figure 2c). The LSVs persist in the presence of the jets and, in the saturated state, move only with the jets. The growth of kinetic energy is eventually saturated on a timescale that is long compared with the dynamical timescale of the convection but short compared with the viscous timescale. 5
7 Figure 3: Snapshots of the horizontal and vertical cross-sections of (a) the axial vorticity and (b) the vertical velocity for L y = Flow in the saturated phase for L y = 1 We now describe the flow in the final saturated state, first recalling results obtained for the standard case of L y = 1. Figure 3 shows snapshots of horizontal and vertical cross-sections of ω z andu z. Thedominantfeatureintheplotofω z isalarge-scalevortexthatispredominantlyzinvariant. The LSV consists of a concentrated cyclone and a more dilute anticyclone, which grow until they reach the box size. LSVs are a robust feature of turbulent, rotationally-constrained convection in Cartesian geometry when L x = L y (Chan, 2007; Käpylä et al., 2011; Favier et al., 2014; Guervilly et al., 2014; Rubio et al., 2014; Stellmach et al., 2014). The asymmetry between cyclone and anticyclone is thought to be due both to the instability of the large-scale anticyclone, whose local vorticity is comparable with the planetary vorticity 2Ω, and to the preference for cyclonic convective vortices (Guervilly et al., 2014). This asymmetry eventually disappears as therossbynumberisloweredviaadecreaseoftheekmannumber(julienet al.,2012;stellmach et al., 2014). In spectral space, the LSV is a superposition of the wavenumbers (k x,k y ) = (0,1) and (k x,k y ) = (1,1)forthevelocitycomponentu x andof(k x,k y ) = (1,0)and(k x,k y ) = (1,1)foru y. The LSV therefore has a visible signature in snapshots of u x xz and u y yz. However, the LSV drifts in random horizontal directions over time, so long time-averages of either u x xz or u y yz produce a zero profile. In this sense, the LSVs do not produce any persistent unidirectional flow. The u z plot of figure 3 shows that small-scale convective flows are embedded within the large-scale cyclone. The convective structures can be elongated in one horizontal direction, especially in the regions of intense shear surrounding the LSV. There is no discernable box-size vertical flow associated with the presence of the LSV, so LSVs do not directly contribute to the outward heat transport, although they might have a feedback effect (Julien et al., 2012; Guervilly et al., 2014) Flow in the saturated phase for L y > 1 Figure 4 shows snapshots of the horizontal and vertical cross-sections of ω z and u z for L y = 8 in the final saturated state. The dominant feature in the plot of ω z is the unidirectional flow that consists of five jets, each with a cyclonic and anti-cyclonic component, directed along x. A multitude of small-scale convective vortices are embedded within the jets. The cyclonic LSVs that are prominent in the plots of ω z z during the initial growing phase shown in figure 2 are also visible in the saturated phase. However, as seen in figure 2c, they are no longer located in the middle of the bands of mean (i.e. x-averaged) cyclonic axial vorticity, but lie towards 6
8 Figure 4: Snapshots of the horizontal and vertical cross-sections of (a) the axial vorticity and (b) the vertical velocity for Ly = 8. 7
9 Figure 5: Snapshots of the depth-averaged axial vorticity (bottom panel) and profile of the mean u x as a function of y (top panel) for (a) L y = 2, (b) L y = 4 and (c) L y = 8. the flanks. The plot of u z shows that the convective structures can be elongated along the x-direction in the shear regions but that there is no discernible box-size vertical flow associated with the presence of the jets. Figure 5 shows snapshots of ω z z and of the profile of u x xz in the saturated stage for L y = 8 and for the intermediate cases L y = 2 and 4. The unidirectional flow consists of one jet for L y = 2 and two jets for L y = 4. No persistent unidirectional flow along y is produced. The presence of a unidirectional flow might be expected in this system by analogy with 2D turbulence, in which the upscale energy transfer feeds the largest mode available (e.g. Bouchet & Simonnet, 2009). However, as noted by (Frishman et al., 2017) in a study of forced 2D turbulence, the largest mode argument does not explain why there is more than one stable jet for L y = 4 and L y = 8. In all the cases studied, LSVs coexist with the jets. The LSVs are advected by the mean flow u x and always remain at the same location in y. At most two cyclones can coexist in the same cyclonic band and they remain in the flanks of their respective jets, with no cyclone sharing its y location. The cyclones are never observed at the location corresponding to u x xz = 0, which is the region of maximum shear y u x xz. Furthermore, no large-scale anticyclone is visible in the bands of anticyclonic vorticity; hence the asymmetry between large-scale cyclones and anticyclones is present for L y > L x, as it is when L y = L x. The mean flows themselves however are not asymmetric, since cyclonic and anticyclonic bands have comparable width (and hence comparable amplitude). 3.2 How much horizontal anisotropy is required to drive the jets? Since large-scale motions consist mostly of horizontal flows, we follow (Guervilly et al., 2014) in usingtheparameterγ = Re 2 /(3Re 2 z)todeterminetheirpresence; Γ > 1indicatesthepresenceof 8
10 α L y Figure 6: Flow anisotropy coefficient α, defined by (4), as a function of L y. either jets or LSVs. Such is the case for all the simulations reported in table 1. To distinguish between those simulations in which only LSVs form and those in which jets are present, we calculate the coefficient α defined by α = u2 x u 2 y u 2 x + u 2 y, (4) which measures the anisotropy of the flow; α is close to unity when jets are present. Figure 6 shows α as a function of L y ; crucially, only a small deviation of L y from L x (namely L y = 1.08 for the parameters used here) is required to observe the spontaneous emergence of jets. 3.3 Long term evolution and merging of the jets TostudythedriftingormergingofthejetsandLSVsovertime, weplotinfigure7thespace-time diagram of the axial vorticity averaged in x and z, ω z xz (the mean vorticity), together with the r.m.s. of the fluctuating component, defined as ω z z ω z xz. The mean vorticity contains information about both the jets and the LSV, while the fluctuating part contains information about the LSV. Over this time series, the jets do not drift or merge; indeed, this is representative of all the cases of table 1, irrespective of the number of jets. The absence of measurable drift is consistent with the fact that the horizontal average of the horizontal velocity is negligible in these simulations. LSVs can persist in the same band of mean cyclonic vorticity for the whole duration of the time integration or they might be intermittent. LSVs on opposite flanks of the same cyclonic band are advected in opposite directions and do not appear to interact with each other. However, in some occurrences, the formation of a new LSV might destabilise the neighbouring LSV, since the disappearance of an LSV sometimes closely follows the appearance of its immediate neighbour. It is plausible that the persistence or intermittency of the LSVs is linked to the width of the band of mean vorticity in which they are located. Indeed, in figure 7, the LSVs within the narrower band located around y = 1 are more intermittent than those within the wider band located around y = 3. The comparison of the behaviour of the LSVs for the cases L y = 2 and 2.5 (for which one jet is present) also supports this argument. Table 1 gives the number of jets, denoted by n, where, by jets, we mean unidirectional flows that persist in time averages. For clarification, we reiterate that we are defining a jet as consisting of both its cyclonic and anticyclonic component (i.e. both a red and a blue region in plots such as figure 2). Note that simulations that produce only LSVs (i.e. L y < 1.08) have n = 0; cases with L y 3 have multiple jets (n > 1). This result is somewhat surprising since, in a system where there is no gradient of the planetary rotation (i.e. no β-effect), one might expect 9
11 Figure7: Space-timediagramofthedepth-averagedaxialvorticityforEk = 10 5, Ra = , L x = 1 and L y = 4: (a) x-averaged, and (b) the r.m.s. of the fluctuating component. the upscale energy transfer to persist until the largest available scale is attained. Over the entire duration of the simulations with Ek = 10 5 there is no evidence of a trend for any slow drift or merger of the jets. That said, for such a small value of Ek, computational limitations restrict the time integration to less than one viscous timescale (which corresponds to Ek 1 in our time unit). In all the simulations reported so far, we have described the large-scale structure that gradually emerges from an initial small-scale perturbation and that persists over the entire time integration. The system spontaneously adopts a number of jets such that L y /(2n) L x when possible (i.e. when L y is large enough). This observation can be explained by the snapshots of the initial growing stage in figure 2: the LSVs form first and grow to the size L x, before extending into jets that retain the initial size of the LSV in the y-direction. It is however plausible that on longer timescales the jets could readjust their size by merging. In order to investigate very long term aspects of the dynamics, it is necessary to consider an increased value of the Ekman number. To this end, we perform a simulation with Ek = 10 4, Ra = (Ra/Ra c = 4.2) and (L x,l y ) = (1,4) for a duration of five viscous timescales. The time series of ˆRe x, ˆRe y and ˆRe z, together with the space-time diagram of ω z xz, are shown in figure 8. The initial phases of growth and saturation, before one viscous timescale has elapsed, are similar to the case already described for Ek = 10 5, with the formation of four jets. Note that no persistent LSV is observed in the jet flanks here because L x is not large enough to permit it at this value of Ek. Over approximately two viscous timescales, one of the mean cyclonic bands slowly drifts, preceding a relatively rapid expansion and merger, with only one jet remaining; this jet merger is accompanied by an increase by a factor two of ˆRe x. On long timescales, the jets therefore evolve by increasing their lengthscale as a result of the merger of bands of mean cyclonic vorticity, leading to an increase in the amplitude of u x. This relation between the amplitude of u x and the jet lengthscale is consistent with the results presented in 10
12 1000 a) ˆRe x ˆRey ˆRez time Figure 8: a) Time series of the r.m.s. values of the three velocity components and b) space-time diagramofthedepth-averagedaxialvorticityforek = 10 4, Ra/Ra c = 4.2and(L x,l y ) = (1,4). 11
13 Rex n = 1 n = n = 4 n = L y /n Figure 9: Re x as a function of L y /n (for cases with n > 0) for Ek = 10 5, Ra = and L x = ˆRe x ˆRey ˆRez time Figure 10: Time series of the r.m.s. values of the three velocity components for the simulation with (L x,l y ) = (1,1) whose initial condition is constructed from (L x,l y ) = (0.5,1) (Ek = 10 5, Ra = ). table 1 (for which Ek = 10 5 ). Figure 9 shows Re x as a function of the jet lengthscale L y /n for the cases with n > 0 of table 1, corroborating the linear dependence of Re x on L y /n. 3.4 Bistability of the large-scale flow In their study of forced 2D turbulence, Bouchet & Simonnet (2009) show that when the aspect ratio of the computational domain slightly exceeds unity, the flow switches randomly between a dipole vortex and a unidirectional flow. For the ratio L y /L x = 1.02, their system spends more time in the dipole vortex configuration, in contrast with the case of L y /L x = 1.04, for which the unidirectional flow configuration is preferred. These topological changes are found to be slow processes that occur over a viscous timescale, and which therefore can be observed only through long time integrations. In all our simulations reported in table 1 (all of which have Ek = 10 5 ), which are started from an initial small-scale perturbation, no random topological changes were observed, even for small deviations of L y /L x from unity. However, the results of (Bouchet & Simonnet, 2009) show that the solution might be dependent on initial conditions and may vary on viscous timescales. To study whether the two configurations LSV-only or jets can both occur for a given 12
14 600 a) ˆRe x ˆRey ˆRez time Figure 11: (a) Time series of the r.m.s. values of the three velocity components for the simulation with Ek = 10 4, Ra/Ra c = 4.2 and (L x,l y ) = (2,2.15). (b) (c) Snapshots of the z-averaged axial vorticity at the times indicated by the dashed vertical lines in (a). aspect ratio L y /L x, we should ideally perform simulations for several viscous timescales. Since computational constraints exclude this possibility for small Ekman numbers, we use two different approaches. First, we force the system to adopt a given large-scale configuration by using an initial condition of finite amplitude. To determine whether the jets can be stable when L x = L y, we use the data from a snapshot of the case (L x,l y ) = (0.5,1) (which is thus in the jet configuration), which we replicate in the x-direction to obtain a box of dimension (L x,l y ) = (1,1). This fabricated data, with the addition of a small amount of noise, is then used as an initial condition for the domain (L x,l y ) = (1,1). Figure 10 shows the time series of ˆRe x, ˆRe y and ˆRe z. The jet configuration is recognisable initially by ˆRe x > ˆRe y ˆRe z. After about 1000 time units, or equivalently about 25 turnover timescales based on the horizontal velocity and the box size, the jet configuration is lost in favour of the LSV configuration (where ˆRe x ˆRe y ). Furthermore, we perform simulations for which a snapshot from a case with (L x,l y ) = (1,1.08)(ajetconfiguration)isusedasaninitialconditionforeither(L x,l y ) = (1,1) or(l x,l y ) = (1.08,1.08). In bothcases, thelarge-scale flow rapidlyevolvesfrom jetstoanlsv, similarly to the case shown in figure 10. Based on these numerical simulations, we conclude that the jets become rapidly unstable when L x = L y. The second approach consists, as in 3.3, in increasing the Ekman number in order to simulate many viscous timescales. We perform a simulation with Ek = 10 4, Ra/Ra c = 4.2 and (L x,l y ) = (2,2.15) for a duration of 10 viscous timescales. The time series of ˆRe x, ˆRe y and ˆRe z are shown in figure 11a. An LSV forms initially, followed by spontaneous topological changes that occur after the elapse of more than one viscous timescale, which corresponds to 13
15 ˆRe x ˆRey ˆRez time Figure 12: Time series of the r.m.s. values of the three velocity components for the simulation with Ek = 10 4, Ra/Ra c = 4.2 and (L x,l y ) = (2,2.2) time units in figure 11a. The LSV configuration is characterised by ˆRe y > ˆRe z and is illustrated in figure 11b, which shows a snapshot of ω z z. The jet configuration is characterised by ˆRe y ˆRe z and is illustrated by figure 11c. For this particular aspect ratio (L y /L x = 1.075), the system spends most of its time in the LSV configuration, i.e. the probability of obtaining the LSV configuration is higher. As L y /L x decreases towards unity, the probability of a jet configuration decreases. Conversely, as L y /L x increases, the probability of achieving the LSV configuration decreases. This is illustrated in figure 12, which shows the times series of ˆRe x, ˆRe y and ˆRe z for (L x,l y ) = (2,2.2). We find that for L y /L x = 1.2, the system remains in the jet configuration during the whole time integration of 10 viscous timescales. For the parameters studied in this paragraph, there is a narrow window, given by L y /L x [1.06,1.15], in which both the jet state and the LSV-only state can exist. 3.5 Domain of existence When L y = L x, LSVs form when the convective vortices are constrained by rotation to be anisotropic (i.e. narrow in the horizontal directions and tall in the vertical) and when they are sufficiently energetic to merge (Guervilly et al., 2014). At least for the Ekman numbers considered in Guervilly et al. (2014) (Ek ), these conditions are quantified by (i) small local Rossby numbers, Ro l z = u z /(2Ωl) 0.1, where l is the typical horizontal lengthscale of the convection, and (ii) Ra/Ra c 3. To meet both criteria, the Ekman number must be suitably small. An additional condition for the formation of LSVs is that the computational domain must be wide enough to accommodate a sufficient number of convective cells in the horizontal direction (typically more than four convective cells along L x or L y ). To determine whether similar constraints apply to the formation of the jets, we first performed simulations at larger Ekman numbers than considered above. Table 2 lists a number of outputs from simulations run at Ek = and Ek = 10 4 for Ra/Ra c 4 and for various box configurations. The ratio Γ is again used to indicate the formation of dominant horizontal flows (related to the formation of large-scale flows) when greater than unity, and the coefficient α is used to indicate the formation of unidirectional flows when close to unity. For Ek = and all the box configurations, the local Rossby number is Ro l z = 0.12 and Γ 1; neither an LSV nor jets form because conditions (i) and (ii) cannot be met simultaneously for this Ekman number. Note that, for this value of Ek, only about three convective cells fit along a unit horizontal length; we thus have to use wide computational domains. When Ek is decreased to 10 4, the local Rossby number decreases to 0.09 for this value of Ra/Ra c. For all box configurations, Γ > 1, i.e. large-scale horizontal flows are predominant. For L x = L y = 2, 14
16 Ek Ra Ra/Ra c L x L y Re x Re y Re z Re Ro l z Γ α N x N y N z Table 2: Simulations performed in order to determine the domain of existence of the large-scale flows. a) L y = 1 L y = 2 b) L y = 1 L y = Γ 6 Γ Ra/Ra c Ro l z Figure 13: Γ as a function of (a) the normalised Rayleigh number and (b) the local Rossby number for Ek = α = 0, i.e. the large-scale flow is an LSV, as expected. For (L x,l y ) = (1,2), α > 0, reflecting the formation of jets. These observations are all confirmed by visual inspection of the axial vorticity (not shown here for brevity). These simulations thus show that the constraint Ek 10 4 is required for the formation of jets, a similar criterion to that for the formation of LSVs. To test whether both of the conditions (i) and (ii) apply to the formation of jets, we perform a number of simulations with increasing Ra for Ek = 10 5 and with (L x,l y ) = (1,2) (listed in table 2). Figure 13a shows the value of Γ as a function of Ra/Ra c for these cases, which is compared with cases with (L x,l y ) = (1,1) at the same Ek as used in Guervilly et al. (2014). Γ increases well above unity for both L y = 1 and L y = 2 when Ra/Ra c > 3. Condition (ii) is therefore quantitatively the same for both jets and LSVs. Figure 13b shows Γ as a function of Ro l z and indicates that condition (i) is quantitatively different for jets (L y = 2) and LSVs (L y = 1): the amplitude of the jets start to decay when the local Rossby number is greater than 0.08, which is smaller than the Rossby number at which Γ starts to decline for the LSV. We thus conclude that the jets are more fragile than the LSVs when the local Rossby number approaches 0.1. Note that for the largest Rayleigh number computed for L y = 2, Γ takes similar values for L y = 1 and L y = 2; however, the large-scale flow for L y = 2 still consists of an n = 1 jet accompanied by LSVs in the mean cyclonic region. 3.6 Effect on the heat transfer We anticipate that convective heat transfer might be impeded by the presence of jets due to either (i) the shear flow itself, or (ii) the local increase of the rotation rate in the regions of large-scale cyclonic vorticity. It should though be noted that the influence of a mean shear flow on convection is complex and may even enhance the efficiency of the heat transport (Clever & Busse, 1992). Effect (i) is expected to be important when most of the convective heat transport 15
17 31 30 Nu 29 n = 0 without LSV n = 0 with LSV 28 n = 1 n = 2 n = 4 n = Re x Figure 14: Nu as a function of the Reynolds number based on u x for the simulations of table 1 (L x = 1 and varying L y ) and for a reference case without a large-scale flow (L x = L y = 0.3). is directed parallel to the mean velocity gradients. Such is the case outside the tangent cylinder in spherical geometry (e.g. Yadav et al., 2016), and in planar geometry when the rotation axis is perpendicular to the direction of gravity (e.g. von Hardenberg et al., 2015). In our system, the effect of the shear is difficult to predict: on the one hand, the convective heat flux is perpendicular to the mean velocity gradients, so the shear might not be particularly disruptive; on the other, the elongation of the vertical flows in the x-direction observed in figure 4b indicates that the vertical velocity is affected by the shear. Effect (ii) is due to the well-known inhibiting effect of rotation on convection (Chandrasekhar, 1961) and has previously been observed in the presence of LSVs when L x = L y ; here the contribution of the LSV vorticity to the background rotation is not negligible locally (Guervilly et al., 2014). The decrease in the heat transport efficiency associated with LSVs depends on the horizontal box size since wider horizontal domains promote the formation of broader LSVs of increasing strength. In the case of jets, the efficiency of the heat transfer may also depend on their amplitude and lengthscale. The global efficiency of the heat transfer is often quantified by the Nusselt number, Nu, which is a measure of the total heat flux through the layer normalised by the heat flux in the absence of convective motions. In our system, in the absence of heat sources and with fixed-temperature boundary conditions, the Nusselt number is defined as Nu = 1+ θ 2, (5) where the angle brackets denote time- and volume-averages. To compare quantitatively the effect of the presence of large-scale flows on Nu, we need to calculate Nu in a reference case without large-scale flows. This can be obtained by considering a small horizontal box size that accommodates only about three convective cells, thereby preventing the formation of a largescale flow. For Ek = 10 5 and Ra = , the reference case is computed for L x = L y = 0.3. Figure 14 shows the Nusselt number as a function of Re x for the simulations of table 1 and for the reference case. Overall, Nu is smaller by approximately 10% in cases with a large-scale flow than it is in the reference case. As we might have expected, there is a decreasing trend of Nu with Re x for the cases with jets (n > 0). However, at similar values of Re x, Nu is relatively unaffected by the number of jets. The cases with only an LSV (n = 0) have smaller values of Nu than cases with jets with similar values of Re x. We now examine the local effect of the large-scale flows on the heat transfer. Figure 15 shows the time-averaged heat flux at the upper boundary, q = z T(z = 1) x, for (L x,l y ) = (1,4). For comparison, we also plot the time averages of u x xz and ω z xz. The heat flux is lower by approximately 25% in the regions of positive ω z xz than in the regions of negative ω z xz. 16
18 0.10 a) 0.05 ux xz ωz xz b) c) q y Figure 15: Profiles of (a) the time-averaged mean flow, u x xz, (b) the mean axial vorticity, ω z xz, and (c) the mean heat flux at the upper boundary, q; Ek = 10 5, Ra = , L x = 1 and L y = 4. The red bands correspond to the regions of positive mean vorticity. 17
19 This main difference is not due to the shearing of the convective flows by the jets because the highs of q are not correlated with regions of smallest shear. The heat flux is particularly low in the locations of the large-scale cyclone, which appear distinctly as peaks in the profile of ω z xz. This detrimental effect of large-scale cyclonic flow on the heat transport is due to the local increase of rotation rate, as in the case of the LSVs when L x = L y. This is consistent with the fact that the mean axial vorticity persistently reaches values of 0.2 (in units of 2Ω). In the reference case L x = L y = 0.3, q 18 on average, so this argument also explains why the anticyclonic regions enhance the heat flux. In cases with LSVs only (L y [1,1.05]), the mean axial vorticity persistently reaches values of 0.5, which explains why LSVs have a more pronounced effect than jets on Nu (figure 14). 4 Conclusion We have studied the formation of large-scale flows in rotating anisotropic convection; this has been achieved through numerical simulations in planar geometry, with horizontal anisotropy introduced by choosing a box aspect ratio with unequal horizontal dimensions (L y > L x ). For the standard case of L y /L x = 1, which has received much attention in the literature (e.g. Chan, 2007; Käpylä et al., 2011; Julien et al., 2012; Favier et al., 2014; Guervilly et al., 2014; Rubio et al., 2014; Stellmach et al., 2014), the stable configuration consists of a depth-invariant large-scale vortex, formed by the clustering of small-scale convective vortices, and which grows to the horizontal box size. When L y /L x is increased above unity, the system undergoes a transition from the LSV-only configuration to one in which depth-invariant unidirectional flows emerge. As for LSVs, the necessary conditions for the formation of these unidirectional flows are small Rossby numbers (based on the convective lengthscale and velocity) and sufficiently large Reynolds numbers. Over less than one global viscous timescale, persistent unidirectional flows with multiple jets are produced when L y 3L x. On longer timescales, the multiple jets merge to reach the largest available lengthscale in the system. Interestingly, we found that LSVs of size comparable with L x systematically coexist with the jets. The LSVs are located in the flanks of the jets and can be persistent or intermittent depending on the lengthscale of the jet. The transition from the LSV-only configuration to a configuration involving both jets and LSVs occurs for a small degree of anisotropy, namely L y /L x 1.08 when Ek = 10 5 and L x = 1. In the vicinity of this transition, the system is bistable; random changes between the two large-scale states are observed when the simulations are integrated over several viscous timescales. This bistable regime is observed only for a narrow interval of the box aspect ratio (L y /L x [1.06,1.15] for Ek = 10 4 and L x = 2). Large-scale flows are known to influence the efficiency of the heat transfer in rotating convection (e.g. von Hardenberg et al., 2015; Yadav et al., 2016). In our model, although there are no sizeable vertical flows associated with the presence of the LSVs or the jets, we find that the large-scale flows modify the transport properties of the convection by changing locally the rotation rate. As the axial vorticity of the large-scale flows is skewed towards cyclonic vorticity in our simulations, these flows reduce the efficiency of the convective heat transport. In more realistic planetary conditions at smaller Ekman and Rossby numbers, the skewness of the axial vorticity is expected to diminish, and hence the effect of the large-scale flows on the efficiency of the heat transfer may be different (Julien et al., 2012). In all of the simulations described above, the domains are characterised by L x = 1, with differing values of L y. Although a systematic exploration of L x -L y space is not computationally feasible when both L x and L y exceed unity, it is nonetheless of interest to investigate the nature of the solution in the regime 1 < L x < L y, given that, for 1 < L x = L y, the attracting solution will consist of only an LSV. To this end, we simulated the case of L x = 2, L y = 4, with the standard parameter values of Ek = 10 5, Ra = and Pr = 1. Following a small perturbation to the static state, the general form of the evolution is as described above for the 18
20 Figure 16: Snapshot of the depth-averaged axial vorticity in the saturated regime for the case of L x = 2, L y = 4. cases with L x = 1, with an initial rapid growth and decline in the kinetic energy, followed by a slow evolution to the large-scale state (as shown, for example, in figure 1). In this extended domain, the slow evolution is accordingly lengthened, with equilibration of the kinetic energy not occurring until t Figure 16 shows a snapshot of the depth-averaged axial vorticity in the long-time stationary regime, with one jet and with LSVs embedded in the flanks of the cyclonic component; the LSVs extend essentially across the entire x-domain. In comparison with the case of L x = 1, L y = 2, which has the same aspect ratio in the horizontal plane, the flows in the extended domain (L x = 2, L y = 4) are, as expected, much more vigorous, with, for example, Re x = 9870, in comparison with Re x = 5548 in the smaller domain; as a result, this leads to a reduction in the heat transfer (Nu = 26.9 in comparison with Nu = 28.6). In spherical shell models of rotating convection in the molecular envelope of gas giant planets, shallow large-scale vortices have been observed to coexist with deep zonal flows (i.e. flows that extend throughout most of the layer depth) (Heimpel et al., 2016). The persistence of deep LSVs when the horizontal symmetry is broken, together with their coexistence with jets as demonstrated here, is a novel feature in rotating convection. We have used a simplified model that does not permit the saturation of the lengthscale of the large-scale flows as a result of either the curvature of the boundaries or the latitudinal variation of the rotation background. Our model is therefore unable to provide predictions for the amplitude and lengthscale of the largescale flows in terms of the properties of the convective flows. However, it does offer a valuable proof of concept of the possible coexistence of deep LSVs and jets in rotating convection. In numerical simulations using periodic boundaries, the dynamics should ideally not be influenced by the choice of the computational domain size. However, large-scale flows can be affected by confinement effects due to the finite domain size, as we find here. Numerical simulations in very wide computational domains are required to study the influence of confinement on the large-scale dynamics, but, unfortunately, such simulations are currently out of reach of 3D models of rapidly rotating convection. The analogy between the large-scale dynamics observed in 3D simulations and in forced 2D turbulence simulations (e.g. Frishman et al. (2017)) supports the study of the evolution of large-scale flows over long timescales using 2D turbulence models in the wide domains that 3D models cannot explore. It would be particularly interesting to perform 2D simulations with a forcing possessing similar statistical properties to the small-scale convective buoyancy driving found in 3D models. An alternative and complementary approach to 3D direct numerical simulations, proposed by Julien & Knobloch (2007), is to model rotating Rayleigh-Bénard convection using an asymptotically reduced set of PDEs valid in the limit of small Rossby number. Certain results obtained with the reduced system are found to be in 19
Generation of magnetic fields by large-scale vortices in rotating convection
Generation of magnetic fields by large-scale vortices in rotating convection Céline Guervilly, David Hughes & Chris Jones School of Mathematics, University of Leeds, UK Generation of the geomagnetic field
More informationarxiv: v1 [physics.flu-dyn] 17 Aug 2017
Multiple zonal jets and convective heat transport barriers in a quasi-geostrophic model of planetary cores Céline Guervilly 1 & Philippe Cardin 2 1 School of Mathematics, Statistics and Physics, Newcastle
More informationarxiv: v1 [physics.flu-dyn] 27 Aug 2014
arxiv:148.6483v1 [physics.flu-dyn] 27 Aug 214 Inverse cascade and symmetry breaking in rapidly-rotating Boussinesq convection B. Favier, 1, a) L.J. Silvers, 1 and M.R.E. Proctor 2 1) Centre for Mathematical
More informationConvection-driven dynamos in the limit of rapid rotation
Convection-driven dynamos in the limit of rapid rotation Michael A. Calkins Jonathan M. Aurnou (UCLA), Keith Julien (CU), Louie Long (CU), Philippe Marti (CU), Steven M. Tobias (Leeds) *Department of Physics,
More informationSimulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
More informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More informationInverse cascade and symmetry breaking in rapidly rotating Boussinesq convection
Inverse cascade and symmetry breaking in rapidly rotating Boussinesq convection Benjamin Favier, Lara Silvers, Michael Proctor To cite this version: Benjamin Favier, Lara Silvers, Michael Proctor. Inverse
More informationGeostrophic turbulence and the formation of large scale structure
Geostrophic turbulence and the formation of large scale structure Edgar Knobloch University of California, Berkeley, CA 9472, USA knobloch@berkeley.edu http://tardis.berkeley.edu Ian Grooms, Keith Julien,
More informationRapidly Rotating Rayleigh-Bénard Convection
Rapidly Rotating Rayleigh-Bénard Convection Edgar Knobloch Department of Physics University of California at Berkeley 27 February 2008 Collaborators: Keith Julien, Applied Mathematics, Univ. of Colorado,
More informationES265 Order of Magnitude Phys & Chem Convection
ES265 Order of Magnitude Phys & Chem Convection Convection deals with moving fluids in which there are spatial variations in temperature or chemical concentration. In forced convection, these variations
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
This is a repository copy of Combining music and life story to enhance participation in family interaction in semantic dementia: a longitudinal study of one family s experience. White Rose Research Online
More informationTheoretical Geomagnetism. Lecture 3. Core Dynamics I: Rotating Convection in Spherical Geometry
Theoretical Geomagnetism Lecture 3 Core Dynamics I: Rotating Convection in Spherical Geometry 1 3.0 Ingredients of core dynamics Rotation places constraints on motions. Thermal (and chemical buoyancy)
More informationA Further Study of Dynamics of the Four-Day Circulation. By Yoshihisa Matsuda
February 1982 Y. Matsuda 245 A Further Study of Dynamics of the Four-Day Circulation in the Venus Atmosphere By Yoshihisa Matsuda Department of Astronomy and Earth Sciences, Tokyo Gakugei University, Koganei
More informationThis is a repository copy of Aggregation of growing crystals in suspension: III. Accounting for adhesion and repulsion.
This is a repository copy of Aggregation of growing crystals in suspension: III. Accounting for adhesion and repulsion. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/88281/
More informationGlobal magnetorotational instability with inflow The non-linear regime
Global magnetorotational instability with inflow The non-linear regime Evy Kersalé PPARC Postdoctoral Research Associate Dept. of Appl. Math. University of Leeds Collaboration: D. Hughes & S. Tobias (Dept.
More informationScaling laws for planetary dynamos driven by helical waves
Scaling laws for planetary dynamos driven by helical waves P. A. Davidson A. Ranjan Cambridge What keeps planetary magnetic fields alive? (Earth, Mercury, Gas giants) Two ingredients of the early theories:
More informationThis is a repository copy of Measure for Measure. Semi-cylindrical illuminance: a semi-conceived measure?.
This is a repository copy of Measure for Measure. Semi-cylindrical illuminance: a semi-conceived measure?. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/112512/ Version:
More informationThis is a repository copy of Conceptualizing a circular framework of supply chain resource sustainability.
This is a repository copy of Conceptualizing a circular framework of supply chain resource sustainability. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk// Version: Accepted
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More informationThis is a repository copy of An in vivo platform for identifying inhibitors of protein aggregation.
This is a repository copy of An in vivo platform for identifying inhibitors of protein aggregation. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/91486/ Version: Supplemental
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationVortex Dynamos. Steve Tobias (University of Leeds) Stefan Llewellyn Smith (UCSD)
Vortex Dynamos Steve Tobias (University of Leeds) Stefan Llewellyn Smith (UCSD) An introduction to vortices Vortices are ubiquitous in geophysical and astrophysical fluid mechanics (stratification & rotation).
More informationThis is a repository copy of The Uptake of HO ₂ on Meteoric Smoke Analogues.
This is a repository copy of The Uptake of HO ₂ on Meteoric Smoke Analogues. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/107542/ Version: Accepted Version Article: James,
More informationThis is a repository copy of Microbiology: Mind the gaps in cellular evolution.
This is a repository copy of Microbiology: Mind the gaps in cellular evolution. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/114978/ Version: Accepted Version Article:
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationMax Planck Institut für Plasmaphysik
ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately
More informationPrototype Instabilities
Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies
More informationAnisotropic turbulence in rotating magnetoconvection
Anisotropic turbulence in rotating magnetoconvection André Giesecke Astrophysikalisches Institut Potsdam An der Sternwarte 16 14482 Potsdam MHD-Group seminar, 2006 André Giesecke (AIP) Anisotropic turbulence
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationAsymptotic theory for torsional convection in rotating fluid spheres
Under consideration for publication in J. Fluid Mech. 1 Asymptotic theory for torsional convection in rotating fluid spheres By KEKE Z H A N G 1, KAMENG L A M A N D DALI K O N G 3 1,3 College of Engineering,
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationPassive Scalars in Stratified Turbulence
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Passive Scalars in Stratified Turbulence G. Brethouwer Linné Flow Centre, KTH Mechanics, SE-100 44 Stockholm, Sweden E. Lindborg Linné Flow Centre,
More informationCoriolis effects on the elliptical instability in cylindrical and spherical rotating containers
18 ème Congrès Français de Mécanique Grenoble, 7-31 août 007 Coriolis effects on the elliptical instability in cylindrical and spherical rotating containers Michael Le Bars, Stéphane Le Dizès & Patrice
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More informationATMOSPHERIC AND OCEANIC FLUID DYNAMICS
ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-scale Circulation G E O F F R E Y K. V A L L I S Princeton University, New Jersey CAMBRIDGE UNIVERSITY PRESS An asterisk indicates more advanced
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationTurbulent three-dimensional MHD dynamo model in spherical shells: Regular oscillations of the dipolar field
Center for Turbulence Research Proceedings of the Summer Program 2010 475 Turbulent three-dimensional MHD dynamo model in spherical shells: Regular oscillations of the dipolar field By R. D. Simitev, F.
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationChapter 9 NATURAL CONVECTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition in SI Units Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 9 NATURAL CONVECTION PM Dr Mazlan Abdul Wahid Universiti Teknologi
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
This is a repository copy of A comparison of the sampling efficiency of bioaerosol samplers and particle counters in natural and controlled environments. White Rose Research Online URL for this paper:
More informationCoherent structures in stably stratified plane Couette flow
Coherent structures in stably stratified plane Couette flow D. Olvera * & R. R. Kerswell School of Mathematics, University of Bristol, Bristol, UK. * do2542@bristol.ac.uk Abstract A large body of recent
More informationThe interaction of vorticity and rate-of-strain in homogeneous sheared turbulence
PHYSICS OF FLUIDS VOLUME 12, NUMBER 4 APRIL 2000 The interaction of vorticity and rate-of-strain in homogeneous sheared turbulence K. K. Nomura a) and P. J. Diamessis Department of Mechanical and Aerospace
More informationConvective Heat and Mass Transfer Prof. A. W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A. W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No.# 01 Lecture No. # 41 Natural Convection BLs So far we have considered
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationThis is a repository copy of An application of Kozeny Carman flow resistivity model to predict the acoustical properties of polyester fibre.
This is a repository copy of An application of Kozeny Carman flow resistivity model to predict the acoustical properties of polyester fibre. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/93426/
More informationScaling properties of convection-driven dynamos in rotating spherical shells and application to planetary magnetic fields
Geophys. J. Int. (006) 66, 97 4 doi: 0./j.365-46X.006.03009.x Scaling properties of convection-driven dynamos in rotating spherical shells and application to planetary magnetic fields U. R. Christensen
More informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationINVESTIGATION OF TRANSITIONAL AND TURBULENT HEAT AND MOMENTUM TRANSPORT IN A ROTATING CAVITY
INVESTIGATION OF TRANSITIONAL AND TURBULENT HEAT AND MOMENTUM TRANSPORT IN A ROTATING CAVITY Ewa Tuliska-Snitko, Wojciech Majchrowski, Kamil Kiełcewski Institute of Thermal Engineering, Ponan University
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationVortices in accretion discs: formation process and dynamical evolution
Vortices in accretion discs: formation process and dynamical evolution Geoffroy Lesur DAMTP (Cambridge UK) LAOG (Grenoble) John Papaloizou Sijme-Jan Paardekooper Giant vortex in Naruto straight (Japan)
More informationDynamics of plumes driven by localized heating in a stably stratified ambient
Dynamics of plumes driven by localized heating in a stably stratified ambient Abstract Juan M. Lopez 1 and Francisco Marques 2 1 School of Mathematical and Statistical Sciences, Arizona State University,
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationBuoyancy Fluxes in a Stratified Fluid
27 Buoyancy Fluxes in a Stratified Fluid G. N. Ivey, J. Imberger and J. R. Koseff Abstract Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed
More informationOn dynamo action produced by boundary thermal coupling
On dynamo action produced by boundary thermal coupling Binod Sreenivasan To cite this version: Binod Sreenivasan. On dynamo action produced by boundary thermal coupling. Physics of the Earth and Planetary
More informationChapter 5. Experimental Results - Dynamics
Chapter 5 Experimental Results - Dynamics 16 5.1 Overview The dynamics of the mean flow fields described in Chapter 4 will be presented in this chapter using mainly two tools: the power spectral density
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationSpherical Shallow Water Turbulence: Cyclone-Anticyclone Asymmetry, Potential Vorticity Homogenisation and Jet Formation
Spherical Shallow Water Turbulence: Cyclone-Anticyclone Asymmetry, Potential Vorticity Homogenisation and Jet Formation Jemma Shipton Department of Atmospheric, Oceanic and Planetary Physics, University
More informationarxiv: v1 [astro-ph.sr] 3 Oct 2011
Small-scale dynamo action in rotating compressible convection B. Favier and P.J. Bushby School of Mathematics and Statistics, Newcastle University, Newcastle upon Tyne NE 7RU, UK arxiv:.374v [astro-ph.sr]
More informationNatalia Tronko S.V.Nazarenko S. Galtier
IPP Garching, ESF Exploratory Workshop Natalia Tronko University of York, York Plasma Institute In collaboration with S.V.Nazarenko University of Warwick S. Galtier University of Paris XI Outline Motivations:
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationSize Scaling and Nondiffusive Features of Electron Heat Transport in Multi-Scale Turbulence
Size Scaling and Nondiffusive Features of Electron Heat Transport in Multi-Scale Turbulence Z. Lin 1, Y. Xiao 1, W. J. Deng 1, I. Holod 1, C. Kamath, S. Klasky 3, Z. X. Wang 1, and H. S. Zhang 4,1 1 University
More informationOPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW
OPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW Shingo Motoki, Genta Kawahara and Masaki Shimizu Graduate School of Engineering Science Osaka University 1-3 Machikaneyama, Toyonaka, Osaka
More informationFlow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers ( )
Advances in Fluid Mechanics VII 391 Flow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers (10 6 10 9 ) R. L. Frederick & S. Courtin
More informationRotating dynamo turbulence: theoretical and numerical insights
Rotating dynamo turbulence: theoretical and numerical insights Nathanaël Schaeffer 1, H.-C. Nataf 1, A. Fournier 2, J. Aubert 2, D. Jault 1, F. Plunian 1, P. Zitzke 1 1: ISTerre / CNRS /Université Joseph
More informationThis is a repository copy of Internal twist drill coolant channel modelling using computational fluid dynamics.
This is a repository copy of Internal twist drill coolant channel modelling using computational fluid dynamics. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/89798/ Version:
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationThis is a repository copy of Myths about HIV and AIDS among serodiscordant couples in Malawi.
This is a repository copy of Myths about HIV and AIDS among serodiscordant couples in Malawi. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk// Version: Accepted Version Article:
More informationDynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit. Transport Physics of the Density Limit
Dynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit Transport Physics of the Density Limit R. Hajjar, P. H. Diamond, M. Malkov This research was supported by the U.S. Department of Energy,
More informationPlanetary core dynamics and convective heat transfer scalingy
Geophysical and Astrophysical Fluid Dynamics, Vol. 101, Nos. 5 6, October December 2007, 327 345 Planetary core dynamics and convective heat transfer scalingy J. M. AURNOU* Department of Earth and Space
More informationApplied Computational Fluid Dynamics
Lecture 9 - Kolmogorov s Theory Applied Computational Fluid Dynamics Instructor: André Bakker André Bakker (2002-2005) Fluent Inc. (2002) 1 Eddy size Kolmogorov s theory describes how energy is transferred
More informationThis is a repository copy of Observation of electrolytic capacitor ageing behaviour for the purpose of prognostics.
This is a repository copy of Observation of electrolytic capacitor ageing behaviour for the purpose of prognostics. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/104114/
More informationThe Magnetorotational Instability
The Magnetorotational Instability Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics March 10, 2014 These slides are based off of Balbus & Hawley (1991), Hawley
More informationNumerical Study of the Moving Boundary Problem During Melting Process in a Rectangular Cavity Heated from Below
American Journal of Applied Sciences 4 (4): 25-256, 2007 ISSN 546-9239 2007 Science Publications Corresponding Author: Numerical Study of the Moving Boundary Problem During Melting Process in a Rectangular
More informationControl of Neo-classical tearing mode (NTM) in advanced scenarios
FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline
More informationSchool and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks
2292-13 School and Conference on Analytical and Computational Astrophysics 14-25 November, 2011 Angular momentum transport in accretion disks Gianluigi Bodo Osservatorio Astronomico, Torino Italy Angular
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationVARYING THE SPHERICAL SHELL GEOMETRY IN ROTATING THERMAL CONVECTION
Geophysical and Astrophysical Fluid Dynamics Vol. 98, No. 2, April 2004, pp. 153 169 VARYING THE SPHERICAL SHELL GEOMETRY IN ROTATING THERMAL CONVECTION F.M. AL-SHAMALI a, *, M.H. HEIMPEL a and J.M. AURNOU
More informationThis is a repository copy of Optimised PEI impregnation of activated carbons - Enhancement of CO2 capture under post-combustion conditions.
This is a repository copy of Optimised PEI impregnation of activated carbons - Enhancement of CO2 capture under post-combustion conditions. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/89498/
More informationPlanetary dynamos: Dipole-multipole transition and dipole reversals
Planetary dynamos: Dipole-multipole transition and dipole reversals Ulrich Christensen Max-Planck-Institute for Solar System Research Katlenburg-Lindau, Germany in collaboration with Hagay Amit, Julien
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationScaling laws for dynamos in rotating spheres: from planets to stars
Scaling laws for dynamos in rotating spheres: from planets to stars Ulrich Christensen Max Planck Institute for Solar System Research, Katlenburg- Lindau, Germany In collaboration with: Julien Aubert,
More informationNonlinear MRI. Jeremy Goodman. CMPD/CMSO Winter School 7-12 January 2008, UCLA
Nonlinear MRI Jeremy Goodman CMPD/CMSO Winter School 7-12 January 2008, UCLA Questions for the nonlinear regime How does the linear instability saturate? i.e., What nonlinear mode-mode interactions brake
More informationHomogeneous Rayleigh-Bénard convection
Slide 1 Homogeneous Rayleigh-Bénard convection scaling, heat transport and structures E. Calzavarini and F. Toschi, D. Lohse, R. Tripiccione, C. R. Doering, J. D. Gibbon, A. Tanabe Euromech Colloquium
More informationSimilarity and incomplete similarity
TIFR, Mumbai, India Refresher Course in Statistical Mechanics HBCSE Mumbai, 12 November, 2013 Copyright statement Copyright for this work remains with. However, teachers are free to use them in this form
More informationAsymptotic analysis of long-time behaviour of zonal flows in two-dimensional turbulence on a β plane. Kiori OBUSE. January 2012
RIMS-1737 Asymptotic analysis of long-time behaviour of zonal flows in two-dimensional turbulence on a β plane By Kiori OBUSE January 212 RESEARCH INSTITUTE FOR MATHEMATICAL SCIENCES KYOTO UNIVERSITY,
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationMagnetic Field Intensification and Small-scale Dynamo Action in Compressible Convection
Magnetic Field Intensification and Small-scale Dynamo Action in Compressible Convection Paul Bushby (Newcastle University) Collaborators: Steve Houghton (Leeds), Nigel Weiss, Mike Proctor (Cambridge) Magnetic
More informationBoussinesq dynamics of an idealized tropopause
Abstract Boussinesq dynamics of an idealized tropopause Olivier Asselin, Peter Bartello and David Straub McGill University olivier.asselin@mail.mcgill.ca The near-tropopause flow has been interpreted using
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationModeling the atmosphere of Jupiter
Modeling the atmosphere of Jupiter Bruce Turkington UMass Amherst Collaborators: Richard S. Ellis (UMass Professor) Andrew Majda (NYU Professor) Mark DiBattista (NYU Postdoc) Kyle Haven (UMass PhD Student)
More informationGeodynamo Simulations and Electromagnetic Induction Studies
Chapter 2 Solid Earth Simulation Geodynamo Simulations and Electromagnetic Induction Studies Project Representative Yozo Hamano Japan Agency for Marine-Earth Science and Technology Authors Yozo Hamano
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationTurbulent Dynamos: The Quick and the Fast
In preparation. 1 Turbulent Dynamos: The Quick and the Fast By S T E V E N M. T O B I A S 1 A N D F A U S T O C A T T A N E O 2 1 Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT,
More information