We continue our investigation into the nonlinear evolution of the Goldreich-Schubert-Fricke (GSF) instability in differentially rotating radiation zones. This instability may be a key player in transporting angular momentum in stars and giant planets, but its nonlinear evolution remains mostly unexplored. In a previous paper we considered the equatorial instability, whereas here we simulate the instability at a general latitude for the first time. We adopt a local Cartesian Boussinesq model in a modified shearing box for most of our simulations, but we also perform some simulations with stress-free, impenetrable, radial boundaries. We first revisit the linear instability and derive some new results, before studying its nonlinear evolution. The instability is found to behave very differently compared with its behaviour at the equator. In particular, here we observe the development of strong zonal jets ("layering" in the angular momentum), which can considerably enhance angular momentum transport, particularly in axisymmetric simulations. The jets are, in general, tilted with respect to the local gravity by an angle that corresponds initially with that of the linear modes, but which evolves with time and depends on the strength of the flow. The instability transports angular momentum much more efficiently (by several orders of magnitude) than it does at the equator, and we estimate that the GSF instability could contribute to the missing angular momentum transport required in both red giant and subgiant stars. It could also play a role in the long-term evolution of the solar tachocline and the atmospheric dynamics of hot Jupiters.

Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thin stably-stratified interfaces, as a result of double-diffusive convection. These layers could have important effects on wave propagation and tidal dissipation that have not been fully explored. We analyse the effects of these layers on the propagation and transmission of internal waves within giant planets, extending prior work in a local Cartesian model. We adopt a simplified global Boussinesq planetary model in which we explore the internal waves in a non-rotating spherical body. We begin by studying the free modes of a region containing semi-convective layers. We then analyse the transmission of internal waves through such a region. The free modes depend strongly on the staircase properties, and consist of modes with both internal and interfacial gravity wave-like behaviour. We determine the frequency shifts of these waves as a function of the number of steps to explore their potential to probe planetary internal structures. We also find that wave transmission is strongly affected by the presence of a staircase. Very large-wavelength waves are transmitted efficiently, but small-scale waves are only transmitted if they are resonant with one of the free modes. The effective size of the core is therefore larger for non-resonant modes.

We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.

Convection is thought to act as a turbulent viscosity in damping tidal flows and in driving spin and orbital evolution in close convective binary systems. This turbulent viscosity should be reduced, compared to mixing-length predictions, when the forcing (tidal) frequency $|{\omega}_t|$ exceeds the turnover frequency ${\omega}_{cv}$ of the dominant convective eddies. However, two contradictory scaling laws have been proposed and this issue remains highly disputed. To revisit this controversy, we conduct the first direct numerical simulations (DNS) of convection interacting with the equilibrium tidal flow in an idealized global model of a low-mass star. We present direct computations of the turbulent effective viscosity, ${\nu}_E$, acting on the equilibrium tidal flow. We unexpectedly report the coexistence of the two disputed scaling laws, which reconciles previous theoretical (and numerical) findings. We recover the universal quadratic scaling ${\nu}_E \propto (|{\omega}_t|/{\omega}_{cv})^{-2}$ in the high-frequency regime $|{\omega}_t|/{\omega}_{cv} \gg 1$. Our results also support the linear scaling ${\nu}_E \propto (|{\omega}_t|/{\omega}_{cv})^{-1}$ in an intermediate regime with $1 \leq |{\omega}_t|/{\omega}_{cv} = \mathcal{O}(10)$. Both regimes may be relevant to explain the observed properties of close binaries, including spin synchronization of solar-type stars and the circularization of low-mass stars. The robustness of these two regimes of tidal dissipation, and the transition between them, should be explored further in more realistic models. A better understanding of the interaction between convection and tidal flows is indeed essential to correctly interpret observations of close binary stars and short-period planetary orbits.

Tidal interactions are important in driving spin and orbital evolution in planetary and stellar binary systems, but the fluid dynamical mechanisms responsible remain incompletely understood. One key mechanism is the interaction between tidal flows and convection. Turbulent convection is thought to act as an effective viscosity in damping large-scale tidal flows, but there is a long-standing controversy over the efficiency of this mechanism when the tidal frequency exceeds the turnover frequency of the dominant convective eddies. This high frequency regime is relevant for many applications, such as for tides in stars hosting hot Jupiters. We explore the interaction between tidal flows and convection using hydrodynamical simulations within a local Cartesian model of a small patch of a convection zone of a star or planet. We adopt the Boussinesq approximation and simulate Rayleigh-B\'enard convection, modelling the tidal flow as a background oscillatory shear flow. We demonstrate that the effective viscosity of both laminar and turbulent convection is approximately frequency-independent for low frequencies. When the forcing frequency exceeds the dominant convective frequency, the effective viscosity scales inversely with the square of the tidal frequency. We also show that negative effective viscosities are possible, particularly for high frequency tidal forcing, suggesting the surprising possibility of tidal anti-dissipation. These results are supported by a complementary high frequency asymptotic analysis that extends prior work by Ogilvie & Lesur. We discuss the implications of these results for interpreting the orbital decay of hot Jupiters, and for several other astrophysical problems.

We present an investigation into the nonlinear evolution of the Goldreich-Schubert-Fricke (GSF) instability using axisymmetric and three-dimensional simulations near the equator of a differentially rotating radiation zone. This instability may provide an important contribution to angular momentum transport in stars and planets. We adopt a local Boussinesq Cartesian shearing box model, which represents a small patch of a differentially rotating stellar radiation zone. Complementary simulations are also performed with stress-free, impenetrable boundaries in the local radial direction. The linear and nonlinear evolution of the equatorial axisymmetric instability is formally equivalent to the salt fingering instability. This is no longer the case in 3D, but we find that the instability behaves nonlinearly in a similar way to salt fingering. Axisymmetric simulations -- and those in 3D with short dimensions along the local azimuthal direction -- quickly develop strong jets along the rotation axis, which inhibit the instability and lead to predator-prey-like temporal dynamics. In 3D, the instability initially produces homogeneous turbulence and enhanced momentum transport, though in some cases jets form on a much longer timescale. We propose and validate numerically a simple theory for nonlinear saturation of the GSF instability and its resulting angular momentum transport. This theory is straightforward to implement in stellar evolution codes incorporating rotation. We estimate that the GSF instability could contribute towards explaining the missing angular momentum transport required in red giant stars, and play a role in the long-term evolution of the solar tachocline.

Recent Juno observations have suggested that the heavy elements in Jupiter could be diluted throughout a large fraction of its gaseous envelope, providing a stabilising compositional gradient over an extended region of the planet. This could trigger layered semi-convection, which, in the context of giant planets more generally, may explain Saturn's luminosity excess and play a role in causing the abnormally large radii of some hot Jupiters. In giant planet interiors, it could take the form of density staircases, which are convective layers separated by thin stably stratified interfaces. In addition, the efficiency of tidal dissipation is known to depend strongly on the planetary internal structure. We aim to study the resulting tidal dissipation when internal waves are excited in a region of layered semi-convection by tidal gravitational forcing due to other bodies (such as moons in giant planet systems, or stars in hot Jupiter systems). We adopt a local Cartesian model with a background layered density profile subjected to an imposed tidal forcing, and we compute the viscous and thermal dissipation rates numerically. We find that the rates of tidal dissipation can be enhanced in a region of layered semi-convection compared to a uniformly convective medium, where the latter corresponds with the usual assumption adopted in giant planet interior models. In particular, a region of layered semi-convection possesses a richer set of resonances, allowing enhanced dissipation for a wider range of tidal frequencies. The details of these results significantly depend on the structural properties of the layered semi-convective regions. Layered semi-convection could contribute towards explaining the high tidal dissipation rates observed in Jupiter and Saturn, which have not yet been fully explained by theory. Further work is required to explore the efficiency of this mechanism in global models.

The Kelvin-Helmholtz instability has been proposed as a mechanism to extract energy from magnetohydrodynamic (MHD) kink waves in flux tubes, and to drive dissipation of this wave energy through turbulence. It is therefore a potentially important process in heating the solar corona. However, it is unclear how the instability is influenced by the oscillatory shear flow associated with an MHD wave. We investigate the linear stability of a discontinuous oscillatory shear flow in the presence of a horizontal magnetic field within a Cartesian framework that captures the essential features of MHD oscillations in flux tubes. We derive a Mathieu equation for the Lagrangian displacement of the interface and analyse its properties, identifying two different instabilities: a Kelvin-Helmholtz instability and a parametric instability involving resonance between the oscillatory shear flow and two surface Alfv\'{e}n waves. The latter occurs when the system is Kelvin-Helmholtz stable, thus favouring modes that vary along the flux tube, and as a consequence provides an important and additional mechanism to extract energy. When applied to flows with the characteristic properties of kink waves in the solar corona, both instabilities can grow, with the parametric instability capable of generating smaller scale disturbances along the magnetic field than possible via the Kelvin-Helmholtz instability. The characteristic time-scale for these instabilities is $\sim 100$ s, for wavelengths of $200$ km. The parametric instability is more likely to occur for smaller density contrasts and larger velocity shears, making its development more likely on coronal loops than on prominence threads.

Layered semi-convection could operate in giant planets, potentially explaining the constraints on the heavy elements distribution in Jupiter deduced recently from Juno observations, and contributing to Saturn's luminosity excess or the abnormally large radius of some hot Jupiters. This is a state consisting of density staircases, in which convective layers are separated by thin stably stratified interfaces. The efficiency of tidal dissipation in a planet depends strongly on its internal structure. It is crucial to improve our understanding of the mechanisms driving this dissipation, since it has important consequences to predict the long-term evolution of any planetary system. In this work, our goal is to study the resulting tidal dissipation when internal waves are excited by other bodies (such as the moons of giant planets) in a region of layered semi-convection. We find that the rates of tidal dissipation can be significantly enhanced in a layered semi-convective medium compared to a uniformly convective medium, especially in the astrophysically relevant sub-inertial frequency range. Thus, layered semi-convection is a possible candidate to explain high tidal dissipation rates recently observed in Jupiter and Saturn.

Layered semi-convection is a possible candidate to explain Saturn's luminosity excess and the abnormally large radius of some hot Jupiters. In giant planet interiors, it could lead to the creation of density staircases, which are convective layers separated by thin stably stratified interfaces. We study the propagation of internal waves in a region of layered semi-convection, with the aim to predict energy transport by internal waves incident upon a density staircase. The goal is then to understand the resulting tidal dissipation when these waves are excited by other bodies such as moons in giant planets systems. We use a local Cartesian analytical model, taking into account the complete Coriolis acceleration at any latitude, thus generalizing previous works. We find transmission of incident internal waves to be strongly affected by the presence of a density staircase, even if these waves are initially pure inertial waves (which are restored by the Coriolis acceleration). In particular, low-frequency waves of all wavelengths are perfectly transmitted near the critical latitude. Otherwise, short-wavelength waves are only efficiently transmitted if they are resonant with a free mode (interfacial gravity wave or short-wavelength inertial mode) of the staircase. In all other cases, waves are primarily reflected unless their wavelengths are longer than the vertical extent of the entire staircase (not just a single step). We expect incident internal waves to be strongly affected by the presence of a density staircase in a frequency-, latitude- and wavelength-dependent manner. First, this could lead to new criteria to probe the interior of giant planets by seismology; and second, this may have important consequences for tidal dissipation and our understanding of the evolution of giant planet systems.

I discuss two related nonlinear mechanisms of tidal dissipation that require finite tidal deformations for their operation: the elliptical instability and the precessional instability. Both are likely to be important for the tidal evolution of short-period extrasolar planets. The elliptical instability is a fluid instability of elliptical streamlines, such as in tidally deformed non-synchronously rotating or non-circularly orbiting planets. I summarise the results of local and global simulations that indicate this mechanism to be important for tidal spin synchronisation, planetary spin-orbit alignment and orbital circularisation for the shortest period hot Jupiters. The precessional instability is a fluid instability that occurs in planets undergoing axial precession, such as those with spin-orbit misalignments (non-zero obliquities). I summarise the outcome of local MHD simulations designed to study the turbulent damping of axial precession, which suggest this mechanism to be important in driving tidal evolution of the spin-orbit angle for hot Jupiters. Avenues for future work are also discussed.

The WASP-18 system, with its massive and extremely close-in planet, WASP-18b (M$_{p}$ = 10.3M$_{J}$, a = 0.02 AU, P = 22.6 hours), is one of the best known exoplanet laboratories to directly measure Q', the modified tidal quality factor and proxy for efficiency of tidal dissipation, of the host star. Previous analysis predicted a rapid orbital decay of the planet toward its host star that should be measurable on the time scale of a few years, if the star is as dissipative as is inferred from the circularization of close-in solar-type binary stars. We have compiled published transit and secondary eclipse timing (as observed by WASP, TRAPPIST, and Spitzer) with more recent unpublished light curves (as observed by TRAPPIST and HST) with coverage spanning nine years. We find no signature of a rapid decay. We conclude that the absence of rapid orbital decay most likely derives from Q' being larger than was inferred from solar-type stars, and find that Q' $\geq $1$\times10^{6}$, at 95\,\% confidence; this supports previous work suggesting that F-stars, with their convective cores and thin convective envelopes, are significantly less tidally dissipative than solar-type stars, with radiative cores and large convective envelopes.

The spin axis of a rotationally deformed planet is forced to precess about its orbital angular momentum vector, due to the tidal gravity of its host star, if these directions are misaligned. This induces internal fluid motions inside the planet that are subject to a hydrodynamic instability. We study the turbulent damping of precessional fluid motions, as a result of this instability, in the simplest local computational model of a giant planet (or star), with and without a weak internal magnetic field. Our aim is to determine the outcome of this instability, and its importance in driving tidal evolution of the spin-orbit angle in precessing planets (and stars). We find that this instability produces turbulent dissipation that is sufficiently strong that it could drive significant tidal evolution of the spin-orbit angle for hot Jupiters with orbital periods shorter than about 10-18 days. If this mechanism acts in isolation, this evolution would be towards alignment or anti-alignment, depending on the initial angle, but the ultimate evolution (if other tidal mechanisms also contribute) is expected to be towards alignment. The turbulent dissipation is proportional to the cube of the precession frequency, so it leads to much slower damping of stellar spin-orbit angles, implying that this instability is unlikely to drive evolution of the spin-orbit angle in stars (either in planetary or close binary systems). We also find that the instability-driven flow can act as a system-scale dynamo, which may play a role in producing magnetic fields in short-period planets.

I present results from the first global hydrodynamical simulations of the elliptical instability in a tidally deformed gaseous planet (or star) with a free surface. The elliptical instability is potentially important for tidal evolution of the shortest-period hot Jupiters. I model the planet as a spin-orbit aligned or anti-aligned, and non-synchronously rotating, tidally deformed, homogeneous fluid body. A companion paper presented an analysis of the global modes and instabilities of such a planet. Here I focus on the nonlinear evolution of the elliptical instability. This is observed to produce bursts of turbulence that drive the planet towards synchronism with its orbit in an erratic manner. If the planetary spin is initially anti-aligned, the elliptical instability also drives spin-orbit alignment on a similar timescale as the spin synchronisation. The instability generates differential rotation inside the planet in the form of zonal flows, which play an important role in the saturation of the instability, and in producing the observed burstiness. These results are broadly consistent with the picture obtained using a local Cartesian model (where columnar vortices played the role of zonal flows). I also simulate the instability in a container that is rigid (but stress-free) rather than free, finding broad quantitative agreement. The dissipation resulting from the elliptical instability could explain why the shortest-period hot Jupiters tend to have circular orbits inside about 2-3 days, and predicts spin-synchronisation (and spin-orbit alignment) out to about 10-15 days. However, other mechanisms must be invoked to explain tidal circularisation for longer orbital periods.

We revisit the global modes and instabilities of homogeneous rotating ellipsoidal fluid masses, which are the simplest global models of rotationally and tidally deformed gaseous planets or stars. The tidal flow in a short-period planet may be unstable to the elliptical instability, a hydrodynamic instability that can drive tidal evolution. We perform a global (and local WKB) analysis to study this instability using the elegant formalism of Lebovitz & Lifschitz. We survey the parameter space of global instabilities with harmonic orders $\ell\leq 5$, for planets with spins that are purely aligned (prograde) or anti-aligned (retrograde) with their orbits. In general, the instability has a much larger growth rate if the planetary spin and orbit are anti-aligned rather than aligned. We have identified a violent instability for anti-aligned spins outside of the usual frequency range for the elliptical instability (when $\frac{n}{\Omega}\lesssim -1$, where $n$ and $\Omega$ are the orbital and spin angular frequencies, respectively) if the tidal amplitude is sufficiently large. We also explore the instability in a rigid ellipsoidal container, which is found to be quantitatively similar to that with a realistic free surface. Finally, we study the effect of rotation and tidal deformation on mode frequencies. We find that larger rotation rates and larger tidal deformations both decrease the frequencies of the prograde sectoral surface gravity modes. This increases the prospect of their tidal excitation, potentially enhancing the tidal response over expectations from linear theory. In a companion paper, we use our results to interpret global simulations of the elliptical instability.

We perform global two-dimensional hydrodynamical simulations of Keplerian discs with free eccentricity over thousands of orbital periods. Our aim is to determine the validity of secular theory in describing the evolution of eccentric discs, and to explore their nonlinear evolution for moderate eccentricities. Linear secular theory is found to correctly predict the structure and precession rates of discs with small eccentricities. However, discs with larger eccentricities (and eccentricity gradients) are observed to precess faster (retrograde relative to the orbital motion), at a rate that depends on their eccentricities (and eccentricity gradients). We derive analytically a nonlinear secular theory for eccentric gas discs, which explains this result as a modification of the pressure forces whenever eccentric orbits in a disc nearly intersect. This effect could be particularly important for highly eccentric discs produced in tidal disruption events, or for narrow gaseous rings; it might also play a role in causing some of the variability in superhump binary systems. In two dimensions, the eccentricity of a moderately eccentric disc is long-lived and persists throughout the duration of our simulations. Eccentric modes are however weakly damped by their interaction with non-axisymmetric spiral density waves (driven by the Papaloizou-Pringle instability, which occurs in our idealised setup with solid walls), as well as numerical diffusion.

We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a "vertical-shear instability", which recent nonlinear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. (2013). This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power law radial entropy profiles ("locally polytropic discs"), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. In all cases the fastest growing modes have very short wavelengths and are localised at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilised classical inertial waves ("r-modes"). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the nonlinear saturation and resulting angular momentum transport.

Eccentric Keplerian discs are believed to be unstable to three-dimensional hydrodynamical instabilities driven by the time-dependence of fluid properties around an orbit. These instabilities could lead to small-scale turbulence, and ultimately modify the global disc properties. We use a local model of an eccentric disc, derived in a companion paper, to compute the nonlinear vertical ("breathing mode") oscillations of the disc. We then analyse their linear stability to locally axisymmetric disturbances for any disc eccentricity and eccentricity gradient using a numerical Floquet method. In the limit of small departures from a circular reference orbit, the instability of an isothermal disc is explained analytically. We also study analytically the small-scale instability of an eccentric neutrally stratified polytropic disc with any polytropic index using a WKB approximation. We find that eccentric discs are generically unstable to the parametric excitation of small-scale inertial waves. The nonlinear evolution of these instabilities should be studied in numerical simulations, where we expect them to lead to a decay of the disc eccentricity and eccentricity gradient as well as to induce additional transport and mixing. Our results highlight that it is essential to consider the three-dimensional structure of eccentric discs, and their resulting vertical oscillatory flows, in order to correctly capture their evolution.

We formulate a local dynamical model of an eccentric disc in which the dominant motion consists of elliptical Keplerian orbits. The model is a generalization of the well known shearing sheet, and is suitable for both analytical and computational studies of the local dynamics of eccentric discs. It is spatially homogeneous in the horizontal dimensions but has a time-dependent geometry that oscillates at the orbital frequency. We show how certain averages of the stress tensor in the local model determine the large-scale evolution of the shape and mass distribution of the disc. The simplest solutions of the local model are laminar flows consisting of a (generally nonlinear) vertical oscillation of the disc. Eccentric discs lack vertical hydrostatic equilibrium because of the variation of the vertical gravitational acceleration around the eccentric orbit, and in some cases because of the divergence of the orbital velocity field associated with an eccentricity gradient. We discuss the properties of the laminar solutions, showing that they can exhibit extreme compressional behaviour for eccentricities greater than about $0.5$, especially in discs that behave isothermally. We also derive the linear evolutionary equations for an eccentric disc that follow from the laminar flows in the absence of a shear viscosity. In a companion paper we show that these solutions are linearly unstable and we determine the associated growth rates and unstable modes.

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield $0.75\pm 0.06$. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.

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.

We investigate whether the elliptical instability is important for tidal dissipation in gaseous planets and stars. In a companion paper, we found that the conventional elliptical instability results in insufficient dissipation because it produces long-lived vortices that then quench further instability. Here, we study whether the addition of a magnetic field prevents those vortices from forming, and hence leads to enhanced dissipation. We present results from magnetohydrodynamical simulations that evolve the elliptical instability in a local patch of a rotating planet or star, in the presence of a weak magnetic field. We find that magnetic fields do indeed prevent vortices from forming, and hence greatly enhance the steady state dissipation rate. In addition, the resulting turbulence acts as a small-scale dynamo, amplifying the initially weak field. The inferred tidal dissipation is potentially important at short orbital periods. For example, it can circularise hot Jupiters with orbital periods shorter than 2.5 days, and synchronise their spins with their orbits out to 6 days. However, it appears unable to account for the hot Jupiters that appear to have been circularised out to six to ten day orbital periods. It also cannot account for the inferred circularisation of many close binary stars.

Tidally distorted rotating stars and gaseous planets are subject to a well-known linear fluid instability -- the elliptical instability. It has been proposed that this instability might drive enough energy dissipation to solve the long-standing problem of the origin of tidal dissipation in stars and planets. But the nonlinear outcome of the elliptical instability has yet to be investigated in the parameter regime of interest, and the resulting turbulent energy dissipation has not yet been quantified. We do so by performing three dimensional hydrodynamical simulations of a small patch of a tidally deformed fluid planet or star subject to the elliptical instability. We show that when the tidal deformation is weak, the nonlinear outcome of the instability leads to the formation of long-lived columnar vortices aligned with the axis of rotation. These vortices shut off the elliptical instability, and the net result is insufficient energy dissipation to account for tidal dissipation. However, further work is required to account for effects neglected here, including magnetic fields, turbulent convection, and realistic boundary conditions.

We perform idealised numerical simulations of magnetic buoyancy instabilities in a model of the solar tachocline. We introduce a simplified model of magnetic flux pumping in an upper layer (the convection zone), and study the effects of its inclusion on the evolution of buoyancy instabilities in a lower layer (the radiative interior). We study its effects on the instability of both a preconceived magnetic slab and of a shear-generated magnetic layer. In the former, we find that in the regime in which the downward pumping velocity is comparable with the Alfven speed of the magnetic layer, flux pumping is able to hold back the bulk of the magnetic field, with only small pockets of strong field able to rise into the upper layer. In simulations in which the magnetic layer is generated by shear, we find that the shear velocity is not necessarily required to exceed that of the pumping (therefore the kinetic energy of the shear is not required to exceed that of the overlying convection), for strong localised pockets of magnetic field to be produced which can rise into the upper layer. This is because magnetic flux pumping acts to store the field below the interface, allowing it to be amplified both by the shear, and by vortical fluid motions, until pockets of field can achieve sufficient strength to rise into the upper layer. In addition, we find that the interface between the two layers is a natural location for the production of strong vertical gradients in the magnetic field. If these gradients are sufficiently strong to allow the development of magnetic buoyancy instabilities, strong shear is not necessarily required to drive them (c.f. previous work by Vasil & Brummell). We find that the addition of magnetic flux pumping appears to be able to assist shear-driven magnetic buoyancy in producing strong flux concentrations that can rise into the convection zone from the radiative interior.

We perform a stability analysis of a tidally excited nonlinear internal gravity wave near the centre of a solar-type star in two-dimensions. The motivation is to understand the tidal interaction between short-period planets and their solar-type host stars, which involves the launching of gravity waves at the top of the radiation zone that propagate towards the stellar centre. Studying the instabilities of these waves near the centre, where nonlinearities are most important, is essential, since it may have implications for the survival of these planets. When the waves have sufficient amplitude to overturn the stratification, they break and form a critical layer, which efficiently absorbs subsequent ingoing wave angular momentum, and can result in the planet spiralling into the star. However, previous simulations do not find the waves to undergo instability for smaller amplitudes. This work has two aims: to determine any instabilities that set in for small-amplitude waves, and to further understand the breaking of large-amplitude waves. Our main result is that the waves undergo parametric instabilities for any amplitude. However, because the nonlinearity is spatially localised in the innermost wavelengths, their growth rates are sufficiently small that they do not result in astrophysically important tidal dissipation. The resulting modified tidal quality factors are estimated to be Q'_star>10^7, and possibly much greater, so the dissipation is much weaker than that which results from critical-layer absorption. These results support our explanation for the survival of all currently observed short-period planets around solar-type main-sequence stars: that planets unable to cause wave breaking at the centre of their host stars are likely to survive against tidal decay. This hypothesis will be tested by ongoing and future observations of transiting planets, such as WASP and Kepler.

Internal gravity waves are excited at the interface of convection and radiation zones of a solar-type star by the tidal forcing of a short-period planet. The fate of these waves as they approach the centre of the star depends on their amplitude. We discuss the results of numerical simulations of these waves approaching the centre of a star, and the resulting evolution of the spin of the central regions of the star, and the orbit of the planet. If the waves break, we find efficient tidal dissipation, which is not present if the waves perfectly reflect from the centre. This highlights an important amplitude dependence of the (stellar) tidal quality factor Q', which has implications for the survival of planets on short-period orbits around solar-type stars, with radiative cores.

We study the fate of internal gravity waves approaching the centre of a non-rotating solar-type star, by performing 3D numerical simulations using a Boussinesq-type model. These waves are excited at the top of the radiation zone by the tidal forcing of a short-period planet on a circular, coplanar orbit. This extends previous work done in 2D by Barker & Ogilvie. We first derive a linear wave solution, which is not exact in 3D; however, the reflection of ingoing waves from the centre is close to perfect for moderate amplitude waves. Waves with sufficient amplitude to cause isentropic overturning break, and deposit their angular momentum near the centre. This forms a critical layer, at which the angular velocity of the flow matches the orbital angular frequency of the planet. This efficiently absorbs ingoing waves, and spins up the star from the inside out, while the planet spirals into the star. We also perform numerical integrations to determine the linearised adiabatic tidal response throughout the star, for solar-type stellar models with masses in the range 0.5 < m_star/M_sun < 1.1, throughout their main sequence lifetimes. The aim is to study the launching region for these waves at the top of the radiation zone in more detail, and to determine the accuracy of a semi-analytic approximation for the tidal torque on the star, that was derived under the assumption that all ingoing wave angular momentum is absorbed in a critical layer. The main conclusions are that this nonlinear mechanism of tidal dissipation could provide an explanation for the survival of all short-period extrasolar planets observed around FGK stars, while it predicts the destruction of more massive planets. This work provides further support for the model outlined in a previous paper by Barker & Ogilvie, and makes predictions that will be tested by ongoing observational studies, such as WASP and Kepler.

We study the fate of internal gravity waves, which are excited by tidal forcing by a short-period planet at the interface of convection and radiation zones, approaching the centre of a solar-type star. We study at what amplitude these wave are subject to instabilities. These instabilities lead to wave breaking whenever the amplitude exceeds a critical value. Below this value, the wave reflects perfectly from the centre of the star. Wave breaking results in spinning up the central regions of the star, and the formation of a critical layer, which acts as an absorbing barrier for ingoing waves. As these waves are absorbed, the star is spun up from the inside out. This results in an important amplitude dependence of the tidal quality factor Q'. If the tidal forcing amplitude exceeds the value required for wave breaking, efficient dissipation results over a continuous range of tidal frequencies, leading to Q' \approx 10^5 (P/1day)^(8/3), for the current Sun. This varies by less than a factor of 5 throughout the range of G and K type main sequence stars, for a given orbit. We predict fewer giant planets with orbital periods of less than about 2 days around such stars, if they cause breaking at the centre, due to the efficiency of this process. This mechanism would, however, be ineffective in stars with a convective core, such as WASP-18, WASP-12 and OGLE-TR-56, perhaps partly explaining the survival of their close planetary companions.

Tidal friction is thought to be important in determining the long-term spin-orbit evolution of short-period extrasolar planetary systems. Using a simple model of the orbit-averaged effects of tidal friction, we study the evolution of close-in planets on inclined orbits, due to tides. We analyse the effects of the inclusion of stellar magnetic braking by performing a phase-plane analysis of a simplified system of equations, including the braking torque. The inclusion of magnetic braking is found to be important, and its neglect can result in a very different system history. We then present the results of numerical integrations of the tidal evolution equations, where we find that it is essential to consider coupled evolution of the orbital and rotational elements, including dissipation in both the star and planet, to accurately model the evolution. The main result of our integrations is that for typical Hot Jupiters, tidal friction aligns the stellar spin with the orbit on a similar time as it causes the orbit to decay. This means that if a planet is observed to be aligned, then it probably formed coplanar. This reinforces the importance of Rossiter-McLaughlin effect observations in determining the degree of spin-orbit alignment in transiting systems. We apply these results to the XO-3 system, and constrain the tidal quality factors Q' in both the star and planet in this system. Using a model in which inertial waves are excited by tidal forcing in the outer convective envelope and dissipated by turbulent viscosity, we calculate Q' for a range of F-star models, and find it to vary considerably within this class of stars. This means that assuming a single Q' applies to all stars is probably incorrect. We propose an explanation for the survival of WASP-12 b & OGLE-TR-56 b, in terms of weak dissipation in the star.

Tidal friction is thought to be important in determining the long-term spin-orbit evolution of short-period extrasolar planetary systems. Using a simple model of the orbit-averaged effects of tidal friction Eggleton, Kiseleva & Hut (1998), we analyse the effects of the inclusion of stellar magnetic braking on the evolution of such systems. A phase-plane analysis of a simplified system of equations, including only the stellar tide together with a model of the braking torque proposed by Verbunt & Zwaan (1981), is presented. The inclusion of stellar magnetic braking is found to be extremely important in determining the secular evolution of such systems, and its neglect results in a very different orbital history. We then show the results of numerical integrations of the full tidal evolution equations, using the misaligned spin and orbit of the XO-3 system as an example, to study the accuracy of simple timescale estimates of tidal evolution. We find that it is essential to consider coupled evolution of the orbit and the stellar spin in order to model the behaviour accurately. In addition, we find that for typical Hot Jupiters the stellar spin-orbit alignment timescale is of the same order as the inspiral time, which tells us that if a planet is observed to be aligned, then it probably formed coplanar. This reinforces the importance of Rossiter-McLaughlin effect observations in determining the degree of spin-orbit alignment in transiting systems.

We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar diffusion then plays an important role in the formation process. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D numerical simulations, we follow the evolution of the fast-mode wave. The simulations are done with a multifluid, adaptive mesh refinement MHD code. Initial perturbations with wavelengths that are 2 orders of magnitude larger than the dissipation length are strongly affected by the ion-neutral drift. Only in situations where there are large variations in the magnetic field corresponding to a highly turbulent gas can fast-mode waves generate dense cores. This means that, in most cores, no substructure can be produced. However, Core D of TMC-1 is an exception to this case. Due to its atypically high ionisation fraction, waves with wavelengths up to 3 orders of magnitude greater than the dissipation length can be present. Such waves are only weakly affected by ambipolar diffusion and can produce dense substructure without large wave-amplitudes. Our results also explain the observed transition from Alfvenic turbulent motion at large scales to subsonic motions at the level of dense cores.