[Mahyar Javidi]

Hi everybody
I am running a new solver that I developed based on interFoam by adding a new field (volume fraction) to it and using a new viscosity model. It workes up to a certain time then I get this error: Floating point exception
I tried to change the solvers in fvSolution and schems in fvSchemes but no thing really changed.
I just copy/paste the last part that I get from solving a case by this solver, which may help to get a clue.



MULES: Solving for alpha1
Phase-1 volume fraction = 0.516683 Min(alpha1) = -6.43313e-18 Max(alpha1) = 1.1454
MULES: Solving for alpha1
Phase-1 volume fraction = 0.516631 Min(alpha1) = -1.54198e-21 Max(alpha1) = 1.1454
DICPCG: Solving for p_rgh, Initial residual = 0.00254906, Final residual = 8.62831e-05, No Iterations 3
time step continuity errors : sum local = 8.62297e-05, global = -2.64316e-08, cumulative = -2.53353e-05
DICPCG: Solving for p_rgh, Initial residual = 0.00031195, Final residual = 1.17906e-05, No Iterations 28
time step continuity errors : sum local = 1.61631e-05, global = 1.28852e-06, cumulative = -2.40467e-05
DICPCG: Solving for p_rgh, Initial residual = 0.000107948, Final residual = 9.00722e-09, No Iterations 82
time step continuity errors : sum local = 1.41569e-08, global = 1.41857e-09, cumulative = -2.40453e-05
GAMG: Solving for Qhi, Initial residual = 0.00241765, Final residual = 8.19049e-05, No Iterations 1
ExecutionTime = 87.7 s ClockTime = 89 s

Courant Number mean: 0.0534621 max: 0.428816
Interface Courant Number mean: 0.00291851 max: 0.357529
deltaT = 0.000111111
Time = 0.00222222

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so"
#6 Foam::viscosityModels:articleVolumFraction::calc Nu() const in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/lib/libuserparticleVolumFraction.so"
#7 Foam::viscosityModels:articleVolumFraction::corr ect() in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/lib/libuserparticleVolumFraction.so"
#8 Foam::twoPhaseMixture::calcNu() in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so"
#9
in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/bin/interFoamParticle"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/bin/interFoamParticle"
Floating point exception (core dumped)


I would appreciate any help
Mahyar

[ssss]

So as you can see in the log, there is an error in the calcNu() function, when you call pow(). Try searching in your code if there can be a divition by zero or if there is a sqrt which can be imaginary

[Mahyar Javidi]

Thanks for the reply. I will check the code and let you know what was the issue.

Regards
Mahyar

[Mahyar Javidi]

Hi

The problem was division by the zero. Solver is running perfectly now. Thanks for the help.

Mahyar

[ssss]

You are welcome

[Kittipos]

Hi

I used rhoPorousSimpleFoam to running simple case but it has an error

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #4 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #5 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" #6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #7 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" Floating point exception (core dumped)

Could you explain me pls

[ltiago94]

Hello everyone,

How do you think I can find if in my case it is dividing by zero?

I'm using interFoam and when I run it, immediatly I get this error.

Do you know how to solve it?

Thanks
Luís