Planetary rotation organizes fluid motions into coherent, long-lived swirls, known as large-scale vortices (LSVs), which play an important role in the dynamics and long-term evolution of geophysical and astrophysical fluids. Here, using direct numerical simulations, we show that LSVs in rapidly rotating mixed convective and stably stratified fluids, which approximates the two-layer, turbulent-stratified dynamics of many geophysical and astrophysical fluids, have a generic shape and that their size can be predicted. We show that LSVs emerge in the convection zone from upscale energy transfers and can penetrate into the stratified layer. At the convective-stratified interface, the LSV cores have a positive buoyancy anomaly. Due to the thermal wind constraint, this buoyancy anomaly leads to winds in the stratified layer that decay over a characteristic vertical length scale. Thus LSVs take the shape of a depth-invariant cylinder with a finite-size radius in the turbulent layer and of a penetrating half dome in the stratified layer. Importantly, we demonstrate that when LSVs penetrate all the way through the stratified layer and reach a boundary that is no-slip, they saturate by boundary friction. We provide a prediction for the penetration depth and maximum radius of LSVs as a function of the LSV vorticity, the stratified layer depth, and the stratification. Our results, which apply for cyclonic LSVs, suggest that LSVs in slowly rotating stars and Earth's liquid core are confined to the convective layer, while in Earth's atmosphere and oceans they can penetrate far into the stratified layer.

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

Complex mixing and magnetic field generation occurs within stellar interiors particularly where there is a strong shear flow. To obtain a comprehensive understanding of these processes, it is necessary to study the complex dynamics of shear regions. Due to current observational limitations, it is necessary to investigate the inevitable small-scale dynamics via numerical calculations. Here, we examine direct numerical calculations of a local model of unstable shear flows in a compressible polytropic fluid primarily in a two-dimensional domain, where we focus on determining how key parameters affect the global properties and characteristics of the resulting saturated turbulent phase. We consider the effect of varying both the viscosity and the thermal diffusivity on the non-linear evolution. Moreover, our main focus is to understand the global properties of the saturated phase, in particular estimating for the first time the spread of the shear region from an initially hyperbolic tangent velocity profile. We find that the vertical extent of the mixing region in the saturated regime is generally determined by the initial Richardson number of the system. Further, the characteristic quantities of the turbulence, i.e. typical length-scale and the root-mean-square velocity are found to depend on both the Richardson number, and the thermal diffusivity. Finally, we present our findings of our investigation into saturated flows of a `secular' shear instability in the low P\'eclet number regime with large Richardson numbers.

Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.

Shear flows are ubiquitous in astrophysical objects including planetary and stellar interiors, where their dynamics can have significant impact on thermo-chemical processes. Investigating the complex dynamics of shear flows requires numerical calculations that provide a long time evolution of the system. To achieve a sufficiently long lifetime in a local numerical model the system has to be forced externally. However, at present, there exist several different forcing methods to sustain large-scale shear flows in local models. In this paper we examine and compare various methods used in the literature in order to resolve their respective applicability and limitations. These techniques are compared during the exponential growth phase of a shear flow instability, such as the Kelvin-Helmholtz (KH) instability, and some are examined during the subsequent non-linear evolution. A linear stability analysis provides reference for the growth rate of the most unstable modes in the system and a detailed analysis of the energetics provides a comprehensive understanding of the energy exchange during the system's evolution. Finally, we discuss the pros and cons of each forcing method and their relation with natural mechanisms generating shear flows.

Shear flows have an important impact on the dynamics in an assortment of different astrophysical objects including accreditation discs and stellar interiors. Investigating shear flow instabilities in a polytropic atmosphere provides a fundamental understanding of the motion in stellar interiors where turbulent motions, mixing processes, as well as magnetic field generation takes place. Here, a linear stability analysis for a fully compressible fluid in a two-dimensional Cartesian geometry is carried out. Our study focuses on determining the critical Richardson number for different Mach numbers and the destabilising effects of high thermal diffusion. We find that there is a deviation of the predicted stability threshold for moderate Mach number flows along with a significant effect on the growth rate of the linear instability for small P\'eclet numbers. We show that in addition to a Kelvin-Helmholtz instability a Holmboe instability can appear and we discuss the implication of this in stellar interiors.

Regions of quiet Sun generally exhibit a complex distribution of small-scale magnetic field structures, which interact with the near-surface turbulent convective motions. Furthermore, it is probable that some of these magnetic fields are generated locally by a convective dynamo mechanism. In addition to the well-known granular and supergranular convective scales, various observations have indicated that there is an intermediate scale of convection, known as mesogranulation, with vertical magnetic flux concentrations accumulating preferentially at mesogranular boundaries. Our aim is to investigate the small-scale dynamo properties of a convective flow that exhibits both granulation and mesogranulation, comparing our findings with solar observations. Adopting an idealised model for a localised region of quiet Sun, we use numerical simulations of compressible magnetohydrodynamics, in a 3D Cartesian domain, to investigate the parametric dependence of this system (focusing particularly upon the effects of varying the aspect ratio and the Reynolds number). In purely hydrodynamic convection, we find that mesogranulation is a robust feature of this system provided that the domain is wide enough to accommodate these large-scale motions. The mesogranular peak in the kinetic energy spectrum is more pronounced in the higher Reynolds number simulations. We investigate the dynamo properties of this system in both the kinematic and the nonlinear regimes and we find that the dynamo is always more efficient in larger domains, when mesogranulation is present. Furthermore, we use a filtering technique in Fourier space to demonstrate that it is indeed the larger scales of motion that are primarily responsible for driving the dynamo. In the nonlinear regime, the magnetic field distribution compares very favourably to observations, both in terms of the spatial distribution and the measured field strengths.

We perform one of the first studies into the nonlinear evolution of tidally excited inertial waves in a uniformly rotating fluid body, exploring a simplified model of the fluid envelope of a planet (or the convective envelope of a solar-type star) subject to the gravitational tidal perturbations of an orbiting companion. Our model contains a perfectly rigid spherical core, which is surrounded by an envelope of incompressible uniform density fluid. The corresponding linear problem was studied in previous papers which this work extends into the nonlinear regime, at moderate Ekman numbers (the ratio of viscous to Coriolis accelerations). By performing high-resolution numerical simulations, using a combination of pseudo-spectral and spectral element methods, we investigate the effects of nonlinearities, which lead to time-dependence of the flow and the corresponding dissipation rate. Angular momentum is deposited non-uniformly, leading to the generation of significant differential rotation in the initially uniformly rotating fluid, i.e. the body does not evolve towards synchronism as a simple solid body rotator. This differential rotation modifies the properties of tidally excited inertial waves, changes the dissipative properties of the flow, and eventually becomes unstable to a secondary shear instability provided that the Ekman number is sufficiently small. Our main result is that the inclusion of nonlinearities eventually modifies the flow and the resulting dissipation from what linear calculations would predict, which has important implications for tidal dissipation in fluid bodies. We finally discuss some limitations of our simplified model, and propose avenues for future research to better understand the tidal evolution of rotating planets and stars.

Mean-field dynamo theory suggests that turbulent convection in a rotating layer of electrically-conducting fluid produces a significant alpha-effect, which is one of the key ingredients in any mean-field dynamo model. Provided that this alpha-effect operates more efficiently than (turbulent) magnetic diffusion, such a system should be capable of sustaining a large-scale dynamo. However, in the Boussinesq model that was considered by Cattaneo&Hughes (2006) the dynamo produced small-scale, intermittent magnetic fields with no significant large-scale component. In this paper, we consider the compressible analogue of the rotating convective layer that was considered by Cattaneo&Hughes (2006). Varying the horizontal scale of the computational domain, we investigate the dependence of the dynamo upon the rotation rate. Our simulations indicate that these turbulent compressible flows can drive a small-scale dynamo but, even when the layer is rotating very rapidly (with a mid-layer Taylor number of Ta=10^8), we find no evidence for the generation of a significant large-scale component of the magnetic field on a dynamical timescale. Like Cattaneo&Hughes(2006), we measure a negligible (time-averaged) alpha-effect when a uniform horizontal magnetic field is imposed across the computational domain. Although the total horizontal magnetic flux is a conserved quantity in these simulations, the (depth-dependent) horizontally-averaged magnetic field always exhibits strong fluctuations. If these fluctuations are artificially suppressed within the code, we measure a significant mean electromotive force that is comparable to that found in related calculations in which the alpha-effect is measured using the test-field method, even though we observe no large-scale dynamo action.

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

We study dynamo action in a convective layer of electrically-conducting, compressible fluid, rotating about the vertical axis. At the upper and lower bounding surfaces, perfectly-conducting boundary conditions are adopted for the magnetic field. Two different levels of thermal stratification are considered. If the magnetic diffusivity is sufficiently small, the convection acts as a small-scale dynamo. Using a definition for the magnetic Reynolds number $R_M$ that is based upon the horizontal integral scale and the horizontally-averaged velocity at the mid-layer of the domain, we find that rotation tends to reduce the critical value of $R_M$ above which dynamo action is observed. Increasing the level of thermal stratification within the layer does not significantly alter the critical value of $R_M$ in the rotating calculations, but it does lead to a reduction in this critical value in the non-rotating cases. At the highest computationally-accessible values of the magnetic Reynolds number, the saturation levels of the dynamo are similar in all cases, with the mean magnetic energy density somewhere between 4 and 9% of the mean kinetic energy density. To gain further insights into the differences between rotating and non-rotating convection, we quantify the stretching properties of each flow by measuring Lyapunov exponents. Away from the boundaries, the rate of stretching due to the flow is much less dependent upon depth in the rotating cases than it is in the corresponding non-rotating calculations. It is also shown that the effects of rotation significantly reduce the magnetic energy dissipation in the lower part of the layer. We also investigate certain aspects of the saturation mechanism of the dynamo.

This article is focused on the dynamics of a rotating electrically conducting fluid in a turbulent state. As inside the Earth's core or in various industrial processes, a flow is altered by the presence of both background rotation and a large scale magnetic field. In this context, we present a set of 3D direct numerical simulations of incompressible decaying turbulence. We focus on parameters similar to the ones encountered in geophysical and astrophysical flows, so that the Rossby number is small, the interaction parameter is large, but the Elsasser number, defining the ratio between Coriolis and Lorentz forces, is about unity. These simulations allow to quantify the effect of rotation and thus inertial waves on the growth of magnetic fluctuations due to Alfv\'en waves. Rotation prevents the occurrence of equipartition between kinetic and magnetic energies, with a reduction of magnetic energy at decreasing Elsasser number {\Lambda}. It also causes a decrease of energy transfer mediated by cubic correlations. In terms of flow structure, a decrease of {\Lambda} corresponds to an increase in the misalignment of velocity and magnetic field.