CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   CVode error while running (https://www.cfd-online.com/Forums/openfoam/87275-cvode-error-while-running.html)

NewKid April 15, 2011 10:25

CVode error while running
 
I got error messgae below while running Cantera solved CH4 combustion case. Anyone can tell me the possible reasons cause it?

Error message:

Time = 3
Solving chemistry
[CVODES ERROR] CVode
At t = 0 and h = 3.55081e-29, the error test failed repeatedly or with |h| = hmin.

************************************************
Cantera Error!
************************************************

Procedure: CVodesIntegrator
Error: CVodes error encountered.

"H2":1e-15 "H":9.97698e-16 "O":1e-15 "O2":0.21 "OH":1e-15 "H2O":1e-15 "HO2":1.07538e-15 "H2O2":1e-15 "C":1e-15 "CH":1e-15 "CH2":1e-15 "CH2(S)":1e-15 "CH3":1e-15 "CH4":1e-15 "CO":1e-15 "CO2":1e-15 "HCO":1e-15 "CH2O":1e-15 "CH2OH":1e-15 "CH3O":1e-15 "CH3OH":1e-15 "C2H":1e-15 "C2H2":1e-15 "C2H3":1e-15 "C2H4":1e-15 "C2H5":1e-15 "C2H6":1e-15 "HCCO":1e-15 "CH2CO":1e-15 "HCCOH":1e-15 "N":1e-15 "NH":1e-15 "NH2":1e-15 "NH3":1e-15 "NNH":1e-15 "NO":1e-15 "NO2":1e-15 "N2O":1e-15 "HNO":1e-15 "CN":1e-15 "HCN":1e-15 "H2CN":1e-15 "HCNN":1e-15 "HCNO":1e-15 "HOCN":1e-15 "HNCO":1e-15 "NCO":1e-15 "N2":0.78 "AR":0.01 "C3H7":1e-15 "C3H8":1e-15 "CH2CHO":1e-15 "CH3CHO":1e-15
With the state: T:298 p:3.14616e+14 and the mixture
53
(
1e-15
1e-15
1e-15
0.21
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
1e-15
0.78
0.01
1e-15
1e-15
1e-15
1e-15
)
Cantera complained in cell 639
#0 Foam::error::printStack(Foam::Ostream&) in "/OpenFOAM/OpenFOAM-1.5-dev/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/OpenFOAM/OpenFOAM-1.5-dev/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) at /OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/errorManip.H:86
#3 Foam::canteraLocalTimeChemistryModel::solve(double , double) at ~/OpenFOAM/slax-1.5-dev/Libraries/canteraThermosChemistry/canteraLocalTimeChemistryModel.C:132
#4 main at /OpenFOAM/OpenFOAM-1.5-dev/applications/solvers/AlternateChemistry/Steady/alternateSteadyReactingFoam/chemistry.H:9
#5 __libc_start_main in "/lib/libc.so.6"
#6 _start in "/home/slax/OpenFOAM/slax-1.5-dev/applications/bin/linuxGccDPOpt/alternateSteadyReactingFoam"

From function canteraLocalTimeChemistryModel::solve
in file canteraLocalTimeChemistryModel.C at line 127.
FOAM aborting (FOAM_ABORT set)
Aborted

gschaider April 18, 2011 09:12

Quote:

Originally Posted by NewKid (Post 303791)
I got error messgae below while running Cantera solved CH4 combustion case. Anyone can tell me the possible reasons cause it?

Error message:

Time = 3
Solving chemistry
[CVODES ERROR] CVode
At t = 0 and h = 3.55081e-29, the error test failed repeatedly or with |h| = hmin.

************************************************
Cantera Error!
************************************************

Procedure: CVodesIntegrator
Error: CVodes error encountered.

"H2":1e-15 "H":9.97698e-16 "O":1e-15 "O2":0.21 "OH":1e-15 "H2O":1e-15 "HO2":1.07538e-15 "H2O2":1e-15 "C":1e-15 "CH":1e-15 "CH2":1e-15 "CH2(S)":1e-15 "CH3":1e-15 "CH4":1e-15 "CO":1e-15 "CO2":1e-15 "HCO":1e-15 "CH2O":1e-15 "CH2OH":1e-15 "CH3O":1e-15 "CH3OH":1e-15 "C2H":1e-15 "C2H2":1e-15 "C2H3":1e-15 "C2H4":1e-15 "C2H5":1e-15 "C2H6":1e-15 "HCCO":1e-15 "CH2CO":1e-15 "HCCOH":1e-15 "N":1e-15 "NH":1e-15 "NH2":1e-15 "NH3":1e-15 "NNH":1e-15 "NO":1e-15 "NO2":1e-15 "N2O":1e-15 "HNO":1e-15 "CN":1e-15 "HCN":1e-15 "H2CN":1e-15 "HCNN":1e-15 "HCNO":1e-15 "HOCN":1e-15 "HNCO":1e-15 "NCO":1e-15 "N2":0.78 "AR":0.01 "C3H7":1e-15 "C3H8":1e-15 "CH2CHO":1e-15 "CH3CHO":1e-15
With the state: T:298 p:3.14616e+14 and the mixture

CVode/Cantera can't satisfy the requested accuracy of the solution. That can have a lot of causes. But I think that a pressure of 3e14 is a bit excessive. Seems to me the solution already blew up at the previous timestep. I'd suggest writting out all timesteps. Having a look where it blew up. Think what could be the cause. Play around with intitial values, relaxation factors, check the chemical system.

NewKid April 21, 2011 04:46

gschaider, seems that I've resolved this problem.
The main reason seems to be the unstructed mesh, now I use one mesh made by blockMesh, everything seems OK.
And I changed one initial condition of p from "totalPressure" to "fixedValue", and T from "totalTemprature" to "fixedValue".


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