Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Agreements”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Instance Segmentation for Direct Measurements of Satellites in Metal Powders and Automated Microstructural Characterization from Image Data

In this work, we propose instance segmentation as a useful tool for image analysis in materials science. Instance segmentation is an advanced technique in computer vision which generates individual segmentation masks for every object of interest that is recognized in an image. Using an out-of-the-box implementation of Mask R-CNN, instance segmentation is applied to images of metal powder particles produced through gas atomization. Leveraging transfer learning allows for the analysis to be conducted with a very small training set of labeled images. As well as providing another method for measuring the particle size distribution, we demonstrate the first direct measurements of the satellite content in powder samples. After analyzing the results for the labeled data dataset, the trained model was used to generate measurements for a much larger set of unlabeled images. The resulting particle size measurements showed reasonable agreement with laser scattering measurements. The satellite measurements were self-consistent and showed good agreement with the expected trends for different samples. Finally, we present a small case study showing how instance segmentation can be used to measure spheroidite content in the UltraHigh Carbon Steel DataBase, demonstrating the flexibility of the technique.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monte Carlo simulations of Cu/Ni–Si–Mn co-precipitation in duplex stainless steels

First-Passage Kinetic Monte Carlo (FPKMC) simulations of species migration in duplex stainless steels were performed in order to establish relationships between alloy content, segregation behavior, and the formation of Ni-Si-Mn rich particles in cast duplex stainless steels during thermal aging. The Ni-Si-Mn-rich second phase forms after extended aging at reactor operating temperatures and degrades alloy properties. Simulations of Ni-Si-Mn cluster formation were validated through comparison with experimental results obtained through Atom Probe Tomography (APT) on similar alloy compositions, identifying several trends. First, Cu promotes the formation of Ni-Si-Mn clusters, but only when Ni or Mn prefer segregation to the surface of Cu particles; without the segregation of these species, the critical composition for the clusters to form was not achieved. Second, the width of precipitate-denuded zones near ?/d interfaces increases with decreasing Cu content. This finding was in strong agreement with APT results, further validating our model. Finally, our model predicts that Ni-Si-Mn cluster formation will be the most extensive when the Si:Mn ratio is approximately 1:1 and the least extensive when one of these key elements is less concentrated (Mn = 0.5 at.% or Si = 0.5 at.%). The implications of these results on how to improve alloy properties are discussed.

Cu co-precipitation, Atom probe tomography (APT), ↗
DOE OSTI2022

Development of grain-scale slip activity and lattice rotation fields in Inconel 718

Using a combination of in-situ high-resolution digital image correlation (HR-DIC), Heaviside-DIC method (H-DIC), and crystal plasticity finite element (CPFE), we investigate the evolution of intragranular lattice rotations and slip activity during monotonic and cyclic loading in a high performance, polycrystalline face centered cubic material. The CPFE employs a quasi-3D model microstructure, which is a highly resolved mirror representation of the experimental in-situ test sample. In agreement, the measurements and calculations reveal that most grains, regardless of their size and lattice orientation, develop intragranular lattice rotation gradients that span the grain. For a small cluster of grains on the deformed material, we perform HR-DIC analysis of slip lines to demonstrate agreement in the active slip systems and changes in this local slip activity across the individual grains. The combined analysis reveals that deforming grains are divided into sub-granular regions of uniform lattice rotation and these regions are most often associated with only one or two active slip systems. The gradient lines that divide them correspond to changes in the predominant slip system. The model is used to examine the evolution of intragranular lattice rotation in a single fully reversed tension-compression cycle. The calculations indicate that intragranular gradients intensify during the reverse loading path as nearest neighboring regions appear to shed lattice rotation, increasing the lattice rotation in some regions, while shutting down rotation in neighboring regions. So these findings provide insight into the irreversible changes that develop within deforming grains at the scale of the grain, particularly the heterogeneous development of intragranular lattice rotation in early stages of deformation, which could serve as precursors to localization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗
DOE OSTI2023

Emergency department overcrowding and its associated factors at HARME medical emergency center in Eastern Ethiopia

Introduction: Emergency department (ED) overcrowding has become a significant concern as it can lead to compromised patient care in emergency settings. Various tools have been used to evaluate overcrowding in ED. However, there is a lack of data regarding this issue in resource-limited countries, including Ethiopia. This study aimed to validate NEDOCS, assess level of ED overcrowding and identify associated factors at HARME Medical Emergency Center, located in Hiwot Fana Comprehensive Specialized Hospital, Harar, Ethiopia. Methods: A cross-sectional study was conducted at the HARME Medical Emergency Center, Hiwot Fana Comprehensive Specialized Hospital, involving a total of 899 patients during 120 sampling intervals. The area under the receiver operating characteristic curves (AUC) was calculated to evaluate the agreement between objective and subjective assessments of ED overcrowding. A multivariable logistic regression analysis was employed to identify factors associated with ED overcrowding and statistically significant association was declared using 95% confidence level and a p-value < 0.05. Results: The interrater agreement showed a strong correlation with a Cohen's kappa (κ) of 0.80. The National Emergency Department Overcrowding Study Score demonstrated a strong association with subjective assessments from residents and case team nurses, with an AUC of 0.81 and 0.79, respectively. According to residents' perceptions, ED were considered overcrowded 65.8% of the time. Factors significantly associated with ED overcrowding included waiting time for triage (AOR: 2.24; 95% CI: 1.54–3.27), working time (AOR: 2.23; 95% CI: 1.52–3.26), length of stay (AOR: 2.40; 95% CI: 1.27–4.54), saturation level (AOR: 2.35; 95% CI: 1.31–4.20), chronic illness (AOR: 2.19; 95% CI: 1.37–3.53), and abnormal pulse rate (AOR: 1.52; 95% CI: 1.06–2.16). Conclusion: The study revealed that ED were overcrowded approximately two-thirds of the time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Serpent transient analyses of GIACINT geometrical change experiments

This paper simulates seven transient experiments with geometrical movements performed at the GIACINT facility of Belarus. These experiments include three slow control rods insertion in more than 25 s, one fast control rod insertion in less than 1 s, and three water moderator draining transients. In the fast control rod insertion experiment, the control rods are first pushed by a mechanical spring and then let fall by gravity in the assembly. In this experiment, the control rods speed varies because of the initial spring force, the gravitational acceleration, and the water buoyancy force. In all other transient experiments, the geometry change occurs with a constant speed. The Serpent computer program was utilized to simulate these transient experiments with geometrical changes by performing two separate computations. One time-independent computation is performed in criticality mode and the other time-dependent computation is performed in dynamic-source mode. The first computation writes the neutron population and delayed neutron precursors for starting the transient simulation. The second computation divides the transient time into several time bins, reads the two files from the first computation, and updates them at the end of each time bin. The Serpent simulation results obtained for this set of experiments are in good agreement with the time dependent experimental measurements. In addition, Serpent and MCNP simulations results show an excellent agreement for the water moderator transient

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

An economic assessment of behind-the-meter photovoltaics paired with batteries on the Hawaiian Islands.

Due to natural variability and uncertainty, the ever-increasing penetration of solar generation in Hawaii presents challenges to power grid operators to maintain reliable system operation. Demand response (DR) has the potential to be a cost-effective tool for Hawaii to reach its aggressive renewable energy goals while maintaining the reliability of power grids. The Hawaii Public Utilities Commission has approved the Hawaiian Electric Company's revised portfolio of DR programs. The companies have released a grid services purchase agreement and subscribed an initial tranche of load into their DR programs. This paper presents innovative analytical methods and comprehensive economic assessment for distributed photovoltaics (PV) paired with battery energy storage systems (BESSs) for two new DR programs, including fast frequency response and capacity grid service. Optimal dispatch and sizing methods are proposed for the paired system considering different tariff schedules and PV compensation programs across five islands. It was found that while the best resource configuration and potential economic benefits vary with tariff structure, a BESS paired with PV can be optimally dispatched to generate multiple value streams simultaneously. Compensation from DR programs is an important value stream to help increase the cost-effectiveness of the integrated system.

Battery energy storage system↗
DOE OSTI2026

Flow pattern and void fraction characterization in nitrogen/water flows through a diamond-type triply periodic minimal surface lattice

Here, this study presents, to the authors knowledge, the first experimental investigation on the void fraction and flow patterns in two-phase flows through a diamond-type Triply Periodic Minimal Surface (TPMS) lattice. An additively manufactured TPMS structure was tested under upward co-current flow of a water/nitrogen mixture. Superficial velocities varied for a total of 42 test conditions (gas: 0.01–2.4 m/s; liquid: 0.01–2.4 m/s; mass flux: 20–2370 kg/m 2 ·s). High-speed video and X-ray imaging enabled time-averaged void fraction measurements and identified six distinct flow regimes which were used to develop a flow pattern map. Comparison of the void fraction data with correlations from literature demonstrated the Rouhani and Axelsson (1970) [49] correlation modified by Steiner (1993) [52] provided the best agreement, which was improved with empirically fit coefficients. This approach predicted the void fraction with errors < ±20% for 64% of the data, with a mean absolute percent deviation (MAPD) of 23%. The measured frictional pressure drop was compared to correlations from literature which captured the observed trends but did not provide good accuracy. The best agreement was found after optimizing the empirical coefficients of the Muller-Steinhagen & Heck (1986) [56] correlation; this approach captured 33% of the data within ±20% with a MAPD of 46% over the full range and captured 71% of the data within ±20% a MAPD of 13.3% at mass fluxes >1100 kg/m 2 -s. These results provide foundational insight into TPMS two-phase flow behavior and inform modeling and design of advanced heat exchange components incorporating TPMS geometries.

42 - ENGINEERING↗
DOE OSTI2023

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Empirical estimation of densities in NaCl-KCl-UCl 3 and NaCl-KCl-YCl 3 molten salts using Redlich-Kister expansion

Densities of molten KCl-NaCl-UCl 3 and KCl-NaCl-YCl 3 ternary systems have been estimated using a multidimensional Redlich-Kister model. Temperature and composition dependent Redlich-Kister functions have been used to generate binary interaction parameters in the outlined ternary salt systems. These binary interactions have been used in the extrapolation to ternary system densities. The results of the density extrapolations by Muggianu interpolation scheme provide agreement within 2–3% for the NaCl-KCl-YCl 3 liquids and 11% in NaCl-KCl-UCl 3 liquids compared to the available experimental data. Modeling NaCl-KCl-UCl 3 molten phase density with a ternary interaction parameter improved the agreement within 4%. Thermophysical modeling used in this study has shown promising results for use in other material properties, such as viscosity, thermal conductivity, and heat capacity of the molten salts. Lastly, the outlined modeling method applied in these specific molten salt ternaries can be used for quaternary or higher multicomponent molten salt systems.

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

Autoignition of CRC diesel surrogates at low temperature combustion conditions: Rapid compression machine experiments and modeling

As federal programs require increasingly stringent engine emissions and fuel economy standards, these ambitions can only be met if next-generation combustion technology is developed focusing on high-efficiency and low-emissions engines. Recent research has indicated the need to operate engines at higher compression ratios and with low temperature combustion (LTC) to achieve the needed gains in engine efficiency and reductions in emissions. Because there is a lack of understanding of the chemistry of diesel fuel components and their mixtures at these LTC conditions, this limits the ability to develop predictive chemical kinetic models that can be used to optimize engine combustion. The current study aims to fill in gaps in fundamental combustion data on surrogate fuel mixtures relevant to diesel fuels. Specifically, four multicomponent diesel surrogates formulated by the Coordinating Research Council (CRC) to emulate an ultra-low-sulfur research-grade #2 certification diesel fuel (CFA), namely V0a (4 components), V0b (5 components), V1 (8 components), and V2 (9 components), have been investigated in a rapid compression machine (RCM) through determination of total and first-stage ignition delay times. Autoignition characteristics of lean to rich fuel/O 2 /N 2 mixtures, for the four CRC surrogates and CFA, have been measured using an RCM at LTC relevant pressures and temperatures, in the ranges of 10–20 bar and 650–1000 K, respectively. The equivalence ratios have been varied by independently changing the oxygen mole fraction and the fuel mole fraction in the test mixtures, thereby illustrating the individual effects of oxygen concentration and fuel loading on diesel autoignition. Autoignition results of these four CRC surrogates are compared among them and with those of CFA. Some degree of agreement in autoignition response between each CRC surrogate and CFA is observed, while discrepancies are also identified and discussed. In addition, a detailed chemical kinetic model for diesel surrogates has been developed and validated against these newly-acquired RCM data. We find that this model shows reasonable agreement with the overall ignition delay time results of the current RCM experiments. Chemical kinetic analyses of the developed model were further conducted to help identify the reactions controlling the autoignition processes and the consumption of fuel components in CRC surrogates.

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

Temporally and spatially resolved X-ray densitometry in a shock tube

Gas densities were measured by x-ray absorption spectroscopy to examine spatial inhomogeneity in density behind incident and reflected shock waves in a miniature (12.7 mm bore) high repetition rate shock tube. Shock waves were generated in argon. The pressure and temperature behind the incident shock were P 2 ~ 1.72 bar and T 2 ~ 800 K, while those behind the reflected shock were P 5 ~ 6.5 bar and T 5 ~ 1480 K. Time-resolved line-of-sight transmission measurements using x-rays at 9 keV photon energy were made along various axial and transverse locations covering the entire shock tube cross section. The x-ray transmission at each location was converted to pathlength-integrated densities, which in some conditions could be further converted to pathlength-averaged densities. The measured gas densities were compared with those calculated by the normal shock wave relations, and good agreement was found in the preshock and immediately behind the incident shock regions. Similar agreement was found for reflected shock conditions very near the endwall. However, increasingly large differences between the measured and calculated gas densities were found as distance from the endwall increased. Furthermore, reconstruction by Abel inversion allowed time and radially resolved densities to be obtained from the transverse measurements. While the profiles measured in this experiment are of insufficient signal-to noise level and resolution to definitively assign boundary layer thicknesses, future improvements may be able to do this. The radially resolved densities reveal growing sidewall thermal boundary layers in gases behind the reflected shock. These reconstructions were compared with a model of the thermal boundary layer and are consistent with values of the thermal boundary layer thickness being between 0.25 mm and 1 mm at 0.7 ms after shock reflection.

42 ENGINEERING↗
DOE OSTI2022

Time-resolved measurements of HO 2 radical in a heated plasma flow reactor

Time-resolved, absolute HO 2 number density in diluted H 2 –O 2 –Ar, CH 4 –O 2 -Ar, and C 2 H 4 –O 2 –Ar mixtures excited by a repetitive ns pulse discharge in a heated plasma flow reactor is measured by Cavity Ringdown Spectroscopy (CRDS). The experimental results are obtained at $\textit{T}$ = 300-600 K and $\textit{P}$ = 130 Torr, both during the discharge pulse burst and in the afterglow. In this work, the HO 2 number density is inferred from the CRDS data using a spectral model exhibiting good agreement with previous measurements of absolute HO 2 absorption cross sections. In the room-temperature H 2 –O 2 mixture, as well as in CH 4 –O 2 and C 2 H 4 –O 2 mixtures over the entire temperature range studied, HO2 is generated only during the discharge burst and decays in the afterglow. However, in the H 2 –O 2 mixture at elevated temperatures, $\textit{T}$ = 400-600 K, HO 2 persists in the afterglow up to 10 ms after the discharge burst, comparable with the flow residence time in the reactor. Comparison with kinetic modeling shows that the sustained reactivity after the source of radicals is turned off is due to a chain propagation / hydrogen oxidation process, which dominates the radical recombination reactions. The kinetic modeling predictions are in good agreement with the relative HO 2 number density measured in all three mixtures, although the model underpredicts the absolute number densities in H 2 –O 2 at $\textit{T}$ = 400-600 K by up to a factor of two. Detection of the sustained low-temperature reactivity in H 2 –O 2 , initiated by the radical generation in the plasma, suggests that the plasma excitation may also affect kinetics of oxidation and reforming of fuels exhibiting low-temperature chemistry below hot ignition point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Non-equilibrium low-temperature plasma-assisted combustion of iso-octane: Perturbing pyrolysis and oxidation kinetics

Here, in this study, a plasma-coupled flow reactor facility is used to examine the effects of non-equilibrium low-temperature plasmas on perturbing the pyrolysis and oxidation kinetics of iso-octane. Experiments were performed in highly dilute reactive mixtures of nitrogen, at near isothermal conditions for temperatures ranging from 523 K to 1203 K. Experiments cumulatively demonstrated enhanced chemical reactivity with the plasma for temperatures below 900 K, and a lowering of the hot-ignition temperature. Detailed kinetic insight was derived from a 0D plasma-coupled kinetic model, utilizing a constructed mechanism that combined both plasma-specific chemistry and the neutral combustion chemistry. For pyrolysis conditions, the model displayed relatively good agreement with fuel consumption and the formation of most intermediates compared to the experimental data, demonstrating the model is able to accurately predict primary radical formation from the plasma directly interacting with the fuel. Enhanced reactivity was attributed to collisional quenching of excited-states of N 2 with fuel, which led to efficient fuel fragmentation and enhancement of the H-radical flux. For oxidation conditions, the model displayed satisfactory agreement with the experiments. Model predictions were able to accurately predict fuel consumption and most intermediate speciation data for T > 800 K, but most discrepancies were towards T < 800 K in particular with oxygenated intermediates. In the presence of oxygen, plasma effects were predominantly spent on efficient enhancement of O- and H-radical fluxes, leading to further fuel fragmentation and initiation of both the OH- and HO 2 -radical pools. Subsequent reactivity of iso-octane was then dictated by the response of the temperature-dependent neutral chemistry. At low-temperatures (T = 643 K), enhanced fuel radicals and O 2 -additon chemistry lead to the formation of oxygenated species, while at intermediate temperatures (T = 843 K) net decrease in OH-radical reactivity led to an increase in hydrocarbon speciation. Near the self-ignition threshold (T = 1163 K), radicals generated by high-temperature branching reactions dominate the oxidation process and effectively ignition. This study ultimately demonstrated that the enhancement of radicals afforded by the plasma causes a deviation in known understanding of iso-octane kinetics in some regards and warrants future studies to reconcile these discrepancies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗