Improved numerical inversion methods for the recovery of bivariate distributions of polymer properties from 2D probability generating function domains

Abstract The 2D probability-generating function technique is a powerful method for modeling bivariate distributions of polymer properties. It is based on the transformation of bivariate population balance equations using 2D probability generating functions (pgf) followed by a recovery of the distributions from the transform domain by numerical inversion. A key step of this method is the inversion of the pgf transforms. Available numerical inversion methods yield excellent results for pgf transforms of distributions with independent dimensions with similar orders of magnitude, for example bivariate molecular weight distributions in copolymerization systems. However, numerical problems are found for 2D distributions in which the independent dimensions have very different ranges of values, such as the molecular weight distribution-branching distribution in branched polymers. In this work, two new 2D pgf inversion methods are developed, which regard the pgf as a complex variable. The superior accuracy of these innovative methods makes them suitable for recovering any type of bivariate distribution. This enhances the capabilities of the 2D pgf modeling technique for simulation and optimization of polymer processes. An application example of the technique in a polymeric system of industrial interest is presented.

[1]  S. Lazzari,et al.  Modeling Multiradicals in Crosslinking MMA/EGDMA Bulk Copolymerization , 2014 .

[2]  Ward Whitt,et al.  Numerical inversion of probability generating functions , 1992, Oper. Res. Lett..

[3]  P. L. Mills,et al.  Determination of polymer chain length distributions by numerical inversion of z-transforms , 1986 .

[4]  A Novel Approach to Population Balance Modeling of Reactive Polymer Modification Leading to Branching , 2013 .

[5]  Costas Kiparissides,et al.  Prediction of the Bivariate Molecular Weight - Long Chain Branching Distribution in Highly Branched Polymerization Systems Using Monte Carlo and Sectional Grid Methods , 2007 .

[6]  Michael Wulkow,et al.  Computer Aided Modeling of Polymer Reaction Engineering—The Status of Predici, I‐Simulation , 2008 .

[7]  A. Brandolin,et al.  Modeling of RAFT Polymerization using Probability Generating Functions. Detailed Prediction of Full Molecular Weight Distributions and Sensitivity Analysis , 2014 .

[8]  J. K. Cavers,et al.  On the Fast Fourier Transform Inversion of Probability Generating Functions , 1978 .

[9]  A. Brandolin,et al.  Recovery of molecular weight distributions from transformed domains. Part I. Application of pgf to mass balances describing reactions involving free radicals , 2002 .

[10]  Jitendra Kumar,et al.  Accurate and efficient solution of bivariate population balance equations using unstructured grids , 2013 .

[11]  Moon Ho Lee,et al.  2-D Laplace-Z Transformation , 2006, IEICE Trans. Fundam. Electron. Commun. Comput. Sci..

[12]  Michael Wulkow,et al.  Modeling of Functional Group Distribution in Copolymerization: A Comparison of Deterministic and Stochastic Approaches , 2014 .

[13]  Catalytic Degradation of Polystyrene: Modeling of Molecular Weight Distribution , 2011 .

[14]  Zhenping Cheng,et al.  Reversible addition‐fragmentation chain transfer polymerization of vinyl acetate under high pressure , 2015 .

[15]  M. Morbidelli,et al.  Evaluation of the Chain Length Distribution in Free-Radical Polymerization, 2. Emulsion Polymerization , 2002 .

[16]  C. Kiparissides,et al.  Prediction of the Viscoelastic Behavior of Low‐Density Polyethylene Produced in High‐Pressure Tubular Reactors , 2015 .

[17]  Piet D. Iedema,et al.  Modeling Molecular Weight and Degree of Branching Distribution of low density Polyethylene , 2000 .

[18]  Peter P. Valko,et al.  Numerical inversion of 2-D Laplace transforms applied to fractional diffusion equations , 2005 .

[19]  A. Brandolin,et al.  Mathematical Modeling of Bivariate Polymer Property Distributions Using 2D Probability Generating Functions, 1 – Numerical Inversion Methods , 2010 .

[20]  C. Kiparissides,et al.  Dynamic prediction of the bivariate molecular weight–copolymer composition distribution using sectional-grid and stochastic numerical methods , 2008 .

[21]  E. Saldívar‐Guerra,et al.  Full Molecular Weight Distribution in RAFT Polymerization. New Mechanistic Insight by Direct Integration of the Equations , 2011 .

[22]  L. Brancik.  Improved method of numerical inversion of two-dimensional Laplace transforms for dynamical systems simulation , 2002, 9th International Conference on Electronics, Circuits and Systems.

[23]  A. Brandolin,et al.  Graft Copolymers for Blend Compatibilization: Mathematical Modeling of the Grafting Process , 2012 .

[24]  J. Vlach,et al.  Numerical inversion of multidimensional Laplace transform , 1975, Proceedings of the IEEE.

[25]  Ward Whitt,et al.  Numerical Inversion of Multidimensional Laplace Transforms by the Laguerre Method , 1998, Perform. Evaluation.

[26]  John N. Daigle,et al.  Queue length distributions from probability generating functions via discrete fourier transforms , 1989 .

[27]  C. Schütte,et al.  A Hybrid Galerkin-Monte-Carlo Approach to Higher-Dimensional Population Balances in Polymerization Kinetics , 2010 .

[28]  A. Hamielec,et al.  The chemical composition component of the distribution of chain length and long chain branching for copolymerization of olefins and polyolefin chains containing terminal double-bonds , 1997 .

[29]  F. Teymour,et al.  Reconstruction of the Chain Length Distribution for Vinyl-Divinyl Copolymerization Using the Numerical Fractionation Technique , 2003 .

[30]  S. Costeux.  Statistical Modeling of Randomly Branched Polymers Produced by Combination of Several Single-Site Catalysts: Toward Optimization of Melt Properties , 2003 .

[31]  Lubomír Brancík,et al.  Technique of 3D NILT based on complex Fourier series and quotient-difference algorithm , 2010, 2010 17th IEEE International Conference on Electronics, Circuits and Systems.

[32]  A. Brandolin,et al.  Mathematical Modeling of Bivariate Distributions of Polymer Properties Using 2D Probability Generating Functions. Part II: Transformation of Population Mass Balances of Polymer Processes , 2013 .

[33]  Piet D. Iedema,et al.  Predicting Molecular Weight and Degree of Branching Distribution of Polyethylene for Mixed Systems with a Constrained Geometry Metallocene Catalyst in Semibatch and Continuous Reactors , 2003 .

[34]  F. Teymour,et al.  Mathematical model for vinyl-divinyl polymerization , 2007 .

[35]  H. Mori,et al.  RAFT Polymerization of Phenyl Vinyl Sulfide Using Trithiocarbonate Mediating Agents and Synthesis of Block Copolymers , 2015 .

[36]  Costas Kiparissides,et al.  A comprehensive model for the calculation of molecular weight–long-chain branching distribution in free-radical polymerizations , 1998 .