Physical flux calculation along x direction for Euler equation.
JZ: here, we divide by a polynomial, we should be careful! We are leaving the polynomial space here!
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | state(:,:) | The state in nodal space. Dimension is the number of vars, i.e. 5 for Euler |
||
real(kind=rk), | intent(in) | :: | isenCoeff | Adiabatice index, also known as isentropic expansion factor. |
||
real(kind=rk), | intent(in) | :: | penalty_char(nPoints) | The value of the characteristic function (stemming from penalization) |
||
real(kind=rk), | intent(in) | :: | porosity | The porosity at the current point |
||
integer, | intent(in) | :: | nPoints | number of points |
||
integer, | intent(in) | :: | rot(5) | rotation |
||
real(kind=rk), | intent(inout) | :: | physFlux(:,:) | The physical flux along the x axis for all variables |
||
real(kind=rk), | intent(in) | :: | U_o(nPoints) | The velocity of the obstacle |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
real(kind=rk), | private | :: | inv_s1 | ||||
real(kind=rk), | private | :: | s2 | ||||
real(kind=rk), | private | :: | s3 | ||||
real(kind=rk), | private | :: | s4 | ||||
real(kind=rk), | private | :: | s5 | ||||
real(kind=rk), | private | :: | s1 | ||||
real(kind=rk), | private | :: | pressure | ||||
real(kind=rk), | private | :: | vel1 | ||||
real(kind=rk), | private | :: | cPor | ||||
integer, | private | :: | iPoint |