BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gordon I. Ogilvie

Gordon I. Ogilvie appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.HEphysics.flu-dynastro-phastro-ph.GAnlin.PS

Catalog footprint

What is connected

33works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

Hydrodynamics of warps in the local model of astrophysical discs

We show how the local approximation of astrophysical discs, which is the basis for the well known model of the shearing box, can be used to study many aspects of the dynamics of warped discs. In the local model, inclination of the orbit of a test particle with respect to the reference orbit corresponds to a vertical oscillation of the particle at the orbital frequency. Warping of a disc corresponds to a locally axisymmetric corrugation of the midplane of the disc that oscillates vertically at the orbital frequency, while evolution of the warp corresponds to a modulation of the complex amplitude of the vertical oscillation. We derive a conservation law for this amplitude that is the local equivalent of the conservation of angular momentum and therefore governs the evolution of the warp. For lengthscales that are long compared to the vertical scaleheight of the disc, the known non-resonant and resonant regimes of warp dynamics, including the diffusive and wavelike regimes of Keplerian discs, occur in the local model in the same way as in a global view of warped discs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Standing solitary waves as transitions to spiral structures in gravitationally unstable accretion disks

Astrophysical disks that are sufficiently cold and dense are linearly unstable to the formation of axisymmetric rings as a result of the disk's gravity. In practice, spiral structures are formed, which may in turn produce bound fragments. We study a nonlinear dynamical path that can explain the development of spirals in a local model of a gaseous disk on the subcritical side of the gravitational instability bifurcation. Axisymmetric equilibria can be radially periodic or localized, in the form of standing solitary waves. The solitary solutions have an energy slightly larger than a smooth disk. They are further unstable to non-axisymmetric perturbations with a wide range of azimuthal wavenumbers. The solitary waves may act as a pathway to spirals and fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Importance of magnetic fields in highly eccentric discs with applications to tidal disruption events

Whether tidal disruption events (TDEs) circularise or accrete directly as a highly eccentric disc is the subject of current research and appears to depend sensitively on the disc thermodynamics. In a previous paper we applied the theory of eccentric discs to TDE discs using an $α-$prescription for the disc stress, which leads to solutions that exhibit extreme, potentially unphysical, behaviour. In this paper we further explore the dynamical vertical structure of highly eccentric discs using alternative stress models that are better motivated by the behaviour of magnetic fields in eccentric discs. We find that the presence of a coherent magnetic field has a stabilising effect on the dynamics and can significantly alter the behaviour of highly eccentric radiation dominated discs. We conclude that magnetic fields are important for the evolution of TDE discs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. I. Hydrodynamic simulations

High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of black-hole X-ray binary systems promise insight into strongly curved spacetime. `Discoseismic' oscillations with frequencies set by the intrinsic properties of the central black hole, in particular `trapped inertial waves' (r-modes), offer an attractive explanation for HFQPOs. To produce an observable signature, however, such oscillations must be excited to sufficiently large amplitudes. Turbulence driven by the magnetorotational instability (MRI) fails to provide the necessary amplification, but r-modes may still be excited via interaction with accretion disc warps or eccentricities. We present 3D global hydrodynamic simulations of relativistic accretion discs, which demonstrate for the first time the excitation of trapped inertial waves by an imposed eccentricity in the flow. While the r-modes' saturated state depends on the vertical boundary conditions used in our unstratified, cylindrical framework, their excitation is unambiguous in all runs with eccentricity >0.005 near the ISCO. These simulations provide a proof of concept, demonstrating the robustness of trapped inertial wave excitation in a non-magnetized context. In a companion paper, we explore the competition between this excitation, and damping by magnetohydrodynamic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. II. Magnetohydrodynamic simulations

Trapped inertial oscillations (r-modes) provide a promising explanation for high-frequency quasi-periodic oscillations (HFQPOs) observed in the emission from black hole X-ray binary systems. An eccentricity (or warp) can excite r-modes to large amplitudes, but concurrently the oscillations are likely damped by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). We force eccentricity in global, unstratified, zero-net flux MHD simulations of relativistic accretion discs, and find that a sufficiently strong disc distortion generates trapped inertial waves despite this damping. In our simulations, eccentricities above ~ 0.03 in the inner disc excite trapped waves. In addition to the competition between r-mode damping and driving, we observe that larger amplitude eccentric structures modify and in some cases suppress MRI turbulence. Given the variety of distortions (warps as well as eccentricities) capable of amplifying r-modes, the robustness of trapped inertial wave excitation in the face of MRI turbulence in our simulations provides support for a discoseismic explanation for HFQPOs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wind-MRI interactions in local models of protoplanetary discs: I. Ohmic resistivity

A magnetic disc wind is an important mechanism that may be responsible for driving accretion and structure formation in protoplanetary discs. Recent numerical simulations have shown that these winds can take either the traditional `hourglass' symmetry about the mid-plane, or a `slanted' symmetry dominated by a mid-plane toroidal field of a single sign. The formation of this slanted symmetry state has not previously been explained. We use radially local 1D vertical shearing box simulations to assess the importance of large-scale MRI channel modes in influencing the formation and morphologies of these wind solutions. We consider only Ohmic resistivity and explore the effect of different magnetisations, with the mid-plane $β$ parameter ranging from $10^5$ to $10^2$. We find that our magnetic winds go through three stages of development: cyclic, transitive and steady, with the steady wind taking a slanted symmetry profile similar to those observed in local and global simulations. We show that the cycles are driven by periodic excitation of the $n=2$ or $3$ MRI channel mode coupled with advective eviction, and that the transition to the steady wind is caused by a much more slowly growing $n=1$ mode altering the wind structure. Saturation is achieved through a combination of advective damping from the strong wind, and suppression of the instability due to a strong toroidal field. A higher disc magnetisation leads to a greater tendency towards, and more rapid settling into the slanted symmetry steady wind, which may have important implications for mass and flux transport processes in protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Growth of eccentric modes in disc-planet interactions

We formulate a set of linear equations that describe the behaviour of small eccentricities in a protoplanetary system consisting of a gaseous disc and a planet. Eccentricity propagates through the disc by means of pressure and self-gravity, and is exchanged with the planet via secular interactions. Excitation and damping of eccentricity can occur through Lindblad and corotation resonances, as well as viscosity. We compute normal modes of the coupled disc-planet system in the case of short-period giant planets orbiting inside an inner cavity, possibly carved by the stellar magnetosphere. Three-dimensional effects allow for a mode to be trapped in the inner parts of the disc. This mode can easily grow within the disc's lifetime. An eccentric mode dominated by the planet can also grow, although less rapidly. We compute the structure and growth rates of these modes and their dependence on the assumed properties of the disc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lecture notes: Astrophysical fluid dynamics

These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes, and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is 'frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as determining their stability and their response to tidal forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear hydrodynamical evolution of eccentric Keplerian discs in two dimensions: validation of secular theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear tides in a homogeneous rotating planet or star: global modes and elliptical instability

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}Ω\lesssim -1$, where $n$ and $Ω$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global evolution of the magnetic field in a thin disc and its consequences for protoplanetary systems

The strength and structure of the large-scale magnetic field in protoplanetary discs are still unknown, although they could have important consequences for the dynamics and evolution of the disc. Using a mean-field approach in which we model the effects of turbulence through enhanced diffusion coefficients, we study the time-evolution of the large-scale poloidal magnetic field in a global model of a thin accretion disc, with particular attention to protoplanetary discs. With the transport coefficients usually assumed, the magnetic field strength does not significantly increase radially inwards, leading to a relatively weak magnetic field in the inner part of the disc. We show that with more realistic transport coefficients that take into account the vertical structure of the disc and the back-reaction of the magnetic field on the flow as obtained by Guilet & Ogilvie (2012), the magnetic field can significantly increase radially inwards. The magnetic-field profile adjusts to reach an equilibrium value of the plasma $β$ parameter (the ratio of midplane thermal pressure to magnetic pressure) in the inner part of the disc. This value of $β$ depends strongly on the aspect ratio of the disc and on the turbulent magnetic Prandtl number, and lies in the range $10^4-10^7$ for protoplanetary discs. Such a magnetic field is expected to affect significantly the dynamics of protoplanetary discs by increasing the strength of MHD turbulence and launching an outflow. We discuss the implications of our results for the evolution of protoplanetary discs and for the formation of powerful jets as observed in T-Tauri star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamic instability in eccentric astrophysical discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local and global dynamics of eccentric astrophysical discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal dissipation in stars and giant planets

Astrophysical fluid bodies that orbit close to one another induce tidal distortions and flows that are subject to dissipative processes. The spin and orbital motions undergo a coupled evolution over astronomical timescales, which is relevant for many types of binary star, short-period extrasolar planetary systems and the satellites of the giant planets in the solar system. I review the principal mechanisms that have been discussed for tidal dissipation in stars and giant planets in both linear and nonlinear regimes. I also compare the expectations based on theoretical models with recent observational findings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal interactions of a Maclaurin spheroid. I: Properties of free oscillation modes

We review the work of Bryan (1889) on the normal modes of a Maclaurin spheroid, carrying out numerical calculations of the frequencies and spatial forms of these modes that have not been previously published. We study all modes of degree $l \le 4$, which includes both inertial modes and surface gravity modes, with the aim of better understanding the effect of rapid rotation on tidal interactions. The inclusion of these higher degree modes greatly increases the number of frequencies at which tidal resonances may occur. We derive an expression for the decay rates of these modes to first order in viscosity and explicitly plot these for modes. We see that the equatorial bulge of the spheroid has a significant effect on the decay rates (changing some of these by a factor of 2 between an eccentricity of $e=0$ and $0.5$), and a more modest effect on the mode frequencies. This suggests that models of tidal interaction between rapidly rotating stars and giant planets that model the Coriolis force while neglecting the centrifugal distortion of the body may be in error by an order unity factor. In a subsequent paper we shall examine the case of a forced flow in this spheroid, and complete the model by considering how the tides raised by the orbiting companion change the orbital elements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal interactions of a Maclaurin spheroid. II: Resonant excitation of modes by a close, misaligned orbit

We model a tidally forced star or giant planet as a Maclaurin spheroid, decomposing the motion into the normal modes found by Bryan (1889). We first describe the general prescription for this decomposition and the computation of the tidal power. Although this formalism is very general, forcing due to a companion on a misaligned, circular orbit is used to illustrate the theory. The tidal power is plotted for a variety of orbital radii, misalignment angles, and spheroid rotation rates. Our calculations are carried out including all modes of degree $l \le 4$, and the same degree of gravitational forcing. Remarkably, we find that for close orbits ($a/R_* \approx 3$) and rotational deformations that are typical of giant planets ($e\approx 0.4$) the $l=4$ component of the gravitational potential may significantly enhance the dissipation through resonance with surface gravity modes. There are also a large number of resonances with inertial modes, with the tidal power being locally enhanced by up to three orders of magnitude. For very close orbits ($a/R_* \approx 3$), the contribution to the power from the $l=4$ modes is roughly the same magnitude as that due to the $l = 3$ modes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamic instability in warped astrophysical discs

Warped astrophysical discs are usually treated as laminar viscous flows, which have anomalous properties when the disc is nearly Keplerian and the viscosity is small: fast horizontal shearing motions and large torques are generated, which cause the warp to evolve rapidly, in some cases at a rate that is inversely proportional to the viscosity. However, these flows are often subject to a linear hydrodynamic instability, which may produce small-scale turbulence and modify the large-scale dynamics of the disc. We use a warped shearing sheet to compute the oscillatory laminar flows in a warped disc and to analyse their linear stability by the Floquet method. We find widespread hydrodynamic instability deriving from the parametric resonance of inertial waves. Even very small, unobservable warps in nearly Keplerian discs of low viscosity can be expected to generate hydrodynamic turbulence, or at least wave activity, by this mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local and global dynamics of warped astrophysical discs

Astrophysical discs are warped whenever a misalignment is present in the system, or when a flat disc is made unstable by external forces. The evolution of the shape and mass distribution of a warped disc is driven not only by external influences but also by an internal torque, which transports angular momentum through the disc. This torque depends on internal flows driven by the oscillating pressure gradient associated with the warp, and on physical processes operating on smaller scales, which may include instability and turbulence. We introduce a local model for the detailed study of warped discs. Starting from the shearing sheet of Goldreich & Lynden-Bell, we impose the oscillating geometry of the orbital plane by means of a coordinate transformation. This warped shearing sheet (or box) is suitable for analytical and computational treatments of fluid dynamics, magnetohydrodynamics, etc., and it can be used to compute the internal torque that drives the large-scale evolution of the disc. The simplest hydrodynamic states in the local model are horizontally uniform laminar flows that oscillate at the orbital frequency. These correspond to the nonlinear solutions for warped discs found in previous work by Ogilvie, and we present an alternative derivation and generalization of that theory. In a companion paper we show that these laminar flows are often linearly unstable, especially if the disc is nearly Keplerian and of low viscosity. The local model can be used in future work to determine the nonlinear outcome of the hydrodynamic instability of warped discs, and its interaction with others such as the magnetorotational instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local outflows from turbulent accretion disks

The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Jet launching from accretion discs in the local approximation

The acceleration of an outflow along inclined magnetic field lines emanating from an accretion disc can be studied in the local approximation, as employed in the computational model known as the shearing box. By including the slow magnetosonic point within the computational domain, the rate of mass loss in the outflow can be calculated. The accretion rates of mass and magnetic flux can also be determined, although some effects of cylindrical geometry are omitted. We formulate a simple model for the study of this problem and present the results of one-dimensional numerical simulations and supporting calculations. Quasi-steady solutions are obtained for relatively strong poloidal magnetic fields for which the magnetorotational instability is suppressed. In this regime the rate of mass loss decreases extremely rapidly with increasing field strength, or with decreasing surface density or temperature. If the poloidal magnetic field in an accretion disc can locally achieve an appropriate strength and inclination then a rapid burst of ejection may occur. For weaker fields it may be possible to study the launching process in parallel with the magnetorotational instability, but this will require three-dimensional simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the interaction between tides and convection

We study the interaction between tides and convection in astrophysical bodies by analysing the effect of a homogeneous oscillatory shear on a fluid flow. This model can be taken to represent the interaction between a large-scale periodic tidal deformation and a smaller-scale convective motion. We first consider analytically the limit in which the shear is of low amplitude and the oscillation period is short compared to the timescales of the unperturbed flow. In this limit there is a viscoelastic response and we obtain expressions for the effective elastic modulus and viscosity coefficient. The effective viscosity is inversely proportional to the square of the oscillation frequency, with a coefficient that can be positive, negative or zero depending on the properties of the unperturbed flow. We also carry out direct numerical simulations of Boussinesq convection in an oscillatory shearing box and measure the time-dependent Reynolds stress. The results indicate that the effective viscosity of turbulent convection falls rapidly as the oscillation frequency is increased, attaining small negative values in the cases we have examined, although significant uncertainties remain because of the turbulent noise. We discuss the implications of this analysis for astrophysical tides.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The gravitational instability of a stream of co-orbital particles

We describe the dynamics of a stream of equally spaced macroscopic particles in orbit around a central body (e.g. a planet or star). A co-orbital configuration of small bodies may be subject to gravitational instability, which takes the system to a spreading, disordered and collisional state. We detail the linear instability's mathematical and physical features using the shearing sheet model and subsequently track its nonlinear evolution with local N-body simulations. This model provides a convenient tool with which to understand the gravitational and collisional dynamics of narrow belts, such as Saturn's F-ring and the streams of material wrenched from tidally disrupted bodies. In particular, we study the tendency of these systems to form long-lived particle aggregates. Finally, we uncover an unexpected connection between the linear dynamics of the gravitational instability and the magnetorotational instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal disruption of satellites and formation of narrow rings

In this paper we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using $N$-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellites and differentiated bodies containing a denser core. We show that the Roche limit for a rubble pile is closer to the planet than for a fluid body of the same mean density. The Roche limit for a differentiated body is also closer to the planet than for a homogeneous satellite of the same mean density. Within its Roche limit, a homogeneous satellite totally disrupts and forms a narrow ring. The initial stages of the disruption are similar to the evolution of a viscous fluid ellipsoid, which can be computed semi-analytically. On the other hand, when a differentiated satellite is just within the Roche limit only the mantle is disrupted. This process is similar to Roche-lobe overflow in interacting binary stars and produces two narrow rings on either side of a remnant satellite. We argue that the Uranian rings, and possibly their shepherd satellites, could have been formed through the tidal disruption of a number of protomoons that were formed inside the corotation radius of Uranus and migrated slowly inwards as a result of tidal interaction with the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tides in rotating barotropic fluid bodies: the contribution of inertial waves and the role of internal structure

We discuss the linear response to low-frequency tidal forcing of fluid bodies that are slowly and uniformly rotating, are neutrally stratified and may contain a solid or fluid core. This problem may be regarded as a simplified model of astrophysical tides in convective regions of stars and giant planets. The response can be separated into non-wavelike and wavelike parts, where the former is related instantaneously to the tidal potential and the latter may involve resonances or other singularities. The imaginary part of the potential Love number of the body, which is directly related to the rates of energy and angular momentum exchange in the tidal interaction and to the rate of dissipation of energy, may have a complicated dependence on the tidal frequency. However, a certain frequency-average of this quantity is independent of the dissipative properties of the fluid and can be determined by means of an impulse calculation. The result is a strongly increasing function of the size of the core when the tidal potential is a sectoral harmonic, especially when the body is not strongly centrally condensed. However, the same is not true for tesseral harmonics, which receive a richer response and may therefore be important in determining tidal evolution even though they are usually subdominant in the expansion of the tidal potential. We also discuss analytically the low-frequency response of a slowly rotating homogeneous fluid body to tidal potentials proportional to spherical harmonics of degrees less than five. Tesseral harmonics of degrees greater than two, such as are present in the case of a spin-orbit misalignment, can resonate with inertial modes of the full sphere, leading to an enhanced tidal interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transport of magnetic flux and the vertical structure of accretion discs: I. Uniform diffusion coefficients

The evolution of a large-scale poloidal magnetic field in accretion discs is an important problem because of its role in the launching of jets and winds and in determining the intensity of turbulence. In this paper, we develop a formalism to calculate the transport magnetic flux in a thin accretion disc, thus determining its evolution on a viscous/resistive timescale. The governing equations are derived by performing an asymptotic expansion in the limit of a thin disc, in the regime where the magnetic field is dominated by its vertical component. Turbulent viscosity and resistivity are included, with an arbitrary vertical profile that can be adjusted to mimic the vertical structure of the turbulence. At a given radius and time, the rates of transport of mass and magnetic flux are determined by a one-dimensional problem in the vertical direction, in which the radial gradients of various quantities appear as source terms. We solve this problem to obtain the transport rates and the vertical structure of the disc. The present paper is then restricted to the idealised case of uniform diffusion coefficients, while a companion paper will study more realistic vertical profiles of these coefficients. We show the advection of weak magnetic fields to be significantly faster than the advection of mass, contrary to what a crude vertical averaging might suggest. This results from the larger radial velocities away from the mid-plane, which barely affect the mass accretion owing to the low density in these regions but do affect the advection of magnetic flux. Possible consequences of this larger accretion velocity include a potentially interesting time-dependence with the magnetic flux distribution evolving faster than the mass distribution. If the disc is not too thin, this fast advection may also partially solve the long-standing problem of too efficient diffusion of an inclined magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A model of the entropy flux and Reynolds stress in turbulent convection

We propose a closure model for the transport of entropy and momentum in astrophysical turbulence, intended for application to rotating stellar convective regions. Our closure model is first presented in the Boussinesq formalism, and compared with laboratory and numerical experimental results on Rayleigh-Benard convection and Homogeneous Rayleigh-Benard convection. The predicted angular momentum transport properties of the turbulence in the slowly rotating case recover the well-known Lambda-effect, with an amplitude uniquely related to the convective heat flux. The model is then extended to the anelastic case as well as the fully compressible case. In the special case of spherical symmetry, the predicted radial heat flux is equivalent to that of mixing-length theory. For rotating stars, our model describes the coupled transport of heat and angular momentum, and provides a unified formalism in which to study both differential rotation and thermal inhomogeneities in stellar convection zones.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hydrodynamical simulations of viscous overstability in Saturn's rings

We perform axisymmetric hydrodynamical simulations that describe the nonlinear outcome of the viscous overstability in dense planetary rings. These simulations are particularly relevant for Cassini observations of fine-scale structure in Saturn's A and B-ring, which take the form of periodic microstructure on the 0.1 km scale, and irregular larger-scale structure on 1-10 km. Nonlinear wavetrains dominate all the simulations, and we associate them with the observed periodic microstructure. The waves can undergo small chaotic fluctuations in their phase and amplitude, and may be punctuated by more formidable `wave defects' distributed on longer scales. It is unclear, however, whether the defects are connected to the irregular larger-scale variations observed by Cassini. The long-term behaviour of the simulations is dominated by the imposed boundary conditions, and more generally by the limitations of the local model we use: the shearing box. When periodic boundary conditions are imposed, the system eventually settles on a uniform travelling wave of a predictable wavelength, while reflecting boundaries, and boundaries with buffer zones, maintain a disordered state. The simulations omit self-gravity, though we examine its influence in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The viscous overstability, nonlinear wavetrains, and finescale structure in dense planetary rings

This paper addresses the fine-scale axisymmetric structure exhibited in Saturn's A and B-rings. We aim to explain both the periodic microstructure on 150-220m, revealed by the Cassini UVIS and RSS instruments, and the irregular variations in brightness on 1-10km, reported by the Cassini ISS. We propose that the former structures correspond to the peaks and troughs of the nonlinear wavetrains that form naturally in a viscously overstable disk. The latter variations on longer scales may correspond to modulations and defects in the wavetrains' amplitudes and wavelength. We explore these ideas using a simple hydrodynamical model which captures the correct qualitative behaviour of a disk of inelastically colliding particles, while also permitting us to make progress with analytic and semi-analytic techniques. Specifically, we calculate a family of travelling nonlinear density waves and determine their stability properties. Detailed numerical simulations that confirm our basic results will appear in a following paper.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tidal dissipation in rotating fluid bodies: a simplified model

We study the tidal forcing, propagation and dissipation of linear inertial waves in a rotating fluid body. The intentionally simplified model involves a perfectly rigid core surrounded by a deep ocean consisting of a homogeneous incompressible fluid. Centrifugal effects are neglected, but the Coriolis force is considered in full, and dissipation occurs through viscous or frictional forces. The dissipation rate exhibits a complicated dependence on the tidal frequency and generally increases with the size of the core. In certain intervals of frequency, efficient dissipation is found to occur even for very small values of the coefficient of viscosity or friction. We discuss the results with reference to wave attractors, critical latitudes and other features of the propagation of inertial waves within the fluid, and comment on their relevance for tidal dissipation in planets and stars.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Three-dimensional eccentric discs around Be stars

One-armed oscillation modes in the circumstellar discs of Be stars may explain the cyclical variations in their emission lines. We show that a three-dimensional effect, involving vertical motion and neglected in previous treatments, profoundly influences the dynamics. Using a secular theory of eccentric discs that reduces the problem to a second-order differential equation, we show that confined prograde modes are obtained for all reasonable disc temperatures and stellar rotation rates. We confirm these results using a numerical analysis of the full set of linearized equations for three-dimensional isothermal discs including viscous terms that couple the horizontal motions at different altitudes. In order to make these modes grow, viscous damping must be overcome by an excitation mechanism such as viscous overstability.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Wave attractors and the asymptotic dissipation rate of tidal disturbances

Linear waves in bounded inviscid fluids do not generally form normal modes with regular eigenfunctions. Examples are provided by inertial waves in a rotating fluid contained in a spherical annulus, and internal gravity waves in a stratified fluid contained in a tank with a non-rectangular cross-section. For wave frequencies in the ranges of interest, the inviscid linearized equations are spatially hyperbolic and their characteristic rays are typically focused on to wave attractors. When these systems experience periodic forcing, for example of tidal origin, the response of the fluid can become localized in the neighbourhood of a wave attractor. In this paper I define a prototypical problem of this form and construct analytically the long-term response to a periodic body force in the asymptotic limit of small viscosity. The vorticity of the fluid is localized in a detached shear layer close to the wave attractor in such a way that the total rate of dissipation of energy is asymptotically independent of the viscosity. I further demonstrate that the same asymptotic dissipation rate is obtained if a non-viscous damping force is substituted for the Navier-Stokes viscosity. I discuss the application of these results to the problem of tidal forcing in giant planets and stars, where the excitation and dissipation of inertial waves may make a dominant, or at least important, contribution to the orbital and spin evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $α_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $α_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $α_{\rm v}$ is approximately equal to $α_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $α_{\rm v}$ and $α_{\rm h}$ we conclude that for $β= p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

On the difficulty of launching an outflow from an accretion disk

We solve for the local vertical structure of a thin accretion disk threaded by a poloidal magnetic field. The angular velocity deviates from the Keplerian value as a result of the radial Lorentz force, but is constant on magnetic surfaces. Angular momentum transport and energy dissipation in the disk are parametrized by an alpha-prescription, and a Kramers opacity law is assumed to hold. We also determine the stability of the equilibria with respect to the magnetorotational (or Balbus-Hawley) instability. If the magnetic field is sufficiently strong, stable equilibria can be found in which the angle of inclination, i, of the magnetic field to the vertical at the surface of the disk has any value in the range [0,90 degrees). By analyzing the dynamics of a transonic outflow in the corona of the disk, we show that a certain potential difference must be overcome even when i > 30 degrees. We determine this potential difference as a function of i for increasing values of the vertical magnetic field strength. For magnetorotationally stable equilibria, the potential difference increases faster than the fourth power of the magnetic field strength, quickly exceeding a value corresponding to the central temperature of the disk, and is minimized with respect to i at i approximately equal to 38 degrees. We show that this property is relatively insensitive to the form of the opacity law. Our results suggest that an additional source of energy, such as coronal heating, may be required for the launching of an outflow from a magnetized disk.

0 citations0 reviewsNo code linkedOpen access