Modeling and Simulation of an Isothermal Suspension Polymerization Reactor for PMMA Production Using Python

This paper presents the modeling and simulation of a suspension polymerization for methyl methacrylate in an isothermal reactor to produce poly methyl methacrylate using Python 3.5. The numeral solution to the stiff ordinary differential equations was performed by building a custom module which was used with the inbuilt NumPy and matplotlib modules that come with the Anaconda python distro. Python was used in order to obtain a realistic solution that considers the gel, glass and cage effects that affect the non-linear polymerization kinetics established in literature. The results showed that a maximum monomer conversion of about 92.8% at a minimum batch time of about 2.2 hours could be achieved at the specified conditions to obtain a polydisperse polymer with an index of 27. It is further concluded that Python can be employed to perform similar studies with equal success as any other programming language.


Introduction
Poly (methyl methacrylate) (PMMA) is a transparent thermoplastic used in a wide range of fields including medicine, art, and aesthetics.It is also used as an alternative to glass due to its shatter-resistance property [1].There has been a widespread need to use PMMA based components in Uganda especially for construction and medicine purposes, however advancing the technology that addresses industrial polymerization processes is often hampered by a lack of fundamental research activities in academic and the few industrial labs available in Industrially, PMMA is obtained from the free-radical polymerization of methyl methacrylate [2] which may be produced via methanolysis of methylcrylamide sulphate [3].The polymerization reaction is commonly initiated by thermal radical forming agents often employing organic peroxy compounds and azo-compounds [4].Suspension polymerization is used for the commercial production of many important polymers other than PMMA such as polyvinyl chloride and polystyrene [5] with similar properties to those produced from bulk polymerization.The limitations of increased viscosity, reaction temperature, and volume contraction in the latter process are always avoided [6].The process may be described as being operated batch-wise however, in essence it is often in a semi-batch mode since some material enters the reactor after the start of polymerization.
During the process, the monomer is stirred with about twice its volume of water and dispersants forming droplet-like distribution of the monomer phase in water.The mechanism of polymerization corresponds to that of bulk polymerization and as such, the droplets become increasingly viscous during the course of polymerization.Distributors such as gelatin stabilize the suspension and prevent the droplets from adhering to one another during collisions.The system is heterogeneous therefore the polymer is collected as granular beads on completion of polymerization and it is washed to remove adhering distributor and salts then dried [7].
Various works have been published on the modeling of methyl methacrylate polymerization processes.A major feature of these studies is the observation in the increase of the mass reaction viscosity with monomer conversion resulting into a deviation from the "normal" kinetic.The severe reduction in the mobility of the macroradicals auto accelerates the reaction causing an increase in the gel effect and a decrease in the termination rate constant.If the polymerization is carried out at lower temperatures, the glass effect occurs when the transition state is reached at a certain conversion causing a decrease in the propagation rate constant.This results in the interruption of the reaction before the monomer is completely consumed.The increased viscosity at high monomer conversion leads to a decrease in the initiator efficiency due to the cage effect [8].In most of the reviewed works, gel and glass effects are modeled according to Chiu et al. [9] and the cage effect is modeled according to Achilias-Kiparissides [10].Most of the models considered are reported to have worked well for MMA polymerization [11], however the general model developed by Seth and Gupta [12] and reworked by Ghosh et al. [13] using free volume theory is noted as being most superior [8].Several simulation packages to assess the performance of the models under consideration have been used such as MATLAB [14], gPROMS [15], FORTRAN [16], which are based on utilization of sequential quadratic programming, ordinary differential equations and fuzzy-neural modeling.
From the above, it can be seen that techniques applied in modeling polymeri- should keep in mind that software cannot be a substitute for engineering judgement and its use without understanding is a dangerous abandonment of professional responsibility [17].
The aim of this paper is to model and simulate the polymerization of methyl methacrylate to PMMA in an isothermal suspension reactor.The reactor performance is studied by modeling and simulation using Python 3.5 instead of using experimental assessment.This is because the procedure requires low computational demand in performing the reactor studies and besides high-quality data is available for validation and comparison with literature values.

Model Formulation
The kinetics of MMA polymerization are described basing on the free radical polymerization mechanism scheme in Table 1 below.The scheme map comprises the initiation step, where the initiator (I) decomposes to form reactive radicals (R), the addition of monomer molecules (M) to the reactive polymer chain formed from the reactive radicals and the termination step where deactivation of polymer radicals occur.The initiator used is Benzoyl Peroxide (BPO).
P and D represent the live and dead polymer with n and m monomeric units respectively; whereas the kinetic constants for the initiator decomposition, initiation, propagation and termination are K d , K i , K p and K t respectively.
In this work, the process model in [12] [13] with slight adjustments was followed.The assumptions used in this work are: 1) The reactor operates isothermally at 70˚C 2) All reaction steps in the system are irreversible Table 1.Kinetic scheme for the free-radical polymerization of MMA. Initiation: The method of moments is invoked in order to reduce the infinite system of molar balance equations required to describe the molecular weight distribution [18] that would have presented roughly 10,000 -50,000 stiff differential equations to be solved simultaneously [19].Therefore the first three radical species and dead polymer moment balances are; where i λ is the momentum i of the radical species and i µ is the momentum i of the dead polymeric species summing to; , A detailed procedure of deriving the moment equations can be found elsewhere in the works of [20] [21].
The initial conditions used in this work: ( ) At any time and point within the system, predictions for the monomer conversion (X m ), number and weight average molecular weights (M n and M w ) can be given as; The enthalpy change Q during the reaction progress can be defined as In order to account for the cage, gel and glass effects, [12] reformulated the model of [10] basing on free volume theory and the equations below were used to define the variation of the initiator efficiency f, the termination and propagation rate constants with temperature ,0 In the above equations, the volume of the mixture V is determined from; With m ρ and p ρ as the densities of the monomer and polymer whereas the volume of the cage, gel and glass effects of the monomer V fm and the corresponding volume of the polymer V fp are calculated from; ( ) 0.149 2.9 10 273.15 ( ) ( ) The ratios of the critical volumes of monomer i = 1 and initiator i = I over the polymer are determined as  2.

Results and Discussion
A program to solve the above equations was developed using the object-oriented language Python 3.5 installed via the Anaconda distro.The NumPy module [22] and a customized module were used to solve for the stiff ODEs and results were ported to matplotlib whose output has a close resemblance to the well-known MATLAB format.In order to simulate the performance of the model, it was assumed that 2 tonnes of monomer were fed to the reactor and the initiator (1% wt. of monomer) was added in 1 minute at constant rate after which the addition rate was stopped.The temperature was fixed at 70˚C.The reaction progressed for the 3 hours and the reactor profiles in Figures 1-8 were obtained.As the initiator is added during the first minute at a constant rate, there's an initial increase at loading time but then the initiator and monomer amounts (see Figure 3) decrease with time when the initiator radicals combine with the monomer particle leading to chain propagation which finally results into formation of the polymer.
The reaction between the monomer and initiator is highly exothermic leading   For the polymer being produced, the molecular weight predictions of the simulated model (Mn, Mw and PDI) are given in Figures 5-8.It can be seen that the there is a slight drop in the polymer molecular weight due to a contraction in volume at low conversions.
However on the onset of the gel effect, larger polymer chains are produced which leads to a rise in the predicted molecular weights.As monomer conversion increases, the viscosity of reacting mixture increases.This causes a cage effect till the polymer chain stops growing.The polydispersity index predicted in Figure 8 is 27, a value greater than one indicating that the produced polymer is polydisperse.

Conclusion
The simulation of suspension polymerization model has been performed in Python 3.5.The prediction results were in good agreement as those reported in literature.A maximum monomer conversion of about 92.8% at a minimum batch time of about 2.2 hours was achieved at the specified conditions.The results of this work have demonstrated that Python can be employed to perform modeling and simulation studies of polymerization based processes effectively with equal success as any other programming language.

A
. Okullo et al.DOI: 10.4236/aces.2017.74029409 Advances in Chemical Engineering and Science the country.

3 )
Chain transfer to the monomer is negligible compared to other reaction steps 4) The dispersion generated in water isolates the monomer thus negligible evaporation rate 5) Each droplet behaves as a batch reactor operating in bulk 6) Termination is due to disproportionation only 7) The volume of the liquid phase is unaffected by the small amount of initiator 8) There is no polymer at the beginning of the process Using the above assumptions, the initiator (I), monomer (M) and radical species (R) balances are;

A.
Okullo et al.DOI: 10.4236/aces.2017.74029412 Advances in Chemical Engineering and Sciencewhere MW m is the molecular weight of the MMA monomer.Then the polydispersity index can be written of the monomer with respect to V, denoted as m φ and the volumetric fraction of the polymer in respect to the same, p

∨
are the critical volumes for the monomer, initiator and polymer respectively.The molecular weights for the monomer, polymer and in-The parameters used in solving the above equations are given in Table

Figure 2 .
Figure 2. Rate of change of initiator as the reaction progresses.

Figure 4 .
Figure 4. Enthalpy change with time.There is a fairly linear increase in the monomer conversion with time until the start of the gel effect as observed in Figure1.The reaction rate is then auto-accelerated thus shifting the reaction to the limiting conversion.It should be noted that maximum monomer conversion and polymer production are limited by short reaction batch times, the existence of diffusional effects in the reaction medium and by the fast decay of the initiator as observed in Figure2.

Figure 5 .
Figure 5. Variation of Mn and Mw with conversion.

Figure 6 .
Figure 6.Weight average molecular weight versus time.

Figure 7 .
Figure 7. Number average molecular weight versus time.