Commit 9d20963f authored by Lluis Jofre Cruanyes's avatar Lluis Jofre Cruanyes
Browse files

Implementation of HES flux scheme II

parent abe3d454

src/FlowSolverRHEA.cpp

+2 −1
Original line number Diff line number Diff line
@@ -4713,7 +4713,8 @@ double HesFluxApproximateRiemannSolver::calculateIntercellFlux( const double &rh 
    double V_3_RRR = ( -1.0 )*rho_RRR*ds_drhoE_RRR;



    /// SOBOLEV-LIKE SENSOR: approximated derivatives

    double eps = 1e-12*pow( max( abs(V_1_L), max( abs(V_2_L), max( abs(V_3_L), max( abs(V_1_R), max( abs(V_2_R), abs(V_3_R) ) ) ) ) ), 2.0 );

    //double eps = 1e-12*pow( max( abs(V_1_L), max( abs(V_2_L), max( abs(V_3_L), max( abs(V_1_R), max( abs(V_2_R), abs(V_3_R) ) ) ) ) ), 2.0 ) + epsilon;

    double eps = 1e-12*pow( max( abs(V_1_L), max( abs(V_2_L), max( abs(V_3_L), max( abs(V_1_R), max( abs(V_2_R), abs(V_3_R) ) ) ) ) ), 2.0 ) + 1.0e-10;	

    double DV1_L = ( -V_1_RR + 8.0*V_1_R - 8.0*V_1_LL + V_1_LLL )/( 12.0*delta );

    double DV2_L = ( -V_2_RR + 8.0*V_2_R - 8.0*V_2_LL + V_2_LLL )/( 12.0*delta );

    double DV3_L = ( -V_3_RR + 8.0*V_3_R - 8.0*V_3_LL + V_3_LLL )/( 12.0*delta );