|
[Sponsors] |
October 15, 2015, 04: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, 04: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, 05: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, 05: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, 05: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, 05: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, 06: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, 06: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, 06: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, 08: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, 08: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, 08: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, 09: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, 09: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, 09: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, 09: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, 09: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, 09: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, 09: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, 09: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 |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
error compiling modified applications | yvyan | OpenFOAM Programming & Development | 21 | March 1, 2016 04:53 |
[OpenFOAM] Native ParaView Reader Bugs | tj22 | ParaView | 270 | January 4, 2016 11:39 |
[General] 2 datas on one plot | Akuji | ParaView | 46 | December 1, 2013 14:06 |
Floating point error and divergence detected | aannjj | FLUENT | 0 | July 2, 2013 03:44 |
[Gmsh] Gmsh and samplesurface | touf | OpenFOAM Meshing & Mesh Conversion | 2 | December 10, 2007 02:27 |