Search NASA⌕ Search

Engineering topics

Faravelli, Tiziano

Publications and source records attributed to Faravelli, Tiziano.

7 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗
DOE OSTI2023

Automated identification and calculation of prompt effects in kinetic mechanisms using statistical models

The kinetics of prompt dissociation involves rovibrationally excited species (generally formed by exothermic reactions) which may dissociate or isomerize prior to thermalization via collisions with the bath gas. Treating such rovibrationally excited species (so-called "hot" species) with standard kinetic phenomenology may result in incorrect macroscopic representation of their reactivity. Here this work presents the first fully automated methodology for the calculation of prompt effects of a chosen species in a kinetic mechanism, including (i) reaction selection; (ii) theoretical calculation of rate constants and prompt branching fractions; and (iii) final rate constant fitting. The energy partition between hot fragments is estimated using a variety of statistical models, including a new physically sound microcanonical statistical model based on the rovibrational density of states of the fragments. The methodology is validated against literature data for the prompt dissociations of HCO and C 3 H 7 radicals. The microcanonical statistical model is in better agreement with trajectory simulations for larger species and is thus applicable for practical systems that typically involve large molecules, for which direct dynamics calculations are impractical. The automated workflow is applied to the evaluation of the effects of prompt dissociation for two isomeric radicals C 4 H 7 1-3 (1-methylallyl) and C 4 H 7 1-4 (3-buten-1-yl). Twelve H-atom abstraction reactions are selected and the corresponding rate constants are computed with first principles theory. The microcanonical statistical model predicts that prompt dissociations of C 4 H 7 1-3 and C 4 H 7 1-4 are already significant at 1000K, resulting in differences of up to an order of magnitude at 2000K with respect to the phenomenological thermal rate constants. To illustrate the effects of prompt dissociation on simulations of experimental data, the calculated prompt rate constants are implemented in both CRECK and C3MechV3.3 kinetic mechanisms. Simulations of experimental flame data illustrate the noticeable impact of prompt dissociation kinetics on the high-temperature combustion reactivity of C 4 H 8 -1 and C 4 H 8 -2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

HȮ 2 + HȮ 2 : High level theory and the role of singlet channels

At low temperatures and high pressures, the HȮ 2 + HȮ 2 reaction is an important reaction in combustion. There are significant unresolved discrepancies between prior theoretical and experimental studies of this reaction. It has generally been presumed to occur as an abstraction on the triplet surface to produce H 2 O 2 + 3 O 2 . Here, we employ a combination of high-level electronic structure theory (the ANL0 composite method or multi-reference methods as appropriate), sophisticated transition state theory (vibrationally adiabatic torsions, variational, and variable reaction coordinate), and master equation analyses to predict the thermal kinetics on the H 2 O 4 surface. Notably, this analysis suggests a significant branching to HȮ 3 + ȮH near 1000 K via reaction on the singlet surface. This channel does not appear to have been considered in prior combustion models. HȮ 3 itself is a metastable complex (bound by only 3 kcal mol -1 ) that rapidly dissociates to ȮH + 3 O 2 . Thus, the net reaction for this channel, HȮ 2 + HȮ 2 → ȮH + ȮH + O 2 , converts two low reactivity HȮ 2 radicals into two highly reactive ȮH radicals. The ramifications of these newly derived rate expressions are highlighted through kinetic modeling studies of H 2 , CH 3 OH, $\textit{n}$-heptane, and isooctane oxidation; all at the high pressures of relevance to combustion devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer

A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.

09 BIOMASS FUELS↗
DOE OSTI2016

A new predictive multi-zone model for HCCI engine combustion

Here, this work introduces a new predictive multi-zone model for the description of combustion in Homogeneous Charge Compression Ignition (HCCI) engines. The model exploits the existing OpenSMOKE++ computational suite to handle detailed kinetic mechanisms, providing reliable predictions of the in-cylinder auto-ignition processes. All the elements with a significant impact on the combustion performances and emissions, like turbulence, heat and mass exchanges, crevices, residual burned gases, thermal and feed stratification are taken into account. Compared to other computational approaches, this model improves the description of mixture stratification phenomena by coupling a wall heat transfer model derived from CFD application with a proper turbulence model. Furthermore, the calibration of this multi-zone model requires only three parameters, which can be derived from a non-reactive CFD simulation: these adaptive variables depend only on the engine geometry and remain fixed across a wide range of operating conditions, allowing the prediction of auto-ignition, pressure traces and pollutants. This computational framework enables the use of detail kinetic mechanisms, as well as Rate of Production Analysis (RoPA) and Sensitivity Analysis (SA) to investigate the complex chemistry involved in the auto-ignition and the pollutants formation processes. In the final sections of the paper, these capabilities are demonstrated through the comparison with experimental data.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2015

Kinetic modeling of particle size distribution of soot in a premixed burner-stabilized stagnation ethylene flame

A detailed model of soot formation is proposed, which consists of a gas-phase kinetic model for the pyrolysis and oxidation of selected hydrocarbon fuels and a kinetic mechanism of soot nucleation and mass/size growth through coagulation and surface reactions. The gas-phase model (Ranzi et al., 2012) was expanded to include the chemistry of Polycyclic Aromatic Hydrocarbons (PAHs) up to four-to-five ring PAHs, with a modular and hierarchical approach. The discrete sectional approach was employed to solve the size evolution of the particle size distribution function (PSDF). Analogy and similarity rules were employed to describe heterogeneous reaction kinetics of soot surface reactions. A variable collision efficiency was assumed for the coalescence of small soot particles. Larger particles were assumed to undergo aggregation. The predicted PSDFs are found to be in reasonably good agreement with the experimental data for nascent soot measured in an atmospheric-pressure premixed ethylene–oxygen–argon flame in the burner-stabilized stagnation flame configuration. Sensitivity analyses of the PSDF, number density, and volume fraction were carried out with respect to the rate parameters of addition reactions of acetylene, PAHs, resonantly stabilized radical reactions, and coalescence and aggregation. The findings demonstrate that the reaction of PAHs and acetylene with soot surfaces and the kinetics of coalescence and aggregation exhibit dominant effects on the detailed and global soot properties for the flame studied, in agreement with conclusions of a large range of previous modeling studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗