We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation and failure flows for dynamical behaviour of regolith. In the context of projected small body (asteroid or small moons) sample return missions, we investigate the efficiency and feasibility of a particular material sampling method: Brushes sweep material from the asteroid's surface into a collecting tray. We analyze the influence of different material parameters of regolith such as cohesion and angle of internal friction on the sampling rate. Furthermore, we study the sampling process in two environments by varying the surface gravity (Earth's and Phobos') and we apply different rotation rates for the brushes. We find good agreement of our sampling simulations on Earth with experiments and provide estimations for the influence of the material properties on the collecting rate.

It has been suggested that the comet-like activity of Main Belt Comets are due to the sublimation of sub-surface water-ice that has been exposed as a result of their surfaces being impacted by m-sized bodies. We have examined the viability of this scenario by simulating impacts between m-sized and km-sized objects using a smooth particle hydrodynamics approach. Simulations have been carried out for different values of the impact velocity and impact angle as well as different target material and water-mass fraction. Results indicate that for the range of impact velocities corresponding to those in the asteroid belt, the depth of an impact crater is slightly larger than 10 m suggesting that if the activation of MBCs is due to the sublimation of sub-surface water-ice, this ice has to exist no deeper than a few meters from the surface. Results also show that ice-exposure occurs in the bottom and on the interior surface of impact craters as well as the surface of the target where some of the ejected icy inclusions are re-accreted. While our results demonstrate that the impact scenario is indeed a viable mechanism to expose ice and trigger the activity of MBCs, they also indicate that the activity of the current MBCs is likely due to ice sublimation from multiple impact sites and/or the water contents of these objects (and other asteroids in the outer asteroid belt) is larger than the 5% that is traditionally considered in models of terrestrial planet formation providing more ice for sublimation. We present details of our simulations and discuss their results and implications.

Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request.

So far, multiple stellar systems harbor more than 130 extra solar planets. Dynamical simulations show that the outcome of planetary formation process can lead to various planetary architecture (i.e. location, size, mass and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (hereafter HZ). In this study, we make a comparison of several binary star systems and their efficiency to move icy asteroids from beyond the snow-line into orbits crossing the HZ. We modeled a belt of 10000 asteroids (remnants from the late phase of planetary formation process) beyond the snow-line. The planetesimals are placed randomly around the primary star and move under the gravitational influence of the two stars and a gas giant. As the planetesimals do not interact with each other, we divided the belt into 100 subrings which were separately integrated. In this statistical study, several double star configurations with a G-type star as primary are investigated. Our results show that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result into a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star can not only favor a faster depletion of our disk of planetesimals but can also bring 4 -- 5 times more water into the whole HZ.

We present preliminary results of terrestrial planet formation using on the one hand classical numerical integration of hundreds of small bodies on CPUs and on the other hand -- for comparison reasons -- the results of our GPU code with thousands of small bodies which then merge to larger ones. To be able to determine the outcome of collision events we use our smooth particle hydrodynamics (SPH) code which tracks how water is lost during such events.

As part of a national scientific network 'Pathways to Habitability' the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize the state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).

We investigate the outcome of collisions of Ceres-sized planetesimals composed of a rocky core and a shell of water ice. These collisions are not only relevant for explaining the formation of planetary embryos in early planetary systems, but also provide insight into the formation of asteroid families and possible water transport via colliding small bodies. Earlier studies show characteristic collision velocities exceeding the bodies' mutual escape velocity which - along with the distribution of the impact angles - cover the collision outcome regimes 'partial accretion', 'erosion', and 'hit-and-run' leading to different expected fragmentation scenarios. Existing collision simulations use bodies composed of strengthless material; we study the distribution of fragments and their water contents considering the full elasto-plastic continuum mechanics equations also including brittle failure and fragmentation.

Context. We investigate the effects of including material strength in multi-material planetesimal collisions. Aims. The differences between strengthless material models and including the full elasto-plastic model for solid bodies with brittle failure and fragmentation when treating collisions of asteroid-sized bodies as they occur frequently in early planetary systems are demonstrated. Methods. We study impacts of bodies of Ceres-mass with a solid rock impactor and a target with 30 weight-% water content as surface ice.} The initial impact velocities and impact parameters are varied between the escape velocity $v_\mathrm{esc}$ to about 6 $v_\mathrm{esc}$ and from head-on collisions to close fly-bys, respectively. We simulate the collisions using our own SPH code utilizing both strengthless material and the full elasto-plastic material model including brittle failure. Results. One of the most prominent differences is the higher degree of fragmentation and shattered debris clouds in the solid model. In most collision scenarios however, the final outcomes are very similar and differ primarily by the about one order of magnitude higher degree of fragmentation in the solid case. Also, the survivors tend to be of less mass in the solid case which also predicts a higher water loss than the strengthless hydro model. This may be an effect of the relatively low-energy impacts that cannot destroy the solid material instantly. As opposed to giant impacts we also observe an indication that some water ice gets transferred between the bodies.

We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.

Context.Within the sequential accretion scenario of planet formation, planets are build up through a sequence sticking collisions. The outcome of collisions between porous dust aggregates is very important for the growth from very small dust particles to planetesimals. In this work we determine the necessary material properties of dust aggregates as a function the porosity. Aims: Continuum models such as SPH that are capable of simulating collisions of macroscopic dust aggregates require a set of material parameters. Some of them such as the tensile and shear strength are difficult to obtain from laboratory experiments. The aim of this work is to determine these parameters from ab-initio molecular dynamics simulations. Methods: We simulate the behavior of porous dust aggregates using a detailed micro-physical model of the interaction of spherical grains that includes adhesion forces, rolling, twisting, and sliding. Using different methods of preparing the samples we study the strength behavior of our samples with varying porosity and coordination number of the material. Results: For the tensile strength, we can reproduce data from laboratory experiments very well. For the shear strength, there are no experimental data available. The results from our simulations differ significantly from previous theoretical models, which indicates that the latter might not be sufficient to describe porous dust aggregates. Conclusions: We have provided functional behavior of tensile and shear strength of porous dust aggregates as a function of the porosity that can be directly applied in continuum simulations of these objects in planet formation scenarios.

We carry out three-dimensional Smoothed Particle Hydrodynamics simulations of spherical homogeneous SiO2 dust aggregates to investigate how the mass and the porosity of the aggregates affects their ability to survive an impact at various different collision velocities (between 1 - 27.5m/s). We explore how the threshold velocities for fragmentation vary with these parameters. Crucially, we find that the porosity plays a part of utmost importance in determining the outcome of collisions. In particular, we find that aggregates with filling factors >37% are significantly weakened and that the velocity regime in which the aggregates grow is reduced or even non-existent (instead, the aggregates either rebound off each other or break apart). At filling factors less than ~37% we find that more porous objects are weaker but not as weak as highly compact objects with filling factors >37%. In addition we find that (for a given aggregate density) collisions between very different mass objects have higher threshold velocities than those between very similar mass objects. We find that fragmentation velocities are higher than the typical values of 1m/s and that growth can even occur for velocities as high as 27.5m/s. Therefore, while the growth of aggregates is more likely if collisions between different sized objects occurs or if the aggregates are porous with filling factor <37%, it may also be hindered if the aggregates become too compact.

We give a brief introduction to smoothed particle hydrodynamics methods for continuum mechanics. Specifically, we present our 3D SPH code to simulate and analyze collisions of asteroids consisting of two types of material: basaltic rock and ice. We consider effects like brittle failure, fragmentation, and merging in different impact scenarios. After validating our code against previously published results we present first collision results based on measured values for the Weibull flaw distribution parameters of basalt.

The early planetesimal growth proceeds through a sequence of sticking collisions of dust agglomerates. Very uncertain is still the relative velocity regime in which growth rather than destruction can take place. The outcome of a collision depends on the bulk properties of the porous dust agglomerates. Continuum models of dust agglomerates require a set of material parameters that are often difficult to obtain from laboratory experiments. Here, we aim at determining those parameters from ab-initio molecular dynamics simulations. Our goal is to improveon the existing model that describe the interaction of individual monomers. We use a molecular dynamics approach featuring a detailed micro-physical model of the interaction of spherical grains. The model includes normal forces, rolling, twisting and sliding between the dust grains. We present a new treatment of wall-particle interaction that allows us to perform customized simulations that directly correspond to laboratory experiments. We find that the existing interaction model by Dominik & Tielens leads to a too soft compressive strength behavior for uni and omni-directional compression. Upon making the rolling and sliding coefficients stiffer we find excellent agreement in both cases. Additionally, we find that the compressive strength curve depends on the velocity with which the sample is compressed. The modified interaction strengths between two individual dust grains will lead to a different behaviour of the whole dust agglomerate. This will influences the sticking probabilities and hence the growth of planetesimals. The new parameter set might possibly lead to an enhanced sticking as more energy can be stored in the system before breakup.

In the framework of the coagulation scenario, kilometre-sized planetesimals form by subsequent collisions of pre-planetesimals of sizes from centimetre to hundreds of metres. Pre-planetesimals are fluffy, porous dust aggregates, which are inhomogeneous owing to their collisional history. Planetesimal growth can be prevented by catastrophic disruption in pre-planetesimal collisions above the destruction velocity threshold. We develop an inhomogeneity model based on the density distribution of dust aggregates, which is assumed to be a Gaussian distribution with a well-defined standard deviation. As a second input parameter, we consider the typical size of an inhomogeneous clump. These input parameters are easily accessible by laboratory experiments. For the simulation of the dust aggregates, we utilise a smoothed particle hydrodynamics (SPH) code with extensions for modelling porous solid bodies. The porosity model was previously calibrated for the simulation of silica dust, which commonly serves as an analogue for pre-planetesimal material. The inhomogeneity is imposed as an initial condition on the SPH particle distribution. We carry out collisions of centimetre-sized dust aggregates of intermediate porosity. We vary the standard deviation of the inhomogeneous distribution at fixed typical clump size. The collision outcome is categorised according to the four-population model. We show that inhomogeneous pre-planetesimals are more prone to destruction than homogeneous aggregates. Even slight inhomogeneities can lower the threshold for catastrophic disruption. For a fixed collision velocity, the sizes of the fragments decrease with increasing inhomogeneity. Pre-planetesimals with an active collisional history tend to be weaker. This is a possible obstacle to collisional growth and needs to be taken into account in future studies of the coagulation scenario.

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.

A highly favoured mechanism of planetesimal formation is collisional growth. Single dust grains, which follow gas flows in the protoplanetary disc, hit each other, stick due to van der Waals forces and form fluffy aggregates up to centimetre size. The mechanism of further growth is unclear since the outcome of aggregate collisions in the relevant velocity and size regime cannot be investigated in the laboratory under protoplanetary disc conditions. Realistic statistics of the result of dust aggregate collisions beyond decimetre size is missing for a deeper understanding of planetary growth. Joining experimental and numerical efforts we want to calibrate and validate a computer program that is capable of a correct simulation of the macroscopic behaviour of highly porous dust aggregates. After testing its numerical limitations thoroughly we will check the program especially for a realistic reproduction of various benchmark experiments. We adopt the smooth particle hydrodynamics (SPH) numerical scheme with extensions for the simulation of solid bodies and a modified version of the Sirono porosity model. Experimentally measured macroscopic material properties of silica dust are implemented. We calibrate and test for the compressive strength relation and the bulk modulus. SPH has already proven to be a suitable tool to simulate collisions at rather high velocities. In this work we demonstrate that its area of application can not only be extended to low-velocity experiments and collisions. It can also be used to simulate the behaviour of highly porous objects in this velocity regime to a very high accuracy.The result of the calibration process in this work is an SPH code that can be utilised to investigate the collisional outcome of porous dust in the low-velocity regime.

Recent years have shown many advances in our knowledge of the collisional evolution of protoplanetary dust. Based on a variety of dust-collision experiments in the laboratory, our view of the growth of dust aggregates in protoplanetary disks is now supported by a deeper understanding of the physics involved in the interaction between dust agglomerates. However, the parameter space, which determines the collisional outcome, is huge and sometimes inaccessible to laboratory experiments. Very large or fluffy dust aggregates and extremely low collision velocities are beyond the boundary of today's laboratories. It is therefore desirable to augment our empirical knowledge of dust-collision physics with a numerical method to treat arbitrary aggregate sizes, porosities and collision velocities. In this article, we implement experimentally-determined material parameters of highly porous dust aggregates into a Smooth Particle Hydrodynamics (SPH) code, in particular an omnidirectional compressive-strength and a tensile-strength relation. We also give a prescription of calibrating the SPH code with compression and low-velocity impact experiments. In the process of calibration, we developed a dynamic compressive-strength relation and estimated a relation for the shear strength. Finally, we defined and performed a series of benchmark tests and found the agreement between experimental results and numerical simulations to be very satisfactory. SPH codes have been used in the past to study collisions at rather high velocities. At the end of this work, we show examples of future applications in the low-velocity regime of collisional evolution.

In the recent years we have performed various experiments on the collision dynamics of highly porous dust aggregates and although we now have a comprehensive picture of the micromechanics of those aggregates, the macroscopic understanding is still lacking. We are therefore developing a mechanical model to describe dust aggregate collisions with macroscopic parameters like tensile strength, compressive strength and shear strength. For one well defined dust sample material, the tensile and compressive strength were measured in a static experiment and implemented in a Smoothed Particle Hydrodynamics (SPH) code. A laboratory experiment was designed to compare the laboratory results with the results of the SPH simulation. In this experiment, a mm-sized glass bead is dropped into a cm-sized dust aggregate with the previously measured strength parameters. We determine the deceleration of the glass bead by high-speed imaging and the compression of the dust aggregate by x-ray micro-tomography. The measured penetration depth, stopping time and compaction under the glass bead are utilized to calibrate and test the SPH code. We find that the statically measured compressive strength curve is only applicable if we adjust it to the dynamic situation with a 'softness' parameter. After determining this parameter, the SPH code is capable of reproducing experimental results, which have not been used for the calibration before.

Following the recent insight in the material structure of comets, protoplanetesimals are assumed to have low densities and to be highly porous agglomerates. It is still unclear if planetesimals can be formed from these objects by collisional growth. Therefore, it is important to study numerically the collisional outcome from low velocity impacts of equal sized porous agglomerates which are too large to be examined in a laboratory experiment. We use the Lagrangian particle method Smooth Particle Hydrodynamics to solve the equations that describe the dynamics of elastic and plastic bodies. Additionally, to account for the influence of porosity, we follow a previous developed equation of state and certain relations between the material strength and the relative density. Collisional growth seems possible for rather low collision velocities and particular material strengths. The remnants of collisions with impact parameters that are larger than 50% of the radius of the colliding objects tend to rotate. For small impact parameters, the colliding objects are effectively slowed down without a prominent compaction of the porous structure, which probably increases the possibility for growth. The protoplanetesimals, however, do not stick together for the most part of the employed material strengths. An important issue in subsequent studies has to be the influence of rotation to collisional growth. Moreover, for realistic simulations of protoplanetesimals it is crucial to know the correct material parameters in more detail.

We argue that the quiescent value of the viscosity parameter of the accretion disc in WZ Sge may be alpha_{cold}~0.01, in agreement with estimates of alpha_{cold} for other dwarf novae. Assuming the white dwarf in WZ Sge to be magnetic, we show that, in quiescence, material close to the white dwarf can be propelled to larger radii, depleting the inner accretion disc. The propeller therefore has the effect of stabilizing the inner disc and allowing the outer disc to accumulate mass. The outbursts of WZ Sge are then regulated by the (magnetically determined) evolution of the surface density of the outer disc at a radius close to the tidal limit. Numerical models confirm that the recurrence time can be significantly extended in this way. The outbursts are expected to be superoutbursts since the outer disc radius is forced to exceed the tidal (3:1 resonance) radius. The large, quiescent disc is expected to be massive, and to be able to supply the observed mass accretion rate during outburst. We predict that the long-term spin evolution of the white dwarf spin will involve a long cycle of spin up and spin down phases.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

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 have performed Smoothed Particle Hydrodynamics (SPH) simulations to study the time evolution of one and two protoplanets embedded in a protoplanetary accretion disc. We investigate accretion and migration rates of a single protoplanet depending on several parameters of the protoplanetary disc, mainly viscosity and scale height. Additionally, we consider the influence of a second protoplanet in a long time simulation and examine the migration of the two planets in the disc, especially the growth of eccentricity and chaotic behaviour. One aim of this work is to establish the feasibility of SPH for such calculations considering that usually only grid-based methods are adopted. To resolve shocks and to prevent particle penetration, we introduce a new approach for an artificial viscosity, which consists of an additional artificial bulk viscosity term in the SPH-representation of the Navier-Stokes equation. This allows for an accurate treatment of the physical kinematic viscosity to describe the shear, without the use of an artificial shear viscosity.

We perform three-dimensional, Newtonian hydrodynamic simulations with a nuclear equation of state to investigate the accretion dynamics in neutron star black hole systems. We find as a general result that non-spinning donor stars yield larger circularization radii than corotating donors. Therefore, the matter from a neutron star without spin will more likely settle into an accretion disk outside the Schwarzschild radius. With the used stiff equation of state we find it hard to form an accretion disk that is promising to launch a gamma-ray burst. In all relevant cases the core of the neutron star survives and keeps orbiting the black hole as a mini neutron star for the rest of the simulation time (up to several hundred dynamical neutron star times scales). The existence of this mini neutron star leaves a clear imprint on the gravitational wave signal which thus can be used to probe the physics at supra-nuclear densities.

We argue that the outbursts of the FU Orionis stars occur on timescales which are much longer than expected from the standard disc instability model with \alpha_{c} \gtrsim 10^{-3}. The outburst, recurrence, and rise times are consistent with the idea that the accretion disc in these objects is truncated at a radius R_{i} \sim 40 \rsun. In agreement with a number of previous authors we suggest that the inner regions of the accretion discs in FU Ori objects are evacuated by the action of a magnetic propeller anchored on the central star. We develop an analytic solution for the steady state structure of an accretion disc in the presence of a central magnetic torque, and present numerical calculations to follow its time evolution. These calculations confirm that a recurrence time that is consistent with observations can be obtained by selecting appropriate values for viscosity and magnetic field strength.

We investigate the spin evolution of the unusual magnetic CV AE Aqr. As a prototype for a potentially large population of CVs subject to a thermally unstable phase of mass transfer, understanding its future is crucial. We present a new definition of the magnetospheric radius in terms of the white dwarf's spin period, and use this along with numerical simulations to follow the spin evolution of AE Aqr. We also present preliminary SPH results suggesting the existence of a stable propeller state. These results highlight the complexity of mCVs and may provide am improved understanding of the evolution of all types of CVs.

We study analytically and numerically the stability of the pressure-less, viscously spreading accretion ring. We show that the ring is unstable to small non-axisymmetric perturbations. To perform the perturbation analysis of the ring we use a stretching transformation of the time coordinate. We find that to 1st order, one-armed spiral structures, and to 2nd order additionally two-armed spiral features may appear. Furthermore, we identify a dispersion relation determining the instability of the ring. The theoretical results are confirmed in several simulations, using two different numerical methods. These computations prove independently the existence of a secular spiral instability driven by viscosity, which evolves into persisting leading and trailing spiral waves. Our results settle the question whether the spiral structures found in earlier simulations of the spreading ring are numerical artifacts or genuine instabilities.