|
[Sponsors] |
Floating point exception with multiphaseEulerFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 14, 2023, 09:06 |
Floating point exception with multiphaseEulerFoam
|
#1 |
New Member
Björn
Join Date: Jun 2023
Posts: 2
Rep Power: 0 |
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; } } |
|
Tags |
multiphase modelling, openfoam, sigfpe |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |