CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development at line 0.

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

LinkBack Thread Tools Display Modes
Old   November 18, 2014, 01:49
Default at line 0.
Join Date: Jul 2013
Posts: 55
Rep Power: 5
Thamali is on a distinguished road
Dear Foamers,

I am solving a case using an Euler-Euler solver developed for combustion of wood chips.
When I run the case using the solver,it stops at calculating "rCH4" and FPE error appears.
#0  Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OSspecific/POSIX/printStack.C:221
#1  Foam::sigFpe::sigHandler(int) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OSspecific/POSIX/signals/sigFpe.C:117
#2   in "/lib/x86_64-linux-gnu/"
#3   in "/lib/x86_64-linux-gnu/"
#4  Foam::pow(double, double) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/doubleFloat.H:78
#5  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/fields/Fields/scalarField/scalarField.C:118 (discriminator 2)
#6  void Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:274
 at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:328
 at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:350 (discriminator 1)
 at ~/OpenFOAM/thamali-2.2.2/applications/solvers/latest/testonlymoistandcharnpyrovolat/UEqn1.H:95 (discriminator 1)
#10  __libc_start_main in "/lib/x86_64-linux-gnu/"
 at ??:?
Floating point exception (core dumped)
But,when I check the parameters YCH4,YO2 or tg of previous time step necessary for calculation,all are positive values.
When I,unset FOAM_SIGFPE,
The "rCH4" become negative at 2 points(1 point in a boundary and 1point in internalField).
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/thamali/OpenFOAM/thamali-2.2.2/run/tutorials/combustion/woodchipcombustion/new/latest/withwallradiation/test/drypyrochacovolat/system/ at line 0.

    From function operator>>(Istream&, Scalar&)
    in file lnInclude/Scalar.C at line 91.

FOAM exiting
I checked case log and for tg;
smoothSolver:  Solving for tg, Initial residual = nan, Final residual = nan, No Iterations 1000
I tried by changing "fvSchemes " and they did not work.
I am attaching the files calculating "rCH4" and file calculating Gasphase mas fractions(Y) and tg.

I highly appreciate your help.
If you want to look at other files of the solver,I am happy put upload.

Thanks in advance.
Attached Files
File Type: h YghgEqn.H (4.1 KB, 1 views)
File Type: h UEqn1.H (10.9 KB, 2 views)
Thamali is offline   Reply With Quote


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
Problem with Gmsh nishant_hull Open Source Meshers: Gmsh, Netgen, CGNS, ... 23 August 5, 2015 02:09
Annoying issue of automatic "Rescale to Data Range " with paraFoam/paraview 3.12 keepfit OpenFOAM Paraview & paraFoam 60 September 18, 2013 03:23
OpenFOAM 1.7.1 installation problem on OpenSUSE 11.3 flakid OpenFOAM Installation 16 December 28, 2010 09:48
Regarding FoamX running Kindly help out hariya03 OpenFOAM Pre-Processing 0 April 18, 2008 04:26
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00

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