Numerical modeling of self-propagating polymerization fronts: The role of kinetics on front stability.

Frontal propagation of a highly exothermic polymerization reaction in a liquid is studied with the goal of developing a mathematical model of the process. As a model case we consider monomers such as methacrylic acid and n-butyl acrylate with peroxide initiators, although the model is not limited to these reactants and can be applied to any system with the similar basic polymerization mechanism. A three-step reaction mechanism, including initiation, propagation and termination steps, as well as a more simple one-step mechanism, were considered. For the one-step mechanism the loss of stability of propagating front was observed as a sequence of period doubling bifurcations of the front velocity. It was shown that the one-step model cannot account for less than 100% conversion and product inhomogeneities as a result of front instability, therefore the three-step mechanism was exploited. The phenomenon of superadiabatic combustion temperature was observed beyond the Hopf bifurcation point for both kinetic schemes and supported by the experimental measurements. One- and two-dimensional numerical simulations were performed to observe various planar and nonplanar periodic modes, and the results for different kinetic schemes were compared. It was found that stability of the frontal mode for a one-step reaction mechanism does not differ for 1-D and 2-D cases. For the three-step reaction mechanism 2-D solutions turned out to be more stable with respect to the appearance of nonplanar periodic modes than corresponding 1-D solutions. Higher Zeldovich numbers (i.e., higher effective activation energies or lower initial temperatures) are necessary for the existence of planar and nonplanar periodic modes in the 2-D reactor with walls than in the 1-D case. (c) 1997 American Institute of Physics.

[1]  G. Odian,et al.  Principles of polymerization , 1981 .

[2]  J. Pojman,et al.  Free-Radical Frontal Polymerization: Self-Propagating Thermal Reaction Waves , 1996 .

[3]  Gregory I. Sivashinsky,et al.  PROPAGATION OF A PULSATING REACTION FRONT IN SOLID FUEL COMBUSTION , 1978 .

[4]  István Nagy,et al.  THERMOCHROMIC COMPOSITE PREPARED VIA A PROPAGATING POLYMERIZATION FRONT , 1995 .

[5]  J. Pojman,et al.  Factors affecting propagating fronts of addition polymerization: Velocity, front curvature, temperatue profile, conversion, and molecular weight distribution , 1995 .

[6]  John A. Pojman,et al.  Spin mode instabilities in propagating fronts of polymerization , 1995 .

[7]  D. Frank-Kamenetskii,et al.  Diffusion and heat exchange in chemical kinetics , 1955 .

[8]  John A. Pojman,et al.  Convective instabilities in traveling fronts of addition polymerization , 1992 .

[9]  J. Pojman,et al.  The true molecular weight distributions of acrylate polymers formed in propagating fronts , 1996 .

[10]  N. S. Enikolopyan,et al.  Theory of polymerization-front propagation , 1974 .

[11]  J. Pojman.  Traveling fronts of methacrylic acid polymerization , 1991 .

[12]  Bernard J. Matkowsky,et al.  Fronts, relaxation oscillations, and period doubling in solid fuel combustion , 1987 .

[13]  S. P. Davtyan,et al.  Two-dimensional combustion modes in condensed flow , 1992 .

[14]  S. D. Ittel,et al.  Dependence of Free-Radical Propagation Rate Constants on the Degree of Polymerization , 1996 .

[15]  Vladimir A. Volpert,et al.  Radially Symmetric Flow of Reacting Liquid with Changing Viscosity , 1993, SIAM J. Appl. Math..