Pavel Ševeček

About me

PhD thesis: Simulations of asteroid collisions using a hybrid SPH/N-body approach

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

Teaching

Numerical solution of partial differential equations (selected chapters, fall/winter 2019)
Lectures cover basics of numerical solution of PDEs, using both Eulerian (FDM, FEM) and Lagrangian (SPH) approaches.
Presentation slides Example codes Videos OpenSPH simulations

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);

                                    }