Godunov flux for the 1D Euler equation.
Implement correct Riemann problem for the equation with porosity. Right now the plain Euler problem is solved and the porosity and penalty_char are just averaged centrally.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(atl_euler_type), | intent(in) | :: | euler | |||
real(kind=rk), | intent(in) | :: | state_left(:,:) | The state on the face from its left limit (in conservative variables). |
||
real(kind=rk), | intent(in) | :: | state_right(:,:) | The state on the face from its right limit (in conservative variables). |
||
real(kind=rk), | intent(in) | :: | material_left(:,:) | The left value of the characteristic function (stemming from penalization) |
||
real(kind=rk), | intent(in) | :: | material_right(:,:) | The right value of the characteristic function (stemming from penalization) |
||
integer, | intent(in) | :: | nPoints | Number of points to evaluate the flux at. |
||
real(kind=rk), | intent(out) | :: | flux(:,:) | Resulting flux for the left element (in conservative variables). |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
real(kind=rk), | private | :: | p_left | ||||
real(kind=rk), | private | :: | p_right | ||||
real(kind=rk), | private | :: | v_left | ||||
real(kind=rk), | private | :: | v_right | ||||
real(kind=rk), | private | :: | estate(3) | ||||
real(kind=rk), | private | :: | pstate(3) | ||||
real(kind=rk), | private | :: | isen_coeff | ||||
real(kind=rk), | private | :: | icm1 | ||||
real(kind=rk), | private | :: | dummy | ||||
real(kind=rk), | private | :: | dummy_2 | ||||
real(kind=rk), | private | :: | material(ubound(material_left,1)) | ||||
integer, | private | :: | iPoint | ||||
integer, | private | :: | matpoint | ||||
integer, | private | :: | mm |