The Visible and Infrared Survey Telescope for Astronomy (VISTA) is the 4-metre wide-field survey telescope at ESO's Paranal Observatory, equipped with the world's largest near-infrared imaging camera (VISTA IR Camera, VIRCAM), with 1.65 degree diameter field of view, and 67 Mpixels giving 0.6 square degrees active pixel area, operating at wavelengths 0.8 - 2.3 microns. We provide a short history of the project, and an overview of the technical details of the full system including the optical design, mirrors, telescope structure, IR camera, active optics, enclosure and software. The system includes several innovative design features such as the f/1 primary mirror, the dichroic cold-baffle camera design and the sophisticated wavefront sensing system delivering closed-loop 5-axis alignment of the secondary mirror. We conclude with a summary of the delivered performance, and a short overview of the six ESO public surveys in progress on VISTA.

We present three dimensional smoothed particle hydrodynamics (SPH) calculations of irradiation-driven warping of accretion discs. Initially unwarped planar discs are unstable to the radiation reaction when the disc is illuminated by a central radiation source. The disc warps and tilts and precesses slowly in a retrograde direction; its shape continuously flexes in response to the changing orientation of the Roche potential. We simulate ten systems: eight X-ray binaries, one cataclysmic variable (CV), and a `generic' low mass X-ray binary (LMXB). We adopt system parameters from observations and tune a single parameter: our model X-ray luminosity ($L_{*}$) to reproduce the observed or inferred super-orbital periods. Without exception, across a wide range of parameter space, we find an astonishingly good match between the observed $L_{X}$ and the model $L_{*}$. We conclude irradiation-driven warping is the mechanism underlying the long periods in X-ray binaries. Our Her X-1 simulation simultaneously reproduces the observed $L_{X}$, the "main-" and "short-high" X-ray states and the orbital inclination. Our simulations of SS 433 give a maximum warp angle of $18.6^{\circ}$, a good match to the cone traced by the jets, but this angle is reached only in the outer disc. In all cases, the overall disc tilt is less than $\degrees{13}$ and the maximum disc warp is less than and or equal to $\degrees{21}$.

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad-hoc manner. Thus no complete and direct comparison with observations of the pulsar population within the Galactic disk has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed \textit{typical} MSP formation process. Results also focus on reproducing the observed $P\dot{P}$ diagram for Galactic pulsars and how this precludes short timescales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.

The majority of pulsar population synthesis studies performed to date have focused on isolated pulsar evolution. Those that have incorporated pulsar evolution within binary systems have tended to either treat binary evolution poorly of evolve the pulsar population in an ad-hoc manner. Here we present the first model of the Galactic field pulsar population that includes a comprehensive treatment of both binary and pulsar evolution. Synthetic observational surveys mimicking a variety of radio telescopes are then performed on this population. As such, a complete and direct comparison of model data with observations of the pulsar population within the Galactic disk is now possible. The tool used for completing this work is a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Here we give a brief overview of the method and assumptions involved with each component. Some preliminary results are also presented as well as plans for future applications of the code.

Aims: We investigate the behaviour of dust in protoplanetary disks under the action of gas drag in the presence of a planet. Our goal is twofold: to determine the spatial distribution of dust depending on grain size and planet mass, and therefore to provide a framework for interpretation of coming observations and future studies of planetesimal growth. Method: We numerically model the evolution of dust in a protoplanetary disk using a two-fluid (gas + dust) Smoothed Particle Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the evolution of a minimum mass solar nebula (MMSN) disk comprising 1% dust by mass in the presence of an embedded planet. We run a series of simulations which vary the grain size and planetary mass to see how they affect the resulting disk structure. Results: We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system with a given stellar, disk and planetary mass will have a completely different appearance depending on the grain size. For low mass planets in our MMSN disk, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates at the external edge of the planetary gap and speculate that the presence of a planet in the disk may enhance the formation of a second planet by facilitating the growth of planetesimals in this high density region.

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

Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems.

We present three dimensional smoothed particle hydrodynamics (SPH) calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and in general the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios q ~ 0.1 a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (c.f. X 1916-053). For less extreme mass ratios the entire disc tilts out of the orbital plane (c.f. Her X-1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction we obtain both positive and negative superhumps simultaneously in the dissipation light curve (c.f. V603 Aql).

We compare the accretion flow onto the neutron star induced by Roche lobe overflow with that by the overflow from the Be disk, in a zero eccentricity, short period binary with the same mass transfer rate, performing three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a persistent accretion disk is formed around the neutron star in both cases. The circularization radius of the material transferred via Roche lobe overflow is larger than that of the material transfered from the Be disk. Thus, the growth of the accretion disk in the former case becomes significantly slower than in the latter case. In both cases, the mass accretion rate is very small and varies little with orbital phase, which is consistent with the observed X-ray behaviour of Be/X-ray binaries with circular orbits (e.g. XTE J1543-568).

We present results from three XMM-Newton observations of the M31 low mass X-ray binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-hr period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anti-correlated with source intensity. However, our new observations do not favour a strict intensity dependance, but rather suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio <~ 0.3 (period <~ 4 hr), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of ~11 P_orb, and prograde disc precession on a period of ~29 P_orb. This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on a 81+/-3 hr cycle.

The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that the dip depth was anticorrelated with the X-ray luminosity. However, we present results from three new XMM-Newton observations that suggest that the evolution of dipping is instead due to precession of the accretion disc. Such precession is expected in neutron star LMXBs with mass ratios <0.3 (i.e. with orbital periods <4 hr), such as the Galactic dipping LMXB 4U 1916-053. We simulated the accretion disc of Bo 158 using cutting-edge 3D smoothed particle hydrodynamics (SPH), and using the observed parameters. Our results show disc variability on two time-scales. The disc precesses in a prograde direction on a period of 81+/-3 hr. Also, a radiatively-driven disc warp is present in the inner disc, which undergoes retrograde precesson on a \~31 hr period. From the system geometry, we conclude that the dipping evolution is driven by the disc precession. Hence we predict that the dipping behaviour repeats on a ~81 hr cycle.

We present the results of a three dimensional, locally isothermal, non-self-gravitating SPH code which models protoplanetary disks with two fluids: gas and dust. We ran simulations of a 1 Msun star surrounded by a 0.01 Msun disk comprising 99% gas and 1% dust in mass and extending from 0.5 to ~300 AU. The grain size ranges from 0.001 mm to 10 m for the low resolution (~25 000 SPH particles) simulations and from 0.1 mm to 10 cm for the high resolution (~160 000 SPH particles) simulations. Dust grains are slowed down by the sub-Keplerian gas and lose angular momentum, forcing them to migrate towards the central star and settle to the midplane. The gas drag efficiency varies according to the grain size, with the larger bodies being weakly influenced and following marginally perturbed Keplerian orbits, while smaller grains are strongly coupled to the gas. For intermediate sized grains, the drag force decouples the dust and gas, allowing the dust to preferentially migrate radially and efficiently settle to the midplane. The resulting dust distributions for each grain size will indicate, when grain growth is added, the regions when planets are likely to form.

We present seven nights' blue spectroscopy of the novalike variable star, V3885 Sgr. The line spectrum shows a typical combination of broad absorption and emission in H gamma and H beta and He I, which is associated with the accretion disk. We also observe anti-phased narrow emission, which we attribute to irradiation of the secondary star. The He II 4686 and N III--C III--C IV emission lines are devoid of structure and are most likely formed in an outflow. We measure radial velocity shifts in the absorption and emission lines, from which we fit an orbital period of 4.97126+/-0.00036h. From the velocity semi-amplitudes of the disk and companion star, we are able to constrain the binary mass ratio to q>0.7. The phase-folded spectra provide dense coverage of the entire orbital cycle. Doppler tomograms of the hydrogen and He I lines reveal spiral structure in the accretion disk and the irradiated donor star. We believe this is the first unambiguous detection of spiral waves in a novalike variable.

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.

Two dimensional SPH simulations of a precessing accretion disc in a q=0.1 binary system (such as XTE J1118+480) reveal complex and continuously varying shape, kinematics, and dissipation. The stream-disc impact region and disc spiral density waves are prominent sources of energy dissipation.The dissipated energy is modulated on the period P_{sh} = ({P_{orb}}^{-1}-{P_{prec}}^{-1}^{-1} with which the orientation of the disc relative to the mass donor repeats. This superhump modulation in dissipation energy has a variation in amplitude of ~10% relative to the total dissipation energy and evolves, repeating exactly only after a full disc precession cycle. A sharp component in the light curve is associated with centrifugally expelled material falling back and impacting the disc. Synthetic trailed spectrograms reveal two distinct "S-wave" features, produced respectively by the stream gas and the disc gas at the stream-disc impact shock. These S-waves are non-sinusoidal, and evolve with disc precession phase. We identify the spiral density wave emission in the trailed spectrogram. Instantaneous Doppler maps show how the stream impact moves in velocity space during an orbit. In our maximum entropy Doppler tomogram the stream impact region emission is distorted, and the spiral density wave emission is uppressed. A significant radial velocity modulation of the whole line profile occurs on the disc precession period. We compare our SPH simulation with a simple 3D model: the former is appropriate for comparison with emission lines while the latter is preferable for skewed absorption lines from precessing discs.

We present the results of three-dimensional numerical simulations that include the effects of hydrodynamical forces and gas drag upon an evolving dusty gas disk. We briefly describe a new parallel, two phase numerical code based upon the smoothed particle hydrodynamics (SPH) technique in which the gas and dust phases are represented by two distinct types of particles. We use the code to follow the dynamical evolution of a population of grains in a gaseous protoplanetary disk in order to understand the distribution of grains of different sizes within the disk. Our ``grains'' range from metre to submillimetre in size.

We have developed a new numerical technique for simulating dusty-gas flows. Our unique code incorporates gas hydrodynamics, self-gravity and dust drag to follow the dynamical evolution of a dusty-gas medium. We have incorporated several descriptions for the drag between gas and dust phases and can model flows with submillimetre, centimetre and metre size "dust". We present calculations run on the APAC* supercomputer following the evolution of the dust distribution in the pre-solar nebula. *Australian Partnership for Advanced Computing

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 demonstrate that measurable vertical structure can be excited in the accretion disc of a close binary system by a dipolar magnetic field centred on the secondary star. We present the first high resolution hydrodynamic simulations to show the initial development of a uniform warp in a tidally truncated accretion disc. The warp precesses retrogradely with respect to the inertial frame. The amplitude depends on the phase of the warp with respect to the binary frame. A warped disc is the best available explanation for negative superhumps.

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 propose a mechanism for the superhump modulations observed in optical photometry of at least two black hole X-ray transients (SXTs). As in extreme mass-ratio cataclysmic variables (CVs), superhumps are assumed to result from the presence of the 3:1 orbital resonance in the accretion disc. This causes the disc to become non-axisymmetric and precess. However the mechanism for superhump luminosity variations in low mass X-ray binaries (LMXBs) must differ from that in CVs, where it is attributed to a tidally-driven modulation of the disc's viscous dissipation, varying on the beat between the orbital and disc precession period. By contrast in LMXBs, tidal dissipation in the outer accretion disc is negligible: the optical emission is overwhelming dominated by reprocessing of intercepted central X-rays. Thus a different origin for the superhump modulation is required. Recent observations and numerical simulations indicate that in an extreme mass-ratio system the disc area changes on the superhump period. We deduce that the superhumps observed in SXTs arise from a modulation of the reprocessed flux by the changing area. Therefore, unlike the situation in CVs, where the superhump amplitude is inclination-independent, superhumps should be best seen in low-inclination LMXBs, whereas an orbital modulation from the heated face of the secondary star should be more prominent at high inclinations. Modulation at the disc precession period (10s of days) may indicate disc asymmetries such as warping. We comment on the orbital period determinations of LMXBs, and the possibility and significance of possible permanent superhump LMXBs.

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.

We investigate the role of the eccentric disc resonance in systems with mass ratios q greater than 1/4, and demonstrate the effects that changes in the mass flux from the secondary star have upon the disc radius and structure. The addition of material with low specific angular momentum to its outer edge restricts a disc radially. Should the mass flux from the secondary be reduced, it is possible for the disc in a system with mass ratio as large as 1/3 to expand to the 3:1 eccentric inner Lindblad resonance and for superhumps to be excited.

We consider the precession rates of eccentric discs in close binaries, and compare theoretical predictions with the results of numerical disc simulations and with observed superhump periods. A simple dynamical model for precession is found to be inadequate. For mass ratios less than approximately 1/4 a linear dynamical model does provide an upper limit for disc precession rates. Theory suggests that pressure forces have a significant retrograde impact upon the precession rate (Lubow 1992). We find that the disc precession rates for three systems with accurately known mass ratios are significantly slower than predicted by the dynamical theory, and we attribute the difference to pressure forces. By assuming that pressure forces of similar magnitude occur in all superhumping systems, we obtain an improved fit to superhump observations.

We describe three presentations relating to the X-ray pulsar GX 1+4 at a workshop on magnetic fields and accretion at the Astrophysical Theory Centre, Australian National University on 1998, November 12-13. Optical and X-ray spectroscopy indicate that GX 1+4 is seen through a cloud of gravitationaly bound matter. We discuss an unstable negative feedback mechanism (originally proposed by Kotani et al, 1999), based on X-ray heating of this matter which controls the accretion rate when the source is in a low X-ray luminosity state. A deep minimum lasting ~6 hours occurred during observations with the RXTE satellite over 1996, July 19-21. The shape of the X-ray pulses changed remarkably from before to after the minimum. These changes may be related to the transition from neutron star spin-down to spin-up which occurred at about the same time. Smoothed particle hydrodynamic simulations of the effect of adding matter with opposite angular momentum to an existing disc, show that it is possible for a number of concentric rings with alternating senses of rotation to co-exist in a disc. This could provide an explanation for the step-like changes in Pdot which are observed in GX 1+4. Changes at the inner boundary of the disc occur at the same timescale as that imposed at the outer boundary. Reversals of material torque on the neutron star occur at a minimum in L_X.

The University of Durham Mark 6 Telescope has been used to make observations of a number of AGNs visible from the Southern Hemisphere. Results are presented for limits to VHE gamma ray emission from 1ES 0323+022, PKS 0829+046, 1ES 1101-232, Cen A, PKS 1514-24, RXJ 10578-275 and 1ES 2316-423, both for steady emission and for outbursts on timescales of ~ 1 day.

We numerically investigate the hydrodynamics of accretion disk reversal and relate our findings to the observed spin-rate changes in the accreting X-ray pulsar GX~1+4. In this system, which accretes from a slow wind, the accretion disk contains two dynamically distinct regions. In the inner part viscous forces are dominant and disk evolution occurs on a viscous timescale. In the outer part dynamical mixing of material with opposite angular momentum is more important, and the externally imposed angular momentum reversal timescale governs the flow. In this outer region the disk is split into concentric rings of material with opposite senses of rotation that do not mix completely but instead remain distinct, with a clear gap between them. We thus predict that torque reversals resulting from accretion disk reversals will be accompanied by minima in accretion luminosity.

We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.

We investigate the effects of tidally induced asymmetric disc structure on accretion onto the white dwarf in intermediate polars. Using numerical simulation, we show that it is possible for tidally induced spiral waves to propagate sufficiently far into the disc of an intermediate polar that accretion onto the central white dwarf could be modulated as a result. We suggest that accretion from the resulting asymmetric inner disc may contribute to the observed X-ray and optical periodicities in the light curves of these systems. In contrast to the stream-fed accretion model for these periodicities, the tidal picture predicts that modulation can exist even for systems with weaker magnetic fields where the magnetospheric radius is smaller than the radius of periastron of the mass transfer stream. We also predict that additional periodic components should exist in the emission from low mass ratio intermediate polars displaying superhumps.

We consider the implications of the detection of spiral structure in the accretion disc of the binary IP Pegasi. We use numerical simulations of the development of a disc outburst to construct predicted Doppler tomograms, which are found to be in close agreement with the observations if the spiral pattern arises as a transient feature when the disc expands viscously at the start of the outburst. The good agreement of such viscous disc simulations with the data is consistent with models in which most of the angular momentum transport in the disc originates in internal stresses rather than globally excited waves or shocks. Future detailed observations of the development of transient spiral features offer the potential to measure the dependence of the disc viscosity on the local physical conditions in the disc.

We numerically study the tidal instability of accretion discs in close binary systems using a two-dimensional SPH code. We find that the precession rate of tidally unstable, eccentric discs does not only depend upon the binary mass ratio q. Although the (prograde) disc precession rate increases with the strength of the tidal potential, it also increases with the shear viscosity. Increasing the disc temperature has a retrograde impact upon the precession rate. We find that motion relative to the binary potential results in superhump-like, periodic luminosity variations in the outer reaches of an eccentric disc. The nature and location of the luminosity modulation is a function of the shear viscosity, nu. Light curves most similar to observations are obtained for nu values appropriate for a dwarf nova in outburst. We investigate the thermal-tidal instability model for superoutburst. A dwarf nova outburst is simulated by instantaneously increasing nu, which causes a rapid radial expansion of the disc. Should the disc encounter the 3:1 eccentric inner Lindblad resonance and become tidally unstable, then tidal torques become much more efficient at removing angular momentum from the disc. The disc then shrinks and mass flux through the disc increases. The resulting increase in disc luminosity is found to be consistent with the excess luminosity of a superoutburst.

In this paper we briefly discuss the ways in which angular momentum transport is included in simulations of non-self-gravitating accretion disks, concentrating on disks in close binaries. Numerical approaches fall in two basic categories; particle based Lagrangian schemes, and grid based Eulerian techniques. Underlying the choice of numerical technique are assumptions that are made about disk physics, in particular about the angular momentum transport mechanism. Grid-based simulations have generally been of hot, relatively inviscid disks whereas particle-based simulations are more commonly of cool, viscous disks. Calculations of the latter type have been instrumental in developing a model for the superhump phenomenon. We describe how we use an artificial viscosity term to introduce angular momentum transport into our smoothed particle hydrodynamics (SPH) disk code.

We numerically study the precessing disk model for superhump in the SU~UMa subclass of cataclysmic variables, using a two dimensional SPH code specifically designed for thin disk problems. Two disk simulations for a binary with mass ratio $q=\frac{3}{17}$ (similar to OY~Car) are performed, in order to investigate the Lubow (1991 a,b) tidal resonance instability mechanism. In the first calculation, a disk evolves under steady mass transfer from $L_1$. In the second simulation, mass is added in Keplerian orbit to the inner disk. The two disks follow similar evolutionary paths. However the $L_1$ stream-disk interaction is found to slow the disk's radial expansion and to circularise gas orbits. The initial eccentricity growth in our simulations is exponential at a rate slightly less than predicted by Lubow (1991a). We do not observe a clearing of material from the resonance region via the disk's tidal response to the $m=2$ component of the binary potential as was described in Lubow (1992). Instead the $m=2$ response weakens as the disk eccentricty increases. Both disks reach an eccentric equilibrium state, in which they undergo prograde precession. The rate of viscous energy dissipation in the disks has a periodic excess with a period matching the disk's rotation. The source is identified as a large region in the outer disk, and the mechanism by which it is produced is identified. The time taken for the periodic excess to develop is consistent with the first appearance of superhumps in a superoutburst.

We model a massless viscous disk using Smoothed Particle Hydrodynamics (SPH) and note that it evolves according to the Lynden-Bell \& Pringle theory (1974) until a non-axisymmetric instability develops at the inner edge of the disk. This instability may have the same origin as the instability of initially axisymmetric viscous disks discussed by Lyubarskij et al. (1994). To clarify the evolution we evolved single and double rings of particles. It is actually inconsistent with the SPH scheme to set up a single ring as an initial condition because SPH assumes a smoothed initial state. As would be expected from an SPH simulation, the ring rapidly breaks up into a band. We analyse the stability of the ring and show that the predictions are confirmed by the simulation.