CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Floating point exception with multiphaseEulerFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 14, 2023, 09:06
Default Floating point exception with multiphaseEulerFoam
  #1
New Member
 
Björn
Join Date: Jun 2023
Posts: 2
Rep Power: 0
Björnpork is on a distinguished road
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:rintStack(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:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
#6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<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:ow<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:haseSyste m>::momentumTransfer() at ??:?
#15 Foam::PhaseTransferPhaseSystem<Foam::OneResistance HeatTransferPhaseSystem<Foam::MomentumTransferPhas eSystem<Foam:haseSystem> > >::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;
}
}
Björnpork is offline   Reply With Quote

Reply

Tags
multiphase modelling, openfoam, sigfpe

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 11:43
A floating point exception has occurred: floating point exception [Overflow]. starlight STAR-CCM+ 4 May 4, 2016 09:08
A floating point exception - SEM Model yansheng STAR-CCM+ 1 April 4, 2016 04:57
Floating point exception from twoPhaseEulerFoam openfoammaofnepo OpenFOAM Running, Solving & CFD 1 March 19, 2016 13:56
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 19:28


All times are GMT -4. The time now is 01:52.