**· Find Similar Abstracts (with default settings below)****· Full Refereed Journal Article (PDF/Postscript)****· Full Refereed Scanned Article (GIF)****· arXiv e-print**(arXiv:astro-ph/9409058)**· References in the article****· Citations to the Article (281)**(Citation History)**· Refereed Citations to the Article****· Also-Read Articles**(Reads History)- ·
**· Translate This Page**

Title: | Hydra: an Adaptive-Mesh Implementation of P 3M-SPH | |

Authors: | Couchman, H. M. P.; Thomas, P. A.; Pearce, F. R. | |

Publication: | Astrophysical Journal v.452, p.797 (ApJ Homepage) | |

Publication Date: | 10/1995 | |

Origin: | APJ; KNUDSEN | |

Astronomy Keywords: | HYDRODYNAMICS, METHODS: NUMERICAL | |

DOI: | 10.1086/176348 | |

Bibliographic Code: | 1995ApJ...452..797C |

The code evolves a mixture of purely gravitational particles and gas particles. SPH gas forces are calculated in the standard way from near neighbors. Gravitational forces are calculated using the mesh refinement scheme described by Couchman (1991). The AP

The choice of time integration scheme is investigated in detail. We find that a simple single-step predictor-corrector type integrator, which is equivalent to Leapfrog for velocity-independent forces, is the most efficient.

A method for generating an initial distribution of particles by allowing a uniform temperature gas of SPH particles to relax within a periodic box is presented. The average SPH density that results varies by ≈ ±1.3%. This is the fluctuation amplitude on roughly the Nyquist frequency; for smaller wavenumbers the fluctuations have lower amplitudes. We present a modified form of the Layzer-Irvine equation which includes the thermal contribution of the gas together with radiative cooling.

The SPH and time integration schemes were tested and compared by running a series of tests of sound waves and shocks. These tests were also used to derive time-step constraints sufficient to ensure both energy and entropy conservation.

We have compared the results of simulations of spherical infall and collapse with varying numbers of particles. We show that many thousands of particles are necessary in a halo to correctly model the collapse.

As a further test, the cluster simulation of Thomas & Couchman (1992) has been rerun with the new code, which includes a number of improvements in the SPH implementation. We find close agreement except in the core properties of the cluster which are strongly affected by entropy scatter in the older simulation. This demonstrates the crucial importance of conserving entropy in SPH simulations.

Bibtex entry for this abstract Preferred format for this abstract (see Preferences) |