[agustinvo]

Hello,


I actually have some simulations around where I have a Reynolds stress model. It is similar to the "Speziale, Sarkar and Gatski Reynolds-stress" already implemented in OpenFOAM. Well, actually is an updated Manceau EB-RSM model (I took as a base the one this repo). These models have several explicit terms, whose sign I don't know a priori,

When running this simulation I face divergece, I see the pressure residuals exploding with the epsilon field.I use these configurations in fvSolution:

    p_rgh
    {
        solver          GAMG;
        smoother        DIC;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator faceAreaPair;
        mergeLevels 1;
        tolerance       1e-15;
        relTol 1e-6;
    }

    "(U|h|T|e)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-15;
        relTol          1e-6;
    }

    "(k|epsilon|kt|epsilont|R)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-15;
        relTol          1e-6;
    }

    f
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-15;
        relTol          1e-6;
    }

Have you ever simulated this kind of models before, or at least with a very explicit equation?


Regards,
Agustín