Particle initialization through a novel packing algorithm

A novel particle packing algorithm has been derived basing on a simplified standard SPH scheme and is used to initialize the particle distribution for generic SPH solvers. Using some intrinsic features of the SPH schemes, the proposed algorithm leads to a final particle distribution that is very stable and is not affected by a further resettlement during the early stages of the evolution. Moreover, the computational costs maintain very low since the final configuration is attained very quickly. Finally, a strong point of the packing algorithm is that it can be easily derived using whatever SPH scheme.