Closed form numerical approach for a kinetic interpretation of high-cis polybutadiene rubber vulcanization with sulphur

A novel mathematical approach to predict the vulcanization degree of high-cis polybutadiene rubber vulcanized with sulphur is presented. The model has kinetic base, it is constituted by four reactions occurring in series and parallel and takes contemporarily into consideration, within a simplified but reliable scheme, the actual reactions occurring during polybutadiene sulphur curing, namely primary crosslinking and possible de-vulcanization. The first order differential equation system obtained is suitably rearranged and a closed form expression for the vulcanization degree is derived, depending the four kinetic constants characterizing the chemistry describing reactions. Instead of using classic least-squares optimization routines to characterize kinetic constants on experimental data, a simplified but reliable approach is proposed, where a system of four non-linear equations is solved with a recursive strategy, allowing estimating kinetic constants that proved to fit well normalized experimental data. The procedure is fast and its reliability is tested on a number of experimental data available, relying into a high-cis polybutadiene rubber cured under different temperatures and accelerators concentrations. Very good approximations of experimental data are obtained, also in comparison with a heuristic numerical approach where optimization is obtained interactively.     

Effective closed form starting point determination for kinetic model interpreting NR vulcanized with sulphur

In this paper, a closed form analytical approach for a recently presented kinetic model proposed in Milani and Milani (Polym Test, 2013, under review) to interpret NR sulphur vulcanization in presence of either experimental or surrogate rheometer curves is proposed. The model has kinetic base and is aimed at predicting, by means of a very refined approach, the vulcanization degree of NR vulcanized with sulphur. It needs as input only rheometer curves to fit and provides as output kinetic constants of the single reactions occurring during the crosslink process. In Milani and Milani (Polym Test, 2013, under review) a cure chemical scheme constituted by five reactions occurring in series and parallel was adopted. The chemical scheme, translated mathematically into a differential equations system, was suitably re-arranged and a single analytical equation was derived, representing rubber crosslink degree evolution upon time. The main drawback of such procedure is that the five kinetic constants corresponding to each reaction were determined through a standard non-linear least squares procedure, trying to minimize the deviation of the analytical cure curve from experimental data. Such a limitation is here superseded and a major improvement is proposed utilizing (1) a closed form solution which does not require any optimization algorithm and (2) finding analytically a starting point for the unknown kinetic constants, very near to the actual solution and thus very convenient for a successive least squares minimization. In the model, it is shown how the analytical condition deduced from the scorch point (second derivative of the rheometer curve equal to zero) and two further conditions, e.g. the time at 90 % of vulcanization and the reversion percentage, allow the simple direct evaluation of kinetic constants, providing a closed form analytical formula to predict well the state of cure of the rubber under consideration. To assess the results obtained with the model proposed, several examples on two different NRs are discussed. The approach proved to be extremely robust and much faster when compared with the model proposed by Milani and Milani (Polym Test, 2013, under review).     

Curing degree prediction for S-TBBS-DPG natural rubber by means of a simple numerical model accounting for reversion and linear interaction

The paper presents a simple numerical model able to provide directly kinetic constants and reliable numerical rheometer curves for S-TBBS-DPG natural rubber. The approach is suitable to calculate the kinetic constants and maximum torque (MH) at any S-TBBS-DPG concentration, following a 3D mathematical interpolation/extrapolation procedure, when kinetic constants on few grid points of S-TBBS-DPG concentrations are available. In particular, the possibility to estimate with sufficient accuracy the behavior of natural rubber at any intermediate concentration of S-TBBS-DPG having engineering relevance has been proved, calibrating the model by means of simple closed form standard best fitting on few experimental data. The model used is a three kinetic parameters one, derived from the well known Han's and co-workers approach, where constants have been evaluated normalizing experimental rheometers curves following the commonly accepted Sun and Isayev method. The procedure has been validated against experimentally obtained rheometer curves by means of inverse analysis, exhibiting excellent prediction capabilities. The approach may be used for practical purposes in order to avoid expensive and cumbersome experimental investigations.     

Parabola-Hyperbola P-H kinetic model for NR sulphur vulcanization

A Parabola-Hyperbola (P-H) kinetic model for NR sulphur vulcanization is presented. The idea originates from the fitting composite Parabola-Parabola-Hyperbola (P-P-H) function used by the authors in [1,2] to approximate experimental rheometer curves with the knowledge of a few key parameters of vulcanization, such as the scorch point, initial vulcanization rate, 90% of vulcanization, maximum point and reversion percentage. After proper normalization of experimental data (i.e. excluding induction and normalizing against maximum torque), the P-P-H model reduces to the discussed P-H composite function, which is linked to the kinetic scheme originally proposed by Han and co-workers [3]. Typically, it is characterized by three kinetic constants, where classically the first two describe incipient curing and stable/instable crosslinks and the last reproduces reversion.The powerfulness of the proposed approach stands into the very reduced number of input parameters required to accurately fit normalized experimental data (i.e. rate of vulcanization at scorch, vulcanization at 90%, maximum point and reversion percentage), and the translation of a mere geometric data-fitting into a kinetic model. Kinetic constants knowledge from simple geometric fitting allows characterizing rubber curing also at temperature different from those experimentally tested.The P-H model can be applied also in the so-called backward direction, i.e. assuming Han's kinetic constants known from other models and deriving the geometric fitting parameters as result.Some existing experimental data available, relying into rheometer curves conducted at 5 different temperatures on the same rubber blend are used to benchmark the P-H kinetic approach proposed, in both backward and forward direction. Very good agreement with previously presented kinetic approaches and experimental data is observed.     

A new simple numerical model based on experimental scorch curve data fitting for the interpretation of sulphur vulcanization

Sulphur was the first agent used to vulcanize commercial elastomers (e.g. natural rubber) and allows meaningful cost reductions during the industrial process (production cost ratio between peroxides and accelerated sulphur is around 5). Therefore, accelerated sulphur vulcanization is the most popular technique for the production of polydiene and EPDM elastomers items. At present, crosslinking mechanisms are not analytically known in detail, therefore reticulation kinetic has to be deduced from mechanical properties obtained during standardized tests, as for instance the oscillating disc rheometer. In the present paper, we propose a numerical model to fit experimental rheometer data based on a simple composite three functions curve, able to describe the increase of the viscosity at successive curing times at different controlled temperature to use during the production of thick items vulcanized with sulphur. It is believed that rheometer curve is able to give an indirect information on the rubber reticulation kinetic at different temperatures, to use in a successive step to establish simplified analytical kinetic formulas to adopt in the accelerated sulphur vulcanization of polydiene and EPDM elastomers. In the model, it is necessary to collect rheometer curves at different specimen temperatures, because vulcanization in industrial practice occurs at variable temperatures during curing, with considerable differences from the core to boundary of the item. Once that rheometer curves are suitably collected in a database, they are used to predict the optimal vulcanization of real items industrially produced. Finally, a so called alternating tangent approach (AT) is implemented to determine optimal input parameters (curing external temperature T _n and rubber exposition time t) to use in the production process. Output mechanical property (objective function) to optimize is represented by the average tensile strength of the item. A meaningful example of engineering interest, consisting of a thick 2D EPDM cylinder is illustrated to validate the model proposed.     

EPDM accelerated sulfur vulcanization: a kinetic model based on a genetic algorithm

A simple closed form equation for the prediction of crosslinking of EPDM during accelerated sulfur vulcanization is presented. Such a closed form solution is derived from a second order non homogeneous differential equation, deduced from a kinetic model. The kinetic model is based on the assumption that, during vulcanization, a number of partial reactions occurs, both in series and in parallel, which determine the formation of intermediate compounds, including activated and matured polymer. Once written standard first order differential equations for each partial reaction, the differential equation system so obtained is rearranged and, after few considerations, a single second order non homogeneous differential equation with constant coefficients is derived, for which a solution may be found in closed form, provided that the non-homogeneous term is approximated with an exponential function. To estimate numerically the degree of crosslinking, kinetic model constants are evaluated through a simple data fitting, performed on experimental rheometer cure curves. The fitting procedure is a new one, and is achieved using an ad-hoc genetic algorithm, provided that a few points, strictly necessary to estimate model unknown constants with sufficient accuracy, are selected from the whole experimental cure curve. To assess the results obtained with the model proposed, a number of different compounds are analyzed, for which experimental or numerical data are available from the literature. The important cases of moderate and strong reversions are also considered, experiencing a convincing convergence of the analytical model proposed. For the single cases analyzed, partial reaction kinetic constants are also provided.     

Natural rubber/styrene-butadiene rubber blends prepared by solution mixing: Influence of vulcanization temperature using a Semi-EV sulfur curing system on the microstructural properties

Blends of natural rubber (NR) and styrene-butadiene rubber (SBR) were prepared by solution mixing and vulcanized with sulfur and accelerator in a Semi-EV system at 433 K and 443 K in order to study the vulcanization kinetic and the influence of vulcanization temperature on final structure of the blends. The vulcanization kinetic studied through the variation in rheometer curves was analyzed using the Ding and Leonov model, which takes into account the reversion effect during the cure process. The average free nanohole volume and the fractional free volume of samples with different NR/SBR ratio were estimated using positron annihilation lifetime spectroscopy (PALS). Also, the crosslink density was determined by means of swelling tests in a solvent. For all the compounds, a correlation between the free nanohole volume and the delta torque obtained from the respective rheometer curves was established.     

Vulcanization kinetics study of natural rubber compounds having different formulation variables

Vulcanization kinetics of natural rubber (NR) compounds with efficient vulcanization system was studied through phenomenological approach using the experimentally cure data obtained from a moving die rheometer. The cure kinetic parameters were defined using the proposed models by Claxton?CLiska and Deng?CIsayev with the support of curve fitting software. The effects of the amount of accelerators, sulfur and silica in the formulations on the cure characteristics and cure kinetic parameters at high cure temperatures were investigated. Kinetic data results showed that the above two models were able to describe the curing behaviour of the studied compounds satisfactorily. It showed that the fitting of the experimental data with Claxton?CLiska and Deng?CIsayev could provide a good platform to investigate the cure kinetics of the prepared NR compounds.     

Closed form analytical approach for a second order non-linear ODE interpreting EPDM vulcanization with peroxides

In this paper, the recently presented kinetic model proposed in Milani and Milani (J Math Chem 51(3):1116–1133, 2013) to interpret EPDM peroxide vulcanization is extensively revised and the resultant second order ODE is solved by means of an approximate but effective closed form analytical approach. The model has kinetic base and it is aimed at predicting, by means of a very refined approach, the vulcanization degree of rubber vulcanized with peroxides. Such a procedure takes contemporarily into consideration, albeit within a simplified scheme, the actual reactions occurring during peroxidic curing, namely initiation, H-abstraction, combination and addition, and supersedes the simplified approach used in practice, which assumes for peroxidic curing a single first order reaction. The main drawback of the overall procedure proposed in Milani and Milani (J Math Chem 51(3):1116–1133, 2013) is that the single second order non-linear differential equation obtained mathematically and representing the crosslink evolution with respect to time, was solved numerically by means of a Runge–Kutta approach. Such a limitation is here superseded and a major improvement is proposed allowing the utilization of an approximate but still effective closed form solution. After some simplifications applied on some parts of the solving function not allowing direct closed form integration, an analytical function is proposed. Kinetic parameters within the analytical model are evaluated through least squares where target data are represented by few experimental normalized rheometer curve values. In order to have an insight into the reliability of the numerical approach proposed, a case of technical interest of an EPDM with low unsaturation and crosslinked with three different peroxides at three increasing temperatures is critically discussed.     

Direct and closed form analytical model for the prediction of reaction kinetic of EPDM accelerated sulphur vulcanization

A novel direct model with kinetic base for the prediction of the final vulcanization level of EPDM cured with sulphur is presented. The model bases on a preliminary characterization of rubber through standard rheometer tests and allows an accurate prediction of the crosslinking degree at both successive curing times and different controlled temperatures. Both the case of indefinite increase of the torque and reversion can be handled. The approach proposed bases on a previously presented exponential model, where a calibration of three kinetic constants at fixed temperature by means of non-linear least square fitting was required. Here the exponential model is superseded and kinetic constants are evaluated through simple closed form formulas. The applicability of the approach is immediate and makes the model extremely appealing when fast and reliable estimates of crosslinking density of cured EPDM are required. To show the capabilities of the approach proposed, a comprehensive comparison with both available experimental data and results obtained numerically with the exponential model for real compounds at different temperatures is finally provided.     

A contribution to the theory of accelerated sulphur vulcanization of natural rubber and polybutadiene BR with tetramethyl thiuram disulphide and bis(2-benzothiazolyl) disulphide

A number of simple natural rubber (NR) and polybutadiene (BR) formulations were vulcanized in a differential scanning calorimeter. From the shapes of the exothermal cure curves the interaction of curatives can be examined. It is suggested that polysulphidic-accelerator-terminated pendent groups could either crosslink, or be chelated with zinc ions. The latter reaction allows interaction with accelerator molecules to reform the unchelated pendent group, and a zinc accelerator complex, that, in turn, can react with sulphur and rubber molecules to give pendent groups. Differences in the behavior of NR and BR are attributed to the relative rates of these two reactions in the elastomers. Bis(2-benzothiazolyl) disulphide does not react with zinc, and it is suggested that it coordinates with the zinc dithiolate accelerator complex, increasing its activity.     

Effective GA approach for a direct evaluation of reaction kinetic within EPDM accelerated sulphur crosslinking

A direct genetic algorithm (GA) approach with kinetic base, to provide effective numerical estimates of vulcanization level for EPDM cross-linked with accelerated sulphur is presented. The model requires a preliminary characterization of rubber through standard rheometer tests. A recently presented kinetic exponential model is used as starting point to develop the algorithm proposed. In such a model, three kinetic constants have to be determined by means of a non-linear least-squares curve fitting. The approach proposed circumvents a sometimes inefficient and not convergent non-linear data fitting, disregarding at a first attempt reversion and finding the local minimum of a suitable two-variable error function, to have an estimate of the first two kinetic constants. A comparison between present GA approach and traditional gradient based algorithms is discussed. The last constant, representing reversion is again evaluated through a minimization performed on a single variable error function. The applicability of the approach is immediate and makes the model extremely appealing when fast and reliable estimates of crosslinking density of cured EPDM are required. To show the capabilities of the approach proposed, a comprehensive comparison with both available experimental data and results obtained numerically with a least square exponential model for a real compound at different temperatures is provided.     

Differential model accounting for reversion for EPDM vulcanized with peroxides

One of the main drawbacks of EPM/EPDM rubber vulcanization by peroxides is the lack of selectivity, which leads to a number of side reactions. The reaction mechanisms at the base of peroxides crosslinking are generally known and include the formation of alkyl and allyl (in the EPDM case) macro-radicals through H-abstraction from the macromolecular chains and the combination of these macro-radicals, which macroscopically is known with the term “vulcanization”. In the paper, a simple but effective mathematical model having kinetic base, to predict the vulcanization degree of rubber vulcanized with peroxides, is presented. The approach takes contemporarily into consideration, albeit within a simplified scheme, the actual reactions occurring during peroxidic curing, namely initiation, H-abstraction, combination and addition, and supersedes the simplified approach used in practice, which assumes for peroxidic curing a single first order reaction. After a suitable re-arrangement of the first order system of differential equations obtained from the actual kinetic system adopted, a single second order non-linear differential equation is obtained and numerically solved by means of a Runge–Kutta approach. Kinetic parameters to set are evaluated by means of a standard least squares procedure where target data are represented by experimental values available, i.e. normalized rheometer curves or percentage crosslink density experimentally evaluated by means of more sophisticated procedures. In order to have an insight into the reliability of the numerical approach proposed, two cases of technical interest are investigated in detail: the first is an EPDM crosslinked with two different peroxides, whereas the second is a compound with high level of unsaturation, showing reversion at medium-high vulcanization temperature ( $175^\circ \text{ C}$ ).     

Analysis of DSC Data Relating to Proteins Undergoing Irreversible Thermal Denaturation

New approaches to the analysis of differential scanning calorimetry (DSC) data relating to proteins undergoing irreversible thermal denaturation have been demonstrated. The experimental approaches include obtaining a set of DSC curves at various scanning rates and protein concentrations, and also reheating experiments. The mathematical methods of analysis include construction of a linear anamorphosis and simultaneous fitting of a theoretical expression for the dependence of the excess heat capacity on temperature to a set of experimental DSC curves. Different kinetic models are discussed: the one-step irreversible model, the model including two consecutive irreversible steps, the Lumry and Eyring model with a fast equilibrating first step, and the whole kinetic Lumry and Eyring model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Heterogeneous chemical kinetics: A two-dimensional insight into sigmoidα-t curve

In discussing the interrelation between the wide variety of heterogeneous chemical reactions and the constant form of conventional kinetic a(t) and â(t) experimental curves, these curves are constructed in terms of 2D Dirichlet tessellations. In mathematical respect the result is the possibility to indicate 11 singular points on each curve, providing a more detailed treatment of experimental data. In chemical respect this enables one to take the diversity of crystal structures and chemical interactions into account.     

Kinetic modelling of the thermal oxidation of polyisoprene elastomers. Part 3: Oxidation induced changes of elastic properties

This paper is the third of a series elaborating a non-empirical kinetic model for the thermal oxidation of a sulfur vulcanized polyisoprene. Here, we try to identify kinetic parameters for post-crosslinking and reversion (“decrosslinking”) from torsion measurements, under nitrogen at temperatures ranging from 100 to 160 °C. The kinetic parameters relative to oxidative reversion (selective scissions on sulfur crosslinks) are also determined. Then a system of 13 differential equations is derived from the mechanistic scheme composed of 15 elementary reactions. Diffusion and reactions are coupled in the balance equation of oxygen in order to establish the degradation thickness profiles from which it is possible to determine the modulus profiles. The latter are used, through a composite mechanics model, to predict the global torsion stiffness of a rubber barrel. The results obtained at 100, 110, 130, 140, 150 and 160 °C are in good agreement with experimental data.     

Mechanistic modeling of reversion phenomenon in sulphur cured natural rubber vulcanization kinetics

A novel kinetic model of natural rubber sulphur vulcanization is proposed. The modeling approach takes into account current knowledge on the different polysulfidic structures present during vulcanization, and the associated individual reactions. A simplified scheme is proposed, giving a mechanistic view of the reversion phenomenon, which results in a decrease of the elastic modulus (related to the sulphur crosslink density) for long vulcanization times at high temperature. The resulting set of differential equations is solved by an appropriate numerical method to predict the evolution of the degree of vulcanization for isothermal cure conditions.     

Kinetic approach to partially overlapped thermal decomposition processes

Practical usefulness of the kinetic deconvolution for partially overlapped thermal decomposition processes of solids was examined by applying to the co-precipitated basic zinc carbonate and zinc carbonate. Comparing with the experimental deconvolutions by thermoanalytical techniques and mathematical deconvolutions using different statistical fitting functions, performance of the kinetic deconvolution based on an accumulative kinetic equation for the independent processes overlapped partially was evaluated in views of the peak deconvolution and kinetic evaluation. Two-independent kinetic processes of thermal decompositions of basic zinc carbonate and zinc carbonate were successfully deconvoluted by means of the thermoanalytical measurements in flowing CO₂ and by applying sample controlled thermal analysis (SCTA). The deconvolutions by the mathematical curve fittings using different fitting functions and subsequent formal kinetic analysis provide acceptable values of the mass-loss fractions and apparent activation energies of the respective reaction processes, but the estimated kinetic model function changes depending on the fitting functions employed for the peak deconvolution. The mass-loss fractions and apparent kinetic parameters of the respective reaction processes can be optimized simultaneously by the kinetic deconvolution based on the kinetic equation through nonlinear least square analysis, where all the parameters indicated acceptable correspondences to those estimated through the experimental and mathematical deconvolutions. As long as the reaction processes overlapped are independent kinetically, the simple and rapid procedure of kinetic deconvolution is useful as a tool for characterizing the partially overlapped kinetic processes of the thermal decomposition of solids.