Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

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

Switching kinetics in ferroelectric zinc magnesium oxide thin films

The switching kinetics of RF-magnetron reactively sputtered ~200 nm thick Zn 1-x Mg x O (ZMO) ferroelectric thin films with x = 0.41 and x = 0.27 prepared on Pt/Ti/SiO 2 /Si substrates were studied at applied fields near the coercive fields, ranging from 3.4 to 5.1 MV cm -1 , and at temperatures ranging from room temperature to 100 °C. Polarization reversal in ZMO followss the Kolmogorov-Avrami-Ishibashi kinetics model for nucleation and growth at applied fields near the coercive field, and obeys an individual column switching (ICS) model at higher applied fields required to switch most of the spontaneous polarization. Switching in the high-applied field regime can be described using either the Gaussian or inverse gamma distribution functions depending on the mole fraction of magnesium in the film. The switching current transients of the high-Mg content ZMO film (x = 0.41) are always bi-modal, whereas the low Mg composition film (x = 0.27) is described by the Gaussian distribution initially following wake-up, and becomes bimodal with continued cycling. Rayleigh-like behavior of the dielectric constant revealed a simultaneous increase of the irreversible and decrease of the reversible contributions to the dielectric constant, which was ascribed to an increase in the density of mobile domain walls with cycling and resulted in faster switching.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1984

Numerical modeling of D-mappings with applications to chemical kinetics

Numerical modeling of D-mappings was studied and applied to solving nonlinear stiff systems. These mappings were locally linearized for convergence analysis, and some applications were made to chemical kinetics. The technique avoids using multistep implicit codes that require inversion of Jacobian matrices, but depends on the Jacobians for its convergence analysis.

Dey, S. K.↗
DOE OSTI2024

Precision Calcination Mechanism of CaCO 3 to High-Porosity Nanoscale CaO CO 2 Sorbent Revealed by Direct In Situ Observations

Deploying energy storage and carbon capture at scale is hindered by the substantial endothermic penalty of decomposing CaCO 3 to CaO and CO 2 , and the rapid loss of CO 2 absorption capacity by CaO sorbent particles due to sintering at the high requisite decomposition temperatures. The decomposition reaction mechanism underlying sorbent deactivation remains unclear at the atomic level and nanoscale due to past reliance on postmortem characterization methods with insufficient spatial and temporal resolution. Thus, elucidating the important CaCO 3 decomposition reaction pathway requires direct observation by time-resolved (sub-)nanoscale methods. Here, chemical and structural dynamics during the decomposition of CaCO 3 nanoparticles to nanoporous CaO particles comprising high-surface-area CaO nanocrystallites are examined. Comparing in situ transmission electron microscopy (TEM) and synchrotron X-ray diffraction experiments gives key insights into the dynamics of nanoparticle calcination, involving anisotropic CaCO 3 thermal distortion before conversion to thermally dilated energetically stable CaO crystallites. Time-resolved TEM uncovered a novel CaO formation mechanism involving heterogeneous nucleation at extended CaCO 3 defects followed by sweeping reaction front motion across the initial CaCO 3 particle. These observations clarify longstanding, yet incomplete, reaction mechanisms and kinetic models lacking accurate information about (sub-)nanoscale dynamics, while also demonstrating calcination of CaCO 3 without sintering through rapid heating and precise temperature control.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiple Roles of Alkanethiolate‐Ligands in Direct Formation of H 2 O 2 over Pd Nanoparticles

Abstract Coadsorbed organic species including thiolates can promote direct synthesis of hydrogen peroxide from H 2 and O 2 over Pd particles. Here, density functional theory based kinetic modeling, augmented with activity measurements and vibrational spectroscopy are used to provide atomistic understanding of direct H 2 O 2 formation over alkylthiolate(RS) Pd. We find that the RS species are oxidized during reaction conditions yielding RSO 2 as the effective ligand. The RSO 2 ligand shows superior ability for proton transfer to the intermediate surface species OOH, which accelerates the formation of H 2 O 2 . The ligands promote the selectivity also by blocking sites for unselective water formation and by modifying the electronic structure of Pd. The work rationalizes observations of enhanced selectivity of direct H 2 O 2 formation over ligand‐funtionalized Pd nanoparticles and shows that engineering of organic surface modifiers can be used to promote desired hydrogen transfer routes.

Chen, Lin↗
DOE OSTI2023

Hydrogen Evolution on Electrode‐Supported Pt n Clusters: Ensemble of Hydride States Governs the Size Dependent Reactivity

Abstract We report the size‐dependent activity and stability of supported Pt 1,4,7,8 for electrocatalytic hydrogen evolution reaction, and show that clusters outperform polycrystalline Pt in activity, with size‐dependent stability. To understand the size effects, we use DFT calculations to study the structural fluxionality under varying potentials. We show that the clusters can reshape under H coverage and populate an ensemble of states with diverse stoichiometry, structure, and thus reactivity. Both experiment and theory suggest that electrocatalytic species are hydridic states of the clusters (≈2 H/Pt). An ensemble‐based kinetic model reproduces the experimental activity trend and reveals the role of metastable states. The stability trend is rationalized by chemical bonding analysis. Our joint study demonstrates the potential‐ and adsorbate‐coverage‐dependent fluxionality of subnano clusters of different sizes and offers a systematic modeling strategy to tackle the complexities.

Zhang, Zisheng↗
DOE OSTI2022

Multiple Roles of Alkanethiolate‐Ligands in Direct Formation of H 2 O 2 over Pd Nanoparticles

Abstract Coadsorbed organic species including thiolates can promote direct synthesis of hydrogen peroxide from H 2 and O 2 over Pd particles. Here, density functional theory based kinetic modeling, augmented with activity measurements and vibrational spectroscopy are used to provide atomistic understanding of direct H 2 O 2 formation over alkylthiolate(RS) Pd. We find that the RS species are oxidized during reaction conditions yielding RSO 2 as the effective ligand. The RSO 2 ligand shows superior ability for proton transfer to the intermediate surface species OOH, which accelerates the formation of H 2 O 2 . The ligands promote the selectivity also by blocking sites for unselective water formation and by modifying the electronic structure of Pd. The work rationalizes observations of enhanced selectivity of direct H 2 O 2 formation over ligand‐funtionalized Pd nanoparticles and shows that engineering of organic surface modifiers can be used to promote desired hydrogen transfer routes.

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

Hydrogen Evolution on Electrode–Supported Ptn Clusters: Ensemble of Hydride States Governs the Size Dependent Reactivity

We report the size-dependent activity and stability of supported Pt1,4,7,8 for electrocatalytic hydrogen evolution reaction, and show that clusters outperform polycrystalline Pt in activity, with sizedependent stability. To understand the size effects, we use DFT calculations to study the structural fluxionality under varying potentials. We show that the clusters can reshape under H coverage and populate an ensemble of states with diverse stoichiometry, structure, and thus reactivity. Both experiment and theory suggest that electrocatalytic species are hydridic states of the clusters (~2 H/Pt). An ensemble-based kinetic model reproduces the experimental activity trend and reveals the role of metastable states. Furthermore, the stability trend is rationalized by chemical bonding analysis. Our joint study demonstrates the potential- and adsorbate-coverage-dependent fluxionality of subnano clusters of different sizes and offers a systematic modeling strategy to tackle the complexities.

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

Response surface methodology guided adsorption and recovery of free fatty acids from oil using resin

The presence of free fatty acids interferes with the conversion of plant oils to biodiesel. In this work, four strong and weak base resins were evaluated for the removal of free fatty acids (FFA) from oil. Amberlite FPA 51 showed the highest adsorption capacity of FFA. A resin concentration above 3% could enable a higher percentage FFA adsorption. The adsorption process fitted a pseudo-first-order kinetic model and achieved equilibrium in approximately 8 h. A full factorial design was used to optimize the resin and FFA concentrations at a fixed temperature (40°C). A ratio of resin to fatty acid concentrations above 1.875 was sufficient for 70% adsorption and the amount adsorbed continued to increase with further added resin. A two-step washing of resin using hexane and ethanol recovered approximately 67.55%±4.05% of the initially added fatty acid. The resin that was used was regenerated with 5% NaOH and re-used for a minimum of three consecutive cycles. However, the adsorption capacity diminished to 75% of the initial cycle in cycles 2 and 3. Thus, the work presents a resin-based process for deacidification of oil to reduce fatty acid content of oil for biodiesel production.

09 BIOMASS FUELS↗
DOE OSTI2024

Carbon‐13 Hyperpolarization of α‐Ketocarboxylates with Parahydrogen in Reversible Exchange

Abstract Signal Amplification by Reversible Exchange (SABRE) is a relatively simple and fast hyperpolarization technique that has been used to hyperpolarize the α‐ketocarboxylate pyruvate, a central metabolite and the leading hyperpolarized MRI contrast agent. In this work, we show that SABRE can readily be extended to hyperpolarize 13 C nuclei at natural abundance on many other α‐ketocarboxylates. Hyperpolarization is observed and optimized on pyruvate (P 13C =17 %) and 2‐oxobutyrate (P 13C =25 %) with alkyl chains in the R‐group, oxaloacetate (P 13C =11 %) and alpha‐ketoglutarate (P 13C =13 %) with carboxylate moieties in the R group, and phenylpyruvate (P 13C =2 %) and phenylglyoxylate (P 13C =2 %) with phenyl rings in the R‐group. New catalytically active SABRE binding motifs of the substrates to the hyperpolarization transfer catalyst – particularly for oxaloacetate – are observed. We experimentally explore the connection between temperature and exchange rates for all of these SABRE systems and develop a theoretical kinetic model, which is used to fit the hyperpolarization build‐up and decay during SABRE activity.

McBride, Stephen J. [Department of Chemistry North↗
DOE OSTI2021

High‐efficiency adsorption of various heavy metals by tea residue biochar loaded with nanoscale zero‐valent iron

Abstract Nanoscale zero‐valent iron (nZVI) loaded on tea residue biochar activated with phosphoric acid composites was synthesized. This adsorbent can efficiently adsorb multiple heavy metal ions such as Pb (II), Cr (III), Cu (II), and Ni (II) simultaneously. The adsorption process was found to follow the pseudo‐second‐order kinetic model. Equilibrium adsorption capacities reached 288.18, 162.34, 242.72, and 267.22 mg·g −1 for Pb(II), Cr(III), Cu(II), and Ni(II), respectively. The common interfering ions exhibited lower influence on its adsorption capacity. The adsorption mechanism includes reduction, ion exchange, surface complexation and precipitation. Furthermore, the waste adsorbent was converted into a catalyst in situ, which could catalyze the degradation of Congo red by establishing an electron transfer pathway. This research provides a new concept for the sustainable utilization of industrial tea residues and the harmless treatment of waste nZVI composite adsorbents.

Liu, Shuang↗
DOE OSTI2021

Green synthesis of sulfur‐doped g‐C 3 N 4 nanosheets for enhanced removal of oxytetracycline under visible‐light irradiation and reduction of its N ‐nitrosodimethylamine formation potential

Abstract BACKGROUND In this study, sulfur‐doped g‐C 3 N 4 (SCN) nanosheets were synthesized using thiourea as a precursor, and then used to catalyze the degradation of oxytetracycline (OTC) under visible‐light irradiation. OTC was chosen as target contaminant because of it being an important precursor of N ‐nitrosodimethylamine (NDMA). RESULTS The removal of OTC was influenced by its initial concentration, initial pH, SCN dosage and water matrices. OTC (≤10 mg L −1 ) was almost completely removed in 40 min under the conditions of an initial pH of 7.0 and a SCN dose of 1.0 g L −1 , and its degradation was well described by the pseudo‐first‐order kinetics model. Water matrices would also influence the degradation of OTC in the system. Electron paramagnetic resonance analysis and trapping experiments confirmed that the function of the radicals for OTC removal followed the order of • O 2 − > • OH > h + . Eight transformation products were determined during the degradation of OTC, it being removed through four pathways: methyl oxidization, demethylation, decarbonylation and secondary alcohol oxidation. SCN was stable in the system, as the removal efficiency of OTC was not significantly reduced even after three recycles. In addition, 55.5% of NDMA formation potential of OTC was reduced after SCN catalytic oxidation. CONCLUSIONS The results imply that the visible‐light‐driven SCN catalytic oxidation is an effective technology not only for the removal of OTC but also for the reduction of NDMA formation potential. © 2021 Society of Chemical Industry

Dou, Yicheng↗
DOE OSTI2022

Climatic controls on soil and saprock nitrogen distribution and persistence in the Sierra Nevada

Nitrogen (N) is an essential nutrient in soil that regulates plant growth, terrestrial sequestration of atmospheric carbon dioxide, and persistence of organic compounds. However, major knowledge gaps remain about how climate change may impact N accumulation and persistence, especially in deep soil and saprock (friable weakly weathered bedrock). Our objective was to understand how climate impacts the accumulation of N in soil and saprock and how climate impacts soil N distribution, persistence, and mechanisms of N persistence. We investigated N concentration in bulk soil and density fractions. We estimated N persistence along a bio-climatic sequence—sites range from a low-elevation oak savannah, mid-elevation pine-oak/mixed-conifer forest, to a high-elevation subalpine forest—in the southern Sierra Nevada in California. A combination of radiocarbon and elemental composition measurements along with a first-order kinetic model was used. The N concentration in the bulk soil and density fractions declined with depth, and there was a relatively greater mineral-associated heavy fraction (HF) N in deeper samples. The cooler/wetter mixed conifer site held 37% of profile N in saprock, which was greater than that of the entire soil profile at the drier/hotter oak savannah. The majority of N in soil, which was in the HF, was not influenced by climate proxies tested. However, both unprotected and occluded fractions of N were strongly influenced by climate. Soil N mean residence time (MRT) showed that drier/hotter climates have a shorter MRT, compared to mid-elevation sites with cooler/wetter climates. The effect of climate on deep saprock N storage might be indirect, primarily through climatic influence on the thickness of saprock. Overall, our findings suggest the mineral-associated HF N pool will not be vulnerable to changes in climate and will continue to contribute to the persistent soil N pool. The amount of topsoil and subsoil unprotected and occluded N can be explained by gross primary productivity and mean annual precipitation indicating that changes in climate can influence N partitioning. N stored in deep soil and saprock may be less vulnerable to climate than N stored in drier/hotter climates with less deeply stored N. It is critical to dig deeper to understand terrestrial ecosystems’ response to climate.

58 GEOSCIENCES↗
DOE OSTI2020

Isopropanol dehydration reaction rate kinetics measurement using H 2 O time histories

H 2 O formation during thermal decomposition of isopropanol was measured behind reflected shock waves at temperatures ranging from 1127 to 1621 K at an average pressure of 1.42 atm using a laser absorption technique. Of the five modern chemical kinetics models used to compare the H 2 O time histories, the model from Li et al (Combust. Flame 2019;207:171-185) showed the best overall agreement. Sensitivity and rate of production analyses using the Li et al model (as well as those from AramcoMech 3.0, CRECK, and Togbe et al (Energy Fuel. 2011;25:676-683)) showed unimolecular dehydration of isopropanol, iC 3 H 7 OH ⇌ H 2 O + C 3 H 6 (R1), is nearly the sole reaction controlling H2O production at early times, allowing for an a priori measurement of the forward rate constant k 1 . The Arrhenius expression k1 (s –1 ) = 2.60 × 10 13 exp(-31 120 K/T) was determined to represent best the data from this study. According to the models and previous experimental investigations, the pressure investigated is well within the high-pressure limit (HPL) for this reaction. Additional, higher-pressure experiments also confirmed the HPL assumption. An uncertainty analysis was performed by varying secondary reactions within their uncertainties and examining their effect on the overall prediction, establishing an uncertainty within ±20% for all but the highest temperature cases, which have a maximum uncertainty of ±40%. Experiments conducted with a radical trapper, toluene, showed little influence from radical chemistry, suggesting this estimated uncertainty is fairly conservative. Experimental data from Heyne et al (Z Phys Chem. 2015;229:881-907) were found to be in good agreement with the rate measurements from this study and, therefore, a second Arrhenius expression, k 1 (s –1 ) = 2.11 × 10 13 exp(-30 820 K/T), was found to represent both datasets well. This second expression has a larger temperature range of 976-1621 K. Here the present study provides the first high-temperature data collected for this reaction, adding to the limited data available for isopropanol in the literature.

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

High‐pressure oxidation of n ‐butane

Abstract The oxidation of n ‐butane at elevated pressures has been investigated by experiments in a laminar flow reactor at 100 bar and temperatures of 450–900 K. The onset temperature for reaction increased from 550 K under oxidizing conditions (Φ = 0.02) to 625 K under reducing conditions (Φ = 13). NTC behavior was observed at 600–650 K (Φ = 0.02) and 625–675 K (Φ = 1.0). A detailed chemical kinetic model for the oxidation of n ‐butane was established. The present model and those suggested in literature were evaluated against the present experimental results and literature data at elevated pressures. None of the tested models could accurately reproduce the NTC behavior of n ‐butane under stoichiometric conditions of the present study, but all evaluated models could reproduce experimental data from literature with different levels of accuracy.

Hashemi, Hamid↗
DOE OSTI2025

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2023

Analysis of degradation products in thermally treated TATB

Delineating the chemical composition of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) residues produced from the exposure to abnormal thermal environments should lead to a better understanding of the decomposition paths. Identifying and quantifying each compound in thermally produced residues, monitors which compounds are degrading or forming along the decomposition route, as well as providing input for the kinetic models of those pathways. Here, in this paper, we report the methodology of isolating, identifying, and where possible, quantifying soluble compounds present in solid residues of thermally treated TATB (330 °C for tens of minutes). Samples were extracted with DMSO, separated using chromatography, and quantified using their absorption at 354 nm. Identification of unknown compounds was accomplished using high resolution mass spectrometry. TATB, F1 (diamino-dinitro-benzofurazan), HO-TATB (2,4,6-triamino-1-hydroxyl-3,5-dinitrobenzene), and T4A (1-chloro-3,5-dinitro-2,4,6-triaminobenzene) were trace compounds detected in the unreacted TATB. Ten more compounds that formed in the residues were structurally identified including F2 (amino-nitro-difurazan). Several more compounds were observed but not completely identified. We propose possible structures for the unknowns. Of the compounds formed, F1 was the most abundant compound reaching 4.5 % by weight of the degraded solid sample. Other degradation compounds were estimated to sum to trace levels, well below 1 %. Most compounds were new, having not been detected and identified in previous studies of production grade and thermally aged TATB. Many compounds only reached detectable concentrations after several min of thermal exposure.

36 MATERIALS SCIENCE↗