The equations of smoothed particle magnetohydrodynamics (SPMHD), even with the various corrections to instabilities so far proposed, have been observed to be unstable when a very steep density gradient is necessarily combined with a variable smoothing length formalism. Here we consider in more detail the modifications made to the SPMHD equations in LBP2015 that resolve this instability by replacing the smoothing length in the induction and anisotropic force equations with an average smoothing length term. We then explore the choice of average used and compare the effects on a test `cylinder-in-a-box' problem and the collapse of a magnetised molecular cloud core. We find that, aside from some benign numerical effects at low resolutions for the quadratic mean, the formalism is robust as to the choice of average but that in complicated models it is essential to apply the average to both equations; in particular, all four averages considered exhibit similar conservation properties. This improved formalism allows for arbitrarily small sink particles and field geometries to be explored, vastly expanding the range of astronomical problems that can be modeled using SPMHD.

In this paper we describe an adaptive softening length formalism for collisionless N-body and self-gravitating Smoothed Particle Hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially variable softening lengths can be used in N-body calculations without secular increases in energy. The formalism requires the calculation of a small additional term to the gravitational force related to the gradient of the softening length. The extra term is similar in form to the usual SPH pressure force (although opposite in direction) and is therefore straightforward to implement in any SPH code at almost no extra cost. For N-body codes some additional cost is involved as the formalism requires the computation of the density via a summation over neighbouring particles using the smoothing kernel. The results of numerical tests demonstrate that, for homogeneous mass distributions, the use of adaptive softening lengths gives a softening which is always close to the `optimal' choice of fixed softening parameter, removing the need for fine-tuning. For a heterogeneous mass distribution (as may be found in any large scale N-body simulation) we find that the errors on the least-dense component are lowered by an order of magnitude compared to the use of a fixed softening length tuned to the densest component. For SPH codes our method presents a natural and elegant choice of softening formalism which makes a small improvement to both the force resolution and the total energy conservation at almost zero additional cost.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

Toy Stars are gas masses where the compressibility is treated without approximations but gravity is replaced by a force which, for any pair of masses, is along their line of centres and proportional to their separation. They provide an invaluable resource for testing the suitability of numerical codes for astrophysical gas dynamics. In this paper we derive the equations for both small amplitude oscillations and non linear solutions for rotating and pulsating Toy Stars in two dimensions, and show that the solutions can be reduced to a small number of ordinary differential equations. We compare the accurate solutions of these equations with Smoothed Particle Hydrodynamics (SPH) simulations. The two dimensional Toy Star solutions are found to provide an excellent benchmark for SPH algorithms, highlighting many of the strengths and also some weaknesses of the method.

We describe a new, faster implicit algorithm for solving the radiation hydrodynamics equations in the flux-limited diffusion approximation for smoothed particle hydrodynamics. This improves on the method elucidated in Whitehouse & Bate by using a Gauss-Seidel iterative method rather than iterating over the exchange of energy between pairs of particles. The new algorithm is typically many thousands of times faster than the old one, which will enable more complex problems to be solved. The new algorithm is tested using the same tests performed by Turner & Stone for ZEUS-2D, and repeated by Whitehouse & Bate.

In two previous papers (Price & Monaghan 2004a,b) (papers I,II) we have described an algorithm for solving the equations of Magnetohydrodynamics (MHD) using the Smoothed Particle Hydrodynamics (SPH) method. The algorithm uses dissipative terms in order to capture shocks and has been tested on a wide range of one dimensional problems in both adiabatic and isothermal MHD. In this paper we investigate multidimensional aspects of the algorithm, refining many of the aspects considered in papers I and II and paying particular attention to the code's ability to maintain the div B = 0 constraint associated with the magnetic field. In particular we implement a hyperbolic divergence cleaning method recently proposed by Dedner et al. (2002) in combination with the consistent formulation of the MHD equations in the presence of non-zero magnetic divergence derived in papers I and II. Various projection methods for maintaining the divergence-free condition are also examined. Finally the algorithm is tested against a wide range of multidimensional problems used to test recent grid-based MHD codes. A particular finding of these tests is that in SPMHD the magnitude of the divergence error is dependent on the number of neighbours used to calculate a particle's properties and only weakly dependent on the total number of particles. Whilst many improvements could still be made to the algorithm, our results suggest that the method is ripe for application to problems of current theoretical interest, such as that of star formation.

In this paper we present the equations and some basic applications of relativistic SPH. The equations are generated under the assumption that there are no interaction terms between the fluid modelled and the background space-time, i.e. we are neglecting the perturbations to the metric produced by the modelled fluid. This corresponds to a stationary metric, which for this work, we further limit to be static as well. The equations use a new signal velocity term and artificial viscosity to smear out the effects of strong shocks and are tested here against Newtonian and relativistic shocks.

In this paper we show how the Smoothed Particle Hydrodynamics (SPH) equations for ideal magnetohydrodynamics (MHD) can be written in conservation form with the positivity of the dissipation guaranteed. We call the resulting algorithm Smoothed Particle Magnetohydrodynamics (SPMHD). The equations appear to be accurate, robust and easy to apply and do not suffer from the instabilities known to exist previously in formulations of the SPMHD equations. In addition we formulate our MHD equations such that errors associated with non-zero divergence of the magnetic field are naturally propagated by the flow and should therefore remain small. In this and a companion paper (Price and Monaghan 2003b) we present a wide range of numerical tests in one dimension to show that the algorithm gives very good results for one dimensional flows in both adiabatic and isothermal MHD. For the one dimensional tests the field structure is either two or three dimensional. The algorithm has many astrophysical applications and is particularly suited to star formation problems.

In this paper we show how a Lagrangian variational principle can be used to derive the SPMHD (smoothed particle magnetohydrodynamics) equations for ideal MHD. We also consider the effect of a variable smoothing length in the SPH kernels after which we demonstrate by numerical tests that the consistent treatment of terms relating to the gradient of the smoothing length in the SPMHD equations significantly improves the accuracy of the algorithm. Our results complement those obtained in a companion paper (Price and Monaghan 2003a, paper I) for non ideal MHD where artificial dissipative terms were included to handle shocks.

There have been some issues in the past in attempts to simulate magnetic fields using the Smoothed Particle Hydrodynamics (SPH) method. SPH is well suited to star formation problems because of its Lagrangian nature. We present new, stable and conservative methods for magnetohydrodynamics (MHD) in SPH and present numerical tests on both waves and shocks in one dimension to show that it gives robust and accurate results.

In this paper an SPH version of the alpha turbulence model devised by Holm and his colleagues is formulated for compressible flow with a resolution that varies in space and time. The alpha model involves two velocity fields. One velocity field is obtained from the momentum equation, the other by averaging this velocity field as in the version of SPH called XSPH. The particles (fluid elements) are moved with the averaged velocity. In analogy to the continuum alpha model we obtain a particle Lagrangian from which the SPH alpha equations can be derived. The system satisfies a discrete Kelvin circulation theorem identical to that obtained with no velocity averaging. In addition, the energy, linear and angular momentum are conserved. We show that the continuum equivalent of the SPH equations are identical to the continuum alpha model, and we conjecture that they will have the same desirable features of the continuum modelincluding the reduction of energy in the high wave number modes even when the dissipation is zero. Regardless of issues concerning turbulence modelling, the SPH alpha model is a powerful extension of the XSPH algorithm which reduces disorder at short length scales and retains the constants of the motion. The SPH alpha model is simple to implement.

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.