Elsevier

Computers & Chemical Engineering

Volume 36, 10 January 2012, Pages 35-47
Computers & Chemical Engineering

Improved single phase modeling of propylene polymerization in a fluidized bed reactor

https://doi.org/10.1016/j.compchemeng.2011.07.015Get rights and content

Abstract

An improved model for the production of polypropylene in a gas phase fluidized bed reactor was developed. Comparative simulation studies were carried out using the well-mixed, constant bubble size and the improved models. The improved model showed different prediction characteristics of polymer production rate as well as heat and mass transfer behavior as compared to other published models. All the three models showed similar dynamic behavior at the startup conditions but the improved model predicted a narrower safe operation window. Furthermore, the safe ranges of variation of the main operating parameters such as catalyst feed rate and superficial gas velocity calculated by the improved and well mixed models are wider than that obtained by the constant bubble size model. The improved model predicts the monomer conversion per pass through the bed which varies from 0.28 to 5.57% within the practical ranges of superficial gas velocity and catalyst feed rate.

Highlights

► The improved model showed different prediction characteristics of polymer production rate as well as heat and mass transfer behavior as compared to other published models. ► The well-mixed, constant bubble size and the improved models showed similar dynamic behavior at the startup conditions but the improved model predicted a narrower safe operation window. ► The improved model predicts the monomer conversion per pass through the bed which varies from 0.28 to 5.57% within the practical ranges of superficial gas velocity and catalyst feed rate.

Introduction

The various advantages of fluidized bed reactors (FBR) such as their ability to carry out a variety of multiphase chemical reactions, good mixing of particles, high rate of mass and heat transfer and their ability to operate in continuous state have made it one of the most widely used reactors for polyolefin production. Consequently, considerable attention has been paid to model propylene polymerization in fluidized bed reactors. Fig. 1 illustrates the schematic of an industrial gas-phase fluidized bed polypropylene reactor. As shown in this figure, small particles of Ziegler–Natta catalyst and triethyl aluminum co-catalyst are charged continuously to the reactor and react with the reactants to produce a broad distribution of polymer particles. The catalyst particles are porous and are composed of small sub-fragments containing titanium active metal. As the monomer diffuses through the porous catalyst, it polymerizes by reaction on the active sites of the catalyst surface. The catalyst fragments become dispersed during the polymerization and the particles grow into the final polymer product (Zacca, Debling, & Ray, 1996). The feed gas, which consists of propylene, hydrogen and nitrogen, provides the fluidization through the distributor, acts as the heat transfer medium and supplies the reactants for the growing polymer particles. The fluidized particles disengage from the unreacted gases in the disengaging zone. The solid-free gas is combined with fresh feed after heat removal and recycled back to the gas distributor. The polypropylene product is continuously withdrawn from near the base of the reactor and above the gas distributor. The unreacted gas is recovered from the product which proceeds to the finishing section of the plant.

In heterogeneous systems, polymerization occurs in the presence of different phases with inter-phase mass and heat transfer and chemical reaction. Phenomena such as complex flow of gas and solids, kinetics of heterogeneous polymerization and various heat and mass transfer mechanisms must be incorporated in a realistic modeling approach. Several different methods for describing the hydrodynamics of the fluidized bed polyolefin reactor have been proposed in the literature. McAuley, Talbot, and Harris (1994) and Xie, McAuley, Hsu, and Bacon (1994)considered the fluidized bed polyolefin reactor as a well mixed reactor. They compared the simple two-phase and the well mixed models at steady state conditions and showed that the well mixed model does not exhibit significant error in the prediction of the temperature and monomer concentration in the reactor as compared with the simple two-phase model at steady state conditions. Choi and Ray (1985) presented a simple two-phase model in which the reactor consists of emulsion and bubble phases. They assumed that the polymerization reaction occurs only in the emulsion since the bubbles are solid-free. Fernandes and Lona (2001) proposed a heterogeneous three-phase model that considers bubble, emulsion and particulate phases with plug flow behavior. Hatzantonis, Yiannoulakis, Yiagopoulos, and Kiparissides (2000) considered the reactor being comprised of perfectly mixed emulsion phase and a bubble phase divided into several solid-free well-mixed compartments in series. Alizadeh, Mostoufi, Pourmahdian, and Sotudeh-Gharebagh (2004) proposed a tanks-in-series model to represent the hydrodynamics of the reactor. Harshe, Utikar, and Ranade (2004) developed a comprehensive mathematical model based on the mixing cell framework to simulate transient behavior in the fluidized bed polypropylene reactors. This dynamic model was coupled with a steady state population balance equation. Rigorous multi-monomer, multisite polymerization kinetics was incorporated in the model. Ibrehem, Hussain, and Ghasem (2008) proposed a fluidized bed comprising of bubble, cloud, emulsion and solid phases and considered the polymerization reactions occurring in the emulsion and solid phases. Their model also accounts for the effect of catalyst particles type and porosity on the rate of reaction.

In the present study, these previous works on the modeling of gas-phase olefin polymerization fluidized-bed reactors were extended to account for the dynamic behavior of propylene polymerization in fluidized bed reactors. To describe the homopolymerization of propylene over a heterogeneous Ziegler–Natta catalyst, a two-site homopolymerization kinetic scheme was employed. Extensive simulations were carried out to determine the influence of key process parameters (such as superficial gas velocity, catalyst feed rate, feed composition and feed temperature) on the dynamic response of the reactor. In order to obtain a better understanding of the reactor performance, a comparative study between the well-mixed, constant bubble size and improved model was carried out.

Section snippets

Polymerization mechanism

Similar methodology as of McAuley, MacGregor, and Hamilec (1990), Kissin (1985), and de Carvalho, Gloor, and Hamielec (1989) were employed and the following mechanism scheme was developed for the homopolymerization of propylene over the Ziegler–Natta catalyst containing multiple active sites. Mass and heat transfer resistances and existence of multiple active sites are two factors responsible for exhibiting broad molecular weight distribution using Ziegler–Natta catalyst. It has been shown that

Kinetic model

A mathematical kinetic model based on the mechanism described in the previous section was derived in the present study. This model consists of mass balances on the species presented in the reactor which are represented by a series of algebraic and differential equations as described in the next section. Characterization of polymer properties was modeled using the population balance and the method of moments. Application of population balance and the method of moments allows for the prediction

Modeling of fluidized bed reactor

Complex mixing and contacting flow patterns, transport phenomena and polymerization reactions make the fluidized bed reactors non-ideal and difficult to characterize. Many studies have attempted to model such non-ideality using various mixing models to describe the behavior of fluidized bed reactors. A combination of kinetics, hydrodynamics and transport phenomena is required for modeling the non ideal fluidized bed reactor.

The simple two-phase flow structure for the gas-phase olefin

Results and discussion

To investigate the effect of different hydrodynamic sub-models on the dynamic response and grade transition, dynamic modeling of the process was carried out using the three hydrodynamic sub-models described above along with the reaction sub-model. Comprehensive simulations were carried out in order to evaluate the effect of the key process parameters such as catalyst feed rate, superficial gas velocity, propylene feed concentrations and carbon monoxide on the polymer production rate, emulsion

Conclusions

An improved dynamic model for the production of polypropylene in a gas phase fluidized bed reactor was developed to describe the behavior of the fluidized bed reactor of polypropylene production. Comparative simulation studies were carried out using the well-mixed, constant bubble size and the improved models in order to investigate the effects of the mixing, operating conditions, kinetic and hydrodynamic parameters on the reactor performance. The improved model showed different prediction

Acknowledgments

The authors would like to thank the support of the Research Council of the University of Malaya under research grant (No. RG054/09AET) and the Petroleum and Chemical Engineering Department at Sultan Qaboos University and the School of Chemical Engineering, College of Engineering, University of Tehran for their support to this research.

References (27)

  • J.J. Zacca et al.

    Reactor residence time distribution effects on the multistage polymerization of olefins-I. Basic principles and illustrative examples, polypropylene

    Chemical Engineering Science

    (1996)
  • K.Y. Choi et al.

    The dynamic behavior of fluidized bed reactors for solid catalyzed gas phase olefin polymerization

    Chemical Engineering Science

    (1985)
  • H.P. Cui et al.

    Characterization of dynamic gas–solid distribution in the fluidized beds

    Chemical Engineering Journal

    (2000)
  • Cited by (32)

    • Modeling of the molecular weight distribution and short chain branching distribution of linear low-density polyethylene from a pilot scale gas phase polymerization process

      2022, Chemical Engineering Science
      Citation Excerpt :

      This reaction scheme (designated as reaction scheme 1) is introduced in order to be compared with a non-standard reaction scheme which considers the “comonomer effect” (designated as reaction scheme 2) and which is adopted in this work. The reactor model is based on the following assumptions: all components are well mixed within the fluidized bed reactor (McAuley et al., 1994; Shamiri et al., 2012); their concentrations and temperature are uniform throughout; the mass and heat transfer resistances between the gas and the solid particles are negligible. The uniform temperature in the reactor is confirmed by the fact that the temperature rise measured inside the reactor is only 2 ∼ 3 K; because the single-pass conversion of ethylene is less than 1 %.

    • Predictive two-phase modeling and improvement for industrial Unipol<sup>Ⓡ</sup> polypropylene process – Development and a case study

      2020, Computers and Chemical Engineering
      Citation Excerpt :

      Thus, these parameters were taken from the reported values in the literature (Luo et al., 2009; Khare et al., 2004; Zheng et al., 2011). On the other hand, the kinetics of propylene homopolymerization using multi-site Z-N catalysts has been extensively studied in the literature (Luo et al., 2009; Shamiri et al., 2010; 2011; 2012; Khare et al., 2004). In this model, six elementary reaction-type categories were considered commensurate with the current case.

    View all citing articles on Scopus
    View full text