Macroscopic Mechanistic Modeling and Optimization of a Self-Initiated High-Temperature Polymerization Reactor

Size: px
Start display at page:

Download "Macroscopic Mechanistic Modeling and Optimization of a Self-Initiated High-Temperature Polymerization Reactor"

Transcription

1 211 American Conol Conference on O'Farrell Seet, San Francisco, CA, USA June 29 - July 1, 211 Macroscopic Mechanistic Moeling an Optimization of a Self-Initiate High-Temperature Polymerization Reactor Thomas Rier, Sriraj Srinivasan, Masou Soroush, George A. Kalfas, Michael C. Gray, an Anrew M. Rappe Absact This paper presents a macroscopic mechanistic mathematical moeling an optimization stuy of a batch polymerization reactor in which self-initiate free-raical polymerization of n-butyl acrylate at 14 an 16 C takes place. The moel is obtaine using a comprehensive free-raical polymerization reaction mechanism. The rate constant of the monomer self-initiation is estimate from monomer conversion measurements. The moel is valiate using a ifferent set of conversion measurements. The valiation results show that the macroscopic mechanistic moel is accurate enough for optimization of the self-initiate polymerization reactor to prouce high quality acrylic resins. The moel is then use to calculate an optimal batch-reactor temperature profile that yiels an en-batch polymer prouct with esire properties (conversion an number-average molecular weight). I. INTRODUCTION Acrylic resins are wiely use in automobile paint an coating formulations as weather-proofing agents an constitute a multi-billion ollar global market that is continuously growing [1], [2]. In the last several ecaes, because of singent environmental regulations in North America, the volatile organic content in paints an coatings has been reuce to less than 3 g/l [3]. Acrylic resins with high molecular weight an high solvent content were replace with low molecular weight, low solvent content resins. Hightemperature (greater than 1 C) free-raical polymerization processes were evelope to prouce these resins in place of conventional low-temperature polymerization processes. Initiators such as organic peroxies an azoniiles have been commonly use to initiate high-temperature polymerization of alkyl acrylates, but their high costs an tenency to prouce coloration when they are left unreacte in the en-prouct have mae them less atactive to use [4], [5]. Spontaneous thermal polymerization of alkyl acrylates was foun to take place in the absence of any known ae external initiators [6]. Such a process coul provie an atactive alternative for the prouction of low molecular weight polymers of alkyl acrylates. Initial stuies [7], [8] T. Rier, S. Srinivasan an M. Soroush are with the Department of Chemical an Biological Engineering, Drexel University, Philaelphia, PA 1914, U.S.A. G.A. Kalfas an M.C. Gray are with DuPont Experimental Station, Wilmington, DE 1988, USA A.M. Rappe is with Makineni Theoretical Laboratories, Department of Chemisy, University of Pennsylvania, Philaelphia, PA, USA Corresponing author: M. Soroush; (phone); (fax); ms1@rexel.eu This work was supporte by the National Science Founation Grant CBET using eleconspray ionization-fourier ansform mass specomey (ESI-FTMS), nuclear magnetic resonance (NMR) specoscopy an etaile analytic macroscopic moeling were unable to reveal the ue nature of the initiating species an the mechanism of initiation. More recently, using quantum chemical calculations [9], [1], it has been shown that monomer self-initiation in alkyl acrylate polymerization is the ominant initiating mechanism. Further evience to the occurrence of the self-initiation mechanism was shown by accurately fitting the compute suctures to maix-assiste laser esorption ionization-time of flight mass specomey [11] an ESI-FTMS specum [12] of poly-alkyl acrylate samples that were generate using the spontaneous thermal polymerization process. Macroscopic kinetic moels are important tools to unerstan complex polymerization processes [12],. The influence of various parameters on the overall kinetics an thermoynamics of the process can be ientifie with the evelopment of reasonably accurate mechanistic moels. The avantage of mechanistic over semi-empirical moels is in their ability to account for the rate constants of iniviual reactions an provie information about the impact of each of these reactions on the overall polymerization kinetics [14]. Also, the kinetics compute using quantum chemical calculations can be use in the mechanistic moels to improve the preiction of variables such as monomer conversion, average molecular weights an chain branches per 1 monomer units [15]. A macroscopic moel can be use to calculate optimal reactor operating conitions, esign moel-base conollers, conuct experimental esign, an esign proucts with improve quality [16]. A macroscopic mechanistic moel can be evelope using several approaches. The black box approach is base upon empirical correlations that are obtaine via well-esigne experiments, but have limite preictive power. On the other han, the etaile analytic approach is capable of escribing the composition of polymer mass at any time instant by using the information about elementary kinetic mechanisms. However, the lack of aequate quantitative knowlege on many iniviual reactions in a polymerization process has limite the preiction capabilities of the etaile analytical moels. A thir moeling metho is that of tenency moeling, in which number of monomer units in polymer chains is not accounte for, an the Schulz- Flory isibution is use to obtain the moments of live an ea polymer chains. While this simplification reuces the complexity involve in solving the moel equations, /11/$ AACC 371

2 the tenency moeling approach provies little information on the iniviual raicals an their chain length uring polymerization [17]. Propagation rate constants in free-raical polymerization of styrene, methyl methacrylate an alkyl acrylates have been measure through pulse laser polymerization-size exclusion chromatography (PLP-SEC) technique [18], [19], [2]. Because this technique is unable to accurately etermine rate constants at high temperatures, alternative methos shoul be use for the estimation of the rate constants in high-temperature polymerization. One of these alternative methos is macroscopic kinetic moeling combine with estimation of the reaction rate constants from measurements obtaine from polymer sample analyses [8]. Previous macroscopic mechanistic moeling stuies of spontaneous thermal polymerization of styrene an methyl methacrylate have provie qualitative evience to the occurrence of selfinitiation [21], [22]. Styrene an methyl methacrylate selfinitiation rate constants have been reporte [22], [8] to have an orer greater than two. Extensive stuies of macroscopic kinetic moeling of high-temperature polymerization of alkyl acrylates, initiate by ae thermal initiators, have reporte rate constants of various reactions [8], [23]. Rantow et al. [8] explore spontaneous polymerization of n-butyl acrylate using a etaile analytic moeling approach an were able to preict the activation energy an rate constants of spontaneous initiation, backbiting an β-scission reactions. However, no funamental unerstaning of the initiation mechanism was obtaine from this stuy. In this paper, a macroscopic mechanistic kinetic moel is evelope for self-initiate polymerization of n-butyl acrylate (nba) in a batch reactor. The rate constant of the selfinitiation reaction is estimate from monomer conversion measurements, an the moel is valiate using a ifferent set of measurements. The moel is then use to calculate an optimal temperature profile for a batch polymerization reactor. The organization of the rest of this paper is as follows. Section II escribes the evelopment of the mathematical moel. Section III iscusses the parameter estimation. Section IV presents the optimization stuy. Finally, concluing remarks are given in Section V. A. Reaction Mechanisms II. MATHEMATICAL MODELING The reaction scheme consiere here for thermally selfinitiate free-raical polymerization of nba involves the following reactions (Table 1): thermal self-initiation, chain propagation, ina-molecular chain ansfer to polymer (backbiting), inter-molecular chain ansfer to polymer, chain ansfer to solvent, an chain ansfer to monomer, β- scission, termination by isproportionation an termination by combination. Here, M an S represent the monomer an solvent, respectively. D n enotes a ea polymer chain with n monomer units but without a terminal ouble bon. U n is a ea polymer chain with n monomer units an a terminal TABLE I REACTION MECHANISMS a. Monomer self-initiation reactions k 3M i R 1 + R2 b. Propagation reactions R n + M k p R n+1 R n + M kt p R n+1 (+LCB) + M kt p +1 (+SCB) + U m k mac c. Backbiting reactions (n > 2) R n+m k bb R n. β scission reactions (n > 3) k β U k β U n 2 + R2 k β U n m + Rm k β U m + m e. Intermolecular chain ansfer to polymer reactions + D k P m Dn + Rm k P + U m Dn + Rm f. Chain ansfer to monomer reactions + M k M Dn + R1 + M k t M Dn + R1 + M k t M Dn + R1 g. Chain ansfer to solvent reactions + S k S Dn + R + S k t S Dn + R + S k t S Dn + R h. Termination by coupling reactions R n + R m R n + R m R n + R m + Rm + Rm k t c k t tc k t tc tc tc + Rm ktc tt i. Termination by isproportionation reactions + Rm k t + k R t t m + Rm kt t + Rm + Rm + Rm t t kt tt 372

3 ouble bon. is a seconary raical with n monomer units. is a tertiary raical with n monomer units forme by the intermolecular chain ansfer to polymer reactions. represents a tertiary raical with n monomer units generate by the backbiting reactions. SCB enotes a short branching point, an LCB a long chain branching point. The self-initiation reaction is consiere to occur via the iraical mechanism of self-initiation reporte in [1]. The overall self-initiation reaction is inclue in the mechanism given above. It follows a secon orer reaction. The seconary raicals forme by the self-initiation reaction react with a monomer, leaing to chain propagation. Tertiary raicals forme via chain ansfer to polymer an backbiting either propagates to form short chain branches or unergo β-scission. A β-scission reaction generates a ea polymer with a terminal ouble bon an a seconary raical. The β- scission reactions atibute to lowering of the polymer chain length an the average molecular weights. B. Rate Equations On the basis of the reaction mechanisms liste in Table 1, assuming elementary reactions one can erive rate expressions for the various species involve. The prouction rates of M, S, R, R 1, R 2, R 3, R, R, R, D an U, respectively, are: r M = 3k i [M] 2 k p [M][R ] k t p[m]([r ] + [R ]) k M [M][R ] k t M [M]([R ] + [R ]) r S = k S [S][R ] k t S [S]([R ] + [R ]) r R = k S [S][R ] + k t S [S]([R ] + [R ]) k p [M][R ] k mac [R ][U].5k t [R ][R ].5k t t([r ] + [R ])[R ] r R 1 = k i [M] 2 + k p [M]([R ] [R 1]) + k M [M]([R ] R 1]) + k t M [M]([R ] + [R ]) k S [S][R 1] k mac [U][R 1].5k t [R ][R 1].5k t t([r ] +[R ])[R 1] k P [R 1]([U] + [D]) r R 2 = k i [M] 2 + k p [M]([R 1] [R 2]) + k β [R ] k mac [U][R 2] k P ([D] + [U])[R 2] (k M [M] + k S [S])[R 2].5k t [R ][R 2].5k t t([r ] + [R ])[R 2] r R 3 = k p [M]([R 2] [R 3]) k mac [U][R 3] k P ([D] +[U])[R 3] (k M [M] + k S [S])[R 3].5k t [R ][R 3].5k t t([r ] + [R ])[R 3] r R = 2k i [M] 2 + k t p[m]([r ] + [R ]) k mac [U][R ] k bb [R ] + 2k β ([R ] + [R ]) k P ([D] +[U])[R ] + (k t M [M] + k t S [S])([R ] +[R ]) k t [R ] 2.5k t t([r ] + [R ])[R ] r R = kp[m][r t ] 2k β [R ] + k P ([D] + [U])[R ] (k t M [M] + k t S [S])[R ].5kt[R t ][R ].5kt tt (2[R ] + [R ])[R ] + k mac [U][R ] r R = kp[m][r t ] + k bb [R ] 2k β [R ] (k t M [M] + k t S [S])[R ].5kt[R t ] [R ].5kt tt ([R ] + 2[R ])[R ] r D = k P [R ][U] + (k M [M] + k S [S])[R ] +(k t M [M] + k t S [S])([R ] + [R ]) +.5k t [R ] 2 +.5kt[R t ]([R ] + [R ]) +.5kt tt ([R ][R ] + [R ] 2 + [R ] 2 ) r U = 2k β ([R ] + [R ]) k P [R ][U] k mac [R ][U] +.5k t [R ] 2 +.5kt t [R ] ([R ] + [R ]) +.5kt tt ([R ][R ] +[R ] 2 + [R ] 2 ) where [z] enotes the molar concenation of species z, R is a seconary raical, R is a tertiary raical forme by the intermolecular chain ansfer to polymer reactions, an R represents a tertiary raical generate by the backbiting reactions. C. Batch Reactor Moel Mole balances for all species involve yiel a mathematical moel of the batch reactor: [M] = r M, [M]() = [M] [S] = r S, [S]() = [S] [R] = r R, [R]() = [R1] = r R 1, [R1]() = [R2] = r R 2, [R2]() = [R3] = r R 3, [R3]() = [R ] = r R, [R ]() = [R ] = r R, [R ]() = [R ] = r R, [R ]() = [D] = r D, [D]() = [U] = r U, [U]() = where [M] an [S] are the initial concenations of M an S, respectively. The number-average molecular weight of the ea polymer chains, M n, an the monomer conversion,, are then obtaine from: M n = [M] [M] M m [D] + [U] = 1 [M] [M] 373

4 TABLE II REACTION RATE CONSTANT VALUES. Parameter Frequency Factor Activation Ref. Energy kj.mol 1 k p L.mol 1.s [28] kp t L.mol 1.s [29] k bb s [3] k M L.mol 1.s [12] k t L.mol 1.s [31] kt tt L.mol 1.s k β s k P L.mol 1.s [32] C S Measurements Moel Preiction Time (min) Fig. 1. Moel preiction fitte to the measurement set 1 at 14 C TABLE III OTHER KINETIC PARAMETER VALUES. Parameter Dimensionless Value δ s.1 δ st.7 δ t.9 γ III. PARAMETER ESTIMATION Parameter estimation was carrie out using least-squareerror minimization an measurements of monomer conversion from isothermal spontaneous thermal polymerization of nba (4 wt% nba an 6 wt% xylene) at 14 an 16 C [8]. The minimization was conucte using the optimization routine fminsearch from the Optimization Toolbox of MATLAB R21b. Except for the frequency factor an activation energy of the self-initiation reaction that were estimate in this stuy, the other reaction kinetic parameter values were taken from the literature (Table 2). Two ifference sets of conversion measurements (set 1 an set 2) were use in this stuy. Each set inclue gravimeic measurements of monomer conversion from isothermal batch operations at 14 an 16 C. First, the self-initiation reaction rate constant at 14 C was estimate from set-1 conversion measurements obtaine at 14 C. Secon, the self-initiation reaction rate constant at 16 C was estimate from set-1 conversion measurements from batch operations at 16 C. Thir, using the Arrhenius plot, the frequency factor an activation energy of the selfinitiation reaction were calculate: Activation Energy = kj.mol 1 Frequency Factor = L.mol 1.s 1 The quality of the moel fitting the conversion ata set 1 can be seen in Figures 1 an 2. The secon set of conversion ata (set 2) at 14 an 16 C was then use for the moel valiation. The valiation results are shown in Figures 3 an 4; as can be seen, the moel preicts the conversion reasonably well. IV. BATCH REACTOR OPTIMIZATION Using the reactor moel presente in Section IIIB, we calculate an optimal temperature profile, T (t), for a batch.2 Measurements Moel Preiction Fig. 2. Moel preiction fitte to the measurement set 1 at 16 C Measurements Moel Preiction Fig. 3. Moel valiation: moel preiction an measurement set 2 at 14 C Measurements Moel Preiction Fig. 4. Moel valiation: moel preiction an measurement set 2 at 16 C 374

5 TABLE IV REACTION RATE CONSTANT DEFINITIONS/RELATIONS. k t = k tc + k t kt t = kt t c + kt t kt tt = kt tt c + kt tt k t = δ s k t kt tt = δ t kt tt kt t = δ st k t kt tt k tc = (1 δ s )k t kt tt c = (1 δ t )k t tt kt t c = (1 δ st ) k t kt tt k S = C S k p k t M = kt p k p k M k t S = kt p k S k p k mac = γk p T (K) Time(min) Fig. 5. Initial guess for the reactor temperature profile 49 nba polymerization reactor with a batch cycle time of 7 hours (incluing the heating-up an cooling-own perios). The batch reactor temperature shoul be K at the beginning an the en of each batch cycle; T () = K an T (7 36) = K. The reactor is initially loae with monomer (4 wt.%) an solvent (6 wt.%). The optimal temperature profile shoul minimize the performance inex: [ [1 (t f )] M ] 2 n(t f ) 5 where t f = 7 36 s. In other wors, the optimal temperature profile shoul yiel a polymer prouct with a monomer conversion an a number-average molecular weight as close as possible to one an 5, kg.kmol 1, respectively. To simplify the optimization problem, we parameterize the reactor temperature profile using a fourth-orer polynomial in time; that is, where T (t) = at 4 + bt 3 + ct e e = K, = (at 3 f + bt 2 f + ct f ) K.s 1 This parameterization of the optimal temperature profile reuce the search for an optimal profile to a search for three constants, a, b an c. The optimization routine fminsearch from the Optimization Toolbox of MATLAB R21b with the initial guesses: a = b = c = was use to calculate the optimal values of a, b an c. These initial guesses correspon to the initial-guess temperature profile shown in Figure T (K) Fig. 6. Time(min) Optimal reactor temperature profile The optimization routine calculate the optimal values: a = b = c = which correspon to the optimal temperature profile shown in Figure 6. Enforcing the optimal temperature profile to the reactor leas to the prouction of a polymer with (t f ) =.814 an M n (t f ) = 4918 kg.kmol 1. The optimization results shown herein were for a batch reactor; they were presente to just show an application of the moel. A realistic case is that of a semi-batch reactor in which the monomer is ae to the reactor over time such that the mass fraction of free (unreacte) monomer in the reactor never excees 1%. This upper limit on the monomer mass fraction is to ensure the reactor safety. In this case, the entire solvent is loae initially, an the optimization variables are the reactor temperature an the monomer fee flow rate. V. CONCLUSIONS A macroscopic mechanistic mathematical moeling an optimization stuy of a batch polymerization reactor in which self-initiate free-raical polymerization of n-butyl acrylate at 14 an 16 C takes place, was presente. 375

6 The moel was obtaine using a comprehensive free-raical polymerization reaction mechanism. The rate constant of the monomer self-initiation was estimate from monomer conversion measurements. The moel valiation results showe that the macroscopic mechanistic moel is accurate enough for optimizing the self-initiate polymerization process to prouce high quality acrylic resins. On the basis of the moel, optimization was conucte, an an optimal batchreactor temperature profile that yiels an en-batch polymer prouct with esire properties (conversion an numberaverage molecular weight) was calculate. Other applications of the moel inclue moel-base conoller esign, reactor scaleup an safety analysis. VI. ACKNOWLEDGEMENT This material is base upon work supporte by the National Science Founation uner Grant No. CBET The authors thank Mona Bavarian for her invaluable help with the simulations, an gratefully acknowlege reviewers comments. REFERENCES [1] Report on Automotive Coatings Ahesives an Sealants, Freeonia Group, May 28; [2] Cunningham, M.F., Hutchinson, R. Inusial applications an processes. in Hanbook of Raical Polymerization, K. Matyjaszewski an T. P. Davis, Eitors, John Wiley an Sons Inc., Hoboken, NJ, 22, p 333. [3] VOC s Directive, EU Committee of the American Chamber of Commerce in Belgium, ASBL/VZw, Brussels, July 8, [4] Cheifari, J., Jeffery, J., Mayaunne, R.T.A., Moa G., Rizzaro, E., Thang, S.G., Chain Transfer to Polymer: A Convenient Route to Macromonomers, Macromolecules, vol. 32, 1999, p 77. [5] Peter, A.C. Moskal, M.G. Arkema Inc., King of Prussia, PA. [6] Gray, M.C., Quan, C., Soroush, M., Thermally initiate polymerization process. US Patent Application Number 6/484,393, file on July 2, 23. [7] Quan, C., Soroush, M., Gray, M.C., Hansen, J.E., Simonsick, W.J., Characterization of Thermally Polymerize n-butyl Acrylate an Ethyl Acrylate, Macromolecules, vol. 38, 25, pp [8] Rantow, F.S., Soroush M., Gray, M.C., Kalfas, G.A., Spontaneous Polymerization an Chain Microsucture Evolution in High- Temperature Solution Polymerization of n-butyl Acrylate, Polymer, vol. 47, 26, pp [9] Srinivasan, S., Lee, M.W, Gray, M.C., Soroush, M., Rappe, A.M., Computational Stuy of Self-initiation Mechanism in Thermal Polymerization of Methyl Acrylate, J. Phys. Chem. A, vol. 113, 29, pp [1] Srinivasan, S., Lee, M.W., Gray, M.C., Soroush, M., Rappe, A.M., Self-Initiation Mechanism in Spontaneous Thermal Polymerization of Ethyl an n-butyl Acrylate: a Theoretical Stuy, J. Phys. Chem. A, vol. 114, 21, pp [11] Srinivasan, S., Kalfas, G.A., Petkovska, V.I., Bruni, C., Gray, M.C., Soroush, M., Experimental stuy of the spontaneous thermal homopolymerization of methyl an n-butyl acrylate, J. Appl. Pol. Sci., Vol. 118, 21, pp [12] Maeer, S., Gilbert, R.G., Measurement of Transfer Constant for Butyl Acrylate Free-Raical Polymerization, Macromolecules, vol. 33, 1998, pp Nikitin, A.N., Hutchinson, R.A., Wang, W., Kalfas, G.A., Richars, J.R., Bruni, C., Effect of Inamolecular Transfer to Polymer on Stationary Free-Raical Polymerization of Alkyl Acrylates, 5 Consieration of Solution Polymerization up to High Temperatures, Macromolecular Reaction Engin., vol. 4, 21, pp [14] Wang, W., Nikitin, A., Hutchinson, R., Consieration of Macromonomer Reactions in n-butyl Acrylate Free Raical Polymerization, Macromol. Rapi Commn., 29, vol. 3, pp [15] Arzameni, G., Plessis, C., Leiza, J., Asua, J., Effect of the Inamolecular Chain Transfer to Polymer on PLP/SEC Experiments of Alkyl Acrylates. Macromolecular Theory an Simulations, vol. 12, 23, pp [16] Hakim, M., Verhoeven, V., McManus, N.T., Dube, M.A., Penliis, A. High-Temperature Solution Polymerization of Butyl Acrylate/Methyl Methacrylate: Reactivity Ratio Estimation, J. Appl. Poly. Sci., vol. 77, 2, pp [17] Villermaux, J.; Blavier, L., Free Raical Polymerization Engineering I: A New Metho for Moeling Free Raical Homogeneous Polymerization Reactions, Chem. Eng. Sci., vol. 39, 1984, pp [18] Gilbert, R.G., Critically-Evaluate Propagation Rate Coefficients in Free Raical Polymerizations I. Styrene an Methyl Methacrylate, Pure & Appl. Chem., vol. 68, 1996, pp [19] Asua, J.M., Beuermann, S., Buback, M., Castignolles, P., Charleux, B., Gilbert, R.G., Hutchinson, R.A., Leiza, J.R., Nikitin, A.N., Vairon, J. P., Van Herk, A.M., Critically Evaluate Rate Coefficients for Free Raical Polymerizations, Propagation Rate Coefficient for Butyl Acrylate, Macromol. Chem., vol. 25, 24, pp [2] Willemse, R..E., Van Herk, A.M., Panchenko, E., Junkers, T., Buback, M., PLP-ESR Monitoring of Michain Raicals in n-butyl Acrylate Polymerization, Macromolecules, 25, 38, pp [21] Pryor, W.A., Lasswell, L.D., Avances in free raical chemisy, Acaemic Press, New York, 1975, Vol. V. [22] Stickler, M., Meyerhoff, G., The Spontaneous Thermal Polymerization of Methyl Methacrylate. 5. Experimental Stuy an Computer Simulation of the High Conversion Reaction at 13 C, Polymer, vol. 22, 1981, pp [23] Blavier, L., Villermaux, J., Free Raical Polymerization Engineering II: Moeling of Homogeneous Polymerization of Styrene in a Batch Reactor, Influence of Initiator. Chem. Eng. Sci. vol. 39, 1984, pp [24] Peck, A.N.F., Hutchinson, R.A., Gray, M.C., Branching an Scission Reactions in High Temperature Acrylate Polymerizations, Polymer Preprints, vol. 43, 22, p 154. [25] Peck, A.N.F., Hutchinson, R.A., Seconary Reactions in the High Temperature Free Raical Polymerization of Butyl Acrylate, Macromolecules, vol. 37, 24, pp [26] Subrahmanyam, B., Baruah, S.; Rahman, M., Laskar, N., Mazumer, R., Stuies on High Conversion Polymerization of n-alkyl Acrylates, Polymer, vol , pp [27] Willemse, R..E, van Herk, A.M., Determination of Propagation Rate Coefficients of a Family of Acrylates with PLP-MALDI-ToF-MS, Macromol. Chem. an Phy., vol. 211, 21, pp [28] Asua, J.M., Beuermann, S., Buback, M., Castignolles, P., Charleux, B., Gilbert, R.G., Hutchinson, R., Leiza, J.R., Nikitin, A.N., Vairon, J.P., van Herk, A.M., Critically Evaluate Rate Coefficients for Free- Raical Polymerization, 5 Propagation Rate Coefficient for Butyl Acrylate, Macromol. Chem. Phys., vol. 25(16), 24, pp [29] Hesse, P., Ph.D. Thesis, University of Gottingen, 28. [3] Nikitin, A.N., Hutchinson, R.A., Wang, W., Kalfas, G.A., Richars, J.R., Bruni, C., The Effect of Inamolecular Transfer to Polymer on Stationary Free-Raical Polymerization of Alkyl Acrylates: 3. Consieration of Solution Polymerization up to High Conversions, Macromolecular Theory an Simulations, vol. 4, 29, pp [31] Nikitin, A.N., Hutchinson, R.A., Polymerization of Alkyl Acrylates: 2. Improve Consieration of Termination, Macromolecular Theory an Simulations, vol. 15, 26, pp [32] Arzameni, G., Plessis, C., Leiza, J.R., Asua, J.M., Effect of the Inamolecular Chain Transfer to Polymer on PLP/SEC Experiments of Alkyl Acrylates, Macromolecular Theory an Simulations, vol. 12, 23, pp

polymerization of n-butyl acrylate

polymerization of n-butyl acrylate SUPPORTING INFORMATION Intermolecular transfer to polymer in the radical polymerization of n-butyl acrylate Nicholas Ballard, 1 Shaghayegh Hamzehlou, 1 José M. Asua 1* 1 POLYMAT and Kimika Aplikatua Saila,

More information

Theoretical Study of Chain Transfer to Solvent Reactions of Alkyl Acrylates

Theoretical Study of Chain Transfer to Solvent Reactions of Alkyl Acrylates pubs.acs.org/jpca Theoretical Study of Chain Transfer to Solvent Reactions of Alkyl Acrylates Nazanin Moghadam, Sriraj Srinivasan, Michael C. Grady, # Andrew M. Rappe, and Masoud Soroush*, Department of

More information

Theoretical Study of Intermolecular Chain Transfer to Polymer Reactions of Alkyl Acrylates

Theoretical Study of Intermolecular Chain Transfer to Polymer Reactions of Alkyl Acrylates pubs.acs.org/iecr Theoretical Study of Intermolecular Chain Transfer to Polymer Reactions of Alkyl Acrylates Nazanin Moghadam, Shi Liu, Sriraj Srinivasan, Michael C. Grady, # Andrew M. Rappe, and Masoud

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process Latest Trens in Circuits, Automatic Control an Signal Processing Optimal operating strategies for semi-batch reactor use for chromium sluge regeneration process NOOSAD DAID, MACKŮ LUBOMÍR Tomas Bata University

More information

Research Article When Inflation Causes No Increase in Claim Amounts

Research Article When Inflation Causes No Increase in Claim Amounts Probability an Statistics Volume 2009, Article ID 943926, 10 pages oi:10.1155/2009/943926 Research Article When Inflation Causes No Increase in Claim Amounts Vytaras Brazauskas, 1 Bruce L. Jones, 2 an

More information

Computational Study of Cyclohexanone Monomer Co-initiation Mechanism in Thermal Homo-polymerization of Methyl Acrylate and Methyl Methacrylate

Computational Study of Cyclohexanone Monomer Co-initiation Mechanism in Thermal Homo-polymerization of Methyl Acrylate and Methyl Methacrylate pubs.acs.org/jpca Computational Study of Cyclohexanone Monomer Co-initiation Mechanism in Thermal Homo-polymerization of Methyl Acrylate and Methyl Methacrylate Shi Liu, Sriraj Srinivasan, Michael C. Grady,

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

Process Design for Operability: Louis P. Russo and B. Wayne Bequette

Process Design for Operability: Louis P. Russo and B. Wayne Bequette Process Design for Operability: A Styrene Polymerization Application Louis P. Russo an B. Wayne Bequette Department of Chemical Engineering, Rensselaer Polytechnic nstitute, Troy, NY USA Abstract: Steay-state

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late

More information

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:

More information

A Molten Solid Approach for Simulating Urea-Water Solution Droplet Depletion

A Molten Solid Approach for Simulating Urea-Water Solution Droplet Depletion ILASS Americas 7th Annual Conference on Liqui Atomization an Spray Systems, Raleigh, NC, May 015 A Molten Soli Approach for Simulating Urea-Water Solution Droplet Depletion Shaoping Quan* 1, Mingjie Wang

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Basic IIR Digital Filter Structures

Basic IIR Digital Filter Structures Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Kinetics. Chemical kinetics contends with a number of issues. Specfically, it seeks to do the following:

Kinetics. Chemical kinetics contends with a number of issues. Specfically, it seeks to do the following: Kinetics Chemical kinetics contens with a number of issues. pecfically, it seeks to o the following: Detere an preict the rate of chemical reactions Detere the reaction mechanism. With the previous goal,

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

PoS(RAD COR 2007)030. Three-jet DIS final states from k -dependent parton showers. F. Hautmann University of Oxford

PoS(RAD COR 2007)030. Three-jet DIS final states from k -dependent parton showers. F. Hautmann University of Oxford OUTP-8--P Three-jet DIS final states from k -epenent parton showers University of Oxfor E-mail: hautmann@thphys.ox.ac.uk H. Jung DESY E-mail: jung@mail.esy.e Experimental measurements of angular correlations

More information

Integrated Data Reconciliation with Generic Model Control for the Steel Pickling Process

Integrated Data Reconciliation with Generic Model Control for the Steel Pickling Process Korean J. Chem. Eng., (6), 985-99 (3) Integrate Data Reconciliation with Generic Moel Control for the Steel Picling Process Paisan Kittisupaorn an Pornsiri Kaewprait Department of Chemical Engineering,

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

Influence of Radiation on Product Yields in a Film Boiling Reactor

Influence of Radiation on Product Yields in a Film Boiling Reactor R&D NOTES Influence of Raiation on Prouct Yiels in a Film Boiling Reactor C. Thomas Aveisian, Wing Tsang, Terence Daviovits, an Jonah R. Allaben Sibley School of Mechanical an Aerospace Engineering, Cornell

More information

Analyzing the discrepancies in the activation energies of the backbiting and β-scission reactions in the radical polymerization of n-butyl acrylate

Analyzing the discrepancies in the activation energies of the backbiting and β-scission reactions in the radical polymerization of n-butyl acrylate Analyzing the discrepancies in the activation energies of the backbiting and β-scission reactions in the radical polymerization of n-butyl acrylate S. Hamzehlou, a N. Ballard, a Y. Reyes, b A. Aguirre,

More information

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS - PAPER NO. 79-2108 PRAIRIE RAINFALL,CHARACTERISTICS G.E. Dyck an D.M. Gray Research Engineer an Chairman Division of Hyrology University of Saskatchewan Saskatoon, Saskatchewan, Canaa For presentation

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Consider for simplicity a 3rd-order IIR filter with a transfer function. where

Consider for simplicity a 3rd-order IIR filter with a transfer function. where Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Deriving ARX Models for Synchronous Generators

Deriving ARX Models for Synchronous Generators Deriving AR Moels for Synchronous Generators Yangkun u, Stuent Member, IEEE, Zhixin Miao, Senior Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract Parameter ientification of a synchronous generator

More information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive Control of a Laboratory Time Delay Process Experiment Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results

More information

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL 386 EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL Lien-Kwei CHIEN 1, Yan-Nam OH 2 An Chih-Hsin CHANG 3 SUMMARY In this stuy, the fille material in Yun-Lin

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Available online at ScienceDirect. Procedia Engineering 100 (2015 )

Available online at   ScienceDirect. Procedia Engineering 100 (2015 ) Available online at www.scienceirect.com ScienceDirect Proceia Engineering 1 (215 ) 36 369 25th DAAAM International Symposium on Intelligent Manufacturing an Automation, DAAAM 214 Enhancement of a Semi-Batch

More information

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution International Journal of Statistics an Systems ISSN 973-675 Volume, Number 3 (7), pp. 475-484 Research Inia Publications http://www.ripublication.com Designing of Acceptance Double Sampling Plan for Life

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

TCP throughput and timeout steady state and time-varying dynamics

TCP throughput and timeout steady state and time-varying dynamics TCP throughput an timeout steay state an time-varying ynamics Stephan Bohacek an Khushboo Shah Dept. of Electrical an Computer Engineering Dept. of Electrical Engineering University of Delaware University

More information

More from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives.

More from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives. Math 1314 ONLINE More from Lesson 6 The Limit Definition of the Derivative an Rules for Fining Derivatives Eample 4: Use the Four-Step Process for fining the erivative of the function Then fin f (1) f(

More information

Quality competition versus price competition goods: An empirical classification

Quality competition versus price competition goods: An empirical classification HEID Working Paper No7/2008 Quality competition versus price competition goos: An empirical classification Richar E. Balwin an Taashi Ito Grauate Institute of International an Development Stuies Abstract

More information

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization Web-Base Technical Appeni: Multi-Prouct Firms an Trae Liberalization Anrew B. Bernar Tuck School of Business at Dartmouth & NBER Stephen J. Reing LSE, Yale School of Management & CEPR Peter K. Schott Yale

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

A Look at the ABC Conjecture via Elliptic Curves

A Look at the ABC Conjecture via Elliptic Curves A Look at the ABC Conjecture via Elliptic Curves Nicole Cleary Brittany DiPietro Alexaner Hill Gerar D.Koffi Beihua Yan Abstract We stuy the connection between elliptic curves an ABC triples. Two important

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlle Inuction Motor Drive V. S. S. Pavan Kumar Hari, Avanish Tripathi 2 an G.Narayanan 3 Department of Electrical Engineering,

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Some New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,

Some New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco, Jeju, Korea, April 16-20, 2017, on USB 2017 Some New Thoughts on the Multipoint Metho for Reactor Physics Applications Sanra Dulla, Piero Ravetto, Paolo Saracco, Politecnico i Torino, Dipartimento Energia,

More information

Topic Modeling: Beyond Bag-of-Words

Topic Modeling: Beyond Bag-of-Words Hanna M. Wallach Cavenish Laboratory, University of Cambrige, Cambrige CB3 0HE, UK hmw26@cam.ac.u Abstract Some moels of textual corpora employ text generation methos involving n-gram statistics, while

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences

More information

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints EUROSTEEL 2, August 3 September 2, 2, Buapest, Hungary A NUMERICAL EVALUATION OF CHS T JOINTS UNDER AXIAL LOADS Raphael S. a Silva a, Luciano R. O. e Lima b, Pero C. G. a S. Vellasco b, José G. S. a Silva

More information

Real-time optical characterization of GaP heterostructures by p-polarized reflectance

Real-time optical characterization of GaP heterostructures by p-polarized reflectance Thin Soli Films 3334 998 6469 Real-time optical characterization of GaP heterostructures by p-polarize reflectance N. Dietz a,, K. Ito b a Department of Physics, North Carolina State Uniersity, Raleigh,

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Critical Size and Particle Growth

Critical Size and Particle Growth Critical Size an article Growth rof. Sotiris E. ratsinis article Technology Laboratory Department of Mechanical an rocess Engineering, ETH Zürich, Switzerlan www.ptl.ethz.ch 1 Nucleation-Conensation A

More information

Real-time economic optimization for a fermentation process using Model Predictive Control

Real-time economic optimization for a fermentation process using Model Predictive Control Downloae from orbit.tu.k on: Nov 5, 218 Real-time economic optimization for a fermentation process using Moel Preictive Control Petersen, Lars Norbert; Jørgensen, John Bagterp Publishe in: Proceeings of

More information

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function

More information

Design and performance of in vitro transcription rate regulatory circuits

Design and performance of in vitro transcription rate regulatory circuits Conference on Decision an Control (CDC), 2008 http://www.cs.caltech.eu/~murray/papers/2008j_fm08-cc.html Design an performance of in vitro transcription rate regulatory circuits Elisa Franco an Richar

More information

Damage detection of shear building structure based on FRF response variation

Damage detection of shear building structure based on FRF response variation , pp.18-5 http://x.oi.org/10.1457/astl.013.3.05 Damage etection of shear builing structure base on FRF response variation Hee-Chang Eun 1,*, Su-Yong Park 1, Rae-Jung im 1 1 Dept. of Architectural Engineering,

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Classical Quantum Theory

Classical Quantum Theory J.P. Wesley Weiherammstrasse 4 78176 Blumberg, Germany Classical Quantum Theory From the extensive observations an the ieas of Newton an from classical physical optics the velocity of a quantum particle

More information

2Algebraic ONLINE PAGE PROOFS. foundations

2Algebraic ONLINE PAGE PROOFS. foundations Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

An Anisotropic Hardening Model for Springback Prediction

An Anisotropic Hardening Model for Springback Prediction An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS

More information

University of Groningen. Rheokinetics Cioffi, Mario

University of Groningen. Rheokinetics Cioffi, Mario University of Groningen Rheokinetics Cioffi, Mario IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document version below.

More information

Text S1: Simulation models and detailed method for early warning signal calculation

Text S1: Simulation models and detailed method for early warning signal calculation 1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Collective optical effect in complement left-handed material

Collective optical effect in complement left-handed material Collective optical effect in complement left-hane material S.-C. Wu, C.-F. Chen, W. C. Chao, W.-J. Huang an H. L. Chen, National Nano Device Laboratories, 1001-1 Ta-Hsueh Roa, Hsinchu, Taiwan R.O.C A.-C.

More information

A check digit system over a group of arbitrary order

A check digit system over a group of arbitrary order 2013 8th International Conference on Communications an Networking in China (CHINACOM) A check igit system over a group of arbitrary orer Yanling Chen Chair of Communication Systems Ruhr University Bochum

More information

Event based Kalman filter observer for rotary high speed on/off valve

Event based Kalman filter observer for rotary high speed on/off valve 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 28 WeC9.6 Event base Kalman filter observer for rotary high spee on/off valve Meng Wang, Perry Y. Li ERC for Compact

More information

18 EVEN MORE CALCULUS

18 EVEN MORE CALCULUS 8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;

More information

New Statistical Test for Quality Control in High Dimension Data Set

New Statistical Test for Quality Control in High Dimension Data Set International Journal of Applie Engineering Research ISSN 973-456 Volume, Number 6 (7) pp. 64-649 New Statistical Test for Quality Control in High Dimension Data Set Shamshuritawati Sharif, Suzilah Ismail

More information