CFD Online Discussion Forums - Solver Error before it even Runs

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Solver Error before it even Runs (https://www.cfd-online.com/Forums/openfoam-solving/77699-solver-error-before-even-runs.html)

Solver Error before it even Runs

Hello Foamers!

I am trying to run a compressible LES case. However when I run the solver, it doesnt even go through a single iteration and it gives me the terrible, floating point exception / SigFpe error. Im pretty sure my BC are set right. But can anyone shed any light on such a problem. Or can someone post up their BC for an compressible LES model so I can compare?

I am trying to simulate a mixer just for your information!

Cheers,
T

im guessing here...p=0?

Quote:

Originally Posted by trex930 (Post 265229)

it gives me the terrible, floating point exception / SigFpe error. Im pretty sure my BC are set right.

Thank you very much for this very detailed error report.
I am pretty sure the SigFpe error is right.

Quote:

Originally Posted by lars. (Post 265394)

Thank you very much for this very detailed error report.
I am pretty sure the SigFpe error is right.

Would you like a copy of the error? I would gladly put it up if you are able to help me out.

that would help, yes

Quote:

Originally Posted by niklas (Post 265397)

that would help, yes

ok here it is! Cheers, thanks in advance

Code:

Creating field DpDt

#0  Foam::error::printStack(Foam::Ostream&) in "/usr/OpenFOAM/OpenFOAM-1.6/lib/l

inux64GccDPOpt/libOpenFOAM.so"

#1  Foam::sigFpe::sigFpeHandler(int) in "/usr/OpenFOAM/OpenFOAM-1.6/lib/linux64G

ccDPOpt/libOpenFOAM.so"

#2  ?? in "/lib64/libc.so.6"

#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<d

ouble> const&) in "/usr/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so

"

#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricFie

ld<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double

, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, F

oam::fvsPatchField, Foam::surfaceMesh> const&) in "/usr/OpenFOAM/OpenFOAM-1.6/ap

plications/bin/linux64GccDPOpt/rhoPisoFoam"

#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMes

h> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricFiel

d<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::Geomet

ricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/usr/Open

FOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoPisoFoam"

#6  main in "/usr/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoPiso

Foam"

#7  __libc_start_main in "/lib64/libc.so.6"

#8  __gxx_personality_v0 in "/usr/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64

GccDPOpt/rhoPisoFoam"

Floating exception

just a heads up. The same boundary condition worked for the same model but different mesh...i know that should help me conclude that its a mesh problem. But if there was a mesh problem why wouldnt ccm26ToFoam catch it?

Thanks again
T

you have a division by zero in your bc or in your field.

It usually is epsilon, or p.
many forget to set p to 1.0e+5 (or the actual pressure) from 0 when going from incompressible to compressible.

thats what I thought as well. But the code ran properly with the other model. all I did was change the mesh. Heres epsilon

Code:

boundaryField

{

    Walls

    {

        type            compressible::epsilonWallFunction;

        value           uniform 200;

    }

    BigInlet

    {

        type            compressible::turbulentMixingLengthDissipationRateInlet;

        mixingLength    0.150;

        value           uniform 200;

    }

    SmallInlet

    {

        type            compressible::turbulentMixingLengthDissipationRateInlet;

        mixingLength    0.0475;

        value           uniform 200;

    }

    Outlet

    {

        type            inletOutlet;

        inletValue      uniform 200;

        value           uniform 200;

    }

}

and p

Code:

internalField   uniform 1.0e5;

boundaryField

{

    Walls

    {

        type            zeroGradient;

    }

    BigInlet

    {

        type            zeroGradient;

    }

    SmallInlet

    {

        type            zeroGradient;

    }

    Outlet

    {

        type            calculated;

    }

}

Any clue?:confused: I've been stuck on this guy for awhile...

You cant have a calculated bc, calculated based on what?

either set outlet to

outlet
{
type fixedValue;
value uniform 1.0e+5;
}

or zeroGradient

aww man... im so dumb..haha

Thanks a bunch!

Hello,

I am trying to implement combustion in OpenFOAM and I get the error below. This happens mid-way during the simulation, after a few time steps.

Any idea why this could be occurring? I am relatively new to using OpenFOAM, so any information would be very much appreciated.

Thanks so much!
gk

[24] #0 Foam::error::printStack(Foam::Ostream&) in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[24] #1 Foam::sigFpe::sigFpeHandler(int) in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[24] #2 in "/lib/libc.so.6"
[24] #3 in "/lib/libm.so.6"
[24] #4 pow in "/lib/libm.so.6"
[24] #5 Foam::ODEChemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::perfectGas> > > >::omega(Foam::Reaction<Foam::sutherlandTranspor t< Foam::specieThermo<Foam::janafThermo<Foam::perfect Gas> > > > const&, Foam::Field<double> const&, double, double, double&, double&, int&, double&, double&, int&) const in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so"
[24] #6 Foam::ODEChemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::perfectGas> > > >::tc() const in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so"
[24] #7
[24] in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/pFoam"
[24] #8 __libc_start_main in "/lib/libc.so.6"
[24] #9
[24] in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/pFoam"
[node69:00818] *** Process received signal ***
[node69:00818] Signal: Floating point exception (8)
[node69:00818] Signal code: (-6)
[node69:00818] Failing at address: 0x58f800000332
[node69:00818] [ 0] /lib/libc.so.6(+0x33af0) [0x2b0cdf0abaf0]
[node69:00818] [ 1] /lib/libc.so.6(gsignal+0x35) [0x2b0cdf0aba75]
[node69:00818] [ 2] /lib/libc.so.6(+0x33af0) [0x2b0cdf0abaf0]
[node69:00818] [ 3] /lib/libm.so.6(+0x13e81) [0x2b0cdebf0e81]
[node69:00818] [ 4] /lib/libm.so.6(pow+0x15) [0x2b0cdec02765]
[node69:00818] [ 5] /home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so(_ZNK4Foam17ODEChemistryModelI NS_17psiChemistryModelENS_19sutherlandTransportINS _12specieThermoINS_11janafThermoINS_10perfectGasEE EEEEEE5omegaERKNS_8ReactionIS8_EERKNS_5FieldIdEEdd RdSI_RiSI_SI_SJ_+0x285) [0x2b0cdd978ff5]
[node69:00818] [ 6] /home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so(_ZNK4Foam17ODEChemistryModelI NS_17psiChemistryModelENS_19sutherlandTransportINS _12specieThermoINS_11janafThermoINS_10perfectGasEE EEEEEE2tcEv+0x57e) [0x2b0cdd98424e]
[node69:00818] [ 7] pFoam() [0x426bf3]
[node69:00818] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x2b0cdf096c4d]
[node69:00818] [ 9] pFoam() [0x421119]
[node69:00818] *** End of error message ***