simpomann |
December 12, 2012 07:39 |
interFoam/kOmegaSST tank filling with printStackError/Mules
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).
Timestep decreases and this happens. I checked velocities and pressure at the last working timestep, nothing out of hand here.
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);
}
}
// ************************************************************************* //
p_rgh
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;
}
}
// ************************************************************************* //
k
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;
}
}
// ************************************************************************* //
Omega
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;
}
}
// ************************************************************************* //
Nut
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;
}
}
I guess inlet is self explaining. The only outlet is "atmosphaere". All the other stuff are pipes etc, so basically walls.
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
|