error in pimpleDyMFoam
IMPLE: iteration 23
smoothSolver: Solving for Ux, Initial residual = 6.6535e-08, Final residual = 3.46032e-09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1.21205e-08, Final residual = 6.81186e-10, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 2.24525e-08, Final residual = 1.12119e-09, No Iterations 1 GAMG: Solving for p, Initial residual = 0.000130821, Final residual = 8.62738e-09, No Iterations 15 GAMG: Solving for p, Initial residual = 6.31422e-05, Final residual = 8.05468e-09, No Iterations 14 time step continuity errors : sum local = 2.16879e-13, global = -4.03075e-14, cumulative = 7.64631e-10 PIMPLE: converged in 23 iterations ExecutionTime = 130.99 s ClockTime = 209 s Courant Number mean: 0.00185858 max: 2.00168 deltaT = 4.79699e-05 Time = 0.00383814 AMI: Creating addressing and weights between 2118 source faces and 2118 target faces AMI: Patch source sum(weights) min/max/average = 0.999962, 1.00553, 1.0006 AMI: Patch target sum(weights) min/max/average = 0.999958, 1.0055, 1.0006 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 9.66287e-09, No Iterations 41 GAMG: Solving for pcorr, Initial residual = 0.34805, Final residual = 8.79421e-09, No Iterations 30 time step continuity errors : sum local = 1.05264e-14, global = -1.74044e-15, cumulative = 7.64629e-10 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.0154663, Final residual = 0.000213435, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00384019, Final residual = 7.15799e-05, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00752693, Final residual = 0.00013007, No Iterations 1 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::DICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) at ??:? #4 Foam::DICPreconditioner::DICPreconditioner(Foam::l duMatrix::solver const&, Foam::dictionary const&) at ??:? #5 Foam::lduMatrix::preconditioner::addsymMatrixConst ructorToTable<Foam::DICPreconditioner>::New(Foam:: lduMatrix::solver const&, Foam::dictionary const&) at ??:? #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:? #7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:? #9 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 ??:? #10 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #11 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? #12 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? #13 ? at ??:? #14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #15 ? at ??:? when I run pimpleDyMFoam, this error occurred.... I don't know why it happens, what made this problem... ---------fvschemes----------------------------- /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.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) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss linearUpwindV grad(U); div((nuEff*dev2(T(grad(U))))) Gauss linear; div(phi,k) Gauss upwind; div(phi,nuTilda) Gauss upwind; } laplacianSchemes { default Gauss linear limited corrected 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited corrected 0.33; } fluxRequired { default no; pcorr ; p ; } wallDist { method meshWave; } // ************************************************** *********************** // ----------- fvsolution /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.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 1e-08; relTol 0; } p { $pcorr; tolerance 1e-08; relTol 0; } pFinal { $p; tolerance 1e-08; relTol 0; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-08; relTol 0.1; } UFinal { $U; tolerance 1e-08; relTol 0; } cellMotionUx { solver PCG; preconditioner DIC; tolerance 1e-08; relTol 0; } "(U|nuTilda)" { solver smoothSolver; smoother GaussSeidel; tolerance 1e-05; relTol 0.1; } "(U|nuTilda)Final" { $U; tolerance 1e-06; relTol 0; } } PIMPLE { correctPhi yes; nOuterCorrectors 150; nCorrectors 1; nNonOrthogonalCorrectors 1; residualControl { U { tolerance 1e-04; relTol 1e-04; } p { tolerance 1e-04; relTol 1e-04; } nuTilda { tolerance 1e-04; relTol 1e-04; } nut { tolerance 1e-04; relTol 1e-04; } k { tolerance 1e-04; relTol 1e-04; } } pRefCell 0; pRefValue 0; } relaxationFactors { equations { "U.*" 0.3; "(U|nuTilda).*" 1; } fields { p 0.3; pFinal 1; } } // ************************************************** *********************** // |
Have you checked your mesh?
I had a similar error. Turns out after performing some transformation the mesh failed 4 check test! |
1. May be you should try more nCorrectors and lesses nOuterCorrectors
Code:
PIMPLE |
All times are GMT -4. The time now is 18:28. |