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.

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.

The diurnal tidal forces can excite a normal mode of the Earth's core, the free inner core nutation (FICN), which is characterized by a tilt of the rotation axis of the inner core with respect to the rotation axis of the outer core. The differential rotation between the inner core and the outer core induces fluid motions in the outer core and gives rise to Ohmic dissipation in the presence of the Earth's internal magnetic field. Nutation measurements can reflect such dissipation if it is sufficiently strong and thus can provide insights into the properties and dynamics of the Earth's core. In this study we perform a set of numerical calculations of the linear perturbations in the outer core induced by the FICN at very low Ekman numbers (as small as $10^{-11}$). Our numerical results show that the back-reaction of the magnetic field notably alters the structure and length scale of the perturbations induced by the FICN, and thus influences the Ohmic dissipation resulting from the perturbations. When the Ekman number is sufficiently small, Ohmic dissipation tends to be insensitive to the fluid viscosity and to the magnetic diffusivity, which allows us to estimate the Ohmic dissipation associated with the FICN without relying on an extrapolation. In contrast to the results of Buffett (2010b), the estimated Ohmic dissipation based on our numerical calculations is too weak to account for the observed damping of the FICN mode. This also implies that nutation measurements cannot provide effective constraints on the strength of the magnetic field inside the Earth's outer core.

We develop a fully nonlinear approximation to the short-wavelength limit of eccentric waves in astrophysical discs, based on the averaged Lagrangian method of Whitham (1965). In this limit there is a separation of scales between the rapidly varying eccentric wave and the background disc. Despite having small eccentricities, such rapidly varying waves can be highly nonlinear, potentially approaching orbital intersection, and this can result in strong pressure gradients in the disc. We derive conditions for the steepening of nonlinearity and eccentricity as the waves propagate in a radially structured disc in this short-wavelength limit and show that the behaviour of the solution can be bounded by the behaviour of the WKB solution to the linearised equations.

The evolution of a large-scale poloidal magnetic field in an accretion disc is an important problem because it determines the launching of winds and the feasibility of the magnetorotational instability to generate turbulence or channel flows. Recent studies, both semi-analytical calculations and numerical simulations, have highlighted the crucial role non-ideal MHD effects (Ohmic resistivity, Hall drift and ambipolar diffusion), relevant in the protoplanetary disc context, might play in magnetic flux evolution in the disc. We investigate the flux transport in discs through the use of two one-dimensional semi-analytic models in the vertical direction, exploring regimes where different physical source terms and effects dominate. The governing equations for both models are derived by performing an asymptotic expansion in the limit of a thin disc, with the different regimes isolated through setting the relative order of the leading terms between variables. Flux transport rates and vertical structure profiles are calculated for a range of diffusivities and disc magnetisations. We found that Ohmic and ambipolar diffusivities drive radially outward flux transport with an outwardly inclined field. A wind outflow drives inward flux transport, which is significantly enhanced in the presence of Hall drift in the positive polarity case, $\eta_H (\boldsymbol{B}_z \cdot \boldsymbol{\Omega}) > 0$, an effect which has only been briefly noted before. Coupled only with outward inclination, the Hall effect reduces the flux transport given by a background Ohmic and/or ambipolar diffusivity, but drives no flux transport when it is the only non-ideal effect present.

We show that the ideal hydrodynamics of an eccentric astrophysical disc can be derived from a variational principle. The nonlinear secular theory describes the slow evolution of a continuous set of nested elliptical orbits as a result of the pressure in a thin disc. In the artificial but widely considered case of a 2D disc, the hydrodynamic Hamiltonian is just the orbit-averaged internal energy of the disc, which can be determined from its eccentricity distribution using the geometry of the elliptical orbits. In the realistic case of a 3D disc, the Hamiltonian needs to be modified to take into account the dynamical vertical structure of the disc. The simplest solutions of the theory are uniformly precessing nonlinear eccentric modes, which make the energy stationary subject to the angular momentum being fixed. We present numerical examples of nonlinear eccentric modes up to their limiting amplitudes. Although it lacks dissipation, which is important in many astrophysical contexts, this formalism allows a simpler theoretical approach to the nonlinear dynamics of eccentric discs than that derived from stress integrals, and also connects better with established methods of celestial mechanics for cases in which the disc interacts gravitationally with one or more orbital companion.

We study the hydrodynamical stability of the laminar flows associated with warped astrophysical discs using numerical simulations of warped shearing boxes. We recover linear growth rates reported previously due to a parametric resonance of inertial waves, and show that the nonlinear saturated state can significantly reduce the laminar flows, meaning that the warp would evolve on much longer time scales than would be concluded from the internal torques due to these laminar flows. Towards larger warp amplitudes, we find first of all a reversal of angular momentum flux, indicating that the mass distribution would evolve in an anti-diffusive manner, and second that the linear growth rates disappear, possibly because of the very strong shear in the laminar flows in this regime. For discs with small enough viscosity, a nonlinear state can still be found when linear growth rates are absent by introducing a large enough perturbation, either by starting from a nonlinear state obtained at smaller warp amplitude, or by starting from a state with no laminar flows.

The excitation of trapped inertial waves (r-modes) by warps and eccentricities in the inner regions of a black hole accretion disc may explain the high-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of Galactic X-ray binaries. However, it has been suggested that strong vertical magnetic fields push the oscillations' trapping region toward the innermost stable circular orbit (ISCO), where conditions could be unfavourable for their excitation. This paper explores the effects of large-scale magnetic fields that exhibit \textit{both} toroidal and vertical components, through local and global linear analyses. We find that a strong toroidal magnetic field can reduce the detrimental effects of a vertical field: in fact, the isolation of the trapping region from the ISCO may be restored by toroidal magnetic fields approaching thermal strengths. The toroidal field couples the r-modes to the disc's magneto-acoustic response and inflates the effective pressure within the oscillations. As a consequence, the restoring force associated with the vertical magnetic field's tension is reduced. Given the analytical and numerical evidence that accretion discs threaded by poloidal magnetic field lines develop a strong toroidal component, our result provides further evidence that the detrimental effects of magnetic fields on trapped inertial modes are not as great as previously thought.

Non-linear evolution of the parametric instability of inertial waves inherent to eccentric discs is studied by way of a new local numerical model. Mode coupling of tidal deformation with the disc eccentricity is known to produce exponentially growing eccentricities at certain mean-motion resonances. However, the details of an efficient saturation mechanism balancing this growth still are not fully understood. This paper develops a local numerical model for an eccentric quasi-axisymmetric shearing box which generalises the often-used cartesian shearing box model. The numerical method is an overall second order well-balanced finite volume method which maintains the stratified and oscillatory steady-state solution by construction. This implementation is employed to study the non-linear outcome of the parametric instability in eccentric discs with vertical structure. Stratification is found to constrain the perturbation energy near the mid-plane and localise the effective region of inertial wave breaking that sources turbulence. A saturated marginally sonic turbulent state results from the non-linear breaking of inertial waves and is subsequently unstable to large-scale axisymmetric zonal flow structures. This resulting limit-cycle behaviour reduces access to the eccentric energy source and prevents substantial transport of angular momentum radially through the disc. Still, the saturation of this parametric instability of inertial waves is shown to damp eccentricity on a time-scale of a thousand orbital periods. It may thus be a promising mechanism for intermittently regaining balance with the exponential growth of eccentricity from the eccentric Lindblad resonances and may also help explain the occurrence of "bursty" dynamics such as the superhump phenomenon.

We derive expressions for the local ideal magnetohydrodynamic (MHD) equations for a warped astrophysical disc using a warped shearing box formalism. A perturbation expansion of these equations to first order in the warping amplitude leads to a linear theory for the internal local structure of magnetised warped discs in the absence of MRI turbulence. In the special case of an external magnetic field oriented normal to the disc surface these equations are solved semi-analytically via a spectral method. The relatively rapid warp propagation of low-viscosity Keplerian hydrodynamic warped discs is diminished by the presence of a magnetic field. The magnetic tension adds a stiffness to the epicyclic oscillations, detuning the natural frequency from the orbital frequency and thereby removing the resonant forcing of epicyclic modes characteristic of hydrodynamic warped discs. In contrast to a single hydrodynamic resonance we find a series of Alfv\'{e}nic-epicyclic modes which may be resonantly forced by the warped geometry at critical values of the orbital shear rate $q$ and magnetic field strength. At these critical points large internal torques are generated and anomalously rapid warp propagation occurs. As our treatment omits MRI turbulence, these results are of greatest applicability to strongly magnetised discs.

Thin astrophysical discs are very often modelled using the equations of two-dimensional hydrodynamics. We derive an extension of this model that describes more accurately the behaviour of a thin disc in the absence of self-gravity, magnetic fields and complex internal motions. The ideal fluid theory is derived directly from Hamilton's Principle for a three-dimensional fluid after making a specific approximation to the deformation gradient tensor. We express the equations in Eulerian form after projection on to a reference plane. The disc is thought of as a set of fluid columns, each of which is capable of a time-dependent affine transformation, consisting of a translation together with a linear transformation in three dimensions. Therefore, in addition to the usual two-dimensional hydrodynamics in the reference plane, the theory allows for a deformation of the midplane (as occurs in warped discs) and for the internal shearing motions that accompany such deformations. It also allows for the vertical expansions driven in non-circular discs by a variation of the vertical gravitational field around the horizontal streamlines, or by a divergence of the horizontal velocity. The equations of the affine model embody conservation laws for energy and potential vorticity, even for non-planar discs. We verify that they reproduce exactly the linear theories of three-dimensional warped and eccentric discs in a secular approximation. However, the affine model does not rely on any secular or small-amplitude assumptions and should be useful in more general circumstances.

The high-frequency quasi-periodic oscillations (HFQPOs) that punctuate the light curves of X-ray binary systems present a window onto the intrinsic properties of stellar-mass black holes and hence a testbed for general relativity. One explanation for these features is that relativistic distortion of the accretion disc's differential rotation creates a trapping region in which inertial waves (r-modes) might grow to observable amplitudes. Local analyses, however, predict that large-scale magnetic fields push this trapping region to the inner disc edge, where conditions may be unfavorable for r-mode growth. We revisit this problem from a pseudo-Newtonian but fully global perspective, deriving linearized equations describing a relativistic, magnetized accretion flow, and calculating normal modes with and without vertical density stratification. In an unstratified model, the choice of vertical wavenumber determines the extent to which vertical magnetic fields drive the r-modes toward the inner edge. In a global model fully incorporating density stratification, we confirm that this susceptibility to magnetic fields depends on disc thickness. Our calculations suggest that in thin discs, r-modes may remain independent of the inner disc edge for vertical magnetic fields with plasma betas as low as $\beta\approx 100-300$. We posit that the appearance of r-modes in observations may be more determined by a competition between excitation and damping mechanisms near the ISCO than the modification of the trapping region by magnetic fields.

We investigate effects of the presence of a magnetic field on tidal dissipation in rotating fluid bodies. We consider a simplified model consisting of a rigid core and a fluid envelope, permeated by a background magnetic field (either a dipolar field or a uniform axial field). The wavelike tidal responses in the fluid layer are in the form of magnetic-Coriolis waves, which are restored by both the Coriolis force and the Lorentz force. Energy dissipation occurs through viscous damping and Ohmic damping of these waves. Our numerical results show that the tidal dissipation can be dominated by Ohmic damping even with a weak magnetic field. The presence of a magnetic field smooths out the complicated frequency-dependence of the dissipation rate, and broadens the frequency spectrum of the dissipation rate, depending on the strength of the background magnetic field. However, the frequency-averaged dissipation is independent of the strength and structure of the magnetic field, and of the dissipative parameters, in the approximation that the wave-like response is driven only by the Coriolis force acting on the non-wavelike tidal flow. Indeed, the frequency-averaged dissipation quantity is in good agreement with previous analytical results in the absence of magnetic fields. Our results suggest that the frequency-averaged tidal dissipation of the wavelike perturbations is insensitive to detailed damping mechanisms and dissipative properties.

We conduct a linear stability calculation of an ideal Keplerian flow on which a sinusoidal zonal flow is imposed. The analysis uses the shearing sheet model and is carried out both in isothermal and adiabatic conditions, with and without self-gravity (SG). In the non-SG regime a structure in the potential vorticity (PV) leads to a non-axisymmetric Kelvin-Helmholtz (KH) instability; in the short-wavelength limit its growth rate agrees with the incompressible calculation by Lithwick (2007), which only considers perturbations elongated in the streamwise direction. The instability's strength is analysed as a function of the structure's properties, and zonal flows are found to be stable if their wavelength is $\gtrsim 8H$, where $H$ is the disc's scale height, regardless of the value of the adiabatic index $\gamma$. The non-axisymmetric KH instability can operate in Rayleigh-stable conditions, and it therefore represents the limiting factor to the structure's properties. Introducing SG triggers a second non-axisymmetric instability, which is found to be located around a PV maximum, while the KH instability is linked to a PV minimum, as expected. In the adiabatic regime, the same gravitational instability is detected even when the structure is present only in the entropy (not in the PV) and the instability spreads to weaker SG conditions as the entropy structure's amplitude is increased. This eventually yields a non-axisymmetric instability in the non-SG regime, albeit of weak strength, localised around an entropy maximum.

Recent results by Martin et al. (2014) showed in 3D SPH simulations that tilted discs in binary systems can be unstable to the development of global, damped Kozai-Lidov (KL) oscillations in which the discs exchange tilt for eccentricity. We investigate the linear stability of KL modes for tilted inviscid discs under the approximations that the disc eccentricity is small and the disc remains flat. By using 1D equations, we are able to probe regimes of large ratios of outer to inner disc edge radii that are realistic for binary systems of hundreds of AU separations and are not easily probed by multi-dimensional simulations. For order unity binary mass ratios, KL instability is possible for a window of disc aspect ratios H/r in the outer parts of a disc that roughly scale as (n_b/n)^2 < H/r < n_b/n, for binary orbital frequency n_b and orbital frequency n at the disc outer edge. We present a framework for understanding the zones of instability based on the determination of branches of marginally unstable modes. In general, multiple growing eccentric KL modes can be present in a disc. Coplanar apsidal-nodal precession resonances delineate instability branches. We determine the range of tilt angles for unstable modes as a function of disc aspect ratio. Unlike the KL instability for free particles that involves a critical (minimum) tilt angle, disc instability is possible for any nonzero tilt angle depending on the disc aspect ratio.

We study the interaction between stellar irradiation and tidal heating in gaseous planets with short orbital periods. The intentionally simplified atmospheric model we employ makes the problem analytically tractable and permits the derivation of useful scaling relations. We show that many tidal models provide thermal feedback, producing interior radiative zones and leading to enhanced g-mode dissipation with a wide spectrum of resonances. These resonances are dynamically tuned by the thermal feedback, and so represent a novel form of thermomechanical feedback, coupling vibrational modes to the very slow thermal evolution of the planet. We then show that stellar irradiation allows the heat produced by these modes to be trapped at depth with high efficiency, leading to entropy increase in the central convective region, as well as expansion of the planet's radius sufficient to match observed swelling. We find that thermally driven winds play an essential role in this process by making the thermal structure of the atmosphere spherically symmetric within a few scale heights of the photosphere. We characterise the relationship between the swelling factor, the orbital period and the host star and determine the timescale for swelling. We show that these g-modes suffice to produce bloating on the order of the radius of the planet over $\mathrm{Gyr}$ timescales when combined with significant insolation and we provide analytic relations for the relative magnitudes of tidal heating and insolation.

Spin-orbit misalignments have been detected in exoplanetary systems and binary star systems. Tidal interactions may have played an important role in the evolution of the spin-orbit angle. In this study, we investigate the tidal interactions in spin-orbit misaligned systems. In particular, we focus on the tidal response of a rotating fluid body to the obliquity tide, which may be important for the evolution of the spin-orbit angle but hardly affects the orbital evolution. The obliquity tide also provides a torque for the mutual precession of the spin and orbital axes around the total angular momentum vector, which has not yet been considered in previous studies on the tidal interactions. In this paper, we first formulate a set of linearized equations describing the tidal response in spin-orbit misaligned systems, taking into account the precessional motion. Numerical solutions in a homogeneous fluid and in a polytrope of index 1 show that dissipative inertial waves can be excited on top of precession by the obliquity tide in the presence of a rigid core. The tidal quality factor associated with the obliquity tide $Q'_{210}$ can be several orders of magnitude smaller than those associated with other tidal components if their frequencies fall outside the frequency range of inertial waves. Therefore, it is possible that the spin-orbit misalignment undergoes much more rapid decay than the orbital decay in hot Jupiter systems owing to the enhanced dissipation of the obliquity tide.

We explore the evolution of the eccentricity of an accretion disc perturbed by an embedded planet whose mass is sufficient to open a large gap in the disc. Various methods for representing the orbit-averaged motion of an eccentric disc are discussed. We characterize the linear instability which leads to the growth of eccentricity by means of hydrodynamical simulations. We numerically recover the known result that eccentricity growth in the disc is possible when the planet-to-star mass ratio exceeds 0.003. For mass ratios larger than this threshold, the precession rates and growth rates derived from simulations, as well as the shape of the eccentric mode, compare well with the predictions of a linear theory of eccentric discs. We study mechanisms by which the eccentricity growth eventually saturates into a non-linear regime.

This chapter explores the physics shared by planetary rings and the various disks that populate the Universe. It begins with an observational overview, ranging from protoplanetary disks to spiral galaxies, and then compares and contrasts these astrophysical disks with the rings of the Solar System. Emphasis is placed on fundamental physics and dynamics, and how research into the two classes of object connects. Topics covered include disk formation, accretion, collisional processes, waves, instabilities, and satellite-disk interactions.

This work presents a linear analytical calculation on the stability and evolution of a compressible, viscous self-gravitating (SG) Keplerian disc with both horizontal thermal diffusion and a constant cooling timescale when an axisymmetric structure is present and freely evolving. The calculation makes use of the shearing sheet model and is carried out for a range of cooling times. Although the solutions to the inviscid problem with no cooling or diffusion are well known, it is non-trivial to predict the effect caused by the introduction of cooling and of small diffusivities; this work focuses on perturbations of intermediate wavelengths, therefore representing an extension to the classical stability analysis on thermal and viscous instabilities. For density wave modes the analysis can be simplified by means of a regular perturbation analysis; considering both shear and thermal diffusivities, the system is found to be overstable for intermediate and long wavelengths for values of the Toomre parameter $Q \lesssim 2$; a non-SG instability is also detected for wavelengths $\gtrsim 18H$, where $H$ is the disc scale height, as long as $\gamma \lesssim 1.305$. The regular perturbation analysis does not however hold for the entropy and potential vorticity slow modes as their ideal growth rates are degenerate. To understand their evolution, equations for the axisymmetric structure's amplitudes in these two quantities are analytically derived and their instability regions obtained. The instability appears boosted by increasing the value of the adiabatic index and of the Prandtl number, while it is quenched by efficient cooling.

Jets, from the protostellar to the AGN context, have been extensively studied but their connection to the turbulent dynamics of the underlying accretion disc is poorly understood. Following a similar approach to Lesur et al. (2013), we examine the role of the magnetorotational instability (MRI) in the production and acceleration of outflows from discs. Via a suite of one-dimensional shearing-box simulations of stratified discs we show that magneto-centrifugal winds exhibit cyclic activity with a period of 10-20 Omega^{-1}, a few times the orbital period. The cycle seems to be more vigorous for strong vertical field; it is robust to the variation of relevant parameters and independent of numerical details. The convergence of these solutions (in particular the mass loss rate) with vertical box size is also studied. By considering a sequence of magnetohydrostatic equilibria and their stability, the periodic activity may be understood as the succession of the following phases: (a) a dominant MRI channel mode, (b) strong magnetic field generation, (c) consequent wind launching, and ultimately (d) vertical expulsion of the excess magnetic field by the expanding and accelerating gas associated with the wind. We discuss potential connections between this behaviour and observed time-variability in disk-jet systems.

In Keplerian accretion disks, turbulence and magnetic fields may be jointly excited through a subcritical dynamo process involving the magnetorotational instability (MRI). High-resolution simulations exhibit a tendency towards statistical self-organization of MRI dynamo turbulence into large-scale cyclic dynamics. Understanding the physical origin of these structures, and whether they can be sustained and transport angular momentum efficiently in astrophysical conditions, represents a significant theoretical challenge. The discovery of simple periodic nonlinear MRI dynamo solutions has recently proven useful in this respect, and has notably served to highlight the role of turbulent magnetic diffusion in the seeming decay of the dynamics at low magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a common regime in accretion disks. The connection between these simple structures and the statistical organization reported in turbulent simulations remained elusive, though. Here, we report the numerical discovery in moderate aspect ratio Keplerian shearing boxes of new periodic, incompressible, three-dimensional nonlinear MRI dynamo solutions with a larger dynamical complexity reminiscent of such simulations. These "chimera" cycles are characterized by multiple MRI-unstable dynamical stages, but their basic physical principles of self-sustainment are nevertheless identical to those of simpler cycles found in azimuthally elongated boxes. In particular, we find that they are not sustained at low Pm either due to subcritical turbulent magnetic diffusion. These solutions offer a new perspective into the transition from laminar to turbulent instability-driven dynamos, and may prove useful to devise improved statistical models of turbulent accretion disk dynamos.

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.

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 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.

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.

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

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.

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.

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.

The ejecta discharged by impacting meteorites can redistribute a planetary ring's mass and angular momentum. This `ballistic transport' of ring properties instigates a linear instability that could generate the 100--1000-km undulations observed in Saturn's inner B-ring and in its C-ring. We present semi-analytic results demonstrating how the instability sustains steadily travelling nonlinear wavetrains. At low optical depths, the instability produces approximately sinusoidal waves of low amplitude, which we identify with those observed between radii 77,000 and 86,000 km in the C-ring. On the other hand, optical depths of 1 or more exhibit hysteresis, whereby the ring falls into multiple stable states: the homogeneous background equilibrium or large-amplitude wave states. Possibly the `flat zones' and `wave zones' between radii 93,000 and 98,000 km in the B-ring correspond to the stable homogeneous and wave states, respectively. In addition, we test the linear stability of the wavetrains and show that only a small subset are stable. In particular, stable solutions all possess wavelengths greater than the lengthscale of fastest linear growth. We supplement our calculations with a weakly nonlinear analysis that suggests the C-ring reproduces some of the dynamics of the complex Ginzburg--Landau equation. In the third paper in the series, these results will be tested and extended with numerical simulations.

Saturn's inner B-ring and its C-ring support wavetrains of contrasting amplitudes but with similar length scales, 100-1000 km. In addition, the inner B-ring is punctuated by two intriguing `flat' regions between radii 93,000 km and 98,000 km in which the waves die out, whereas the C-ring waves coexist with a forest of plateaus, narrow ringlets, and gaps. In both regions the waves are probably generated by a large-scale linear instability whose origin lies in the meteoritic bombardment of the rings: the ballistic transport instability. In this paper, the third in a series, we numerically simulate the long-term nonlinear evolution of this instability in a convenient local model. Our C-ring simulations confirm that the unstable system forms low-amplitude wavetrains possessing a preferred band of wavelengths. B-ring simulations, on the other hand, exhibit localised nonlinear wave `packets' separated by linearly stable flat zones. Wave packets travel slowly while spreading in time, a result that suggests the observed flat regions in Saturn's B-ring are shrinking. Finally, we present exploratory runs of the inner B-ring edge which reproduce earlier numerical results: ballistic transport can maintain the sharpness of a spreading edge while building a `ramp' structure at its base. Moreover, the ballistic transport instability can afflict the ramp region, but only in low-viscosity runs.

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.

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 $\beta$ parameter (the ratio of midplane thermal pressure to magnetic pressure) in the inner part of the disc. This value of $\beta$ 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.

In a uniformly rotating fluid, inertial waves propagate along rays that are inclined to the rotation axis by an angle that depends on the wave frequency. In closed domains, multiple reflections from the boundaries may cause inertial waves to focus on to particular structures known as wave attractors. Such structures have previously been studied from a theoretical point of view, in laboratory experiments, in linear numerical calculations and in some recent numerical simulations. In the present paper, two-dimensional direct numerical simulations of an inertial wave attractor are presented. In the linear regime, we first recover the results of the linear calculations and asymptotic theory of Ogilvie (2005) who considered a prototypical problem involving the focusing of linear internal waves into a narrow beam centred on a wave attractor in a steady state. The velocity profile of the beam and its scalings with the Ekman number, as well as the asymptotic value of the dissipation rate, are found to be in agreement with the linear theory. We also find that, as the beam builds up around the wave attractor, the power input by the applied force reaches its limiting value more rapidly than the dissipation rate, which saturates only when the beam has reached its final thickness. In the nonlinear regime, the beam is strongly affected and becomes unstable to a subharmonic instability. This instability transfers energy to secondary waves possessing shorter wavelengths and lower frequencies. The onset of the instability of a narrow inertial wave beam is investigated by means of a separate linear analysis and the results, such as the onset of the instability, are found to be consistent with the global simulations of the wave attractor. The excitation of such secondary waves described theoretically in this work has also been seen in recent laboratory experiments on internal gravity waves.

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.

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.

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.

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, nonlinear magnetohydrodynamic process whose study is relevant to the understanding of accretion and magnetic field generation in astrophysics. Transition to this form of dynamo is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. This suggests that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in various astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows.

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.

We investigate the radial transport of magnetic flux in a thin accretion disc, the turbulence being modelled by effective diffusion coefficients (viscosity and resistivity). Both turbulent diffusion and advection by the accretion flow contribute to flux transport, and they are likely to act in opposition. We study the consequences of the vertical variation of the diffusion coefficients, due to a varying strength of the turbulence. For this purpose, we consider three different vertical profiles of these coefficients. The first one is aimed at mimicking the turbulent stress profile observed in numerical simulations of MHD turbulence in stratified discs. This enables us to confirm the robustness of the main result of Paper I obtained for uniform diffusion coefficients that, for weak magnetic fields, the contribution of the accretion flow to the transport velocity of magnetic flux is much larger than the transport velocity of mass. We then consider the presence of a dead zone around the equatorial plane, where the physical resistivity is high while the turbulent viscosity is low. We find that it amplifies the previous effect: weak magnetic fields can be advected orders of magnitude faster than mass, for dead zones with a large vertical extension. The ratio of advection to diffusion, determining the maximum inclination of the field at the surface of the disc, is however not much affected. Finally, we study the effect of a non-turbulent layer at the surface of the disc, which has been suggested as a way to reduce the diffusion of the magnetic flux. We find that the reduction of the diffusion requires the conducting layer to extend below the height at which the magnetic pressure equals the thermal pressure. As a consequence, if the absence of turbulence is caused by the large-scale magnetic field, the highly conducting layer is inefficient at reducing the diffusion.

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.

Magnetorotational turbulence and magnetically driven disc winds are often considered as separate processes. However, realistic astrophysical discs are expected to be subject to both effects, although possibly at different times and locations. We investigate here the potential link between these two phenomena using a mixed numerical and analytical approach. We show in particular that large-scale MRI modes which dominate strongly magnetised discs (plasma beta~10) naturally produce magnetically driven outflows in the nonlinear regime. We show that these outflows share many similarities with local and global disc wind solutions found in the literature. We also investigate the 3D stability of these outflows and show that they are unstable on dynamical timescales. The implications of these results for the transition between a jet-emitting disc and a standard "viscous" disc are discussed.

Planetary rings sustain a continual bombardment of hypervelocity meteoroids that erode the surfaces of ring particles on time scales of 10^5 - 10^7 years. The debris ejected from such impacts re-accretes on to the ring, though often at a slightly different orbital radius from the point of emission. This `ballistic transport' leads to a rearrangement of the disk's mass and angular momentum, and gives rise to a linear instability that generates structure on relatively large scales. It is likely that the 100-km undulations in Saturn's inner B-ring and the plateaus and 1000-km waves in Saturn's C-ring are connected to the nonlinear saturation of the instability. In this paper the physical problem is reformulated so as to apply to a local patch of disk (the shearing sheet). This new streamlined model helps facilitate our physical understanding of the instability, and also makes more tractable the analysis of its nonlinear dynamics. We concentrate on the linear theory in this paper, showing that the instability is restricted to a preferred range of intermediate wavenumbers and optical depths. We subsequently apply these general results to the inner B-ring and the C-ring and find that in both regions the ballistic transport instability should be near marginality, a fact that may have important consequences for its prevalence and nonlinear development. Owing to damping via self-gravity wakes, the instability should not be present in the A-ring. A following paper will explore the instability's nonlinear saturation and how it connects to the observed large-scale structure.

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.

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.

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.

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.

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.

The nature of dynamo action in shear flows prone to magnetohydrodynamic instabilities is investigated using the magnetorotational dynamo in Keplerian shear flow as a prototype problem. Using direct numerical simulations and Newton's method, we compute an exact time-periodic magnetorotational dynamo solution to the three-dimensional dissipative incompressible magnetohydrodynamic equations with rotation and shear. We discuss the physical mechanism behind the cycle and show that it results from a combination of linear and nonlinear interactions between a large-scale axisymmetric toroidal magnetic field and non-axisymmetric perturbations amplified by the magnetorotational instability. We demonstrate that this large scale dynamo mechanism is overall intrinsically nonlinear and not reducible to the standard mean-field dynamo formalism. Our results therefore provide clear evidence for a generic nonlinear generation mechanism of time-dependent coherent large-scale magnetic fields in shear flows and call for new theoretical dynamo models. These findings may offer important clues to understand the transitional and statistical properties of subcritical magnetorotational turbulence.

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 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.

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.

The mechanism of angular momentum transport in accretion discs has long been debated. Although the magnetorotational instability appears to be a promising process, poorly ionized regions of accretion discs may not undergo this instability. In this letter, we revisit the possibility of transporting angular momentum by turbulent thermal convection. Using high-resolution spectral methods, we show that strongly turbulent convection can drive outward angular momentum transport at a rate that is, under certain conditions, compatible with observations of discs. We find however that the angular momentum transport is always much weaker than the vertical heat transport. These results indicate that convection might be another way to explain global disc evolution, provided that a sufficiently unstable vertical temperature profile can be maintained.

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.

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.

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.

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.

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.

We perform a linear analysis to investigate the dynamical response of a non-synchronized hot Jupiter to stellar irradiation. In this work, we consider the diurnal Fourier harmonic of the stellar irradiation acting at the top of a radiative layer of a hot Jupiter with no clouds and winds. In the absence of the Coriolis force, the diurnal thermal forcing can excite internal waves propagating into the planet's interior when the thermal forcing period is longer than the sound crossing time of the planet's surface. When the Coriolis effect is taken into consideration, the latitude-dependent stellar heating can excite weak internal waves (g modes) and/or strong baroclinic Rossby waves (buoyant r modes) depending on the asynchrony of the planet. When the planet spins faster than its orbital motion (i.e. retrograde thermal forcing), these waves carry negative angular momentum and are damped by radiative loss as they propagate downwards from the upper layer of the radiative zone. As a result, angular momentum is transferred from the lower layer of the radiative zone to the upper layer and generates a vertical shear. We estimate the resulting internal torques for different rotation periods based on the parameters of HD 209458b.

We consider the inward propagation of warping and eccentric disturbances in discs around black holes under a wide variety of conditions. In our calculations we use secular theories of warped and eccentric discs and assume the deformations to be stationary and propagating in a disc model similar to regions (a) and (b) of Shakura and Sunyaev discs. We find that the propagation of deformations to the innermost regions of the disc is facilitated for low viscous damping and high accretion rate. We relate our results to the possible excitation of trapped inertial modes, and to the observations of high-frequency quasi-periodic oscillations (QPOs) in black hole systems in the very high spectral state.

Accretion discs around black holes in X-ray binary stars are warped if the spin axis of the black hole is not perpendicular to the binary orbital plane. They can also become eccentric through an instability involving a resonance with the binary orbit. Depending on the thickness of the disc and the efficiency of dissipative processes, these global deformations may be able to propagate into the innermost part of the disc in the form of stationary bending or density waves. We describe the solutions in the linear regime and discuss the conditions under which a warp or eccentricity is likely to produce significant activity in the inner region, which may include the excitation of quasi-periodic oscillations.

High-frequency quasi-periodic oscillations detected in the light curves of black hole candidates can, according to one model, be identified with hydrodynamic oscillations of the accretion disc. We describe a non-linear coupling mechanism, suggested by Kato, through which inertial waves trapped in the inner regions of accretion discs around black holes are excited. Global warping and/or eccentricity of the disc have a fundamental role in this coupling: they combine with trapped modes, generating negative energy waves, that are damped as they approach the inner edge of the disc or their corotation resonance. As a result of this damping, inertial oscillations are amplified. We calculate the resulting eigenfunctions and their growth rates.

(Abriged) The magnetorotational instability (MRI) is believed to be an efficient way to transport angular momentum in accretion discs. It has also been suggested as a way to amplify magnetic fields in discs, the instability acting as a nonlinear dynamo. Recent numerical work has shown that a large-scale magnetic field, which is predominantly azimuthal, can be sustained by motions driven by the MRI of this same field. Following this idea, we present an analytical calculation of the MRI in the presence of an azimuthal field with a non-trivial vertical structure. We find that the mean radial EMF associated to MRI modes tends to reduce the magnetic energy, acting like a turbulent resistivity by mixing the non-uniform azimuthal field. Meanwhile, the azimuthal EMF generates a radial field that, in combination with the Keplerian shear, tends to amplify the azimuthal field and can therefore assist in the dynamo process. This effect, however, is reversed for sufficiently strong azimuthal fields, naturally leading to a saturation of the dynamo and possibly to a cyclic behaviour of the magnetic field, as found in previous numerical works.

(abridged) MHD turbulence is known to exist in shearing boxes with either zero or nonzero net magnetic flux. However, the way turbulence survives in the zero-net-flux case is not explained by linear theory and appears as a purely numerical result. Aims: We look for a nonlinear mechanism able to explain the persistence of MHD turbulence in shearing boxes with zero net magnetic flux, and potentially leading to large-scale dynamo action. Method: Spectral nonlinear simulations of the magnetorotational instability are shown to exhibit a large-scale axisymmetric magnetic field, maintained for a few orbits. The generation process of this field is investigated using the results of the simulations and an inhomogeneous linear approach. Results: The mechanism by which turbulence is sustained in zero-net-flux shearing boxes is shown to be related to the existence of a large-scale azimuthal field, surviving for several orbits. In particular, it is shown that MHD turbulence in shearing boxes can be seen as a dynamo process coupled to a magnetorotational-type instability.

We study the development of finite eccentricity in accretion disks in close binary systems using a two-dimensional grid-based numerical scheme. We perform detailed parameter studies to explore the dependence on viscosity, disk aspect ratio, the inclusion of a mass-transfer stream and the role of the boundary conditions. We consider mass ratios 0.05<q<0.3 appropriate to superoutbursting cataclysmic binary systems. Instability to the formation of a precessing eccentric disk that attains a quasi-steady state with mean eccentricity in the range 0.3-0.5 occurs readily. The shortest growth times are ~15 binary orbits for the largest viscosities and the instability mechanism is for the most part consistent with the mode-coupling mechanism associated with the 3:1 resonance proposed by Lubow. However, the results are sensitive to the treatment of the inner boundary and to the incorporation of the mass-transfer stream. In the presence of a stream we found a critical viscosity below which the disk remains circular. Incorporation of a mass-transfer stream tends to impart stability for small enough viscosity (or, equivalently, mass-transfer rate through the disk) and does assist in obtaining a prograde precession rate that is in agreement with observations. For the larger q the location of the 3:1 resonance is pushed outwards towards the Roche lobe where higher-order mode couplings and nonlinearity occur. It is likely that three-dimensional simulations that properly resolve the disk's vertical structure are required to make significant progress in this case.

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.

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non-rotating wall-bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non-axisymmetric instabilities of shear-induced axisymmetric toroidal velocity or magnetic fields, such as Kelvin-Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high-Reynolds number shear flows could act as a large-scale driving mechanism of turbulent flows that would in turn generate an independent small-scale dynamo.

This paper examines the onset of the viscous overstability in dense particulate rings. First, we formulate a dense gas kinetic theory that is applicable to the Saturnian system. Our model is essentially that of Araki and Tremaine (1986), which we show can be both simplified and generalised. Second, we put this model to work computing the equilibrium properties of dense planetary rings, which we subsequently compare with the results of N-body simulations, namely those of Salo (1991). Finally, we present the linear stability analyses of these equilibrium states, and derive criteria for the onset of viscous overstability in the self-gravitating and non-self-gravitating cases. These are framed in terms of particle size, orbital frequency, optical depth, and the parameters of the collision law. Our results compare favourably with the simulations of Salo et al. (2001). The accuracy and practicality of the continuum model we develop encourages its general use in future investigations of nonlinear phenomena.

According to one model, high-frequency quasi-periodic oscillations (QPOs) can be identified with inertial waves, trapped in the inner regions of accretion discs around black holes due to relativistic effects. In order to be detected, their amplitudes need to reach large enough values via some excitation mechanism. We work out in detail a non-linear coupling mechanism suggested by Kato, in which a global warping or eccentricity of the disc has a fundamental role. These large-scale deformations combine with trapped modes to generate `intermediate' waves of negative energy that are damped as they approach either their corotation resonance or the inner edge of the disc, resulting in amplification of the trapped waves. We determine the growth rates of the inertial modes, as well as their dependence on the spin of the black hole and the properties of the disc. Our results indicate that this coupling mechanism can provide an efficient excitation of trapped inertial waves, provided the global deformations reach the inner part of the disc with non-negligible amplitude.

Maxwell's investigations into the stability of Saturn's rings provide one of the earliest analyses of the dynamics of astrophysical discs. Current research in planetary rings extends Maxwell's kinetic theory to treat dense granular gases of particles undergoing moderately frequent inelastic collisions. Rather than disrupting the rings, local instabilities may be responsible for generating their irregular radial structure. Accretion discs around black holes or compact stars consist of a plasma permeated by a tangled magnetic field and may be compared with laboratory fluids through an analogy that connects Maxwell's researches in electromagnetism and viscoelasticity. A common theme in this work is the appearance of a complex fluid with a dynamical constitutive equation relating the stress in the medium to the history of its deformation.

We describe an instability of viscoelastic Couette-Taylor flow that is directly analogous to the magnetorotational instability (MRI) in astrophysical magnetohydrodynamics, with polymer molecules playing the role of magnetic field lines. By determining the conditions required for the onset of instability and the properties of the preferred modes, we distinguish it from the centrifugal and elastic instabilities studied previously. Experimental demonstration and investigation should be much easier for the viscoelastic instability than for the MRI in a liquid metal. The analogy holds with the case of a predominantly toroidal magnetic field such as is expected in an accretion disk and it may be possible to access a turbulent regime in which many modes are unstable.

A three-dimensional nonlinear dynamo process is identified in rotating plane Couette flow in the Keplerian regime. It is analogous to the hydrodynamic self-sustaining process in non-rotating shear flows and relies on the magneto-rotational instability of a toroidal magnetic field. Steady nonlinear solutions are computed numerically for a wide range of magnetic Reynolds numbers but are restricted to low Reynolds numbers. This process may be important to explain the sustenance of coherent fields and turbulent motions in Keplerian accretion disks, where all its basic ingredients are present.

We calculate the excitation and dissipation of low-frequency tidal oscillations in uniformly rotating solar-type stars. For tidal frequencies smaller than twice the spin frequency, inertial waves are excited in the convective envelope and are dissipated by turbulent viscosity. Enhanced dissipation occurs over the entire frequency range rather than in a series of very narrow resonant peaks, and is relatively insensitive to the effective viscosity. Hough waves are excited at the base of the convective zone and propagate into the radiative interior. We calculate the associated dissipation rate under the assumption that they do not reflect coherently from the center of the star. Tidal dissipation in a model based on the present Sun is significantly enhanced through the inclusion of the Coriolis force but may still fall short of that required to explain the circularization of close binary stars. However, the dependence of the results on the spin frequency, tidal frequency, and stellar model indicate that a more detailed evolutionary study including inertial and Hough waves is required. We also discuss the case of higher tidal frequencies appropriate to stars with very close planetary companions. The survival of even the closest hot Jupiters can be plausibly explained provided that the Hough waves they generate are not damped at the center of the star. We argue that this is the case because the tide excited by a hot Jupiter in the present Sun would marginally fail to achieve nonlinearity. As conditions at the center of the star evolve, nonlinearity may set in at a critical age, resulting in a relatively rapid inspiral of the hot Jupiter.

Subcritical transition to turbulence in Keplerian accretion disks is still a controversial issue and some theoretical progress is required in order to determine whether or not this scenario provides a plausible explanation for the origin of angular momentum transport in non-magnetized accretion disks. Motivated by the recent discoveries of exact nonlinear steady self-sustaining solutions in linearly stable non-rotating shear flows, we attempt to compute similar solutions in Rayleigh-stable rotating plane Couette flows and to identify transition mechanisms in such flows by combining nonlinear continuation methods and asymptotic theory. We obtain exact nonlinear solutions for Rayleigh-stable cyclonic regimes but show that it is not possible to compute solutions for Rayleigh-stable anticyclonic regimes, including Keplerian flow, using similar techniques. We also present asymptotic descriptions of these various problems at large Reynolds numbers that provide some insight into the differences between the non-rotating and Rayleigh-stable anticyclonic regimes and derive some necessary conditions for mechanisms analogous to the non-rotating self-sustaining process to be present in flows on the Rayleigh line. Our results demonstrate that subcritical transition mechanisms cannot be identified in wall-bounded Rayleigh-stable anticyclonic shear flows by transposing directly the phenomenology of subcritical transition in cyclonic and non-rotating wall-bounded shear flows. Asymptotic developments, however, leave open the possibility that nonlinear self-sustaining solutions may exist in unbounded or periodic flows on the Rayleigh line. These could serve as a starting point to discover solutions in Rayleigh-stable flows, but the nonlinear stability of Keplerian accretion disks remains to be determined.

Mean-motion resonances between a Keplerian disc and an orbiting companion are analysed within a Hamiltonian formulation using complex canonical Poincare variables, which are ideally suited to the description of eccentricity and inclination dynamics. Irreversibility is introduced by allowing for dissipation within the disc. A method is given for determining the rates of change of eccentricity and inclination variables of the disc and companion associated with resonances of various orders, including both reversible and irreversible effects, which extend and generalize previous results. Preliminary applications to protoplanetary systems and close binary stars are discussed.

We investigate the growth or decay rate of the fundamental mode of even symmetry in a viscous accretion disc. This mode occurs in eccentric discs and is known to be potentially overstable. We determine the vertical structure of the disc and its modes, treating radiative energy transport in the diffusion approximation. In the limit of very long radial wavelength, an analytical criterion for viscous overstability is obtained, which involves the effective shear and bulk viscosity, the adiabatic exponent and the opacity law of the disc. This differs from the prediction of a two-dimensional model. On shorter wavelengths (a few times the disc thickness), the criterion for overstability is more difficult to satisfy because of the different vertical structure of the mode. In a low-viscosity disc a third regime of intermediate wavelengths appears, in which the overstability is suppressed as the horizontal velocity perturbations develop significant vertical shear. We suggest that this effect determines the damping rate of eccentricity in protoplanetary discs, for which the long-wavelength analysis is inapplicable and overstability is unlikely to occur on any scale. In thinner accretion discs and in decretion discs around Be stars overstability may occur only on the longest wavelengths, leading to the preferential excitation of global eccentric modes.

Irregular structure in planetary rings is often attributed to the intrinsic instabilities of a homogeneous state undergoing Keplerian shear. Previously these have been analysed with simple hydrodynamic models. We instead employ a kinetic theory, in which we solve the linearised moment equations derived in Shu and Stewart 1985 for a dilute ring. This facilitates an examination of velocity anisotropy and non-Newtonian stress, and their effects on the viscous and viscous/gravitational instabilities thought to occur in Saturn's rings. Because we adopt a dilute gas model, the applicability of our results to the actual dense rings of Saturn are significantly curtailled. Nevertheless this study is a necessary preliminary before an attack on the difficult problem of dense ring dynamics. We find the Shu and Stewart formalism admits analytic stability criteria for the viscous overstability, viscous instability, and thermal instability. These criteria are compared with those of a hydrodynamic model incorporating the effective viscosity and cooling function computed from the kinetic steady state. We find the two agree in the `hydrodynamic limit' (i.e. many collisions per orbit) but disagree when collisions are less frequent, when we expect the viscous stress to be increasingly non-Newtonian and the velocity distribution increasingly anisotropic. In particular, hydrodynamics predicts viscous overstability for a larger portion of parameter space. We also numerically solve the linearised equations of the more accurate Goldreich and Tremaine 1978 kinetic model and discover its linear stability to be qualitatively the same as that of Shu and Stewart's. Thus the simple collision operator adopted in the latter would appear to be an adequate approximation for dilute rings, at least in the linear regime.

The migration of a planet through a gaseous disc causes the locations of their resonant interactions to drift and can alter the torques exerted between the planet and the disc. We analyse the time-dependent dynamics of a non-coorbital corotation resonance under these circumstances. The ratio of the resonant torque in a steady state to the value given by Goldreich & Tremaine (1979) depends essentially on two dimensionless quantities: a dimensionless turbulent diffusion time-scale and a dimensionless radial drift speed. When the drift speed is comparable to the libration speed and the viscosity is small, the torque can become much larger than the unsaturated value in the absence of migration, but is still proportional to the large-scale vortensity gradient in the disc. Fluid that is trapped in the resonance and drifts with it acquires a vortensity anomaly relative to its surroundings. If the anomaly is limited by viscous diffusion in a steady state, the resulting torque is inversely proportional to the viscosity, although a long time may be required to achieve this state. A further, viscosity-independent, contribution to the torque comes from fluid that streams through the resonant region. In other cases, torque oscillations occur before the steady value is achieved. We discuss the significance of these results for the evolution of eccentricity in protoplanetary systems. We also describe the possible application of these findings to the coorbital region and the concept of runaway (or type III) migration. [Abridged]

We have applied an eccentric accretion disc theory in simplified form to the case of an accretion disc in a binary system, where the disc contains the 3:1 Lindblad resonance. This is relevant to the case of superhumps in SU Ursae Majoris cataclysmic variables and other systems, where it is thought that this resonance leads to growth of eccentricity and a modulation in the light curve due to the interaction of a precessing eccentric disc with tidal stresses. A single differential equation is formulated which describes the propagation, resonant excitation and viscous damping of eccentricity. The theory is first worked out in the simple case of a narrow ring and leads to the conclusion that the eccentricity distribution is locally suppressed by the presence of the resonance, creating a dip in the eccentricity at the resonant radius. Application of this theory to the superhump case confirms this conclusion and produces a more accurate expression for the precession rate of the disc than has been previously accomplished with simple dynamical estimates.

The nonlinear dynamics of a warped accretion disc is investigated in the important case of a thin Keplerian disc with negligible viscosity and self-gravity. A one-dimensional evolutionary equation is formally derived that describes the primary nonlinear and dispersive effects on propagating bending waves other than parametric instabilities. It has the form of a derivative nonlinear Schroedinger equation with coefficients that are obtained explicitly for a particular model of a disc. The properties of this equation are analysed in some detail and illustrative numerical solutions are presented. The nonlinear and dispersive effects both depend on the compressibility of the gas through its adiabatic index Gamma. In the physically realistic case Gamma<3, nonlinearity does not lead to the steepening of bending waves but instead enhances their linear dispersion. In the opposite case Gamma>3, nonlinearity leads to wave steepening and solitary waves are supported. The effects of a small effective viscosity, which may suppress parametric instabilities, are also considered. This analysis may provide a useful point of comparison between theory and numerical simulations of warped accretion discs.

We consider the alignment torque between a spinning black hole and an accretion disc whose angular momenta are misaligned. This situation must hold initially in almost all gas accretion events on to supermassive black holes, and may occur in binaries where the black hole receives a natal supernova kick. We show that the torque always acts to align the hole's spin with the total angular momentum without changing its magnitude. The torque acts dissipatively on the disc, reducing its angular momentum, and aligning it with the hole if and only if the angle theta between the angular momenta J_d of the disc and J_h of the hole satisfies the inequality cos theta > -J_d / 2 J_h. If this condition fails, which requires both theta > pi/2 and J_d < 2 J_h, the disc counteraligns.

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.

We investigate the nonlinear dynamics of turbulent shear flows, with and without rotation, in the context of a simple but physically motivated closure of the equation governing the evolution of the Reynolds stress tensor. We show that the model naturally accounts for some familiar phenomena in parallel shear flows, such as the subcritical transition to turbulence at a finite Reynolds number and the occurrence of a universal velocity profile close to a wall at large Reynolds number. For rotating shear flows we find that, depending on the Rayleigh discriminant of the system, the model predicts either linear instability or nonlinear instability or complete stability as the Reynolds number is increased to large values. We investigate the properties of Couette-Taylor flows for varying inner and outer cylinder rotation rates and identify the region of linear instability (similar to Taylor's), as well as regions of finite-amplitude instability qualitatively compatible with recent experiments. We also discuss quantitative predictions of the model in comparison with a range of experimental torque measurements. Finally, we consider the relevance of this work to the question of the hydrodynamic stability of astrophysical accretion discs.

[Abridged] The torque exerted by an external potential on a two-dimensional gaseous disk at non-co-orbital corotation resonances is studied by means of numerical simulations. The degree of saturation of these resonances is important in determining whether an eccentric giant planet embedded in a protoplanetary disk experiences an eccentricity excitation or damping. We perform calculations restricted to one or two resonances to investigate the properties of two neighboring corotation resonances, as well as the properties of a corotation resonance that overlaps a Lindblad resonance. We find that these properties hardly differ from the case of an isolated corotation resonance. In particular, although the torque of two neighboring corotation resonances may differ from the sum of the torques of the corresponding resonances considered as isolated, it never exceeds the sum of the fully unsaturated isolated corotation resonances, and it saturates in a fashion similar to an isolated resonance. Similarly, the presence of an underlying Lindblad resonance hardly affects the corotation torque, even if that resonance implies a torque strong enough to significantly redistribute the azimuthally averaged surface density profile, in which case the corotation torque scales with the resulting vortensity gradient. This set of numerical experiments thus essentially validates previous analytic studies. As a side result, we show that corotation libration islands misrepresented by a mesh of too low resolution can lead to a strongly overestimated corotation torque. This may be an explanation why the eccentricity of embedded Jupiter-sized planets was never observed to undergo an excitation in the numerical simulations performed so far.

[Abridged] Tides may play an important role in determining the observed distributions of mass, orbital period, and eccentricity of the extrasolar planets. In addition, tidal interactions between giant planets in the solar system and their moons are thought to be responsible for the orbital migration of the satellites, leading to their capture into resonant configurations. We treat the underlying fluid dynamical problem with the aim of determining the efficiency of tidal dissipation in gaseous giant planets. In cases of interest, the tidal forcing frequencies are comparable to the spin frequency of the planet but small compared to its dynamical frequency. We therefore study the linearized response of a slowly and possibly differentially rotating planet to low-frequency tidal forcing. Convective regions of the planet support inertial waves, while any radiative regions support generalized Hough waves. We present illustrative numerical calculations of the tidal dissipation rate and argue that inertial waves provide a natural avenue for efficient tidal dissipation in most cases of interest. The resulting value of Q depends in a highly erratic way on the forcing frequency, but we provide evidence that the relevant frequency-averaged dissipation rate may be asymptotically independent of the viscosity in the limit of small Ekman number. In short-period extrasolar planets, if the stellar irradiation of the planet leads to the formation of a radiative outer layer that supports generalized Hough modes, the tidal dissipation rate can be enhanced through the excitation and damping of these waves. These dissipative mechanisms offer a promising explanation of the historical evolution and current state of the Galilean satellites as well as the observed circularization of the orbits of short-period extrasolar planets.

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from one Earth mass (1 M$_e$) to one Jupiter mass (1 M$_J$) by using the ZEUS hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses $M \gsim 0.1$ M$_J$ produce significant perturbations in the disc's surface density. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc midplane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the planet's orbit. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M$_J$ and is highly efficient, occurring at the local viscous rate. The migration timescales for planets of mass less than 0.1 M$_J$, based on torques from disc material outside the planets' Roche lobes, are in excellent agreement with the linear theory of Type I (non-gap) migration for three-dimensional discs. The transition from Type I to Type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M$_J$ with little migration. Planets with final masses of order 10 M$_J$ would undergo large migration, which makes formation and survival difficult.

We determine the torque exerted in a steady state by an external potential on a three-dimensional gaseous disk at a non-coorbital corotation resonance. Our model accounts for the feedback of the torque on the surface density and vorticity in the corotation region, and assumes that the disk has a barotropic equation of state and a nonzero effective viscosity. The ratio of the torque to the value given by the formula of Goldreich & Tremaine depends essentially on a single dimensionless parameter, which quantifies the extent to which the resonance is saturated. We discuss the implications for the eccentricity evolution of young planets.

We demonstrate a close analogy between a viscoelastic medium and an electrically conducting fluid containing a magnetic field. Specifically, the dynamics of the Oldroyd-B fluid in the limit of large Deborah number corresponds to that of a magnetohydrodynamic (MHD) fluid in the limit of large magnetic Reynolds number. As a definite example of this analogy, we compare the stability properties of differentially rotating viscoelastic and MHD flows. We show that there is an instability of the Oldroyd-B fluid that is physically distinct from both the inertial and elastic instabilities described previously in the literature, but is directly equivalent to the magnetorotational instability in MHD. It occurs even when the specific angular momentum increases outwards, provided that the angular velocity decreases outwards; it derives from the kinetic energy of the shear flow and does not depend on the curvature of the streamlines. However, we argue that the elastic instability of viscoelastic Couette flow has no direct equivalent in MHD.

[Abridged] The turbulent stresses that lead to angular momentum transport in accretion discs have often been treated as resulting from an isotropic effective viscosity, related to the pressure through the alpha parametrization of Shakura & Sunyaev. This simple approach may be adequate for the simplest aspects of accretion disc theory, and was necessitated historically by an incomplete understanding of the origin of the turbulence. More recently, Balbus & Hawley have shown that the magnetorotational instability provides a robust mechanism of generating turbulent Reynolds and Maxwell stresses in sufficiently ionized discs. The alpha viscosity model fails to describe numerous aspects of this process. This paper introduces a new analytical model that aims to represent more faithfully the dynamics of magnetorotational turbulent stresses and bridge the gap between analytical studies and numerical simulations. Covariant evolutionary equations for the mean Reynolds and Maxwell tensors are presented, which correctly include the linear interaction with the mean flow. Non-linear and dissipative effects, in the absence of an imposed magnetic flux and in the limit of large Reynolds number and magnetic Reynolds number, are modelled through five non-linear terms that represent known physical processes and are strongly constrained by symmetry properties and dimensional considerations. The resulting model should be particularly useful in understanding the dynamics of warped, eccentric and tidally distorted discs, non-Keplerian accretion flows close to black holes, and a variety of time-dependent accretion phenomena.

We present three-dimensional numerical simulations of the interaction of a circular-orbit planet with a protoplanetary disk. We calculate the flow pattern, the accretion rate, and torques on the planet. We consider planet masses ranging from 1 Earth mass to 1 Jupiter mass.

We study the viscous effects on the interaction between the coplanar Be-star disc and the neutron star in Be/X-ray binaries, using a three-dimensional, smoothed particle hydrodynamics code. For simplicity, we assume the Be disc to be isothermal at the temperature of half the stellar effective temperature. In order to mimic the gas ejection process from the Be star, we inject particles with the Keplerian rotation velocity at a radius just outside the star. Both Be star and neutron star are treated as point masses. We find that the Be-star disc is effectively truncated if the Shakura-Sunyaev viscosity parameter alpha_SS >> 1, which confirms the previous semi-analytical result. In the truncated disc, the material decreted from the Be star accumulates, so that the disc becomes denser more rapidly than if around an isolated Be star. The resonant truncation of the Be disc results in a significant reduction of the amount of gas captured by the neutron star and a strong dependence of the mass capture rate on the orbital phase. We also find that an eccentric mode is excited in the Be disc through direct driving due to a one-armed bar potential of the binary. The strength of the mode becomes greater in the case of a smaller viscosity. In a high-resolution simulation with alpha_SS=0.1, the eccentric mode is found to precess in a prograde sense. The mass capture rate by the neutron star modulates as the mode precesses.

We consider the evolution of a warped disc around a Kerr black hole, under conditions such that the warp propagates in a wavelike manner. This occurs when the dimensionless effective viscosity, alpha, that damps the warp is less than the characteristic angular semi-thickness, H/R, of the disc. We adopt linearized equations that are valid for warps of sufficiently small amplitude in a Newtonian disc, but also account for the apsidal and nodal precession that occur in the Kerr metric. Through analytical and time-dependent studies, we confirm the results of Demianski & Ivanov, and of Ivanov & Illarionov, that such a disc takes on a characteristic warped shape. The inner part of the disc is not necessarily aligned with the equator of the hole, even in the presence of dissipation. We draw attention to the fact that this might have important implications for the directionality of jets emanating from discs around rotating black holes.

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic ($\gamma = 5/3$). Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance ${\cal M}_{\rm r} \ga 0.01$ undergo shocks within a distance of order the resonance radius.

An important and widely neglected aspect of the interaction between an accretion disc and a massive companion with a coplanar orbit is the vertical component of the tidal force. As shown by Lubow, the response of the disc to vertical forcing is resonant at certain radii, at which a localized torque is exerted, and from which a compressive wave (p mode) may be emitted. Although these vertical resonances are weaker than the corresponding Lindblad resonances, the m=2 inner vertical resonance in a binary star is typically located within the tidal truncation radius of a circumstellar disc. In this paper I develop a general theory of vertical resonances, allowing for non-linearity of the response, and dissipation by radiative damping and turbulent viscosity. The problem is reduced to a universal, non-linear ordinary differential equation with two real parameters. Solutions of the complex non-linear Airy equation are presented to illustrate the non-linear saturation of the resonance and the effects of dissipation. It is argued that the m=2 inner vertical resonance is unlikely to truncate the disc in cataclysmic variable stars, but contributes to angular momentum transport and produces a potentially observable non-axisymmetric structure.

A planet of low mass orbiting in a two-dimensional gaseous disc generates a one-armed spiral wake. We explain this phenomenon as the result of constructive interference between wave modes in the disc, somewhat similar to the Kelvin wedge produced in the wake of a ship. The same feature is not expected in a three-dimensional disc with thermal stratification.

The non-axisymmetric features observed in the discs of dwarf novae in outburst are usually considered to be spiral shocks, which are the non-linear relatives of tidally excited waves. This interpretation suffers from a number of problems. For example, the natural site of wave excitation lies outside the Roche lobe, the disc must be especially hot, and most treatments of wave propagation do not take into account the vertical structure of the disc. In this paper I construct a detailed semi-analytical model of the non-linear tidal distortion of a thin, three-dimensional accretion disc by a binary companion on a circular orbit. The analysis presented here allows for vertical motion and radiative energy transport, and introduces a simple model for the turbulent magnetic stress. The m=2 inner vertical resonance has an important influence on the amplitude and phase of the tidal distortion. I show that the observed patterns find a natural explanation if the emission is associated with the tidally thickened sectors of the outer disc, which may be irradiated from the centre. According to this hypothesis, it may be possible to constrain the physical parameters of the disc through future observations.

In a planetary system, a secular particle resonance occurs at a location where the precession rate of a test particle (e.g. an asteroid) matches the frequency of one of the precessional modes of the planetary system. We investigate the secular interactions of a system of mutually inclined planets with a gaseous protostellar disk that may contain a secular nodal particle resonance. We determine the normal modes of some mutually inclined planet-disk systems. The planets and disk interact gravitationally, and the disk is internally subject to the effects of gas pressure, self-gravity, and turbulent viscosity. The behavior of the disk at a secular resonance is radically different from that of a particle, owing mainly to the effects of gas pressure. The resonance is typically broadened by gas pressure to the extent that global effects, including large-scale warps, dominate. The standard resonant torque formula is invalid in this regime. Secular interactions cause a decay of the inclination at a rate that depends on the disk properties, including its mass, turbulent viscosity, and sound speed. For a Jupiter-mass planet embedded within a minimum-mass solar nebula having typical parameters, dissipation within the disk is sufficient to stabilize the system against tilt growth caused by mean-motion resonances.

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.

A new theory of eccentric accretion discs is presented. Starting from the basic fluid-dynamical equations in three dimensions, I derive the fundamental set of one-dimensional equations that describe how the mass, angular momentum and eccentricity vector of a thin disc evolve as a result of internal stresses and external forcing. The analysis is asymptotically exact in the limit of a thin disc, and allows for slowly varying eccentricities of arbitrary magnitude. The theory is worked out in detail for a Maxwellian viscoelastic model of the turbulent stress in an accretion disc. This generalizes the conventional alpha viscosity model to account for the non-zero relaxation time of the turbulence, and is physically motivated by a consideration of the nature of magnetohydrodynamic turbulence. It is confirmed that circular discs are typically viscously unstable to eccentric perturbations, as found by Lyubarskij, Postnov & Prokhorov, if the conventional alpha viscosity model is adopted. However, the instability can usually be suppressed by introducing a sufficient relaxation time and/or bulk viscosity. It is then shown that an initially uniformly eccentric disc does not retain its eccentricity as had been suggested by previous analyses. The evolutionary equations should be useful in many applications, including understanding the origin of planetary eccentricities and testing theories of quasi-periodic oscillations in X-ray binaries.

The launching of magnetohydrodynamic outflows from accretion disks is considered. We formulate a model for the local vertical structure of a thin disk threaded by a poloidal magnetic field of dipolar symmetry. The model consists of an optically thick disk matched to an isothermal atmosphere. The disk is supposed to be turbulent and possesses an effective viscosity and an effective magnetic diffusivity. In the atmosphere, if the magnetic field lines are inclined sufficiently to the vertical, a magnetocentrifugal outflow is driven and passes through a slow magnetosonic point close to the surface. We determine how the rate of mass loss varies with the strength and inclination of the magnetic field. In particular, we find that for disks in which the mean poloidal field is sufficiently strong to stabilize the disk against the magnetorotational instability, the mass loss rate decreases extremely rapidly with increasing field strength, and is maximal at an inclination angle of 40-50 degrees. For turbulent disks with weaker mean fields, the mass loss rate increases monotonically with increasing strength and inclination of the field, but the solution branch terminates before achieving excessive mass loss rates. Our results suggest that efficient jet launching occurs for a limited range of field strengths, and a limited range of inclination angles in excess of 30 degrees. In addition, we determine the direction and rate of radial migration of the poloidal magnetic flux, and discuss whether configurations suitable for jet launching can be maintained against dissipation.

We study the radiation-driven warping of accretion discs in the context of X-ray binaries. The latest evolutionary equations are adopted, which extend the classical alpha theory to time-dependent thin discs with non-linear warps. We also develop accurate, analytical expressions for the tidal torque and the radiation torque, including self-shadowing. We investigate the possible non-linear dynamics of the system within the framework of bifurcation theory. First, we re-examine the stability of an initially flat disc to the Pringle instability. Then we compute directly the branches of non-linear solutions representing steadily precessing discs. Finally, we determine the stability of the non-linear solutions. Each problem involves only ordinary differential equations, allowing a rapid, accurate and well resolved solution. We find that radiation-driven warping is probably not a common occurrence in low-mass X-ray binaries. We also find that stable, steadily precessing discs exist for a narrow range of parameters close to the stability limit. This could explain why so few systems show clear, repeatable `super-orbital' variations. The best examples of such systems, Her X-1, SS 433 and LMC X-4, all lie close to the stability limit for a reasonable choice of parameters. Systems far from the stability limit, including Cyg X-2, Cen X-3 and SMC X-1, probably experience quasi-periodic or chaotic variability as first noticed by Wijers & Pringle. We show that radiation-driven warping provides a coherent and persuasive framework but that it does not provide a generic explanation for the long-term variabilities in all X-ray binaries.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

The non-linear fluid dynamics of a warped accretion disc was investigated in an earlier paper by developing a theory of fully non-linear bending waves in a thin, viscous disc. That analysis is here extended to take proper account of thermal and radiative effects by solving an energy equation that includes viscous dissipation and radiative transport. The problem is reduced to simple one-dimensional evolutionary equations for mass and angular momentum, expressed in physical units and suitable for direct application. This result constitutes a logical generalization of the alpha theory of Shakura & Sunyaev to the case of a time-dependent warped accretion disc. The local thermal-viscous stability of such a disc is also investigated.

We consider the stability of warping modes in Keplerian discs. We find them to be parametrically unstable using two lines of attack, one based on three-mode couplings and the other on Floquet theory. We confirm the existence of the instability, and investigate its nonlinear development in three dimensions, via numerical experiment. The most rapidly growing non-axisymmetric disturbances are the most nearly axisymmetric (low m) ones. Finally, we offer a simple, somewhat speculative model for the interaction of the parametric instability with the warp. We apply this model to the masing disc in NGC 4258 and show that, provided the warp is not forced too strongly, parametric instability can fix the amplitude of the warp.

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, $\alpha_{\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, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\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 $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = 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.

We consider the dynamics of a protostellar disk surrounding a star in a circular-orbit binary system. Our aim is to determine whether, if the disk is initially tilted with respect to the plane of the binary orbit, the inclination of the system will increase or decrease with time. The problem is formulated in the binary frame in which the tidal potential of the companion star is static. We consider a steady, flat disk that is aligned with the binary plane and investigate its linear stability with respect to tilting or warping perturbations. The dynamics is controlled by the competing effects of the m=0 and m=2 azimuthal Fourier components of the tidal potential. In the presence of dissipation, the m=0 component causes alignment of the system, while the m=2 component has the opposite tendency. We find that disks that are sufficiently large, in particular those that extend to their tidal truncation radii, are generally stable and will therefore tend to alignment with the binary plane on a time-scale comparable to that found in previous studies. However, the effect of the m=2 component is enhanced in the vicinity of resonances where the outer radius of the disk is such that the natural frequency of a global bending mode of the disk is equal to twice the binary orbital frequency. Under such circumstances, the disk can be unstable to tilting and acquire a warped shape, even in the absence of dissipation. The outer radius corresponding to the primary resonance is always smaller than the tidal truncation radius. For disks smaller than the primary resonance, the m=2 component may be able to cause a very slow growth of inclination through the effect of a near resonance that occurs close to the disk center. We discuss these results in the light of recent observations of protostellar disks in binary systems.

Differentially rotating, "advection-dominated" accretion flows are considered in which the heat generated by viscous dissipation is retained in the fluid. The equations of time-dependent quasi-spherical accretion are solved in a simplified one-dimensional model that neglects the latitudinal dependence of the flow. A self-similar solution is presented that has finite size, mass, angular momentum and energy. This may be expected to be an attractor for the initial-value problem in which a cool and narrow ring of fluid orbiting around a central mass heats up, spreads radially and is accreted. The solution provides some insight into the dynamics of quasi-spherical accretion and avoids many of the strictures of the steady self-similar solution of Narayan & Yi. Special attention is given to the astrophysically important case in which the adiabatic exponent gamma=5/3; even in this case, the flow is found to be differentially rotating and bound to the central object, and accretion can occur without the need for powerful outflows.

The dynamics of a viscous accretion disc subject to a slowly varying warp of large amplitude is considered. Attention is restricted to discs in which self-gravitation is negligible, and to the generic case in which the resonant wave propagation found in inviscid Keplerian discs does not occur. The equations of fluid dynamics are derived in a coordinate system that follows the principal warping motion of the disc. They are reduced using asymptotic methods for thin discs, and solved to extract the equation governing the warp. In general, this is a wave equation of parabolic type with non-linear dispersion and diffusion, which describes fully non-linear bending waves. This method generalizes the linear theory of Papaloizou & Pringle (1983) to allow for an arbitrary rotation law, and extends it into the non-linear domain, where it connects with a generalized version of the theory of Pringle (1992). The astrophysical implications of this analysis are discussed briefly.

A new description of the dynamics of warped accretion discs is presented. A theory of fully nonlinear, slowly varying bending waves is developed, involving a proper treatment of viscous fluid dynamics but neglecting self-gravitation. The special case of resonant wave propagation found in inviscid Keplerian discs is not considered. The resulting description goes beyond existing linear and nonlinear theories, and is suitable for numerical implementation.

We critically assess the role of the Blandford-Znajek mechanism in the powering of outflows from accretion disk-fed black holes. We argue that there is no reason to suppose that the magnetic field threading the central spinning black hole differs significantly in strength from that threading the central regions of the disk. In this case, we show that the electromagnetic output from the inner disk regions is expected in general to dominate over that from the hole. Thus the spin (or not) of the hole is probably irrelevant to the expected electromagnetic power output from the system. We also point out that the strength of the poloidal field in the center of a standard accretion disk has been generally overestimated, and discuss scenarios which might give rise to more significant central poloidal fields.

We extend our analysis of the three-dimensional response of a vertically polytropic disk to tidal forcing at Lindblad resonances by including the effects of a disk atmosphere. The atmosphere is modeled as an isothermal layer that joins smoothly on to an underlying polytropic layer. The launched wave progressively enters the atmosphere as it propagates away from the resonance. The wave never propagates vertically, however, and the wave energy rises to a (finite) characteristic height in the atmosphere. The increase of wave amplitude associated with this process of wave channeling is reduced by the effect of the atmosphere. For waves of large azimuthal mode number m generated by giant planets embedded in a disk, the increase in wave amplitude is still substantial enough to be likely to dissipate the wave energy by shocks for even modest optical depths (tau greater than about 10) over a radial distance of a few times the disk thickness. For low-m waves generated in circumstellar disks in binary stars, the effects of wave channeling are less important and the level of wave nonlinearity increases by less than a factor of 10 in going from the disk edge to the disk center. For circumbinary disks, the effects of wave channeling remain important, even for modest values of optical depth.

We analyze the linear, 3D response to tidal forcing of a disk that is thin and thermally stratified in the direction normal to the disk plane. We model the vertical disk structure locally as a polytrope which represents a disk of high optical depth. We solve the 3D gas-dynamic equations semi-analytically in the neighborhood of a Lindblad resonance. These solutions match asymptotically on to those valid away from resonances and provide solutions valid at all radii. We obtain the following results. 1) A variety of waves are launched at resonance. However, the f mode carries more than 95% of the torque exerted at the resonance. 2) These 3D waves collectively transport exactly the amount of angular momentum predicted by the 2D torque formula. 3) Near resonance, the f mode occupies the full vertical extent of the disk. Away from resonance, the f mode becomes confined near the surface of the disk, and, in the absence of other dissipation mechanisms, damps via shocks. The radial length scale for this process is roughly r_L/m (for resonant radius r_L and azimuthal wavenumber m), independent of the disk thickness H. This wave channeling process is due to the variations of physical quantities in r and is not due to wave refraction. 4) However, the inwardly propagating f mode launched from an m=2 inner Lindblad resonance experiences relatively minor channeling. We conclude that for binary stars, tidally generated waves in highly optically thick circumbinary disks are subject to strong nonlinear damping by the channeling mechanism, while those in circumstellar accretion disks are subject to weaker nonlinear effects. We also apply our results to waves excited by young planets for which m is approximately r/H and conclude that the waves are damped on the scale of a few H.

The theory of waves and instabilities in a differentially rotating disc containing a poloidal magnetic field is developed within the framework of ideal magnetohydrodynamics. A continuous spectrum, for which the eigenfunctions are localized on individual magnetic surfaces, is identified but is found not to contain any instabilities associated with differential rotation. The normal modes of a weakly magnetized thin disc are studied by extending the asymptotic methods used previously to describe the equilibria. Waves propagate radially in the disc according to a dispersion relation which is determined by solving an eigenvalue problem at each radius. The dispersion relation for a hydrodynamic disc is re-examined and the modes are classified according to their behaviour in the limit of large wavenumber. The addition of a magnetic field introduces new, potentially unstable, modes and also breaks up the dispersion diagram by causing avoided crossings. The stability boundary to the magnetorotational instability in the parameter space of polytropic equilibria is located by solving directly for marginally stable equilibria. For a given vertical magnetic field in the disc, bending of the field lines has a stabilizing effect and it is shown that stable equilibria exist which are capable of launching a predominantly centrifugally driven wind.

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.