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

FOAM FATAL ERROR: Maximum number of iterations exceeded: 100

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 10, 2019, 11:15
Default FOAM FATAL ERROR: Maximum number of iterations exceeded: 100
  #1
New Member
 
Antonio
Join Date: Mar 2019
Posts: 6
Rep Power: 2
antoniomollo is on a distinguished road
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"
Attached Files
File Type: zip constant.zip (1.8 KB, 0 views)
File Type: zip system.zip (8.6 KB, 0 views)
antoniomollo is offline   Reply With Quote

Old   October 10, 2019, 14:44
Default
  #2
Senior Member
 
Join Date: Aug 2015
Posts: 313
Rep Power: 9
clapointe is on a distinguished road
You have massive min/max pressures (and negative pressure!) from the beginning. Suggests boundary/initial condition problem.

Caelan
clapointe is offline   Reply With Quote

Old   October 10, 2019, 15:08
Default
  #3
New Member
 
Antonio
Join Date: Mar 2019
Posts: 6
Rep Power: 2
antoniomollo is on a distinguished road
Quote:
Originally Posted by clapointe View Post
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 is offline   Reply With Quote

Old   October 10, 2019, 15:10
Default
  #4
New Member
 
Antonio
Join Date: Mar 2019
Posts: 6
Rep Power: 2
antoniomollo is on a distinguished road
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;
}
}


// ************************************************** *********************** //
antoniomollo is offline   Reply With Quote

Old   October 11, 2019, 00:22
Default
  #5
Senior Member
 
Join Date: Aug 2015
Posts: 313
Rep Power: 9
clapointe is on a distinguished road
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
clapointe is offline   Reply With Quote

Reply

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 On
Pingbacks are On
Refbacks are On


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 468 September 17, 2019 21:28
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


All times are GMT -4. The time now is 08:28.