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 |