Next Article in Journal
Fermentative Conversion of Two-Step Pre-Treated Lignocellulosic Biomass to Hydrogen
Next Article in Special Issue
Fischer-Tropsch Synthesis: Cd, In and Sn Effects on a 15%Co/Al2O3 Catalyst
Previous Article in Journal
Development of Nickel-BTC-MOF-Derived Nanocomposites with rGO Towards Electrocatalytic Oxidation of Methanol and Its Product Analysis
Previous Article in Special Issue
The Evolution of Catalysis for Alkyd Coatings: Responding to Impending Cobalt Reclassification with Very Active Iron and Manganese Catalysts, Using Polydentate Nitrogen Donor Ligands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fischer–Tropsch Synthesis: Computational Sensitivity Modeling for Series of Cobalt Catalysts

1
Asbury University, One Macklem Drive, Wilmore, KY 40390, USA
2
Center for Applied Energy Research, University of Kentucky, 2540 Research Park Dr., Lexington, KY 40511, USA
3
Department of Biomedical Engineering and Chemical Engineering, University of Texas at San Antonio/Department of Mechanical Engineering, 1 UTSA Circle, San Antonio, TX 78249, USA
*
Author to whom correspondence should be addressed.
Catalysts 2019, 9(10), 857; https://doi.org/10.3390/catal9100857
Submission received: 22 July 2019 / Revised: 24 September 2019 / Accepted: 5 October 2019 / Published: 15 October 2019
(This article belongs to the Special Issue Iron and Cobalt Catalysts)

Abstract

:
Nearly a century ago, Fischer and Tropsch discovered a means of synthesizing organic compounds ranging from C1 to C70 by reacting carbon monoxide and hydrogen on a catalyst. Fischer–Tropsch synthesis (FTS) is now known as a pseudo-polymerization process taking a mixture of CO as H2 (also known as syngas) to produce a vast array of hydrocarbons, along with various small amounts of oxygenated materials. Despite the decades spent studying this process, it is still considered a black-box reaction with a mechanism that is still under debate. This investigation sought to improve our understanding by taking data from a series of experimental Fischer–Tropsch synthesis runs to build a computational model. The experimental runs were completed in an isothermal continuous stirred-tank reactor, allowing for comparison across a series of completed catalyst tests. Similar catalytic recipes were chosen so that conditional comparisons of pressure, temperature, SV, and CO/H2 could be made. Further, results from the output of the reactor that included the deviations in product selectivity, especially that of methane and CO2, were considered. Cobalt was chosen for these exams for its industrial relevance and respectfully clean process as it does not intrinsically undergo the water–gas shift (WGS). The primary focus of this manuscript was to compare runs using cobalt-based catalysts that varied in two oxide catalyst supports. The results were obtained by creating two differential equations, one for H2 and one for CO, in terms of products or groups of products. These were analyzed using sensitivity analysis (SA) to determine the products or groups that impact the model the most. The results revealed a significant difference in sensitivity between the two catalyst–support combinations. When the model equations for H2 and CO were split, the results indicated that the CO equation was significantly more sensitive to CO2 production than the H2 equation.

1. Introduction

Fischer–Tropsch synthesis (FTS) is a process used to produce a vast range of organic compounds. Discovered nearly a century ago, FTS reacts H2 and CO over a metal catalyst to produce C1-C70-chain compounds, including gasoline and alternate materials [1,2]. Thus, it is often used to generate fuels, monomers for common polymers, rubber, etc. Despite the long history of research and industrial use of FTS, it is still unknown how this process works on a microscopic level [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. There are at least three proposed mechanisms, namely, the carbide originally proposed by Fischer and Tropsch, CO insertion, and formate mechanisms. Yet, none of these can fully describe the entire FTS process, mainly because of the diversity of materials generated by different active metals [1].
To further comprehend FTS, many have attempted to use mathematical modeling as a useful means of describing the process [23,24,25,26,27,28]. Density functional theory (DFT) calculations, paired with experimental procedures, have provided further insight into the FTS mechanism [18,19,20,21]. These models generally fall into one of two categories—probabilistic models and kinetic models. Probabilistic models attempt to describe the propagation of the carbon chain to generate the products observed [23,24,25]. Kinetic models focus more specifically on the rate of the reaction, determining a mechanism to describe the reaction as a whole [26,27,28]. While these models have significant application for understanding FTS, they do not exist without limitation. Kinetic models depend on the mechanism of the reaction, which, as stated earlier, is not fully understood due to the diverse array of products. Furthermore, they become quite mathematically complex. Probabilistic models focus on the likelihood of a specific carbon chain propagating or terminating. This type of modeling is limited by the monomer used for propagation, which is not fully understood.
This research sought to create a model that describes the rate of change of hydrogen and carbon monoxide as the reaction progresses. Using a series of industrially relevant cobalt-catalyst runs with varying supports, a system of differential equations was created to model the unreacted syngas in terms of the products generated. Each chemical product output term (such as moles of methane, water, carbon dioxide, etc.) was scaled by an impact parameter, a value in units of hr-1. The model was then tested using sensitivity analysis (SA) to determine which impact parameter and, consequently, which associated chemical product term altered the solution of the model the most. From the results of this analysis, conclusions were drawn about the chemical implications, which may then be used to further examine catalytic runs.

2. Computational Modeling

The computational model was established using chemical engineering process principles from the stirred-tank slurry reactor and basic chemical kinetics. The rate of H2 and CO coming out of the reactor (unreacted syngas) was modeled using ordinary first-order differential equations (ODEs) [29]. More specifically, the rate of change of H2 in the reaction can be described by the hydrogen flow into the reactor minus everything coming out of the reactor (H2, water, methane, hydrocarbon products). All the products (excluding unreacted hydrogen) were scaled by an “impact parameter”. This model was applied to CO as well to produce two ODEs as follows:
d ( M H 2 ( t ) ) d t = Q H 2 q H 2 M H 2 ( t ) Z 1 M H 2 O ( t ) Z 2 M C H 4 ( t ) Z 3 M l i q ( t ) Z 9 M C 2 C 4 ( t ) [ Z 10 M C O 2 ( t ) ] *
d ( M C O ( t ) ) d t = Q C O q C O M C O ( t ) Z 4 M H 2 O ( t ) Z 5 M C O 2 ( t ) Z 6 M C H 4 ( t ) Z 7 M l i q ( t ) Z 8 M C 2 C 4 ( t )
where QH2 and QCO are input rates of H2 and CO in moles per hour, respectively; qH2 and qCO are the output rates for H2 and CO in units of hr-1, respectively; MH2(t), MCO(t), MH2O(t), MCO2(t), MCH4(t), MC2C4(t), and Mliq(t) are the time-dependent functions for the moles of H2, CO, CO2, H2O, CH4, C2–C4 products, and C5+ (or liquid) products, respectively; and Z1–10 are the impact parameters in units of hr-1. It is important to note that MH2(t) and MCO(t) are the unknown functions for this system. All other mole functions were interpolated from the experimental data (see Supplementary Materials for data), but these two are the functions for which the system is solved. The MH2(t) and MCO(t) solutions of the ODE system were compared to the experimentally observed values of these functions (see Supplementary Materials) and the resulting sum of squared error was the “response” for the sensitivity analysis. The * symbol at the end of the hydrogen equation is to note that the bracketed term was added later, thus explaining the variance in sequence of impact parameters. Originally, the carbon dioxide term was not included in the hydrogen equation but was later added to explore the possibility of a water–gas shift (WGS).
Due to the innovative nature of the differential equations modeling the reactor, the solver used to evaluate the solutions to this model was the Livermore Solver for Ordinary Differential Equations (LSODA). LSODA performs a test for numerical stiffness as the time steps evolve. A stiff differential equation is one which exhibits solutions with different time scales. For example, stiffness may occur in chemical modeling when simulating reactions with very different reaction rate constants. To solve stiff equations numerically, one must use an implicit solver that requires time-consuming nonlinear equation solutions (e.g., Newton’s method for nonlinear systems) for each time step. The LSODA solver primarily uses a non-stiff solver for speed but dynamically checks the solution to determine if a stiff method is required and applies an implicit solver for those time steps. Since the values of the impact parameters were unknown, this approach provided the most flexibility in the numerical solutions. Initial testing of the solution showed good qualitative agreement with the data, indicating that the LSODA solver worked properly with this model. Once the solution was established, sensitivity analysis (SA) was performed. SA is a branch of computational methods that provide insight into how much the total change in the output of a model can be divided up and allocated to the “input” variables/parameters of the model [30]. In other words, SA measures the effect on the output of the model that results from changing the input factors. While there are two types of SA, local and global, local was not used in this work since it involves only measures the change in output at a specified point. Global SA measures the output across a specified range for the input parameters, which was more practical in analyzing this model [30].
While there are many methods for global SA, the Sobol method and Morris method [31,32,33] were of primary interest for application to the model. These are commonly used for such systems and they are both well validated in multiple applications [32,34,35]. The Sobol method is a variance-based global SA method that generates sensitivity indices from ratios of variance terms [33]. The limitation of Sobol is that calculating the variances is very complex mathematically, and the method requires a large number of samples (often more than 10,000) to counteract any error from negative indices. Furthermore, the primary function of the Sobol method is to analyze the interactions by looking at the second- and higher-order indices, but since the primary focus of this work was the impact of a factor on the model output, the Morris method was used as the SA tool for this preliminary study.
The Morris method is known as a one-at-a-time (OAT) global SA method [31]. This method measures the impact of each input factor on the model by changing only one input factor at a time and observing the change in the model output. Specifically, Morris generates a set of elementary effects (EEs) for each input factor across the given range [31]. An EE is generated by starting with the model evaluated at a specific set of input factors. Then, the model is evaluated again with only one factor increased or decreased by a specified amount, D. The unchanged model output is subtracted from the altered model output, the absolute value is taken for the difference, and this absolute value is divided by D, as follows:
E E i = | Y ( X 1 , X 2 , , X i + D ) Y ( X 1 , X 2 , , X i ) | D
In the equation above, Y represents the model, X1…i is the set of parameters, and i = {1, 2, …, n}, n being the number of parameters in the model. A set of these EEs using multiple values for D is generated for each input parameter across the given range. Since each set, known as the distribution of EEs, can be quite large, they are sampled. The original Morris sampling technique was a type of random sampling that did not use a high number of samples (10–50), but Campolongo et al. discovered a better sampling technique [31]. In this new technique, 500–1000 original samples are taken from the distribution of EEs. These are then analyzed to determine the optimal 10–50 samples (those that most closely describe the set as a whole). The optimal samples, or optimal trajectories, create a new set that is much smaller and easier to analyze statistically. The average, µ*, and the standard deviation, σ, are taken for each optimal sample set. µ* represents the impact of the specific parameter on the model, and σ represents the interaction between a parameter and the other input factors [31]. The greater the magnitude for these values, the greater the impact (µ*) and interaction (σ) for the corresponding input parameter. Since Morris is a relatively simple method that focuses more on the impact and less on the interactions between parameters, it was chosen as the global SA method for this work. It is important to note that, when performing SA on this model, the impact parameters were the input parameters. While each impact parameter has a corresponding output product (e.g., k1 corresponds to the hydrogen term for H2O output), there is not a necessary relationship between the degree of impact and the product itself. However, this work hypothesizes that there is at least a small direct correlation between the magnitude of impact for an impact parameter and the amount of corresponding product produced.

3. Results and Discussion

To ensure the effectiveness of this modeling, the mass closure for these runs had to be accounted for. Since this work was to understand, specifically, how the reactants affected each product through the sensitivity analysis, our entire process would collapse and be non-functional if mass closure was not accounted for. To do this, the mass flow controllers were calibrated beforehand, the flow in and out was measured regularly, and all liquids pulled from the system were weighed and accounted for. The difference between output and input displayed almost 96% closure for mass balance for all our work.
Using the data collected from the experimental runs, the mole functions were interpolated, the differential equations were solved, and the solution was then tested using Morris method SA. Initially, the DEs were loosely coupled, meaning that the solutions for MH2(t) and MCO(t) could be somewhat dependent on one another and the input parameters (Z1–Z9) could potentially interact. Each set of experimental data was analyzed nine times with Morris, on varying ranges from 0.1–50 to 0.1–500. The number of samples generated varied between 500 and 1000, and the optimal trajectories varied between 0 and 50, as described by Campolongo et al. [31]. The µ* and σ values were collected from the Morris output and the average for all nine analyses was taken. The µ* and σ averages for each data set were compared to the other data sets with the same catalyst–support combination and then compared to the values for the other category of data to discern any variance between the two types of supports.
Next, the two DEs were uncoupled. The CO2 term for the hydrogen equation, along with k10, was added to explore the possibility of WGS occurring in the 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst. Both equations were analyzed using Morris, following the same procedure as described above. The µ* and σ values were recorded for both the H2 and the CO equations. These were then compared, primarily the µ* values for k5 and k10, the CO2 impact parameters for the CO and H2, respectively. In theory, should µ* for k10 consistently be greater than that of k5, it could be indicative of hydrogen having a more substantial role in producing CO2. This could only happen through WGS, and thus if Z10 > Z5, it could indicate the possibility of WGS. Otherwise, if Z5 > Z10, this would most likely perpetuate the prevailing theory, primarily for cobalt catalyst FTS, that the Boudouard reaction is producing the CO2 [36,37].
Listed below are the data collected from the outlined modeling and sensitivity analysis procedure. As a key, Table 1 is the list of the five product categories and the corresponding impact parameters for the two ODEs in the model. For example, Z4 is the impact parameter for water in the carbon monoxide equation.
Recall that the µ* value for each parameter indicates how much that parameter influences the model. In other words, changing the impact parameter with the highest µ* value will alter the output of the model more than any of the other parameters. The µ* value does not indicate any level of interaction among the parameters, only how significantly that parameter affects the model. To examine the interactions between the parameters effectively, Sobol method SA should be used. Table 2 and Table 3 show the data collected from Morris method SA performed on the coupled DEs.
From these data, the order of impact for the nine impact parameters can be observed for both supports. For 20 wt% Co/SiO2, the order of impact is Z7 > Z3 > Z4 > Z1 > Z6 > Z2 > Z8 > Z9 > Z5. In contrast, the order for 0.5 wt% Pt, 25 wt% Co/Al2O3 is Z7 > Z3 > Z4 > Z1 > Z6 > Z2 > Z8 > Z5 > Z9. Although a subtle difference, the order of k9 and k5 is inverted, indicating that k5 had a greater impact on the model with the cobalt–alumina catalyst than the cobalt–silica. These orders of impact held constant across all three sets of data for each type of support, as can be seen in Table 2 and Table 3.
Table 4 lists the averages for each category side-by-side in order to compare the two supports more closely. From the hypothesis, this variation in the order of the parameter impact indicates that CO2, the corresponding product for Z5, is more impactful on the alumina support than the silica.
Furthermore, based on the proposed correlation between the products and corresponding parameters, alumina produces a greater amount of carbon dioxide than does silica. To verify this hypothesis, the original data from the data sets were examined, and the theory was confirmed. As Table 5 shows, the cobalt–alumina catalysts did in fact produce consistently more carbon dioxide than the cobalt–silica.
Due to this CO2 production difference, the two DEs were separated and examined individually using the Morris method. Only the Category 2 data were explored since the CO2 parameter was more impactful for the cobalt–alumina support. The additional carbon dioxide must come from an additional source. According to the hypothesis, that source could be either water–gas shift (WGS) or the Boudouard reaction. In order to account for any CO2 produced as a result of hydrogen interaction, thus allowing for the possibility of WGS, the parameter Z10 was added to the hydrogen equation along with the CO2 mole function. Table 6 and Table 7 show the data collected from Morris method SA for the separated equations.
The parameters Z5 and Z10 were of primary interest. As is evident from the data, the impact of k5 was invariably greater in magnitude for the CO equation than the impact of Z10 for the H2 equation across all three data sets in Category 2. Since Z5 > Z10, the likelihood of carbon dioxide being produced from WGS is very minimal. The results indicate that the extra CO2 is a product of the Boudouard reaction.

4. Materials and Methods

4.1. Catalyst Preparation

The cobalt and platinum starting materials were cobalt and platinum nitrate purchased from Alfa Aesar (Haverhill, MA, USA). A 25% cobalt alumina catalyst was prepared through a slurry impregnation procedure using Catalox 150 (high purity γ-alumina, ≈150 m2/g, Sasol, Hamburg, Germany) as a support. Next, the calcination of the catalysts was performed under air at 623 K for 4 h using a tube furnace (Lindberg/MPH, Riverside Road Riverside, MI 49084). The alumina was calcined for 10 h before impregnation and then cooled under an inert gas to room temperature. Co(NO3)2∙6H2O (99.9% purity) was used as the precursor for Co. In this method, which follows a Sasol patent [38], the ratio of the volume of solution used to the weight of alumina was 1:1, such that approximately 2.5 times the pore volume of solution was used to prepare the loading solution. A two-step impregnation method was used allowing a loading of two portions of 12.5% Co by weight. After the impregnation, the catalyst was then dried under vacuum by means of a rotary evaporator (Buchi R-100, Flawil, Switzerland) at 353 K, then slowly increased to 368 K. After the second impregnation/drying step, the catalyst was calcined under air flow at 623 K for 4 h.
The 20% Co/SiO2 catalyst was also prepared using the aqueous slurry phase impregnation method, with cobalt nitrate as the cobalt precursor. The support was PQ-SiO2 CS-2133, surface area about 352 m2/g. The catalyst was calcined in flowing air or flowing ~5% nitric oxide [39] in nitrogen at a rate of 1 L/min for 4 h at 623 K.

4.2. Characterization

4.2.1. BET Surface Area and Porosity Measurements

Surface area, pore volume, and average pore radius of the catalyst calcined at 623 K was evaluated through Brunauer-Emmett-Teller (BET) by using a Micromeritics Tri-Star 3000 gas adsorption analyzer system (Norcross, GA, USA) (Table 8). About 0.35 g of the catalyst was first weighed out, then added to a 3/8 in. o.d. sample tube. The adsorption gas was N2; sample analysis was performed at the boiling temperature of liquid nitrogen. Prior to the measurement, the chamber temperature was gradually increased to 433 K, then evacuated for several hours to approximately 6.7 Pa.

4.2.2. Hydrogen Chemisorption

Hydrogen chemisorption was conducted using temperature-programmed desorption (TPD), Table 8, with a Zeton-Altamira AMI-200 instrument (Pittsburgh, PA, USA). The Co/Al2O3 (≈0.22 g) was activated in a flow of 10 cm3/min of H2 mixed with 20 cm3/min of argon at 553 K for 10 h and then cooled under flowing H2 to 373 K. The sample was held at 373 K under flowing argon to remove and/or prevent adsorption of weakly bound species prior to increasing the temperature slowly to 623 K, the temperature at which oxidation of the catalyst was carried out.
The TPD spectrum was integrated and the number of moles of desorbed hydrogen determined by comparing its area to the areas of calibrated hydrogen pulses. The loop volume was first determined by establishing a calibration curve with syringe injections of hydrogen into a helium flow. Dispersion calculations were based on the assumption of a 1:1 H:CO stoichiometric ratio and a spherical cobalt cluster morphology. After TPD of hydrogen, the sample was reoxidized at 623 K using pulses of oxygen. The percentage of reduction was calculated by assuming that the metal reoxidized to Co3O4. Results of chemisorption are reported in Table 9.

4.2.3. Catalyst Activity Testing

Reaction experiments were conducted using a 1 L stirred-tank slurry reactor (STSR; Fort Worth, TX, USA) equipped with a magnetically driven stirrer with turbine impeller, a gas inlet line, and a vapor outlet line with a stainless steel (SS) fritted filter (2 mm) placed external to the reactor. A tube fitted with a SS fritted filter (0.5 mm opening) extending below the liquid level of the reactor was used to withdraw reactor wax (i.e., wax that is solid at room temperature), thereby maintaining a relatively constant liquid level in the reactor. Separate Brooks Instrument mass flow controllers (Hatfield, PA, USA) were used to control the flow rates of hydrogen and carbon monoxide. Carbon monoxide, prior to use, was passed through a vessel containing lead oxide on alumina to remove traces of iron carbonyls. The gases were premixed in an equalization vessel and fed to the STSR below the stirrer, which was operated at 750 rpm. The reactor temperature was maintained constant (≈274 K) using a temperature controller.
A 12 g amount of the oxide cobalt catalyst (sieved 63–106 µm) was loaded into a fixed-bed reactor for 12 h of ex situ reduction at 623 K and atmospheric pressure using a gas mixture of H2/He (60 L/h) with a molar ratio of 1:3. The reduced catalyst was transferred to an already capped 1 L STSR containing 310 g of melted Polywax 3000 (Baker Hughes, Houston, TX) under the protection of inert nitrogen gas. The reactor used for the reduction was weighed three separate times, before and after reduction, and after catalyst transfer, in order to obtain an accurate amount of reduced catalyst that was added. The transferred catalyst was again reduced in situ at 503 K at atmospheric pressure using pure hydrogen (20 L/h) for another 24 h before starting the FT reaction. Each run was held at 493 K and 18.7 atm (275 psig). The CO:H2 ratio was 1:2 for each experiment. The weight hourly space velocity (WHSV) varied from 3.0 to 6.0 slph/g catalyst, depending on the purpose of the experiment.
Gas, water, oil, light wax, and heavy wax samples were collected daily and analyzed. Heavy wax samples were collected in a 473 K hot trap connected to the filter-containing dip tube. The vapor phase in the region above the reactor slurry passed continuously to the warm (373 K) and then the cold (273 K) traps located external to the reactor. The light wax and water mixture were collected daily from the warm trap and an oil plus water sample from the cold trap. Tail gas from the cold trap was analyzed with an online HP Quad Series Micro gas chromatograph (GC) (Palo Alto, CA, USA), providing molar compositions of C1–C5 olefins and paraffins, as well as H2, CO and CO2. The analysis of the aqueous phase used an SRI 8610C GC (20720 Earl St. Torrance, CA 90503, USA) with a thermal conductivity detector (TCD).
Products were analyzed with an Agilent 6890 GC (Santa Clara, CA, USA) with a flame ionization detector (FID) and a 60 m DB-5 column. Hydrogen and carbon monoxide conversions were calculated based on GC analysis of the gas products, the gas feed rate, and the gas flow that was measured at the outlet of the reactor.

5. Conclusions

In conclusion, the mathematical model created in this work accurately described the FTS process in terms of rate of change of syngas. When the model was analyzed using Morris method SA, a significant difference was observed between cobalt–silica and cobalt–platinum–alumina catalysts, namely a difference in the effect of the impact parameters k5 and k9. Based on the initial hypothesis, this difference indicated that the cobalt–platinum–alumina generated more carbon dioxide than did the cobalt–silica. This was verified by the raw data as well. These conclusions are highly specific to the two catalyst–support types in question, and thus more general conclusions must come from further research. Thus, the equations were separated, and the cobalt–platinum–alumina catalyst data were reexamined, now with an additional CO2 term on the hydrogen equation to test for the possibility of WGS. The observed results showed that the carbon monoxide contribution to CO2 was consistently higher than the hydrogen contribution. Since WGS is the only pathway in which hydrogen produces CO2 during FTS, WGS was excluded as a source of CO2 production. Hence, the Boudouard reaction was determined to be the more probable source of additional carbon dioxide. Though the Boudouard reaction yields a high amount of carbon on the surface, to the best of our knowledge, there is no direct evidence that this reaction leads directly and solely to inert graphitic carbon (coke). Moreover, one of the main mechanistic routes proposed for FT is the carbide mechanism, where CO completely dissociates before reacting with hydrogen. In order to understand this rate of formation for C, evidence of more than one type of carbon on the catalyst surface, an active phase and an unreactive phase, has been observed on multiple catalysts such as cobalt [40] and Fe [41].
Further indications using sensitivity analysis show that both the H2 and CO equations displaying Z3 and Z7 obtained the highest value, indicating the liquids are most sensitive to any deviations brought on by H2 and CO. Though expected, these details could further indicate the importance for a specific balance in the FTS system between the two reactants, suggesting the difficulty in tuning a specific catalyst for C5+ products as a whole. Additionally, the C5+ products are more sensitive in the CO reaction than the H2, which implies that the balance for the FTS system lies more heavily in the CO than for H2, primarily describing the importance of CO, possibly due to the metal’s ability to back donation [42] in the FTS process.
While this research is preliminary, the model and analysis produced some excellent descriptive results for FTS. Although this work only examined cobalt catalyst data, future work could include the incorporation of multiple catalysts and supports. Impact parameters could potentially be optimized in future work to find the ideal values for a given catalyst. The model can be adapted to include terms for olefins and paraffins, and the liquid (C5+) term could be broken down into individual terms for each number of carbons per chain (e.g., C5, C6, C7, etc.). The equations could be changed to predict specific product amounts based on initial conditions such as catalyst, support, temperature, pressure, flowrate, etc. Ultimately, this type of FTS modeling not only describes the process in various and versatile ways but can also be adapted and altered to potentially produce a predictive model.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4344/9/10/857/s1.

Author Contributions

H.W. gathered information and wrote most of the paper. W.D.S. provided most of the analytical results and oversaw part of the process, G.J. provided the catalyst characterization and finalized edits, M.K.G. provided experimental oversight, data collection/analysis and editing. D.C. designed and coded the algorithm and oversaw the entire project.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fischer, F.; Tropsch, H. Synthesis of Petroleum at Atmospheric Pressure from Gasification Products of Coal. Brennstoff-Chemie 1926, 7, 97–104. [Google Scholar]
  2. Fischer, F.; Tropsch, H. Development of the Benzene Synthesis from Carbon Monoxide and Hydrogen at Atmospheric Pressure. Brennstoff-Chemie 1930, 11, 489–500. [Google Scholar]
  3. Schulz, H. Short history and present trends of Fischer–Tropsch synthesis. Appl. Catal. 1999, 186, 3–12. [Google Scholar] [CrossRef]
  4. Davis, B.H. Fischer–Tropsch Synthesis: Reaction mechanisms for iron catalysts. Catal. Today 2009, 141, 25–33. [Google Scholar] [CrossRef]
  5. Navarro, V.; van Spronsen, M.A.; Frenken, J.W.M. In situ observation of self-assembled hydrocarbon Fischer–Tropsch products on a cobalt catalyst. Nat. Chem. 2016, 8, 929–934. [Google Scholar] [CrossRef]
  6. Chakrabarti, D.; Gnanamani, M.K.; Shafer, W.D.; Ribeiro, M.C.; Sparks, D.E.; Prasad, V.; de Klerk, A.; Davis, B.H. Fischer–Tropsch Mechanism: 13C18O Tracer Studies on a Ceria–Silica Supported Cobalt Catalyst and a Doubly Promoted Iron Catalyst. Ind. Eng. Chem. Res. 2015, 54, 6438–6453. [Google Scholar] [CrossRef]
  7. Tuxen, A.; Carenco, S.; Chintapalli, M.; Chuang, C.; Escudero, C.; Pach, E.; Jiang, P.; Borondics, F.; Beberwyck, B.; Alivisatos, A.P.; et al. Size-Dependent Dissociation of Carbon Monoxide on Cobalt Nanoparticles. J. Am. Chem. Soc. 2013, 135, 2273–2278. [Google Scholar] [CrossRef]
  8. Qi, Y.; Yang, J.; Chen, D.; Holmen, A. Recent Progresses in Understanding of Co-Based Fischer–Tropsch Catalysis by Means of Transient Kinetic Studies and Theoretical Analysis. Catal. Lett. 2015, 145, 145–161. [Google Scholar] [CrossRef]
  9. Van Santen, R.A.; Markvoort, A.J.; Filot, I.A.W.; Ghouriab, M.M.; Hensen, E.J.M. Mechanism and Microkinetics of the Fischer–Tropsch reaction. Phys. Chem. Chem. Phys. 2013, 15, 17038–17063. [Google Scholar] [CrossRef]
  10. Banerjee, A.; van Bavel, A.V.; Kuipers, H.P.C.E.; Saeys, M. CO Activation on Realistic Cobalt Surfaces: Kinetic Role of Hydrogen. ACS Catal. 2017, 7, 5289–5293. [Google Scholar] [CrossRef]
  11. Todic, B.; Bhatelia, T.; Froment, G.F.; Ma, W.; Jacobs, G.; Davis, B.H.; Bukur, D.B. Kinetic Model of Fischer-Tropsch Synthesis in a Slurry Reactor on Co-Re/Al2O3 Catalyst. Ind. Eng. Chem. Res. 2013, 52, 669–679. [Google Scholar] [CrossRef]
  12. Keyvanloo, K.; Fisher, M.J.; Hecker, W.C.; Lancee, R.J.; Jacobs, G.; Bartholomew, C.H. Kinetics of deactivation by carbon of a cobalt Fischer–Tropsch catalyst: Effects of CO and H2 partial pressures. J. Catal. 2015, 327, 33–47. [Google Scholar] [CrossRef]
  13. Ma, W.; Jacobs, G.; Sparks, D.E.; Spicer, R.L.; Davis, B.H.; Klettlinger, J.L.S.; Yen, C.H. Fischer-Tropsch synthesis: Kinetics and water effect study over 25%Co/Al2O3 catalysts. Catal. Today 2014, 228, 158–166. [Google Scholar] [CrossRef]
  14. Botes, F.G.; van Dyk, B.; McGregor, C. The Development of a Macro Kinetic Model for a Commercial Co/Pt/Al2O3 Fischer-Tropsch Catalyst. Ind. Eng. Chem. Res. 2009, 48, 10439–10447. [Google Scholar] [CrossRef]
  15. Brady, R.C.; Pettit, R. Mechanism of the Fischer-Tropsch reaction. The chain propagation step. J. Am. Chem. Soc. 1981, 103, 1287–1289. [Google Scholar] [CrossRef]
  16. Yang, J.; Shafer, W.D.; Pendyala, V.R.R.; Jacobs, G.; Chen, D.; Holmen, A.; Davis, B.H. Fischer–Tropsch Synthesis: Using Deuterium as a Tool to Investigate Primary Product Distribution. Catal. Lett. 2014, 144, 524–530. [Google Scholar] [CrossRef]
  17. Shafer, W.D.; Jacobs, G.; Davis, B.H. Fischer–Tropsch Synthesis: Investigation of the Partitioning of Dissociated H2 and D2 on Activated Cobalt Catalyst. ACS Catal. 2012, 2, 1452–1456. [Google Scholar] [CrossRef]
  18. Yang, J.; Qia, Y.; Zhu, J.; Zhu, Y.-A.; Chen, D.; Holmen, A. Reaction mechanism of CO activation and methane formation on Co Fischer–Tropsch catalyst: A combined DFT, transient, and steady-state kinetic modeling. J. Catal. 2013, 308, 37–49. [Google Scholar] [CrossRef]
  19. Ojeda, M.; Nabar, R.; Nilekar, A.U.; Ishikawa, A.; Mavrikakis, M.; Iglesia, E. CO activation pathways and the mechanism of Fischer–Tropsch synthesis. J. Catal. 2010, 272, 287–297. [Google Scholar] [CrossRef]
  20. Ojeda, M.; Li, A.; Nabar, R.; Nilekar, A.U.; Mavrikakis, M.; Iglesia, E. Kinetically Relevant Steps and H2/D2 Isotope Effects in Fischer−Tropsch Synthesis on Fe and Co Catalysts. J. Phys. Chem. C 2010, 114, 19761–19770. [Google Scholar] [CrossRef]
  21. Gracia, J.M.; Prinsloo, F.F.; Niemantsverdriet, J.W. Mars-van Krevelen-like Mechanism of CO Hydrogenation on an Iron Carbide Surface. Catal. Lett. 2009, 133, 257–261. [Google Scholar] [CrossRef] [Green Version]
  22. Davis, B.H. Fischer-Tropsch synthesis: Current mechanism and futuristic needs. Fuel Process. Technol. 2001, 71, 157–166. [Google Scholar] [CrossRef]
  23. Filip, L.; Zámostný, P.; Rauch, R. Mathematical model of Fischer-Tropsch synthesis using variable alpha-parameter to predict product distribution. Fuel 2019, 243, 603–609. [Google Scholar] [CrossRef]
  24. Ahón, V.R.; Costa, E.F.; Monteagudo, J.E.P.; Fontes, C.E.; Biscaia, E.C.; Lage, P.L.C. A comprehensive mathematical model for the Fischer–Tropsch synthesis in well-mixed slurry reactors. Chem. Eng. Sci. 2005, 60, 677–694. [Google Scholar] [CrossRef]
  25. Yu, F.; Lin, T.; Wang, X.; Li, S.; Lu, Y.; Wang, H.; Zhong, L.; Sun, Y. Highly selective production of olefins from syngas with modified ASF distribution model. Appl. Catal. A Gen. 2018, 563, 146–153. [Google Scholar] [CrossRef]
  26. Nikbakht, N.; Mirzaei, A.A.; Atashi, H. Kinetic modeling of the Fischer-Tropsch reaction over a zeolite supported Fe-Co-Ce catalyst prepared using impregnation procedure. Fuel 2018, 229, 209–216. [Google Scholar] [CrossRef]
  27. Ostadi, M.; Rytter, E.; Hillestad, M. Evaluation of kinetic models for Fischer–Tropsch cobalt catalysts in a plug flow reactor. Chem. Eng. Res. Des. Trans. Inst Chem. Eng. Part A 2016, 114, 236–246. [Google Scholar] [CrossRef]
  28. Nabipoor Hassankiadeh, M.; Haghtalab, A. Product Distribution of Fischer-Tropsch Synthesis in a Slurry Bubble Column Reactor Based on Langmuir-Freundlich Isotherm. Chem. Eng. Commun. 2013, 200, 1170–1186. [Google Scholar] [CrossRef]
  29. Seborg, D.E.; Mellichamp, D.A.; Thomas, F. Process Dynamics and Control; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2011. [Google Scholar]
  30. Papraćanin, E. Local and Global Sensitivity Analysis of Modeling Parameters for Composting Process. Technol. Acta 2019, 11, 9–16. [Google Scholar]
  31. Campolongo, F.; Cariboni, J.; Saltelli, A. An effective screening design for large models. Environ. Model. Softw. 2007, 22, 1509–1518. [Google Scholar] [CrossRef]
  32. Zhang, X.; Trame, M.; Lesko, L.; Schmidt, S. Sobol Sensitivity Analysis: A Tool to Guide the Development and Evaluation of Systems Pharmacology Models. CPT Pharmacomet. Syst. Pharmacol. 2015, 4, 69–79. [Google Scholar] [CrossRef]
  33. Sobol, I. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math. Comput. Simul. 2001, 55, 271–280. [Google Scholar] [CrossRef]
  34. Cariboni, J.; Gatelli, D.; Liska, R.; Saltelli, A. The role of sensitivity analysis in ecological modeling. Ecol. Model. 2007, 203, 167–182. [Google Scholar] [CrossRef]
  35. Nossent, J.; Elsen, P.; Bauwens, W. Sobol’ sensitivity analysis of a complex environmental model. Environ. Model. Softw. 2011, 26, 1515–1525. [Google Scholar] [CrossRef]
  36. Hunt, J.; Ferrari, A.; Lita, A.; Crosswhite, M.; Ashley, B.; Stiegman, A.E. Microwave-Specific Enhancement of the Carbon-Carbon Dioxide (Boudouard) Reaction. J. Phys. Chem. 2013, 117, 26871–26880. [Google Scholar] [CrossRef]
  37. Bremmer, G.M.; Zacharaki, E.; Sjȁstad, A.O.; Navarro, V.; Frenken, J.W.M.; Kooyman, P.J. In situ TEM observation of the Boudouard Reaction: Multi-layered graphene formation from CO on cobalt nanoparticles at atmospheric pressure. Faraday Discuss. 2017, 197, 337–351. [Google Scholar] [CrossRef]
  38. Espinoza, R.L.; Visagie, J.L.; van Berge, P.J.; Bolder, F.H. Catalysts. US Patent 5,733,839, 31 March 1998. [Google Scholar]
  39. Jacobs, G.; Ma, W.; Davis, B.H.; Cronauer, D.C.; Kropf, J.; Marshall, C.L. Fischer–Tropsch Synthesis: TPR-XAFS Analysis of Co/Silica and Co/Alumina Catalysts Comparing a Novel NO Calcination Method with Conventional Air Calcination. Catal. Lett. 2010, 140, 106–115. [Google Scholar] [CrossRef]
  40. Chen, W.; Kimpel, T.F.; Song, Y.; Chiang, K.-F.; Zijlstra, B.; Pestman, R.; Wang, P.; Hensen, E.J.M. Influence of Carbon Deposits on the Cobalt-Catalyzed Fischer−Tropsch Reaction: Evidence of a Two-Site Reaction Model. ACS Catal. 2018, 8, 1580–1590. [Google Scholar] [CrossRef]
  41. Zhang, C.; Zhao, G.; Liu, K.; Yang, Y.; Xiang, H.; Li, Y. Adsorption and reaction of CO and hydrogen on iron-based Fischer–Tropsch synthesis catalysts. J. Mol. Catal. A 2010, 328, 35–43. [Google Scholar] [CrossRef]
  42. Shafer, W.D.; Gnanamani, M.K.; Graham, U.M.; Yang, J.; Masuku, C.M.; Jacobs, G.; Davis, B.H. Fischer–Tropsch: Product Selectivity–The Fingerprint of Synthetic Fuels. Catalysts 2019, 9, 259. [Google Scholar] [CrossRef]
Table 1. Key for understanding the impact parameters and their corresponding equation (columns 2 and 3) and product (column 1).
Table 1. Key for understanding the impact parameters and their corresponding equation (columns 2 and 3) and product (column 1).
ProductCOH2
H2OZ4Z1
CH4Z6Z2
CO2Z5Z10*
C2–C4Z8Z9
C5+Z7Z3
The * is to note that Z10 was added during the second set of sensitivity analysis runs with the split equations. It does not appear in the first set of data.
Table 2. Averages of the µ* values across the nine sensitivity analysis (SA) runs of each data set in Category 1, the data from experiments using 20 wt% Co/SiO2 catalyst.
Table 2. Averages of the µ* values across the nine sensitivity analysis (SA) runs of each data set in Category 1, the data from experiments using 20 wt% Co/SiO2 catalyst.
ParametersData Set 1AData Set 1BData Set 1CAverage
Z11633.961866.032117.571872.52
Z2137.167125.673135.359132.733
Z321,848.627,452.032,795.527,365.4
Z43216.693743.594242.443734.24
Z525.488421.085144.291030.2882
Z6276.976253.349274.167268.164
Z745,393.156,067.567,022.856,161.1
Z877.275082.030195.384084.8964
Z938.478140.221047.240841.9800
The last column is an average of µ* across all three runs for each impact parameter.
Table 3. Averages of the µ* values across the nine SA runs of each data set in Category 2, the data from experiments using 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
Table 3. Averages of the µ* values across the nine SA runs of each data set in Category 2, the data from experiments using 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
ParametersData Set 2AData Set 2BData Set 2CAverage
Z12628.022147.252356.762377.34
Z2196.136166.964179.690180.930
Z339,224.130,725.434,196.234,715.2
Z45269.894190.984644.674701.85
Z594.454864.402165.795674.8841
Z6401.960338.096364.823368.293
Z780,809.962,518.770,123.671,150.8
Z8129.969315.635129.089191.564
Z962.252257.509063.468761.0766
The last column is an average of µ* across all three runs for each impact parameter.
Table 4. Comparison of the average µ* values between the two supports.
Table 4. Comparison of the average µ* values between the two supports.
ParametersData Category 1Data Category 2
Z11872.522377.34
Z2132.733180.930
Z327,365.434,715.2
Z43734.244701.85
Z530.288274.8841
Z6268.164368.293
Z756,161.171,150.8
Z884.8964191.564
Z941.980061.0766
Category 1 is 20 wt% Co/SiO2, Category 2 is 0.5 wt% Pt, 25 wt% Co/Al2O3.
Table 5. Average moles of CO2 produced per sample period for each data set.
Table 5. Average moles of CO2 produced per sample period for each data set.
Data SetCategory 1Category 2
A0.03171080.0791931
B0.03159820.0824212
C0.06428150.0675372
Average0.04253020.0763838
Recall that Category 1 contains data sets with 20 wt% Co/SiO2 and Category 2 contains data sets with 0.5 wt% Pt, 25 wt% Co/Al2O3. The bold values are the averages for each category.
Table 6. Average across the nine SA runs for the µ* values of each data set with the 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
Table 6. Average across the nine SA runs for the µ* values of each data set with the 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
H2 ParametersData Set 2AData Set 2BData Set 2CAverage
Z12415.851932.342146.772164.99
Z2181.013154.389165.666167.023
Z335,685.728,044.2231,052.331,594.1
Z957.391953.270458.446356.3696
Z1041.227629.266029.332633.2754
Hydrogen equation only. The last column is an average of µ* across all three runs for each parameter. The row in bold is the parameter of interest (i.e., the parameter for CO2).
Table 7. Average across the nine SA runs for the µ* values of each data set with the 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
Table 7. Average across the nine SA runs for the µ* values of each data set with the 0.5 wt% Pt, 25 wt% Co/Al2O3 catalyst.
CO ParametersData Set 2AData Set 2BData Set 2CAverage
Z44867.073805.304251.974308.11
Z586.307159.412160.416068.7117
Z6371.787310.962339.218340.656
Z774,015.757,129.364,064.865,069.9
Z8118.735107.272119.497115.168
Carbon monoxide equation only. The last column is an average of µ* across all three runs for each parameter. The row in bold is the parameter of interest (i.e., the parameter for CO2).
Table 8. Surface area and pore size distribution of silica-supported cobalt catalysts.
Table 8. Surface area and pore size distribution of silica-supported cobalt catalysts.
Support/CatalystsCo, wt%Reduction Temperature (K)SBET, m2/gTotal Pore Volume, cm3/gAverage Pore Radius (nm)
20% Co/SiO2206232480.8565.6
0.5%Pt-25%Co/Al2O32562398.50.2184.4
Table 9. H2 chemisorption (temperature-programmed desorption (TPD)) and pulse reoxidation of metallic phases of cobalt-supported catalysts (cat.).
Table 9. H2 chemisorption (temperature-programmed desorption (TPD)) and pulse reoxidation of metallic phases of cobalt-supported catalysts (cat.).
Support/Catalystsµmol H2 Desorbed Per g Cat.Uncorrected% DispersionUncorrected Diameter (nm)µmol O2 Consumed Per g Cat.% ReductionCorrected% DispersionCorrected Diameter (nm)
20% Co/SiO241.22.4342.51086485.0620.4
0.5%Pt-25%Co/Al2O3132.96.216.5173361.310.310

Share and Cite

MDPI and ACS Style

Williams, H.; Gnanamani, M.K.; Jacobs, G.; Shafer, W.D.; Coulliette, D. Fischer–Tropsch Synthesis: Computational Sensitivity Modeling for Series of Cobalt Catalysts. Catalysts 2019, 9, 857. https://doi.org/10.3390/catal9100857

AMA Style

Williams H, Gnanamani MK, Jacobs G, Shafer WD, Coulliette D. Fischer–Tropsch Synthesis: Computational Sensitivity Modeling for Series of Cobalt Catalysts. Catalysts. 2019; 9(10):857. https://doi.org/10.3390/catal9100857

Chicago/Turabian Style

Williams, Harrison, Muthu K. Gnanamani, Gary Jacobs, Wilson D. Shafer, and David Coulliette. 2019. "Fischer–Tropsch Synthesis: Computational Sensitivity Modeling for Series of Cobalt Catalysts" Catalysts 9, no. 10: 857. https://doi.org/10.3390/catal9100857

APA Style

Williams, H., Gnanamani, M. K., Jacobs, G., Shafer, W. D., & Coulliette, D. (2019). Fischer–Tropsch Synthesis: Computational Sensitivity Modeling for Series of Cobalt Catalysts. Catalysts, 9(10), 857. https://doi.org/10.3390/catal9100857

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop