(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause converging flows which lead to strong spiral shocks once each superhump cycle. The dissipation associated with these shocks powers the superhump. We compare 2D and 3D results, and conclude that 3D simulations are necessary to faithfully simulate the disc dynamics. We ran our simulations for unprecedented durations, so that an eccentric equilibrium is established except at high mass ratios where the growth rate of the instability is very low. Our improved simulations give a closer match to the observed relationship between superhump period excess and binary mass ratio than previous numerical work. The observed black hole X-ray transient superhumpers appear to have systematically lower disc precession rates than the cataclysmic variables. This could be due to higher disc temperatures and thicknesses. The modulation in total viscous dissipation on the superhump period is overwhelmingly from the region of the disc within the 3:1 resonance radius. As the eccentric instability develops, the viscous torques are enhanced, and the disc consequently adjusts to a new equilibrium state, as suggested in the thermal-tidal instability model. We quantify this enhancement in the viscosity, which is ~10 per cent for q=0.08. We characterise the eccentricity distributions in our accretion discs, and show that the entire body of the disc partakes in the eccentricity.

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios q < 0.1, applicable to ultra-compact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For q ~ 0.1, the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q. At q = 0.01, the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com (Morales-Rueda et al. 2003), which has q = 0.02. The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultra-compact binaries.

We show that the combination of a weak magnetic propeller and accretion disc resonances can effectively halt accretion in short period cataclysmic variables for large fractions of their lifetimes. This may help to explain the discrepancy between the observed and predicted orbital period distributions of cataclysmic variables at short periods. Orbital resonances cause the disc to become eccentric, allowing material to fall back onto the donor star or out of the system. A weak magnetic field on a rapidly spinning primary star propels disc material outwards, allowing it to access these resonances. Numerical and analytic calculations show that this state can be long lived (~10^11 yr). This is because the magnetic propeller is required only to maintain access to the resonances, and not to push matter out of the Roche lobe, so that the spin down time-scale is much longer than for a classical propeller model.

The galactic microquasar GRS 1915+105 has been in a continuous state of outburst since 1992, over 20 times longer than any other black hole X-ray transient. Assuming that the outburst is powered via accretion of an irradiated gaseous disc, we calculate how the predicted outburst duration varies according to the efficiency of the self-irradiation mechanism. At least one current model leads to the conclusion that the end of the outburst is imminent. The timing of the decline of GRS 1915+105, whenever it arrives, will be an excellent discriminator of the self-irradiation mechanism in X-ray transients, allowing us to infer the fraction of the disc that is heated by the incident X-rays and the magnitude of the mass loss rate in the form of a wind.

We present Monte Carlo scattered light models of a warped disk that reproduce the observed photopolarimetric variability of the classical T Tauri star AA Tau. For a system inclination of 75 deg and using an analytic description for a warped inner disk, we find that the shape and amplitude of the photopolarimetric variability are reproduced with a warp which occults the star, located at 0.07 au, amplitude 0.016 au, extending over radial and azimuthal ranges 0.0084 au and 145 deg. We also show a time sequence of high spatial resolution scattered light images, showing a dark shadow cast by the warp sweeping round the disk. Using a modified smooth particle hydrodynamics code, we find that a stellar dipole magnetic field of strength 5.2 kG, inclined at 30 deg to the stellar rotation axis can reproduce the required disk warping to explain AA Tau's photopolarimetric variability.

The 2001 outburst of WZ Sagittae has shown the most compelling evidence yet for an enhancement of the mass transfer rate from the donor star during a dwarf nova outburst in the form of hot-spot brightening. I show that even in this extreme case, the brightening can be attributed to tidal heating near the interaction point of an accretion stream with the expanding edge of an eccentric accretion disc, with no need at all for an increase in the mass transfer rate. Furthermore, I confirm previous suggestions that an increase in mass transfer rate through the stream damps any eccentricity in an accretion disc and suppresses the appearance of superhumps, in contradiction to observations. Tidal heating is expected to be most significant in systems with small mass ratios. It follows that systems like WZ Sagittae - which has a tiny mass ratio - are those most likely to show a brightening in the hot-spot region.

We develop a new analytic solution for the steady-state structure of a thin accretion disc under the influence of a magnetic field that is anchored to the central star. The solution takes a form similar to that of Shakura and Sunyaev and tends to their solution as the magnetic moment of the star tends to zero. As well as the Kramer's law case, we obtain a solution for a general opacity. The effects of varying the mass transfer rate, spin period and magnetic field of the star as well as the opacity model applied to the disc are explored for a range of objects. The solution depends on the position of the magnetic truncation radius. We propose a new approach for the identification of the truncation radius and present an analytic expression for its position.

The outbursts of the dwarf nova WZ Sagittae occur on timescales which are much longer than expected from the standard thermal viscous disc instability model, which correctly reproduces the outburst behaviour of many other dwarf novae. In addition short `echo' outbursts have sometimes been observed after the main outburst. The increased timescales could be explained by a magnetic propeller, anchored to the primary star, which clears the inner disc . We argue that the propeller, combined with the resonances of the binary orbit, may also be the cause of the echo outbursts. We present numerical models which reproduce the echo outbursts.

We present very high resolution hydrodynamical simulations of accretion discs in black hole X-ray binaries accreting near the Eddington limit. The results show that mass loss, irradiation and tidal interactions all have a profound effect on the observed behaviour of long period X-ray transients. In particular, the interplay of all of these effects in the outer regions of the accretion disc is able to drive long time-scale (weeks to years) variability is these objects, and is a possible origin for some of the extreme variability of GRS1915+105.

In the proceedings of this, and of several recent close binary conferences, there have been several contributions describing smoothed particle hydrodynamics simulations of accretion disks. It is apposite therefore to review the numerical scheme itself with emphasis on its advantages for disk modelling, and the methods used for modelling viscous processes.

We propose a simple explanation for the constant mean brightness observed between outbursts in dwarf novae. Secular changes in the total energy dissipation rate of the accretion disc brought about by variations in surface density, temperature and disc radius can be regulated by the gradual cooling of a small, critically-stable hot inner region. The hypothesis is supported by two-dimensional time-dependent numerical models of dwarf nova accretion discs.

We present a model of an outburst of the soft X-ray transient A0620-003. A two-dimensional time-dependent smoothed particle hydrodynamics scheme is used to simulate the evolution of the accretion disc through a complete outburst. The scheme includes the full tidal potential of the binary and a simple treatment of the thermal-viscous disc instability. In the case where the mass accretion rate onto the primary determines the fraction of the disc that can be kept in a hot, high viscosity state by the resulting X-ray emission, we find that the shape of the X-ray light curve is ultimately determined by the relative sizes of the irradiated and unirradiated parts of the disc and the growth time-scale of the tidal instability. The model accounts for the rebrightening that has been observed in the light curves of A0620-003 and several other transients. The primary maximum and subsequent decline are due to the accretion of gas within the irradiated portion of the disc, while the secondary maximum is caused by the accretion of gas in the outer part of the disc that is initially shadowed from the central X-rays, but subject to tidal forces. We propose that tidal effects at the disc edge can be sufficient to drive accretion on a time-scale shorter than that expected for a standard alpha-viscosity disc. The final decay is subsequently controlled by the gradual retreat of the irradiated portion of the disc. If the entire disc is kept in the high-viscosity state by the irradiation, no rebrightening is possible.

We present the first detailed hydrodynamic simulation of a superoutburst to incorporate the full tidal potential of a binary system. A two-dimensional smoothed particle hydrodynamics code is used to simulate a superoutburst in a binary with the parameters of the SU UMa system Z Chamaeleontis. The simulated light curves shows all the features observed in such systems. Analysis of the mass flux through the disc and the growth rate of the superhumps and disc eccentricity show that the superoutburst-superhump phenomenon is a direct result of tidal instability. No enhanced mass transfer from the secondary is required to initiate or sustain these phenomena. Comparisons of superoutbursts with normal outbursts are made and we show that the model can be reconciled with the behavior of U Geminorum type dwarf novae, which show no superoutbursts.

We present a numerical scheme for the evolution of an accretion disc through a dwarf nova outburst. We introduce a time-varying artificial viscosity into an existing smoothed particle hydrodynamics code optimised for two and three-dimensional simulations of accretion discs. The technique gives rise to coherent outbursts and can easily be adapted to include a complete treatment of thermodynamics. We apply a two-dimensional isothermal scheme to the system SS Cygni and present a wide range of observationally testable results.