CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Unable to set PIMPLE residual control data after upgrading from v 4.1 to v6.0 (https://www.cfd-online.com/Forums/openfoam-solving/225779-unable-set-pimple-residual-control-data-after-upgrading-v-4-1-v6-0-a.html)

alexfells April 7, 2020 10:17

Unable to set PIMPLE residual control data after upgrading from v 4.1 to v6.0
 
Hi,

I am upgrading my simulation from OpenFOAM version 4.1 to version 6.0 and after doing so, I am now unable to control my residuals in fvSolution. I am aware that PIMPLE was upgraded for v6 and that the "residualControl" variable name has been changed to "outerCorrectorResidualControl" and I have updated my fvSolution file as per https://www.cfd-online.com/Forums/op...kin-6-0-a.html

Would anyone be able to point me in the direction of where I am going wrong? I'm sure it's something fairly simple but I would appreciate the help.

My fvSolution file is below:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "alpha.H2O.*"
    {
        nAlphaCorr      2;
        nAlphaSubCycles 4;
        cAlpha          4;

        MULESCorr      yes;
        nLimiterIter    3;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance      1e-8;
        relTol          0;
    }

    "(pcorr|pcorrFinal)"
    {
        solver          GAMG;
        smoother            DICGaussSeidel;
        tolerance      1e-5;
        relTol          0;
    }

    p_rgh
    {
        solver          GAMG;
        tolerance          1e-6;
        relTol              0;

        smoother            DICGaussSeidel;
        nPreSweeps          0;
        nPostSweeps        2;
            nFinestSweeps      2;
            minIter1                1;

        cacheAgglomeration      true;
        nCellsInCoarsestLevel  2588;
        agglomerator            faceAreaPair;
        mergeLevels            1;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance      1e-05;
        relTol          0.1;
    }

    UFinal
    {
        solver                smoothSolver;
        smoother        symGaussSeidel;
        tolerance        1e-06;
        relTol          0;
    }

    "(flm|fmm)"
    {
        solver PBiCG;
        preconditioner DILU;
        tolerance 1e-07;
        relTol 0;
    }

    "(flm|fmm)Final"
    {
        solver PBiCG;
        preconditioner DILU;
        tolerance 1e-08;
        relTol 0;
    }
}

PIMPLE
{
    momentumPredictor                no;
    nOuterCorrectors            10;
    nCorrectors                4;
    nNonOrthogonalCorrectors    2;
    //pRefCell                    0;
    //pRefValue                  0;

    outerCorrectorResidualControl
    {
        "(U|p_rgh)"
        {
            relTol          0;
            tolerance      0.001;
        }

        "(flm|fmm)"
        {
            relTol          0;
            tolerance      0.0001;
        }
    }
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Output from my parallel simulation:

Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: no residual control data found. Calculations will employ 10 corrector loops

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties


alexfells April 7, 2020 13:03

Issue solved and due to stupidity on my part.

I hadn't updated the HPC job submission script to use OFv6 and therefore, the simulation was being performed using OFv4.1. As I had changed the "residualControl" to "outerCorrectorResidualControl", PIMPLE was defaulting to 10 iterations.

Cheers,

Alex


All times are GMT -4. The time now is 13:59.