CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

error in pimpleDyMFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 8, 2017, 04:07
Default error in pimpleDyMFoam
  #1
Member
 
bye bye my blue
Join Date: Sep 2016
Posts: 37
Rep Power: 2
bye bye my blue is an unknown quantity at this point
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:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:ICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam:ICPreconditioner:ICPreconditioner(Foam::l duMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix:reconditioner::addsymMatrixConst ructorToTable<Foam:ICPreconditioner>::New(Foam:: lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix:reconditioner::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;
}

















}

// ************************************************** *********************** //
bye bye my blue is offline   Reply With Quote

Old   February 2, 2017, 01:49
Default
  #2
New Member
 
Andre C
Join Date: Nov 2013
Posts: 3
Rep Power: 5
spartan42 is on a distinguished road
Have you checked your mesh?

I had a similar error. Turns out after performing some transformation the mesh failed 4 check test!
spartan42 is offline   Reply With Quote

Old   February 2, 2017, 11:34
Default
  #3
Senior Member
 
khedar
Join Date: Oct 2016
Posts: 110
Rep Power: 2
khedar is on a distinguished road
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.
khedar is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Access pointDisplacement field from within pimpleDyMFoam paul b OpenFOAM Programming & Development 0 December 2, 2016 04:16
pimpleDyMFoam error message laurentb OpenFOAM Running, Solving & CFD 7 May 13, 2015 05:48
Savonius pimpleDyMFoam vitokad OpenFOAM Pre-Processing 10 September 16, 2014 09:30
Problem with transient simulation of pump inner flow by pimpleDyMFoam on 1.6-ext renyun0511 OpenFOAM Running, Solving & CFD 0 April 10, 2013 03:26
pimpleDyMFoam samiam1000 OpenFOAM 2 September 19, 2012 10:11


All times are GMT -4. The time now is 12:25.