CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Methane combustion

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 2 Post By Micael
  • 1 Post By sberg

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 5, 2005, 06:51
Default Methane combustion
  #1
Steve
Guest
 
Posts: n/a
Has anyone tried the Jones and Lindstedt four equation mechanism for methane combustion in Fluent?

I have implemented the four steps in the "propdb.scm" file. It works fine for the EDM, but when EDC is activated Fluent gives the error message:

********** ISAT_ABORT ************** routine = phi2y location = 1 message = non-positive species sum

****** END ISAT_ABORT **************

Error: isat terminated. Error Object: #f

Does anyone know how to implement this mechanism in Fluent so the EDC works?

Steve
  Reply With Quote

Old   November 17, 2011, 22:26
Default Hi
  #2
New Member
 
San Siro
Join Date: Nov 2011
Posts: 7
Rep Power: 14
San Siro is on a distinguished road
Quote:
Originally Posted by Steve
;123005
Has anyone tried the Jones and Lindstedt four equation mechanism for methane combustion in Fluent?

I have implemented the four steps in the "propdb.scm" file. It works fine for the EDM, but when EDC is activated Fluent gives the error message:

********** ISAT_ABORT ************** routine = phi2y location = 1 message = non-positive species sum

****** END ISAT_ABORT **************

Error: isat terminated. Error Object: #f

Does anyone know how to implement this mechanism in Fluent so the EDC works?

Steve
Hi Steve,I encountered the same question just like you did.I'm wondering if you have solved this problem.Looking forward to you reply.
San Siro is offline   Reply With Quote

Old   November 18, 2011, 08:07
Default Hi San Siro
  #3
New Member
 
Steve Berg
Join Date: Mar 2009
Posts: 9
Rep Power: 17
sberg is on a distinguished road
At that time I contacted Fluent support and they said I had to use a Chemkin file. They even provided me the file. It's named 4_reactions.che and I attached it below. As far as I remembered the Chemkin file ran smoothly with the EDC model in Fluent 5.3 and 6. I hope it will work for you as well.

Kind regards
Steve
----------------

ELEMENTS
H O C N
END
SPECIES
CH4 O2 CO2 CO H2O H2 N2
END
THERMO
END
REACTIONS KJOULES/MOLE
!Units are cm, mol, s, kJ/mol and K.
2CH4+O2=>2CO+4H2 7.824435194e13 0.0 0.12552
FORD/CH4 0.5/O2 1.25/
CH4+H2O=>CO+3H2 3.1E11 0.0 0.12552
CO+H2O=>CO2+H2 2.75E12 0.0 0.08368
2H2+O2<=>2H2O 4.44566E21 -1.0 0.16736
FORD/H2 0.5/O2 2.25/
RORD/H2O -1.0/
END
sberg is offline   Reply With Quote

Old   November 22, 2011, 23:03
Default HI
  #4
New Member
 
San Siro
Join Date: Nov 2011
Posts: 7
Rep Power: 14
San Siro is on a distinguished road
When I used the laminal finite-rate model to simulate a reaction involves both volumetric reactions and wall surface reactions,I selected the Two-Step Solution Procedure.I begin by solving the flow, energy, and species equations with reactions disabled.When the basic flow pattern has been established,I enable the reactions and continue the calculation.
Then fluent gives the error message:

Error: Non-physical Specific Heat = 0.000000e+000.
Error Object: ()

I ignore the message and continue to calculate,fluent gives the following error message:

Error: isat terminated.
Error Object: ().
This problem has troubled me for a long time.Is there any suggestions how can I solve this problem and continue my calculation.
Thank
San Siro
San Siro is offline   Reply With Quote

Old   March 9, 2012, 17:04
Default
  #5
Senior Member
 
Micael
Join Date: Mar 2009
Location: Canada
Posts: 157
Rep Power: 18
Micael is on a distinguished road
Quote:
Originally Posted by sberg View Post
At that time I contacted Fluent support and they said I had to use a Chemkin file. They even provided me the file. It's named 4_reactions.che and I attached it below. As far as I remembered the Chemkin file ran smoothly with the EDC model in Fluent 5.3 and 6. I hope it will work for you as well.

Kind regards
Steve
----------------

ELEMENTS
H O C N
END
SPECIES
CH4 O2 CO2 CO H2O H2 N2
END
THERMO
END
REACTIONS KJOULES/MOLE
!Units are cm, mol, s, kJ/mol and K.
2CH4+O2=>2CO+4H2 7.824435194e13 0.0 0.12552
FORD/CH4 0.5/O2 1.25/
CH4+H2O=>CO+3H2 3.1E11 0.0 0.12552
CO+H2O=>CO2+H2 2.75E12 0.0 0.08368
2H2+O2<=>2H2O 4.44566E21 -1.0 0.16736
FORD/H2 0.5/O2 2.25/
RORD/H2O -1.0/
END
How reverse order for H2O should be interpreted? It is written as -1 (minus one). What the minus sign means? I mean, it is unphysical to have a negative order, that would mean the less concentration, the fastest kinetic (and infinite kinetic at zero concentration!)

---
Edit: My guess so far is that adding a minus or plus sign mean a correction on the order. Anybody can confirm?
Micael is offline   Reply With Quote

Old   March 9, 2012, 20:05
Default
  #6
New Member
 
Steve Berg
Join Date: Mar 2009
Posts: 9
Rep Power: 17
sberg is on a distinguished road
The H2 + 0.5 O2 --> H2O reaction is a global one and it includes several elementary steps with many reactive intermediates like H, O, OH. Such global reactions have complicated rate laws. The reaction orders are usually not integers and they can be negative.

You can probably find some more details, if you are interested, in Warnatz, Maas and Dibble's textbook "Combustion."
sberg is offline   Reply With Quote

Old   March 9, 2012, 20:28
Default
  #7
New Member
 
Steve Berg
Join Date: Mar 2009
Posts: 9
Rep Power: 17
sberg is on a distinguished road
And of course, a negative reaction order means the rate decreases when the concentration of that component increases

r = k(T) [H2]^0.5 [O2]^2.25 [H2O]^(-1) = k(T) [H2]^0.5 [O2]^2.25 / [H2O]
amin_jalalian likes this.
sberg is offline   Reply With Quote

Old   March 10, 2012, 18:25
Default
  #8
Senior Member
 
Micael
Join Date: Mar 2009
Location: Canada
Posts: 157
Rep Power: 18
Micael is on a distinguished road
Thanks for your reply.

It is still not clear to me. I thought the chemkin format means:

r = kf(T) [H2]^0.5 [O2]^2.25 - kb(T) [H2O]^(-1)

kf: forward constant
kb: reverse constant

Of course, put this way, it does not make sense (as the product decrease, the reverse kinetic increase).

Looking at the file, it seems that the reaction is considered reversible. The way you put it, I don't see how it can be reversible.
Micael is offline   Reply With Quote

Old   March 10, 2012, 20:04
Default
  #9
New Member
 
Steve Berg
Join Date: Mar 2009
Posts: 9
Rep Power: 17
sberg is on a distinguished road
Yes, you are correct.

2H2+O2<=>2H2O is often simplified to be non-reversible and represented by r = k(T) [H2]^0.5 [O2]^2.25 [H2O]^(-1). But according to the Chemkin output, as you noticed, it isn't the case. The reaction seems to be reversible with a reaction order of -1.
sberg is offline   Reply With Quote

Old   April 1, 2014, 18:33
Default Different version of chemkin
  #10
New Member
 
Kevin
Join Date: Jun 2013
Posts: 4
Rep Power: 13
leexit is on a distinguished road
Hi
I just wonder if the chemkin file that is available in Fluent can be older version of chemkin. I use CHEMKIN-II Version 3.1 and do not think it allows me to enter the reaction order. Has anyone ever tried to use the Jones and Lindstedt mechanism with my version of chemkin?
Thank you in advance
leexit is offline   Reply With Quote

Old   March 17, 2021, 08:40
Default
  #11
New Member
 
Prashant Gahlot
Join Date: Mar 2021
Posts: 1
Rep Power: 0
Prashantgahlot is on a distinguished road
Quote:
Originally Posted by sberg View Post
At that time I contacted Fluent support and they said I had to use a Chemkin file. They even provided me the file. It's named 4_reactions.che and I attached it below. As far as I remembered the Chemkin file ran smoothly with the EDC model in Fluent 5.3 and 6. I hope it will work for you as well.

Kind regards
Steve
----------------

ELEMENTS
H O C N
END
SPECIES
CH4 O2 CO2 CO H2O H2 N2
END
THERMO
END
REACTIONS KJOULES/MOLE
!Units are cm, mol, s, kJ/mol and K.
2CH4+O2=>2CO+4H2 7.824435194e13 0.0 0.12552
FORD/CH4 0.5/O2 1.25/
CH4+H2O=>CO+3H2 3.1E11 0.0 0.12552
CO+H2O=>CO2+H2 2.75E12 0.0 0.08368
2H2+O2<=>2H2O 4.44566E21 -1.0 0.16736
FORD/H2 0.5/O2 2.25/
RORD/H2O -1.0/
END



Tried this in my fluent, but its showing floating point error. Can Somebody please guide me through this?
Prashantgahlot is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Temperature drops in Methane combustion CPRL Siemens 0 July 24, 2007 13:56
How to model combustion with only methane stream? ligang zheng FLUENT 2 May 1, 2007 10:16
Methane combustion in gas burner Stefano CFX 1 June 7, 2005 03:12
Presumed PDF for methane combustion Andreas Siemens 0 May 21, 2003 10:07
Methane Combustion Lars FLUENT 4 March 5, 2003 09:24


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