Adrian J Barker
University of Leeds, Department of Applied Mathematics, Faculty Member
Research Interests:
In close exoplanetary systems, tidal interactions drive orbital and spin evolution of planets and stars over long time-scales. Tidally forced inertial waves (restored by the Coriolis acceleration) in the convective envelopes of low-mass... more
In close exoplanetary systems, tidal interactions drive orbital and spin evolution of planets and stars over long time-scales. Tidally forced inertial waves (restored by the Coriolis acceleration) in the convective envelopes of low-mass stars and giant gaseous planets contribute greatly to the tidal dissipation when they are excited and subsequently damped (e.g. through viscous friction), especially early in the life of a system. These waves are known to be subject to non-linear effects, including triggering differential rotation in the form of zonal flows. In this study, we use a realistic tidal body forcing to excite inertial waves through the residual action of the equilibrium tide in the momentum equation for the waves. By performing 3D non-linear hydrodynamical simulations in adiabatic and incompressible convective shells, we investigate how the addition of non-linear terms affects the tidal flow properties, and the energy and angular momentum redistribution. In particular, we ...
Research Interests:
Tidal dissipation is responsible for circularizing the orbits and synchronizing the spins of solar-type close binary stars, but the mechanisms responsible are not fully understood. Previous work has indicated that significant enhancements... more
Tidal dissipation is responsible for circularizing the orbits and synchronizing the spins of solar-type close binary stars, but the mechanisms responsible are not fully understood. Previous work has indicated that significant enhancements to the theoretically predicted tidal dissipation rates are required to explain the observed circularization periods (P circ) in various stellar populations and their evolution with age. This was based partly on the common belief that the dominant mechanism of tidal dissipation in solar-type stars is turbulent viscosity acting on equilibrium tides in convective envelopes. In this paper, we study tidal dissipation in both convection and radiation zones of rotating solar-type stars following their evolution. We study equilibrium tide dissipation, incorporating a frequency-dependent effective viscosity motivated by the latest hydrodynamical simulations, and inertial wave (dynamical tide) dissipation, adopting a frequency-averaged formalism that account...
Research Interests: Physics and Dissipation
Research Interests:
We study thermal convection in a rotating fluid, with the ultimate goal of explaining the structure of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the... more
We study thermal convection in a rotating fluid, with the ultimate goal of explaining the structure 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. 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 ± 0.06. We conclude that the mixing...
Research Interests:
Abstract. 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... more
Abstract. 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.
Research Interests:
Research Interests:
Research Interests:
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... more
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,...
Research Interests:
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... more
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 sat...
Research Interests:
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... more
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' ≥ 1×10^6, at 95 % confidence;...
Research Interests:
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... more
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 ver...
Research Interests:
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... more
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.7...
Research Interests:
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... more
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 mi...
Research Interests:
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... more
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 pr...
Research Interests:
We study the fate of internal gravity waves approaching the centre of an initially non-rotating solar-type star, primarily using two-dimensional numerical simulations based on a cylindrical model. A train of internal gravity waves is... more
We study the fate of internal gravity waves approaching the centre of an initially non-rotating solar-type star, primarily using two-dimensional numerical simulations based on a cylindrical model. A train of internal gravity waves is excited by tidal forcing at the interface between the convection and radiation zones of such a star. We derive a Boussinesq-type model of the central region of a star and find a nonlinear wave solution that is steady in the frame rotating with the angular pattern speed of the tidal forcing. We then use spectral methods to integrate the equations numerically, with the aim of studying at what amplitude the wave is subject to instabilities. These instabilities are found to lead to wave breaking whenever the amplitude exceeds a critical value. Below this critical value, the wave reflects perfectly from the centre of the star. Wave breaking leads to mean flow acceleration, which corresponds to a spin up of the central region of the star, and the formation of...
Eccentric Keplerian discs are believed to be unstable to three-dimensional hydro-dynamical instabilities driven by the time-dependence of fluid properties around an orbit. These instabilities could lead to small-scale turbulence, and... more
Eccentric Keplerian discs are believed to be unstable to three-dimensional hydro-dynamical 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 in-dex using a WKB approximation. We find that eccentric discs are generically unstable to the parametric excitation of small-scale inertial wav...
Research Interests:
Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilize ordinary convection and produce a stably stratified layer near the core... more
Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilize ordinary convection and produce a stably stratified layer near the core of these planets. This region could consist of semiconvective 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 semiconvective layers. We then analyse the transmis...
Research Interests:
We study tidal dissipation in stars with masses in the range 0.1–1.6 M⊙ throughout their evolution, including turbulent effective viscosity acting on equilibrium tides and inertial waves (IWs) in convection zones, and internal gravity... more
We study tidal dissipation in stars with masses in the range 0.1–1.6 M⊙ throughout their evolution, including turbulent effective viscosity acting on equilibrium tides and inertial waves (IWs) in convection zones, and internal gravity waves in radiation zones. We consider a range of stellar evolutionary models and incorporate the frequency-dependent effective viscosity acting on equilibrium tides based on the latest simulations. We compare the tidal flow and dissipation obtained with the conventional equilibrium tide, which is strictly invalid in convection zones, finding that the latter typically overpredicts the dissipation by a factor of 2–3. Dissipation of IWs is computed using a frequency-averaged formalism accounting for realistic stellar structure for the first time, and is the dominant mechanism for binary circularization and synchronization on the main sequence. Dissipation of gravity waves in the radiation zone assumes these waves to be fully damped (e.g. by wave breaking)...
Research Interests:
We continue our investigation into the non-linear 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... more
We continue our investigation into the non-linear 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 non-linear 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 non-linear 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 mom...
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... more
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 ℓ≤ 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 n/Ω≲ -1, where n and Ω are the orbital and spin angular frequencies, respectiv...
Research Interests:
The interaction between equilibrium tides and convection in stellar envelopes is often considered important for tidal evolution in close binary and extrasolar planetary systems. Its efficiency for fast tides has however long been... more
The interaction between equilibrium tides and convection in stellar envelopes is often considered important for tidal evolution in close binary and extrasolar planetary systems. Its efficiency for fast tides has however long been controversial, when the tidal frequency exceeds the turnover frequency of convective eddies. Recent numerical simulations indicate that convection can act like an effective viscosity which decays quadratically with tidal frequency for fast tides, resulting in inefficient dissipation in many applications involving pre- and main-sequence stars and giant planets. A new idea was however recently proposed by Terquem (2021), who suggested Reynolds stresses involving correlations between tidal flow components dominate the interaction instead of correlations between convective flow components as usually assumed. They further showed that this can potentially significantly enhance tidal dissipation for fast tides in many applications. Motivated by the importance of t...
Research Interests:
Research Interests:
Research Interests:
Research Interests:
Research Interests:
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... more
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 impor...
Research Interests:
Research Interests:
Research Interests:
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... more
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énard convection, modelling the tidal flow as a background oscillatory shear flow. We demonstrate that th...
Turbulent convection is thought to act as an effective viscosity in damping equilibrium tidal flows, driving spin and orbital evolution in close convective binary systems. Compared to mixing-length predictions, this viscosity ought to be... more
Turbulent convection is thought to act as an effective viscosity in damping equilibrium tidal flows, driving spin and orbital evolution in close convective binary systems. Compared to mixing-length predictions, this viscosity ought to be reduced when the tidal frequency |ωt| exceeds the turnover frequency ωcv of the dominant convective eddies, but the efficiency of this reduction has been disputed. We re-examine this long-standing controversy using direct numerical simulations of an idealized global model. We simulate thermal convection in a full sphere, and externally forced by the equilibrium tidal flow, to measure the effective viscosity νE acting on the tidal flow when |ωt|/ωcv ≳ 1. We demonstrate that the frequency reduction of νE is correlated with the frequency spectrum of the (unperturbed) convection. For intermediate frequencies below those in the turbulent cascade (|ωt|/ωcv ∼ 1−5), the frequency spectrum displays an anomalous 1/ωα power law that is responsible for the freq...
Turbulent convection is thought to act as an effective viscosity (νE) in damping tidal flows in stars and giant planets. However, the efficiency of this mechanism has long been debated, particularly in the regime of fast tides, when the... more
Turbulent convection is thought to act as an effective viscosity (νE) in damping tidal flows in stars and giant planets. However, the efficiency of this mechanism has long been debated, particularly in the regime of fast tides, when the tidal frequency (ω) exceeds the turnover frequency of the dominant convective eddies (ωc). We present the results of hydrodynamical simulations to study the interaction between tidal flows and convection in a small patch of a convection zone. These simulations build upon our prior work by simulating more turbulent convection in larger horizontal boxes, and here we explore a wider range of parameters. We obtain several new results: (1) νE is frequency dependent, scaling as ω−0.5 when ω/ωc ≲ 1, and appears to attain its maximum constant value only for very small frequencies (ω/ωc ≲ 10−2). This frequency reduction for low-frequency tidal forcing has never been observed previously. (2) The frequency dependence of νE appears to follow the same scaling as ...
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... more
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 account...
Research Interests:
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... more
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.
Research Interests:
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... more
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.