About me

My papers

Vernazza, P.; Jorda, L.; Ševeček, P.; Brož, M.; et al. A basin-free spherical shape as an outcome of a giant impact on asteroid Hygiea

Ševeček, P.; Brož, M.; Jutzi, M. Impacts into rotating targets: angular momentum draining and efficient formation of synthetic families

Ševeček, P.; Brož, M.; Nesvorný, D.; Enke, B.; et al. SPH/N-Body simulations of small (D = 10km) asteroidal breakups and improved parametric relations for Monte-Carlo collisional models

Ševeček, P.; Golubov, O.; Scheeres, D. J.; Krugly, Yu. N. Obliquity dependence of the tangential YORP

Ševeček, P.; Brož, M.; Čapek, D.; Ďurech, J. The thermal emission from boulders on (25143) Itokawa and general implications for the YORP effect


OpenSPH code

OpenSPH runs hydrodynamical and N-body simulations and was written for asteroid collisions and subsequent gravitational evolution. The code offers SPH and N-body solvers with several different equations of state and material rheologies. It is written in C++14 with a modular object-oriented design, focused on extensibility and maintainability, and it can be used either as a library or as a standalone graphical program that allows to set up the problem in a convenient graphical node editor. The graphical program further allows real-time visualization of the simulation, diagnostics and tools for analysis of the results.

Code snippet

AsymmetricSolver solver(*scheduler, settings, equations);
for (Size matId = 0; matId < particles->getMaterialCnt(); ++matId) {
  solver.create(*particles, particles->getMaterial(matId));
PredictorCorrector stepper(particles, settings);

Statistics stats;
for (Float t = 0; t < t_end; t += stepper.getTimeStep()) {
  stepper.step(*scheduler, solver, stats);