Christoph's photo Christoph Pfrommer
Back to Lectures overview

    Modern Computational Astrophysics: Concepts and Applications

    The objectives of this course are to endow students with the capacity to identify and classify common numerical problems in modern astrophysics. The course aims at an active understanding of numerical methods and algorithms as well as their ranges of applicability. Solving basic astrophysical problems with adequate numerical techniques and determining the range of validity is an essential part of the course. This course (PHY-755) presents theoretical concepts which are then practiced in course PHY-765. A participation of both courses is strongly recommended.


    Prof. Christoph Pfrommer & Dr. Philipp Girichidis, tutorials by Timon Thomas

    Place and Dates:

    The lectures take place every Tuesday, at 10:15-11:45 in room The exercise classes take place every Thursday, at 10:15-11:45 in room

    Not that this course is complemeted by a Seminar on Computational Astrophysics that takes place immediately after the lectures every Tuesday, at 12:15-13:45 in room

    The lectures are based on the revised lecture notes originally written by Prof. Springel (Max-Planck Institute for Astrophysics) and are available as a PDF file (password protected).


    • Reviewing basic concepts of numerical simulations (preparatory course)
      • Integer arithmetic
      • Floating point arithmetic

    • Integration of ordinary differential equations
      • The Euler method: explicit and implicit
      • Runge-Kutta methods and adaptive step sizes
      • The leapfrog and symplectic integrators

    • Basic gas dynamics
      • Basic conservation laws
      • Euler and Navier-Stokes equations
      • Shocks, fluid instabilities and turbulence

    • Eulerian hydrodynamics
      • Types and solution schemes of partial differential equations
      • Advection and the Riemann problem
      • Finite volume discretization
      • Godunov's method and Riemann solvers
      • Extensions to multiple dimensions and high-order accuracy

    • Smoothed particle hydrodynamics (SPH)
      • Kernel interpolants
      • Variational derivation of SPH
      • Artificial viscosity
      • Advantages and disadvantages of SPH

    • Diffusion equation
      • Explicit method
      • Implicit methods

    • Fourier transform techniques
      • Convolution problems
      • The discrete Fourier transform
      • Non-periodic problems with 'zero padding'

    • Collisionless particle systems
      • N-particle ensembles
      • Uncorrelated (collisionless) systems
      • N-body models of collisionless systems

    • Tree algorithms
      • Multipole expansion
      • Hierarchical grouping
      • Tree walk

    • The particle-mesh technique
      • Mass/charge assignment
      • Solving for the potential and force calculation
      • Interpolating from the mesh to the particles

    • Paralellization techniques
      • Hardware overview
      • Shared memory parallelization with OpenMP
      • Distributed memory parallelization with MPI

    Lecture slides on Parallel Computing

    Homework Assignments

    All homework assignments and sample code can be obtained via this link (password protected).

    Credit Points:

    Credit points for this lectures can be used for the Master Science Astrophysics Modules PHY-755 and PHY-765.

    Grades and credit points are given on the basis of an oral exam. To be eligible for this exam, you need to obtain at least 50% of the points of your homework assignments.