FOAM FATAL ERROR: Maximum number of iterations exceeded: 100

antoniomollo · October 10, 2019, 11:15

I am a new user of OpenFoam. I have this error, I tried different ways to avoid this but the result is the same. Can anyone help me?
I understood that there is a problem with energy equation...but I don't know.
Thanks.

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : rhoPimpleFoam
Date : Oct 10 2019
Time : 17:05:43
Host : "antoniom91-HP-Pavilion-x360-Convertible"
PID : 5472
I/O : uncollated
Case : /home/antoniom91/OpenFOAM/antoniom91-6/run/VFP_simulation2_rho_simple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: Corrector convergence criteria found
"(U|k|epsilon)": tolerance 0.0001, relTol 0
Calclations will do 50 corrections if the convergence criteria are not met

Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

pressureControl
pMax 1500000
pMin 100000

Creating turbulence model

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

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.0014163604 max: 5.065195

Starting time loop

Courant Number mean: 0.0013968052 max: 4.9952613
deltaT = 9.8619329e-06
Time = 9.86193e-06

PIMPLE: Iteration 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.01541721, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.023338802, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.021696637, No Iterations 1
smoothSolver: Solving for h, Initial residual = 1, Final residual = 0.002021401, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0082120498, No Iterations 2
GAMG: Solving for p, Initial residual = 0.031349804, Final residual = 0.001157596, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00086029742, global = -6.1886942e-05, cumulative = -6.1886942e-05
pressureControl: p max 2048360.3
GAMG: Solving for p, Initial residual = 0.35939903, Final residual = 0.0023679751, No Iterations 2
GAMG: Solving for p, Initial residual = 0.019107717, Final residual = 0.0011200636, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0019050911, global = 0.0014242547, cumulative = 0.0013623678
smoothSolver: Solving for epsilon, Initial residual = 0.7155894, Final residual = 0.031617312, No Iterations 1
bounding epsilon, min: -32.865838 max: 210676.04 average: 329.31726
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.050063332, No Iterations 1
PIMPLE: Iteration 2
smoothSolver: Solving for Ux, Initial residual = 0.1885659, Final residual = 0.00552603, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.43538988, Final residual = 0.0041218538, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.43413677, Final residual = 0.0042201254, No Iterations 1
smoothSolver: Solving for h, Initial residual = 0.74736897, Final residual = 0.0061303172, No Iterations 4
GAMG: Solving for p, Initial residual = 0.75047041, Final residual = 0.0038408842, No Iterations 1
GAMG: Solving for p, Initial residual = 0.001399869, Final residual = 0.001399869, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0023359741, global = -0.0012983479, cumulative = 6.4019828e-05
pressureControl: p max 37145843
pressureControl: p min -1102504.2
GAMG: Solving for p, Initial residual = 0.54157365, Final residual = 0.0045490989, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0031541605, Final residual = 0.0031541605, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0055287247, global = -0.0020528002, cumulative = -0.0019887804
pressureControl: p max 48090802
pressureControl: p min -34043479
smoothSolver: Solving for epsilon, Initial residual = 0.10537911, Final residual = 0.0010463339, No Iterations 1
bounding epsilon, min: -118502.74 max: 1.1528773e+09 average: 32663.937
smoothSolver: Solving for k, Initial residual = 0.85080698, Final residual = 0.018905995, No Iterations 2
PIMPLE: Iteration 3
smoothSolver: Solving for Ux, Initial residual = 0.14157198, Final residual = 0.0061435985, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.24628262, Final residual = 0.013605369, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.24371158, Final residual = 0.013440339, No Iterations 1
smoothSolver: Solving for h, Initial residual = 0.16678901, Final residual = 0.0014007156, No Iterations 3
GAMG: Solving for p, Initial residual = 0.43836041, Final residual = 0.0015686454, No Iterations 1
GAMG: Solving for p, Initial residual = 0.00058446854, Final residual = 0.00058446854, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00089785328, global = -0.00065931356, cumulative = -0.0026480939
pressureControl: p max 45989243
pressureControl: p min -17902714
GAMG: Solving for p, Initial residual = 0.38769917, Final residual = 0.0017128437, No Iterations 1
GAMG: Solving for p, Initial residual = 0.00071476724, Final residual = 0.00071476724, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0011236984, global = -0.00081905666, cumulative = -0.0034671506
pressureControl: p max 87261603
pressureControl: p min -31612707
smoothSolver: Solving for epsilon, Initial residual = 0.10539707, Final residual = 0.0041859928, No Iterations 1
bounding epsilon, min: -6.4329239e+08 max: 3.2489285e+11 average: 28509511
smoothSolver: Solving for k, Initial residual = 0.24679149, Final residual = 0.0048989557, No Iterations 2
PIMPLE: Iteration 4
smoothSolver: Solving for Ux, Initial residual = 0.99892642, Final residual = 0.020548135, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.99985405, Final residual = 0.088095772, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.9998084, Final residual = 0.015464463, No Iterations 2
smoothSolver: Solving for h, Initial residual = 1, Final residual = 0.0032430201, No Iterations 4

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded: 100

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam:

erfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:

er fectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 73.

FOAM aborting

#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hePsiThermo<Foam:

siThermo, Foam:

ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:

erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#3 Foam::hePsiThermo<Foam:

siThermo, Foam:

ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:

erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#4 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"

clapointe · October 10, 2019, 14:44

You have massive min/max pressures (and negative pressure!) from the beginning. Suggests boundary/initial condition problem.

Caelan

antoniomollo · October 10, 2019, 15:08

Quote:

Originally Posted by clapointe

You have massive min/max pressures (and negative pressure!) from the beginning. Suggests boundary/initial condition problem.

Caelan

I have corrected the boundary. But the error is the same.

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : rhoPimpleFoam
Date : Oct 10 2019
Time : 21:05:53
Host : "antoniom91-HP-Pavilion-x360-Convertible"
PID : 5296
I/O : uncollated
Case : /home/antoniom91/OpenFOAM/antoniom91-6/run/VFP_simulation2_rho_simple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: Corrector convergence criteria found
"(U|k|epsilon)": tolerance 0.0001, relTol 0
Calclations will do 50 corrections if the convergence criteria are not met

Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

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

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.0014163604 max: 5.065195

Starting time loop

Courant Number mean: 0.0013968052 max: 4.9952613
deltaT = 9.8619329e-06
Time = 9.86193e-06

PIMPLE: Iteration 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0067606, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.023338802, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.021696637, No Iterations 1
smoothSolver: Solving for h, Initial residual = 1, Final residual = 0.002021401, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0087385529, No Iterations 3
GAMG: Solving for p, Initial residual = 0.04134101, Final residual = 0.0024747188, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.8173259e-05, global = -3.2038111e-06, cumulative = -3.2038111e-06
GAMG: Solving for p, Initial residual = 0.43076641, Final residual = 0.004968437, No Iterations 3
GAMG: Solving for p, Initial residual = 0.034180016, Final residual = 0.0020686053, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.92591e-05, global = -2.765535e-06, cumulative = -5.9693461e-06
smoothSolver: Solving for epsilon, Initial residual = 0.02463015, Final residual = 0.00069668943, No Iterations 1
bounding epsilon, min: -32.879214 max: 25055.446 average: 267.28537
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.043250315, No Iterations 1
PIMPLE: Iteration 2
smoothSolver: Solving for Ux, Initial residual = 0.72034491, Final residual = 7.5697148e-05, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.41292388, Final residual = 4.2416415e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.40803199, Final residual = 4.3650929e-05, No Iterations 1
smoothSolver: Solving for h, Initial residual = 0.1605451, Final residual = 0.00087375001, No Iterations 3
GAMG: Solving for p, Initial residual = 0.3454394, Final residual = 0.0042882372, No Iterations 2
GAMG: Solving for p, Initial residual = 0.003520348, Final residual = 0.003520348, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.685274e-05, global = 1.2610556e-05, cumulative = 6.64121e-06
GAMG: Solving for p, Initial residual = 0.28431697, Final residual = 0.0037378878, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0029046327, Final residual = 0.0029046327, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0001278336, global = 1.9418362e-05, cumulative = 2.6059572e-05
smoothSolver: Solving for epsilon, Initial residual = 0.02724506, Final residual = 4.5704602e-06, No Iterations 1
bounding epsilon, min: -64190.597 max: 6224817.7 average: 1134.893
smoothSolver: Solving for k, Initial residual = 0.2591801, Final residual = 0.00015350397, No Iterations 1
PIMPLE: Iteration 3
smoothSolver: Solving for Ux, Initial residual = 0.99999963, Final residual = 0.032678231, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.99997201, Final residual = 0.022285855, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.99999693, Final residual = 0.058668532, No Iterations 1
smoothSolver: Solving for h, Initial residual = 1, Final residual = 0.0032136039, No Iterations 3

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded: 100

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam:

erfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:

er fectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 73.

FOAM aborting

#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hePsiThermo<Foam:

siThermo, Foam:

ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:

erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#3 Foam::hePsiThermo<Foam:

siThermo, Foam:

ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:

erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#4 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"

antoniomollo · October 10, 2019, 15:10

These are the boundary:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
iniettore_1
{
type calculated;
value uniform 0;
}
iniettore_3
{
type calculated;
value uniform 0;
}
lato_N2
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
iniettore_2
{
type calculated;
value uniform 0;
}
iniettore_4
{
type calculated;
value uniform 0;
}
braccio_iniettori
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
lato_ugello
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
ugello
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 2;

boundaryField
{
iniettore_1
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.02;
value uniform 2;
}
iniettore_3
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.02;
value uniform 2;
}
lato_N2
{
type epsilonWallFunction;
value uniform 2;
}
outlet
{
type zeroGradient;
}
iniettore_2
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.02;
value uniform 2;
}
iniettore_4
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.02;
value uniform 2;
}
braccio_iniettori
{
type epsilonWallFunction;
value uniform 2;
}
lato_ugello
{
type epsilonWallFunction;
value uniform 2;
}
ugello
{
type epsilonWallFunction;
value uniform 2;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.375;

boundaryField
{
iniettore_1
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.375;
}
iniettore_3
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.375;
}
lato_N2
{
type kqRWallFunction;
value uniform 0.375;
}
outlet
{
type zeroGradient;
}
iniettore_2
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.375;
}
iniettore_4
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.375;
}
braccio_iniettori
{
type kqRWallFunction;
value uniform 0.375;
}
lato_ugello
{
type kqRWallFunction;
value uniform 0.375;
}
ugello
{
type kqRWallFunction;
value uniform 0.375;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
iniettore_1
{
type calculated;
value uniform 0;
}
iniettore_3
{
type calculated;
value uniform 0;
}
lato_N2
{
type nutkWallFunction;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
iniettore_2
{
type calculated;
value uniform 0;
}
iniettore_4
{
type calculated;
value uniform 0;
}
braccio_iniettori
{
type nutkWallFunction;
value uniform 0;
}
lato_ugello
{
type nutkWallFunction;
value uniform 0;
}
ugello
{
type nutkWallFunction;
value uniform 0;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
iniettore_1
{
type zeroGradient;
}
iniettore_3
{
type zeroGradient;

}
lato_N2
{
type zeroGradient;
}
outlet
{
type waveTransmissive;
rho rho;
gamma 1.4;
fieldInf 101325;
lInf 10;
value uniform 101325;

}
iniettore_2
{
type zeroGradient;

}
iniettore_4
{
type zeroGradient;

}
braccio_iniettori
{
type zeroGradient;
}
lato_ugello
{
type zeroGradient;
}
ugello
{
type zeroGradient;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
iniettore_1
{
type zeroGradient;

}
iniettore_3
{
type zeroGradient;

}
lato_N2
{
type fixedValue;
value uniform 820;
}
outlet
{
type zeroGradient;

}
iniettore_2
{
type zeroGradient;
}
iniettore_4
{
type zeroGradient;

}
braccio_iniettori
{
type zeroGradient;
}
lato_ugello
{
type fixedValue;
value uniform 820;
}
ugello
{
type zeroGradient;
}
}

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

/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
iniettore_1
{
type fixedValue;
value uniform (0 -300 0);
}
iniettore_3
{
type fixedValue;
value uniform (0 0 300);
}
lato_N2
{
type noSlip;
}
outlet
{
type zeroGradient;
}
iniettore_2
{
type fixedValue;
value uniform (0 300 0);
}
iniettore_4
{
type fixedValue;
value uniform (0 0 -300);
}
braccio_iniettori
{
type noSlip;
}
lato_ugello
{
type noSlip;
}
ugello
{
type noSlip;
}
}

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

clapointe · October 11, 2019, 00:22

Just a quick glance, but it looks like you tried to specify velocity at 4 patches, which I presume are inlets/outlets, and pressure at none of them. This could be why the simulation is blowing up. I don't know which solver you're using, but I would take a look at the tutorials for examples of how to set pressure/velocity at inlets/outlets.

Caelan

HeloïseM · March 2, 2023, 06:13

Hello, I have a similar problem i think but all seems okay regarding boundary conditions and setFieldsDict. I am running a specific mesh with the ddtFoam solver to modelize deflagration to detonation combustion in a rotational detonation engine.
The checkMesh indicates that everything is okay with it.

I have this error:

Acoustic Courant Number mean: 1.12966e-05 max: 0.317766
deltaT = 4.29456e-11
Time = 0.000055444344, time step = 680
Quenching 0.170678 % of the turbulent reaction
Calculating Riemann fluxes, second order = true
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for Ux, Initial residual = 3.43464e-05, Final residual = 2.94524e-15, No Iterations 1
DICPCG: Solving for Uy, Initial residual = 0.00309516, Final residual = 1.94156e-15, No Iterations 1
DICPCG: Solving for Uz, Initial residual = 0.00869599, Final residual = 5.231e-17, No Iterations 1
DICPCG: Solving for fH, Initial residual = 1.51402e-07, Final residual = 1.51402e-07, No Iterations 0
diagonal: Solving for tauHI, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for tauLO, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for c, Initial residual = 6.34503e-06, Final residual = 1.37472e-16, No Iterations 1
Combustion progress = 0.784048%
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoEu, Initial residual = 0, Final residual = 0, No Iterations 0
[8]
[8]
[8] --> FOAM FATAL ERROR:
[8] Maximum number of iterations exceeded
Test = 10280.7 K,
eu_ = 1.19965e+07 J/kg, e(Test,fH) = 1.18617e+07 J/kg
[8]
[8] From function eUnburnedThermo::correct(const volScalarField& fH)
[8] in file eUnburnedThermo/eUnburnedThermo.C at line 413.
[8]
FOAM parallel run aborting
[8]
[8] #0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
[8] #1 Foam::error::abort() at ??:?
[8] #2 Foam::eUnburnedThermo::correct() at ??:?
[8] #3
[8] at ??:?
[8] #4 __libc_start_main in "/lib64/libc.so.6"
[8] #5
[8] at ??:?
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 8 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

Could someone try to help me please ?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend"	bigphil	OpenFOAM CC Toolkits for Fluid-Structure Interaction	686	December 22, 2022 09:10
Free surface issues with interDyMFoam for hydroturbine	oumnion	OpenFOAM Running, Solving & CFD	0	October 6, 2017 14:05
SLTS+rhoPisoFoam: what is rDeltaT???	nileshjrane	OpenFOAM Running, Solving & CFD	4	February 25, 2013 04:13
pisoFoam with k-epsilon turb blows up - Some questions	Heroic	OpenFOAM Running, Solving & CFD	26	December 17, 2012 03:34
formatted point data- icem	HSK	CFX	12	August 11, 2011 21:25

October 10, 2019, 14:44		#2
clapointe Senior Member Join Date: Aug 2015 Posts: 494 Rep Power: 14	You have massive min/max pressures (and negative pressure!) from the beginning. Suggests boundary/initial condition problem. Caelan

October 10, 2019, 15:10		#4
antoniomollo New Member Antonio Join Date: Mar 2019 Posts: 8 Rep Power: 7	These are the boundary: /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object alphat; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { iniettore_1 { type calculated; value uniform 0; } iniettore_3 { type calculated; value uniform 0; } lato_N2 { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } outlet { type calculated; value uniform 0; } iniettore_2 { type calculated; value uniform 0; } iniettore_4 { type calculated; value uniform 0; } braccio_iniettori { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } lato_ugello { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } ugello { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object epsilon; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 2; boundaryField { iniettore_1 { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.02; value uniform 2; } iniettore_3 { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.02; value uniform 2; } lato_N2 { type epsilonWallFunction; value uniform 2; } outlet { type zeroGradient; } iniettore_2 { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.02; value uniform 2; } iniettore_4 { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.02; value uniform 2; } braccio_iniettori { type epsilonWallFunction; value uniform 2; } lato_ugello { type epsilonWallFunction; value uniform 2; } ugello { type epsilonWallFunction; value uniform 2; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.375; boundaryField { iniettore_1 { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.375; } iniettore_3 { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.375; } lato_N2 { type kqRWallFunction; value uniform 0.375; } outlet { type zeroGradient; } iniettore_2 { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.375; } iniettore_4 { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.375; } braccio_iniettori { type kqRWallFunction; value uniform 0.375; } lato_ugello { type kqRWallFunction; value uniform 0.375; } ugello { type kqRWallFunction; value uniform 0.375; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { iniettore_1 { type calculated; value uniform 0; } iniettore_3 { type calculated; value uniform 0; } lato_N2 { type nutkWallFunction; value uniform 0; } outlet { type calculated; value uniform 0; } iniettore_2 { type calculated; value uniform 0; } iniettore_4 { type calculated; value uniform 0; } braccio_iniettori { type nutkWallFunction; value uniform 0; } lato_ugello { type nutkWallFunction; value uniform 0; } ugello { type nutkWallFunction; value uniform 0; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 101325; boundaryField { iniettore_1 { type zeroGradient; } iniettore_3 { type zeroGradient; } lato_N2 { type zeroGradient; } outlet { type waveTransmissive; rho rho; gamma 1.4; fieldInf 101325; lInf 10; value uniform 101325; } iniettore_2 { type zeroGradient; } iniettore_4 { type zeroGradient; } braccio_iniettori { type zeroGradient; } lato_ugello { type zeroGradient; } ugello { type zeroGradient; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { iniettore_1 { type zeroGradient; } iniettore_3 { type zeroGradient; } lato_N2 { type fixedValue; value uniform 820; } outlet { type zeroGradient; } iniettore_2 { type zeroGradient; } iniettore_4 { type zeroGradient; } braccio_iniettori { type zeroGradient; } lato_ugello { type fixedValue; value uniform 820; } ugello { type zeroGradient; } } // ************************************************ ********************* // /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 6 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { iniettore_1 { type fixedValue; value uniform (0 -300 0); } iniettore_3 { type fixedValue; value uniform (0 0 300); } lato_N2 { type noSlip; } outlet { type zeroGradient; } iniettore_2 { type fixedValue; value uniform (0 300 0); } iniettore_4 { type fixedValue; value uniform (0 0 -300); } braccio_iniettori { type noSlip; } lato_ugello { type noSlip; } ugello { type noSlip; } } // ************************************************ ********************* //

October 11, 2019, 00:22		#5
clapointe Senior Member Join Date: Aug 2015 Posts: 494 Rep Power: 14	Just a quick glance, but it looks like you tried to specify velocity at 4 patches, which I presume are inlets/outlets, and pressure at none of them. This could be why the simulation is blowing up. I don't know which solver you're using, but I would take a look at the tutorials for examples of how to set pressure/velocity at inlets/outlets. Caelan

March 2, 2023, 06:13		#6
HeloïseM New Member Héloïse Magliano Join Date: Feb 2023 Posts: 4 Rep Power: 3	Hello, I have a similar problem i think but all seems okay regarding boundary conditions and setFieldsDict. I am running a specific mesh with the ddtFoam solver to modelize deflagration to detonation combustion in a rotational detonation engine. The checkMesh indicates that everything is okay with it. I have this error: Acoustic Courant Number mean: 1.12966e-05 max: 0.317766 deltaT = 4.29456e-11 Time = 0.000055444344, time step = 680 Quenching 0.170678 % of the turbulent reaction Calculating Riemann fluxes, second order = true diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DICPCG: Solving for Ux, Initial residual = 3.43464e-05, Final residual = 2.94524e-15, No Iterations 1 DICPCG: Solving for Uy, Initial residual = 0.00309516, Final residual = 1.94156e-15, No Iterations 1 DICPCG: Solving for Uz, Initial residual = 0.00869599, Final residual = 5.231e-17, No Iterations 1 DICPCG: Solving for fH, Initial residual = 1.51402e-07, Final residual = 1.51402e-07, No Iterations 0 diagonal: Solving for tauHI, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for tauLO, Initial residual = 0, Final residual = 0, No Iterations 0 DICPCG: Solving for c, Initial residual = 6.34503e-06, Final residual = 1.37472e-16, No Iterations 1 Combustion progress = 0.784048% diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoEu, Initial residual = 0, Final residual = 0, No Iterations 0 [8] [8] [8] --> FOAM FATAL ERROR: [8] Maximum number of iterations exceeded Test = 10280.7 K, eu_ = 1.19965e+07 J/kg, e(Test,fH) = 1.18617e+07 J/kg [8] [8] From function eUnburnedThermo::correct(const volScalarField& fH) [8] in file eUnburnedThermo/eUnburnedThermo.C at line 413. [8] FOAM parallel run aborting [8] [8] #0 Foam::error:rintStack(Foam::Ostream&) at ??:? [8] #1 Foam::error::abort() at ??:? [8] #2 Foam::eUnburnedThermo::correct() at ??:? [8] #3 [8] at ??:? [8] #4 __libc_start_main in "/lib64/libc.so.6" [8] #5 [8] at ??:? -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 8 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. Could someone try to help me please ?