Generating Optimal Initial Conditions for Smoothed Particle   Hydrodynamics Simulations
  We review existing SPH setup methods and outline their advantages, limitations and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics (SPH) simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tesselations. Particles are placed and iteratively moved based on their proximity to neighboring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements. 
Steven Diehl, Gabriel Rockefeller, Christopher L. Fryer, David   Riethmiller, Thomas S. Statler
astro-ph.IM astro-ph.CO astro-ph.GA