General Relativity (GR)
Nothing too exciting. Starting hard into curved spacetime.
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
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.
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
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