1. Introduction
Polyolefins are the most important plastics produced in the world. Although polyolefins have the largest market share in the plastics industry, competition in this market is hardened with the regular construction of new plants around the world, prompting old plants to reduce their operating costs in order to remain competitive. The demand for polyolefins is growing, taking into account specific needs. Polyolefin manufacturers that can respond to this demand guarantee a market share with higher profits.
Polypropylene production is the largest petrochemical production at present. The demand for this high-quality and inexpensive polymer is growing, so it is gradually replacing polyethylene. It is indispensable for creating rigid packaging, plates, films, automotive parts, synthetic paper, ropes, carpets, as well as for creating a variety of domestic equipment. At the beginning of the twenty-first century, the production of polypropylene occupied the second place in the polymer industry. Taking into account the demands of various industries, it can be concluded that in the near future the trend of large-scale production of propylene and ethylene will continue.
Currently, the bulk of polypropylene in the world is produced using titanium-containing Ziegler–Natta catalysts. Titanium–magnesium catalysts of the IV generation provide a polymer yield during the polymerization in a liquid monomer of 30–60 kg per 1 kg of catalyst with polypropylene isotacticity of more than 96 wt. %. However, in Russia today there is no industrial production of catalysts of the IV generation. All domestic polypropylene producers traditionally use imported catalysts. Within the framework of import substitution, the task is to produce and use domestic olefin polymerization catalysts.
Usually, the production of polypropylene is carried out using mathematical models of the production process [
1], presented in the form of simulators. Scientific research continues on modeling olefin polymerization reactions using domestic titanium–magnesium catalysts. In addition to well-founded reactor models, reliable kinetic parameters are needed for the correct use of these models and simulators. If reliable estimates of kinetic parameters are not available, engineers are forced to rely on rough assumptions, which can lead to undesirable uncertainty when making decisions in the production process.
Polymerization of propylene in the liquid phase is one of the most important industrial processes in the production of polypropylene. This process in liquid monomer accounts for about 35% of the total polypropylene production.
Polymerization process: a catalyst (titanium–magnesium Ziegler–Natta catalyst of the fourth generation), pre-activated with alkyl aluminum and treated with an external electron donor, is continuously fed into the first reactor together with a monomer and a chain transfer agent (hydrogen). The flow coming out of the first reactor, after mixing with the side feed flow containing a monomer and a chain transfer agent, is fed to the second loop reactor. The flow coming out of the second reactor is processed to separate the solid-phase polymer from the unreacted monomer. The solid phase of the polymer particles is deactivated by steam and then extruded. In the case of heterogeneous (impact-resistant) copolymers, polymerization is continued in a gas-phase reactor, where the addition of ethylene-propylene rubber provides better impact resistance. The loop reactor operates completely filled with the suspension. The axial pump located in the lower part of the reactor promotes high-speed (5–7 m/s) recirculation of the reaction mixture. The resulting turbulent flow ensures efficient heat transfer between the cooling water of the jacket and the reaction medium.
Although the process can operate with a single loop reactor, it is usually designed using two loop reactors. One or two additional gas-phase reactors can be placed in series to produce heterogeneous copolymers. The two-loop reactor system has the ability to produce bimodal polymers with enhanced mechanical properties.
In the works of Alshaiban A. et al. [
2], the kinetic scheme of the polymerization process using Ziegler–Natta catalysts is most fully presented. It can be concluded that when creating a model of the industrial polymerization process of propylene on a titanium–magnesium catalyst, a simplified kinetics model is usually used.
2. Features of the Process of Polymerization of Propylene on Titanium–Magnesium Catalysts
Despite the fact that heterogeneous Ziegler–Natta catalysts have been used in industry for almost half a century, understanding the differences between the types of active centers that are presented in heterogeneous titanium–magnesium catalysts is a difficult task [
3,
4,
5,
6]. Due to their interaction with co-catalysts, internal and external donors, resistance to the transfer of mass and heat between and inside the particles during the polymerization, the effect of hydrogen on the polymerization rate, as well as sensitivity to poisons because of this complexity, it is difficult to obtain polyolefins with a controlled microstructure using titanium–magnesium catalysts [
7,
8,
9].
When the same heterogeneous Ziegler–Natta catalyst is combined with various external donors, it can produce polypropylene with different average molecular weights, polydispersity and tacticity under the same polymerization conditions [
10,
11,
12]. For example, the addition of an external silane donor, such as dicyclopentyldimethoxysilane, can double the average molecular weight of polypropylene obtained using TiCl
4/dialkyl phthalate/MgCl
2 [
13]. Various internal/external donor systems also affect the stereoregularity of polypropylene [
14,
15,
16].
Hydrogen increases the polymerization rate of propylene when using highly active catalysts based on TiCl
4 [
17,
18,
19]. Studies with TiCl
4/dialkyl phthalate/MgCl
2 have shown a significant increase in the polymerization rate with the addition of hydrogen during polymerization [
20]. Busico et al. [
21] explained this effect of increasing the rate by releasing 2,1-terminal dormant centers, which is further confirmed by analyses of the polymer terminal group [
22,
23]. In publications [
24,
25], it is also noted that hydrogen increases the tacticity of polypropylene.
Temperature is an important variable affecting the polymerization rate of propylene. The polymerization rate with most TiCl
4/MgCl
2 catalysts is the highest, from 60 to 70 °C, starting to decrease above 70 °C [
26,
27,
28]. The polymerization temperature also affects the stereoregularity of polypropylene [
29]. Polypropylene chains can be classified according to their crystallizability into three main fractions: low, medium and high-crystalline. Kissin et al. It has been shown that when the polymerization temperature increases, the fraction with a high degree of crystallinity increases, but the average number of mmmm pentads for the entire polymer sample decreases [
30,
31], which reflects the complex nature of the active centers of heterogeneous Ziegler–Natta catalysts.
The article [
32,
33] quantified the effect of an electron and hydrogen donor on the kinetics of propylene polymerization for a fourth-generation commercial Ziegler–Natta catalyst modified by an external electron donor. A quantitative calculation model has been developed, taking into account the influence of hydrogen and an external electron donor on the polymerization rate of propylene for this type of catalytic system (TiCl
4/MgCI
2 with a Ti content of 1.7% by weight, the co-catalyst is triethylaluminium, the electron donor is dicyclopentyldimethoxysilane).
Despite the fact that heterogeneous Ziegler–Natta catalysts have been used in industry for almost half a century, understanding the differences between the types of active centers that are in heterogeneous titanium–magnesium catalysts is a difficult task. That is, because of this complexity, it is difficult to obtain polyolefins with a controlled microstructure using titanium–magnesium catalysts [
32].
When the heterogeneous Ziegler–Natta catalyst is combined with different external donors, it can produce polypropylenes with different average molecular weights, polydispersity and tacticity under the same polymerization conditions [
33]. For example, the addition of an external silane donor such as dicyclopentyldimethoxysilane can double the weight average molecular weight of polypropylene produced with TiCl
4/dialkyl phthalate/MgCl
2 [
16]. Various internal/external donor systems also affect the stereoregularity of polypropylene [
17].
Hydrogen increases the rate of propylene polymerization when highly active catalysts based on TiCl
4 are used [
18,
19]. Researches about TiCl
4/dialkyl phthalate/MgCl
2 have shown a significant increase in the polymerization rate when hydrogen is added during the polymerization. Busico et al. [
20] explained this rate increase effect by releasing the 2,1-terminal dormant centers, which was additionally supported by polymer end group analyses. It was also noted that hydrogen increases the tacticity of polypropylene [
2,
21,
22].
Temperature is an important parameter affecting the rate of polymerization of propylene. The polymerization rate with TiCl
4/MgCl
2 catalysts has the highest values from 60 to 70 °C, starting to decrease above 70 °C [
23,
24,
25]. The polymerization temperature also affects the stereoregularity of polypropylene. Kissin et al. [
26] showed that when the polymerization temperature is increased, the fraction yield with a high degree of crystallinity increases, but the average number of pentads for the entire polymer sample decreases, which reflects the complex nature of the active sites of heterogeneous Ziegler–Natta catalysts.
Keii et al. [
23] presented the influence of an electron donor and hydrogen on the kinetics of propylene polymerization for a commercial fourth-generation Ziegler–Natta catalyst modified with an external electron donor. A quantitative calculation model has been developed, taking into account the effect of hydrogen and an external electron donor on the propylene polymerization rate for a given type of catalytic system (TiCl
4/MgCI
2 with a Ti content of 1.7 wt %, triethylaluminum as a co-catalyst, and dicyclopentyldimethoxysilane as an electron donor).
The highest activity values in the polymerization system are observed for temperatures of 70–90 °C, monomer concentration of 403 g/L in the presence of hydrogen (0.15 mol under the studied conditions). Under the specified polymerization conditions, polypropylene powder is characterized by optimal characteristics in terms of bulk density and morphology, i.e., high value of the average particle size and the absence of a dusty fraction. The isotacticity of PP slightly increases with increasing temperature, hydrogen concentration and the duration of the polymerization process. As for the molecular characteristics of the resulting PP, a decrease in the temperature and hydrogen content, an increase in the duration of polymerization and monomer content leads to an increase in average molecular weights and decrease the melt flow index. The molecular weight distribution (MWD) of PP samples slightly narrows with increasing temperature, duration of process and hydrogen concentration.
The results obtained in this work make it possible to choose the optimal conditions for the synthesis of propylene in liquid monomer for the catalytic system under study, and the identified process parameters effectively control the molecular and rheological characteristics of PP.
Ti-Mg catalysts (TMC) are of great importance in the production of polypropylene. Their typical composition is as follows: TiCl4/D1/MgCl2-AlEt3/D2, where the electron donor compounds D1 (internal donor) and D2 (external donor) provide high stereospecificity of the catalytic system.
Tsutsui T. et al. [
21] studied the chain transfer reaction with hydrogen during the polymerization of propylene on Ti-Mg catalysts of the TiCl
4/D1/MgCl
2-AlEt3/D2 composition in a wide range of hydrogen concentrations. A two-step reaction mechanism was proposed. This mechanism takes into account the fractional order of the reaction with respect to the hydrogen concentration. The reaction constants for the chain transfer reaction with hydrogen are determined for TMC with various donors: 1,3-diester or dibutyl phthalate as D1 and tetra-ethoxysilane or dicyclopentyldimethoxysilane as D2. In the polymerization of propylene with the TMC, the length of the polymer chain is mainly determined by the ratio of propylene and hydrogen concentrations, since the rate constants of chain propagation and transfer are comparable. The rate constant of chain transfer with hydrogen in the polymerization of ethylene is much lower, and higher concentrations of hydrogen are required to achieve the same degree of polymerization.
Despite the fact that the model is idealized and ignores the presence of different types of centers and mass transfer resistances inside the particles, it is able to estimate the apparent kinetic constants of polymerization, which adequately describe the kinetics of polymerization, and it can be included in the process simulator for better prediction of polymer properties, and used to develop more economical transitions between brands.
The need to create a mathematical model of the process consists in the high demand for obtaining high-quality plastic with specified physical and mechanical properties (melt flow rate, Izod impact strength, elastic modulus, frost resistance). Consequently, it is reasonable, on the basis of mathematical modelling of the process, to analyse the kinetic laws and molecular weights of the resulting polymer with the aim of further practical application, in particular, to identify a rational thermal regime for conducting the copolymerization process. A mathematical model of the copolymerization of propylene and ethylene in a reactor with a fluidized bed was created, taking into account the main reactions, the identification of the kinetic constants of the reaction rates, numerical experiments based on the obtained models of chemical kinetics. Research has been carried out and mathematical models have been developed for the chemical kinetics of propylene homopolymerization processes in a liquid monomer medium on titanium–magnesium catalysts.
3. Kinetic Scheme of the Polymerization Process
Experimental data show that the introduction of hydrogen makes it possible to increase the activity by nearly 40%. This increase in activity is the result of the transition of centers that were deactivated by regioirregular insertion into an active state. After regular insertion of propylene through the Ti–C bond the centers remain in the active state, while upon regioirregular insertion, the active centers pass into an inactive “dormant” state due to steric hindrance between propylene and the methyl group of the growing end of the chain at the titanium atom. After the chain transfer with such centers, hydrogen forms Ti–H bonds, which form active Ti–C bonds upon subsequent insertion of propylene. As a result, in the presence of hydrogen, the center passes from the dormant state to the active state (
Figure 1).
An increase in the hydrogen concentration in the system leads to the expected decrease in the molecular weight of polyplopylene and to an increase in the MFI (melt flow index) values due to the chain transfer reaction with hydrogen.
Table 1 presents the mechanism of the polymerization reaction of polypropylene on a titanium–magnesium catalyst.
The following notation is used:
—a living chain of polymerization degree r with an active centre of type j;
—a dead chain of polymerization degree r synthesized on an active centre of type j;
—a potential active center of type j;
A—alkyl aluminum co-catalyst;
Don—electron donor.
Speed constants are as follows:
—activation;
—growth;
—transfer of the chain to hydrogen;
—transfer of the chain to the monomer;
—deactivation;
, —transformations of active centers.
The kinetic scheme does not explicitly identify the initiation stage, since calculations usually consider the initiation rate constant to be equal to the growth rate constant, so these stages are combined.
The mathematical model being compiled includes a stage of reversible transformation of a stereospecific active centre of the first type into a non-specific active centre of the second type in the presence of an electron donor (Don). When the type of centre changes from 2 to 1 (coordination with the electron donor molecule) or from 1 to 2 (when the electron donor molecule is released), the length of the polymer chain does not change (r remains unchanged), but the number of stereo blocks increases by 1 (i + 1).
Due to the complexity of describing heterogeneous systems, due to the need to take into account many parameters, the mathematical model for describing the ethylene polymerization process will be supplemented by a number of equations taking into account the effect of the growth of a polymer particle on the nature of the flow structure, the effect of temperature on the growth and transfer constants of the chain, and a change in catalyst activity.
Molar balances of reagents entering the reactor are as follows:
Molar balances of vacant AC of two types are as follows:
To reduce the number of equations included in the polymerization process model, the method of moments was used. The moments of the distributions of living and dead chains were entirely determined by the formulas:
where n is the order of moments, j is the type of active centres.
Equations of the moments of the distribution of live chains with two types of AC are as follows:
Equations of the moments of the distribution of dead chains of two types are as follows:
Formulas for calculating average molecular weights and polydispersity coefficient are as follows:
where
- the molecular weight of the monomer link (
= 42 g/mol).
Equations of the moments of the distribution of isotactic chains are as follows:
The molar and mass fractions of purely isotactic chains are calculated by the formulas:
Equations of the moments of the distribution of atactic chains are as follows:
The molar and mass fractions of purely atactic chains are calculated by the formulas:
Equations of the moments of the distribution of stereoblock chains are as follows:
The molar and mass fractions of block chains are calculated using the formulas:
According to the compiled model (1)–(42), calculations were carried out for the following basic values of the model parameters:
T = 70 °C, kp1 = 6316 L/(mol min), kp2 = 631 L/(mol min),
ka1 = 2829 L/(mol min), ka2 = 2829 L/(mol min), km1 = 0 L/(mol min),
km2 = 0 L/(mol min), k21 = 150 L/(mol min), k12 = 0.01 min−1,
kH1 = 110 L/(mol min), kH2 = 110 L/(mol min), kd1 = 0.023 min−1,
kd2 = 0.023 min−1, = 0.000021 mol/(L min), = 0.000021 mol/(L min), = 1.96 mol/(L min), = 0.00144 mol/(L min), = 0.000042 mol/(L min), = 0.0000173 mol/(L min).
The concentration and type of donor play an important role in the polymerization of propylene, modeling their effect on the properties of PP is important for controlling transitions between brands to ensure that the required properties of the polymer are achieved in a short transition time. Controlling the concentration of electron donors is the most common way to control the tacticity of polypropylene synthesized at TMC.
Figure 2 and
Figure 3 show the effect of the initial electron donor concentration on the average molecular weights M
n, M
w and the polydispersity coefficient K
wn.
It can be seen from calculations based on the model that in the absence of an electron donor, stereoblock chains are not formed, only purely isotactic or atactic chains are synthesized. When an electron donor is introduced into the reactor, a significant part of the second type centres (non-specific) is transformed into the first type (isospecific) centres, reducing the proportion of atactic polypropylene, followed by an increase in the proportion of isotactic and stereoblock chains, as is shown in
Figure 6,
Figure 8,
Figure 9 and
Figure 10. According to the calculations, most stereoblock chains are diblocks with a very small proportion of multiblock chains (
Figure 10).
Figure 6 shows that the proportion of purely isotactic chains increases with increasing donor concentration. At the same time, the average molecular weight of isotactic chains decreases, which can be explained by the action of the donor as a pseudo-chain transfer agent, reducing the average length of the stereoblock.
When the initial donor concentration is reduced by half compared to the base value, the mass fraction of atactic polypropylene increases, and if the donor concentration doubles, the mass fraction of atactic polypropylene decreases. The mass fraction of stereoblock chains is also affected by the concentration of the electron donor, but diblock chains continue to be the dominant population among stereoblock chains for all calculations, according to the model.
Calculations show that the average values of molecular weights increase, and the polydispersity index decreases with an increase in the number of blocks in the chain. Both trends are expected, since longer chains will have a higher probability of the active centre transformation stage than shorter chains. The effect on Kwn is a simple consequence of sampling a narrow fraction of polymers: single-block chains follow Flory statistics with Kwn = 2, and chains with two, three or more blocks will have narrower MWDs, since they are selected from a narrower fraction with increasing molecular weight.
The influence of the donor type can also be studied by changing the values of the coefficients of the rate of transformation of active centers by the donor k12 and k21. The use of the best donors reduces the mass fraction of stereo block chains, since the transformation of isotactic centers into atactic centers is less likely during the lifetime of the polypropylene chain.
Similarly, the effect of changes in hydrogen concentration on the microstructure of polypropylene is studied. The molecular weight decreases with increasing hydrogen concentration. A change in the hydrogen concentration also affects the proportion of stereoblock chains: with an increase in the hydrogen concentration, polymer chains become shorter on average and, consequently, the probability of transformation of active centers occurring during the growth time of polymer chains decreases. Therefore, when the concentration of hydrogen increases, the proportion of stereoblock chains decreases.
4. Conclusions
To obtain stereoregular polypropylene with a high index of isotactic polypropylene (iPP), it is essential to carry out the polymerization by the ion-coordination mechanism on the Ziegler–Natta catalysts. The most effective catalysts are TiCl3-Al(C2H5)3 and TiCl3-Al(C2H5)2Cl.
Linear polymer can be formed both in a homogeneous and in a heterogeneous phase during the propylene polymerization by unclaimed isomers. This feature is important for PP when using a catalyst on a solid surface. Additionally, the size of the Ti ions and the difference between them are significant.
Determination of the monomer sorption rate on the catalyst surface results in the formation of a π-complex of the monomer with the TiCl3 surface and, as a result, the orientation of the monomer molecules before the incorporation into the catalytic complex. • TiCl3 exists in the form of various crystalline modifications, which differ from each other by different alternation in the crystal lattice of Ti and Cl, therefore, it has different accessibility of the Ti surface for the sorption of propylene molecules.
An increase in the pressure at T = const leads to an increase in the rate of polymerization. However, the solubility of propylene in the solvent will sharply increase with a high load and the consumption of propylene, which has not entered into thorough polymerization and processing of solvent vapors, will increase. The optimum pressure of the process is from 0.5–1.2 MPa.
When increasing the polymerization time, the polymer yield will increase and it results in worsening heat transfer and the solubility of PP in the solvent. The optimal duration of the process is from 4 to 6 h.
In practice, the molecular weight of the PP is regulated by the consumption of the hydrogen, by the electrolysis of H2O. Evaluation of efficiency of the chain length is carried out according to the melt flow index (MFI). The speed of the polymerization process increases by 30–40% when changing a hydrogen concentration of 0.5 to 2.5%. This can be explained by the maximum saturation of the TiCl3 surface and the formation of active centers.
Due to the presence of more than one type of AC in heterogeneous titanium–magnesium catalysts used in industry for the PP polymerization, the model can be expanded to include several types of centers. Different types of AC will be characterized by different kinetic parameters of polymerization, which can be estimated by deconvolution of experimental MWD.
Modeling of PP polymerization at TMC develops in two directions: the kinetics of the process is detailed and the imperfection of mixing in the polymerizer reactor is taken into account.