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 vertical temperature structure of a protoplanetary disk bears on several processes relevant to planet formation, such as gas and dust grain chemistry, ice lines and convection. The temperature profile is controlled by irradiation from the central star and by any internal source of heat such as might arise from gas accretion. We investigate the heat and angular momentum transport generated by the resistive dissipation of magnetic fields in laminar disks. We use local one-dimensional simulations to obtain vertical temperature profiles for typical conditions in the inner disk (0.5 to 4 au). Using simple assumptions for the gas ionization and opacity, the heating and cooling rates are computed self-consistently in the framework of radiative non-ideal magnetohydrodynamics. We characterize steady solutions that are symmetric about the midplane and which may be associated with saturated Hall-shear unstable modes. We also examine the dissipation of electric currents driven by global accretion-ejection structures. In both cases we obtain significant heating for a sufficiently high opacity. Strong magnetic fields can induce an order-unity temperature increase in the disk midplane, a convectively unstable entropy profile, and a surface emissivity equivalent to a viscous heating of $\alpha \sim 10^{-2}$. These results show how magnetic fields may drive efficient accretion and heating in weakly ionized disks where turbulence might be inefficient, at least for a range of radii and ages of the disk.

Gravitational instability (GI) controls the dynamics of young massive protoplanetary discs. Apart from facilitating gas accretion on to the central protostar, it must also impact on the process of planet formation: directly through fragmentation, and indirectly through the turbulent concentration of small solids. To understand the latter process, it is essential to determine the dust dynamics in such a turbulent flow. For that purpose, we conduct a series of 3D shearing box simulations of coupled gas and dust, including the gas's self-gravity and scanning a range of Stokes numbers, from 0.001 to ~0.2. First, we show that the vertical settling of dust in the midplane is significantly impeded by gravitoturbulence, with the dust scale-height roughly 0.6 times the gas scale height for centimetre grains. This is a result of the strong vertical diffusion issuing from (a) small-scale inertial-wave turbulence feeding off the GI spiral waves and (b) the larger-scale vertical circulations that naturally accompany the spirals. Second, we show that at R=50 AU concentration events involving sub-metre particles and yielding order 1 dust to gas ratios are rare and last for less than an orbit. Moreover, dust concentration is less efficient in 3D than in 2D simulations. We conclude that GI is not especially prone to the turbulent accumulation of dust grains. Finally, the large dust scale-height measured in simulations could be, in the future, compared with that of edge-on discs seen by ALMA, thus aiding detection and characterisation of GI in real systems.

In the early stages of a protoplanetary disk, when its mass is a significant fraction of its star's, turbulence generated by gravitational instability (GI) should feature significantly in the disk's evolution. At the same time, the disk may be sufficiently ionised for magnetic fields to play some role in the dynamics. Though usually neglected, the impact of magnetism on the GI may be critical, with consequences for several processes: the efficiency of accretion, spiral structure formation, fragmentation, and the dynamics of solids. In this paper, we report on global three-dimensional magnetohydrodynamical simulations of a self-gravitating protoplanetary disk using the meshless finite mass (MFM) Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta < 10), an order of magnitude greater than that generated by the magneto-rotational instability alone. We also determine the dynamo's nonlinear back reaction on the gravitoturbulent flow: the saturated state is substantially hotter, with an associated larger Toomre parameter and weaker, more 'flocculent' spirals. But perhaps of greater import is the dynamo's boosting of accretion via a significant Maxwell stress; mass accretion is enhanced by factors of several relative to either pure GI or pure MRI. Our simulations use ideal MHD, an admittedly poor approximation in protoplanetary disks, and thus future studies should explore the full gamut of non-ideal MHD. In preparation for that, we exhibit a small number of Ohmic runs that reveal that the dynamo, if anything, is stronger in a non-ideal environment. This work confirms that magnetic fields are a potentially critical ingredient in gravitoturbulent young disks, possibly controlling their evolution, especially via their enhancement of (potentially episodic) accretion.

The magneto-rotational instability (MRI) is one of the most important processes in sufficiently ionized astrophysical disks. Grid-based simulations, especially those using the local shearing box approximation, provide a powerful tool to study the ensuing nonlinear turbulence. On the other hand, while meshless methods have been widely used in both cosmology, galactic dynamics, and planet formation they have not been fully deployed on the MRI problem. We present local unstratified and vertically stratified MRI simulations with two meshless MHD schemes: a recent implementation of SPH MHD (Price2012), and a MFM MHD scheme with a constrained gradient divergence cleaning scheme, as implemented in the GIZMO code \citep{Hopkins2017}. Concerning variants of the SPH hydro force formulation we consider both the "vanilla" SPH and the PSPH variant included in GIZMO. We find, as expected, that the numerical noise inherent in these schemes affects turbulence significantly. A high order kernel, free of the pairing instability, is necessary. Both schemes can adequately simulate MRI turbulence in unstratified shearing boxes with net vertical flux. The turbulence, however, dies out in zero-net-flux unstratified boxes, probably due to excessive and numerical dissipation. In zero-net-flux vertically stratified simulations, MFM can reproduce the MRI dynamo and its characteristic butterfly diagram for several tens of orbits before ultimately decaying. In contrast, extremely strong toroidal fields, as opposed to sustained turbulence, develop in equivalent simulations using SPH MHD. This unphysical state in SPH MHD is likely caused by a combination of excessive artificial viscosity, numerical resistivity, and the relatively large residual errors in the divergence of the magnetic field remaining even after cleaning procedures are applied.

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.

The origin of large-scale and coherent magnetic fields in astrophysical discs is an important and long standing problem. Researchers commonly appeal to a turbulent dynamo, sustained by the magneto-rotational instability (MRI), to supply the large-scale field. But research over the last decade in particular has demonstrated that various non-ideal MHD effects can impede or extinguish the MRI, especially in protoplanetary disks. In this paper we propose a new scenario, by which the magnetic field is generated and sustained via the gravitational instability (GI). We use 3D stratified shearing box simulations to characterise the dynamo and find that it works at low magnetic Reynolds number (from unity to ~100) for a wide range of cooling times and boundary conditions. The process is kinematic, with a relatively fast growth rate ($< 0.1 \Omega$), and it shares some properties of mean field dynamos. The magnetic field is generated via the combination of differential rotation and spiral density waves, the latter providing compressible horizontal motions and large-scale vertical rolls. At greater magnetic Reynolds numbers the build up of large-scale field is diminished and instead small-scale structures emerge from the breakdown of twisted flux ropes. We propose that GI may be key to the dynamo engine not only in young protoplanetary discs but also in some AGN and galaxies.

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

The prevalence and consequences of convection perpendicular to the plane of accretion discs have been discussed for several decades. Recent simulations combining convection and the magnetorotational instability have given fresh impetus to the debate, as the interplay of the two processes can enhance angular momentum transport, at least in the optically thick outburst stage of dwarf novae. In this paper we seek to isolate and understand the most generic features of disc convection, and so undertake its study in purely hydrodynamical models. First, we investigate the linear phase of the instability, obtaining estimates of the growth rates both semi-analytically, using one-dimensional spectral computations, as well as analytically, using WKBJ methods. Next we perform three-dimensional, vertically stratified, shearing box simulations with the conservative, finite-volume code PLUTO, both with and without explicit diffusion coefficients. We find that hydrodynamic convection can, in general, drive outward angular momentum transport, a result that we confirm with ATHENA, an alternative finite-volume code. Moreover, we establish that the sign of the angular momentum flux is sensitive to the diffusivity of the numerical scheme. Finally, in sustained convection, whereby the system is continuously forced to an unstable state, we observe the formation of various coherent structures, including large- scale and oscillatory convective cells, zonal flows, and small vortices.

Via the process of accretion, magnetorotational turbulence removes energy from a disk's orbital motion and transforms it into heat. Turbulent heating is far from uniform and is usually concentrated in small regions of intense dissipation, characterised by abrupt magnetic reconnection and higher temperatures. These regions are of interest because they might generate non-thermal emission, in the form of flares and energetic particles, or thermally process solids in protoplanetary disks. Moreover, the nature of the dissipation bears on the fundamental dynamics of the magnetorotational instability (MRI) itself: local simulations indicate that the large-scale properties of the turbulence (e.g. saturation levels, the stress-pressure relationship) depend on the short dissipative scales. In this paper we undertake a numerical study of how the MRI dissipates and the small-scale dissipative structures it employs to do so. We use the Godunov code RAMSES and unstratified compressible shearing boxes. Our simulations reveal that dissipation is concentrated in ribbons of strong magnetic reconnection that are significantly elongated in azimuth, up to a scale height. Dissipative structures are hence meso-scale objects, and potentially provide a route by which large scales and small scales interact. We go on to show how these ribbons evolve over time --- forming, merging, breaking apart, and disappearing. Finally, we reveal important couplings between the large-scale density waves generated by the MRI and the small-scale structures, which may illuminate the stress-pressure relationship in MRI turbulence.

Spiral density waves dominate several facets of accretion disc dynamics --- planet-disc interactions and gravitational instability (GI) most prominently. Though they have been examined thoroughly in two-dimensional simulations, their vertical structures in the non-linear regime are somewhat unexplored. This neglect is unwarranted given that any strong vertical motions associated with these waves could profoundly impact dust dynamics, dust sedimentation, planet formation, and the emissivity of the disc surface. In this paper we combine linear calculations and shearing box simulations in order to investigate the vertical structure of spiral waves for various polytropic stratifications and wave amplitudes. For sub-adiabatic profiles we find that spiral waves develop a pair of counter-rotating poloidal rolls. Particularly strong in the nonlinear regime, these vortical structures issue from the baroclinicity supported by the background vertical entropy gradient. They are also intimately connected to the disk's g-modes which appear to interact nonlinearly with the density waves. Furthermore, we demonstrate that the poloidal rolls are ubiquitous in gravitoturbulence, emerging in the vicinity of GI spiral wakes, and potentially transporting grains off the disk midplane. Other than hindering sedimentation and planet formation, this phenomena may bear on observations of the disk's scattered infrared luminosity. The vortical features could also impact on the turbulent dynamo operating in young protoplanetary discs subject to GI, or possibly even galactic discs.

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.

If sufficiently irradiated by its central star, a protoplanetary disks falls into an equilibrium state exhibiting vertical shear. This state may be subject to a hydrodynamical instability, the `vertical shear instability' (VSI), whose breakdown into turbulence transports a moderate amount of angular momentum while also facilitating planet formation, possibly via the production of small-scale vortices. In this paper, we show that VSI modes (a) exhibit arbitrary spatial profiles and (b) remain nonlinear solutions to the incompressible local equations, no matter their amplitude. The modes are themselves subject to parasitic Kelvin-Helmholtz instability, though the disk rotation significantly impedes the parasites and permits the VSI to attain large amplitudes (fluid velocities < 10% the sound speed). This `delay' in saturation probably explains the prominence of the VSI linear modes in global simulations. More generally, the parasites may set the amplitude of VSI turbulence in strongly irradiated disks. They are also important in breaking the axisymmetry of the flow, via the unavoidable formation of vortices. The vortices, however, are not aligned with the orbital plane and thus express a pronounced $z$-dependence. We also briefly demonstrate that the vertical shear has little effect on the magnetorotational instability, whereas magnetic fields easily quench the VSI, a potential issue in the ionised surface regions of the disk and also at larger radii.

Though usually treated in isolation, the magnetorotational and gravitational instabilities (MRI and GI) may coincide at certain radii and evolutionary stages of protoplanetary discs and active galactic nuclei. Their mutual interactions could profoundly influence several important processes, such as accretion variability and outbursts, fragmentation and disc truncation, or large-scale magnetic field production. Direct numerical simulations of both instabilities are computationally challenging and remain relatively unexplored. In this paper, we aim to redress this neglect via a set of 3D vertically stratified shearing-box simulations, combining self-gravity and magnetic fields. We show that gravito-turbulence greatly weakens the zero-net-flux MRI. In the limit of efficient cooling (and thus enhanced GI), the MRI is completely suppressed, and yet strong magnetic fields are sustained by the gravitoturbulence. This turbulent `spiral wave' dynamo may have widespread application, especially in galactic discs. Finally, we present preliminary work showing that a strong net-vertical-flux revives the MRI and supports a magnetically dominated state, in which the GI is secondary.

Local models of gaseous accretion discs have been successfully employed for decades to describe an assortment of small-scale phenomena, from instabilities and turbulence, to dust dynamics and planet formation. For the most part, they have been derived in a physically motivated but essentially ad hoc fashion, with some of the mathematical assumptions never made explicit nor checked for consistency. This approach is susceptible to error, and it is easy to derive local models that support spurious instabilities or fail to conserve key quantities. In this paper we present rigorous derivations, based on an asympototic ordering, and formulate a hierarchy of local models (incompressible, Boussinesq, and compressible), making clear which is best suited for a particular flow or phenomenon while spelling out explicitly the assumptions and approximations of each. We also discuss the merits of the anelastic approximation, emphasising that anelastic systems struggle to conserve energy unless strong restrictions are imposed on the flow. The problems encountered by the anelastic approximation are exacerbated by the disk's differential rotation, but also attend non-rotating systems such as stellar interiors. We conclude with a defence of local models and their continued utility in astrophysical research.

Young protoplanetary discs and the outer radii of active galactic nucleii may be subject to gravitational instability and, as a consequence, fall into a `gravitoturbulent' state. While in this state, appreciable angular momentum can be transported. Alternatively, the gas may collapse into bound clumps, the progenitors of planets or stars. In this paper, we numerically characterize the properties of 3D gravitoturbulence, focussing especially on its dependence on numerical parameters (resolution, domain size) and its excitation of small-scale dynamics. Via a survey of vertically stratified shearing box simulations with PLUTO and RODEO, we find (a) evidence that certain gravitoturbulent properties are independent of horizontal box size only when the box is larger than $\simeq 40 H_0$, where $H_0$ is the height scale, (b) at high resolution, small-scale isotropic turbulence appears off the midplane around $z\simeq 0.5 -1 H_0$, and (c) this small-scale dynamics results from a parametric instability, involving the coupling of inertial waves with a large-scale axisymmetric epicyclic mode. This mode oscillates at a frequency close to $\Omega$ and is naturally excited by gravito-turbulence, via a nonlinear process to be determined. The small-scale turbulence we uncover has potential implications for a wide range of disc physics, e.g. turbulent saturation levels, fragmentation, turbulent mixing, and dust settling.

Z Cam variables are a subclass of dwarf nova that lie near a global bifurcation between outbursting ("limit cycle") and non-outbursting ("standstill") states. It is believed that variations in the secondary star's mass-injection rate instigate transitions between the two regimes. In this paper we explore an alternative trigger for these transitions: stochastic fluctuations in the disk's turbulent viscosity. We employ simple one-zone and global viscous models which, though inappropriate for detailed matching to observed light curves, clearly indicate that turbulent disk fluctuations induce outbursts when the system is sufficiently close to the global bifurcation point. While the models easily produce the observed "outburst dip" pairs exhibited by Z Cam and Nova-like variables, they struggle to generate long trains of outbursts. We conclude that mass transfer variability is the dominant physical process determining the overall Z Cam standstill outburst pattern, but that viscous stochasticity provides an additional ingredient explaining some of the secondary features observed.

A long-standing puzzle in the study of black-hole accretion concerns the presence or not of thermal instability. Classical theory predicts the encircling accretion disk is unstable, as do some self-consistent MHD simulations of the flow. Yet observations of strongly accreting sources generally fail to exhibit cyclic or unstable dynamics on the expected timescales. This paper checks whether turbulent fluctuations impede thermal instability. It also asks if it makes sense to conduct linear stability analyses on a turbulent background. These issues are explored with a set of MRI simulations in thermally unstable local boxes in combination with stochastic equations that approximate the disk energetics. These models show that the disk's thermal behaviour deviates significantly from laminar theory, though ultimately a thermal runaway does occur. We find that the disk temperature evolves as a biased random walk, rather than increasing exponentially, and thus generates a broad spread of outcomes, with instability often delayed for several thermal times. We construct a statistical theory that describes some of this behaviour, emphasising the importance of the `escape time' and its associated probability distribution. In conclusion, turbulent fluctuations on there own cannot stabilise a disk, but they can weaken and delay thermal instability.

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.

In protoplanetary disks the aerodynamical friction between particles and gas induces a variety of instabilities that facilitate planet formation. Of these we examine the so-called `secular gravitational instability' (SGI) in the two-fluid approximation, deriving analytical expressions for its stability criteria and growth rates. Concurrently, we present a physical explanation of the instability that shows how it manifests upon an intermediate range of lengthscales exhibiting geostrophic balance in the gas component. In contrast to a single-fluid treatment, the SGI is quenched within a critical disk radius, as large as 10 AU and 30 AU for cm and mm sized particles respectively, although establishing robust estimates is hampered by uncertainties in the parameters (especially the strength of turbulence) and deficiencies in the razor-thin disk model we employ. It is unlikely, however, that the SGI is relevant for well-coupled dust. We conclude by applying these results to the question of planetesimal formation and the provenance of large-scale dust rings.

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.

Debris discs are commonly swathed in gas which can be observed in UV, in fine structure lines in FIR, and in resolved maps of CO emission. Carbon and oxygen are overabundant in such gas, but it is severely depleted in hydrogen. As a consequence, its ionisation fraction is remarkably high, suggesting magnetohydrodynamic (MHD) processes may be important. In particular, the gas may be subject to the magnetorotational instability (MRI), and indeed recent modelling of $\beta$ Pictoris requires an anomalous viscosity to explain the gas's observed radial structure. In this paper we explore the possibility that the MRI is active in debris-disc gas and responsible for the observed mass transport. We find that non-ideal MHD and dust-gas interactions play a subdominant role, and that linear instability is viable at certain radii. However, owing to low gas densities, the outer parts of the disc could be stabilised by a weak ambient magnetic field, though it is difficult to constrain such a field. Even if the MRI is stabilised by too strong a field, a magnetocentrifugal wind may be launched in its place and this could lead to equivalent (non-turbulent) transport. Numerical simulations of the vertically stratified MRI in conditions appropriate to the debris disc gas should be able to determine the nature of the characteristic behaviour at different radii, and decide on the importance of the MRI (and MHD more generally) on the evolution of these discs.

Recently it has been proposed that the zombie vortex instability (ZVI) could precipitate hydrodynamical activity and angular momentum transport in unmagnetised regions of protoplanetary discs, also known as "dead zones". In this letter we scrutinise, with high resolution 3D spectral simulations, the onset and survival of this instability in the presence of viscous and thermal physics. First, we find that the ZVI is strongly dependent on the nature of the viscous operator. Although the ZVI is easily obtained with hyper-diffusion, it is difficult to sustain with physical (second order) diffusion operators up to Reynolds numbers as high as 10^7. This sensitivity is probably due to the ZVI's reliance on critical layers, whose characteristic lengthscale, structure, and dynamics are controlled by viscous diffusion. Second, we observe that the ZVI is sensitive to radiative processes, and indeed only operates when the Peclet number is greater than a critical value ~10^4, or when the cooling time is longer than ~10 Omega^{-1}. As a consequence, the ZVI struggles to appear at R>0.3 AU in standard 0.01 M_sun T Tauri disc models, though younger more massive disks provide a more hospitable environment. Together these results question the prevalence of the ZVI in protoplanetary discs.

Gravitational instability (GI) features in several aspects of protostellar disk evolution, most notably in angular momentum transport, fragmentation, and the outbursts exemplified by FU Ori and EX Lupi systems. The outer regions of protostellar discs may also be coupled to magnetic fields, which could then modify the development of GI. To understand the basic elements of their interaction, we perform local 2D ideal and resistive MHD simulations with an imposed toroidal field. In the regime of moderate plasma beta, we find that the system supports a hot gravito-turbulent state, characterised by considerable magnetic energy and stress and a surprisingly large Toomre parameter $Q~10$. This result has potential implications for disk structure, vertical thickness, ionisation, etc. Our simulations also reveal the existence of long-lived and dense `magnetic islands' or plasmoids. Lastly, we find that the presence of a magnetic field has little impact on the fragmentation criterion of the disk. Though our focus is on protostellar disks, some of our results may be relevant for the outer radii of AGN.

This paper explores the driving of low-level hydrodynamical activity in protoplanetary-disc dead zones. A small adverse radial entropy gradient, ordinarily stabilised by rotation, excites oscillatory convection (`convective overstability') when thermal diffusion, or cooling, is neither too strong nor too weak. I revisit the linear theory of the instability, discuss its prevalence in protoplanetary discs, and show that unstable modes are exact nonlinear solutions in the local Boussinesq limit. Overstable modes cannot grow indefinitely, however, as they are subject to a secondary parametric instability that limits their amplitudes to relatively low levels. If parasites set the saturation level of the ensuing turbulence then the convective overstability is probably too weak to drive significant angular momentum transport or to generate vortices. But I also discuss an alternative, and far more vigorous, saturation route that generates radial `layers' or `zonal flows' (witnessed also in semiconvection). Numerical simulations are required to determine which outcome is favoured in realistic discs, and consequently how important the instability is for disc dynamics.

We determine how MRI-turbulent stresses depend on gas pressure via a suite of unstratified shearing box simulations. Earlier numerical work reported only a very weak dependence at best, results that call into question the canonical alpha-disk model and the thermal stability results that follow from it. Our simulations, in contrast, exhibit a stronger relationship, and show that previous work was box-size limited: turbulent `eddies' were artificially restricted by the numerical domain rather than by the scale height. Zero-net-flux runs without physical diffusion coefficients yield a stress proportional to $P^{0.5}$, where P is pressure. The stresses are also proportional to the grid length and hence remain numerically unconverged. The same runs with physical diffusivities, however, give a result closer to an alpha-disk: the stress is proportional to $P^{0.9}$. Net-flux simulations without explicit diffusion exhibit stresses proportional to $P^{0.5}$, but stronger imposed fields weaken this correlation. In summary, compressibility is important for the saturation of the MRI, but the exact stress-pressure relationship is difficult to ascertain in local simulations because of numerical convergence issues and the influence of any imposed flux. As a consequence, the interpretation of thermal stability behaviour in local simulations is a problematic enterprise.

We revisit the linear MRI in a cylindrical model of an accretion disk and uncover a number of attractive results overlooked in previous treatments. In particular, we elucidate the connection between local axisymmetric modes and global modes, and show that a local channel flow corresponds to the evanescent part of a global mode. In addition, we find that the global problem reproduces the local dispersion relation without approximation, a result that helps explain the success the local analysis enjoys in predicting global growth rates. MRI channel flows are nonlinear solutions to the governing equations in the local shearing box. However, only a small subset of MRI modes share the same property in global disk models, providing further evidence that the prominence of channels in local boxes is artificial. Finally, we verify our results via direct numerical simulations with the Godunov code RAMSES.

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

In protoplanetary disks, the inner boundary between the turbulent and laminar regions is a promising site for planet formation because solids may become trapped at the interface itself or in vortices generated by the Rossby wave instability. The disk thermodynamics and the turbulent dynamics at that location are entwined because of the importance of turbulent dissipation on thermal ionization and, conversely, of thermal ionisation on the turbulence. However, most previous work has neglected this dynamical coupling and have thus missed a key element of the physics in this region. In this paper, we aim to determine how the the interplay between ionization and turbulence impacts on the formation and evolution of vortices at the interface between the active and the dead zones. Using the Godunov code RAMSES, we have performed a 3D magnetohydrodynamic global numerical simulation of a cylindrical model of an MRI--turbulent protoplanetary disk, including thermodynamical effects as well as a temperature-dependant resistivity. The comparison with an analogous 2D viscous simulation has been extensively used to help identify the relevant physical processes and the disk's long-term evolution. We find that a vortex formed at the interface, due to Rossby wave instability, migrates inward and penetrates the active zone where it is destroyed by turbulent motions. Subsequently, a new vortex emerges a few tens of orbits later at the interface, and the new vortex migrates inward too. The sequence repeats itself, resulting in cycles of vortex formation, migration, and disruption. This behavior is successfully reproduced using two different codes. In this paper, we characterize this vortex life cycle and discuss its implications for planet formation at the dead/active interface. Our simulations highlight the importance of thermodynamical processes for the vortex evolution at the dead zone inner edge.

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.

In protoplanetary disks, the inner boundary between the turbulent and laminar regions could be a promising site for planet formation, thanks to the trapping of solids at the boundary itself or in vortices generated by the Rossby wave instability. At the interface, the disk thermodynamics and the turbulent dynamics are entwined because of the importance of turbulent dissipation and thermal ionization. Numerical models of the boundary, however, have neglected the thermodynamics, and thus miss a part of the physics. The aim of this paper is to numerically investigate the interplay between thermodynamics and dynamics in the inner regions of protoplanetary disks by properly accounting for turbulent heating and the dependence of the resistivity on the local temperature. Using the Godunov code RAMSES, we performed a series of 3D global numerical simulations of protoplanetary disks in the cylindrical limit, including turbulent heating and a simple prescription for radiative cooling. We find that waves excited by the turbulence significantly heat the dead zone, and we subsequently provide a simple theoretical framework for estimating the wave heating and consequent temperature profile. In addition, our simulations reveal that the dead-zone inner edge can propagate outward into the dead zone, before staling at a critical radius that can be estimated from a mean-field model. The engine driving the propagation is in fact density wave heating close to the interface. A pressure maximum appears at the interface in all simulations, and we note the emergence of the Rossby wave instability in simulations with extended azimuth. Our simulations illustrate the complex interplay between thermodynamics and turbulent dynamics in the inner regions of protoplanetary disks. They also reveal how important activity at the dead-zone interface can be for the dead-zone thermodynamic structure.

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.

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.

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 present results from large-scale particle simulations of the viscous overstability in Saturn's rings. The overstability generates a variety of structure on scales covering a few hundred metres to several kilometres, including axisymmetric wavetrains and larger-scale modulations. Such patterns have been observed in Saturn's rings by the Cassini spacecraft. Our simulations model the collisional evolution of particles in a co-rotating patch of the disk. These are the largest N-body simulations of the viscous overstability yet performed. The radial box size is five orders of magnitude larger than a typical particle radius, and so describes a 20-50 km radial portion of the rings. Its evolution is tracked for more than 10,000 orbits. In agreement with hydrodynamics, our N-body simulations reveal that the viscous overstability exhibits a rich set of dynamics characterised by nonlinear travelling waves with wavelengths of a few hundred meters. In addition, wave defects, such as sources and shocks, punctuate this bed of waves and break them up into large-scale divisions of radial width ~5 km. We find that the wavelength of the travelling waves is positively correlated with the mean optical depth. In order to assess the role of the numerical boundary conditions and also background ring structure, we include simulations of broad spreading rings and simulations with a gradient in the background surface density. Overall, our numerical results and approach provide a tool with which to interpret Cassini occultation observations of microstructure in Saturn's rings. We present an example of such a synthetic occultation observation and discuss what features to expect. We also make the entire source code freely available.

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 protoplanetary disks, the inner radial boundary between the MRI turbulent (`active') and MRI quiescent (`dead') zones plays an important role in models of the disk evolution and in some planet formation scenarios. In reality, this boundary is not well-defined: thermal heating from the star in a passive disk yields a transition radius close to the star (<0.1 au), whereas if the disk is already MRI active, it can self-consistently maintain the requisite temperatures out to a transition radius of roughly 1 au. Moreover, the interface may not be static; it may be highly fluctuating or else unstable. In this paper, we study a reduced model of the dynamics of the active/dead zone interface that mimics several important aspects of a real disk system. We find that MRI-transition fronts propagate inward (a `dead front' suppressing the MRI) if they are initially at the larger transition radius, or propagate outward (an `active front' igniting the MRI) if starting from the smaller transition radius. In both cases, the front stalls at a well-defined intermediate radius, where it remains in a quasi-static equilibrium. We propose that it is this new, intermediate stalling radius that functions as the true boundary between the active and dead zones in protoplanetary disks. These dynamics are likely implicated in observations of variable accretion, such as FU Ori outbursts, as well as in those planet formation theories that require the accumulation of solid material at the dead/active interface.

The recurrent outbursts that characterise low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionises near 5000 K. This physics leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low accreting states and a family of hot and efficiently accreting states. Previous models have parametrised the disc turbulence via an alpha (or `eddy') viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other on orbital times; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full MHD turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries, and possibly young stellar objects.

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.

In this paper we investigate how convective instabilities influence heat conduction in the intracluster medium (ICM) of cool-core galaxy clusters. The ICM is a high-beta, weakly collisional plasma in which the transport of momentum and heat is aligned with the magnetic field. The anisotropy of heat conduction, in particular, gives rise to instabilities that can access energy stored in a temperature gradient of either sign. We focus on the heat-flux buoyancy-driven instability (HBI), which feeds on the outwardly increasing temperature profile of cluster cool cores. Our aim is to elucidate how the global structure of a cluster impacts on the growth and morphology of the linear HBI modes when in the presence of Braginskii viscosity, and ultimately on the ability of the HBI to thermally insulate cores. We employ an idealised quasi-global model, the plane-parallel atmosphere, which captures the essential physics -- e.g. the global radial profile of the cluster -- while letting the problem remain analytically tractable. Our main result is that the dominant HBI modes are localised to the the innermost (~<20%) regions of cool cores. It is then probable that, in the nonlinear regime, appreciable field-line insulation will be similarly localised. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, heat conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Finally, our linear solutions provide a convenient numerical test for the nonlinear codes that tackle the saturation of such convective instabilities in the presence of anisotropic transport.

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.

The isorotation contours of the solar convective zone (SCZ) show three distinct morphologies, corresponding to two boundary layers (inner and outer), and the bulk of the interior. Previous work has shown that the thermal wind equation together with informal arguments on the nature of convection in a rotating fluid could be used to deduce the shape of the isorotation surfaces in the bulk of the SCZ with great fidelity, and that the tachocline contours could also be described by relatively simple phenomenology. In this paper, we show that the form of these surfaces can be understood more broadly as a mathematical consequence of the thermal wind equation and a narrow convective shell. The analysis does not yield the angular velocity function directly, an additional surface boundary condition is required. But much can already be deduced without constructing the entire rotation profile. The mathematics may be combined with dynamical arguments put forth in previous works to the mutual benefit of each. An important element of our approach is to regard the constant angular velocity surfaces as an independent coordinate variable for what is termed the "residual entropy," a quantity that plays a key role in the equation of thermal wind balance. The difference between the dynamics of the bulk of the SCZ and the tachocline is due to a different functional form of the residual entropy in each region. We develop a unified theory for the rotational behavior of both the SCZ and the tachocline, using the solutions for the characteristics of the thermal wind equation. These characteristics are identical to the isorotation contours in the bulk of the SCZ, but the two deviate in the tachocline. The outer layer may be treated, at least descriptively, by similar mathematical techniques, but this region probably does not obey thermal wind balance.

We revisit, via a very simplified set of equations, a linear streaming instability (technically an overstability), which is present in, and potentially important for, dusty protoplanetary disks (Youdin & Goodman 2005). The goal is a better understanding of the physical origin of such instabilities, which are notoriously subtle. Rotational dynamics seem to be essential to this type of instability, which cannot be captured by one-dimensional Cartesian models. Dust `pileups' in moving pressure maxima are an important triggering mechanism of the instability, and drag feedback of dust upon the gas allows these maxima to be strengthened. Coriolis forces and the background drift counteract the effects of the pressure force.

This paper concerns kinematic helical dynamos in a spherical fluid body surrounded by an insulator. In particular, we examine their behaviour in the regime of large magnetic Reynolds number $\Rm$, for which dynamo action is usually concentrated upon a simple resonant stream-surface. The dynamo eigensolutions are computed numerically for two representative single-roll flows using a compact spherical harmonic decomposition and fourth-order finite-differences in radius. These solutions are then compared with the growth rates and eigenfunctions of the Gilbert and Ponty (2000) large $\Rm$ asymptotic theory. We find good agreement between the growth rates when $\Rm>10^4$, and between the eigenfunctions when $\Rm>10^5$.

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.

In this paper, we present a model for the effects of the tachocline on the differential rotation in the solar convection zone. The mathematical technique relies on the assumption that entropy is nearly constant ("well-mixed") in isorotation surfaces both outside and within the tachocline. The resulting solutions exhibit nontrivial features that strikingly resemble the true tachocline isorotation contours in unexpected detail. This strengthens the mathematical premises of the theory. The observed rotation pattern in the tachocline shows strong quadrupolar structure, an important feature that is explicitly used in constructing our solutions. The tachocline is treated locally as an interior boundary layer of small but finite thickness, and an explicit global solution is then constructed. A dynamical link can thus be established between the internal jump in the angular velocity at the tachocline and the spread of angular velocities observed near the solar surface. In general, our results suggest that the bulk of the solar convection zone is in thermal wind balance, and that simple quadrupolar stresses, local in radius, mediate the tachocline transition from differential rotation to uniform rotation in the radiative interior.

Simulations of the magnetorotational instability (MRI) in 'unstratified' shearing boxes exhibit powerful coherent flows, whereby the fluid vertically splits into countermoving planar jets or `channels'. Channel flows correspond to certain axisymmetric linear MRI modes, and their preponderance follows from the remarkable fact that they are approximate nonlinear solutions of the MHD equations in the limit of weak magnetic fields. We show in this paper, analytically and with one-dimensional numerical simulations, that this property is also shared by certain axisymmetric MRI modes in vertically-stratified shearing boxes. These channel flows rapidly capture significant amounts of magnetic and kinetic energy, and thus are vulnerable to secondary shear instabilities. We examine these parasites in the vertically stratified context, and estimate the maximum amplitudes that channels attain before they are destroyed. These estimates suggest that a dominant channel flow will usually drive the disk's magnetic field to thermal strengths. The prominence of these flows and their destruction place enormous demands on simulations, but channels in their initial stages also offer a useful check on numerical codes. These benchmarks are especially valuable given the increasing interest in the saturation of the stratified MRI. Lastly we speculate on the potential connection between 'run-away' channel flows and outburst behaviour in protostellar and dwarf nova disks.

We outline a novel linear instability that may arise in the dead-zones of protostellar disks, and possibly the fluid interiors of planets and protoplanets. In essence it is an axisymmetric buoyancy instability, but one that would not be present in a purely hydrodynamical gas. The necessary ingredients for growth include a negative radial entropy gradient (of any magnitude), weak magnetic fields, and efficient resistive diffusion (in comparison with thermal diffusion). The character of the instability is local, axisymmetric, and double-diffusive, and it attacks lengths much shorter than the resistive scale. Like the axisymmetric convective instability, it draws its energy from the negative radial entropy gradient; but by utilising the diffusing magnetic field, it can negate the stabilising influence of rotation. Its nonlinear saturated state, while not transporting appreciable angular momentum, could drive radial and vertical mixing, which may influence the temperature structure of the disk, dust dynamics and, potentially, planet formation.

This paper concerns the interaction between non-axisymmetric inertial waves and their corotation resonances in a hydrodynamical disk. Inertial waves are of interest because they can localise in resonant cavities circumscribed by Lindblad radii, and as a consequence exhibit discrete oscillation frequencies that may be observed. It is often hypothesised that these trapped eigenmodes are affiliated with the poorly understood QPO phenomenon. We demonstrate that a large class of non-axisymmetric 3D inertial waves cannot manifest as trapped normal modes. This class includes any inertial wave whose resonant cavity contains a corotation singularity. Instead, these `singular' modes constitute a continuous spectrum and, as an ensemble, are convected with the flow, giving rise to shearing waves. Lastly, we present a simple demonstration of how the corotation singularity stabilizes three-dimensional perturbations in a slender torus.

We show that the differential rotation profile of the solar convection zone, apart from inner and outer boundary layers, can be reproduced with great accu- racy if the isorotation contours correspond to characteristics of the thermal wind equation. This requires that there be a formal quantitative relationship involving the entropy and the angular velocity. Earlier work has suggested that this could arise from magnetohydrodynamic stability constraints; here we argue that purely hydrodynamical processes could also lead to such a result. Of special importance to the hydrodynamical solution is the fact that the thermal wind equation is insensitive to radial entropy gradients. This allows a much more general class of solutions to fit the solar isorotation contours, beyond just those in which the entropy itself must be a function of the angular velocity. In particular, for this expanded class, the thermal wind solution of the solar rotation profile remains valid even when large radial entropy gradients are present. A clear and explicit example of this class of solution appears to be present in published numerical simulations of the solar convective zone. Though hydrodynamical in character, the theory is not sensitive to the presence of weak magnetic fields. Thus, the identification of solar isorotation contours with the characteristics of the thermal wind equation appears to be robust, accommodating, but by no means requiring, magnetic field dynamics.

We run mean-field shearing-box numerical simulations with a temperature-dependent resistivity and compare them to a reduced dynamical model. Our simulations reveal the co-existence of two quasi-steady states, a `quiet' state and an `active' turbulent state, confirming the predictions of the reduced model. The initial conditions determine on which state the simulation ultimately settles. The active state is strongly influenced by the geometry of the computational box and the thermal properties of the gas. Cubic domains support permanent channel flows, bar-shaped domains exhibit eruptive behaviour, and horizontal slabs give rise to infrequent channels. Meanwhile, longer cooling time-scales lead to higher saturation amplitudes.

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.

Local simulations of the magnetorotational instability (MRI) in accretion disks can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin-Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are the clean isolated channels which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures inevitably achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.

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.

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.

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.