CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   error in pimpleDyMFoam (https://www.cfd-online.com/Forums/openfoam/182359-error-pimpledymfoam.html)

bye bye my blue January 8, 2017 03:07

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;
}

















}

// ************************************************** *********************** //

spartan42 February 2, 2017 00:49

Have you checked your mesh?

I had a similar error. Turns out after performing some transformation the mesh failed 4 check test!

khedar February 2, 2017 10:34

1. May be you should try more nCorrectors and lesses nOuterCorrectors

Code:

PIMPLE
{
    correctPhi          yes;
    nOuterCorrectors    1;
    nCorrectors        10;  // or 20 or 50
    nNonOrthogonalCorrectors 1;}

2. Try using PCG solver for p equation at the start until it runs fine and then switch to GAMG.


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