 |
 |
 |
|
|
|
[Sponsors] | ||||
October 15, 2015, 05:25
|
Floating point exception error
|
#1 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11  |
Hi everyone! I'm tring to simulate the flowfield on a stationary VAWT for some debugging. In the complete simulation (VAWT rotating) iterations start, but do not make sense. So I'm trying the same dictionaries with no rotations. Using pimpleFoam unstead of pimpleDyMFoam.
After running the first iteration for U, this appears: Code:
#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib64/libc.so.6" #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:? #6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? #10 at ??:? #11 __libc_start_main in "/lib64/libc.so.6" #12 at /home/abuild/rpmbuild/BUILD/glibc-2.19/csu/../sysdeps/x86_64/start.S:125 Floating point exception |
|
|
|
|
|
October 15, 2015, 05:40
|
|
#2 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
We need more information.
1) How does the rest of the log file look like? 2) Boundary conditions 3) fvSchemes 4) fvSolution Cheers, Philipp.
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 06:19
|
|
#3 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
rest of the log file
Code:
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare> of231
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare> pimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.1-262087cdf8db
Exec : pimpleFoam
Date : Oct 15 2015
Time : 11:26:00
Host : "ime042"
PID : 27789
Case : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}
No finite volume options present
PIMPLE: max iterations = 100
field "(p|U|k|omega)" : relTol 0, tolerance 0.01
Starting time loop
forces forces:
Not including porosity effects
Courant Number mean: 2.72708e-09 max: 2.97915e-07
Time = 1.09083e-05
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 8.13576e-14, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 6.4221e-14, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 6.59122e-14, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
...ETC...
Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; interpolate(HbyA) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; pcorr ; } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { pcorr { solver GAMG; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration off; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; tolerance 0.02; relTol 0; } p { $pcorr; tolerance 1e-07; relTol 0.01; } pFinal { $p; tolerance 1e-07; relTol 0; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; } UFinal { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; } "(k|omega)" { solver smoothSolver; smoother GaussSeidel; tolerance 1e-7; relTol 0.1; } "(k|omega)Final" { $U; tolerance 1e-7; relTol 0; } } PIMPLE { nNonOrthogonalCorrectors 0; correctPhi no; nOuterCorrectors 100; nCorrectors 2; pRefCell 0; pRefValue 0; residualControl { "(p|U|k|omega)" { tolerance 1e-2; relTol 0; } } } relaxationFactors { fields { } equations { "p.*" 0.3; "U.*" 0.7; "k.*" 0.7; "omega.*" 0.7; } } // ************************************************************************* // |
|
|
|
|
|
|
October 15, 2015, 06:22
|
|
#4 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
boundary
Code:
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
15
(
cubodx
{
type patch;
nFaces 592;
startFace 7406234;
}
cubosopra
{
type patch;
nFaces 588;
startFace 7406826;
}
cubosotto
{
type patch;
nFaces 588;
startFace 7407414;
}
cubosx
{
type patch;
nFaces 590;
startFace 7408002;
}
inlet
{
type patch;
nFaces 478;
startFace 7408592;
}
outlet
{
type patch;
nFaces 476;
startFace 7409070;
}
pala1
{
type wall;
nFaces 9506;
startFace 7409546;
}
pala1sopra
{
type wall;
nFaces 301;
startFace 7419052;
}
pala1sotto
{
type wall;
nFaces 301;
startFace 7419353;
}
pala2
{
type wall;
nFaces 9506;
startFace 7419654;
}
pala2sopra
{
type wall;
nFaces 301;
startFace 7429160;
}
pala2sotto
{
type wall;
nFaces 301;
startFace 7429461;
}
pala3
{
type wall;
nFaces 9506;
startFace 7429762;
}
pala3sopra
{
type wall;
nFaces 301;
startFace 7439268;
}
pala3sotto
{
type wall;
nFaces 301;
startFace 7439569;
}
)
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ 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
{
cubodx
{
type zeroGradient;
}
cubosopra
{
type zeroGradient;
}
cubosotto
{
type zeroGradient;
}
cubosx
{
type zeroGradient;
}
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.375;
}
outlet
{
type inletOutlet;
inletValue uniform 0.375;
value uniform 0.375;
}
pala1
{
type kqRWallFunction;
value uniform 0.375;
}
pala1sopra
{
type kqRWallFunction;
value uniform 0.375;
}
pala1sotto
{
type kqRWallFunction;
value uniform 0.375;
}
pala2
{
type kqRWallFunction;
value uniform 0.375;
}
pala2sopra
{
type kqRWallFunction;
value uniform 0.375;
}
pala2sotto
{
type kqRWallFunction;
value uniform 0.375;
}
pala3
{
type kqRWallFunction;
value uniform 0.375;
}
pala3sopra
{
type kqRWallFunction;
value uniform 0.375;
}
pala3sotto
{
type kqRWallFunction;
value uniform 0.375;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 1118.03;
boundaryField
{
cubodx
{
type zeroGradient;
}
cubosopra
{
type zeroGradient;
}
cubosotto
{
type zeroGradient;
}
cubosx
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform 1118.03;
}
outlet
{
type inletOutlet;
inletValue uniform 1118.03;
value uniform 1118.03;
}
pala1
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala1sopra
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala1sotto
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala2
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala2sopra
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala2sotto
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala3
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala3sopra
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
pala3sotto
{
type omegaWallFunction;
value uniform 1118.03;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
cubodx
{
type zeroGradient;
}
cubosopra
{
type zeroGradient;
}
cubosotto
{
type zeroGradient;
}
cubosx
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value $internalField;
}
pala1
{
type zeroGradient;
}
pala1sopra
{
type zeroGradient;
}
pala1sotto
{
type zeroGradient;
}
pala2
{
type zeroGradient;
}
pala2sopra
{
type zeroGradient;
}
pala2sotto
{
type zeroGradient;
}
pala3
{
type zeroGradient;
}
pala3sopra
{
type zeroGradient;
}
pala3sotto
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
cubodx
{
type slip;
}
cubosopra
{
type slip;
}
cubosotto
{
type slip;
}
cubosx
{
type slip;
}
inlet
{
type fixedValue;
value uniform (10.0 0 0);
}
outlet
{
type zeroGradient;
}
pala1
{
type fixedValue;
value $internalField;
}
pala1sopra
{
type fixedValue;
value $internalField;
}
pala1sotto
{
type fixedValue;
value $internalField;
}
pala2
{
type fixedValue;
value $internalField;
}
pala2sopra
{
type fixedValue;
value $internalField;
}
pala2sotto
{
type fixedValue;
value $internalField;
}
pala3
{
type fixedValue ;
value $internalField;
}
pala3sopra
{
type fixedValue;
value $internalField;
}
pala3sotto
{
type fixedValue;
value $internalField;
}
}
// ************************************************************************* //
|
|
|
|
|
|
|
October 15, 2015, 06:30
|
|
#5 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Michele, you use a wrong syntax for the relaxation factors in your fvSolution file.
Firstly, "p" must be written under "fields" instead of "equations" as far as I know. Thus, I guess your pressure field is not under-relaxated at all. Secondly, using the ".*" notation relaxes all iterations, also the last one, as someone stated to one of my posts here: http://www.cfd-online.com/Forums/ope...tml#post481037 So you should not use the ".*". I will write some more things if you have posted the rest of the log file.
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 06:32
|
|
#6 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
I have already posted the rest of the log file.
Thank you for the help. 
|
|
|
|
|
|
|
October 15, 2015, 07:11
|
|
#7 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Great, didn't see it.
The pressure solver doesn't work. Could you try to run it with some other setting? If the error remains we need to search somewhere else... Code:
"(p|pFinal)"
{
solver GAMG;
tolerance 1e-12;
relTol 0.1;
maxIter 20;
smoother DIC;
nPreSweeps 1;
nPostSweeps 1;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 07:38
|
|
#8 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
I've replaced the code you pasted in the precedent post with the precedent one. And changed all the settings in relaxation factors. My fvSolution file now looks like:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"(p|pFinal)"
{
solver GAMG;
tolerance 1e-12;
relTol 0.1;
maxIter 20;
smoother DIC;
nPreSweeps 1;
nPostSweeps 1;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}
p
{
$pcorr;
tolerance 1e-07;
relTol 0.01;
}
pFinal
{
$p;
tolerance 1e-07;
relTol 0;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
}
UFinal
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
}
"(k|omega)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-7;
relTol 0.1;
}
"(k|omega)Final"
{
$U;
tolerance 1e-7;
relTol 0;
}
}
PIMPLE
{
correctPhi no;
nOuterCorrectors 100;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
"(p|U|k|omega)"
{
tolerance 1e-2;
relTol 0;
}
}
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
omega 0.7;
}
}
// ************************************************************************* //
Trying to launch but it gives me Code:
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 8.13576e-14, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 6.4221e-14, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 6.59122e-14, No Iterations 1
--> FOAM FATAL IO ERROR:
keyword solver is undefined in dictionary "/home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum/system/fvSolution.solvers.p"
file: /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum/system/fvSolution.solvers.p from line 46 to line 47.
From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 437.
FOAM exiting
|
|
|
|
|
|
|
October 15, 2015, 07:43
|
|
#9 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Yes. You need to delete the two old entries "p{ ... }" and "pFinal{ ... }" directly below the new one. Otherwise "my" code is overwritten by your old entries and the solver keyword is missing.
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 09:53
|
|
#10 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
Well, this time it gets 7 iterations to create problems
error Code:
#0 Foam::error::printStack(Foam::Ostream&) at ??:? [3] #1 Foam::sigFpe::sigHandler(int) at ??:? [3] #2 in "/lib64/libc.so.6" [3] #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [3] #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [3] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [3] #6 [3] at ??:? [3] #7 [3] at ??:? [3] #8 [3] at ??:? [3] #9 [3] at ??:? [3] #10 __libc_start_main in "/lib64/libc.so.6" [3] #11 [3] at /home/abuild/rpmbuild/BUILD/glibc-2.19/csu/../sysdeps/x86_64/start.S:125 [ime042:32119] *** Process received signal *** [ime042:32119] Signal: Floating point exception (8) [ime042:32119] Signal code: (-6) [ime042:32119] Failing at address: 0x3e800007d77 [ime042:32119] [ 0] /lib64/libc.so.6(+0x35200) [0x7fe8b18e0200] [ime042:32119] [ 1] /lib64/libc.so.6(gsignal+0x37) [0x7fe8b18e0187] [ime042:32119] [ 2] /lib64/libc.so.6(+0x35200) [0x7fe8b18e0200] [ime042:32119] [ 3] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x369) [0x7fe8b2a17719] [ime042:32119] [ 4] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x2d) [0x7fe8b2a1795d] [ime042:32119] [ 5] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4e7) [0x7fe8b2a11397] [ime042:32119] [ 6] pimpleFoam() [0x439c5a] [ime042:32119] [ 7] pimpleFoam() [0x44e78f] [ime042:32119] [ 8] pimpleFoam() [0x44ea7e] [ime042:32119] [ 9] pimpleFoam() [0x42020d] [ime042:32119] [10] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7fe8b18ccb05] [ime042:32119] [11] pimpleFoam() [0x421071] [ime042:32119] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 3 with PID 32119 on node ime042 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.1-262087cdf8db
Exec : pimpleFoam -parallel
Date : Oct 15 2015
Time : 14:50:22
Host : "ime042"
PID : 32116
Case : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum
nProcs : 4
Slaves :
3
(
"ime042.32117"
"ime042.32118"
"ime042.32119"
)
Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}
No finite volume options present
PIMPLE: max iterations = 100
field "(p|U|k|omega)" : relTol 0, tolerance 0.01
Starting time loop
forces forces:
Not including porosity effects
Courant Number mean: 2.72708e-09 max: 2.97915e-07
Time = 1.09083e-05
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 8.10494e-14, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 6.42711e-14, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 6.55968e-14, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0941324, No Iterations 6
time step continuity errors : sum local = 5.13413e-10, global = -6.965e-12, cumulative = -6.965e-12
GAMG: Solving for p, Initial residual = 0.18443, Final residual = 0.0162853, No Iterations 2
time step continuity errors : sum local = 4.27511e-09, global = -4.4206e-10, cumulative = -4.49025e-10
PIMPLE: iteration 2
smoothSolver: Solving for Ux, Initial residual = 0.551614, Final residual = 0.000245009, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.117547, Final residual = 6.71312e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.130338, Final residual = 3.37683e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.155643, Final residual = 0.0111007, No Iterations 3
time step continuity errors : sum local = 2.67407e-09, global = -1.67831e-10, cumulative = -6.16856e-10
GAMG: Solving for p, Initial residual = 0.131887, Final residual = 0.00901081, No Iterations 2
time step continuity errors : sum local = 3.61971e-09, global = -3.24685e-10, cumulative = -9.41541e-10
PIMPLE: iteration 3
smoothSolver: Solving for Ux, Initial residual = 0.293694, Final residual = 0.000130662, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.0586409, Final residual = 3.56826e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.0660638, Final residual = 1.2281e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.119831, Final residual = 0.00790474, No Iterations 2
time step continuity errors : sum local = 3.18576e-09, global = -2.72153e-10, cumulative = -1.21369e-09
GAMG: Solving for p, Initial residual = 0.0971771, Final residual = 0.00894128, No Iterations 1
time step continuity errors : sum local = 4.45609e-09, global = -5.33984e-10, cumulative = -1.74768e-09
PIMPLE: iteration 4
smoothSolver: Solving for Ux, Initial residual = 0.0963746, Final residual = 3.0991e-05, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.0294886, Final residual = 1.47947e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.0324897, Final residual = 4.06102e-06, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0889663, Final residual = 0.00849765, No Iterations 1
time step continuity errors : sum local = 4.15247e-09, global = -4.6771e-10, cumulative = -2.21539e-09
GAMG: Solving for p, Initial residual = 0.0697397, Final residual = 0.00692439, No Iterations 1
time step continuity errors : sum local = 3.73189e-09, global = -4.04551e-10, cumulative = -2.61994e-09
PIMPLE: iteration 5
smoothSolver: Solving for Ux, Initial residual = 0.0491208, Final residual = 3.19173e-05, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.0162711, Final residual = 1.53472e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.0181073, Final residual = 4.73704e-06, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0667263, Final residual = 0.00600909, No Iterations 1
time step continuity errors : sum local = 3.28066e-09, global = -3.66243e-10, cumulative = -2.98618e-09
GAMG: Solving for p, Initial residual = 0.0517721, Final residual = 0.00505175, No Iterations 1
time step continuity errors : sum local = 2.93846e-09, global = -3.62277e-10, cumulative = -3.34846e-09
PIMPLE: iteration 6
smoothSolver: Solving for Ux, Initial residual = 0.015568, Final residual = 0.000222185, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00836338, Final residual = 4.93298e-05, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.0108277, Final residual = 7.2597e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0534809, Final residual = 0.00328122, No Iterations 3
time step continuity errors : sum local = 1.92493e-09, global = 6.20145e-11, cumulative = -3.28644e-09
GAMG: Solving for p, Initial residual = 0.0597616, Final residual = 0.00585955, No Iterations 1
time step continuity errors : sum local = 3.77776e-09, global = -2.90223e-11, cumulative = -3.31547e-09
PIMPLE: iteration 7
|
|
|
|
|
|
|
October 15, 2015, 09:57
|
|
#11 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Looks like it's the velocity solver this time, right?
I am not a solver expert, so I can't tell you why the smooth solver creates these problems. Usually, I use this setting: Code:
solvers
{
"(p|pFinal)"
{
solver GAMG;
tolerance 1e-12;
relTol 0.1;
maxIter 20;
smoother DIC;
nPreSweeps 1;
nPostSweeps 1;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 1.0e-2;
maxIter 100;
};
}
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 09:59
|
|
#12 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
Ok, so actually my fvSolution looks like
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"(p|pFinal)"
{
solver GAMG;
tolerance 1e-12;
relTol 0.1;
maxIter 20;
smoother DIC;
nPreSweeps 1;
nPostSweeps 1;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 1.0e-2;
maxIter 100;
};
}
PIMPLE
{
correctPhi no;
nOuterCorrectors 100;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
"(p|U|k|omega)"
{
tolerance 1e-2;
relTol 0;
}
}
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
omega 0.7;
}
}
// ************************************************************************* //
|
|
|
|
|
|
|
October 15, 2015, 10:02
|
|
#13 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
I would recommend a different residualControl setting. I usually use the relTol rather than the tolerance, because you don't have to concern yourself with the absolute values of the residuals. But maybe others think differently about that !?
Then, I would also lower the nOuterCorrectors to something like 20. If you can't get residuals down in 20 iterations, 100 won't work, too.
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 10:10
|
|
#14 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
Well, now it keeps going on.
First time step converged in 12 iterations. But numbers don't match with physics. Second time step starts with a max Courant number of 76, further more time step continuity errors keep growing. Code:
forces forces output:
sum of forces:
pressure : (3.20219e+13 6.22034e+11 -4.0116e+11)
viscous : (226981 3372.33 7927.69)
porous : (0 0 0)
sum of moments:
pressure : (-1.15048e+14 4.54921e+13 -3.75908e+14)
viscous : (-1.4342e+06 -2.72523e+06 -1.60083e+07)
porous : (0 0 0)
Courant Number mean: 4.59332e-07 max: 76.5398
Time = 2.18166e-05
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.866664, Final residual = 3.15156e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.257985, Final residual = 6.31871e-06, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.25924, Final residual = 4.27647e-07, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0638671, Final residual = 0.00387404, No Iterations 3
time step continuity errors : sum local = 2.82518e-09, global = 2.29057e-10, cumulative = -4.71539e-09
GAMG: Solving for p, Initial residual = 0.0570386, Final residual = 0.00406766, No Iterations 3
time step continuity errors : sum local = 2.34492e-09, global = 2.2376e-10, cumulative = -4.49163e-09
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.260287, Final residual = 2.3509e-05, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0472878, Final residual = 6.97051e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0472358, Final residual = 6.89489e-06, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0639765, Final residual = 0.00506699, No Iterations 3
time step continuity errors : sum local = 2.88782e-09, global = 2.81694e-10, cumulative = -4.20994e-09
GAMG: Solving for p, Initial residual = 0.0966392, Final residual = 0.00895581, No Iterations 2
time step continuity errors : sum local = 3.7894e-09, global = 4.5091e-10, cumulative = -3.75903e-09
PIMPLE: iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 0.0180231, Final residual = 4.31519e-05, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.00996661, Final residual = 1.28832e-05, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00989444, Final residual = 3.77719e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0819241, Final residual = 0.00574263, No Iterations 3
time step continuity errors : sum local = 2.42303e-09, global = 2.66958e-10, cumulative = -3.49207e-09
GAMG: Solving for p, Initial residual = 0.132618, Final residual = 0.0101418, No Iterations 2
time step continuity errors : sum local = 3.15306e-09, global = 4.19991e-10, cumulative = -3.07208e-09
PIMPLE: iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 0.0849585, Final residual = 0.000508934, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0167581, Final residual = 1.4555e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.0169335, Final residual = 8.71885e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.115883, Final residual = 0.0115837, No Iterations 1
time step continuity errors : sum local = 3.59584e-09, global = 1.59479e-09, cumulative = -1.47728e-09
GAMG: Solving for p, Initial residual = 0.0900661, Final residual = 0.00743689, No Iterations 3
time step continuity errors : sum local = 2.21426e-09, global = 2.69514e-10, cumulative = -1.20777e-09
PIMPLE: iteration 5
DILUPBiCG: Solving for Ux, Initial residual = 0.078644, Final residual = 0.000613499, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.0161631, Final residual = 0.00011581, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.016418, Final residual = 2.26945e-05, No Iterations 3
GAMG: Solving for p, Initial residual = 0.163533, Final residual = 0.015967, No Iterations 1
time step continuity errors : sum local = 3.72005e-09, global = 1.51401e-09, cumulative = 3.06234e-10
GAMG: Solving for p, Initial residual = 0.0970837, Final residual = 0.00890793, No Iterations 3
time step continuity errors : sum local = 2.17655e-09, global = 3.06431e-10, cumulative = 6.12665e-10
PIMPLE: iteration 6
DILUPBiCG: Solving for Ux, Initial residual = 0.0531039, Final residual = 2.68521e-05, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.00897067, Final residual = 7.4357e-05, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.0113846, Final residual = 2.54919e-05, No Iterations 3
GAMG: Solving for p, Initial residual = 0.18581, Final residual = 0.0142419, No Iterations 1
time step continuity errors : sum local = 2.69513e-09, global = 9.0678e-10, cumulative = 1.51945e-09
GAMG: Solving for p, Initial residual = 0.122813, Final residual = 0.0107259, No Iterations 2
time step continuity errors : sum local = 2.05469e-09, global = 5.09273e-10, cumulative = 2.02872e-09
PIMPLE: iteration 7
DILUPBiCG: Solving for Ux, Initial residual = 0.0142203, Final residual = 1.11577e-05, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00513482, Final residual = 7.10856e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00550852, Final residual = 5.4645e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 0.155244, Final residual = 0.0127537, No Iterations 1
time step continuity errors : sum local = 2.21653e-09, global = 7.59961e-10, cumulative = 2.78868e-09
GAMG: Solving for p, Initial residual = 0.105718, Final residual = 0.00990871, No Iterations 2
time step continuity errors : sum local = 1.74794e-09, global = 3.39822e-10, cumulative = 3.1285e-09
PIMPLE: iteration 8
DILUPBiCG: Solving for Ux, Initial residual = 0.00490551, Final residual = 1.22639e-05, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00136004, Final residual = 7.06376e-07, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00267391, Final residual = 5.58931e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.142433, Final residual = 0.00811321, No Iterations 3
time step continuity errors : sum local = 1.27769e-09, global = 4.30803e-12, cumulative = 3.13281e-09
GAMG: Solving for p, Initial residual = 0.169433, Final residual = 0.016526, No Iterations 2
time step continuity errors : sum local = 2.1341e-09, global = 2.21144e-10, cumulative = 3.35395e-09
PIMPLE: iteration 9
DILUPBiCG: Solving for Ux, Initial residual = 0.00287841, Final residual = 1.57161e-05, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.000483634, Final residual = 4.30534e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.00230007, Final residual = 1.06611e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.157906, Final residual = 0.0155495, No Iterations 3
time step continuity errors : sum local = 2.3094e-09, global = 5.37005e-10, cumulative = 3.89096e-09
GAMG: Solving for p, Initial residual = 0.141294, Final residual = 0.0131548, No Iterations 4
time step continuity errors : sum local = 2.03133e-09, global = 4.90192e-10, cumulative = 4.38115e-09
PIMPLE: iteration 10
DILUPBiCG: Solving for Ux, Initial residual = 0.000192267, Final residual = 2.26198e-07, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000125139, Final residual = 1.72955e-07, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.000595267, Final residual = 2.00749e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.165359, Final residual = 0.0152966, No Iterations 4
time step continuity errors : sum local = 2.71149e-09, global = 3.18975e-11, cumulative = 4.41305e-09
GAMG: Solving for p, Initial residual = 0.300512, Final residual = 0.0231589, No Iterations 3
time step continuity errors : sum local = 4.86893e-09, global = -7.42521e-11, cumulative = 4.3388e-09
PIMPLE: iteration 11
DILUPBiCG: Solving for Ux, Initial residual = 0.000350208, Final residual = 1.94364e-06, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00032242, Final residual = 2.32274e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00115122, Final residual = 2.25507e-06, No Iterations 4
GAMG: Solving for p, Initial residual = 0.348331, Final residual = 0.0293308, No Iterations 4
time step continuity errors : sum local = 7.32133e-09, global = -2.00533e-10, cumulative = 4.13826e-09
GAMG: Solving for p, Initial residual = 0.17185, Final residual = 0.0124721, No Iterations 3
time step continuity errors : sum local = 8.7529e-09, global = -1.76329e-10, cumulative = 3.96193e-09
PIMPLE: iteration 12
DILUPBiCG: Solving for Ux, Initial residual = 0.000558672, Final residual = 4.85112e-06, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.000562371, Final residual = 5.38573e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00113069, Final residual = 6.76897e-06, No Iterations 4
GAMG: Solving for p, Initial residual = 0.386677, Final residual = 0.034756, No Iterations 5
time step continuity errors : sum local = 1.65232e-08, global = 1.55496e-10, cumulative = 4.11743e-09
GAMG: Solving for p, Initial residual = 0.0460313, Final residual = 0.00430096, No Iterations 2
time step continuity errors : sum local = 5.36873e-08, global = -5.97013e-10, cumulative = 3.52042e-09
...
PIMPLE: iteration 25
DILUPBiCG: Solving for Ux, Initial residual = 0.45703, Final residual = 0.00364494, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.466817, Final residual = 0.0032379, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.448066, Final residual = 0.00264534, No Iterations 4
GAMG: Solving for p, Initial residual = 4.89754e-07, Final residual = 4.45106e-08, No Iterations 17
time step continuity errors : sum local = 0.0126292, global = -0.000375018, cumulative = -0.000298839
GAMG: Solving for p, Initial residual = 1.2146e-07, Final residual = 7.14817e-08, No Iterations 20
time step continuity errors : sum local = 0.128873, global = 0.00011411, cumulative = -0.00018473
PIMPLE: iteration 26
DILUPBiCG: Solving for Ux, Initial residual = 0.488462, Final residual = 0.00410571, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.489245, Final residual = 0.00398358, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.468887, Final residual = 0.00105651, No Iterations 4
GAMG: Solving for p, Initial residual = 1.55761e-07, Final residual = 1.28471e-08, No Iterations 6
time step continuity errors : sum local = 0.0261039, global = 0.000406571, cumulative = 0.000221841
GAMG: Solving for p, Initial residual = 4.17219e-08, Final residual = 3.32099e-09, No Iterations 12
time step continuity errors : sum local = 0.0568839, global = -0.00190224, cumulative = -0.0016804
PIMPLE: iteration 27
DILUPBiCG: Solving for Ux, Initial residual = 0.485182, Final residual = 0.00444568, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.477219, Final residual = 0.00236323, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.452158, Final residual = 0.000454943, No Iterations 5
GAMG: Solving for p, Initial residual = 9.29377e-08, Final residual = 7.82375e-09, No Iterations 7
time step continuity errors : sum local = 0.08866, global = 0.00500549, cumulative = 0.00332509
GAMG: Solving for p, Initial residual = 6.04943e-08, Final residual = 5.50121e-09, No Iterations 2
time step continuity errors : sum local = 0.142202, global = 0.00511167, cumulative = 0.00843676
PIMPLE: iteration 28
DILUPBiCG: Solving for Ux, Initial residual = 0.5799, Final residual = 0.0032216, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.56742, Final residual = 0.00354331, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.464652, Final residual = 0.00279237, No Iterations 3
GAMG: Solving for p, Initial residual = 1.15256e-07, Final residual = 1.73304e-08, No Iterations 20
time step continuity errors : sum local = 0.296138, global = -0.000460135, cumulative = 0.00797662
GAMG: Solving for p, Initial residual = 1.43725e-08, Final residual = 1.9104e-08, No Iterations 20
time step continuity errors : sum local = 3.256, global = -4.9317e-05, cumulative = 0.0079273
PIMPLE: iteration 29
DILUPBiCG: Solving for Ux, Initial residual = 0.475133, Final residual = 0.00420981, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.514787, Final residual = 0.00477123, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.387523, Final residual = 0.00241095, No Iterations 4
GAMG: Solving for p, Initial residual = 7.66632e-08, Final residual = 4.95751e-09, No Iterations 5
time step continuity errors : sum local = 0.351647, global = 0.00446064, cumulative = 0.0123879
GAMG: Solving for p, Initial residual = 1.11533e-08, Final residual = 8.46322e-10, No Iterations 4
time step continuity errors : sum local = 0.669685, global = -0.0141515, cumulative = -0.00176359
PIMPLE: iteration 30
DILUPBiCG: Solving for Ux, Initial residual = 0.482939, Final residual = 0.00169732, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.492853, Final residual = 0.00167008, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.423742, Final residual = 0.0033366, No Iterations 3
GAMG: Solving for p, Initial residual = 7.33705e-08, Final residual = 7.02836e-09, No Iterations 4
time step continuity errors : sum local = 1.208, global = 0.0817156, cumulative = 0.079952
GAMG: Solving for p, Initial residual = 9.10694e-09, Final residual = 7.2643e-10, No Iterations 8
time step continuity errors : sum local = 1.11649, global = 0.0696683, cumulative = 0.14962
PIMPLE: iteration 31
DILUPBiCG: Solving for Ux, Initial residual = 0.434401, Final residual = 0.00276104, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.438609, Final residual = 0.0035577, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.428942, Final residual = 0.00169325, No Iterations 5
GAMG: Solving for p, Initial residual = 6.13821e-06, Final residual = 8.82587e-07, No Iterations 20
time step continuity errors : sum local = 3.98246, global = 0.277155, cumulative = 0.426775
GAMG: Solving for p, Initial residual = 2.69328e-08, Final residual = 2.1687e-09, No Iterations 1
time step continuity errors : sum local = 4.15637, global = 0.0797123, cumulative = 0.506488
PIMPLE: iteration 32
DILUPBiCG: Solving for Ux, Initial residual = 0.728414, Final residual = 0.00296967, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.541076, Final residual = 0.00394594, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.79722, Final residual = 0.00145519, No Iterations 3
GAMG: Solving for p, Initial residual = 7.0245e-08, Final residual = 6.54451e-09, No Iterations 9
time step continuity errors : sum local = 5.78199, global = 0.245586, cumulative = 0.752073
GAMG: Solving for p, Initial residual = 5.14576e-09, Final residual = 3.40823e-10, No Iterations 8
time step continuity errors : sum local = 8.29289, global = -0.468393, cumulative = 0.283681
PIMPLE: iteration 33
DILUPBiCG: Solving for Ux, Initial residual = 0.431882, Final residual = 0.000428693, No Iterations 6
DILUPBiCG: Solving for Uy, Initial residual = 0.451638, Final residual = 0.00125437, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.500052, Final residual = 0.00311681, No Iterations 4
GAMG: Solving for p, Initial residual = 9.34601e-08, Final residual = 1.17991e-08, No Iterations 20
time step continuity errors : sum local = 23.4012, global = -15.7156, cumulative = -15.4319
GAMG: Solving for p, Initial residual = 4.27502e-08, Final residual = 1.25665e-08, No Iterations 20
time step continuity errors : sum local = 85.7935, global = -34.1987, cumulative = -49.6306
PIMPLE: iteration 34
DILUPBiCG: Solving for Ux, Initial residual = 0.40978, Final residual = 0.00235957, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.451869, Final residual = 0.001927, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.419105, Final residual = 0.00210072, No Iterations 5
GAMG: Solving for p, Initial residual = 1.92655e-06, Final residual = 8.91561e-07, No Iterations 20
time step continuity errors : sum local = 208.798, global = -51.2509, cumulative = -100.881
GAMG: Solving for p, Initial residual = 7.43198e-07, Final residual = 2.23723e-06, No Iterations 20
time step continuity errors : sum local = 1943.26, global = 18.5426, cumulative = -82.3388
PIMPLE: iteration 35
DILUPBiCG: Solving for Ux, Initial residual = 0.330767, Final residual = 0.0021958, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.343671, Final residual = 0.00117926, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.322184, Final residual = 0.00102286, No Iterations 4
GAMG: Solving for p, Initial residual = 2.22497e-07, Final residual = 2.12035e-08, No Iterations 6
time step continuity errors : sum local = 118.654, global = -61.5736, cumulative = -143.912
GAMG: Solving for p, Initial residual = 7.57818e-08, Final residual = 6.6613e-09, No Iterations 10
time step continuity errors : sum local = 186.871, global = 43.3242, cumulative = -100.588
PIMPLE: iteration 36
DILUPBiCG: Solving for Ux, Initial residual = 0.557896, Final residual = 0.00349763, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.566518, Final residual = 0.000803787, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.484968, Final residual = 0.00214561, No Iterations 4
GAMG: Solving for p, Initial residual = 3.62432e-07, Final residual = 3.45964e-08, No Iterations 13
time step continuity errors : sum local = 388.037, global = -4.52998, cumulative = -105.118
GAMG: Solving for p, Initial residual = 5.9511e-08, Final residual = 1.23765e-08, No Iterations 20
time step continuity errors : sum local = 1321.47, global = -0.315589, cumulative = -105.434
PIMPLE: iteration 37
DILUPBiCG: Solving for Ux, Initial residual = 0.486504, Final residual = 0.00191651, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.48672, Final residual = 0.00204066, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.427991, Final residual = 0.000811272, No Iterations 5
GAMG: Solving for p, Initial residual = 4.65018e-08, Final residual = 4.45178e-09, No Iterations 20
time step continuity errors : sum local = 731.468, global = 198.436, cumulative = 93.0026
GAMG: Solving for p, Initial residual = 9.32307e-09, Final residual = 4.49918e-09, No Iterations 20
time step continuity errors : sum local = 4870.59, global = -1325.02, cumulative = -1232.01
PIMPLE: iteration 38
DILUPBiCG: Solving for Ux, Initial residual = 0.479185, Final residual = 0.00165469, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.476365, Final residual = 0.000516798, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.489416, Final residual = 0.00472329, No Iterations 3
GAMG: Solving for p, Initial residual = 1.9266e-06, Final residual = 5.08132e-07, No Iterations 20
time step continuity errors : sum local = 7383.56, global = 3325.75, cumulative = 2093.74
GAMG: Solving for p, Initial residual = 4.86745e-07, Final residual = 4.49227e-07, No Iterations 20
time step continuity errors : sum local = 57164.3, global = -1297.7, cumulative = 796.039
PIMPLE: iteration 39
DILUPBiCG: Solving for Ux, Initial residual = 0.469658, Final residual = 0.00286888, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.486166, Final residual = 0.000287571, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.420831, Final residual = 0.00308034, No Iterations 3
|
|
|
|
|
|
|
October 15, 2015, 10:11
|
|
#15 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
Trying different settings like you suggested. I'll write on it.
|
|
|
|
|
|
|
October 15, 2015, 10:14
|
|
#16 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Don't cut the log file, I miss the most intersting parts
 What about the k and omaga solvers?
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 10:21
|
|
#17 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
this one, as you suggested, right?
Code:
"(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 1.0e-2;
maxIter 100;
};
|
|
|
|
|
|
|
October 15, 2015, 10:23
|
|
#18 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
No, I mean the log file is missing the information about the turbulence.
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
|
October 15, 2015, 10:27
|
|
#19 |
|
Member
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 11 |
As far as I remember, all the log files I've seen did not write nothing about k-omega. Did I forgot to write something?
Anyway, I tried with the following fvSolution Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"(p|pFinal)"
{
solver GAMG;
tolerance 0;
relTol 1e-5;
maxIter 20;
smoother DIC;
nPreSweeps 1;
nPostSweeps 1;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"
{
solver PBiCG;
preconditioner DILU;
tolerance 0;
relTol 1.0e-5;
maxIter 100;
};
}
PIMPLE
{
correctPhi no;
nOuterCorrectors 20;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
"(p|U|k|omega)"
{
tolerance 0;
relTol 1e-5;
}
}
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
omega 0.7;
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}
No finite volume options present
PIMPLE: max iterations = 20
field "(p|U|k|omega)" : relTol 1e-05, tolerance 0
Starting time loop
forces forces:
Not including porosity effects
Courant Number mean: 2.72708e-09 max: 2.97915e-07
Time = 1.09083e-05
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 5.38097e-17, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 2.19952e-16, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.98804e-17, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.000850538, No Iterations 20
time step continuity errors : sum local = 4.63897e-12, global = -2.99077e-14, cumulative = -2.99077e-14
GAMG: Solving for p, Initial residual = 0.182951, Final residual = 1.48404e-05, No Iterations 20
time step continuity errors : sum local = 3.92995e-12, global = -4.25516e-14, cumulative = -7.24593e-14
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.575869, Final residual = 5.20509e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.130436, Final residual = 5.97277e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.130436, Final residual = 6.17317e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 0.14681, Final residual = 1.26348e-05, No Iterations 20
time step continuity errors : sum local = 3.58e-12, global = -4.61894e-14, cumulative = -1.18649e-13
GAMG: Solving for p, Initial residual = 0.123814, Final residual = 6.48693e-06, No Iterations 20
time step continuity errors : sum local = 2.91471e-12, global = -4.44487e-14, cumulative = -1.63097e-13
PIMPLE: iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 0.328115, Final residual = 2.9441e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0652315, Final residual = 2.92304e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.0652764, Final residual = 3.10868e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 0.113768, Final residual = 5.22568e-06, No Iterations 20
time step continuity errors : sum local = 2.41215e-12, global = -3.75052e-14, cumulative = -2.00603e-13
GAMG: Solving for p, Initial residual = 0.091437, Final residual = 3.3779e-06, No Iterations 20
time step continuity errors : sum local = 1.89695e-12, global = -3.15321e-14, cumulative = -2.32135e-13
PIMPLE: iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 0.17629, Final residual = 1.56232e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0338782, Final residual = 1.51358e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.0339012, Final residual = 1.6177e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0863256, Final residual = 2.64111e-06, No Iterations 20
time step continuity errors : sum local = 1.50007e-12, global = -2.45599e-14, cumulative = -2.56695e-13
GAMG: Solving for p, Initial residual = 0.0670754, Final residual = 1.84298e-06, No Iterations 20
time step continuity errors : sum local = 1.15507e-12, global = -1.90806e-14, cumulative = -2.75775e-13
PIMPLE: iteration 5
DILUPBiCG: Solving for Ux, Initial residual = 0.0931013, Final residual = 8.08966e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0182765, Final residual = 7.92359e-11, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.0182847, Final residual = 8.65505e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0635045, Final residual = 1.43359e-06, No Iterations 20
time step continuity errors : sum local = 9.02351e-13, global = -1.44497e-14, cumulative = -2.90225e-13
GAMG: Solving for p, Initial residual = 0.048217, Final residual = 1.04543e-06, No Iterations 20
time step continuity errors : sum local = 6.94057e-13, global = -1.09189e-14, cumulative = -3.01144e-13
PIMPLE: iteration 6
DILUPBiCG: Solving for Ux, Initial residual = 0.048723, Final residual = 4.10896e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0105615, Final residual = 4.24163e-11, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.0105538, Final residual = 4.84504e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0454965, Final residual = 8.2809e-07, No Iterations 20
time step continuity errors : sum local = 5.50428e-13, global = -8.40717e-15, cumulative = -3.09551e-13
GAMG: Solving for p, Initial residual = 0.0339978, Final residual = 6.22983e-07, No Iterations 20
time step continuity errors : sum local = 4.26835e-13, global = -6.33486e-15, cumulative = -3.15886e-13
PIMPLE: iteration 7
DILUPBiCG: Solving for Ux, Initial residual = 0.025365, Final residual = 2.14062e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00664658, Final residual = 2.96478e-11, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00662996, Final residual = 2.91403e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0319445, Final residual = 5.04916e-07, No Iterations 20
time step continuity errors : sum local = 3.45932e-13, global = -4.91469e-15, cumulative = -3.208e-13
GAMG: Solving for p, Initial residual = 0.0236159, Final residual = 3.87828e-07, No Iterations 20
time step continuity errors : sum local = 2.70513e-13, global = -3.69509e-15, cumulative = -3.24496e-13
PIMPLE: iteration 8
DILUPBiCG: Solving for Ux, Initial residual = 0.0132529, Final residual = 3.12824e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00443574, Final residual = 1.94406e-11, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00442741, Final residual = 2.19916e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0221197, Final residual = 3.17427e-07, No Iterations 20
time step continuity errors : sum local = 2.21317e-13, global = -2.41478e-15, cumulative = -3.2691e-13
GAMG: Solving for p, Initial residual = 0.016238, Final residual = 2.47359e-07, No Iterations 20
time step continuity errors : sum local = 1.74362e-13, global = -1.75134e-15, cumulative = -3.28662e-13
PIMPLE: iteration 9
DILUPBiCG: Solving for Ux, Initial residual = 0.00708739, Final residual = 1.48741e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00301968, Final residual = 8.01464e-09, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.0030182, Final residual = 1.53244e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0151832, Final residual = 2.32595e-07, No Iterations 20
time step continuity errors : sum local = 1.6388e-13, global = -6.4121e-16, cumulative = -3.29303e-13
GAMG: Solving for p, Initial residual = 0.0110994, Final residual = 1.86506e-07, No Iterations 20
time step continuity errors : sum local = 1.32314e-13, global = -3.65334e-16, cumulative = -3.29668e-13
PIMPLE: iteration 10
DILUPBiCG: Solving for Ux, Initial residual = 0.00396473, Final residual = 9.72602e-12, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00205423, Final residual = 3.17736e-11, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.00205473, Final residual = 2.35815e-09, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0103786, Final residual = 2.70309e-07, No Iterations 20
time step continuity errors : sum local = 1.91695e-13, global = 6.29131e-17, cumulative = -3.29605e-13
GAMG: Solving for p, Initial residual = 0.00757311, Final residual = 2.43872e-07, No Iterations 20
time step continuity errors : sum local = 1.73711e-13, global = 1.07779e-16, cumulative = -3.29498e-13
PIMPLE: iteration 11
DILUPBiCG: Solving for Ux, Initial residual = 0.00234174, Final residual = 2.28446e-09, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00138961, Final residual = 1.31493e-09, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.00139063, Final residual = 2.06147e-09, No Iterations 3
GAMG: Solving for p, Initial residual = 0.00709553, Final residual = 5.57821e-07, No Iterations 20
time step continuity errors : sum local = 3.97233e-13, global = 2.16375e-16, cumulative = -3.29281e-13
GAMG: Solving for p, Initial residual = 0.00518617, Final residual = 5.17207e-07, No Iterations 20
time step continuity errors : sum local = 3.69376e-13, global = 1.81566e-16, cumulative = -3.291e-13
PIMPLE: iteration 12
DILUPBiCG: Solving for Ux, Initial residual = 0.001446, Final residual = 1.07055e-08, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.000934765, Final residual = 2.94952e-09, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.000936384, Final residual = 1.07785e-09, No Iterations 6
GAMG: Solving for p, Initial residual = 0.00488315, Final residual = 1.22038e-06, No Iterations 20
time step continuity errors : sum local = 8.71415e-13, global = 3.63095e-17, cumulative = -3.29063e-13
GAMG: Solving for p, Initial residual = 0.00361049, Final residual = 1.10981e-06, No Iterations 20
time step continuity errors : sum local = 7.93999e-13, global = 2.2947e-17, cumulative = -3.2904e-13
PIMPLE: iteration 13
DILUPBiCG: Solving for Ux, Initial residual = 0.000913855, Final residual = 8.92874e-10, No Iterations 7
DILUPBiCG: Solving for Uy, Initial residual = 0.000624009, Final residual = 6.93498e-10, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.000628673, Final residual = 2.40093e-09, No Iterations 7
GAMG: Solving for p, Initial residual = 0.00344592, Final residual = 2.92139e-06, No Iterations 20
time step continuity errors : sum local = 2.08991e-12, global = 1.37988e-17, cumulative = -3.29027e-13
GAMG: Solving for p, Initial residual = 0.00265492, Final residual = 2.466e-06, No Iterations 20
time step continuity errors : sum local = 1.76652e-12, global = -4.18087e-18, cumulative = -3.29031e-13
PIMPLE: iteration 14
DILUPBiCG: Solving for Ux, Initial residual = 0.000574472, Final residual = 2.83792e-09, No Iterations 7
DILUPBiCG: Solving for Uy, Initial residual = 0.000391953, Final residual = 8.47879e-10, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.000421002, Final residual = 1.50038e-09, No Iterations 8
GAMG: Solving for p, Initial residual = 0.00263674, Final residual = 5.59006e-06, No Iterations 20
time step continuity errors : sum local = 4.00474e-12, global = 2.33869e-17, cumulative = -3.29007e-13
GAMG: Solving for p, Initial residual = 0.00233292, Final residual = 3.80849e-06, No Iterations 20
time step continuity errors : sum local = 2.73149e-12, global = -1.58426e-17, cumulative = -3.29023e-13
PIMPLE: iteration 15
DILUPBiCG: Solving for Ux, Initial residual = 0.000263722, Final residual = 4.01094e-10, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 0.000157811, Final residual = 1.16427e-09, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.000278886, Final residual = 2.53124e-09, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00328214, Final residual = 1.48079e-05, No Iterations 20
time step continuity errors : sum local = 1.06313e-11, global = 1.10572e-16, cumulative = -3.28913e-13
GAMG: Solving for p, Initial residual = 0.00393972, Final residual = 1.76433e-05, No Iterations 20
time step continuity errors : sum local = 1.26889e-11, global = -9.98943e-17, cumulative = -3.29012e-13
PIMPLE: iteration 16
DILUPBiCG: Solving for Ux, Initial residual = 3.70449e-05, Final residual = 3.67244e-10, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 2.11178e-05, Final residual = 1.22973e-10, No Iterations 9
DILUPBiCG: Solving for Uz, Initial residual = 0.0001729, Final residual = 8.05336e-10, No Iterations 8
GAMG: Solving for p, Initial residual = 0.00536874, Final residual = 9.37444e-05, No Iterations 20
time step continuity errors : sum local = 6.75098e-11, global = -7.25852e-16, cumulative = -3.29738e-13
GAMG: Solving for p, Initial residual = 0.0132991, Final residual = 0.000106057, No Iterations 20
time step continuity errors : sum local = 7.74716e-11, global = 9.05626e-16, cumulative = -3.28833e-13
PIMPLE: iteration 17
DILUPBiCG: Solving for Ux, Initial residual = 0.000318289, Final residual = 3.27937e-10, No Iterations 12
DILUPBiCG: Solving for Uy, Initial residual = 0.000168046, Final residual = 1.47304e-09, No Iterations 8
DILUPBiCG: Solving for Uz, Initial residual = 2.59715e-05, Final residual = 8.5052e-11, No Iterations 10
GAMG: Solving for p, Initial residual = 0.0187933, Final residual = 0.00361403, No Iterations 20
time step continuity errors : sum local = 2.62536e-09, global = -2.1167e-13, cumulative = -5.40503e-13
GAMG: Solving for p, Initial residual = 0.0783631, Final residual = 0.00585178, No Iterations 20
time step continuity errors : sum local = 5.31193e-09, global = -4.36381e-13, cumulative = -9.76884e-13
PIMPLE: iteration 18
DILUPBiCG: Solving for Ux, Initial residual = 0.000485379, Final residual = 1.97349e-09, No Iterations 10
DILUPBiCG: Solving for Uy, Initial residual = 0.000632127, Final residual = 2.95446e-09, No Iterations 10
DILUPBiCG: Solving for Uz, Initial residual = 0.000310312, Final residual = 1.76645e-09, No Iterations 9
GAMG: Solving for p, Initial residual = 0.0984537, Final residual = 0.00645455, No Iterations 20
time step continuity errors : sum local = 4.64713e-09, global = 3.52369e-14, cumulative = -9.41647e-13
GAMG: Solving for p, Initial residual = 0.0469483, Final residual = 0.00256194, No Iterations 20
time step continuity errors : sum local = 6.73626e-09, global = -8.1758e-14, cumulative = -1.02341e-12
PIMPLE: iteration 19
DILUPBiCG: Solving for Ux, Initial residual = 0.000620827, Final residual = 1.78317e-09, No Iterations 11
DILUPBiCG: Solving for Uy, Initial residual = 0.000546878, Final residual = 1.73103e-09, No Iterations 10
DILUPBiCG: Solving for Uz, Initial residual = 0.000262496, Final residual = 6.11569e-10, No Iterations 11
GAMG: Solving for p, Initial residual = 0.158795, Final residual = 0.0197291, No Iterations 20
time step continuity errors : sum local = 3.89702e-08, global = -5.62443e-13, cumulative = -1.58585e-12
GAMG: Solving for p, Initial residual = 0.0769585, Final residual = 0.00629057, No Iterations 20
time step continuity errors : sum local = 4.59415e-08, global = -5.9324e-13, cumulative = -2.17909e-12
PIMPLE: iteration 20
DILUPBiCG: Solving for Ux, Initial residual = 0.0033839, Final residual = 2.3056e-08, No Iterations 65
DILUPBiCG: Solving for Uy, Initial residual = 0.00399156, Final residual = 2.71118e-08, No Iterations 56
DILUPBiCG: Solving for Uz, Initial residual = 0.00267406, Final residual = 1.59844e-08, No Iterations 65
GAMG: Solving for p, Initial residual = 0.191895, Final residual = 0.00318539, No Iterations 20
time step continuity errors : sum local = 3.28746e-08, global = -6.25726e-12, cumulative = -8.43635e-12
GAMG: Solving for p, Initial residual = 0.0450605, Final residual = 0.00172445, No Iterations 20
time step continuity errors : sum local = 3.66081e-07, global = -5.65159e-11, cumulative = -6.49523e-11
DILUPBiCG: Solving for omega, Initial residual = 0.0624003, Final residual = 4.45631e-07, No Iterations 49
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 8.95703e-06, No Iterations 34
PIMPLE: not converged within 20 iterations
ExecutionTime = 374.04 s ClockTime = 374 s
forces forces output:
sum of forces:
pressure : (-2.28408e+17 5.32264e+17 6.93677e+17)
viscous : (-2.95957e+09 9.28817e+07 -1.33636e+10)
porous : (0 0 0)
sum of moments:
pressure : (4.85726e+20 2.13491e+20 1.42405e+19)
viscous : (2.56719e+11 2.87508e+11 1.30512e+11)
porous : (0 0 0)
Courant Number mean: 5.74294e-06 max: 1.27126e+06
Time = 2.18166e-05
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0575416, Final residual = 2.68272e-07, No Iterations 10
DILUPBiCG: Solving for Uy, Initial residual = 0.0395307, Final residual = 2.18704e-07, No Iterations 10
DILUPBiCG: Solving for Uz, Initial residual = 0.0887286, Final residual = 6.59463e-07, No Iterations 11
GAMG: Solving for p, Initial residual = 0.0210242, Final residual = 0.000440192, No Iterations 20
time step continuity errors : sum local = 6.20467e-07, global = -1.58818e-09, cumulative = -1.65313e-09
GAMG: Solving for p, Initial residual = 0.00116681, Final residual = 1.03143e-05, No Iterations 20
time step continuity errors : sum local = 4.39509e-07, global = -2.29711e-10, cumulative = -1.88284e-09
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.104692, Final residual = 7.13485e-07, No Iterations 10
DILUPBiCG: Solving for Uy, Initial residual = 0.293159, Final residual = 1.94508e-06, No Iterations 9
DILUPBiCG: Solving for Uz, Initial residual = 0.0575179, Final residual = 1.97897e-07, No Iterations 10
GAMG: Solving for p, Initial residual = 0.000959239, Final residual = 0.000823077, No Iterations 20
time step continuity errors : sum local = 3.69348e-05, global = -2.23494e-08, cumulative = -2.42322e-08
|
|
|
|
|
|
|
October 15, 2015, 10:43
|
|
#20 |
|
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
Firstly, this is not the reason for your problems but: In my view the relTol of 1e-5 is insane. You see it in the later iterations that the linear solvers take 10 or even 65 iterations to converge. Also pressure takes all 20 iterations. In some ordinary simulations all solvers should take something like less than 5 iterations, often just 1 or two. And the outer algorithm still remains stable. There is no increase in overall accuracy by setting relTol to such a low value.
Now, to the actual problem. PIMPLE doesn't converge. You can see the residuals increasing after some initial decrease. You can try to use more stable numerics in the fvSchemes file. Use (1st order) upwind also for the velocity and uncorrected scheme for laplacian instead of corrected. Did you try to initialize the flow field with potentialFoam?
__________________
The skeleton ran out of shampoo in the shower. |
|
|
|
|
|
 |
| Tags |
| debug information, floating point exception |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| error compiling modified applications | yvyan | OpenFOAM Programming & Development | 21 | March 1, 2016 05:53 |
| [OpenFOAM] Native ParaView Reader Bugs | tj22 | ParaView | 270 | January 4, 2016 12:39 |
| [General] 2 datas on one plot | Akuji | ParaView | 46 | December 1, 2013 15:06 |
| Floating point error and divergence detected | aannjj | FLUENT | 0 | July 2, 2013 04:44 |
| [Gmsh] Gmsh and samplesurface | touf | OpenFOAM Meshing & Mesh Conversion | 2 | December 10, 2007 03:27 |
2Likes
Linear Mode
