OPTIMUM REACTOR DESIGN IN METHANATION PROCESSES WITH NONUNIFORM CATALYSTS

A generic methodology is developed to design a heterogeneous catalytic reactor for methanation processes. For the optimization of a heterogeneous catalytic reactor, nonuniform catalyst pellets such as a layered catalyst are considered with respect to reaction type, reactor performance, and component distribution inside the catalyst. Heterogeneous uniform and nonuniform catalyst models were developed to analyze the effect of mass and heat transfer between both bulk phase and catalyst surface and inside a catalyst pellet. Then, concentration profiles of hydrogen and carbon monoxide in the catalyst pellet and along the reactor axis were obtained by analyzing simulation results. It was shown that the application of different types of nonuniform catalyst pellets at a certain number of separate zones within a reactor could produce higher catalyst performance than a reactor with uniform catalyst. Furthermore, it proved a significant decrease of catalyst deactivation behavior such as coking and sintering. Layered catalysts were optimized to maximize an overall reactor performance over the catalyst lifetime, achieving capital cost reduction characterized by reactor size, catalyst amount, and degree of catalyst deactivation. Last, temperature control throughout the reactor operating periods was strategically planned for a reactor operation with distribution of nonuniform catalyst pellets. This methodology can also be usefully applied to the design of heterogeneous catalytic reactors for other processes such as hydro-treating process and cracking process.     

Optimization of Fischer‐Tropsch Process in a Fixed‐Bed Reactor Using Non‐uniform Catalysts

Optimization of Fischer‐Tropsch (FT) process in a fixed‐bed reactor is carried out using non‐uniform catalysts. The C₅⁺ yield of the reactions is maximized utilizing a combination of non‐uniform catalysts across the bed. A 1D heterogeneous model is developed to simulate the bed containing uniform and non‐uniform catalysts. It is found that the egg‐shell and surface‐layered catalysts result in higher C₅⁺ yield. Moreover, effects of cooling temperature are studied. Genetic Algorithm (GA) and Successive Quadratic Programming (SQP) methods are applied. Feed and cooling temperature are selected as decision variables together with distribution of non‐uniform catalysts along the bed. The optimization result shows 14.47 % increase in the C₅⁺ yield with respect to the base condition.     

Production and Isolation of n‐Butyl Acrylate Using Pervaporation‐Aided Esterification Reaction: Kinetics and Optimization

Synthesis of n‐butyl acrylate by esterification of acrylic acid with n‐butyl alcohol was carried out in a batch membrane reactor. Optimization and design of the experiment was accomplished by response surface methodology with Box‐Behnken experimental design. The effects of different parameters like reaction temperature, catalyst concentration, molar ratio of alcohol to acid, and ratio of membrane surface to initial volume on water flux and conversion of acrylic acid were evaluated. A kinetic model for the esterification‐coupled pervaporation process was developed. Kinetic parameters were estimated by a nonlinear optimization technique in the MATLAB optimization toolbox. The experimental and simulation results were applied for developing a concept to effectively conduct a pilot‐scale esterification‐pervaporation experiment.     

Advanced‐TEMKIN Reaktor: Testung von industriellen Schalenkatalysatoren im Labormaßstab

In a variety of reactions in the chemical industry egg‐shell catalysts with a thin active layer are applied, which are often crushed for laboratory testing. The destruction of the shell can be avoided by a special reactor design. The presented Advanced TEMKIN reactor is a further development of the reactor system for testing egg‐shell catalysts in laboratory scale published by Temkin et al. in 1969. It is suitable for kinetic studies as well as for the detailed research of deactivation processes, as is shown on the example of the selective hydrogenation of acetylene.     

Heterogeneous catalytic reactor design with optimum temperature profile II: application of non-uniform catalyst

A new methodology has been developed to design non-isothermal, non-adiabatic heterogeneous catalytic fixed bed and tubular reactors with optimal temperature profiles inside a reactor. Catalyst characteristics such as pellet diameter, shape and activity distributions inside a pellet are considered simultaneously for reactor design. Various types of non-uniform activity distributions inside a pellet are modelled and optimised for the maximisation of an objective such as yield or selectivity. Dirac-δ, layered and general non-uniform distribution profiles such as egg-shell, egg-yolk and middle peak distributions are applied for the reactor design. The research demonstrates that different catalyst distribution profiles can approach the optimum performance. Whilst it is known that the Dirac-δ profile (and its step-function equivalent) always gives the best performance for clean catalyst, other profiles can approach this performance and might offer advantages in catalyst manufacture and under degraded conditions. A profile-based synthesis approach is applied to generate various shapes of activity profiles for multiple sections along the reactor during the optimisation of non-uniform catalyst pellets. A case study with the ethylene oxidation process illustrates that the catalyst characteristics, such as activity distribution profiles inside a pellet, sizes and shapes can be manipulated to control the temperature through the reactor very effectively, leading to significant improvements in selectivity or yield. The non-uniform catalyst pellet is further applied to various reactor configurations such as inert mixing and side stream distributions. This work is the first to consider all of these effects simultaneously.     

Fischer‐Tropsch Synthesis: Modeling and Performance Study for Fe‐HZSM5 Bifunctional Catalyst

A water‐cooled fixed bed Fischer‐Tropsch reactor packed with Fe‐HZSM5 catalyst has been modeled in two dimensions (radial and axial) using the intrinsic reaction rates previously developed at RIPI. The reactor is used for production of high‐octane gasoline from synthesis gas. The Fischer‐Tropsch synthesis reactor was a shell and tube type with high pressure boiling water circulating on the shell side. By the use of a two‐dimensional model, the effects of some important operating parameters such as cooling temperature, H₂/CO ratio in syngas and reactor tube diameter on the performance capability of the reactor were investigated. Based on these results, the optimum operating conditions and the tube specification were determined. The model has been used to estimate the optimum operating conditions for the pilot plant to be operated in RIPI.     

An optimal catalyst activation policy for poisoning problems

Pontryagins minimum principle is used to calculate the optimum distribution of active material throughout a single pellet that is uniformly experiencing activity decay. The definition of optimality is based upon balancing catalyst costs against the net return from reactant conversion. The optimal policy is full activation to a fractional depth with an inert core, the depth depending upon the Thiele parameter, poisoning time constant, operating time, and an economic parameter. An isothermal, first order, irreversible reaction is considered.Auxiliary calculations of the effectiveness factor are given for the reaction rate in the nonisothermal pellet with the catalyst distribution found to be optimal in the isothermal case.The results of the pellet calculations are utilized in numerically calculating a uniform activation policy that is optimal for a homogeneously poisoned bed. The definition of optimality is on the same basis as that for the single pellet and the results depend upon the same physical parameters in addition to the time constant for convection in the tube relative to that for diffusion in the pellets. Four regimes of behavior arise depending upon catalyst costs: (i) the reactor cannot be operated economically, (ii) only partial activation policies are economical, a single one being optimal, (iii) both partial and full activation schemes are economical, with one of the former being optimal, and (iv) full activation is optimal but the reactor can be operated economically with partially activated pellets.The results for both the single pellet and the packed tube are viewed as the homogeneous poisoning limit of many practical problems and are believed to reflect the major characteristics of more complicated poisoning mechanisms.     

CO2 methanation: Optimal start‐up control of a fixed‐bed reactor for power‐to‐gas applications

Utilizing volatile renewable energy sources (e.g., solar, wind) for chemical production systems requires a deeper understanding of their dynamic operation modes. Taking the example of a methanation reactor in the context of power‐to‐gas applications, a dynamic optimization approach is used to identify control trajectories for a time optimal reactor start‐up avoiding distinct hot spot formation. For the optimization, we develop a dynamic, two‐dimensional model of a fixed‐bed tube reactor for carbon dioxide methanation which is based on the reaction scheme of the underlying exothermic Sabatier reaction mechanism. While controlling dynamic hot spot formation inside the catalyst bed, we prove the applicability of our methodology and investigate the feasibility of dynamic carbon dioxide methanation. © 2016 American Institute of Chemical Engineers AIChE J, 63: 23–31, 2017     

Selective Hydrogenation of 1‐Butene Rich Cuts: the Impact of the Intraparticle Diffusion Limitations on the Selectivity

The impact of intraparticle diffusion limitations on the selectivity of an industrial reactor for selective hydrogenation of 1‐butyne and 1,3‐butadiene contained in 1‐butene rich cuts was evaluated. To this end, a simple model of a trickle‐bed reactor was employed and actual process operating conditions were chosen. A kinetic model was chosen whose parameters correspond to a commercial catalyst. These parameters were calculated from experiments conducted under industrial operating conditions. The complex diffusion and reaction phenomena occurring inside catalyst pellets placed at different depths of the reactor are comprehensively described. 1‐Butene losses in the range 20–30 %, which are usual in commercial plants, were predicted. It was concluded that the operating pressure is crucial for enhancing process selectivity.     

Numerical Investigation of Resonance Behaviour of a Tubular Packed‐Bed Reactor with Non‐Uniform Activity

The effects of periodic temperature perturbations at the inlet of a packed‐bed catalytic reactor were investigated for a system of non‐uniform catalyst activity that could arise through catalyst deactivation. A simplified model consisting of zones of active catalyst and inert packing of similar thermal properties was assumed. Simulations showed that large amplification, expressed as gain, can occur. The largest gain is observed at a resonance frequency. Both gain and resonance frequency depend on the number of layers of active catalyst and inert packing, the depth of the inert packing and catalyst layers.     

Simulation of fixed bed reactor for dimethyl ether synthesis

Dimethyl Ether (DME) is considered as one of the most promising candidates for a substitute for LPG and diesel fuel. We analyzed one-step DME synthesis from syngas in a shell and tube type fixed bed reactor with consideration of the heat and mass transfer between catalyst pellet and reactants gas and effectiveness factor of catalysts together with reactor cooling through reactor wall. Simulation results showed strong effects of pore diffusion. We compared two different arrangements of catalysts, mixture of catalyst pellets (methanol synthesis catalyst and methanol dehydration catalyst) and hybrid catalyst. Hybrid catalyst gave better performance than a mixture of pellets in terms of CO conversion and DME productivity, but more difficulties with reactor temperature control. Use of inert pellets and inter-cooling was also simulated as a means of controlling maximum reactor temperature.     

Industrial Xylene/Ethylbenzene Isomerization Unit Using a Radial‐Flow Reactor and EUO‐Type Zeolite

A mathematical model based on 34 days continuous operation of an industrial isomerization unit was developed. The unit involves a radial‐flow reactor with a catalyst capable of converting xylenes and ethylbenzene to mixed xylenes. The catalyst contains EU‐1 zeolite, platinum, and alumina as binder. Two reactions are considered, i.e., ethylbenzene isomerization and xylene isomerization. The rates are based on the Hougen‐Watson model according to the literature. An optimization procedure based on the trust‐region‐reflective algorithm was carried out in order to obtain new kinetic constants that minimize the difference between the actual and the calculated values. The standard error of the parameters estimated was calculated through the deleted‐one Jackknife method.     

Analysis of the fluid–structure interaction in the optimization‐based design of polymer sheeting dies

A polymer‐sheeting‐die‐design methodology is presented that integrates a simulation of the polymer melt flow and die‐cavity deformation with numerical optimization to compute a die‐cavity geometry capable of giving a nearly uniform exit flow rate. Both the polymer melt flow and sheeting‐die deformation are analyzed with a general‐purpose finite‐element program. The approach includes a user‐defined element that is used to evaluate the purely viscous non‐Newtonian flow in a flat die. The flow analysis, which is simplified with the Hele–Shaw approximation, is coupled to a three‐dimensional finite‐element simulation for die deformation. In addition, shape optimization of a polymer sheeting die is performed by the incorporation of the coupled analyses in our constrained optimization algorithm. A sample problem is discussed to illustrate the die‐design methodology. © 2006 Wiley Periodicals, Inc. J Appl Polym Sci 103: 3994–4004, 2007     

Multiscale modeling of methane catalytic partial oxidation: From the mesopore to the full‐scale reactor operation

A multiscale methodology combining three different reactor length‐scales is presented to investigate the role of the catalyst internal pore structure and metal loading and dispersion on the catalyst layer and full‐scale reactor performances. At the catalyst level, the methodology involves pore‐scale simulations in the three‐dimensional mesopore and macropore space. The information gathered at the catalyst level is delivered to the full‐scale reactor model. The methodology is applied to a honeycomb reactor performing methane partial oxidation considering reaction kinetics described through a detailed multistep reaction mechanism. Realistic mesopore and macropore structures were reconstructed and combined to form specific bidisperse porous washcoat layers. The study shows that species effective diffusivities vary significantly but not in the same proportion for different structures. For structures featuring poor transport characteristics, the integral methane conversion and hydrogen selectivity are strongly affected while the reactor temperatures increase substantially. © 2017 American Institute of Chemical Engineers AIChE J, 64: 578–594, 2018     

Dynamic Simulation and Optimization of a Dual‐Type Methanol Reactor Using Genetic Algorithms

In this investigation, a dynamic simulation and optimization for an auto‐thermal dual‐type methanol synthesis reactor was developed in the presence of catalyst deactivation. Theoretical investigation was performed in order to evaluate the performance, optimal operating conditions, and enhancement of methanol production in an auto‐thermal dual‐type methanol reactor. The proposed reactor model was used to simulate, optimize, and compare the performance of a dual‐type methanol reactor with a conventional methanol reactor. An auto‐thermal dual‐type methanol reactor is a shell‐and‐tube heat exchanger reactor in which the first reactor is cooled with cooling water and the second one is cooled with synthesis gas. The proposed model was validated against daily process data measured of a methanol plant recorded for a period of 4 years. Good agreement was achieved. The optimization was achieve by use of genetic algorithms in two steps and the results show there is a favorable profile of methanol production rate along the dual‐type reactor relative to the conventional‐type reactor. Initially, the optimal ratio of reactor lengths and temperature profiles along the reactor were obtained. Then, the approach was followed to get an optimal temperature profile at three periods of operation to maximize production rate. These optimization approaches increased by 4.7 % and 5.8 % additional yield, respectively, throughout 4 years, as catalyst lifetime. Therefore, the performance of the methanol reactor system improves using optimized dual‐type methanol reactor.     

Catalyst performance for noble metal catalysed alcohol oxidation: reaction-engineering modelling and experiments

A reaction-engineering model is presented, which describes catalyst performance as a function of the catalyst activity profile, the reaction kinetics, and the degree of catalyst deactivation. With this model, the catalyst activity profile can be optimised for Pt catalysed methyl -D-glucopyranoside (slowly-reactive) and glucose (highly-reactive) oxidations. This is done by comparing modelling results with experimentally obtained data for catalysts of different activity distributions. Experiments in a semi-batch stirred reactor showed that for methyl -D-glucopyranoside (MGP) oxidation at oxygen partial pressures below 40 kPa, egg shell catalytic activity distribution gives a higher rate of oxidation than a uniform distribution. It was also observed that with increase in oxygen concentration from 10 to 40 kPa, the rate of deactivation due to catalyst over-oxidation increased dramatically. For glucose oxidation, both catalyst activity distributions give the same oxidation rate for all investigated oxygen partial pressures (5–100 kPa). The developed model adequately describes the observed experimental results of both reactions. It was found that the active metal particle size has a significant influence on the catalyst deactivation for MGP oxidation; the uniform catalyst with higher dispersion shows a higher deactivation rate than the egg shell catalyst. For modelling glucose oxidation, the effect of catalyst particle-to-bubble adhesion and higher diffusivity or partition coefficient for oxygen have to be taken into account.     

Isothermal effectiveness factor—II: Analytical expression for single reaction with arbitrary kinetics,geometry and activity distribution

A previously presented method[1], to predict isothermal effectiveness factor with a single complex reaction in isothermal slab pellets is extended to encompass much more complex and general situations. In this work the geometry of the pellet can be arbitrary and a non uniform distribution of the catalyst is considered.Though the previously presented method[1] has had to be slightly modified to predict with great accuracy the effectiveness factor, (with less than 3% deviation from the exact values), an almost general and very simple algebraic expression is deduced to predict effectiveness factor values within 10% of their respective exact values. Thus for many applications in engineering design and catalytic reactor simulation, this simple general expression can be extremely useful since only one easily generated parameter is needed, as shown throughout the present contribution.     

High‐Pressure High‐Temperature Transparent Fixed‐Bed Reactor for Operando Gas‐Liquid Reaction Follow‐up

A 3‐MPa, 350 °C fixed‐bed reactor was designed to follow‐up gas‐liquid‐solid reactions on a millimetric size heterogeneous catalyst with Raman spectroscopy. The transparent reactor is a quartz cylinder enclosed in a Joule effect heated stainless‐steel tube. A methodology to determine how to focus the microscope for liquid and solid phase characterization is presented. The setup was validated by performing diesel hydrodesulfurization on a CoMo/alumina extrudate catalyst with a conversion very close to expected values along with the acquisition of Raman spectra of the solid catalyst showing an evolution of the catalyst phase during sulfidation.     

Coupling non-isothermal trickle-bed reactor with catalyst pellet models to understand the reaction and diffusion in gas oil hydrodesulfurization

In this work, a trickle-bed reactor coupled with catalyst pellet model is employed to understand the effects of the temperature and catalyst pellet structures on the reaction–diffusion behaviors in gas oil hydrodesulfurization(HDS). The non-isothermal reactor model is determined to be reasonable due to non-negligible temperature variation caused by the reaction heat. The reaction rate along the reactor is mainly dominated by the temperature,and the sulfur concentration gradient in the catalyst pellet decreases gradually along the reactor, leading to the increased internal effectiveness factor. For the fixed catalyst bed volume, there exists a compromise between the catalyst reaction rate and effectiveness factor. Under commonly studied catalyst pellet size of 0.8–3 mm and porosity of 0.4–0.8, an optimization of the temperature and catalyst pellet structures is carried out, and the optimized outlet sulfur content decreases to 7.6 wppm better than the commercial level at 0.96 mm of the catalyst pellet size and 0.40 of the catalyst porosity.     

A novel approach to the design and operation scheduling of heterogeneous catalytic reactors

A number of studies have been conducted to reduce the overall level of catalyst deactivation in heterogeneous catalytic reactors, and improve the performance of reactors, such as yield, conversion or selectivity. The methodology generally includes optimization of the following: (1) operating conditions of the reaction system, such as feed temperature, normal operating temperature, pressure, and composition of feed streams; (2) reactor design parameters, such as dimension of the reactor, side stream distribution along the axis of the reactor beds, the mixing ratio of inert catalyst at each bed; and (3) catalyst design parameters, such as the pore size distribution across the pellet, active material distribution, size and shape of the catalyst, etc. Few studies have examined optimization of the overall catalyst reactor performance throughout the catalyst lifetime, considering catalyst deactivation. Furthermore, little attention has been given to the impact of various configurations of reactor networks and scheduling of the reactor operation (i.e., online and offline-regeneration) on the overall reactor performance throughout the catalyst lifetime. Therefore, we developed a range of feasible sequences of reactors and scheduling of reactors for operation and regeneration, and compared the overall reactor performance of multiple cases. Furthermore, a superstructure of reactor networks was developed and optimized to determine the optimum reactor network that shows the maximum overall reactor performance. The operating schedule of each reactor in the network was considered further. Lastly, the methodology was illustrated using a case study of the MTO (methanol to olefin) process.