[trex930]

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

[niklas]

im guessing here...p=0?

[lars.]

Originally Posted by trex930

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.

[trex930]

Originally Posted by lars.

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.

[niklas]

that would help, yes

[trex930]

Originally Posted by niklas

that would help, yes

ok here it is! Cheers, thanks in advance

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

[niklas]

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.

[trex930]

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

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

internalField   uniform 1.0e5;
boundaryField
{
    Walls
    {
        type            zeroGradient;
    }
    BigInlet
    {
        type            zeroGradient;
    }
    SmallInlet
    {
        type            zeroGradient;
    }
    Outlet
    {
        type            calculated;
    }
}

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

[niklas]

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

either set outlet to

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

or zeroGradient

[trex930]

aww man... im so dumb..haha

Thanks a bunch!

[usergk]

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:rintStack(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:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >:mega(Foam::Reaction<Foam::sutherlandTranspor t< Foam::specieThermo<Foam::janafThermo<Foam:erfect 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:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >::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 ***