Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Low pressure microplasmas enabled by field ionization: Kinetic modeling

A principle of microplasma generation that utilizes field emission of electrons at the cathode and field ionization producing ions at the anode, both processes relying on nanorods or nanotubes, is explored theoretically. In this plasma generation concept, collisional ionization of atoms and molecules by electron impact would play a negligible role. Analytical estimates as well as plasma kinetic modeling by particle-in-cell method with Monte Carlo collisions in argon confirm that this principle can enable substantial plasma densities at near-collisionless microgaps, while requiring relatively low voltages, less than 100 V. An order of magnitude increase in electron number density can be achieved due to enhancement of field emission at the cathode by positive space charge at high field ionization ion current densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Probing the Effects of NO x and SO x Impurities on Oxy-Fuel Combustion in Supercritical CO 2 : Shock Tube Experiments and Chemical Kinetic Modeling

The direct-fired supercritical carbon dioxide cycles are one of the most promising power generation methods in terms of their efficiency and environmental friendliness. Two important challenges in implementing these cycles are the high pressure (300 bar) and high CO 2 dilution (>80%) in the combustor. The design and development of supercritical oxy-combustors for natural gas require accurate reaction kinetic models to predict the combustion outcomes. The presence of a small amount of impurities in natural gas and other feed streams to oxy-combustors makes these predictions even more complex. During oxy-combustion, trace amounts of nitrogen present in the oxidizer is converted to NO x and gets into the combustion chamber along with the recirculated CO 2 . Similarly, natural gas can contain a trace amount of ammonia and sulfurous impurities that get converted to NO x and SO x and get back into the combustion chamber with recirculated CO 2 . In this work, a reaction model is developed for predicting the effect of impurities such as NO x and SO x on supercritical methane combustion. The base mechanism used in this work is GRI Mech 3.0. H 2 S combustion chemistry is obtained from Bongartz et al. while NOx chemistry is from Konnov. The reaction model is then optimized for a pressure range of 30–300 bar using high-pressure shock tube data from the literature. It is then validated with data obtained from the literature for methane combustion, H 2 S oxidation, and NO x effects on ignition delay. The effect of impurities on CH 4 combustion up to 16 atm is validated using NO x -doped methane studies obtained from the literature. In order to validate the model for high-pressure conditions, experiments are conducted at the UCF shock tube facility using natural gas identical mixtures with N 2 O as an impurity at ~100 bar. Current results show that there is a significant change in ignition delay with the presence of impurities. A comparison is made with experimental data using the developed model and predictions are found to be in good agreement. The model developed was used to study the effect of impurities on CO formation from sCO 2 combustors. It was found that NO x helps in reducing CO formation while the presence of H 2 S results in the formation of more CO. The reaction mechanism developed herein can also be used as a base mechanism to develop reduced mechanisms for use in CFD simulations.

Energy & Fuels↗
DOE OSTI2018

Modified Point-Kinetics Model for Neutron Precursors and Fission Product Behavior for Fluid-Fueled Molten Salt Reactors

Fluid-fueled nuclear reactors, such as molten salt reactors (MSRs), have recently gained significant interest. These advanced reactors represent a potential revolutionary shift in the implementation of nuclear power, and as a broad class of reactors, they have the potential to directly address many U.S. energy policy objectives. Fuel that is dissolved in the coolant requires methods to account for the birth, decay, and transport of fission products not only in the core but also throughout the loop and any auxiliary systems, such as off-gas, to which liquid fuel flows, gaseous products are carried, or solid particulates plate out. System models are particularly well suited to explore the wide range of phenomena that are associated with fluid-fueled systems, especially for safeguards analysis. However, before system dynamics can be explored, the compositions of fission products of the salt throughout the loop must be determined as they drive the dynamic behavior of a reactor.This paper describes the derivation of a modified point-kinetics model for obtaining a first-order approximation of the behavior of a salt-fueled system in which neutron precursors and fission products are born in the fuel-salt and transported outside the core. This paper also provides verification of the model using a steady-state analytic solution and provides additional cases exploring the response under transient cases. This model establishes a baseline model that can be used to explore the dynamic response of fluid-fueled reactors and to investigate important safeguards issues such as mass accountability of source terms. Here, the model is implemented in the Oak Ridge National Laboratory–developed, Modelica-based TRANSFORM library that was developed to investigate various aspects of advanced energy systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Development of a kinetic model to describe six types of symbiotic interactions in a formate utilizing microalgae-bacteria cultivation system

This study investigated an algae/formate-utilizing-bacteria system that has been developed for carbon capture. This photomixotrophic consortium consumed formate to support bacterial growth so that the resulting respiration CO 2 could be simultaneously used by algae to avoid CO 2 gas-lipid mass transfer limitation. To understand biomass growth and population interactions in this unique system, a kinetic model has been developed to describe algae and bacteria multiplication, formate and nitrogen utilizations, CO 2 mass transfer, O 2 generation and consumption, lighting condition and shading effect. The simulation indicated that this ecosystem could form six types of interactions (mutualism, commensalism, parasitism, neutralism, amensalism, and competition) depending on the light intensity and nutrient availability. Furthermore, the simulation of both batch culture and chemostat, along with parameter sensitivity tests and experimental observations, offered insights into optimal applications of algae-bacteria consortium based photobiorefinery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas

X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kinetics Modeling for Design of Continuous Enzymatic Hydrolysis

Enzymatic hydrolysis of cellulose to monomeric sugars continues to be a limiting step in cost-effectively producing sugar and fermentation-based biofuels from biomass. In particular, the high cost of enzymes coupled with the long time-scale of reaction pose challenges to economic viability. Performing enzymatic hydrolysis in a continuous mode with enzyme recycle may provide a path towards substantially reduced costs for sugar production, but design and analysis are complicated by a lack of suitable kinetics models. Computationally attainable models, such as fractal-based models, require knowledge of the reaction-history of the biomass, and are thus only suitable for describing the batch reactions from which they are derived. Fundamental models, while potentially more generalizable, are often too computationally intensive to use in process or reactor modeling. In this work, a phenomenological rate model is proposed based on a two-phase substrate representation. Good agreement is seen between batch and continuous enzymatic hydrolysis (CEH) experiment data, which validates the model and enables us to solve for reactor design parameters, such as CEH reactor size and stream flow rates, based on process variables like yield. This model is integrated with techno-economic analysis software to explore economic sensitivities. Important design optimizations and tradeoffs are identified and quantified, including the relative cost imposed by rate slowdown from sugar inhibition versus the cost to remove and concentrate sugars at a lower concentration. It also identifies, high-leverage avenues for further exploration, such as increasing the maximum feasible solids concentration, and sustaining high membrane flux and reliability.

09 BIOMASS FUELS↗
DOE OSTI2023

Insights into the thermal history of carbon fibers using Raman spectroscopy and a novel kinetic model

Carbon fibers and carbon fiber composites are applied in high-performance applications, but a key consideration for application is their relative sensitivity to oxidative environments. To enable in-situ characterization of carbon fibers exposed to oxidative conditions, the Raman spectral response of T700 carbon fibers that have been exposed to a variety of dwell temperatures is reported herein with dwell times reaching up to 1 month. Here, we evaluate the spectra holistically by using integrated absolute difference analysis. By combining this analysis with straightforward kinetic models, we connect the total Raman spectral response to the temperature-time curve that could yield such a shift in spectral parameters. Our work connects the Raman spectral response of carbon fibers to their thermal history and can easily be extended to other graphitic materials, such as nuclear graphite.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Kinetic models of magnetic flux ropes observed in the Earth magnetosphere

Magnetic flux ropes (MFR) are universal magnetoplasma structures (similar to cylindrical screw pinches) formed in reconnecting current sheets. In particular, MFR with scales from about the ion inertial length to MHD range are widely observed in the Earth magnetosphere. Typical MFR have force-free configuration with the axial magnetic field peaking on the MFR axis, whereas bifurcated MFR with an off-axis peak of the axial magnetic field are observed as well. In the present paper, we develop kinetic models of force-free and bifurcated MFR and determine consistent ion and electron distribution functions. The magnetic field configuration of the force-free MFR represents well-known Gold-Hoyle MFR (uniformly twisted MFR). We show that bifurcated MFR are characterized by the presence of cold and hot current-carrying electrons. The developed models are capable to describe MFR observed in the Earth magnetotail as well as MFR recently observed by Magnetospheric Multiscale Mission at the Earth magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Autoignition Delay Time Measurements and Chemical Kinetic Modeling of Hydrogen/Ammonia/Natural Gas Mixtures

Abstract In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels’ characteristics for use in high-pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study CO formation profiles as well as the ignition delay times behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida’s shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles. Experimental results were then compared with the GRI 3.0 mechanism, as well as the proprietary UCF 2022 mechanism utilizing CHEMKIN-Pro software. In general, both models were able to capture the trend in autoignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia-hydrogen mixtures, GRI 3.0 failed to predict ignition delay times whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Finally, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.

Baker, Jessica B.↗
DOE OSTI2022

Experimental Ignition Delay Time Measurements and Chemical Kinetics Modeling of Hydrogen/Ammonia/Natural Gas Fuels

In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels' characteristics for use in high pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study carbon monoxide (CO) formation profiles as well as the ignition delay times (IDTs) behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida's shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles over a temperature range between 1200 K and 1800 K with an average pressure of 2.2 atm. Experimental mixtures included both neat and combination natural gas/hydrogen with ammonia addition, with all mixtures except one having an equivalence ratio of 1. Results were then compared with the GRI 3.0 mechanism, as well as the newly developed UCF 2022 mechanism utilizing chemkin-pro software. In general, both models were able to capture the trend in auto-ignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia–hydrogen mixtures, GRI 3.0 failed to predict ignition delay times, whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Lastly, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

VALIDATION OF A RELAP5 POINT KINETICS MODEL OF TREAT

Development of advanced nuclear fuels requires testing of fuel specimens under various power-cooling mismatch conditions. The Transient Reactor Test (TREAT) facility at the Idaho National Laboratory will be used to provide transient nuclear heating of fuel specimens. Experiments will be conducted to support assessment of accident tolerant fuel designs being proposed for light water reactor applications. The experiments to be conducted at TREAT will provide data for the assessment of advanced, multi-physics computer codes. However, before the experiments can be performed, calculations are required to demonstrate that they will meet program objectives and that they can be performed safely in TREAT in the event of an unplanned reactivity insertion accident that occurs during the experiment. The safety calculations for the first experiments need to be performed in the next year and therefore will be performed with an existing computer code, RELAP5-3D. A point kinetics model of TREAT was developed for RELAP5-3D. The purpose of the model is to calculate the reactor power during experiments and to support the accident analyses of those experiments. The RELAP5-3D model was validated using data from nine TREAT experiments conducted previously. The RELAP5-3D model generates results that are in reasonable agreement with measured values of maximum core power, energy deposition, and maximum fuel temperature for a wide range of reactivity insertions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Elucidating reactivity regimes in cyclopentane oxidation: Jet stirred reactor experiments, computational chemistry, and kinetic modeling

This study is concerned with the identification and quantification of species generated during the combustion of cyclopentane in a jet stirred reactor (JSR). Experiments were carried out for temperatures between 740 and 1250 K, equivalence ratios from 0.5 to 3.0, and at an operating pressure of 10 atm. The fuel concentration was kept at 0.1% and the residence time of the fuel/O 2 /N 2 mixture was maintained at 0.7 s. The reactant, product, and intermediate species concentration profiles were measured using gas chromatography and Fourier transform infrared spectroscopy. The concentration profiles of cyclopentane indicate inhibition of reactivity between 850-1000 K for φ=2.0 and φ=3.0. This behavior is interesting, as it has not been observed previously for other fuel molecules, cyclic or non-cyclic. A kinetic model including both low- and high-temperature reaction pathways was developed and used to simulate the JSR experiments. The pressure-dependent rate coefficients of all relevant reactions lying on the PES of cyclopentyl + O 2 , as well as the C-C and C-H scission reactions of the cyclopentyl radical were calculated at the UCCSD(T)-F12b/cc-pVTZ-F12//M06-2X/6-311++G(d,p) level of theory. The simulations reproduced the unique reactivity trend of cyclopentane and the measured concentration profiles of intermediate and product species. Furthermore, sensitivity and reaction path analyses indicate that this reactivity trend may be attributed to differences in the reactivity of allyl radical at different conditions, and it is highly sensitive to the C-C/C-H scission branching ratio of the cyclopentyl radical decomposition.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Shock‐tube spectroscopic water measurements and detailed kinetics modeling of 1‐pentene and 3‐methyl‐1‐butene

Abstract To understand the effects of the chemical structure of two C 5 alkene isomers on their combustion properties, and to highlight the major chemical reactions occurring during their high‐temperature oxidation, water time histories were measured behind reflected shock waves for the oxidation of 1‐pentene (C 5 H 10 ‐1) and 3‐methyl‐1‐butene (3M1B) in 99.5% Ar. The experiments were carried out at three different equivalence ratios ( φ = 0.5, 1.0, and 2.0) at pressures and temperatures ranging from 1.29 to 1.47 atm and 1 331 to 1 877 K, respectively. The H 2 O quantification extends the database for 1‐pentene and provides new insights for 3M1B. These unique results were used to validate and to develop a new detailed kinetics model. Numerical predictions are presented, and the new model was able to capture the results with suitable accuracy, with 3M1B being notably more reactive than C 5 H 10 ‐1. Sensitivity and rate‐of‐production analyses were performed to help explain the results. Under the present conditions, the reactivity is rapidly initiated by molecular dissociation of a fraction of the pentene isomers. The initiation phase then induces H‐atom abstraction by active radicals (H, OH, O, HO 2 , and CH 3 ) to first produce alkenyl C 5 H 9 radicals (or an alkyl radical and an alkenyl radical by breaking a C─C bond) and subsequent, smaller fragments. The difference in terms of reactivity between the isomers is essentially due to the fact that 3M1B has one particularly weak tertiary allylic C─H bond, which allows for fast H‐atom abstraction compared with 1‐pentene.

Grégoire, Claire M.↗
DOE OSTI2016

Compositional effects on the ignition of FACE gasolines [Compositional effects on the ignition of FACE gasoline fuels: experiments, surrogate fuel formulation, and chemical kinetic modeling]

As regulatory measures for improved fuel economy and decreased emissions are pushing gasoline engine combustion technologies towards extreme conditions (i.e., boosted and intercooled intake with exhaust gas recirculation), fuel ignition characteristics become increasingly important for enabling stable operation. Here, this study explores the effects of chemical composition on the fundamental ignition behavior of gasoline fuels. Two well-characterized, high-octane, non-oxygenated FACE (Fuels for Advanced Combustion Engines) gasolines, FACE F and FACE G, having similar antiknock indices but different octane sensitivities and chemical compositions are studied. Ignition experiments were conducted in shock tubes and a rapid compression machine (RCM) at nominal pressures of 20 and 40 atm, equivalence ratios of 0.5 and 1.0, and temperatures ranging from 650 to 1270 K. Results at temperatures above 900 K indicate that ignition delay time is similar for these fuels. However, RCM measurements below 900 K demonstrate a stronger negative temperature coefficient behavior for FACE F gasoline having lower octane sensitivity. In addition, RCM pressure profiles under two-stage ignition conditions illustrate that the magnitude of low-temperature heat release (LTHR) increases with decreasing fuel octane sensitivity. However, intermediate-temperature heat release is shown to increase as fuel octane sensitivity increases. Various surrogate fuel mixtures were formulated to conduct chemical kinetic modeling, and complex multicomponent surrogate mixtures were shown to reproduce experimentally observed trends better than simpler two- and three-component mixtures composed of n-heptane, iso-octane, and toluene. Measurements in a Cooperative Fuels Research (CFR) engine demonstrated that the multicomponent surrogates accurately captured the antiknock quality of the FACE gasolines. Simulations were performed using multicomponent surrogates for FACE F and G to reveal the underlying chemical kinetics linking fuel composition with ignition characteristics. Finally, a key discovery of this work is the kinetic coupling between aromatics and naphthenes, which affects the radical pool population and thereby controls ignition.

02 PETROLEUM↗
DOE OSTI2016

WE-H-BRA-09: Application of a Modified Microdosimetric-Kinetic Model to Analyze Relative Biological Effectiveness of Ions Relevant to Light Ion Therapy Using the Particle Heavy Ion Transport System

Purpose: To evaluate the dose and biological effectiveness of various ions that could potentially be used for actively scanned particle therapy. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary particles for 1H beams and ten million particles for 4He, 7Li, 10B, 12C, 14N, 16O, and 20Ne were simulated for 0.6cm diameter pencil beams. Beam energies corresponding to Bragg peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially in annuli with outer radius of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model for five different cell types to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. The product of the calculated RBEs and the simulated physical dose was taken to create biological dose and comparisons were then made between the various ions. Results: Transversely, the 10B beam was seen to minimize relative biological dose in both the constant and accelerated dose change regions, proximal to the Bragg Peak, for all beams traveling greater than 50mm. For the 50mm beam, 7Li was seen to provide the most optimal biological dose profile. Radially small fluctuations (<4.2%) were seen in RBE while physical dose was greater than 1% for all beams. Conclusion: Even with the growing usage of 12C, it may not be the most optimal ion in all clinical situations. Boron was calculated to have slightly enhanced RBE characteristics, leading to lower relative biological doses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Severity factor kinetic model as a strategic parameter of hydrothermal processing (steam explosion and liquid hot water) for biomass fractionation under biorefinery concept

Hydrothermal processes are an attractive clean technology and cost-effective engineering platform for bio-refineries based in the conversion of biomass to biofuels and high-value bioproducts under the basis of sustainability and circular bioeconomy. The deep and detailed knowledge of the structural changes by the severity of biomasses hydrothermal fractionation is scientifically and technological needed in order to improve processes effectiveness, reactors designs, and industrial application of the multi-scale target compounds obtained by stream explosion and liquid hot water systems. The concept of the severity factor [log 10 (R 0 )] established > 30 years ago, continue to be a useful index that can provide a simple description of the relationship between the operational conditions for biomass fractionation in second generation of biorefineries. Furthermore, this review develops a deep explanation of the hydrothermal severity factor based in lignocellulosic biomass fractionation with emphasis in research advances, pretreatment operations and the applications of severity factor kinetic model.

09 BIOMASS FUELS↗
DOE OSTI2017

A kinetic modeling study on octane rating and fuel sensitivity under HCCI conditions

With the development of advanced compression ignition engines such as the homogeneous charge compression ignition (HCCI) and the reactivity controlled compression ignition (RCCI), the long existing research and motor octane rating for fuels, i.e., RON and MON, might be insufficient to describe the fuel reactivity and performance with these new combustion strategies due to their distinctive engine operation conditions. It is also noted that with the stringent regulation for CO2 emission, a detailed knowledge of fuel reactivity is urgently needed for the selection of novel alternative bio-fuels for transportation, and optimization of their engine performance in the future. It is therefore very important to understand the reactivity and sensitivity for general classes of fuels as a function of the operating conditions in these advanced combustion engines. In this modeling study, using validated kinetic mechanism for ignition, we have simulated the performance of different classes of fuels including hydrogen, alkanes (straight and branched), alcohols and potential bio-fuels under typical HCCI engine conditions using detailed chemical kinetics. The crank angle corresponding to 50% total heat release (CA50) is utilized as an indicator for the fuel reactivity. Iso-contours of CA50 are then plotted in the engine operating parameter coordinate, such as the intake temperature, and pressure. It is shown that CA50 value can vary significantly with different engine operating conditions, and that fuels with the same RON/MON even exhibit qualitatively different CA50 iso-contours in specific operating regimes, implying differences in the intrinsic reactivity of different classes of fuels, as well as the insufficiency of RON/MON as comprehensive metrics of fuel reactivity in advanced compression ignition (ACI) strategies. In addition, the operating regime with high and low sensitivity, as well as potential non-monotonic change in fuel reactivity, has been identified by the CA50 iso-contour map for all studied fuels.

Wu, Tong↗
DOE OSTI2021

Kinetic Monte Carlo Tool for Kinetic Modeling of Linear Step-Growth Polymerization: Insight into Recycling of Polyurethanes

A kinetic Monte Carlo model of polyurethane polymerization which explicitly tracks the polymer sequences is developed and shared. Here this model is benchmarked against theoretical and experimental polyurethane data and used to investigate the effect on oligomer distributions of unequal reactivity of the first and second isocyanate to react. The reverse reactions using thermodynamic consistency are then added to the framework, and analogous to the addition polymerization concept of ceiling temperature, equilibrium chain length distributions at various temperatures are calculated. For a mixture of three monomers AA, BB, and CC, where BB and CC do not react with one another, are present in stoichiometric proportions, and have different enthalpies of reaction with AA, an odd-even effect emerges. Odd length chains are more likely than even length chains for temperatures at which BB and CC have significantly different equilibrium conversions. The concept of ceiling temperature that is typically cited for addition polymers is extended here to provide a measure of conditions under which depolymerization for recycling is favored.

condensation polymerization↗