CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

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

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 7, 2020, 11:17
Default Unable to set PIMPLE residual control data after upgrading from v 4.1 to v6.0
  #1
New Member
 
Alex Fells
Join Date: Sep 2019
Posts: 5
Rep Power: 7
alexfells is on a distinguished road
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 PIMPLE residual control not workin in OF. 6.0

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 is offline   Reply With Quote

Old   April 7, 2020, 14:03
Default
  #2
New Member
 
Alex Fells
Join Date: Sep 2019
Posts: 5
Rep Power: 7
alexfells is on a distinguished road
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
alexfells is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
laplacianFoam with source term Herwig OpenFOAM Running, Solving & CFD 17 November 19, 2019 14:47
HeatSource BC to the whole region in chtMultiRegionHeater xsa OpenFOAM Running, Solving & CFD 3 November 7, 2016 06:07
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13


All times are GMT -4. The time now is 02:36.