A radially symmetric Gaussian profile about r=0 spreads out radially.
For comparison, the 1d equations with radial source terms are solved in subdirectory 1drad on a fine grid.
Left column: Comparison of q(x,y=0,z=0) (blue line) and 1d solution (red line)
for the unsplit method (method(3) = 22)
Right column: Comparison of q(x,y=0,z=0) (blue line) and 1d solution (red line)
for the split method (method(3) = -22)
Several frames are shown around time t=1.3 when the ingoing shock bounces at the origin. This is not well resolved in this 3d computation.
The results above were obtained using CFL = 0.9 and show results only along one slice.
Below are some comparisons with CFL = 0.7 showing scatter plots of all the cell values of pressure vs. distance from origin. Again the left column is unsplit and the right column is split.
The unsplit method takes about 4 times longer to run (both under an hour on RJL's laptop).
