[Raji]

Hi to everyone,

I am trying to simulate a spouted bed (CFD) using pisoFoam. I am using kepslion turbulence model. When i run the case it is always giving following floating point error:

Time = 0.0773

Courant Number mean: 0.056338759 max: 1.4369808
DILUPBiCG: Solving for Ux, Initial residual = 0.0045467778, Final residual = 4.1594073e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0011529659, Final residual = 1.1293885e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0058914324, Final residual = 3.5906389e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.016142261, Final residual = 0.001555332, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.0015991186, Final residual = 0.00015579393, No Iterations 25
time step continuity errors : sum local = 2.4741439e-07, global = -8.0117863e-09, cumulative = 3.0139362e-06
DICPCG: Solving for p, Initial residual = 0.0026050586, Final residual = 0.00021336197, No Iterations 4
DICPCG: Solving for p, Initial residual = 0.00026164929, Final residual = 2.3073015e-05, No Iterations 25
time step continuity errors : sum local = 3.661006e-08, global = -1.1797049e-08, cumulative = 3.0021391e-06
DICPCG: Solving for p, Initial residual = 0.00074208285, Final residual = 5.9053786e-05, No Iterations 2
DICPCG: Solving for p, Initial residual = 6.4034064e-05, Final residual = 6.1044294e-06, No Iterations 47
time step continuity errors : sum local = 9.6839112e-09, global = 1.1177187e-09, cumulative = 3.0032569e-06
DICPCG: Solving for p, Initial residual = 0.00027902805, Final residual = 1.5131582e-05, No Iterations 2
DICPCG: Solving for p, Initial residual = 1.8245823e-05, Final residual = 9.8413354e-07, No Iterations 45
time step continuity errors : sum local = 1.5590186e-09, global = 2.1035233e-10, cumulative = 3.0034672e-06
DILUPBiCG: Solving for epsilon, Initial residual = 0.00010341292, Final residual = 7.0807648e-07, No Iterations 1
bounding epsilon, min: -55.138134 max: 10498.612 average: 75.832938
DILUPBiCG: Solving for k, Initial residual = 2.3604908e-06, Final residual = 2.3604908e-06, No Iterations 0
bounding k, min: 0 max: 10.185168 average: 0.31627826
[1] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[1] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[1] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[1] #6 Foam::incompressible::RASModels::realizableKE::rCm u(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[1] #7 Foam::incompressible::RASModels::realizableKE::cor rect() at ??:?
[1] #8
[1] at ??:?
[1] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #10
[1] at ??:?
[pk:13432] *** Process received signal ***
[pk:13432] Signal: Floating point exception (8)
[pk:13432] Signal code: (-6)
[pk:13432] Failing at address: 0x3e900003478
[pk:13432] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x352f0) [0x7f0f993322f0]
[pk:13432] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f0f99332267]
[pk:13432] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x352f0) [0x7f0f993322f0]
[pk:13432] [ 3] /opt/openfoam231/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_ 5UListIdEES6_+0xd5) [0x7f0f9a605c05]
[pk:13432] [ 4] /opt/openfoam231/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam6divideINS_1 2fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldI dT_T0_EERKS6_S9_+0x7f) [0x7f0f9cb58b1f]
[pk:13432] [ 5] /opt/openfoam231/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4FoamdvINS_12fvPa tchFieldENS_7volMeshEEENS_3tmpINS_14GeometricField IdT_T0_EEEERKS8_RKS7_+0x14f) [0x7f0f9cb69daf]
[pk:13432] [ 6] /opt/openfoam231/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible9RASModels12realizableKE4rCmuERKNS_14Geometric FieldINS_6TensorIdEENS_12fvPatchFieldENS_7volMeshE EERKNS3_IdS6_S7_EESD_+0x4ac) [0x7f0f9cb78ccc]
[pk:13432] [ 7] /opt/openfoam231/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible9RASModels12realizableKE7correctEv+0xb97) [0x7f0f9cb7b7d7]
[pk:13432] [ 8] pisoFoam() [0x41b8f1]
[pk:13432] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0) [0x7f0f9931da40]
[pk:13432] [10] pisoFoam() [0x41c805]
[pk:13432] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 13432 on node pk exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


My initial condition for k and epsilon are as follows:

-> k:
dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.00325;

boundaryField
{
casing
{
type kqRWallFunction;
value $internalField;
}
steg
{
type kqRWallFunction;
value $internalField;
}
inlet_3
{
type fixedValue;
value $internalField;
}
inlet_1
{
type fixedValue;
value $internalField;
}
inlet_2
{
type fixedValue;
value $internalField;
}
inlet_4
{
type fixedValue;
value $internalField;
}
inlet_5
{
type fixedValue;
value $internalField;
}
inlet_6
{
type fixedValue;
value $internalField;
}
inlet_7
{
type fixedValue;
value $internalField;
}
inlet_8
{
type fixedValue;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue uniform 0;
value $internalField;
}
}

-> epsilon:

dimensions [0 2 -3 0 0 0 0];

internalField uniform 0.000765;

boundaryField
{
casing
{
type epsilonWallFunction;
value $internalField ;
}

steg
{
type epsilonWallFunction;
value $internalField ;
}
outlet
{
type inletOutlet;
inletValue uniform 0;
value $internalField;
}
inlet_1
{
type fixedValue;
value $internalField;
}
inlet_2
{
type fixedValue;
value $internalField;
}
inlet_3
{
type fixedValue;
value $internalField;
}
inlet_4
{
type fixedValue;
value $internalField;
}
inlet_5
{
type fixedValue;
value $internalField;
}
inlet_6
{
type fixedValue;
value $internalField;
}
inlet_7
{
type fixedValue;
value $internalField;
}
inlet_8
{
type fixedValue;
value $internalField;
}

}

Any comment and help will be really appreciated.

Thanks in advance.

[DancingButterfly]

Hey,

I have the same problem.. Did you figure out how to solve it?

Thank you for your help!