Collin
Wittenstein

Geothermal Well Simulations

Deep borehole heat exchangers extract geothermal energy by circulating fluid through closed-loop pipes that extend several kilometers into the Earth. This sounds nice, but simulating them is unfortunately not easy: you're dealing with multi-scale geometry where boreholes are centimeters wide but kilometers deep, and you need to simulate years to decades of operation to understand thermal performance. For arrays of multiple wells interacting thermally over these timescales, it was traditionally thought that a fully three-dimensional simulation of the arrays would be computationally unfeasible. Luckily for me I didn't know this before starting to work on the project.

So, I developed a GPU-accelerated three-dimensional model simulating geothermal well arrays. The key innovation is an operator splitting strategy tailored to the problem's physics. Vertical diffusion gets stabilized explicit Runge-Kutta-Chebyshev methods, horizontal diffusion gets ADI schemes, advection gets semi-Lagrangian treatment. Each operator handled by whatever works best for it - sounds obvious in hindsight, but figuring out this combination took a while.

I released this as the open-source Julia package GeothermalWells.jl, using vendor-agnostic GPU kernels that run on both NVIDIA and AMD hardware. Full three-dimensional simulations of multi-well arrays over decades can now run on a single GPU, enabling systematic design space exploration and optimization of geothermal well systems. Supervised by Prof. Dr. Alan Edelman, Prof. Dr. Robert Metcalfe, and Prof. Dr. Hendrik Ranocha.

Water Wave Simulations

I worked on dispersive shallow water equations. The following visualization shows the evolution of a wave over Gaussian bathymetry, using the hyperbolic approximation of the Serre–Green–Naghdi equations in two dimensions. The simulation employs a new energy-conserving semi-discretization that maintains physical properties over long time integrations. I also worked quite a bit on DispersiveShallowWater.jl, a Julia package that implements structure-preserving numerical methods for one-dimensional dispersive shallow water models. Supervised by Prof. Dr. Hendrik Ranocha.

Rattleback Simulation

For my (maths) bachelor thesis, I simulated the movement of the rattleback (German: Wackelstein), a semi-ellipsoid with an inhomogeneous mass distribution that exhibits quite counterintuitive behavior: it can spontaneously reverse its spinning direction. This phenomenon defies intuition but follows directly from the solutions of its equations of motion. Supervised by Prof. Dr. Hendrik Ranocha and Univ.-Prof. Dr. Elmar Schömer.

Image: Ed Keath, CC BY-SA 3.0, via Wikimedia Commons

For this work, I received the best bachelor's thesis in mathematics of the year award.