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

outerCorrectorResidualControl OpenFOAM v1812 vs. OpenFOAM 6

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 27, 2019, 07:53
Default outerCorrectorResidualControl OpenFOAM v1812 vs. OpenFOAM 6
  #1
New Member
 
Korbinian Faust
Join Date: Feb 2019
Posts: 7
Rep Power: 7
Kobi_Faust is on a distinguished road
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:
Code:
PIMPLE
{
    nOuterCorrectors 50;
    nCorrectors     1;
    consistent      true;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
}
system/fluid1/fvSolution:
Code:
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:
Code:
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
Kobi_Faust is offline   Reply With Quote

Old   February 28, 2019, 04:33
Default
  #2
Senior Member
 
anonymous
Join Date: Jan 2016
Posts: 416
Rep Power: 14
simrego is on a distinguished road
Hi!


It is only residualControl, not outerCorrectorResidualControl
simrego is offline   Reply With Quote

Old   February 28, 2019, 05:47
Default
  #3
New Member
 
Korbinian Faust
Join Date: Feb 2019
Posts: 7
Rep Power: 7
Kobi_Faust is on a distinguished road
Quote:
Originally Posted by simrego View Post
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:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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?
Kobi_Faust is offline   Reply With Quote

Old   February 28, 2019, 06:11
Default
  #4
Senior Member
 
anonymous
Join Date: Jan 2016
Posts: 416
Rep Power: 14
simrego is on a distinguished road
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 is offline   Reply With Quote

Old   February 28, 2019, 06:19
Default
  #5
Senior Member
 
anonymous
Join Date: Jan 2016
Posts: 416
Rep Power: 14
simrego is on a distinguished road
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)
simrego is offline   Reply With Quote

Old   February 28, 2019, 08:51
Default
  #6
New Member
 
Korbinian Faust
Join Date: Feb 2019
Posts: 7
Rep Power: 7
Kobi_Faust is on a distinguished road
Quote:
Originally Posted by simrego View Post
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!!
Kobi_Faust is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, residual control

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Frequently Asked Questions about Installing OpenFOAM wyldckat OpenFOAM Installation 3 November 14, 2023 12:58
OpenFOAM Foundation releases OpenFOAMŪ 3.0.0 CFDFoundation OpenFOAM Announcements from OpenFOAM Foundation 1 November 7, 2015 16:16
OpenFOAM Foundation Releases OpenFOAM v2.3.0 opencfd OpenFOAM Announcements from OpenFOAM Foundation 3 December 23, 2014 04:43
[Gmsh] 2D Mesh Generation Tutorial for GMSH aeroslacker OpenFOAM Meshing & Mesh Conversion 12 January 19, 2012 04:52
64bitrhel5 OF installation instructions mirko OpenFOAM Installation 2 August 12, 2008 19:07


All times are GMT -4. The time now is 22:09.