The 1D convergence tests consist of single Riemann problem with a known solution, the standard Sod shock. The solution on 100, 200, 400 and 800 cell resolution. The L1 density norm over the whole grid and in 3 defined sub-regions are compared.
The code and configuration files plus output for Fyris to run the 1D Riemann Sod problem will be available soon.
The Sod test
Specified as follows:
- adiabatic index: Gamma = 1.4
- grid domain: 0.0 < x < 1.0
- all cells with central coordinates < x0 are set to the Left State, all cells with central coordinates > x0 are set to the Right State (x0 lies on an exact cell boundary)
- left and right boundaries are free or natural boundaries that replicate the adjacent grid values
- tests are run from time = 0 until time = T = 0.2
The models are run and compared with the exact solution.
The exact solutions are integrated over each cell to get cell average values at the same resolution as the test models.
For the exact solutions the cells were subsampled by a factor of 5 and averaged.
L1 norms between the models and the integrated exact solutions are computed as percentages.
All L1 norms were performed on the density.
The exact solutions are included the riemann.c download as a set of text files.
All Fyris Alpha 1D tests were run with a Courant number of 0.8, with an initial step with a Courant number of 0.64 (0.8x0.8).
Four regions are defined for L1 norm analysis:
Region 0 is the whole 0.0 < x < 1.0
Region 1 is the interval 0.20 < x < 0.30 the head of the R wave
Region 2 is 0.49 < x < 0.68, between the CD and the R wave
Region 3 is 0.69 < x < 0.85, between the CD and the S wave
Plots comparing the exact and computed models follow. Finite models are shown as yellow dots, exact solutions as red line and red square dots.
Symbols used in the region plots, sod100 green, sod800 yellow, exact solution red line only.