Commit 23f4fa6a authored by Lluis Jofre Cruanyes's avatar Lluis Jofre Cruanyes
Browse files

Implementation of HES flux scheme III

parent 9d20963f

src/FlowSolverRHEA.cpp

+4 −4
Original line number Diff line number Diff line
@@ -4609,8 +4609,8 @@ double EckepFluxApproximateRiemannSolver::calculateIntercellFlux( const double & 
        double psi_L = u_L*P_L/T_L;

        double psi_R = u_R*P_R/T_R;

        double deltaPsi = psi_R - psi_L;

        //double alpha_3  = ( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + epsilon );

        double alpha_3  = ( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + 1.0e-10 );	// ... modified for OpenACC

        //double alpha_3  = 2.0*( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + epsilon );

        double alpha_3  = 2.0*( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + 1.0e-10 );	// ... modified for OpenACC

        F = bar_F_3 - ( 1.0/2.0 )*alpha_3*deltaV_3;

    }
@@ -4809,8 +4809,8 @@ double HesFluxApproximateRiemannSolver::calculateIntercellFlux( const double &rh 
        double psi_L = u_L*P_L/T_L;

        double psi_R = u_R*P_R/T_R;

        double deltaPsi = psi_R - psi_L;

        //double alpha_3  = ( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + epsilon );

        double alpha_3  = ( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + 1.0e-10 );	// ... modified for OpenACC

        //double alpha_3  = 2.0*( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + epsilon );

        double alpha_3  = 2.0*( bar_F_1*deltaV_1 + bar_F_2*deltaV_2 + bar_F_3*deltaV_3 - deltaPsi )/( deltaV_3*deltaV_3 + 1.0e-10 );	// ... modified for OpenACC

        F = bar_F_3 - ( 1.0/2.0 )*alpha_3*deltaV_3;

        F -= ( 1.0/2.0 )*alpha_S*deltaU_3;

    }