Tags

Updated figures with the new DT computation as well as the correct axis for the skewed shock problem

Now that DT is computed correctly, I realize that for the gresho vortex, we need CFL of 0.5

Usual cfl condition is computed with

 1/dt = cx/dx + cy/dy + cz/dz (loop over directions)

Here, we did

1/dt = cx/dx + cy/dy + cz/dz+cx/dx + cy/dy + cz/dz+cx/dx + cy/dy + cz/dz (loop over faces)

Resulting in dt twice too small

Test case paper result might change from now on. Recover by dividing cfl by two

see the corresponding folder in the repo for the ini files

This version of the code was used to generate the plots for the mhd paper.

Find the setting files and plotting  routines in test_case_for_mhd_paper/

Warning, to switch from standard blast to low beta blast, you have to comment/uncomment the corresponding lines in src/problems/BlastInitKernel.cpp