Selectivity of the Fischer – Tropsch process : deviations from single alpha product distribution explained by gradients in process conditions

The Fischer–Tropsch (FT) process, in which syngas is converted to hydrocarbons, is a key conversion step in one of the most important routes to alternative fuels and chemicals. Various models exist for the selectivity of this reaction. The simplest and most widely used is the Anderson–Schulz–Flory (ASF) model, where the FT reaction is modeled as an addition polymerization reaction with chain growth probability a. The resulting product distribution, depicted as a plot of the logarithmic molar fraction versus the carbon number (‘ASF-plot’), is a straight line. However, deviations from this ASF product distribution are regularly observed. Commonly reported deviations from the ASF model are the C1 (methane) and C2 (ethane/ethene) values. Typically, C1 is formed in a greater quantity than that predicted by the ASF distribution, while a much lower amount of C2 is usually measured. Another deviation from the ASF product distribution shape that is encountered in practice is a chainlength-dependent chain growth probability. Instead of a linear trend in an ‘ASF-plot’ (cf. Fig. 1c), the product distribution shows a double chain growth probability or a curved distribution. The fractions of heavier hydrocarbons are in some cases lower than expected based on the ASF model but in most cases they are higher. The cause of these positive deviations from ASF product distributions is a disputed topic. It has been attributed to flaws in experimentation, the existence of multiple chain growth mechanisms or multiple different active sites on the catalyst. However, more recent work assumes that olefin readsorption and incorporation is responsible. In this communication we propose an alternative explanation for non-ASF product distributions which is nevertheless based on a single-chain-growth-probability model. Our explanation is based on the process condition dependence of a, where we assume a variable-a model that takes into account the local reactor temperature and the syngas ratio. It has previously been demonstrated that gradients in the process conditions may lead to selectivity gradients. Temperature and concentration gradients are the most likely variables that cause variations in selectivity. These gradients can take place on the particle scale due to diffusion limitations inside the catalyst particle, as well as on the reactor scale, e.g. due to heat transport limitations in a packed bed reactor. We will show that these gradients in process variables can result in an apparent non-ASF Fischer– Tropsch product distribution. Gradients in FT process conditions and resulting chain growth probability and product distributions were modeled on the particle scale and the reactor scale. For both cases, a process-condition-dependent a-model proposed by Vervloet et al. was used, which captured the experimental trends well:

The Fischer-Tropsch (FT) process, in which syngas is converted to hydrocarbons, is a key conversion step in one of the most important routes to alternative fuels and chemicals.Various models exist for the selectivity of this reaction.The simplest and most widely used is the Anderson-Schulz-Flory (ASF) model, 1 where the FT reaction is modeled as an addition polymerization reaction with chain growth probability a.The resulting product distribution, depicted as a plot of the logarithmic molar fraction versus the carbon number ('ASF-plot'), is a straight line.
However, deviations from this ASF product distribution are regularly observed.Commonly reported deviations from the ASF model are the C1 (methane) and C2 (ethane/ethene) values.Typically, C1 is formed in a greater quantity than that predicted by the ASF distribution, while a much lower amount of C2 is usually measured. 2Another deviation from the ASF product distribution shape that is encountered in practice is a chainlength-dependent chain growth probability.Instead of a linear trend in an 'ASF-plot' (cf.Fig. 1c), the product distribution shows a double chain growth probability or a curved distribution. 3,4[8][9][10][11][12][13] The cause of these positive deviations from ASF product distributions is a disputed topic.It has been attributed to flaws in experimentation, [8][9][10] the existence of multiple chain growth mechanisms 11 or multiple different active sites on the catalyst. 4,12wever, more recent work assumes that olefin readsorption and incorporation is responsible. 2,6,7,13n this communication we propose an alternative explanation for non-ASF product distributions which is nevertheless based on a single-chain-growth-probability model.Our explanation is based on the process condition dependence of a, where we assume a variable-a model that takes into account the local reactor temperature and the syngas ratio. 14It has previously been demonstrated that gradients in the process conditions may lead to selectivity gradients. 3,8,14Temperature and concentration gradients are the most likely variables that cause variations in selectivity.These gradients can take place on the particle scale due to diffusion limitations inside the catalyst particle, as well as on the reactor scale, e.g.due to heat transport limitations in a packed bed reactor. 8,15We will show that these gradients in process variables can result in an apparent non-ASF Fischer-Tropsch product distribution.
Gradients in FT process conditions and resulting chain growth probability and product distributions were modeled on the particle scale and the reactor scale.For both cases, a process-condition-dependent a-model proposed by Vervloet et al. 14 was used, which captured the experimental trends well: This model has three parameters: k a is the ratio of rate constants for the termination over propagation reactions, b is the syngas ratio power constant, and DE a is the difference in activation energies for the propagation and termination reactions.
The model is based on the assumptions that a is independent of chain length, and is derived from the ratio of propagation and termination reaction rates, which follow an Arrhenius-type temperature dependence. 14More detailed approaches, based on micro-kinetics, 6,16 are also proposed in the literature.These complex models have the drawback that a large number of parameters are unknown and must be estimated theoretically or from experimental data.In this communication, we have

View Article Online
View Journal | View Issue stuck to the simpler a-based model to simply provide insight into the FT process and predict the overall selectivity using only these three parameters.
To demonstrate the effect of reactant concentration profiles at the particle scale, a steady state reaction-diffusion model was applied for a spherical catalyst particle.This model, described in detail by Vervloet et al., 14 is based on dimensionless steady state reaction-diffusion mass balances in a catalyst particle.The relevant process variables considered in the particle model are the particle diameter d cat , the syngas ratio H 2 /CO, the temperature T and the catalyst activity multiplication factor F YS , relative to the Yates & Satterfield kinetics. 18o describe the behavior at the reactor scale, a simple 1-D reactor model was set up for a single tube.Since the FT process is strongly exothermal, often axial temperature gradients will be found.Such gradients can also lead to a deviation from a single a product distribution.We illustrated this by imposing a simple linear temperature profile which can also be encountered in industrial practice. 17In this steady state reactor model, mass transfer limitations and radial gradients were assumed to be absent and the syngas ratio was kept constant at 2. The gas phase, which was modeled to only include H 2 and CO, was therefore in instant equilibrium with the liquid phase according to Henry's law.The linear temperature increase was from T = 210 to 240 1C.Other process conditions and assumptions can be found in Appendix A.
At both the particle and the reactor scale, the reaction rate was assumed to follow Langmuir-Hinshelwood kinetics, proposed by Yates and Satterfield, 18 where the catalyst activity multiplication factor F YS accounts for the increase in catalyst activity since the publication of the original values.The reaction rate and adsorption constant reported by Maretto and Krishna 19 are used.Finally, the overall Fischer-Tropsch product distribution was calculated in two different ways: firstly, using a single average a, based on the reaction rate weighted integrated a over the catalyst or reactor volume; and secondly, by summating local (ASF) product distributions weighed with the reaction rate and incremental volume.
At the particle scale, we demonstrated the effect of diffusion limitations on the chain growth probability, product distribution and C 5+ selectivity of the FT reaction.A system with a moderate diffusion length, representative of packed bed reactors, was simulated.This example is illustrated in Fig. 1, where the reactant concentrations, reaction rate, a gradients and the resulting product distributions are shown.In this diffusionlimited catalyst particle, the calculation of the product distribution from the average a gave a straight distribution with a C 5+ selectivity of 54%, while the weighted summation of local distributions gave a curved plot with a positive deviation in the higher carbon numbers and a C 5+ selectivity of 79%.There are thus large deviations in the shape of the product distribution and in the C 5+ selectivity between these two calculation methods.
The cause of the curved product distribution is the occurrence of gradients in the local syngas ratio in the catalyst particle due to diffusion limitations.Towards the centre of the particle, CO becomes depleted, causing an increase in the H 2 /CO ratio, resulting in a decrease in the local a.Correspondingly, the local productivities of the components are different, resulting in a curved total product distribution when summated.In contrast to the chain growth probability, a maximum reaction rate can occur inside the particle at lower CO concentrations due to the kinetics.
The prominence of this CO depletion depends strongly on the diffusion length, temperature, syngas ratio and catalyst activity.At H 2 /CO = 2, a common choice for experimental conditions since it is close to the overall stoichiometric Fig. 1 Results of the reaction-diffusion model for the spherical catalyst particle at T = 488 K, the gas phase molar ratio H 2 /CO in = 1.8, d cat = 0.5 mm, F YS = 2, P = 30 bar. z = 0 at the centre and z = 1 at the surface of the particle.(A) Liquid H 2 and CO concentration profile in the catalyst particle normalized to CO concentration at the external particle surface, (B) profile of chain growth probability and CO conversion rate, (C) resulting product distributions calculated from local product distributions and from average a.
Table 1 Overview of reaction-diffusion model results.Significant difference (>30% in S C 5+ ) between product distribution calculated by average a and summation of local distributions indicated by ''!''.No significant difference (o1%) indicated by ''=''.Intermediate cases were not observed under these conditions.In all cases, P = 30 bar consumption of the reactants, even catalyst particles as small as 200 mm can become depleted in CO at enhanced catalyst activity.This potentially leads to non-linear ASF product distributions when the concentration dependency of a is accounted for.A high temperature and/or catalyst activity amplifies this effect, because CO is depleted more rapidly at higher reaction rates.Different conditions and diffusion lengths were simulated to represent various reactor systems (i.e.short diffusion lengths to represent the slurry reactor and thin coated reactor internals, and moderate and large diffusion lengths to represent packed bed reactors).We modeled a catalyst with a conservative activity (F YS = 2, representing a factor 2 improvement in catalyst activity from the original values 18 ), and a high catalyst activity (F YS = 10, representing the best state-of-the-art catalysts).The product distribution was calculated in two different ways: using the average chain growth probability and by summating local product distributions.In strongly diffusion-limited systems a significant difference between these two distributions was found (Table 1).
To compare this modeling with experimental results, we used the reaction-diffusion model to simulate the process conditions used by Visconti et al. 6 (p = 20 bar, T = 230 1C, H 2 /CO = 2.1, d cat = 100 mm).For this purpose, the catalyst dependent parameters in the variable-a model (k a , b and DE a , see eqn (1)) were tuned to fit the data.The resulting product distribution is compared to the experimental data in Fig. 2. The modeled product distribution captures the experimental data well.It is especially noteworthy that the CO conversion level is reproduced, as well as the curved product distribution, including the value for C 1 .This result implies that diffusion limitations in the particle could cause a non-ASF product distribution such as that found experimentally.
At the reactor scale, the effect of temperature gradients on the product distribution is demonstrated by imposing a linear temperature increase of 30 K over the length of the reactor (see Fig. 3).Generally, and captured in eqn (1), a decreases with increasing temperature along the reactor.The effect of this a-gradient is a non-linear ASF product distribution, similar to the shape found in diffusion-limited catalyst particles.
It was found that the effect of the imposed temperature profile on the apparent chain growth, and thus on the product distribution, is smaller than the effect of the reactant ratio at the catalyst particle scale.This is due to the parameters in the variable-a model.Although the used reactor model is a gross simplification of a real FT reactor, this example clearly shows that typical axial temperature gradients encountered in FT reactors can cause a non-ASF product distribution too.

Conclusions
Non-ASF product distributions in Fischer-Tropsch synthesis can occur when there is a gradient in process variables (temperature, concentration) leading to a gradient in the chain growth probability a.The origin of this effect can be found on the particle scale, i.e. concentration gradients caused by diffusion limitations in the catalyst, or on the reactor scale, i.e. gradients in local conditions such as temperature.When such gradients are modeled, using the average chain growth probability to predict the product composition can lead to erroneous results.The correct product distribution should be obtained by volume based summation of local product distributions, weighted with the local reaction rate.

Table 2
Values of parameters used in the 1-D reactor model