Computational modeling polymerization

Computer modelling provides powerful and convenient tools for the quantitative analysis of fluid dynamics and heat transfer in non-Newtonian polymer flow systems. Therefore these techniques arc routmely used in the modern polymer industry to design and develop better and more efficient process equipment and operations. The main steps in the development of a computer model for a physical process, such as the flow and deformation of polymeric materials, can be summarized as ... 

Process Modeling. The complexity of emulsion polymerization makes rehable computer models valuable. Many attempts have been made to simulate the emulsion polymerization process for different monomer systems (76—78). 

Fig. 7 gives an example of such a comparison between a number of different polymer simulations and an experiment. The data contain a variety of Monte Carlo simulations employing different models, molecular dynamics simulations, as well as experimental results for polyethylene. Within the error bars this universal analysis of the diffusion constant is independent of the chemical species, be they simple computer models or real chemical materials. Thus, on this level, the simplified models are the most suitable models for investigating polymer materials. (For polymers with side branches or more complicated monomers, the situation is not that clear cut.) It also shows that the so-called entanglement length or entanglement molecular mass Mg is the universal scaling variable which allows one to compare different polymeric melts in order to interpret their viscoelastic behavior. 

Chains with uttdesired functionality from termination by combination or disproportionation cannot be totally avoided. Tn attempts to prepare a monofunctional polymer, any termination by combination will give rise to a difunctional impurity. Similarly, when a difunctional polymer is required, termination by disproportionation will yield a monofunctional impurity. The amount of termination by radical-radical reactions can be minimized by using the lowest practical rate of initiation (and of polymerization). Computer modeling has been used as a means of predicting the sources of chain ends during polymerization and examining their dependence on reaction conditions (Section 7.5.612 0 J The main limitations on accuracy are the precision of rate constants which characterize the polymerization. 

The computer model used for this analysis is based on a plug flow tubular reactor operating under restraints of the commonly accepted kinetic mechanism for polymerization reactions ( 5 ) ... 

In order to test this computer model, we conducted experiments on thermally Initiated styrene polymerization In sealed pressure vessels. We only measured pressures and temperatures In these experiments. We conducted our tests in two phases. 

Advanced computational models are also developed to understand the formation of polymer microstructure and polymer morphology. Nonuniform compositional distribution in olefin copolymers can affect the chain solubility of highly crystalline polymers. When such compositional nonuniformity is present, hydrodynamic volume distribution measured by size exclusion chromatography does not match the exact copolymer molecular weight distribution. Therefore, it is necessary to calculate the hydrodynamic volume distribution from a copolymer kinetic model and to relate it to the copolymer molecular weight distribution. The finite molecular weight moment techniques that were developed for free radical homo- and co-polymerization processes can be used for such calculations [1,14,15]. 

Although the basic mechanisms are generally agreed on, the difficult part of the model development is to provide the model with the rate constants, physical properties and other model parameters needed for computation. For copolymerizations, there is only meager data available, particularly for cross-termination rate constants and Trommsdorff effects. In the development of our computer model, the considerable data available on relative homopolymerization rates of various monomers, relative propagation rates in copolymerization, and decomposition rates of many initiators were used. They were combined with various assumptions regarding Trommsdorff effects, cross termination constants and initiator efficiencies, to come up with a computer model flexible enough to treat quantitatively the polymerization processes of interest to us. 

As a simple computational model for the catalysis of alkene polymerization, let us consider some aspects of the general chain-propagation reaction... 

The last results described are a strong indication that any computer modeling of the activity of early transition metal catalysts for the polymerization of olefins probably requires the inclusion of the counterion in the simulations. 

Weisel, J. W., and Nagaswami, C. (1992). Computer modeling of fibrin polymerization kinetics correlated with electron microscope and turbidity observations Clot structure and assembly are kinetically controlled. Biophys. J. 63, 111-128. 

In the last few years rapid advances have been made in the field of computational crystallography, so that it is now possible to produce highly refined computer models of a wide variety of polymeric materials using X-ray diffraction data. Unfortunately, these achievements have been negated to some extent because the techniques used to collect the data for such refinement programs have not advanced at a comparable rate. In this contribution we describe a computer program which facilitates the reduction of intensity and d-spacing data obtained by the multiple film-pack method, and attempts to quantify the errors associated with such measurements. 

We approached the problem of establishing a structure of the "random coil" conformation by first establishing the limits of the present computational model for a known polymeric structure, the a-helix. Coordinates, created by program MacroModel [22] for the carbonyl groups of a-helical oligomers were used, along with published and experimental dipole transition moments, to compute the VCD and absorption spectra of the a-helical conformer. We found that VCD spectra, independent of chain length, can be calculated for octamers, and that the choice of side chain residues is immaterial for the computed spectra. Both calculated and experimental data were normalized to one residue, to permit a comparison between computed and observed spectra. 

The first discovered solid phase of fullerenes C6o represents typical molecular crystal. Later it was established that high pressure applied to solid C6o at high temperature induces polymerization of C6o [1-2]. Using the computer modeling methods allows confirming the existence of at least three different planar polymerized structures of fullerene Cgo with coordination numbers 2, 4, 6, and besides the values 4 and 6 are more probable ones. 

Michalak A, Ziegler T, The Key Steps in Olefin Polymerization Catalyzed by Late Transition Metals. In Computational Modeling of Homogeneous Catalysis, edited by F Maseras, A Lledos (Kluwer Academic Publishers, 2002)... 

In recent years, computer simulation see Solids Computer Modeling) methods have been used successfidly to examine the relationship between structural properties and transport mechanism in crystalline amorphous and polymeric materials. ... 

Chapter 10 is a modest attempt to introduce polymerization reaction engineering. 1 would hope that this subject will be interesting and useful not just to engineers but to scientists as well, because it is always informative to see how basic concepts are applied in practice. In this connection, the student and practitioner should realize that there are lessons to be learned not only in industrial versions of laboratory-scale reactions but also in how and why some processes are not used. Computer modeling of polymerization processes has not been included in this chapter because of space limitations. 

Segmental Behavior. The understanding of segmental behavior is clearly a critical factor in many of the frontiers discussed, both because synthetic capability has broadened the potential array of new materials and because this understanding is a vital key in predicting the physical behavior and uses of polymeric materials. With computer modeling, significant improvements in this area can be achieved. 

These models require information about mean velocity and the turbulence field within the stirred vessels. Computational flow models can be developed to provide such fluid dynamic information required by the reactor models. Although in principle, it is possible to solve the population balance model equations within the CFM framework, a simplified compartment-mixing model may be adequate to simulate an industrial reactor. In this approach, a CFD model is developed to establish the relationship between reactor hardware and the resulting fluid dynamics. This information is used by a relatively simple, compartment-mixing model coupled with a population balance model (Vivaldo-Lima et al., 1998). The approach is shown schematically in Fig. 9.2. Detailed polymerization kinetics can be included. Vivaldo-Lima et a/. (1998) have successfully used such an approach to predict particle size distribution (PSD) of the product polymer. Their two-compartment model was able to capture the bi-modal behavior observed in the experimental PSD data. After adequate validation, such a computational model can be used to optimize reactor configuration and operation to enhance reactor performance. 

The main objectives in modeling polymerization reactions are to compute polymerization rate and polymer properties for various reaction conditions. These two types of model outputs are not separate but they are usually very closely related. For example, an increase in reaction temperature raises polymerization rate... 

A new approach to calculating viscosities of concentrated polymer solutions has been presented. It consists of the derivation of a semi-empirical computational model containing three parameters characteristic of a particular polymer. Once these parameters have been established, the viscosity of any solution of the polymeric material in a solvent or solvent blend may be calculated. The method should be of particular interest to the coatings industry, where they often require a screening estimate of the potential viscosity-reducing power of a new solvent blend. 

---