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 |