Source author record

J. J. Monaghan

J. J. Monaghan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph physics.flu-dyn

Catalog footprint

What is connected

3works

2topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2009arXiv

On the dynamics of swimming linked bodies

In this paper we study the motion of three linked ellipses moving through a viscous fluid in two dimensions. The angles between the ellipses change with time in a specified manner (the gait) and the resulting time varying configuration is similar to the appearance of a swimming leech. We simulate the motion using the particle method Smoothed Particle Hydrodynamics (SPH) which we test by convergence studies and by comparison with the inviscid results of Kanso et al. (2005) and the viscous results of Eldredge (2006, 2007, 2008). We determine how the average speed and power output depends on the amplitude and oscillation frequency of the gait. We find that the results fit simple scaling rules which can related to the analytical results of G. I. Taylor for the swimming of long narrow animals (1952). We apply our results to estimate the speed of a swimming leech with reasonable accuracy, and we determine the minimum power required to propel the bodies at a specified average speed.

preprint2006arXiv

An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes

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.

preprint2003arXiv

Smoothed Particle Magnetohydrodynamics I. Algorithms and tests in one dimension

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.