Modelling Burner Combustion in a Reactor Network - How to integrate a PFR? (in Python)

601 views
Skip to first unread message

T Berry

unread,
Oct 25, 2016, 10:22:52 AM10/25/16
to Cantera Users' Group
Hi everyone,

I am trying to model the combustion of methane in a shell boiler. I have successfully modeled the combustion using the combustor.py example initially (I thought this to be more suitable than the flame examples as I can input fuel and air separately). I have then added a re-circulation from exhaust back into flame to model the recirculation zones in the furnace (10%).


MY QUESTION is how can I add a plug flow reactor (PFR) onto the exhaust reservoir in order to model the post flame zone, where the rest of the CO is oxidised and NOx continues to be formed. 

I have tried using the FlowReactor, however this cannot be used with other reactors and I have tried incorporating the PFR model in surf_pfr.py, however I cannot figure out which bits I need to include as its modeled around a catalyst surface.

Please find code attached.

Please note: I'm very new to cantera and python, and only just getting used to the way the language works.

Many thanks

Tom
CH4_combustor (with recirc).py

walterm...@me.com

unread,
Oct 28, 2016, 8:24:33 AM10/28/16
to Cantera Users' Group
Dear Tom,

During my PhD I used extensively Cantera for several purposes, including the simulation of  PFR in homogeneous media.
I adapted the example from surface_pfr.py into a class which you may find in my scripts stored at
https://github.com/waltermateriais/xutils/blob/master/xpython/xchemistry/xcantera/flowreactor.py

Please notice that this script uses some functions of my personal package "xutils" and are not fully documented yet. Please find attached an example application as I used in my thesis (I tried to comment, read the code above or contact me for additional details/suggestions). As soon as I pass my defense I will document the whole package and 'tell the community' about its existence. I checked that the current version under development of Cantera will be issued with a advance_to_steady_state function, which will simplify my routine.

Walter
example-pfr.py

mactone

unread,
Jun 7, 2022, 9:58:36 AM6/7/22
to Cantera Users' Group
This thread is long time ago. I can't find any code from github. I wonder if Walter can share it again?!

Best, 
Weidong

walterm...@me.com 在 2016年10月28日 星期五下午8:24:33 [UTC+8] 的信中寫道:
Reply all
Reply to author
Forward
0 new messages