The kinetic modeling of the pyrolysis and combustion of liquid transportation fuels is a very complex task for two different reasons: the challenging characterization of the complex mixture of several hydrocarbon isomers and the complexity of the oxidation mechanisms of large hydrocarbon and oxygenated molecules. While surrogate mixtures of reference components allow to tackle the first difficulty, the complex behavior of the oxidation mechanisms is mostly overcome by reducing the total number of involved species by adopting a lumping approach. After a first investigation of the different liquid fuels (gasoline, kerosene, and diesel fuels), a short discussion on the lumping techniques allows to highlight the advantages of this approach. The lumped POLIMI pyrolysis and oxidation mechanism of hydrocarbon and oxygenated fuels is then used for generating several skeletal mechanisms for typical surrogate mixtures, moving from pure n-heptane up to heavy diesel fuels. These skeletal models are simply reduced with a reaction flux analysis, and they involve between 100 and 200 species. While these sizes already allow detailed computational fluid dynamics (CFD) calculations in internal combustion engines, further reduction phases are necessary when the interest is toward more complex CFD computations. To maintain the standard structure of the skeletal mechanisms, successive reduction phases are not considered. Moreover, new regulations pushed toward a greater use of renewable fuels. For these reasons, the skeletal models are also extended to biogasolines including methanol, ethanol, and n-butanol. Similarly, skeletal models of diesel and biodiesel fuels, including methyl esters, are also provided. Several comparisons with experimental data and complete validations in the operating range of internal combustion engines are also reported. The whole set of comparisons with experimental data obtained in a wide range of conditions not only validate the reduced models of specific transportation fuels but also the complete kinetic scheme POLIMI 1404

No abstract

No abstract

The low- and high-temperature oxidation mechanisms of n-heptane have been extensively studied in recent and past literature because of its importance as a primary reference fuel. Recent advanced analytical methods allowed for the identification of several intermediate oxygenated species at very low-temperature conditions in jet-stirred reactors. On these bases, new classes of successive reactions involving hydroperoxide species, already discussed for propane and n-butane oxidation, were included in the low-temperature oxidation mechanism of n-heptane. These new reactions allowed for the improvement of the overall mechanism, not only obtaining a satisfactorily agreement with reaction products, such as organic acids, diones, and ketones, but also in terms of system reactivity. Moreover, deeper attention was also paid to the formation of ketohydroperoxides, rarely experimentally measured. Because of n-heptane importance as a primary reference fuel, the overall POLIMI kinetic mechanism is validated in a wide range of conditions, in both the high- and low-temperature regimes. Moreover, the reliability of the updated oxidation mechanism is further proven in a couple of more complex applications, such as the autoignition of nheptane droplets in microgravity conditions and the oxidation of lean n-heptane/air mixtures in homogeneous charge compression ignition (HCCI) engines

No abstract

Alkylated aromatics constitute a significant fraction of the components commonly found in commercial fuels. Toluene is typically considered as a reference fuel. Together with n-heptane and iso-octane, it allows for realistic emulations of the behavior of real fuels by the means of surrogate mixture formulations. Moreover, it is a key precursor for the formation of poly-aromatic hydrocarbons, which are of relevance to understanding soot growth and oxidation mechanisms. In this study the POLIMI kinetic model is first updated based on the literature and on recent kinetic modelling studies of toluene pyrolysis and oxidation. Then, important reaction pathways are investigated by means of high-level theoretical methods, thereby advancing the present knowledge on toluene oxidation. H-Abstraction reactions by OH, HO2, O and O2, and the reactivity on the multi well benzyl-oxygen (C6H5CH2 + O2) potential energy surface (PES) were investigated using electronic structure calculations, transition state theory in its conventional, variational, and variable reaction coordinate forms (VRC-TST), and master equation calculations. Exploration of the effect on POLIMI model performance of literature rate constants and of the present calculations provides valuable guidelines for implementation of the new rate parameters in existing toluene kinetic models.

A priori rate predictions for gas phase reactions have undergone a gradual but dramatic transformation, with current predictions often rivaling the accuracy of the best available experimental data. The utility of such kinetic predictions would be greatly magnified if they could more readily be implemented for large numbers of systems. Here, we report the development of a new computational environment, namely, EStokTP, that reduces the human effort involved in the rate prediction for single channel reactions essentially to the specification of the methodology to be employed. The code can also be used to obtain all the necessary master equation building blocks for more complex reactions. In general, the prediction of rate constants involves two steps, with the first consisting of a set of electronic structure calculations and the second in the application of some form of kinetic solver, such as a transition state theory (TST)-based master equation solver. EStokTP provides a fully integrated treatment of both steps through calls to external codes to perform first the electronic structure and then the master equation calculations. It focuses on generating, extracting, and organizing the necessary structural properties from a sequence of calls to electronic structure codes, with robust automatic failure recovery options to limit human intervention. The code implements one or multidimensional hindered rotor treatments of internal torsional modes (with automated projection from the Hessian and with optional vibrationally adiabatic corrections), Eckart and multidimensional tunneling models (such as small curvature theory), and variational treatments (based on intrinsic reaction coordinate following). This focus on a robust implementation of high-level TST methods allows the code to be used in high accuracy studies of large sets of reactions, as illustrated here through sample studies of a few hundred reactions. At present, the following reaction types are implemented in EStokTP: abstraction, addition, isomerization, and beta-decomposition. Preliminary protocols for treating barrierless reactions and multiple-well and/ or multiple-channel potential energy surfaces are also illustrated.

No abstract