atl_physFluxEuler_1d Function

public function atl_physFluxEuler_1d(state, isenCoeff, penalty_scaling, U_o) result(physFlux)

Physical flux calculation along x direction for Euler equation.

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in) :: state(:)

The state in nodal space. Dimension is the number of vars, i.e. 3 for Euler 1d

real(kind=rk), intent(in) :: isenCoeff

Adiabatic index, also known as isentropic expansion factor.

real(kind=rk), intent(in) :: penalty_scaling

The scaling of the mass flux due to the penality in the porous material.

It is given by $1 + (\frac{1}{\phi} - 1) \Chi$ Thus, if there is no porous media, this factor has to be 1.

real(kind=rk), intent(in) :: U_o

Obstacle velocity

Return Value real(kind=rk), (3)

The physical flux along the x axis for all variables