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

reactingFoam: ODE Solver KRR4 "Too many integration steps"

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 6, 2011, 06:27
Default reactingFoam: ODE Solver KRR4 "Too many integration steps"
  #1
New Member
 
Paul
Join Date: Mar 2009
Posts: 19
Rep Power: 9
pajofego is on a distinguished road
Hello all,

I've been trying to do some simulations with reactingFoam. I've also used different ODE Solver (SIBS, RK, KRR4) to see if it's possible to get a stable simulation.

Both Solver (RK, KRR4) stops with following error message:

Code:
Courant Number mean: 0.0104708 max: 0.394021
deltaT = 3.1361e-05
Time = 0.00293728

Solving chemistry


--> FOAM FATAL ERROR: 
Too many integration steps

    From function ODESolver::solve(const ODE& ode, const scalar  xStart, const scalar xEnd,scalarField& yStart, const scalar eps,  scalar& hEst) const
    in file ODESolvers/ODESolver/ODESolver.C at line 107.

FOAM exiting
The maximum integration steps is hard coded as a constant variable in the source file. Could it make sense to keep it free for the user?

The SIBS Solver stops with following error message:
Code:
Solving chemistry


--> FOAM FATAL ERROR: 
attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000;  T = 187.14

    From function janafThermo<equationOfState>::checkT(const scalar T) const
    in file /home/opencfd/OpenFOAM/OpenFOAM-1.7.x/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2  Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> >::H(double) const in "/opt/openfoam171/lib/linuxGccDPOpt/libreactionThermophysicalModels.so"
#3  Foam::ODEChemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > >::solve(double, double) in "/opt/openfoam171/lib/linuxGccDPOpt/libchemistryModel.so"
#4  
 in "/opt/openfoam171/applications/bin/linuxGccDPOpt/reactingFoam"
#5  __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6  
 in "/opt/openfoam171/applications/bin/linuxGccDPOpt/reactingFoam"
Aborted
With my setup I wasn't able to get my Simulation go forward. Maybe someone has some marks or suggestions how I can get a stable ODE Solver for my reactinFoam simulation.

Thanks in advance.
pajofego

*********************************chemestryProperti es**********************************
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format binary;
class dictionary;
location "constant";
object chemistryProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

psiChemistryModel ODEChemistryModel<gasThermoPhysics>;

chemistry on;

chemistrySolver ode;
//chemistrySolver EulerImplicit;
//chemistrySolver sequential;

initialChemicalTimeStep 1e-07;

turbulentReaction on;

sequentialCoeffs
{
cTauChem 0.001;
equilibriumRateLimiter off;
}

EulerImplicitCoeffs
{
cTauChem 0.05;
equilibriumRateLimiter off;
}

odeCoeffs
{
ODESolver KRR4; //RK; //KRR4;
eps 0.05;
scale 1;
}

Cmix Cmix [ 0 0 0 0 0 0 0 ] 0.1;

// ************************************************** *********************** //
pajofego is offline   Reply With Quote

Old   January 7, 2011, 07:24
Default
  #2
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 8
nakul is on a distinguished road
Hi,

These error are coming because your solution begins to diverge. I don't think its a good idea to change the maximum number of integration steps without properly understanding its theoretical background.

Rather its better if you make changes in your test case and try to run it completely with SIBS and then see if you still get any error with KRR4.

You may check your BC as they are a major reason for causing the JANAF error you are getting with SIBS.
nakul is offline   Reply With Quote

Old   February 18, 2013, 22:40
Default
  #3
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 6
conceptone is on a distinguished road
Hi Nakul, I also meet the "too many integration step" problem, can you specify the possibility of the BC or is there any previous post discussing this?
Thanks
conceptone is offline   Reply With Quote

Old   May 18, 2013, 15:15
Default
  #4
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 6
conceptone is on a distinguished road
Hi,pajofego
I confront this problem too, did you find the reason eventually?
Thanks
conceptone 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
ODE solver coefficients mrangitschdowcom OpenFOAM Running, Solving & CFD 2 August 18, 2013 22:16
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
ODE solver extension to RKtime integration for CFD luca_g OpenFOAM Running, Solving & CFD 4 September 26, 2007 06:44
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
Need ODE solver Levi Main CFD Forum 3 April 9, 2003 22:38


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