This is a code that I wrote for fun. I performed a few standard simulations with it. The parameter files for these simulations are provided in the repo.

Running the Code

The code can be run with ./setup [arguments].
For help use ./setup --help

Sedov-Blast wave problem

Simulation is performed using a 400 x 600 cell grid with 2 different reconstruction method: WENO3 and WENO5. Need to fix the scale on python plot.

https://github.com/user-attachments/assets/4454bd73-7b9f-4a44-978a-450a5eedff3b

Kelvin-Helmholtz instability problem

Performed using a 256 x 256 grid.

https://github.com/user-attachments/assets/7ac78e27-1890-4721-a03f-916765af4de7

Rayleigh-Taylor instability problem

Performed using a 100 x 300 grid.

https://github.com/user-attachments/assets/fccd9944-18d1-4832-b176-2fd8b6b9648d

2D Riemann problem

https://github.com/user-attachments/assets/940f2f6d-cf63-4eb4-8496-b0084b3e8611

2D MHD

Ideal MHD equations are being added to the code.

2D Orszag-Tang Vortex problem

Perfomed using a 512 x 512 grid.

https://github.com/user-attachments/assets/9074f9f5-c4a3-4380-b70e-d01abafd4e05