(20g) Catalyst-Scale Modeling of Fischer-Tropsch Synthesis from Syngas Using Comsol Multiphysics
AIChE Annual Meeting
2014
2014 AIChE Annual Meeting
Catalysis and Reaction Engineering Division
Multiscale Modeling
Sunday, November 16, 2014 - 5:30pm to 5:50pm
CATALYST-SCALE MODELING OF FISCHER-TROPSCH SYNTHESIS FROM SYNGAS USING COMSOL MULTIPHYSICS
Arvind Nanduri
Department of Natural Gas Engineering
Texas A&M University-Kingsville
Patrick L. Mills
Department of Chemical and Natural Gas Engineering
Texas A&M University-Kingsville Introduction. In the early 1920â??s, Gas-To-Liquids (GTL) and Coal-To-Liquids (CTL)
technologies were developed to account for the depleting crude oil resources. During this period,
it was Franz Fischer and Hans Tropsch who developed a process to convert synthesis gas (syn gas), derived from coal gasification, to a wide range of high value- addedproducts. This process later came to be known as Fischer-Tropsch (F-T) synthesis. F-T synthesis was an experimental success but its economic viability became a topic of concern as refining of crude oil was a more developed and an economically attractive option [1]. The syn gas used in F-T synthesis can be derived from various feed stocks like coal, natural gas, biomass, and waste. This indirect liquefaction process, also termed as feed-to-liquids, is often referred to as XTLs (Where X: C=coal, G=natural gas, B=biomass, W=waste) [2].
Multi-Tubular Fixed Bed Reactors (MTFBR) and Slurry Bubble Colum Reactors (SBCR) are widely employed for FTS. An MTFBR is similar to a shell and tube heat exchanger with a catalytic reaction taking place on the tube-side and a typical MTFBR contains about 10 to 50,000 tubes. A coolant , generally water, flows on the shell-side to maintain isothermal conditions in the reactor. To model such a system, detailed knowledge about shell-side fluid solid interactions
coupled with tube-side fluid-solid transport kinetics is required. In this study, attention is focused
on modeling of diffusion and nonisothermal reaction in FT catalyst particles since this analysis is a fundamental underpinning for any reactor-scale model. The isothermal 1-D pellet model will serve as the starting point to understand the computational limitations that could be encountered while simulating the non-isothermal case.
The primary objective of this study is three-fold: (1) to simulate an isothermal 1-D catalyst pellet model; (2) use the isothermal pellet model results as a starting point for the non-isothermal case ; and (3) use extrusion coupling in COMSOL Mmultiphysics to link the catalyst pellet (2-D
domain) to a single reactor tube (1-D domain)
Methods. A micro kinetic olefin readsorption model proposed by Wang et al. (2003) for an iron-based catalyst and a macro kinetic model used by Pangarkar et al. (2009) for a cobalt-based catalyst are used to study the coupled diffusion, reaction and transport effects using COMSOL Mmultiphysics. The kinetic equations and parameters are given in Table 1. Transport of diluted species model in COMSOL is used to define the physics involved in thecoupled to catalytic
reaction in the catalyst particle.and transport of the specie.
Table 1: Kinetic Equations and Parameters
Micro Kinetic Olefin Readsorption Model on Fe-based
Catalyst (Wang et al. 2003)
Macro Kinetic model on a Co-based
Catalyst (Pangarkar et al. 2009)
(n =1)
(n â?¥ 2)
(n â?¥ 2)
Where: k0 = 22.1 and Î?E = 19.822 J/mol
Table 1: Continued
Where: n â?¥ 1
Results and Discussion. The concentration profiles of the species in a typical Fe-based and Co- based catalyst pellet are shown in Figure 1 and and Figure 2. Additional results for higher pellet radius under typical process conditions will be discussed and presented. The vapor Liquid equilibrium (VLE) calculations using Soave-Redlich-Kwong (SRK) equation of state will also be presented to identify the phase of the species.
Figure 1: Concentration Profiles of Species in a Typical Fe-based Catalyst Pellet.
Rp (mm) Rp (mm)
Rp (mm)
Figure 2: Concentration Profiles of Species in a Typical Co-based Catalyst Pellet.
It is shown that to properly model the FT system and capture the catalyst performance, detailed models for the transport kinetic interactions that define the yield of various hydrocarbon products are required.
References
1. D. A. Wood, C. Nwaoha and B. F. Towler, "Gas-to-Liquids (GTL): A Review of an Industry Offering Several Routes for monitizing Natural Gas," Journal of Natural Gas Science and Engineering, vol. 9, 2012.
Formatted: Line spacing: 1.5 lines
Formatted: Font: (Default) Times New
Roman, 12 pt
Formatted: Justified, Indent: First line: 0.25", Line spacing: 1.5 lines
2. O. O. James, B. Chowdhury, M. Adediran Mesubi and S. Maity, "Reflections on the
Chemistry of the Fischer-Tropsch Synthesis," RSC Advances, vol. 2, no. 19, 2012.
3. Y. N. Wang, W. P. Ma, Y. J. Lu, J. Yang, Y. Y. Xu, H. W. Xiang, Y. W. Li, Y. L. Zhao and B. J. Zhang, â?? Kinetic Modelling of Fischer-Tropsch Synthesis over an Industrial Fe-Cu-K catalyst,â? Fuel, 2003.
4. K. Pangarkar, T. J. Schildhauer, J. Ruud Van Ommen, J. Nijenhuis, J. A. Moulijn and F.
Kapteijn, â??Experimental and Numerical Comparison of Structured Packings with a
Randomly Packed Bed Reactor for Fischer-Tropsch Synthesis,â? Catalysis Today, 2009.
Topics
Checkout
This paper has an Extended Abstract file available; you must purchase the conference proceedings to access it.
Do you already own this?
Log In for instructions on accessing this content.
Pricing
Individuals
AIChE Pro Members | $150.00 |
AIChE Graduate Student Members | Free |
AIChE Undergraduate Student Members | Free |
AIChE Explorer Members | $225.00 |
Non-Members | $225.00 |