A numerical advection solver for power laws
(Supervision: David Huber & Ralf Kissmann)


In our group, we study particle acceleration at shocks in different astrophysical systems. The shape of the energy spectrum resulting from this acceleration process is determined by the balance of acceleration and energy losses. In this thesis, we focus on efficient numerical modelling of the energy spectrum. Here, we use a semi-Lagrangian code to efficiently describe the energy-loss processes. Most of such codes, however, are optimsed for a smooth variation of the dynamical variables, wheras the acceleration process leads to a power-law in energy. In the bachelor thesis, an alternative implementation of the solver will be implemented and tested, that aims at a better representation of such power-laws. For this, a first implementation of the solver together with simple test cases will be provided.


Back to topics

Nach oben scrollen