November 17, 2013, 01:15
|
LTSInterFoam Crashing
|
#1
|
Member
John Anonymous
Join Date: Jan 2011
Location: Melbourne Australia
Posts: 37
Rep Power: 16
|
Good day All
I am attempting to create a LTSInterFoam model as the interFoam models I currently have developed are taking an inordinate amount of time to finish (two months so far on 8/16 nodes with approximately 500,000 cells). Anyways the error output I am getting is:
Code:
Time = 1
Flow time scale min/max = 0.0223363, 1
Smoothed flow time scale min/max = 0.0223363, 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191047 Min(alpha1) = -93.8154 Max(alpha1) = 26.8419
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191234 Min(alpha1) = -2410 Max(alpha1) = 8705.97
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191417 Min(alpha1) = -808793 Max(alpha1) = 223717
[0] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4 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::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #8
[0] in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #10
[0] in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[Ubuntu:29885] *** Process received signal ***
[Ubuntu:29885] Signal: Floating point exception (8)
[Ubuntu:29885] Signal code: (-6)
[Ubuntu:29885] Failing at address: 0x3e8000074bd
[Ubuntu:29885] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fc8032b0425]
[Ubuntu:29885] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 3] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xd2) [0x7fc8043953c2]
[Ubuntu:29885] [ 4] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1_IS8_EESD_h+0xf1d) [0x7fc804397cfd]
[Ubuntu:29885] [ 5] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4d3) [0x7fc804398c43]
[Ubuntu:29885] [ 6] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7fc805e86e17]
[Ubuntu:29885] [ 7] LTSInterFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x152) [0x4562a2]
[Ubuntu:29885] [ 8] LTSInterFoam() [0x42c93a]
[Ubuntu:29885] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fc80329b76d]
[Ubuntu:29885] [10] LTSInterFoam() [0x43145d]
[Ubuntu:29885] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 29885 on node Ubuntu exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
I have been using the following fvSolutions file:
Code:
solvers
{
pcorr
{
solver PCG;
preconditioner
{
preconditioner GAMG;
smoother DICGaussSeidel;
agglomerator faceAreaPair;
mergeLevels 1;
nCellsInCoarsestLevel 10;
cacheAgglomeration true;
tolerance 1e-6;
relTol 0;
};
tolerance 1e-6;
relTol 0;
};
p_rgh
{
solver GAMG;
smoother GaussSeidel;
agglomerator faceAreaPair;
mergeLevels 1;
nCellsInCoarsestLevel 10;
cacheAgglomeration true;
tolerance 1e-7;
relTol 0.01;
};
p_rghFinal
{
$p_rgh;
tolerance 1e-7;
relTol 0;
}
U
{
solver smoothSolver; //PBiCG;
smoother GaussSeidel; //preconditioner DILU;
tolerance 1e-06;
relTol 0;
nSweeps 1;
}
UFinal
{
solver smoothSolver; //$U;
smoother GaussSeidel;
tolerance 1e-6;
relTol 0;
nSweeps 1;
}
/* "U.*"
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-8;
relTol 0;
};
*/
}
PIMPLE
{
momentumPredictor no; // yes;
nCorrectors 2; // 1;
nNonOrthogonalCorrectors 0;
nAlphaCorr 1;
nAlphaSubCycles 3; //1;
cAlpha 1;
correctPhi no;
maxCo 0.75;
maxAlphaCo 0.75;
nAlphaSweepIter 1;
rDeltaTSmoothingCoeff 0.1;
rDeltaTDampingCoeff 1;
maxDeltaT 1;
}
relaxationFactors
{
fields
{
}
equations
{
"U.*" 1;
}
}
// ************************************************************************* //
And fvSchemes file:
Code:
ddtSchemes
{
default localEuler rDeltaT;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
div(rho*phi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression;
// div(phi,k) Gauss upwind;
// div(phi,omega) Gauss upwind;
div((muEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p_rgh;
pcorr;
alpha1;
}
// ************************************************************************* //
I would appreciate if someone could have a look at these and let me know what I am missing. I suspect this is an issue with alpha?
Regards
JFM
|
|
|