MHD Implemention Prep

General Relativity (GR)

Nothing too exciting. Starting hard into curved spacetime.

Dotfiles Presentation

On Wednesday I gave a presentation on how I organize dotfiles to the computational astrophysics seminar. It went well and now everybody has a scalable and portable way of organizing and moving their dotfiles from machine to machine. The skeleton of my dotfiles can be found here: dotfiles repo

HLLC Solver

I updated my HLLC solver to take a positional argument so that I can test the solver on it’s own the same way I can with Toro’s exact Riemann solver. As you can see in the below image it doesn’t exactly reconstruct all the states on it’s own, though it does work perfectly in an actual simulation. It does show that the estimate of the shock and contact wave speeds are a little low. This behavior seems to persist no matter what choice of wave speeds I use.

exact vs HLLC Top: The output of the exact Riemann solver. Note the perfect reconstruction of the rarefaction fan and the position of the shock wave. Bottom: The output of the HLLC Riemann solver. Note how the fan is not resolved at all and the shock speed is a little underestimated

MHD Solver

I started implementing my MHD solver this week. So far I’ve just been working on setting the ground work for the HLLD solver and VL+CT integrator; there are lots of lots of small changes to make to the code to accommodate the additional complexities of magnetohydrodynamics and keeping the documentation updated requires care. Progress is good though so I am hopeful.

Original HLLD paper: Miyoshi & Kusano 2005

Blog post on the HLLD Algorithm: HLLD Algorithm

Original VL+CT paper: Stone & Gardiner 2009

This post is licensed under CC BY 4.0 by the author.