# RiemannSolver

**Repository Path**: li--ze/RiemannSolver

## Basic Information

- **Project Name**: RiemannSolver
- **Description**: No description available
- **Primary Language**: Unknown
- **License**: Not specified
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2024-12-05
- **Last Updated**: 2024-12-05

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README


# One-dimensional HLLC Riemann solver 

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. 

# Structure of the code
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.

# Compilation and running 
`cd Exec/Sod`   
`sh run_RiemannSolver.sh`

# Visualization
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.

# Results for the test cases
## 1. Sod shock tube test case

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   

$$
(\rho,u,p)=
\begin{cases}
(1.0, 0.0, 100000.0)~x\le0\\
(0.125, 0.0, 10000.0)~x\gt0\\
\end{cases}
$$

<img src="Exec/Sod/Results/SodShockTubeProblem.gif?raw=true&v=100" alt="your_alternative_text" width="100%" height="100%" loop="true" autoplay="true">

## Comparison with exact solution
<img src="Exec/Sod/Results/Comparison_Density.png" alt="Image" width="400"><img src="Exec/Sod/Results/Comparison_Velocity.png" alt="Image" width="400">  
<img src="Exec/Sod/Results/Comparison_Pressure.png" alt="Image" width="400">