Sfx_TestRoeFlux.h

Go to the documentation of this file.

 
#pragma once
 
namespace dbpp 
{
    // E. McNeil test (didn't test yet)
//  class TestRoeFlux : public BaseNumTreatmemt 
//  {
//  public:
        
double CalculMatriceRoe1D (double FR, double FL, double UR, double UL)
{
    double A;
 
    if (UR =! UL)
        A = (FR-FL)/(UR-UL);
    else if (UR == UL)
        A = 1.;    //  A vérifier
 
    return (A);
}
 
 
void TraitementTermesFluxConvectionPressionRoe1D (double *F1, double *F2,
                                                  double *U1, double *U2,
                                                  double *dU1, double  *dU2,
                                                  int NbSections, double B)
{
    //int j;
 
    /*double UL1, UL2, UR1, UR2;
    double FL1, FL2, FR1, FR2, FL2C, FR2C; 
    double A1, A2;*/
 
        //  Calculs préalables: évaluation des dU
    dU1[0] = minmod (U1[1]-U1[0], 0.);
    dU2[0] = minmod (U2[1]-U2[0], 0.);
 
    for (auto j = 1; j < NbSections; j++)
    {
        dU1[j] = minmod (U1[j+1]-U1[j], U1[j]-U1[j-1]);
        dU2[j] = minmod (U2[j+1]-U2[j], U2[j]-U2[j-1]);
    }
 
    for (auto j = 0; j < NbSections-1; j++)
    {
        //  Calcul des Fj+1/2 
        //  (Fj-1/2 représentant simplement le Fj+1/2 associé à la section précédente)
        //  Calcul basé sur Numerical scheme II de Nujic, 1995
        //  Fj+1/2 = 0.5 * (FR + FL - |Aj+1/2|(UR - UL)
 
        auto UL1 = U1[j] + 0.5*dU1[j];
        auto UR1 = U1[j+1] - 0.5*dU1[j+1];
 
        auto UL2 = U2[j] + 0.5*dU2[j];
        auto UR2 = U2[j+1] - 0.5*dU2[j+1];
 
        auto FL1 = UL2;
        auto FR1 = UR2;
 
        auto FL2 = EvaluationF2_C_1D (UL1, UL2, B);
        auto FR2 = EvaluationF2_C_1D (UR1, UR2, B);
 
        auto FL2C = EvaluationF2_I_1D (UL1, UL2, B);
        auto FR2C = EvaluationF2_I_1D (UR1, UR2, B);
 
        //fprintf (FichierDEBUG, "U1[%d]: %lf    U2[%d]: %lf\n", j, U1[1], j, U2[j]);
 
        //fprintf (FichierDEBUG, "FL1[%d]: %lf    FL2[%d]: %lf\n", j, FL1, j, FL2);
        //fprintf (FichierDEBUG, "FR1[%d]: %lf    FR2[%d]: %lf\n", j, FR1, j, FR2);
 
        //  Calcul de la matrice de Roe
        //  basé sur Alcrudo et Garcia-Navarro (1993)
 
        auto A1 = CalculMatriceRoe1D (FR1, FL1, UR1, UL1);
        auto A2 = CalculMatriceRoe1D (FR2C, FL2C, UR2, UL2);
 
        F1[j] = 0.5 * (FR1 + FL1 - (A1 * (UR1 - UL1)));
        F2[j] = 0.5 * (FR2 + FL2 - (A2 * (UR2 - UL2)));
 
        //fprintf (FichierDEBUG, "FF1[%d]: %lf    FF2[%d]: %lf\n\n", j, F1[j], j, F2[j]);
    }
}
 
//  };
} // End of namespace