Scientific discovery via numerical simulations is important in modern astrophysics. This relatively new branch of astrophysics has become possible due to the development of reliable numerical... Show moreScientific discovery via numerical simulations is important in modern astrophysics. This relatively new branch of astrophysics has become possible due to the development of reliable numerical algorithms and the high performance of modern computing technologies. These enable the analysis of large collections of observational data and the acquisition of new data via simulations at unprecedented accuracy and resolution. Ideally, simulations run until they reach some pre-determined termination condition, but often other factors cause extensive numerical approaches to break down at an earlier stage. In those cases, processes tend to be interrupted due to unexpected events in the software or the hardware. In those cases, the scientist handles the interrupt manually, which is time-consuming and prone to errors. We present the Simulation Monitor (SiMon) to automatize the farming of large and extensive simulation processes. Our method is light-weight, it fully automates the entire workflow management, operates concurrently across multiple platforms and can be installed in user space. Inspired by the process of crop farming, we perceive each simulation as a crop in the field and running simulation becomes analogous to growing crops. With the development of SiMon we relax the technical aspects of simulation management. The initial package was developed for extensive parameter searchers in numerical simulations, but it turns out to work equally well for automating the computational processing and reduction of observational data reduction. Show less
Voort, F. van de; Bahé, Y.M.; Bower, R.G.; Correa, C.A.; Crain, R.A.; Schaye, J.; Theuns, T. 2017
We study the evolution of star clusters on circular and eccentric orbits using direct N-body simulations. We model clusters with initially N = 8k and 16k single stars of the same mass, orbiting... Show moreWe study the evolution of star clusters on circular and eccentric orbits using direct N-body simulations. We model clusters with initially N = 8k and 16k single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same N on a circular orbit. We show that then, the evolution of bound particle number and half-mass radius is approximately independent of eccentricity. Secondly, when we scale our results to orbits with the same semimajor axis, we find that the lifetimes are, to first order, independent of eccentricity. When the results of Baumgardt and Makino for a singular isothermal halo are scaled in the same way, the lifetime is again independent of eccentricity to first order, suggesting that this result is independent of the galactic mass profile. From both sets of simulations, we empirically derive the higher order dependence of the lifetime on eccentricity. Our results serve as benchmark for theoretical studies of the escape rate from clusters on eccentric orbits. Finally, our results can be useful for generative models for cold streams and cluster evolution models that are confined to spherical symmetry and/or time-independent tides, such as Fokker-Planck models, Monte Carlo models, and (fast) semi-analytic models. Show less
Astrophysical research in recent decades has made significant progress thanks to the availability of various N-body simulation techniques. With the rapid development of high-performance computing... Show moreAstrophysical research in recent decades has made significant progress thanks to the availability of various N-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to use the computing power of massively parallel clusters with more than 105 GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. In this paper, we propose an adaptive storage scheme for simulation data, inspired by the block time step (BTS) integration scheme found in a number of direct N-body integrators available nowadays, as an urgent response to these challenges. The proposed scheme, namely, the BTS storage scheme, works by minimizing the data redundancy by assigning individual output frequencies to the data as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct N-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used. Show less
Ploeckinger, S.; Recchi, S.; Hensler, G.; Kroupa, P. 2015