CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Solver Error before it even Runs (http://www.cfd-online.com/Forums/openfoam-solving/77699-solver-error-before-even-runs.html)

trex930 July 1, 2010 02:45

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

niklas July 1, 2010 04:07

im guessing here...p=0?

lars. July 2, 2010 04:24

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.

trex930 July 2, 2010 04:30

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.

niklas July 2, 2010 04:32

that would help, yes

trex930 July 2, 2010 04:39

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

niklas July 2, 2010 04:44

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 July 2, 2010 04:52

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...

niklas July 2, 2010 04:58

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 July 2, 2010 05:04

aww man... im so dumb..haha

Thanks a bunch!

usergk April 25, 2011 23:23

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 ***


All times are GMT -4. The time now is 18:32.