Current features in SEREN
- 3 'flavours' of SPH :
- 1, 2 or 3 dimensions
- Gas, inter-cloud particles or boundary particles
(static or frog-marched)
- Barnes-Hut Octal tree for efficient neighbour finding
and calculation of gravitational forces
(Barnes & Hut 1986)
- Periodic boundary conditions
(independent for each dimension)
- Ewald periodic gravity
(Hernquist, Bouchet & Suto 1991,
Klessen 1997)
- Various integration schemes :
- 2nd order Leapfrog kick-drift-kick and drift-kick-drift
- 2nd order Runge-Kutta
- Block timesteps and neighbour-timestep checking
for safety
(Saitoh & Makino 2009)
- Various flavours of artificial dissipation :
- Equations of state :
- Ideal-gas EoS with compressional heating
- Isothermal equation of state
- Polytropic equation of state
- Barotropic equation of state
- Radiative cooling approximation
(Stamatellos et al. 2007)
- Hybrid radiative-cooling/flux-limited diffusion approximation
(Forgan et al. 2009)
- HEALPix ionizing radiation
(Bisbas et al. 2009)
- Sink particles
- 4th-order Hermite integrator for N-body evolution
of star orbital properties (Makino & Aarseth 1992) and automatic calculation
of binary properties
- Hybrid SPH and 4th-order Hermite integrator for simulating gaseous clusters (Hubber, Allison, Smith & Goodwin 2013)
- Parallelized using OpenMP
- Parallelized using MPI (limited beta-implementation)
- Produces output compatable with the
SPLASH visualisation program.