Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion dynamics”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Computational Study of Mixing Solid Materials for CO2 Capture Technology

CO2 is one of the major combustion products that, once released into the air, can contribute to global climate change. Solid sorbents have been reported as promising candidates for CO2 sorbent applications due to their high CO2 absorption capacities at moderate working temperatures. By combining thermodynamic database mining with first-principles density functional theory and phonon lattice dynamics calculations, the National Energy Technology Laboratory has proposed and validated a theoretical screening methodology to identify the most promising CO2 sorbent candidates from a vast array of possible solid materials. The advantage of this method is that it identifies the thermodynamic properties of the CO2 capture reaction as a function of temperature and gas pressure without any experimental input beyond crystallographic structural information of the solid phases involved. The calculated thermodynamic properties of solid materials versus temperature and pressure changes were further used to evaluate the equilibrium properties for the CO2 adsorption/desorption cycles. The selected CO2 sorbent candidates were further considered for experimental validations. In this presentation, the results of M2O/Al2O3 (M=Li, Na, K) mixtures (MAlO2 and M5AlO4) capturing CO2 will be demonstrated in detail.

Duan, Yuhua↗
DOE OSTI2021

Dynamic life-cycle carbon analysis for fast pyrolysis biofuel produced from pine residues: implications of carbon temporal effects

Abstract Background Woody biomass has been considered as a promising feedstock for biofuel production via thermochemical conversion technologies such as fast pyrolysis. Extensive Life Cycle Assessment studies have been completed to evaluate the carbon intensity of woody biomass-derived biofuels via fast pyrolysis. However, most studies assumed that woody biomass such as forest residues is a carbon–neutral feedstock like annual crops, despite a distinctive timeframe it takes to grow woody biomass. Besides, few studies have investigated the impacts of forest dynamics and the temporal effects of carbon on the overall carbon intensity of woody-derived biofuels. This study addressed such gaps by developing a life-cycle carbon analysis framework integrating dynamic modeling for forest and biorefinery systems with a time-based discounted Global Warming Potential (GWP) method developed in this work. The framework analyzed dynamic carbon and energy flows of a supply chain for biofuel production from pine residues via fast pyrolysis. Results The mean carbon intensity of biofuel given by Monte Carlo simulation across three pine growth cases ranges from 40.8–41.2 g CO 2 e MJ −1 (static method) to 51.0–65.2 g CO 2 e MJ −1 (using the time-based discounted GWP method) when combusting biochar for energy recovery. If biochar is utilized as soil amendment, the carbon intensity reduces to 19.0–19.7 g CO 2 e MJ −1 (static method) and 29.6–43.4 g CO 2 e MJ −1 in the time-based method. Forest growth and yields (controlled by forest management strategies) show more significant impacts on biofuel carbon intensity when the temporal effect of carbon is taken into consideration. Variation in forest operations and management (e.g., energy consumption of thinning and harvesting), on the other hand, has little impact on the biofuel carbon intensity. Conclusions The carbon temporal effect, particularly the time lag of carbon sequestration during pine growth, has direct impacts on the carbon intensity of biofuels produced from pine residues from a stand-level pine growth and management point of view. The carbon implications are also significantly impacted by the assumptions of biochar end-of-life cases and forest management strategies.

09 BIOMASS FUELS↗
DOE OSTI2021

Molecular design and fabrication of PIM-1/polyphosphazene blend membranes with high performance for CO 2 /N 2 separation

New polymeric blend membranes for CO 2 separation were synthesized based on insights from molecular dynamics simulation. A molecular-level structure-property relationship in polymers of intrinsic microporosity (PIM) based blend membranes, was investigated in detail computationally. Calculated local density profiles and energy of interaction of the blend membranes, composed of PIM-1 and various polyphosphazenes, showed that using the polyphosphazene with a higher concentration of ether side chains can improve the compatibility with PIM-1. Furthermore, based on the findings of computational studies, blend membranes were experimentally fabricated from PIM-1 and polyphosphazenes with various polyether side chain concentrations. Polyether concentration in polyphosphazenes was correlated with the film properties and gas transport performance of the blend membranes. Blend membranes showed very high CO 2 permeability (3100-5300 barrer) and improved CO 2 /N 2 selectivity (24-28), outperforming all other PIM-based blend membranes reported to date. Moreover, the CO 2 permeability performance of the blend membranes was tested 566 hours under real post-combustion flue gas from a coal-fired power plant, including CO 2 , N 2 , H 2 O, O 2 , SO x and NO x .

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

In-Situ Species Concentration Measurements In Ammonia-Mix Flames Using Ftir Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH 3 /H 2 ) and ammonia-natural gas (NH 3 /CH 4 ) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H 2 O absorption.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of Premixed and Partially Premixed Jet-in-Crossflow Flames at High Pressure

In this paper, we explore the operational map of a lean axial-staged combustor of premixed and partially premixed reacting jet-in-crossflow flames at high -pressure (5 atm). This study attempts to expand the data to relatively high pressure and could significantly aid scaling to real gas turbine engine conditions at 20–30 atm. High-speed camera, particle image velocimetry (PIV), CH* chemiluminescence, temperature, and pressure measurements were taken and processed to allow accurate reconstruction of six operating points relative to computational fluid dynamics (CFD) simulations under minimal adjustments. Variation of lean main stage (φ = 0.575 and 0.73) and rich jet (φ = 1.1, 4, and 8) equivalence ratio has been investigated for a four mm axial jet. The fully premixed flames were found to be controlled by the crossflow temperature before ignition and the crossflow oxygen content during combustion. Analysis of flame shape and position for the partially premixed operating points describes a lee stabilized as well as a more unsteady windward flame branch. Adjustment of added jet fuel and crossflow temperature along with its corresponding oxygen level is required to attain a compact flame body. Here, the risk of delaying combustion progress is significantly increased at a richer jet φ = 8 and an overshooting, spatially divided flame was attained with a main stage φ = 0.73. Control toward a compact flame body is critical to allow combustion at reasonable reaction rate.

03 NATURAL GAS↗
DOE OSTI2023

Laser heating and evaporation of a single droplet

The laser technology is being abundantly studied for controlled energy deposition for a range of applications in aerodynamic flow control, material processing, ignition, and combustion. The absorption of laser radiation by liquid droplets affects further propagation of laser in the atmosphere and causes bleaching of suspended droplets while the ignition and combustion characteristics in combustors are influenced by the evaporation rate of the sprayed fuel. In this work, we present a multi-dimensional mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion dominated regime taking into account absorption of the laser radiation, evaporation process and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for the ethanol and water droplet heating and evaporation. For continuous heating the peak temperature is established by the balance of cooling, evaporation and heating and results in high temperature for larger droplets. It has been shown that for heating by a single laser pulse the maximum temperature of droplets depends only on the peak intensity of the laser radiation. Furthermore, for the peak irradiance close to the transition to the boiling regime, temporal dynamics of the droplet temperature is independent of the droplet size. With proper normalization of time, the dynamics of the droplet shrinkage and cooling is shown to be independent of droplet sizes and peak laser intensities. The influence of cooling and evaporation processes on droplet heating was found to be controlled by the pulse repetition rate for repeated pulse operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An experimental and computational analysis of combustion heat release transformation in dual fuel combustion

Dual fuel (DF) diesel-methane combustion, which employs a high-reactivity fuel (diesel) to ignite a low-reactivity fuel (methane), is a widely studied combustion strategy for internal combustion engines, with significant potential for engine-out emissions reductions without the need for major hardware modifications. A phenomenon, which has been reported in the DF literature, but not explained fully, is the transformation of the shape of the apparent heat release rate (AHRR) curve as the start of injection (SOI) of diesel is advanced beyond a certain threshold; coincidentally, this AHRR transformation is usually accompanied by a sharp decrease in engine-out emissions of oxides of nitrogen (NOx). The goal of the present work is to establish the underlying physical reason(s) that cause the AHRR transformation. The AHRR transformation was observed on a single cylinder research engine (SCRE) at an indicated mean effective pressure (IMEP) of 5 bar at a speed of 1500 rev/min. The transformation occurred over a range of SOIs from 330 to 320 crank angle degrees (CAD). While the 330 CAD SOI exhibited a typical two-stage AHRR curve, with a clearly definable first-stage peak followed by a second-stage AHRR with little-to-no low temperature heat release (LTHR) present and high engine-out NOx, the 320 CAD SOI exhibited a single-stage, Gaussian-like AHRR curve, with noticeable LTHR and at least one order-of-magnitude lower NOx emissions. Here, leveraging analysis of experimental data and three-dimensional computational fluid dynamic simulations, the authors show that the AHRR transformation is impacted mainly by differences in local equivalence ratio distributions within the cylinder at ignition onset for different diesel SOIs.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Advances in imaging of chemically reacting flows

Many important chemically reacting systems are inherently multi-dimensional with spatial and temporal variations in the thermochemical state, which can be strongly coupled to interactions with transport processes. Fundamental insights into these systems require multi-dimensional measurements of the thermochemical state as well as fluid dynamics quantities. Laser-based imaging diagnostics provide spatially and temporally resolved measurements that help address this need. The state of the art in imaging diagnostics is continually progressing with the goal of attaining simultaneous multi-parameter measurements that capture transient processes, particularly those that lead to stochastic events, such as localized extinction in turbulent combustion. Development efforts in imaging diagnostics benefit from advances in laser and detector technology. This article provides a perspective on the progression of increasing dimensionality of laser-based imaging diagnostics and highlights the evolution from single-point measurements to 1D and 2D multi-parameter imaging and 3D high-speed imaging. This evolution is demonstrated using highlights of laser-based imaging techniques in combustion science research as an exemplar of a complex multi-dimensional chemically reacting system with chemistry–transport coupling. Imaging diagnostics impact basic research in other chemically reacting systems as well, such as measurements of near-surface gases in heterogeneous catalysis. The expanding dimensionality of imaging diagnostics leads to larger and more complex datasets that require increasingly demanding approaches to data analysis and provide opportunities for increased collaboration between experimental and computational researchers in tackling these challenges.

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

DNS of Hydrodynamic Instabilities of Laminar H2/O2/N2 Flames at Elevated Pressure

The atypical combustion properties of hydrogen have long provided an interest for hydrogen-enriched combustion, emphasized in recent years by increasingly stringent emission regulations. Significant work has been focused on the small hydrocarbon(C1-C3)/hydrogen blends, but much less is known about the effect of hydrogen on heavier hydrocarbon combustion, of interest for the aviation and automotive industry. In this work, we perform direct numerical simulations (DNS) of spherically expanding laminar n-dodecane/H2 flames in a constant-volume vessel at elevated pressure and a range of hydrogen seeding levels. A low Mach number model is used to represent the gas and flame dynamics allowing for a temporally varying, spatially homogeneous pressure field, while also enabling numerical evolution of the system numerically at time step sizes governed by the advective CFL limit rather than acoustic processes. We use an adaptive mesh refinement (AMR) approach to tackle the large separation of scales associated with high pressure premixed flame surfaces propagating in a laboratory-scale closed domain. We incorporate detailed chemistry and transport models for the reacting flow and focus the analysis on the effects of the pressure rise on the flame propagation characteristics and morphology as the hydrogen content is increased. The presence of fast diffusive H2 and the pressure rise results in dramatic thinning of the flame and the potential onset of hydrodyanmic and theromdiffusive instabilities. We assess the potential role of these instabilities on the ability to accurately measure laminar burning speed experimentally based on the mean flame surface propagation speed.

high pressure↗
DOE OSTI2022

Plasma-assisted chemical-looping combustion: Mechanistic insights into low temperature methane oxidation with CuO

The low-temperature oxidation of CH 4 by CuO in a coaxial, fixed bed, double dielectric barrier discharge (DBD) reactor was investigated with time-dependent species measurements by an electron-ionization molecular beam mass spectrometer (EI-MBMS). In the experiment, 10% methane carried by noble gasses was flown at 50 sccm through 1 g CuO dispersed in quartz wool both under plasma and non-plasma conditions, while time-dependent gas-phase species profiles were collected. Plasma conditions were explored from 300 to 600 °C while the non-plasma conditions were set from 300 to 900 °C. Mechanistic insights into the oxidation of CH 4 by CuO with plasma discharge at lower temperatures (≤ 600 °C) were obtained by quantifying the fuel oxidation, intermediate species, and CO 2 production in comparison to the non-plasma conditions. Here, we observed significant enhancement of fuel oxidation from the plasma discharge between 400 and 500 °C. The CO 2 production at 500 °C with plasma was greater than that at 700 °C without plasma, reducing fuel oxidation temperature by 200+ °C. During tests, three distinct reaction stages were observed: a gas-phase transport limited stage, a surface reaction limited kinetic stage, and an oxygen ion diffusion limited stage. It was observed that plasma greatly improved the reactivity of the second stage at low temperature. In addition, no carbon deposits were observed on the resultant particles, even under the presence of plasma. M. species such as C 4 H 2 and C 6 H 6 not previously observed or predicted in CuO/CH 4 chemical looping were observed, with some species such as CH 3 OH only becoming detectable as total flowrate was increased from 50 to 1500sccm. A non-plasma reaction pathway for CH 4 based the observed species from the MBMS spectrum and previous predictions from reactive molecular dynamics simulations was created, providing a framework from which more complex plasma CuO mechanisms can be crafted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Advanced Structured Adsorbent Architectures for Transformative Carbon Dioxide Capture Performance (Final Report)

Svante is a world leader at using solid sorbents for low-cost Carbon Dioxide (CO 2 ) capture, a technology which is recognized as critical in meeting the dual mandates of energy security/reliability and the mitigation of man-made CO 2 emissions. Svante has been developing proprietary adsorbent material compositions, forming them into structured laminates, developing and optimizing process cycles, and system design for efficient capture of CO 2 from post-combustion flue gases of thermal power plants and industrial facilities. The deployment of first-generation CO 2 capture technology has been significantly hampered by high costs and energy penalties, among other barriers. Second generation CO 2 capture technologies (including the Mark I variant of Svante’s Veloxotherm™ adsorption-based technology), utilizing single adsorbent architecture, show promise for reducing the barriers to deploying CO 2 capture plants in commercially meaningful numbers. The objective of this project was to evaluate the Recipient’s transformational (Mark-II) VeloxoTherm™ Technology via the development and bench-scale testing of an advanced structured adsorbent, including novel Bi-layer, laminated adsorbent structures and segmented beds. Svante selected, synthesized, and characterized tailored solid adsorbents for computational modeling, advanced structured adsorbent development, process simulations, and dynamic bench scale (~1-10 kg/day CO 2 captured) testing using an existing single-bed VeloxoTherm™ Station (VTS) coupled with a natural gas-fired boiler. Segmented beds used the in-house, multi-bed Process Demonstration Unit (PDU) to demonstrate key performance indicators (KPIs), such as recovery, product purity, regeneration energy, and the integrated system's productivity in lifetime analysis. Segmented beds were used at a 1 tonne per day (TPD) unit at an industrial site to provide bench-scale validation of performance in an industrial setting. Svante was developing and optimizing the post-combustion CO 2 adsorption technology architectures, including the Bi-layer and segmented laminated adsorbent structure design, integrated rapid cycle temperature swing adsorption (RC-TSA) cycle, flow path architecture, and adsorbent bed construction and packaging (including gas porting) to progress towards achievement of DOE’s Transformational CO 2 Capture goals of 95% CO 2 purity and a cost of electricity at least 30% lower than a supercritical Pulverized Coal (PC) power plant with CO 2 capture, or approximately $30 per tonne of CO 2 captured ready for demonstration by 2030. The main requirements to reach the DOE target cost of carbon capture below $30/MT using Rapid-Cycling Temperature Swing Adsorption (RC-TSA) are as follows: (1) Increased capacity at different CO 2 concentrations, (2) Increased sorbents cycle life, (3) Increased O 2 resistance, and (4) Decreased steam requirement to extract 1 kg of CO 2 .

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Using the Principals of Electrochemistry to Understand and Overcome the Complicated Degradation Mechanisms of Silicon Anodes in Lithium-Ion Batteries [Slides]

Battery technology is the most significant problem facing widespread market adoption of battery electric vehicles (BEVs). The low energy density of state-of-the-art lithium-ion electrode materials leaves BEV owners and prospective buyers with lower ranges than a comparable internal combustion engine (ICE) vehicle, and vulnerable to a nascent fast charging network. Silicon has the potential to increase the anode energy density by nearly ten times compared to the incumbent material (graphite) and make BEVs a more competitive transportation option. However, lithiated silicon is extremely reactive towards components of the electrolyte and forms a heterogeneous, complicated, and dynamic solid at its surface known as the solid electrolyte interphase (SEI). Ideally, the SEI would passivate the silicon surface, but continuous chemical degradation persists even when the battery is not operating. This reactivity reduces silicon anode lifetimes well below the necessary standards for BEVs. The NREL-led Silicon Consortium Project is dedicated to understanding and solving these mechanisms of degradation. Here, I will discuss an electrochemical method that provides deep insights into the silicon interface during battery operation. I will link these observations to fundamental electrochemical principals and how they translate into actionable strategies that extend the lifetime of silicon anodes.

25 ENERGY STORAGE↗
DOE OSTI2020

A transported Livengood-Wu integral model for knock prediction in CFD simulation

This work describes the development of a transported Livengood-Wu (L-W) integral model for computational fluid dynamics (CFD) simulation to predict auto-ignition and engine knock tendency. The currently employed L-W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition combustion for three different fuels, showing good accuracy in prediction of auto-ignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L-W integral model is coupled with G-equation model to indicate end-gas auto-ignition and knock tendency in CFD simulations of a direct injection spark ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely Co-Optima Alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multi-variable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency indicated by the fuel energy contained in the auto-ignited region of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

Yue, Zongyu↗
DOE OSTI2023

One System, Many Models: Designing a Surrogate Model for Sulfur Thermal Energy Storage: Preprint

Industrial process heating (IPH) relies primarily on thermal energy generated by fossil fuel combustion to produce, treat, and alter manufactured goods. Thermal energy storage (TES) helps reduce the carbon footprint of IPH systems by facilitating the utilization of renewable and waste heat sources. A promising new TES technology uses elemental sulfur as the heat-storage medium. The design of sulfur TES systems can be evaluated with the aid of computational fluid dynamics (CFD). However, the computational cost of such CFD efforts is prohibitive to comprehensive optimizations over design parameters. To reduce this computational cost, machine learning (ML) models can be developed to act as surrogates for CFD. In this paper, we describe the process of building and evaluating surrogate ML models for facilitating optimization of sulfur TES systems for IPH. To enforce the thermodynamic relationship between the two modeled quantities, we develop a hybrid model for sulfur temperature using both direct predictions of temperature and calculations of temperature from predictions of the heat transfer coefficient. The hybrid model enforces this constraint at the expense of the slightly reduced accuracy compared to two disjoint models. The overall high accuracy observed in our model evaluation demonstrates the usefulness of such surrogate modeling for studying TES systems. This work contributes to the field of TES surrogate modeling by offering a novel accurate hybrid approach to predicting simultaneously the heat transfer coefficient and temperature of the heat storage medium.

computational fluid dynamics↗
DOE OSTI2020

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗
DOE OSTI2020

ReaxFF-based molecular dynamics study of bio-derived polycyclic alkanes as potential alternative jet fuels

This work investigates the initial stages of the pyrolysis of HtH-1 (C 18 H 32 ; 2,2,7,7,8a,8b-hexamethyl-dodecahydrobiphenylene) and HtH-2 (C 18 H 34 ; 1,1',3,3,3',3'-hexamethyl-1,1'-bi(cyclohexane)), which are bio-derived polycyclic alkanes and potential jet fuels, using ReaxFF force field based molecular dynamics (MD) simulations. Global Arrhenius parameters, such as activation energies and pre-exponential factors, are calculated and used to analyze the overall decomposition kinetics of the fuels. HtH-1 decomposes faster than HtH-2 at the same temperature and density conditions, and they have a faster decomposition rate compared to some existing jet-fuels, such as JP-10. A systematic reaction analysis framework developed in this work is applied to determine a temperature-dependent decomposition mechanism. At lower temperature, the central C-C bond connecting the two cyclohexane rings is dominantly broken in both HtH-1 and HtH-2. However, C-CH 3 bond breaking becomes dominant with increasing temperature due to the large increase in entropy during this reaction. Major products from HtH-1 are C 5 H 8 and C 4 H 8 , and those from HtH-2 are C 4 H 8 and C 2 H 4 . The major products predict that HtH-1 has a higher sooting tendency than HtH-2, which is consistent with measurements. The impact of HtH-2 on the pyrolysis of HtH-1 is also investigated in their binary mixtures. HtH-1 and HtH-2 decompose by unimolecular reactions, and they rarely interact with each other during the pyrolysis of the mixtures. Furthermore, this work demonstrates that ReaxFF can be used to investigate pyrolysis and combustion chemistry of existing or future fuels and to contribute to the development of their chemical kinetic models without any a priori input and chemical intuition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗