Floating point exception with multiphaseEulerFoam
 

Hello everyone :)
First, I'm pretty new to OpenFoam.



I try to simulate particle flow (air as fluid, solids as particles) in a cyclone seperator. Not with the model from the OpenFoam tutorials but one that I downloaded as a step file and created the mesh in Salome (using ideasUnvToFoam). I'm using OpenFoam10.



When I try to run the case with multiphaseEulerFoam I get the following error message:


#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 powf64 in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
#6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
#7 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, double const&) at ??:?
#8 Foam::dragModels::WenYu::CdRe() const at ??:?
#9 Foam::dragModels::GidaspowErgunWenYu::CdRe() const at ??:?
#10 Foam::dragModels::dispersedDragModel::Ki() const at ??:?
#11 Foam::dragModels::dispersedDragModel::K() const at ??:?
#12 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModel<Foam::dragModel>::ev aluate<double, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > (Foam::dragModel::*)() const, Foam::word const&, Foam::dimensionSet const&, bool) const at ??:?
#13 Foam::blendedDragModel::K() const at ??:?
#14 Foam::MomentumTransferPhaseSystem<Foam::phaseSyste m>::momentumTransfer() at ??:?
#15 Foam::PhaseTransferPhaseSystem<Foam::OneResistance HeatTransferPhaseSystem<Foam::MomentumTransferPhas eSystem<Foam::phaseSystem> > >::momentumTransfer() at ??:?
#16 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#17 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#18 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#19 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
Gleitkomma-Ausnahme (Speicherabzug geschrieben)




I tried a lot but what helped me bit was to create a really coarse mesh and set down the timestep. But event then I get the error message after a few time steps. At least I figured out that the temperature of the gas is negative in the last timestep (which is not really physical). The infomation of the last time step is as follows:


Courant Number mean: 0.00119300503 max: 11.6762713
deltaT = 3.74431908e-08
Time = 0.028879063s

PIMPLE: Iteration 1
MULES: Solving for alpha.solids
MULES: Solving for alpha.gas
smoothSolver: Solving for alpha.solids, Initial residual = 8.69713291e-11, Final residual = 2.39318033e-17, No Iterations 1
smoothSolver: Solving for alpha.gas, Initial residual = 8.69715171e-11, Final residual = 2.55978747e-16, No Iterations 1
solids fraction, min, max = 0.00177253757 -1.27363306e-15 0.311543533
gas fraction, min, max = 0.998227462 0.688456467 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.solids
MULES: Solving for alpha.gas
smoothSolver: Solving for alpha.solids, Initial residual = 8.69713036e-11, Final residual = 2.70647871e-17, No Iterations 1
smoothSolver: Solving for alpha.gas, Initial residual = 8.69715194e-11, Final residual = 2.63513492e-16, No Iterations 1
solids fraction, min, max = 0.00177253757 -1.27363272e-15 0.311543541
gas fraction, min, max = 0.998227462 0.688456459 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing face momentum equations
smoothSolver: Solving for h.solids, Initial residual = 4.61536121e-06, Final residual = 1.74615236e-14, No Iterations 1
smoothSolver: Solving for h.gas, Initial residual = 0.990169816, Final residual = 0.00181652713, No Iterations 10
solids min/max T 113.672551 - 764.376261
gas min/max T -340.372963 - 756.472606
GAMG: Solving for p_rgh, Initial residual = 0.0158926039, Final residual = 1.0533868e-14, No Iterations 1
limitPressure: min -14913.6638
PIMPLE: Iteration 2
MULES: Solving for alpha.solids
MULES: Solving for alpha.gas
smoothSolver: Solving for alpha.solids, Initial residual = 3.73381758e-12, Final residual = 2.34474767e-17, No Iterations 1
smoothSolver: Solving for alpha.gas, Initial residual = 3.73408379e-12, Final residual = 2.29119313e-16, No Iterations 1
solids fraction, min, max = 0.00177253757 -1.17839012e-15 0.311543541
gas fraction, min, max = 0.998227462 0.688456459 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.solids
MULES: Solving for alpha.gas
smoothSolver: Solving for alpha.solids, Initial residual = 3.73382093e-12, Final residual = 2.39683063e-17, No Iterations 1
smoothSolver: Solving for alpha.gas, Initial residual = 3.73409853e-12, Final residual = 2.19498013e-16, No Iterations 1
solids fraction, min, max = 0.00177253757 -1.18956538e-15 0.311543549
gas fraction, min, max = 0.998227462 0.688456451 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing face momentum equations



So the error occurs when it ries to construct the face momentum equations I assume? By what is this caused?
Help will be much appreciated!



The following stuff are the files from the constant folder:
momentumTransport.gas:
simulationType RAS;

RAS
{
model kEpsilon;

turbulence on;
printCoeffs on;
}










momentumTransport.solids:
simulationType RAS;

RAS
{
model kineticTheory;

turbulence on;
printCoeffs on;

kineticTheoryCoeffs
{
equilibrium off;

e 0.8;
alphaMinFriction 0.5;
residualAlpha 1e-5;

granularViscosityModel Syamlal;
granularConductivityModel Syamlal;
granularPressureModel SyamlalRogersOBrien;
frictionalStressModel Schaeffer;
radialModel CarnahanStarling;

SchaefferCoeffs
{
phi 36;
}
}
}








phaseProperties:
type basicMultiphaseSystem;

phases (solids gas);

solids
{
type purePhaseModel;

diameterModel constant;

constantCoeffs
{
d 462e-6;
}

alphaMax 0.65;
residualAlpha 1e-6;
}

gas
{
type purePhaseModel;

diameterModel none;

residualAlpha 1e-5;
}

blending
{
default
{
type continuous;
phase gas;
}
}

surfaceTension
{
gas_solids
{
type constant;
sigma 0;
}
}

drag
{
solids_dispersedIn_gas
{
type GidaspowErgunWenYu;
residualAlpha 1e-5;
residualRe 1e-5;
}
}

virtualMass
{
solids_dispersedIn_gas
{
type constantCoefficient;
Cvm 0;
}
}

heatTransfer
{
solids_dispersedIn_gas
{
type RanzMarshall;
residualAlpha 1e-4;
}
}

phaseTransfer
{}

lift
{}

wallLubrication
{}

turbulentDispersion
{}

interfaceCompression
{}








physicalProperties.gas:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

mixture
{
specie
{
molWeight 28.9;
}
thermodynamics
{
Cp 1007;
Hf 0;
}
equationOfState
{
rho 1.2;
}
transport
{
mu 1.84e-05;
Pr 0.7;
}
}







physicalProperties.solids:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

mixture
{
specie
{
molWeight 100;
}
equationOfState
{
rho 2480;
}
thermodynamics
{
Cp 6000;
Hf 0;
}
transport
{
mu 0;
Pr 24.47;
}
}