Acessibilidade / Reportar erro

MODELING PARTICLE SIZE DISTRIBUTION IN HETEROGENEOUS POLYMERIZATION SYSTEMS USING MULTIMODAL LOGNORMAL FUNCTION

Abstract

This work evaluates the usage of the multimodal lognormal function to describe Particle Size Distributions (PSD) of emulsion and suspension polymerization processes, including continuous reactions with particle re-nucleation leading to complex multimodal PSDs. A global optimization algorithm, namely Particle Swarm Optimization (PSO), was used for parameter estimation of the proposed model, minimizing the objective function defined by the mean squared errors. Statistical evaluation of the results indicated that the multimodal lognormal function could describe distinctive features of different types of PSDs with accuracy and consistency.

Keywords:
Particle Size Distribution; Lognormal function; Parameter estimation; Particle Swarm Optimization; Heterogeneous polymerization

INTRODUCTION

The particle size distribution (PSD) is one of the most important characteristics of polymer latexes/resins, since properties such as viscosity, maximum solids content, adhesion and drying time depend on the profile of this distribution (Vale and McKenna, 2005Vale, H. M., McKenna, T. F., Modeling particle size distribution in emulsion polymerization reactors. Progress in Polymer Science, 30, 1019-1048 (2005). http://dx.doi.org/10.1016/j.progpolymsci.2005.06.006
http://dx.doi.org/10.1016/j.progpolymsci...
). Furthermore, the understanding of the dynamics of PSD is very important for process monitoring and final quality control of heterogeneous polymerization processes (Hosseini et al., 2013Hosseini, A., Milad, O., Engell, S., Mid-course control of particle size distribution in emulsion poly-merization using a hybrid model. IEEE International Conference on Control Applications, Hyderabad, India (2013). DOI:10.1109/CCA.2013.6662836
https://doi.org/10.1109/CCA.2013.6662836...
). Therefore, PSD modeling has received great attention with the proposal of models such as Population Balance Equations (PBE) (Machado et al., 2000Machado, R. A. F., Pinto, J. C., Araújo, P. H. H., Bolzan, A., Mathematical modeling of polystyrene particle size distribution produced by suspension polymerization. Brazilian Journal of Chemical Engineering, 17, 395-407 (2000). http://dx.doi.org/10.1590/S0104-66322000000400004
http://dx.doi.org/10.1590/S0104-66322000...
; Araújo et al., 2001Araújo, P. H. H., De la Cal, J. C., Asua, J. M., Pinto, J. C., Modeling particle size distribution (PSD) in emulsion copolymerization reactions in a continuous loop reactor. Macromol. Theory Simul., 10, 769-79 (2001). DOI:10.1002/1521-3919(20011001)10:8<769::AID-MATS769>3.0.CO;2-Z
https://doi.org/10.1002/1521-3919(200110...
; Kiparissides et al., 2002Kiparissides, C., Achilias, D. S., Frantzikinakis, C. E., The effect of oxygen on the kinetics and particle size distribution in vinylchloride emulsion polymerization. Ind. Eng. Chem. Res., 41, 3097-3109 (2002). DOI:10.1021/ie010928f
https://doi.org/10.1021/ie010928f...
; Immanuel et al., 2002Immanuel, C. D., Cordeiro, C. F., Sundaram, S. S., Meadows, E. S., Crowley, T. J., Doyle III, F. J., Modeling of particle size distribution in emulsion co-polymerization: Comparison with experimental data and parametric sensitivity studies. Comput. Chem. Eng., 26, 1133-1152 (2002). http://dx.doi.org/10.1016/S0098-1354(02)00031-5
http://dx.doi.org/10.1016/S0098-1354(02)...
; Coen et al., 2004Coen, E. M., Peach, S., Morrison B. R., Gilbert R. G., First principles calculation of particle formation in emulsion polymerization: pseudo-bulk systems. Polymer, 45, 595-608 (2004). http://dx.doi.org/10.1016/j.polymer.2004.03.084
http://dx.doi.org/10.1016/j.polymer.2004...
). However, according to Vale and McKenna (2005)Vale, H. M., McKenna, T. F., Modeling particle size distribution in emulsion polymerization reactors. Progress in Polymer Science, 30, 1019-1048 (2005). http://dx.doi.org/10.1016/j.progpolymsci.2005.06.006
http://dx.doi.org/10.1016/j.progpolymsci...
, even when these difficulties in coagulation modeling are overcome, the PBEs are difficult to solve if they include kinetic and/or hydrodynamic complete models. Despite recent advances in techniques for experimental determination, difficulties in determining important variables related to the quality and productivity of polymers still remain (Machado et al., 2007Machado, F., Lima, E. L., Pinto, J. C., Uma Revisão Sobre os Processos de Polimerização em Suspensão. Polímeros: Ciência e Tecnologia, 17, 166-179 (2007). (In Portuguese). http://dx.doi.org/10.1590/S0104-14282007000200016.
http://dx.doi.org/10.1590/S0104-14282007...
). In the case of PSD, these limitations are associated with the existence of polydisperse distributions and the long times often required for sample preparation and analysis.

In this way, simpler and thus less computationally expensive alternatives for the description of complex PSDs are very useful for different types of applications as, for instance, the development of soft-sensors for on-line PSD monitoring (Clementi et al., 2013Clementi, L. A., Vega, J. R., Gugliotta, L. M., Soft sensor for on-line monitoring of particle sizes in miniemulsion polymerization. Mecánica Computacional, 32 (2013). http://www.cimec.org.ar/ojs/index.php/mc/article/viewFile/4494/4424
http://www.cimec.org.ar/ojs/index.php/mc...
), the implementation of optimization and control strategies or the simulation of heterogeneous polymerize- tion reactions using CFD tools (Elgebrandt et al., 2006Elgebrandt, R. C., Fletcher, D. F., Gomes, V. G., Romagnoli, J. A., A framework for modeling particle size effects in emulsion polymerization systems using computational fluid dynamics linked to a detailed population balance model. Computer Aided Chemical Engineering, 21 (2006). DOI:10.1016/S1570-7946(06)80103-3
https://doi.org/10.1016/S1570-7946(06)80...
; Pohn et al., 2013Pohn, J., Cunningham, M., McKenna, T. F. L., Scale-up of emulsion polymerization reactors Part I - development of a model framework. Macromolecular Reaction Engineering, 7, (2013). DOI:10.1002/mren.201300010
https://doi.org/10.1002/mren.201300010...
). These alternatives require a model describing the PSD.

According to Limbert et al. (2001)Limpert, E., Stahel, W., Abbt, M., Lognormal distribution across the sciences: Keys and clues. BioScience, 51, 341-352 (2001). DOI:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2
https://doi.org/10.1641/0006-3568(2001)0...
the multimodal lognormal function has proved to be efficient and flexible in describing wide or narrow, unimodal and multimodal distributions in many particulate systems applications.

This diversity of applications can be observed in works such as Zhao et al. (2003)Zhao, B., Yang, Z., Johnston, M. V., Wang, H., Wexler, A. S., Balthasar, M., Kraft, M., Measurement and numerical simulation of soot particle size distribution functions in a laminar premixed ethylene-oxygen-argon flame. Combustion and Flame, 133, 173-188 (2003). http://dx.doi.org/10.1016/S0010-2180(02)00574-6
http://dx.doi.org/10.1016/S0010-2180(02)...
who investigated the dynamics of soot particle size distributions. The numerical simulation using a kinetic model proposed previously and a stochastic approach to solve aerosol dynamics equations showed a bimodal soot particle size distribution function, the consequence of the interplay between particle-particle coagulation and particle nucleation. In this same respect, Johnsson et al. (2010)Johnsson, J., Bladh, H., Bengtsson, P. E., On the influence of bimodal size distributions in particle sizing using laser-induced incandescence. Applied Physics, B, Lasers e Optics, 99, 817-823 (2010). DOI:10.1007/s00340-010-3986-x
https://doi.org/10.1007/s00340-010-3986-...
claimed that soot particle distributions generally consist of a range of different particle sizes, and the lognormal size distribution has been shown to be a good approximation for various flame conditions; Hwang and Choi (2006)Hwang, S., Choi, S., Use of a lognormal distribution model for estimating soil water retention curves from particle-size distribution data. Journal of Hydrology, 323, 325-334 (2006). DOI:10.1016/j.jhydrol.2005.09.005
https://doi.org/10.1016/j.jhydrol.2005.0...
investigated the predictive potential of the unimodal lognormal distribution model for estimating the water retention curves, and this model was evaluated for a broader range of soil; Yuan et al. (2011)Yuan, Y., Yi, H., Shuai, Y., Liu, B., Tan, H., Inverse problem for aerosol particle size distribution using SPSO associated with multi-lognormal distribution model. Atmospheric Environment, 45, 4892-4897 (2011). http://dx.doi.org/10.1016/j.atmosenv.2011.06.010
http://dx.doi.org/10.1016/j.atmosenv.201...
described the characteristics of atmospheric aerosols using a multi-lognormal distribution model with parameters estimated by Particle Swarm Optimization. The model was validated and analyzed by comparing with the measured optical properties of several bands; Pujol and Pinto (2011)Pujol, J. C. F., Pinto, J. M. A., A neural network approach to fatigue life prediction. International Journal of Fatigue, 33, 313-322 (2011). http://dx.doi.org/10.1016/j.ijfatigue.2010.09.003
http://dx.doi.org/10.1016/j.ijfatigue.20...
approached fatigue life prediction under step-stress conditions by comparing a standard approach based on the lognormal distribution function implemented and fit to experimental data, with an approach for fatigue life prediction from neural networks. Both models were optimized by differential evolution associated to the Newton-Raphson-Method, using a maximum likelihood estimator.

Given the above, this paper aims to evaluate the use of a multimodal lognormal function to describe the behavior of particle size distributions in heterogeneous polymerization systems, using a Particle Swarm Optimization algorithm for parameter estimation.

MATERIALS AND METHODS

Experimental Particle Size Distributions

In order to evaluate whether the proposed multimodal lognormal function was able to describe real PSDs, experimental PSDs of emulsion and suspension polymerization reactions were used for the parameter estimation of the multimodal lognormal function. The choice of the distributions used in this work was based on their different characteristics and on the existence of multimodalities.

The reaction techniques, monomer system, reactor type and operation mode of the experimental PSD used in this work are summarized in Table 1. Four PSDs (PSD01 to PSD04), among which there were unimodal and multimodal distributions, were obtained from the work of Araújo (1999)Araújo, P. H. H., Particle size distribution in heterogeneous polymerization systems. D.Sc. Thesis, COPPE/Universidade Federal do Rio de Janeiro. Rio de Janeiro, Brazil (1999).. In short, this author investigated the effects of operational conditions (temperature, initiator and surfactant concentrations) on the PSD in an emulsion copolymerization of vinyl acetate/Veova10 in a continuous loop reactor with high solids content. The reactor pre-feeding with an aqueous (surfactant, protective colloid and Na2S2O5) and an organic phase (monomer mixture with vinyl acetate to Veova 10 weight ratio of 75/25) caused a fast increase of the particle number at the reaction beginning. As a consequence, the distributions usually show a right skewness. The higher solubility of the oligomeric radical in the aqueous phase, due to the vinyl acetate, affects the main "locus" of the polymerization, promoting re-nucleations.

Table 1
Heterogeneous polymerization systems data.

Another PSD (PSD05) was obtained from the work of Zubitur and Asua (2001)Zubitur, M., Asua, J. M., Factors affecting kinetics and coagulum formation during the emulsion co-polymerization of styrene/butyl acrylate. Polymer, 42, 5979-5985 (2001). http://dx.doi.org/10.1016/S0032-3861(01)00039-8
http://dx.doi.org/10.1016/S0032-3861(01)...
which investigated the factors affecting kinetics and coagulation during the emulsion polymerization of styrene and butyl acrylate, in a semi-batch stirred tank reactor. The PSD showed some bimodality, a peak of small particles, probably resulting from secondary nucleation in the presence of seed particles (shoulder of bigger particles).

Mallikarjunan et al. (2010)Mallikarjunan, V., Pushpavanam, S., Immanuel, C. D., Parameter estimation strategies in batch emulsion polymerization. Chemical Engineering Science, 65, 4967-4982 (2010). http://dx.doi.org/10.1016/j.ces.2010.05.027
http://dx.doi.org/10.1016/j.ces.2010.05....
studied the emulsion copolymerization of vinyl acetate and butyl acrylate, conducted in a semi-batch stirred tank reactor, and another PSD (PSD06) was utilized from the work of these authors. They report a bimodal PSD determined with CHDF, since the system was fed with surfactant in excess.

In addition to the PSDs of emulsion polymerization reactions, PSDs of suspension polymerization reactions of methyl methacrylate were also evaluated. These PSDs are available in the work of Jahanzad et al. (2004)Jahanzad. F., Evolution of particle size distribution in suspension polymerization reactions. Submitted in partial fulfillment of the requirements for the award of the Doctor of Philosophy of the Loughborough University, Leicestershire, UK (2004). https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
https://dspace.lboro.ac.uk/dspace-jspui/...
, and they were included in this work with the objective of investigating the adaptability of the multimodal lognormal function to other heterogeneous polymerization techniques. In this paper, the authors carried out extensive research on the influence of process variables (temperature, inhibition, stabilizer and initiator concentrations, agitation) on the PSD in a batch stirred tank reactor. The distributions obtained by Jahanzad et al. (2004)Jahanzad. F., Evolution of particle size distribution in suspension polymerization reactions. Submitted in partial fulfillment of the requirements for the award of the Doctor of Philosophy of the Loughborough University, Leicestershire, UK (2004). https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
https://dspace.lboro.ac.uk/dspace-jspui/...
are basically unimodal. In general, some factors promoted the growth of the mean particle size and the broadening of the distributions, such as the increase of the reaction temperature, decrease of the stirring rate, decrease of the stabilizer concentration and absence of inhibitor.

In order to normalize the distributions range, a numerical integration (trapezoid method) was used to calculate the area under each PSD. The different types of distributions (number/weight/volume density distributions) of these reference works were normalized by the division of the PSD by the calculated area. In this way, the sum of the relative frequencies was equal to one for each PSD.

Multimodal Lognormal Function

The lognormal distribution is a model with two parameters, whose density function is expressed by Equation (1) (Crow and Shimizu, 1988Crow, E. L., Shimizu, K., Log-normal Distributions: Theory and Application. New York, Dekker (1988).):

(1)

where x is the independent variable, µg σg represent the geometric mean and geometric standard deviation of the distribution, respectively.

For frequency distributions, xk values are weighted by their relative frequencies fk with k = 1, 2 ..., n. In this case, the geometric mean and geometric standard deviation can be defined by Equations (2) and (3), respectively (Crow and Shimizu, 1988Crow, E. L., Shimizu, K., Log-normal Distributions: Theory and Application. New York, Dekker (1988).):

(2)

(3)

where is the absolute frequency.

However, in real situations, these distributions can be composed of more than one mode, being bimodal, trimodal or even n-modal. In such cases, the function can be described by a weighted sum of lognormal functions (Zender, 2010). Equation (4) describes the multimodal lognormal function:

(4)

where MLF(x, σgi, μgi, wi) is the Multimodal Lognormal Function; f(x, σgi, μgi) is the lognormal function mode i; n number of modes and wi is the fraction in mode i assuming .

Taking into account the number of modes of the experimental distributions used in this work, the Bimodal Lognormal Function (BLF) and Trimodal Lognormal Function (TLF) were chosen for modeling experimental PSDs. Equation (5) provides a representation of the function used, with the particle diameter d as independent variable:

(5)

Therefore, there were always five or eight parameters to be estimated for BLF and TLF respectively, since there were three parameters to be estimated for each mode (In μgi, In σgi and wi).

Parameter Estimation

The parameter estimation procedure of the MLFs was supported by Particle Swarm Optimization (Kennnedy and Eberhart, 1995). This algorithm uses its intelligent strategy to randomly minimize the objective function defined by the sum squared errors [SSE], according to Equation (6):

(6)

where fexpk is the experimental frequency value in k; MLF(dk) is the frequency calculated by the multimodal lognormal function in k and n is the number of experimental points.

In short, according to Yang (2014)Yang, X. S., Nature-Inspired Optimization Algorithms. Chapter 7, Particle Swarm Optimization, Elsevier Inc. (2014). DOI:10.1016/B978-0-12-416743-8.00007-5
https://doi.org/10.1016/B978-0-12-416743...
, the PSO strategy considers that a particle i in an interaction k, moves through the search space with two attributes: position (X) and velocity (V) vectors, represented by equations (7) and (8), respectively:

(7)

(8)

where, Pik is the best previous position along the dth dimension of particle i in iteration kth, while Pgk is the best previous position among all the particles along the dth dimension in iteration kth.

According to Blum and Merkle (2008)Blum, C., Merkle, D., Swarm Intelligence - Introduction and Applications. Springer-Verlag Berlin Heidelberg (2008). http://www.springer.com/computer/ai/book/978-3-540-74088-9
http://www.springer.com/computer/ai/book...
the velocity vector is determined by three terms: the "momentum", the "cognitive", and the "social" term. The "momentum" term w.vi represents the previous velocity which is used to carry the particle in the direction it has travelled so far; the "cognitive" term c1r1(Plk -Xik), represents the tendency of the particle to return to the best position it has visited so far; the "social" term c2r2(Pgk -Xik), represents the tendency of the particle to be attracted towards the position of the best position found by the entire swarm. Parameters c1 and c2 are acceleration coefficients, while r1 and r2 are two independent random numbers uniformly distributed in [0,1].

The procedure proposed by Jiao et al. (2008)Jiao, B., Lian, Z., Gu, X., A dynamic inertia weight particle swarm optimization algorithm. Chaos Solitons & Fractals, 37, 698-705 (2008). http://dx.doi.org/10.1016/j.chaos.2006.09.063
http://dx.doi.org/10.1016/j.chaos.2006.0...
, which considers the implementation of the dynamic inertial weight (Equation (9)) as a way to control the particle momentum, weighing the contribution of the previous velocity and thus, basically controlling how much memory of the previous flight direction will influence the new velocity, was used in this work. In this way, algorithm parameters follow the references of Jiao et al. (2008)Jiao, B., Lian, Z., Gu, X., A dynamic inertia weight particle swarm optimization algorithm. Chaos Solitons & Fractals, 37, 698-705 (2008). http://dx.doi.org/10.1016/j.chaos.2006.09.063
http://dx.doi.org/10.1016/j.chaos.2006.0...
and Van den Bergh (2002)van den Berg, F., An Analysis of Particle Swarm Optimizers. PhD Thesis, Department of Computer Science, University of Pretoria, Pretoria, South Africa (2002). who investigated the convergence criteria of the algorithm.

(9)

As this is a strategy that performs a random search, it is necessary to set upper and lower limits of the parameters to be adjusted. So, the limits of the parameters lnµg were defined from the limits of the diameter variable on a logarithm scale; in turn the limits of parameters lnσg were defined based on previous works on modeling of particle size distribution using the lognormal function (Limpert et al., 2001Limpert, E., Stahel, W., Abbt, M., Lognormal distribution across the sciences: Keys and clues. BioScience, 51, 341-352 (2001). DOI:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2
https://doi.org/10.1641/0006-3568(2001)0...
). Table 2 shows the employed search ranges, and the initial parameter estimates generated randomly from a uniform distribution, on the PSO.

Table 2
Search range for the MLFs parameters.

Two stopping criteria of the algorithm were used: the maximum number of iterations (8000) or evaluation of Equation (10):

(10)

where Yg and Yz are, respectively, the best result found and the best z last results.

The parameter estimation strategy was implemented in Fortran 90 language and parameter estimations were carried out in triplicate. Reported optimal parameters are arithmetic mean values of the three estimations. Standard deviations of each parameter are also reported, representing the algorithm variability due to the random feature of the search. Therefore, it is important to point out that standard deviation values have no relationship with the confidence interval of the optimal parameter.

Statistical interpretation of the results was performed through the correlation coefficient [r] between experimental and calculated frequency distributions. In addition, the Kolmogorov-Smirnov test was used to compare experimental frequency distributions with estimated frequency distributions. The Kolmogorov-Smirnov test and Chi-square test are widely used in the scientific literature and the first one is more efficient for continuous frequency distributions (Evren and Tuna, 2012Evren, A., Tuna, E., On some properties of goodness of fit measures based on statistical entropy. International Journal of Research and Reviews in Applied Sciences, 13 (2012). http://www.arpapress.com/Volumes/Vol13Issue1/IJRRAS_13_1_22.pdf
http://www.arpapress.com/Volumes/Vol13Is...
).

RESULTS AND DISCUSSION

The first set of PSDs, PSD01 to PSD04, corresponds to the evolution during a continuous vinyl acetate/Veova 10 emulsion copolymerization reaction conducted in a tubular loop reactor. Due to successive nucleations this reaction presents broad and often multimodal PSDs (Rawlings and Ray, 1988Rawlings, J. B., Ray, W. H., The modeling of batch and continuous emulsion polymerization reactors: I. Model formulation and sensitivity to parameters. Polymer Engineering and Science, 28, 237-256 (1988). http://onlinelibrary.wiley.com/doi/10.1002/pen.760280502/abstract
http://onlinelibrary.wiley.com/doi/10.10...
). This feature generates consequences for the determination of the search interval of the parameter geometric mean. Besides, the positively skewed modes affect strongly the geometric standard deviation, minimizing it in this case.

Figure 1 compares the experimental PSDs with and PSDs predicted using a BLF. It may be observed that the BLF could fit nicely the broad monomodal and positively skewed distribution PSD01, as well as different bimodal distributions with different mode locations and relative fractions (PSD02 and PSD04). As expected, since the proposed model provides the generation of two modes, it was not able to represent all the three modes of PSD03. Nevertheless, even in this case the BLF was able to predict correctly the location, width and relative intensity of the two bigger particle populations, while the less significant mode of the distribution was ignored.

Figure 1
Experimental (Araújo, 1999Araújo, P. H. H., Particle size distribution in heterogeneous polymerization systems. D.Sc. Thesis, COPPE/Universidade Federal do Rio de Janeiro. Rio de Janeiro, Brazil (1999).) and predicted PSDs of a continuous vinyl acetate Veova 10 emulsion copolymerization reaction in a loop reactor using a bimodal lognormal function (BLF).

When a trimodal lognormal function (TLF) is used to represent the evolution of the PSD during this continuous emulsion polymerization reaction, Figure 2, the TLF fits perfectly the experimental distributions, no matter if mono-, bi- or trimodal, showing better agreement than with BLF for all PSDs. This result is due to the higher number of parameters that ensures the representation of trimodal distributions and favors data fitting when the distribution is simpler and, thus, with overlapping modes, with very close parameters, favoring the parameter estimation procedure. The higher number of parameters, though, can lead to a high correlation among the TLF parameters.

Figure 2
Experimental (Araújo, 1999Araújo, P. H. H., Particle size distribution in heterogeneous polymerization systems. D.Sc. Thesis, COPPE/Universidade Federal do Rio de Janeiro. Rio de Janeiro, Brazil (1999).) and predicted PSDs of a continuous vinyl acetate Veova 10 emulsion copolymerization reaction in a loop reactor using a trimodal lognormal function (TLF).

Figure 3 compares the predictions of a TLF with the final experimental PSDs of semi-batch styrene/butyl acrylate (PSD05, Zubitur and Asua, 2001Zubitur, M., Asua, J. M., Factors affecting kinetics and coagulum formation during the emulsion co-polymerization of styrene/butyl acrylate. Polymer, 42, 5979-5985 (2001). http://dx.doi.org/10.1016/S0032-3861(01)00039-8
http://dx.doi.org/10.1016/S0032-3861(01)...
) and vinyl acetate/butyl acrylate (PSD06, Mallikarjunan et al., 2010Mallikarjunan, V., Pushpavanam, S., Immanuel, C. D., Parameter estimation strategies in batch emulsion polymerization. Chemical Engineering Science, 65, 4967-4982 (2010). http://dx.doi.org/10.1016/j.ces.2010.05.027
http://dx.doi.org/10.1016/j.ces.2010.05....
) emulsion copolymerization reactions and the TLF was able to predict correctly both bimodal PSDs, shoulder of bigger particles and two separate populations.

Figure 3
Experimental (Zubitur and Asua, 2001Zubitur, M., Asua, J. M., Factors affecting kinetics and coagulum formation during the emulsion co-polymerization of styrene/butyl acrylate. Polymer, 42, 5979-5985 (2001). http://dx.doi.org/10.1016/S0032-3861(01)00039-8
http://dx.doi.org/10.1016/S0032-3861(01)...
and Mallikarjunan et al., 2010Mallikarjunan, V., Pushpavanam, S., Immanuel, C. D., Parameter estimation strategies in batch emulsion polymerization. Chemical Engineering Science, 65, 4967-4982 (2010). http://dx.doi.org/10.1016/j.ces.2010.05.027
http://dx.doi.org/10.1016/j.ces.2010.05....
) and predicted PSDs of semi-batch emulsion copolymerization reactions in a stirred tank reactor using a trimodal lognormal function (TLF).

In order to verify the adaptability of the TLF to other heterogeneous polymerization techniques, this function was used to predict the PSDs of methyl methacrylate of suspension polymerization reactions (PSD07 and PSD08, Jahanzad et al., 2004Jahanzad. F., Evolution of particle size distribution in suspension polymerization reactions. Submitted in partial fulfillment of the requirements for the award of the Doctor of Philosophy of the Loughborough University, Leicestershire, UK (2004). https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
https://dspace.lboro.ac.uk/dspace-jspui/...
). Again the TLF predicted correctly both mono- as well as multimodal PSDs, despite the different order of magnitude of the particle diameters and forms of the PSDs, when compared to those of the previously evaluated emulsion polymerization reactions.

Table 3 and 4 show the average and standard deviation values for the parameters set for BLF and TLF, respectively. These tables also show the results of the objective function value (SSE) for each case. The error in the parameter estimation was quantified from the results obtained in minimization triplicates.

Table 3
Estimated parameters for BLF model.
Table 4
Estimated parameters for TLF model.

It can be verified that all standard deviation and SSE values are low, confirming the ability of the algorithm to reproduce consistently the PSDs, thus, corroborating PSO robustness. Finally, a comparison between standard deviations presented in Table 3 and 4 shows that they are higher when using TLF, suggesting a greater correlation between the parameters of TLF, since in this function there exists a greater number of possible parameter combinations that result from the minimization of the function.

It can be pointed out that scale differences in the particle diameters between emulsion polymerization [nm] and suspension polymerization [µm] did not affect the proposed model for PSDs, and also that the TFL model was very sensitive to small variations in the distribution like the presence of a very small population of small particles in PSD08 in Figure 4.

Figure 4
Experimental (Jahanzad, 2004Jahanzad. F., Evolution of particle size distribution in suspension polymerization reactions. Submitted in partial fulfillment of the requirements for the award of the Doctor of Philosophy of the Loughborough University, Leicestershire, UK (2004). https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
https://dspace.lboro.ac.uk/dspace-jspui/...
) and predicted PSDs of batch suspension methyl methacrylate polymerization reaction in a stirred tank reactor using a trimodal lognormal function (TLF).

Evaluation of the correlation coefficient (r) was consistent with the results shown in Figures 1 to Figure 4, with r ≥ 0.98 in all cases using BLF and r ≥ 0.99 in all cases using TLF, evidencing a strong correlation between experimental and predicted frequency distributions. The Kolmogorov-Smirnov test and the established critical values resulted in a dcritical value equal to 0.210 for a significance level (α) of 5% and for a number of observed points (n) equal to 40, used in this work. As the largest value dmáx evaluation for PSDs predicted is equal to 0.0147, the hypothesis that the data refer to the same distribution should be accepted, increasing the confidence in the performed parameter estimation.

Finally, it can be pointed out that all variability in the evaluated PSDs could be perfectly assimilated by TLF, and, except for the smaller third population of large particles of PSD03, also by the BLF. The proposal did not present inconsistencies even for the different particle diameter ranges commonly provided by different polymerizations techniques. This difficulty was overcome by using normalized relative frequencies.

CONCLUSIONS

A multimodal lognormal function was proposed to describe the particle size distributions in heterogeneous polymerization systems, with parameters estimated by the global optimization algorithm Particle Swarm Optimization. Distributions with different characteristics, unimodal or multimodal, wide or narrow, but all originating from polymerization reactions in emulsion or suspension were modeled.

High correlations between experimental and predicted frequency distributions using bimodal and trimodal lognormal functions and the Kolmogorov-Smirnov test were obtained for all particle size distributions, confirming that the proposed function can describe successfully particle size distributions with different characteristics.

Even when different particle diameter ranges were used, the proposed model could perfectly represent the distribution profile. Furthermore, the accuracy of the parameter estimation was evidenced by low standard deviation values, attesting to the ability of Particle Swarm Optimization to reproduce consistent results.

Comparing the trimodal lognormal function with the bimodal lognormal function, the use of the former is justified only when the investigated particle size distributions are very polydisperse with three relevant modes.

NOMENCLATURE

BLF Bimodal Lognormal Function CFD Computational Fluid Dynamics CHDF Capillary Hydrodynamic Fractionation d Particle diameter exp experimental f relative frequencies i number of mode MLF Multimodal Lognormal Function N absolute frequency n number of experimental points. Pg previous position among all the particles Pl previous position along the dth dimension of particle PSD Particle Size Distributions PSO Particle Swarm Optimization r correlation coefficient SSE sum squared errors TLF Trimodal Lognormal Function u control parameter V velocity vector w Inertial weight wi fraction of mode i wo Initial inertial weight x independent variable X position vector Yg best result found Yz best z last results Greek Symbols µ g geometric mean σ g geometric standard deviation
  • This is an extended version of the work presented at the 20th Brazilian Congress of Chemical Engineering, COBEQ-2014, Florianópolis, Brazil.

ACKNOWLEDGMENTS

The authors thank the financial support from CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico.

REFERENCES

  • Araújo, P. H. H., Particle size distribution in heterogeneous polymerization systems. D.Sc. Thesis, COPPE/Universidade Federal do Rio de Janeiro. Rio de Janeiro, Brazil (1999).
  • Araújo, P. H. H., De la Cal, J. C., Asua, J. M., Pinto, J. C., Modeling particle size distribution (PSD) in emulsion copolymerization reactions in a continuous loop reactor. Macromol. Theory Simul., 10, 769-79 (2001). DOI:10.1002/1521-3919(20011001)10:8<769::AID-MATS769>3.0.CO;2-Z
    » https://doi.org/10.1002/1521-3919(20011001)10:8<769::AID-MATS769>3.0.CO;2-Z
  • Blum, C., Merkle, D., Swarm Intelligence - Introduction and Applications. Springer-Verlag Berlin Heidelberg (2008). http://www.springer.com/computer/ai/book/978-3-540-74088-9
    » http://www.springer.com/computer/ai/book/978-3-540-74088-9
  • Clementi, L. A., Vega, J. R., Gugliotta, L. M., Soft sensor for on-line monitoring of particle sizes in miniemulsion polymerization. Mecánica Computacional, 32 (2013). http://www.cimec.org.ar/ojs/index.php/mc/article/viewFile/4494/4424
    » http://www.cimec.org.ar/ojs/index.php/mc/article/viewFile/4494/4424
  • Coen, E. M., Peach, S., Morrison B. R., Gilbert R. G., First principles calculation of particle formation in emulsion polymerization: pseudo-bulk systems. Polymer, 45, 595-608 (2004). http://dx.doi.org/10.1016/j.polymer.2004.03.084
    » http://dx.doi.org/10.1016/j.polymer.2004.03.084
  • Crow, E. L., Shimizu, K., Log-normal Distributions: Theory and Application. New York, Dekker (1988).
  • Elgebrandt, R. C., Fletcher, D. F., Gomes, V. G., Romagnoli, J. A., A framework for modeling particle size effects in emulsion polymerization systems using computational fluid dynamics linked to a detailed population balance model. Computer Aided Chemical Engineering, 21 (2006). DOI:10.1016/S1570-7946(06)80103-3
    » https://doi.org/10.1016/S1570-7946(06)80103-3
  • Evren, A., Tuna, E., On some properties of goodness of fit measures based on statistical entropy. International Journal of Research and Reviews in Applied Sciences, 13 (2012). http://www.arpapress.com/Volumes/Vol13Issue1/IJRRAS_13_1_22.pdf
    » http://www.arpapress.com/Volumes/Vol13Issue1/IJRRAS_13_1_22.pdf
  • Hosseini, A., Milad, O., Engell, S., Mid-course control of particle size distribution in emulsion poly-merization using a hybrid model. IEEE International Conference on Control Applications, Hyderabad, India (2013). DOI:10.1109/CCA.2013.6662836
    » https://doi.org/10.1109/CCA.2013.6662836
  • Hwang, S., Choi, S., Use of a lognormal distribution model for estimating soil water retention curves from particle-size distribution data. Journal of Hydrology, 323, 325-334 (2006). DOI:10.1016/j.jhydrol.2005.09.005
    » https://doi.org/10.1016/j.jhydrol.2005.09.005
  • Immanuel, C. D., Cordeiro, C. F., Sundaram, S. S., Meadows, E. S., Crowley, T. J., Doyle III, F. J., Modeling of particle size distribution in emulsion co-polymerization: Comparison with experimental data and parametric sensitivity studies. Comput. Chem. Eng., 26, 1133-1152 (2002). http://dx.doi.org/10.1016/S0098-1354(02)00031-5
    » http://dx.doi.org/10.1016/S0098-1354(02)00031-5
  • Jahanzad. F., Evolution of particle size distribution in suspension polymerization reactions. Submitted in partial fulfillment of the requirements for the award of the Doctor of Philosophy of the Loughborough University, Leicestershire, UK (2004). https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
    » https://dspace.lboro.ac.uk/dspace-jspui/handle/2134/10300
  • Jiao, B., Lian, Z., Gu, X., A dynamic inertia weight particle swarm optimization algorithm. Chaos Solitons & Fractals, 37, 698-705 (2008). http://dx.doi.org/10.1016/j.chaos.2006.09.063
    » http://dx.doi.org/10.1016/j.chaos.2006.09.063
  • Johnsson, J., Bladh, H., Bengtsson, P. E., On the influence of bimodal size distributions in particle sizing using laser-induced incandescence. Applied Physics, B, Lasers e Optics, 99, 817-823 (2010). DOI:10.1007/s00340-010-3986-x
    » https://doi.org/10.1007/s00340-010-3986-x
  • Kennedy, J., Eberhart, R. C., Particle swarm optimization. IEEE International Conference on Neural Networks - Conference Proceedings, 4, 1942-1948 (1995). DOI:10.1109/ICNN.1995.488968
    » https://doi.org/10.1109/ICNN.1995.488968
  • Kiparissides, C., Achilias, D. S., Frantzikinakis, C. E., The effect of oxygen on the kinetics and particle size distribution in vinylchloride emulsion polymerization. Ind. Eng. Chem. Res., 41, 3097-3109 (2002). DOI:10.1021/ie010928f
    » https://doi.org/10.1021/ie010928f
  • Limpert, E., Stahel, W., Abbt, M., Lognormal distribution across the sciences: Keys and clues. BioScience, 51, 341-352 (2001). DOI:10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2
    » https://doi.org/10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2
  • Machado, R. A. F., Pinto, J. C., Araújo, P. H. H., Bolzan, A., Mathematical modeling of polystyrene particle size distribution produced by suspension polymerization. Brazilian Journal of Chemical Engineering, 17, 395-407 (2000). http://dx.doi.org/10.1590/S0104-66322000000400004
    » http://dx.doi.org/10.1590/S0104-66322000000400004
  • Machado, F., Lima, E. L., Pinto, J. C., Uma Revisão Sobre os Processos de Polimerização em Suspensão. Polímeros: Ciência e Tecnologia, 17, 166-179 (2007). (In Portuguese). http://dx.doi.org/10.1590/S0104-14282007000200016
    » http://dx.doi.org/10.1590/S0104-14282007000200016
  • Mallikarjunan, V., Pushpavanam, S., Immanuel, C. D., Parameter estimation strategies in batch emulsion polymerization. Chemical Engineering Science, 65, 4967-4982 (2010). http://dx.doi.org/10.1016/j.ces.2010.05.027
    » http://dx.doi.org/10.1016/j.ces.2010.05.027
  • Pohn, J., Cunningham, M., McKenna, T. F. L., Scale-up of emulsion polymerization reactors Part I - development of a model framework. Macromolecular Reaction Engineering, 7, (2013). DOI:10.1002/mren.201300010
    » https://doi.org/10.1002/mren.201300010
  • Pujol, J. C. F., Pinto, J. M. A., A neural network approach to fatigue life prediction. International Journal of Fatigue, 33, 313-322 (2011). http://dx.doi.org/10.1016/j.ijfatigue.2010.09.003
    » http://dx.doi.org/10.1016/j.ijfatigue.2010.09.003
  • Rawlings, J. B., Ray, W. H., The modeling of batch and continuous emulsion polymerization reactors: I. Model formulation and sensitivity to parameters. Polymer Engineering and Science, 28, 237-256 (1988). http://onlinelibrary.wiley.com/doi/10.1002/pen.760280502/abstract
    » http://onlinelibrary.wiley.com/doi/10.1002/pen.760280502/abstract
  • Vale, H. M., McKenna, T. F., Modeling particle size distribution in emulsion polymerization reactors. Progress in Polymer Science, 30, 1019-1048 (2005). http://dx.doi.org/10.1016/j.progpolymsci.2005.06.006
    » http://dx.doi.org/10.1016/j.progpolymsci.2005.06.006
  • van den Berg, F., An Analysis of Particle Swarm Optimizers. PhD Thesis, Department of Computer Science, University of Pretoria, Pretoria, South Africa (2002).
  • Yang, X. S., Nature-Inspired Optimization Algorithms. Chapter 7, Particle Swarm Optimization, Elsevier Inc. (2014). DOI:10.1016/B978-0-12-416743-8.00007-5
    » https://doi.org/10.1016/B978-0-12-416743-8.00007-5
  • Yuan, Y., Yi, H., Shuai, Y., Liu, B., Tan, H., Inverse problem for aerosol particle size distribution using SPSO associated with multi-lognormal distribution model. Atmospheric Environment, 45, 4892-4897 (2011). http://dx.doi.org/10.1016/j.atmosenv.2011.06.010
    » http://dx.doi.org/10.1016/j.atmosenv.2011.06.010
  • Zhao, B., Yang, Z., Johnston, M. V., Wang, H., Wexler, A. S., Balthasar, M., Kraft, M., Measurement and numerical simulation of soot particle size distribution functions in a laminar premixed ethylene-oxygen-argon flame. Combustion and Flame, 133, 173-188 (2003). http://dx.doi.org/10.1016/S0010-2180(02)00574-6
    » http://dx.doi.org/10.1016/S0010-2180(02)00574-6
  • Zubitur, M., Asua, J. M., Factors affecting kinetics and coagulum formation during the emulsion co-polymerization of styrene/butyl acrylate. Polymer, 42, 5979-5985 (2001). http://dx.doi.org/10.1016/S0032-3861(01)00039-8
    » http://dx.doi.org/10.1016/S0032-3861(01)00039-8

Publication Dates

  • Publication in this collection
    Jul-Sep 2016

History

  • Received
    25 Feb 2015
  • Reviewed
    19 May 2015
  • Accepted
    27 May 2015
Brazilian Society of Chemical Engineering Rua Líbero Badaró, 152 , 11. and., 01008-903 São Paulo SP Brazil, Tel.: +55 11 3107-8747, Fax.: +55 11 3104-4649, Fax: +55 11 3104-4649 - São Paulo - SP - Brazil
E-mail: rgiudici@usp.br