This is a one-dimensional Riemann solver based on the Harten-Lax-Van Leer-Contact (HLLC) scheme with a minmod limiter for limiting the slopes of the primitive variables.
The source code is in Source
, and the test cases are in Exec
. In each of the test cases in
Exec
- for eg. Exec/Sod
, there is a file Initialize.cpp
in which the input parameters, and
initial condition for the simulation can be specified. There is a Makefile in each of the test case
directories.
cd Exec/Sod
sh run_RiemannSolver.sh
There is a Results
directory in each test case directory. It contains a Python script
PlotSolution.py
. python PlotSolution.py
will display a movie of the density as a
function of time.
The domain is -10 < x < 10
with 401 cells in the domain, and a time step of dt=1.25e-5 s
. The
end time is t=0.01 s
. The initial condition is