Switch to English
Cursus: NS-256B
Numerieke methoden voor fysici en astronomen
Cursus informatieRooster
Studiepunten (ECTS)7,5
Categorie / Niveau2 (Bachelor Verdiepend)
CursustypeCursorisch onderwijs
Aangeboden doorFaculteit Betawetenschappen; Undergraduate School Bètawetenschappen;
Contactpersoondr. W.J. van de Berg
Contactpersoon van de cursus
dr. W.J. van de Berg
Overige cursussen docent
dr. W.J. van de Berg
Overige cursussen docent
prof. dr. ir. M. Dijkstra
Overige cursussen docent
dr. F.M.L.C. Freire
Overige cursussen docent
2  (12-11-2018 t/m 01-02-2019)
TimeslotB: DI-ochtend, DO-middag, DO-namiddag
Cursusinschrijving geopendvanaf 17-09-2018 t/m 30-09-2018
Inschrijven via OSIRISJa
Inschrijven voor bijvakkersJa
Na-inschrijving geopendvanaf 22-10-2018 t/m 23-10-2018

On completion of the course, a student:
•          Knows the relevance of numerical methods for physics and astronomy.
•          Can use basic numerical methods for integration, differentiation and finding roots.
•          Can use Monte Carlo methods for evaluating multi-dimensional integrals.
•          Can integrate numerically systems of first order ordinary differential equations for both initial value problems and boundary value problems.
•          Can numerically solve partial differential equations of intermediate complexity using elementary techniques.
•          Can assess the effect of the numerical method applied on numerical stability and on the validity of the results.
•          Is able to analyse and report results of numerical experiments.
•          Will have had direct experience in programming in Python.
The course starts with an introduction to Python, a computer language widely used in scientific research. Next, the three main groups of mathematical problems are dealt in three separate projects. The setup of the course is very “hands-on”: in the introductory lecture for each project the main concepts are explained, the remainder of the time is reserved for computer practica. For each project, you will write your own numerical program, analyze and visualise the results and write a brief report.
Project 1: In the first project, we will use basic numerical methods like integration, differentiation, and finding roots to solve simple physics problems. We will use a root-finding method to determine the angle with which an object should be launched in space in order to reach a certain position. Subsequently, we will use numerical integration of an ordinary differential equation to determine the motion (or phase space trajectory) of a chaotic pendulum and of a classical multi-particle system. We will compare the time averages of an observable with ensemble averages by using Monte Carlo simulations on the same system.     
Project 2: Numerical integration of ordinary differential equations (ODEs)
The motion of particles that interact with each other, for example, through gravitation, can be described with ODE’s if relativistic effects are neglected. In this second project, the orbits in a three-body system is analysed and, time permitting, next the structure of stellar interiors. Each case illustrates one of the two standard ways of implementing the numerical integration of ODEs, initial value problem and boundary value problem. In the first problem, we use the restricted three-body system model to determine the orbits of the astronomical objects using initial conditions, location of objects at initial time. In the second problem, we use the polytrope equation, an approximation to the equation of state for compact stellar objects, such as white dwarfs and neutron stars, to determine the pressure and mass density in their interior as a function of the distance to their center using boundary conditions, on their surface or center.
Project 3: Numerical integration of partial differential equations (PDEs)
PDEs describe the evolution of every continuous quantity, e.g., fluid flow, magnetic fields, probability density functions. In the last project, PDEs are numerically solved by discretization of space and time. We discuss:
·         The difference between implicit and explicit time discretization methods.
·         The accuracy of the applied time and space discretisation.
·         An analytical method to assess the stability of a discretisation.
·         The impact of time step length, grid resolution and applied discretisation on the quality of the model results.
This theory will be applied on two simple problems and one more complex problem, for example, the shallow water equations.

Infinitesimaalrekening, lineaire algebra
Verplicht materiaal
Aanbevolen materiaal
Python (software alleen in CLZ)


Minimum cijfer-

The final grade is the mean of the three project reports. However, each individual project report should not be graded lower than a 5 for a final grade of 6 or higher. Students may improve one project report as retake.

Switch to English