 |
 |
 |
|
|
|
[Sponsors] | ||||
OpenFOAM (interFoam) crashing while solving k-epsilon model. |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 1, 2021, 11:00
|
OpenFOAM (interFoam) crashing while solving k-epsilon model.
|
#1 |
|
New Member
Join Date: Jul 2021
Posts: 16
Rep Power: 5  |
Hello! I am new to OpenFOAM so this question might seem dumb. I want to simulate turbulent flow in an open channel. interFoam solves the case without any problems in laminar case but in turbulent condition it diverges. (high Courant Numbers). I get FPE error as something goes to infinity or division by zero.
I think it has something to do with boundary conditions in k, and epsilon but I could not figure it out. fvSchemes: Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,epsilon) Gauss linear;
div(phi,k) Gauss linear;
}
laplacianSchemes
{
default Gauss linear uncorrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"alpha.water.*"
{
nAlphaCorr 2;
nAlphaSubCycles 1;
cAlpha 1;
MULESCorr yes;
nLimiterIter 5;
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;
}
"pcorr.*"
{
solver PCG;
preconditioner DIC;
tolerance 1e-5;
relTol 0;
}
p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}
p_rghFinal
{
$p_rgh;
relTol 0;
}
U
{
solver smoothSolver; //
smoother symGaussSeidel;
relTol 0;
}
"(k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 1.0e-2;
maxIter 100;
};
}
PIMPLE
{
momentumPredictor no;
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 1;
}
relaxationFactors
{
fields
{
"(p)" 0.3;
}
equations
{
"(U)" 0.7;
"(k)" 0.3;
"(epsilon)" 0.8;
".*" 1;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
scale 0.1;
x1 4; // mainChannelLength
y1 2; // mainChannelWidth
z1 1; // mainChannelHeight
y2 3; // secondaryChannelLength
x2 1; // secondaryChannelWidth
z2 0.5; // secondaryChannelHeight
// Calculating edge coordinates
xx0 #eval{ 0.5*$x1+0.5*$x2 };
yy0 #eval{ $y1+$y2 };
xx1 #eval{ 0.5*$x1-0.5*$x2 };
zz0 #eval{ $z1-$z2 };
//MeshSiz 0.1; // Mesh Size in 2nd channel (spanwise)
//NcellS1 #eval{ std::floor( $x1/$MeshSiz )}; // Channel 1 number of cells (streamwise)
//NcellS2 #eval{ std::floor ($x2/$MeshSiz) }; // Channel 2 number of cells (spanwise)
vertices
(
//Block1
(0 0 0) // 0
($x1 0 0) // 1
($x1 $y1 0) // 2
(0 $y1 0) // 3
(0 0 $z1) // 4
($x1 0 $z1) // 5
($x1 $y1 $z1) // 6
(0 $y1 $z1) // 7
//Block2
($xx0 $y1 $zz0) //8
($xx0 $yy0 $zz0) //9
($xx1 $yy0 $zz0) //10
($xx1 $y1 $zz0) //11
($xx0 $y1 $z1) //12
($xx0 $yy0 $z1) //13
($xx1 $yy0 $z1) //14
($xx1 $y1 $z1) //15
);
blocks
(
hex (0 1 2 3 4 5 6 7) (320 160 80) simpleGrading (1 1 1) //Channel 1
hex (8 9 10 11 12 13 14 15) (400 80 40) simpleGrading (1 1 1) //Channel 2
);
edges
(
);
boundary
(
outlet
{
type patch;
faces
(
(6 2 1 5)
);
}
inlet
{
type patch;
faces
(
(0 4 7 3)
);
}
bottomwall
{
type wall;
faces
(
(0 1 2 3)
);
}
leftwall2
{
type wall;
faces
(
(8 12 15 11)
);
}
atmosphere
{
type patch;
faces
(
(4 5 6 7)
);
}
rightwall
{
type wall;
faces
(
(0 1 5 4)
);
}
leftwall
{
type wall;
faces
(
(2 6 7 3)
);
}
// Box2
c2bottom
{
type wall;
faces
(
(8 9 10 11)
);
}
c2atmosphere
{
type patch;
faces
(
(12 13 14 15)
);
}
c2outlet
{
type patch;
faces
(
(9 13 14 10)
);
}
c2rightwall
{
type wall;
faces
(
(8 9 13 12)
);
}
c2leftwall
{
type wall;
faces
(
(11 10 14 15)
);
}
);
mergePatchPairs
(
(leftwall leftwall2)
);
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application interFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 5;
deltaT 0.001;
writeControl adjustable;
writeInterval 0.05;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
adjustTimeStep yes;
maxCo 1;
maxAlphaCo 1;
maxDeltaT 1;
#sinclude "sampling"
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.1;
boundaryField
{
leftwall
{
type kqRWallFunction;
value uniform 0.1;
}
rightwall
{
type kqRWallFunction;
value uniform 0.1;
}
bottomwall
{
type kqRWallFunction;
value uniform 0.1;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
outlet
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.01;
value uniform 0.06;
}
c2bottom
{
type kqRWallFunction;
value uniform 0.1;
}
c2atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
c2outlet
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
c2rightwall
{
type kqRWallFunction;
value uniform 0.1;
}
c2leftwall
{
type kqRWallFunction;
value uniform 0.1;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 0.1;
boundaryField
{
leftwall
{
type epsilonWallFunction;
value uniform 0.1;
}
rightwall
{
type epsilonWallFunction;
value uniform 0.1;
}
bottomwall
{
type epsilonWallFunction;
value uniform 0.1;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
outlet
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.005;
value uniform 0.000055;
}
c2bottom
{
type epsilonWallFunction;
value uniform 0.1;
}
c2atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
c2outlet
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
c2rightwall
{
type epsilonWallFunction;
value uniform 0.1;
}
c2leftwall
{
type epsilonWallFunction;
value uniform 0.1;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
|
|
|
|
|
|
October 1, 2021, 11:10
|
|
#2 |
|
New Member
Join Date: Jul 2021
Posts: 16
Rep Power: 5 |
When I use adjustable time step, the time step goes to zero: Output here
Code:
PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.576479, Final residual = 1.69443e-09, No Iterations 3 Phase-1 volume fraction = 0.999932 Min(alpha.water) = 0 Max(alpha.water) = 1.1816 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.999932 Min(alpha.water) = 0 Max(alpha.water) = 1.1686 DICPCG: Solving for p_rgh, Initial residual = 0.99993, Final residual = 0.0219691, No Iterations 2 DICPCG: Solving for p_rgh, Initial residual = 2.78473e-38, Final residual = 2.78473e-38, No Iterations 0 time step continuity errors : sum local = 2.34799e+44, global = -8.45321e+43, cumulative = -8.46412e+43 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib64/libc.so.6 #3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #5 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:? #8 ? at ??:? #9 __libc_start_main in /lib64/libc.so.6 #10 ? at ??:? Code:
Starting time loop
--> FOAM Warning :
From function void Foam::particle::locate(const Foam::Vector<double> &, const Foam::Vector<double> *, int, bool, const Foam::string &)
in file particle/particle.C at line 484
Particle initialised with a location outside of the mesh when tracking from centre (0.0225 0.1975 0.0075) of cell 6324 to position (0.02 0.2 0.005)
Reading set description:
gauge_1
gauge_2
Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
Time = 0.001
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0480845, No Iterations 38
DICPCG: Solving for p_rgh, Initial residual = 0.001251, Final residual = 5.07562e-05, No Iterations 24
time step continuity errors : sum local = 8.21433e-05, global = -1.06508e-05, cumulative = -1.06508e-05
DICPCG: Solving for p_rgh, Initial residual = 0.0181754, Final residual = 0.000821557, No Iterations 38
DICPCG: Solving for p_rgh, Initial residual = 0.0023616, Final residual = 9.93711e-05, No Iterations 24
time step continuity errors : sum local = 5.54347e-05, global = 7.24914e-06, cumulative = -3.4017e-06
DICPCG: Solving for p_rgh, Initial residual = 0.00272227, Final residual = 0.00012654, No Iterations 17
DICPCG: Solving for p_rgh, Initial residual = 0.000127637, Final residual = 8.15941e-08, No Iterations 44
time step continuity errors : sum local = 4.43924e-08, global = -1.19515e-09, cumulative = -3.4029e-06
DILUPBiCG: Solving for epsilon, Initial residual = 0.432547, Final residual = 0.00151274, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.00501427, No Iterations 2
ExecutionTime = 2 s ClockTime = 2 s
Courant Number mean: 0.933855 max: 10.6067
Interface Courant Number mean: 0 max: 0
Time = 0.002
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1, Final residual = 6.31571e-09, No Iterations 3
Phase-1 volume fraction = 1 Min(alpha.water) = 0.999999 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1 Min(alpha.water) = 0.999999 Max(alpha.water) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.0257119, Final residual = 0.00127956, No Iterations 41
DICPCG: Solving for p_rgh, Initial residual = 0.00222235, Final residual = 9.88605e-05, No Iterations 29
time step continuity errors : sum local = 2.89055e-05, global = 4.45382e-06, cumulative = 1.05093e-06
DICPCG: Solving for p_rgh, Initial residual = 0.0199581, Final residual = 0.00090856, No Iterations 38
DICPCG: Solving for p_rgh, Initial residual = 0.00437537, Final residual = 0.000210576, No Iterations 28
time step continuity errors : sum local = 1.2806e-05, global = -1.03409e-06, cumulative = 1.68356e-08
DICPCG: Solving for p_rgh, Initial residual = 0.0306757, Final residual = 0.00131807, No Iterations 16
DICPCG: Solving for p_rgh, Initial residual = 0.00133985, Final residual = 9.60899e-08, No Iterations 71
time step continuity errors : sum local = 5.49774e-09, global = 2.29314e-10, cumulative = 1.70649e-08
DILUPBiCG: Solving for epsilon, Initial residual = 0.709524, Final residual = 0.00573007, No Iterations 6
bounding epsilon, min: -0.207824 max: 500.884 average: 9.90857
DILUPBiCG: Solving for k, Initial residual = 0.796152, Final residual = 0.00556545, No Iterations 6
bounding k, min: -58.4318 max: 204.832 average: 4.33867
ExecutionTime = 2.78 s ClockTime = 3 s
Courant Number mean: 0.927321 max: 8.9843
Interface Courant Number mean: 0 max: 0
Time = 0.003
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.99898, Final residual = 4.47903e-11, No Iterations 6
Phase-1 volume fraction = 0.999898 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.999887 Min(alpha.water) = -0.0125567 Max(alpha.water) = 1.00597
DICPCG: Solving for p_rgh, Initial residual = 0.106187, Final residual = 0.00498852, No Iterations 9
DICPCG: Solving for p_rgh, Initial residual = 0.00227057, Final residual = 9.66437e-05, No Iterations 40
time step continuity errors : sum local = 1.0627e-05, global = 1.00609e-06, cumulative = 1.02316e-06
DICPCG: Solving for p_rgh, Initial residual = 0.0188761, Final residual = 0.000923598, No Iterations 34
DICPCG: Solving for p_rgh, Initial residual = 0.000584912, Final residual = 2.85582e-05, No Iterations 27
time step continuity errors : sum local = 5.54231e-06, global = 4.1449e-07, cumulative = 1.43765e-06
DICPCG: Solving for p_rgh, Initial residual = 0.00464214, Final residual = 0.000228561, No Iterations 8
DICPCG: Solving for p_rgh, Initial residual = 0.000243947, Final residual = 9.92472e-08, No Iterations 64
time step continuity errors : sum local = 1.81879e-08, global = 3.21452e-10, cumulative = 1.43797e-06
DILUPBiCG: Solving for epsilon, Initial residual = 0.0321632, Final residual = 0.000107652, No Iterations 1
bounding epsilon, min: -5.00812 max: 311199 average: 135.697
DILUPBiCG: Solving for k, Initial residual = 0.997182, Final residual = 0.00755244, No Iterations 9
bounding k, min: -3.71915 max: 4283.07 average: 155.914
ExecutionTime = 3.49 s ClockTime = 3 s
Courant Number mean: 0.923043 max: 9.75704
Interface Courant Number mean: 0.00260689 max: 9.35186
Time = 0.004
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.491544, Final residual = 6.08724e-10, No Iterations 4
Phase-1 volume fraction = 0.999951 Min(alpha.water) = 0.582409 Max(alpha.water) = 1.00286
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.999956 Min(alpha.water) = 0.582409 Max(alpha.water) = 1.03078
DICPCG: Solving for p_rgh, Initial residual = 0.872382, Final residual = 0.0131008, No Iterations 2
DICPCG: Solving for p_rgh, Initial residual = 0.00083191, Final residual = 3.21955e-05, No Iterations 7
time step continuity errors : sum local = 0.000684782, global = 0.000135832, cumulative = 0.00013727
DICPCG: Solving for p_rgh, Initial residual = 2.76153e-05, Final residual = 1.21254e-06, No Iterations 43
DICPCG: Solving for p_rgh, Initial residual = 9.4024e-07, Final residual = 9.41337e-08, No Iterations 24
time step continuity errors : sum local = 1.07063e-05, global = 1.06688e-06, cumulative = 0.000138337
DICPCG: Solving for p_rgh, Initial residual = 1.0033e-05, Final residual = 4.62191e-07, No Iterations 43
DICPCG: Solving for p_rgh, Initial residual = 4.45772e-07, Final residual = 9.3462e-08, No Iterations 18
time step continuity errors : sum local = 1.10211e-05, global = -1.91586e-06, cumulative = 0.000136421
DILUPBiCG: Solving for epsilon, Initial residual = 0.129346, Final residual = 0.00122242, No Iterations 32
bounding epsilon, min: -731.285 max: 2.60303e+07 average: 947816
DILUPBiCG: Solving for k, Initial residual = 0.785097, Final residual = 0.00508704, No Iterations 26
bounding k, min: -1975.57 max: 20642.1 average: 3522.04
ExecutionTime = 4.21 s ClockTime = 4 s
Courant Number mean: 0.912524 max: 13.9602
Interface Courant Number mean: 0.00281132 max: 13.9602
Time = 0.005
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.538678, Final residual = 2.54492e-09, No Iterations 6
Phase-1 volume fraction = 0.999988 Min(alpha.water) = 0.81784 Max(alpha.water) = 1.0082
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.999989 Min(alpha.water) = 0.797185 Max(alpha.water) = 1.04048
DICPCG: Solving for p_rgh, Initial residual = 0.999924, Final residual = 0.0407437, No Iterations 4
DICPCG: Solving for p_rgh, Initial residual = 0.000732188, Final residual = 2.99942e-05, No Iterations 31
time step continuity errors : sum local = 3171.82, global = 98.9218, cumulative = 98.9219
DICPCG: Solving for p_rgh, Initial residual = 0.999999, Final residual = 0.0496443, No Iterations 20
DICPCG: Solving for p_rgh, Initial residual = 0.0137463, Final residual = 0.000649829, No Iterations 26
time step continuity errors : sum local = 1.77798e+11, global = -2.06326e+09, cumulative = -2.06326e+09
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0484649, No Iterations 20
DICPCG: Solving for p_rgh, Initial residual = 0.013821, Final residual = 9.27114e-08, No Iterations 80
time step continuity errors : sum local = 3.35538e+22, global = -1.09243e+21, cumulative = -1.09243e+21
DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 0.00460737, No Iterations 2
bounding epsilon, min: -3.82187e+43 max: 1.414e+44 average: 9.1937e+39
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.000613373, No Iterations 3
bounding k, min: -1.30274e+34 max: 3.10203e+34 average: 8.83677e+29
ExecutionTime = 4.84 s ClockTime = 5 s
Courant Number mean: 2.83494e+29 max: 3.80432e+32
Interface Courant Number mean: 6.81277e+22 max: 1.43842e+27
Time = 0.006
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.999997, Final residual = 2.32384e-09, No Iterations 22
Phase-1 volume fraction = 238364 Min(alpha.water) = 0 Max(alpha.water) = 1.81157e+10
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 4.1998e+46 Min(alpha.water) = -8.49254e+66 Max(alpha.water) = 1.07193e+67
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0340628, No Iterations 4
DICPCG: Solving for p_rgh, Initial residual = 3.21847e-84, Final residual = 3.21847e-84, No Iterations 0
time step continuity errors : sum local = 2.69341e+27, global = 1.21014e+27, cumulative = 1.21014e+27
DICPCG: Solving for p_rgh, Initial residual = 9.47547e-83, Final residual = 9.47547e-83, No Iterations 0
DICPCG: Solving for p_rgh, Initial residual = 9.47547e-83, Final residual = 9.47547e-83, No Iterations 0
time step continuity errors : sum local = 7.92964e+28, global = -3.68394e+27, cumulative = -2.47381e+27
DICPCG: Solving for p_rgh, Initial residual = 1.01454e-82, Final residual = 1.01454e-82, No Iterations 0
DICPCG: Solving for p_rgh, Initial residual = 1.01454e-82, Final residual = 1.01454e-82, No Iterations 0
time step continuity errors : sum local = 8.49029e+28, global = -6.65261e+27, cumulative = -9.12642e+27
DILUPBiCG: Solving for epsilon, Initial residual = 0.0817333, Final residual = 0.000213674, No Iterations 26
bounding epsilon, min: -2.62723e+66 max: 9.34649e+66 average: 1.28677e+65
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 801874, No Iterations 100
bounding k, min: -2.08023e+79 max: 2.32014e+80 average: 3.38051e+77
ExecutionTime = 5.6 s ClockTime = 6 s
Courant Number mean: 3.94164e+29 max: 2.82322e+32
Interface Courant Number mean: 2.39755e+27 max: 7.8528e+30
Time = 0.007
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1, Final residual = 8.55191e-09, No Iterations 30
Phase-1 volume fraction = -6.03994e+74 Min(alpha.water) = -4.49584e+79 Max(alpha.water) = 6.07195e+72
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib64/libc.so.6
#3 ? at ??:?
#4 ? at ??:?
#5 ? at ??:?
#6 __libc_start_main in /lib64/libc.so.6
#7 ? at ??:?
Floating point exception
|
|
|
|
|
|
|
| Tags |
| interfoam, kepsilonmodel, openfoam, turbulence |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| chtMultiRegionSimpleFoam: maximum number of iterations excedeed. | Nkl | OpenFOAM Running, Solving & CFD | 19 | October 10, 2019 03:42 |
| Segmentation fault when using reactingFOAM for Fluids | Tommy Floessner | OpenFOAM Running, Solving & CFD | 4 | April 22, 2018 13:30 |
| Cannot run the code properly: very large time step continuity error | crst15 | OpenFOAM Running, Solving & CFD | 9 | December 14, 2014 19:17 |
| Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 07:20 |
| pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 09:35 |
Linear Mode
