[vinod956]

Hi friends,
I am very new to Openfoam. I copied a tutorial from combustion XiFoam then modified its geometry as per my project. After that I had run blockMesh. It was successful and polymesh was created. Then tried to run XiFoam solver on it and ended up with the error that no residual control data found. This is how it is shown
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2006 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _b45f8f6f58-20200629 OPENFOAM=2006
Arch : "LSB;label=32;scalar=64"
Exec : XiFoam
Date : Nov 19 2020
Time : 11:21:17
Host : LAPTOP-PE3G8S0K
PID : 434
I/O : uncollated
Case : /home/chitti/OpenFOAM_files/sim3b
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


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

Reading combustion properties

Found ignition cells:

14
(
0
1
2
40
41
42
80
81
82
120
121
160
161
200
)


Ignition on
Reading thermophysical properties

Selecting thermodynamics package
{
type heheuPsiThermo;
mixture homogeneousMixture;
transport const;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy absoluteEnthalpy;
}

min(b) = 1

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model LaunderSharmaKE
RAS
{
RASModel LaunderSharmaKE;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating field Xi

Creating the unstrained laminar flame speed

Selecting laminar flame speed correlation Gulders
Reading strained laminar flame speed field Su

Calculating turbulent flame speed field St

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Time = 5e-06

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.999996, Final residual = 0.000344645, No Iterations 1
Max St-Courant Number = 0.00432165
Igniting cell 0 state : 1 1 0.434 1e-15
Igniting cell 1 state : 1 1 0.434 1e-15
Igniting cell 2 state : 1 1 0.434 1e-15
Igniting cell 40 state : 1 1 0.434 1e-15
Igniting cell 41 state : 1 1 0.434 1e-15
Igniting cell 42 state : 1 1 0.434 1e-15
Igniting cell 80 state : 1 1 0.434 1e-15
Igniting cell 81 state : 1 1 0.434 1e-15
Igniting cell 82 state : 1 1 0.434 1e-15
Igniting cell 120 state : 1 1 0.434 1e-15
Igniting cell 121 state : 1 1 0.434 1e-15
Igniting cell 160 state : 1 1 0.434 4.24516e-13
Igniting cell 161 state : 1 1 0.434 1e-15
Igniting cell 200 state : 1 1 0.434 4.24516e-13
DILUPBiCGStab: Solving for b, Initial residual = 1, Final residual = 0.000382049, No Iterations 1
min(b) = 0.985393
I am also pasting my fvsolution file here
solvers
{
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}

"(p|rho)Final"
{
$p;
tolerance 1e-06;
relTol 0;
}

"(U|b|Su|Xi|ha|hau|k|epsilon)"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}

"(U|b|Su|Xi|ha|hau|k|epsilon)Final"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
}

PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
Can someone please help me in finding the error

[tas38]

You are missing a residual control block from within the PIMPLE block. You need something like the following if you plan to set "nOuterCorrectors" to a large value

PIMPLE
{
    nOuterCorrectors    50;    
    nCorrectors       2;
    nNonOrthogonalCorrectors 0;


    residualControl
    {
        U
        {
                tolerance  1e-6;
                relTol      0;
        }
        p
        {
                tolerance  1e-6;
                relTol      0;
        }
    }
    relaxationFactors
    {
        fields
        {
              p        0.3;
              pFinal   1;
        }
        equations
        {
              "U|k|epsilon"          0.3;
              "(U|k|epsilon)Final"   1;
        }
    
    }


} // END PIMPLE

[vinod956]

Originally Posted by tas38

You are missing a residual control block from within the PIMPLE block. You need something like the following if you plan to set "nOuterCorrectors" to a large value

PIMPLE
{
    nOuterCorrectors    50;    
    nCorrectors       2;
    nNonOrthogonalCorrectors 0;


    residualControl
    {
        U
        {
                tolerance  1e-6;
                relTol      0;
        }
        p
        {
                tolerance  1e-6;
                relTol      0;
        }
    }
    relaxationFactors
    {
        fields
        {
              p        0.3;
              pFinal   1;
        }
        equations
        {
              "U|k|epsilon"          0.3;
              "(U|k|epsilon)Final"   1;
        }
    
    }


} // END PIMPLE

Thank you for your reply. I will try

[Yann]

Hello,

This message is not an error:

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

It is just here to tell you there is no residual control defined, hence the solver will just run for 2 loops no matter what is the residual.

AFAIK this is not a problem. If the solver crashes this is certainly not because you did not define the residual controls. The problem is somewhere else.

Regards,
Yann

[vinod956]

Originally Posted by Yann

Hello,

This message is not an error:

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

It is just here to tell you there is no residual control defined, hence the solver will just run for 2 loops no matter what is the residual.

AFAIK this is not a problem. If the solver crashes this is certainly not because you did not define the residual controls. The problem is somewhere else.

Regards,
Yann

Thanks for your reply
Can u please figure out where the problem is?

[Yann]

I cannot really help, I have never used XiFoam and there is no error displayed in the log on your first post.


Is there anything else written in your terminal when you run the solver?


Yann

[vinod956]

Originally Posted by Yann

I cannot really help, I have never used XiFoam and there is no error displayed in the log on your first post.


Is there anything else written in your terminal when you run the solver?


Yann

The problem is iterations are not running.
Do you know how to run iterations.

[Yann]

Does the solver just hang forever?

Can you post a full log of your run?


I have no experience with XiFoam so if the solver is just hanging I am not sure I can help you. Maybe XiFoam users can help.



Regards,
Yann

[vinod956]

Originally Posted by Yann

Does the solver just hang forever?

Can you post a full log of your run?


I have no experience with XiFoam so if the solver is just hanging I am not sure I can help you. Maybe XiFoam users can help.



Regards,
Yann

chitti@LAPTOP-PE3G8S0K:~/OpenFOAM_files/sim3b$ XiFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2006 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _b45f8f6f58-20200629 OPENFOAM=2006
Arch : "LSB;label=32;scalar=64"
Exec : XiFoam
Date : Nov 21 2020
Time : 08:33:11
Host : LAPTOP-PE3G8S0K
PID : 473
I/O : uncollated
Case : /home/chitti/OpenFOAM_files/sim3b
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


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

Reading combustion properties

Found ignition cells:

14
(
0
1
2
40
41
42
80
81
82
120
121
160
161
200
)


Ignition on
Reading thermophysical properties

Selecting thermodynamics package
{
type heheuPsiThermo;
mixture homogeneousMixture;
transport const;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy absoluteEnthalpy;
}

min(b) = 1

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model LaunderSharmaKE
RAS
{
RASModel LaunderSharmaKE;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating field Xi

Creating the unstrained laminar flame speed

Selecting laminar flame speed correlation Gulders
Reading strained laminar flame speed field Su

Calculating turbulent flame speed field St

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Time = 5e-06

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.999996, Final residual = 0.000344645, No Iterations 1
Max St-Courant Number = 0.00432165
Igniting cell 0 state : 1 1 0.434 1e-15
Igniting cell 1 state : 1 1 0.434 1e-15
Igniting cell 2 state : 1 1 0.434 1e-15
Igniting cell 40 state : 1 1 0.434 1e-15
Igniting cell 41 state : 1 1 0.434 1e-15
Igniting cell 42 state : 1 1 0.434 1e-15
Igniting cell 80 state : 1 1 0.434 1e-15
Igniting cell 81 state : 1 1 0.434 1e-15
Igniting cell 82 state : 1 1 0.434 1e-15
Igniting cell 120 state : 1 1 0.434 1e-15
Igniting cell 121 state : 1 1 0.434 1e-15
Igniting cell 160 state : 1 1 0.434 4.24516e-13
Igniting cell 161 state : 1 1 0.434 1e-15
Igniting cell 200 state : 1 1 0.434 4.24516e-13
DILUPBiCGStab: Solving for b, Initial residual = 1, Final residual = 0.000382049, No Iterations 1
min(b) = 0.985393
chitti@LAPTOP-PE3G8S0K:~/OpenFOAM_files/sim3b$
This is the full log of my run. it is not hanging.