 |
 |
 |
|
|
|
[Sponsors] | ||||
interFoam/kOmegaSST tank filling with printStackError/Mules |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 12, 2012, 08:39
|
interFoam/kOmegaSST tank filling with printStackError/Mules
|
#1 |
|
Member
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 14  |
Hey there,
I am simulating the filling of a tank with interfoam and try to use k-Omega-SST model for turbulence. The mesh has around 400.000 cells, checkMesh says okay. Unluckily I encounter these errors a lot (sometimes after 0.25s, with different transport properties after 0.4s): Code:
MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.524686, Final residual = 0.000681733, No Iterations 2 time step continuity errors : sum local = 1.32975e-10, global = 2.51813e-15, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0920533, Final residual = 7.70789e-08, No Iterations 5 time step continuity errors : sum local = 2.09622e-14, global = -1.45744e-18, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 0.000158889, Final residual = 7.80272e-08, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.00244614, Final residual = 5.717e-10, No Iterations 2 bounding k, min: -0.000200721 max: 3.00901e+20 average: 1.31669e+15 ExecutionTime = 24866.4 s ClockTime = 24916 s Courant Number mean: 2.01358e-08 max: 0.712635 Interface Courant Number mean: 8.29938e-09 max: 0.282305 deltaT = 2.2569e-14 Time = 0.2550889169617 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.408204, Final residual = 0.000561119, No Iterations 2 time step continuity errors : sum local = 2.68439e-10, global = -5.36369e-15, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0698315, Final residual = 5.11837e-08, No Iterations 5 time step continuity errors : sum local = 3.27413e-14, global = -7.29092e-18, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 0.00329008, Final residual = 2.10387e-10, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 0.0124975, Final residual = 7.00019e-09, No Iterations 2 ExecutionTime = 24868.6 s ClockTime = 24919 s Courant Number mean: 3.59009e-08 max: 3.49031 Interface Courant Number mean: 1.43023e-08 max: 0.550178 deltaT = 5.17294e-15 Time = 0.255088916961705 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28037e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.587722, Final residual = 0.000830391, No Iterations 2 time step continuity errors : sum local = 7.64239e-11, global = 4.21113e-16, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0960292, Final residual = 3.83386e-08, No Iterations 6 time step continuity errors : sum local = 5.13861e-15, global = -9.16091e-19, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 0.000484156, Final residual = 1.42813e-08, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.00071846, Final residual = 1.26128e-09, No Iterations 2 ExecutionTime = 24870.9 s ClockTime = 24921 s Courant Number mean: 1.32675e-08 max: 0.702258 Interface Courant Number mean: 5.5049e-09 max: 0.202586 deltaT = 5.76223e-15 Time = 0.25508891696171 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.143881, Final residual = 0.000365565, No Iterations 2 time step continuity errors : sum local = 4.46931e-11, global = -4.75218e-15, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0402696, Final residual = 3.67927e-08, No Iterations 6 time step continuity errors : sum local = 4.85566e-15, global = -1.63317e-18, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 8.34399e-05, Final residual = 9.95434e-12, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 0.00146217, Final residual = 2.85304e-10, No Iterations 2 ExecutionTime = 24873.2 s ClockTime = 24923 s Courant Number mean: 1.83689e-08 max: 0.546694 Interface Courant Number mean: 7.7516e-09 max: 0.262968 deltaT = 6.60544e-15 Time = 0.255088916961717 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.0613232, Final residual = 0.000199917, No Iterations 2 time step continuity errors : sum local = 3.64953e-11, global = 1.40235e-15, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.00862543, Final residual = 7.72936e-08, No Iterations 5 time step continuity errors : sum local = 1.45349e-14, global = 4.97913e-17, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 0.000305148, Final residual = 9.20208e-08, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.00171966, Final residual = 6.46387e-10, No Iterations 2 bounding k, min: -2.15624e+09 max: 2.68344e+21 average: 1.08957e+16 ExecutionTime = 24875.4 s ClockTime = 24925 s Courant Number mean: 2.361e-08 max: 0.507769 Interface Courant Number mean: 1.01369e-08 max: 0.318716 deltaT = 7.64614e-15 Time = 0.255088916961725 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00200676 Min(alpha1) = -6.28036e-20 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.24767, Final residual = 0.000361503, No Iterations 2 time step continuity errors : sum local = 1.34031e-10, global = 3.46138e-15, cumulative = -1.01128e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0585152, Final residual = 3.09369e-08, No Iterations 5 time step continuity errors : sum local = 1.38894e-14, global = 8.65978e-17, cumulative = -1.01128e-05 DILUPBiCG: Solving for omega, Initial residual = 0.00217707, Final residual = 1.56868e-10, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 0.00609363, Final residual = 2.66106e-09, No Iterations 2 ExecutionTime = 24877.6 s ClockTime = 24928 s Courant Number mean: 3.25704e-08 max: 1.92022 Interface Courant Number mean: 1.36638e-08 max: 0.455131 deltaT = 3.18553e-15 Time = 0.255088916961728 [3] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #2 in "/lib64/libc.so.6" [3] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [3] #4 at MULES.C:0 [3] #5 Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [3] #6 [3] in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam" [3] #7 __libc_start_main in "/lib64/libc.so.6" [3] #8 [3] at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116 [c103-06:14726] *** Process received signal *** [c103-06:14726] Signal: Floating point exception (8) [c103-06:14726] Signal code: (-6) [c103-06:14726] Failing at address: 0x35c600003986 [c103-06:14726] [ 0] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10] [c103-06:14726] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7fd91ae43d95] [c103-06:14726] [ 2] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10] [c103-06:14726] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7fd91ca71fa8] [c103-06:14726] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7fd91ca6507e] [c103-06:14726] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7fd91ca65310] [c103-06:14726] [ 6] interFoam() [0x42d35e] [c103-06:14726] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7fd91ae3023d] [c103-06:14726] [ 8] interFoam() [0x425119] [c103-06:14726] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 3 with PID 14726 on node c103-06 exited on signal 8 (Floating point exception). I guess I have to supply my fvSchemes, fvSolutions and my initial conditions to trace this. My guess is a connection to the k-Omega-SST model. Probably my initial settings are bad, but I calculated them several times following the wiki here and can't manage to get it work. Case is currently starting laminar, if this works I will update you asap. U Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type fixedValue;
value uniform (0 0 0);
}
//Durchströmt
inlet
{
type flowRateInletVelocity;
flowRate constant 0.000666667;//incompressible solver: [m³/s]
value uniform (0 0 0); //placeholder
//type fixedValue; //Alternative, Sascha
//value (1.36 -2.56 0);
}
atmosphaere_CATIASTL
{
type zeroGradient;
}
//Alles was umströmt wird
tank_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
efr_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
entl_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
pistole_vorn_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
pistole_hinten_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
innenteile_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
flv_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
einsatz_efr_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
verschlusskopf_CATIASTL
{
type fixedValue;
value uniform (0 0 0);
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type buoyantPressure;
value uniform 0;
}
//Durchströmt
inlet
{
type buoyantPressure;
value uniform 0;
}
atmosphaere_CATIASTL
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;
}
//Alles was umströmt wird
tank_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
efr_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
entl_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
pistole_vorn_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
// sd_CATIASTL
//{
// type buoyantPressure;
// value uniform 0;
//}
pistole_hinten_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
innenteile_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
flv_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
einsatz_efr_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
verschlusskopf_CATIASTL
{
type buoyantPressure;
value uniform 0;
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0.0286;
boundaryField
{
walls
{
type kqRWallFunction;
value $internalField;
}
//Durchströmt
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.041; // 5% turbulent intensity
value $internalField;
}
atmosphaere_CATIASTL
{
type inletOutlet;
inletValue $internalField;
}
//Alles was umströmt wird
tank_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
efr_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
entl_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
pistole_vorn_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
sd_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
pistole_hinten_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
innenteile_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
flv_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
einsatz_efr_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
verschlusskopf_CATIASTL
{
type kqRWallFunction;
value $internalField;
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 1643.71;
boundaryField
{
walls
{
type omegaWallFunction;
value $internalField;
}
//Durchströmt
inlet
{
type fixedValue;
value uniform 280.78;
}
atmosphaere_CATIASTL
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}
//Alles was umströmt wird
tank_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
efr_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
entl_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
pistole_vorn_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
sd_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
pistole_hinten_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
innenteile_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
flv_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
einsatz_efr_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
verschlusskopf_CATIASTL
{
type omegaWallFunction;
value $internalField;
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -1 0 0 0 0 ];
internalField uniform 0.00010198;
boundaryField
{
walls
{
type nutkWallFunction;
value $internalField;
}
//Durchströmt
inlet
{
type calculated;
value uniform 0.00010198;
}
atmosphaere_CATIASTL
{
type calculated;
value $internalField;
}
//Alles was umströmt wird
tank_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
efr_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
entl_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
pistole_vorn_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
sd_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
pistole_hinten_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
innenteile_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
flv_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
einsatz_efr_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
verschlusskopf_CATIASTL
{
type nutkWallFunction;
value $internalField;
}
}
If somebody spots a mistake in one of the files please tell me! I am super lost with this. All help is highly appreciated! Big thanks and best regards! I hope somebody with more experience can help me to fix this, unluckily it's a bit urgent for my thesis. Greetings! Simon |
|
|
|
|
|
December 12, 2012, 08:40
|
|
#2 |
|
Member
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 14 |
fvSchemes
Code:
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear; //Alternative:default cellMDLimited Gauss linear 1
}
divSchemes
{
div(rho*phi,U) Gauss limitedLinearV 1;
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression;
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;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
pcorr;
alpha1;
p_rgh;
}
Code:
solvers
{
p_rgh //Lösung der Druckgleichung
{
solver GAMG;
tolerance 1e-07;
relTol 0.01;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
p_rghFinal //finaler Loop
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-07;
relTol 0;
nVcycles 2;
smoother DICGaussSeidel;
nPreSweeps 2;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
tolerance 1e-07;
relTol 0;
maxIter 10;
}
pcorr //vor dem Loop gelöste Druckvergleichung für Startwerte des Feldes innerhalb des Iterationsschritts
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-05;
relTol 0;
smoother DICGaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration false;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
tolerance 1e-05;
relTol 0;
maxIter 10;
}
"(U|k|omega)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-06;
relTol 0.1;
nSweeps 1;
}
"(U|k|omega)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
}
}
PIMPLE
{
momentumPredictor no;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
nAlphaCorr 1;
nAlphaSubCycles 2;
cAlpha 1; // Alternative 0
}
And alpha (dont expect problems here) Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type zeroGradient;
}
//Durchströmt
inlet
{
type fixedValue;
value uniform 1;
}
atmosphaere_CATIASTL
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
//Alles was umströmt wird
tank_CATIASTL
{
type zeroGradient;
}
efr_CATIASTL
{
type zeroGradient;
}
entl_CATIASTL
{
type zeroGradient;
}
pistole_vorn_CATIASTL
{
type zeroGradient;
}
sd_CATIASTL
{
type zeroGradient;
}
pistole_hinten_CATIASTL
{
type zeroGradient;
}
innenteile_CATIASTL
{
type zeroGradient;
}
flv_CATIASTL
{
type zeroGradient;
}
einsatz_efr_CATIASTL
{
type zeroGradient;
}
verschlusskopf_CATIASTL
{
type zeroGradient;
}
}
// ************************************************************************* //
|
|
|
|
|
|
|
December 12, 2012, 13:11
|
|
#3 |
|
Member
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 14 |
Hey,
I checked a lamniar case. Same error but at a different timestep Code:
MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62932e-21 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62932e-21 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.065834, Final residual = 0.000307154, No Iterations 2 time step continuity errors : sum local = 2.85125e-12, global = -6.3837e-15, cumulative = -3.90402e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.02062, Final residual = 3.97336e-08, No Iterations 7 time step continuity errors : sum local = 3.9867e-16, global = -5.09275e-19, cumulative = -3.90402e-05 ExecutionTime = 15548.4 s ClockTime = 15612 s Courant Number mean: 2.84517e-08 max: 0.867998 Interface Courant Number mean: 5.5157e-09 max: 0.0944208 deltaT = 4.14157e-15 Time = 0.202400597327766 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.0662251, Final residual = 0.000298758, No Iterations 2 time step continuity errors : sum local = 2.77135e-12, global = -4.86472e-15, cumulative = -3.90402e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0204396, Final residual = 3.2508e-08, No Iterations 7 time step continuity errors : sum local = 3.2613e-16, global = -4.25822e-19, cumulative = -3.90402e-05 ExecutionTime = 15550.6 s ClockTime = 15614 s Courant Number mean: 2.84501e-08 max: 0.867991 Interface Courant Number mean: 5.51484e-09 max: 0.0943805 deltaT = 3.81716e-15 Time = 0.202400597327769 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.0657749, Final residual = 0.00031065, No Iterations 2 time step continuity errors : sum local = 2.88026e-12, global = -6.24636e-15, cumulative = -3.90402e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0206025, Final residual = 3.60495e-08, No Iterations 7 time step continuity errors : sum local = 3.61229e-16, global = -3.50359e-19, cumulative = -3.90402e-05 ExecutionTime = 15552.7 s ClockTime = 15616 s Courant Number mean: 2.84444e-08 max: 0.867923 Interface Courant Number mean: 5.6614e-09 max: 0.094348 deltaT = 3.51843e-15 Time = 0.202400597327773 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.0661369, Final residual = 0.000300055, No Iterations 2 time step continuity errors : sum local = 2.77963e-12, global = -4.82001e-15, cumulative = -3.90402e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0204266, Final residual = 3.17341e-08, No Iterations 7 time step continuity errors : sum local = 3.17909e-16, global = -3.68428e-19, cumulative = -3.90402e-05 ExecutionTime = 15554.9 s ClockTime = 15619 s Courant Number mean: 2.8439e-08 max: 0.867925 Interface Courant Number mean: 5.317e-09 max: 0.094304 deltaT = 3.24307e-15 Time = 0.202400597327776 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 MULES: Solving for alpha1 Liquid phase volume fraction = 0.00157538 Min(alpha1) = -1.62931e-21 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.0657622, Final residual = 0.00031015, No Iterations 2 time step continuity errors : sum local = 2.87193e-12, global = -6.3359e-15, cumulative = -3.90402e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.0205844, Final residual = 3.79572e-08, No Iterations 7 time step continuity errors : sum local = 3.79847e-16, global = -3.64393e-19, cumulative = -3.90402e-05 ExecutionTime = 15557.1 s ClockTime = 15621 s Courant Number mean: 2.84311e-08 max: 0.867857 Interface Courant Number mean: 5.31414e-09 max: 0.0942674 deltaT = 2.9895e-15 Time = 0.202400597327779 [3] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #2 in "/lib64/libc.so.6" [3] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [3] #4 at MULES.C:0 [3] #5 Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [3] #6 [3] in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam" [3] #7 __libc_start_main in "/lib64/libc.so.6" [3] #8 [3] at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116 [c103-06:28794] *** Process received signal *** [c103-06:28794] Signal: Floating point exception (8) [c103-06:28794] Signal code: (-6) [c103-06:28794] Failing at address: 0x35c60000707a [c103-06:28794] [ 0] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10] [c103-06:28794] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7f3d0d63ad95] [c103-06:28794] [ 2] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10] [c103-06:28794] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7f3d0f268fa8] [c103-06:28794] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7f3d0f25c07e] [c103-06:28794] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7f3d0f25c310] [c103-06:28794] [ 6] interFoam() [0x42d35e] [c103-06:28794] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7f3d0d62723d] [c103-06:28794] [ 8] interFoam() [0x425119] [c103-06:28794] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 3 with PID 28794 on node c103-06 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Could it be this part from my fvSolution? PIMPLE { momentumPredictor no; nCorrectors 2; nNonOrthogonalCorrectors 0; nAlphaCorr 1; nAlphaSubCycles 2; cAlpha 1; // Alternative 0 } Last edited by simpomann; December 12, 2012 at 13:43. |
|
|
|
|
|
|
February 17, 2014, 18:06
|
|
#4 |
|
New Member
Matej Muller
Join Date: Oct 2011
Location: Slovenia
Posts: 25
Rep Power: 15 |
Hey Simon!
Have you perhaps solved this issue? I'm also trying to run interFoam with the kOmega model and also don't know which solvers and schemes to use, since the tutorial of kOmegaSST exists only for the LTSInterFoam. Have you figured out the right solvers and schemes? thanks for the help, Matej |
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Filling Tank with Water | leff | CFX | 7 | August 21, 2017 08:47 |
| BCs to the outlet patch (filling a simple tank) | Ternox | OpenFOAM Pre-Processing | 4 | August 21, 2017 08:19 |
| tank filling (water/air), I think I messed up the BCs | simpomann | OpenFOAM Running, Solving & CFD | 2 | August 21, 2017 08:16 |
| free-surface liquid filling a fuel tank | lee | CFX | 2 | June 4, 2008 22:01 |
| how to set boundary condition for a tank filling w | namrata | FLUENT | 6 | May 2, 2007 11:23 |
Linear Mode
