Thank you for taking the time to visit my site. I am interested in Computers old and new, programming and software development on HPC with Fortran, MPI, and python, numerical algorithms and simulations, electronics and soldering, and generally anything related to technology. A couple of projects I have worked on are presented below.

Spherical Shell Convection

The molten iron in the Earth’s outer core undergoes turbulent convection. Much of what we know about physical systems such as the Earth’s core comes from computer simulations. To understand more about the proceses happening under our planet’s surface, we can numerically solve the equations governing fluid flow in Earth-like geometry.

The video above shows convection in an Earth-like system, but at much less extreme conditions, since it is currently numerically impossible to simulate the full Earth-system. Rayleigh Number, which determines how much thermal driving is in the system: 2x10⁷, Ekman Number, which shows how fast the Earth is spinning (smaller means faster): 1.3x10⁻⁴, the magnetic Prandtl number Pm is 11, while the thermal Prandtl number Pr is 1. Solved using 120 points in radius, 144 points in co-latitude, and 288 points in azimuth.

Fluid Dynamics

The Navier-Stokes equations describe the way that fluids behave. In all but the most simple of cases, analytical solutions are impossilble, and numerical solutions are required. One particularly well studied and interesting problem is that of Rayleigh-Bénard convection (RBC). RBC describes the buoyancy driven flow of a fluid heated from below and cooled from above. It can be uses as a simplified analogue for understanding a range of astrophysical and geophysical systems, such as the Sun, stars, planetary interiors and the atmosphere.

The above video shows Boussinesq Rayleigh-Bénard convection with fixed flux boundary conditions at a rather high Rayleigh number (5e10). The equations were solved using pseudospectral spatial discretisation using Fourier series in the horizontal direction and Chebyshev polynomials in the vertical. Time stepping was achieved using a Runge-Kutta IMplicit-EXplicit (IMEX) scheme of 3rd order (RK443). The code was written in python for the open-source spectral solver Dedalus.

GameBoy game

I wrote a GameBoy game! Here is a quick teaser video of “SUPER SHOOTER”:

I started this project to satisfy my curiosity in three main areas: old 8-bit computers, experience in C-programming language, and a nostalgia hit for the original GameBoy. The code was developed using the GameBoy Development Kit: gbdk, with testing and debugging taken care of by the bgb emulator. I was able to capture footage using the mGBA Emulator. I aim in future to port this game/write a new game in the updated GBDK-2020.