CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solver Error before it even Runs

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 1, 2010, 02:45
Default Solver Error before it even Runs
  #1
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 7
trex930 is on a distinguished road
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
trex930 is offline   Reply With Quote

Old   July 1, 2010, 04:07
Default
  #2
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
im guessing here...p=0?
niklas is offline   Reply With Quote

Old   July 2, 2010, 04:24
Default
  #3
New Member
 
Lasse
Join Date: Jun 2010
Posts: 5
Rep Power: 7
lars. is on a distinguished road
Quote:
Originally Posted by trex930 View Post
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.
lars. is offline   Reply With Quote

Old   July 2, 2010, 04:30
Default
  #4
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 7
trex930 is on a distinguished road
Quote:
Originally Posted by lars. View Post
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.
trex930 is offline   Reply With Quote

Old   July 2, 2010, 04:32
Default
  #5
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
that would help, yes
niklas is offline   Reply With Quote

Old   July 2, 2010, 04:39
Default
  #6
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 7
trex930 is on a distinguished road
Quote:
Originally Posted by niklas View Post
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
trex930 is offline   Reply With Quote

Old   July 2, 2010, 04:44
Default
  #7
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
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.
niklas is offline   Reply With Quote

Old   July 2, 2010, 04:52
Default
  #8
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 7
trex930 is on a distinguished road
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? I've been stuck on this guy for awhile...
trex930 is offline   Reply With Quote

Old   July 2, 2010, 04:58
Default
  #9
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
You cant have a calculated bc, calculated based on what?

either set outlet to

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

or zeroGradient
niklas is offline   Reply With Quote

Old   July 2, 2010, 05:04
Default
  #10
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 7
trex930 is on a distinguished road
aww man... im so dumb..haha

Thanks a bunch!
trex930 is offline   Reply With Quote

Old   April 25, 2011, 23:23
Default
  #11
Member
 
Join Date: Nov 2010
Posts: 50
Rep Power: 6
usergk is on a distinguished road
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 ***
usergk is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25


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