[Kobi_Faust]

Hi,
I am quite new to OpenFOAM so if any important Information is missing, just ask for it.

Solver: chtMultiRegionFoam
OF-Version: v1812

I have a Problem with the outerCorrectorResidualControl in PIMPLE Mode (Solver: chtMultiRegionFoam) when using OpenFOAM v1812 while on OpenFOAM 6 it works fine. (I use OpenFOAM 6 on my Computer but on our Cluster OpenFOAM v1812 is installed, so I can not change anything on this.)

I'm defining the general PIMPLE-loop in the fvSolution-file placed in the system folder. The outerCorrectorResidualConrol is given in the system/<material>/fvSolution files.

system/fvSolution:

PIMPLE
{
    nOuterCorrectors 50;
    nCorrectors     1;
    consistent      true;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
}

system/fluid1/fvSolution:

solvers
{
    rho
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0.1;
    }

    rhoFinal
    {
        $rho;
        tolerance       1e-7;
        relTol          0;
    }

    p_rgh
    {
        solver           GAMG;
        tolerance        1e-7;
        //relTol           0.01;
        smoother         DIC;
	maxIter		 20;
    }

    p_rghFinal
    {
        $p_rgh;
        tolerance        1e-7;
        relTol           0;
    }

    "(U|h|k|epsilon|R)"
    {
        solver           PBiCGStab;
        preconditioner   DILU;
        tolerance        1e-7;
        relTol           0.1;
    }

    "(U|h|k|epsilon|R)Final"
    {
        $U;
        tolerance        1e-7;
        relTol           0;
    }
}

PIMPLE
{
    
    momentumPredictor   true;
    nOuterCorrectors 50;
    nCorrectors     1;
    consistent      true;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
    outerCorrectorResidualControl
    {
        p_rgh
        {
            tolerance 1e-3;
            relTol 0;
        }
    }
}

relaxationFactors
{
    equations
    {
	"U.*"		0.5;
	"(h|e).*"	0.3;
	"(k|epsilon).*"	0.2;
    }
    fields
    {
    	"p_rgh.*"       0.7;
    }
}

system/solid1/fvSolution:

solvers
{
    h
    {
        solver           PCG;
        preconditioner   DIC;
        tolerance        1e-06;
        relTol           0.1;
    }

    hFinal
    {
        $h;
        tolerance        1e-06;
        relTol           0;
    }
}

PIMPLE
{
    
    momentumPredictor   true;
    nOuterCorrectors 50;
    nCorrectors     1;
    consistent      true;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
    outerCorrectorResidualControl
    {
        h
        {
            tolerance 1e-2;
            relTol 0;
        }
    }
}

When using OF 6 the outerCorrectors loop is left if all ResidualControl criteria are fulfilled (solid1 & fluid1). But when using OF v1812 the solver allways run the 50 OuterCorrector loops (probably without reading the ResidualControls).

What am I doing wrong, or how can I use ResiualControls in chtMultiRegionFoam with OF v1812? Any idea would be helpfull.

Thanks for your help
Kobi

[simrego]

Hi!


It is only residualControl, not outerCorrectorResidualControl

[Kobi_Faust]

Originally Posted by simrego

Hi!


It is only residualControl, not outerCorrectorResidualControl

Thanks for your answer.

I now changed the outerCorrectorResidualControl to residualControl in both regions (thats what you meant, or?), but it is still not used at all. It's again doing all 50 iterations at every timestep (even if all residuals are fulfilled).
When using OF 6 it already shows at the beginning of the output that residuals are given, but in OF v1812 I can see anything simillar.

Here the first lines of the output:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : b97cdf2-20190212 OPENFOAM=1812 patch=190129
Arch   : "LSB;label=32;scalar=64"
Exec   : chtMultiRegionFoam
Date   : Feb 28 2019
Time   : 09:47:47
Host   : xxx
PID    : 25553
I/O    : uncollated
Case   : /localhome/kfaust/Coupling_V3.0_steadyState
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region wasser for time = 0

Create solid mesh for region uran for time = 0

*** Reading fluid mesh thermophysical properties for region wasser

    Adding to thermoFluid

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to rhoFluid

    Adding to UFluid

    Adding to phiFluid

    Adding to hRefFluid

    Adding to ghFluid

    Adding to ghfFluid

    Adding to turbulenceFluid

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
    Adding to reactionFluid

Combustion model not active: combustionProperties not found
Selecting combustion model none
    Adding to radiationFluid

Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding to fieldsFluid

    Adding to QdotFluid

    Adding MRF

No MRF models present

    Adding fvOptions

*** Reading solid mesh thermophysical properties for region uran

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel none
    Adding fvOptions

Creating finite volume options from "system/fvOptions"

Selecting finite volume options type scalarSemiImplicitSource
    Source: energySource
    - selecting all cells
    - selected 13962 cell(s) with volume 0.094188321
Region: wasser Courant Number mean: 2906.2923 max: 9714.7911
Region: uran Diffusion Number mean: 0.0037731489 max: 0.035248857
Region: wasser Courant Number mean: 2906.2923 max: 9714.7911
Region: uran Diffusion Number mean: 0.0037731489 max: 0.035248857
Time = 1


Solving for fluid region wasser
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 0.053847496, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 1, Final residual = 0.05530738, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 1, Final residual = 0.054526391, No Iterations 1
DILUPBiCGStab:  Solving for h, Initial residual = 1, Final residual = 0.028893648, No Iterations 1
Min/max T:299.87069 321.74707
GAMG:  Solving for p_rgh, Initial residual = 0.99321625, Final residual = 7.8937817e-08, No Iterations 16
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (wasser): sum local = 2.5016758e-06, global = -5.4603265e-07, cumulative = -5.4603265e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.063409095, Final residual = 0.00058680414, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 0.020356198, No Iterations 1

Solving for solid region uran
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.00015761022, No Iterations 1
Min/max T:305.4462 400

Solving for fluid region wasser
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.42728706, Final residual = 0.028924828, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.42007585, Final residual = 0.027622375, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.58499379, Final residual = 0.054530451, No Iterations 1
DILUPBiCGStab:  Solving for h, Initial residual = 0.16005171, Final residual = 0.0038650269, No Iterations 1
Min/max T:299.88459 320.56578
GAMG:  Solving for p_rgh, Initial residual = 0.42074818, Final residual = 5.0234018e-08, No Iterations 16
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (wasser): sum local = 2.3008834e-06, global = 4.555841e-07, cumulative = -9.0448551e-08
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.2314172, Final residual = 0.0026828539, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.5568054, Final residual = 0.011547106, No Iterations 1

Solving for solid region uran
DICPCG:  Solving for h, Initial residual = 0.12704141, Final residual = 1.2759894e-05, No Iterations 1
Min/max T:302.69934 400.50223

<etc.>

Any further ideas, what I am doing wrong?

[simrego]

Hmm...


Do you have the original solver? Without modifications?
Because when the solver constructs the PIMPLE class it must tell you the settings.
From line 137:
https://openfoam.com/documentation/g...8C_source.html
You have to see that there is no residual control, OR the settings of the residual control OR you are using the PISO mode. But in your log I can't find any of it.
Not even then pimple loop number.

[simrego]

Oookay. So I checked the source code of the chtMultiRegionFoam, and they are not the same in the different versions.
If i'm correct in the v1812 version you always perform all of the outerCorrector loops. I don't see if the solver creates the PIMPLE class, it is only reads some information from the dictionary (max iter, flow only, turbOnFinalIt.., etc.).
Create a bug report (or maybe it is a feature request), for the esi group, or correct it yourself. (If you can on the cluster)

[Kobi_Faust]

Originally Posted by simrego

If i'm correct in the v1812 version you always perform all of the outerCorrector loops.

Yes, thats how it runs at the moment.


So I think I will write a bug report.

Thanks for your help!!