Compute the exact solution of the Riemann problem for the Euler equation.
For the given state in primitive variables left and right, the solution is computed and stored in me. A Newton-Raphson iteration is used to find the solution in this non-linear system. Necessary parameters besides left and right state are the isentropic coefficient gam, and optionally a target tolerance to reach in the iterative scheme and the maximal number of iterations to perform.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(atl_ere_solState1D_type), | intent(out) | :: | me |
Description of the Riemann solution. If the input results in a critical state, me%critical will be set to true. nIter will contain the actual number of Newton-Raphson iterations, that were performed. |
||
| real(kind=rk), | intent(in) | :: | prim_left(3) |
Primitive variables left of discontinuity
|
||
| real(kind=rk), | intent(in) | :: | prim_right(3) |
Primitive variables right of discontinuity |
||
| real(kind=rk), | intent(in) | :: | gam |
Isentropic expansion coefficient of the gas |
||
| real(kind=rk), | intent(in), | optional | :: | tolerance |
Tolerated variation in the Newton-Raphson iterations, default: 8 epsilon |
|
| integer, | intent(in), | optional | :: | nMaxIter |
Maximal number of iterations to do, default: 100000 |