Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity Coefficients”

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

Modernizing the active neutron collar poison rod correction using simulations to enhance PWR fresh fuel verification

The declared inventory of 235 U in fresh, low enriched uranium, nuclear fuel assemblies is routinely verified by nuclear safeguards inspectorates using the UNCL (Uranium Neutron Collar – Light Water Reactor Fuel) nondestructive assay instrument. The trend in modern fuels is towards higher initial enrichment, which in turn requires a larger number of burnable poison pins with higher Gadolinia concentration to hold down the initial reactivity. UNCL assay error for modern fuel assemblies typically exceeds 10%. The traditional algorithm used to correct the response for burnable poison content needs revision to achieve accuracy comparable to the $\sigma_\text{R}$ of 4.5% expected for poison-free assemblies. Here we review the problem and use a large set of 287 Monte Carlo simulations based on an experimentally benchmarked model corresponding to the 16×16 PWR array used in Brazilian Angra type II and III fuel. Simulated relative responses are used to evaluate the functional form of the poison correction and determine parameters optimal for 16×16 PWR assemblies. This update is conducted using between 4–24 Gadolinia pins ranging 2–11 wt% Gd 2 O 3 with assemblies having mean enrichments ranging 2.5%–5%. Experimental validation of the updated coefficients using seven measurements shows bias is reduced by an order of magnitude and $\sigma_\text{R}$R reduced from 10% to less than 2%. Two updated sets of coefficients are provided, the most accurate of which is directly useable in the existing analysis code (INCC) in use by inspectorate at fuel fabrication facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗
DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗
DOE OSTI2021

Influence of NOx Chemistry on the Prediction of Natural Gas End-Gas Autoignition

Research on the influence of NOx (NO/NO2) chemistry on HC oxidation has shown that NOx plays an important role in the promotion/inhibition of HC autoignition [1-6]. Current data for n-heptane, iso-octane, ethanol, and toluene [3, 4] have shown that NO chemistry inhibits autoignition at low temperatures and concentrations. However, at medium and high temperatures and low concentrations, autoignition is promoted, decreasing the Negative Temperature Coefficient (NTC) effect. This autoignition promotion, however, becomes weaker as NO concentrations are increased. Additionally, fundamental kinetics studies have suggested that NOx chemistry also promotes the autoignition of smaller HC, such as CH4 [7-10], C2H6 [11-14], and C3H8 [15]. Additionally to the chemical kinetics studies, a growing body of evidence supports the importance of NOx chemistry on SI engine simulations [16, 17]. Foong et al. [18] suggested that the addition of NO in their models significantly improved the agreement between simulations and experiments where the onset of knock was advanced due to the presence of NO, and Morganti et al. [19] showed that NO addition in the residual-gas was necessary for their models to obtain good agreement with experimental data. Lastly, Mohr et al. [20, 21] investigated the effect of Natural Gas (NG) reactivity and EGR substitution rate and composition on homogeneous ignition delay, flame speeds, and EGAI for stoichiometric NG/oxidizer/EGR blends. They observed that the addition of EGR composed of only inert species (CO2 and Ar) suppressed homogeneous autoignition in all NG fuels tested under all conditions. However, for all NG fuels tested, increasing EGR rates with reactive species (Ar, CO2, CO, and NO) promoted homogeneous autoignition under all conditions, thus reproducing the same trends for various NG fuel compositions that were observed in other studies and fuels. For this reason, a study is being conducted to analyze the effects of NOx chemistry on the prediction of NG/air/EGR homogeneous autoignition. This study has been divided into two tasks. First, the addition of NOx chemistry to the ARIES82 mechanism [20] is being validated using homogeneous ignition delay data collected by Mohr [20, 21]. Second, once the new modified mechanism with NOx reactions has been validated, it will be employed on engine multi-dimensional modeling presented in [22, 23] to analyze the influence of NOx chemistry on NG EGAI. Results thus far have shown that homogeneous ignition delay calculations are sensitive to NOx chemistry, where the modified mechanisms captured well all trends and closely matched the homogeneous ignition delays observed by Mohr et al. Also, the results in this work suggest that NOx chemistry is necessary to correctly capture the trends in NG homogeneous autoignition when EGR with reactive species is used since the opposite behavior is observed when only inert species are considered for the EGR composition.

03 NATURAL GAS↗
DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗
NASA NTRS2010

Experimental Studies of Liquefaction and Densification of Liquid Oxygen

The propellant combination that offers optimum performance is very reactive with a low average molecular weight of the resulting combustion products. Propellant combinations such as oxygen and hydrogen meet the above criteria, however, the propellants in gaseous form require large propellant tanks due to the low density of gas. Thus, rocketry employs cryogenic refrigeration to provide a more dense propellant stored as a liquid. In addition to propellant liquefaction, cryogenic refrigeration can also conserve propellant and provide propellant subcooling and propellant densification. Previous studies analyzed vapor conditioning of a cryogenic propellant, with the vapor conditioning by either a heat exchanger position in the vapor or by using the vapor in a refrigeration cycle as the working fluid. This study analyzes the effects of refrigeration heat exchanger located in the liquid of the common propellant oxidizer, liquid oxygen. This study predicted and determined the mass condensation rate and heat transfer coefficient for liquid oxygen.

Partridge, Jonathan Koert↗
DOE OSTI2024

The photoionization of methoxymethanol: Fingerprinting a reactive C 2 oxygenate in a complex reactive mixture

Methoxymethanol (CH 3 OCH 2 OH) is a reactive C 2 ether-alcohol that is formed by coupling events in both heterogeneous and homogeneous systems. It is found in complex reactive environments—for example those associated with catalytic reactors, combustion systems, and liquid-phase mixtures of oxygenates. Using tunable synchrotron-generated vacuum-ultraviolet photons between 10.0 and 11.5 eV, we report on the photoionization spectroscopy of methoxymethanol. We determine that the lowest-energy photoionization process is the dissociative ionization of methoxymethanol via H-atom loss to produce [C 2 H 5 O 2 ] + , a fragment cation with a mass-to-charge ratio (m/z) = 61.029. Here, we measure the appearance energy of this fragment ion to be 10.24 ± 0.05 eV. The parent cation is not detected in the energy range examined. To elucidate the origin of the m/z = 61.029 (C 2 H 5 O 2 ) fragment, we used automated electronic structure calculations to identify key stationary points on the cation potential energy surface and compute conformer-specific microcanonical rate coefficients for the important unimolecular processes. The calculated H-atom dissociation pathway results in a [C 2 H 5 O 2 ] + fragment appearance at 10.21 eV, in excellent agreement with experimental results.

Hansen, Niko A.↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Optimum power analysis of a self-reactive wave energy point absorber with mechanically-driven power take-offs

This paper studies a self-reactive ocean wave energy converter (WEC) where energy is absorbed through the relative motion between a floating buoy on the ocean surface and a submerged body. Two types of direct-drive power take off (PTO) systems are examined. One adopts a ball screw system in inverse driving mode with mechanical motion rectifier (MMR) that converts the oscillating bi-directional input heave motion into unidirectional rotation to drive an electromagnetic generator, and another one uses a ball screw system to directly drive the generator (called non-MMR). Dynamic models for both types of PTOs are established and integrated with the overall WEC system model in both the time and frequency domains. Numerical simulation is used to investigate the dynamic performance of both PTOs. The influences of the PTO inerter, which is from the rotational inertia of the flywheel and generator, and the WEC drag damping coefficient are discussed and explored. The analytical closed-form solution for the optimum condition of the non-MMR system is derived based on the criterion of maximum power absorption, and numerical simulation yields the optimum condition for the MMR system. Overall, the results under regular and irregular waves show that the PTO inerter in the MMR system can improve the power absorption for small wave periods and maintain the same performance for large wave periods. The PTO inerter in the non-MMR system shifts the peak frequency. The WEC drag damping negatively influences the power absorption in both PTO systems, especially at large wave periods.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

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

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
NASA NTRS2023

Kankakee Water Resources: Monitoring Temperature and Vegetation to Detect River Flow Impediments at Energy Intake Structures

In recent years, unpredictable grassing events have occurred at the Dresden Generating Station, located on the Kankakee River in northern Illinois. Grassing events are characterized by large mats of aquatic vegetation that accumulate downstream, resulting in the clogging of water intake structures and leading to major disruptions in power generation. Currently, employees at the Dresden Generating Station are responsible for reactively responding to each grassing event individually. This project, in partnership with Constellation Nuclear and the United States Geological Survey (USGS), assessed the feasibility of using Earth observations (Landsat 9 OLI-2, Landsat 8 OLI, Sentinel-2 MSI, DOVE PlanetScope, WorldView-3, and GPM IMERG) to detect floating aquatic vegetation within the Kankakee River and identify predictive factors that trigger grassing events, as doing so will provide the Dresden Generating Station the ability to anticipate future grassing events and enhance general hydrologic modeling efforts held by the USGS. The results of this study illustrated that, while aquatic vegetation can be detected by satellites with up to moderate spatial resolution (30 m), temporal resolution is a major limiting factor for tracking movements in floating aquatic vegetation and identifying predictive measures for these events. In addition, correlation results suggest a possible negative relationship between grassing events and river discharge (-0.875 correlation coefficient). In the future, pairing these results with ground control surveys and sensors with higher temporal capabilities would allow our project partners to predict and proactively address future grassing events, ensuring the reliable operation of the Dresden Generating Station.

Marisa Smedsrud↗
DOE OSTI2020

Radiolytic degradation of formic acid and formate in aqueous solution: modeling the final stages of organic mineralization under advanced oxidation process conditions

The successful use of advanced oxidation processes to treat aqueous solutions containing undesirable organic species requires the degradation of these species to lower molecular weight, lower hazard compounds. Safe application of this technology requires a thorough understanding of the mechanisms of degradation. These oxidative transformations are mainly initiated by the reactions of reactive oxygen species, particularly hydroxyl radicals. These react with organic molecules to generate carbon-centered radicals. In the presence of dissolved oxygen, the carbon-centered radicals are next converted to peroxyl radicals, which then decay to lower molecular weight species by multiple mechanistic pathways. Formic acid and its conjugate base formate are the last stable chemical species produced immediately before the complete mineralization of any organic molecule undergoing oxidative degradation in aqueous solution. Once understood, the radical-induced chemistry of formic acid/formate under these conditions has wide applicability in all advanced oxidation technologies. To develop this quantitative knowledge, here we have performed a series of 60 Co gamma irradiation studies on aqueous formic acid/formate over different pH and solution conditions. The measured species concentration changes, as a function of applied dose, are compared with the predictions of a kinetic computer model constructed from literature reactions and reported rate coefficients. The excellent agreement found between the results and modeling gives confidence in the mechanism presented here and provide the first complete computer model for the radiolytic degradation of formic acid in water.

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

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗
DOE OSTI2022

Understanding the low-temperature chemistry of 1,2,4-trimethylbenzene

1,2,4-trimethylbenzene is an important representative aromatic component of gasoline/diesel/jet fuels and thus it is necessary to understand its low-temperature chemistry. In this paper, ignition delay times (IDTs) of both 1,2,4-trimethylbenzene (124TMB) and its blends with n -heptane were measured at engine-like conditions using both a high-pressure shock tube and a rapid compression machine for fuel in ‘air’ mixtures at pressures of 10 and 30 atm and at temperatures in the range 600 – 1100K. The experiments in this study show for the first time that 124TMB presents a two-stage ignition behavior at engine relevant conditions. Blending n -heptane with 124TMB can significantly increase mixture reactivity at temperatures below 1000K. A new detailed mechanism has been developed to simulate the experimentally measured IDT data. The mechanism can capture well the two-stage ignition behavior as well as the ignition delays at different pressures, equivalence ratios over a wide temperature range, for both pure fuels and their blended mixtures. Flux analyses show that the benzylic radicals (formed via H-atom abstraction from the methyl groups ortho-sites on 124TMB) can add to O 2 forming $\mathrm{R\dot{O}_2}$ radicals, which can isomerize to $\dot{Q}$OOH by intramolecular H-atom transfer from the ortho- methyl group and these $\dot{Q}$OOH radicals undergo a second addition to O 2 . This is analogous to the chain branching reaction pathways of alkanes. The chain branching reaction pathways are responsible for the first-stage heat release of 124TMB. In conclusion, the competitions between chain branching and both chain propagating and chain termination reaction pathways lead to a less pronounced negative temperature coefficient (NTC) behavior for 124TMB oxidation, compared to two-stage ignition behavior observed for alkanes and other fuels.

1,2,4-trimethylbenzene↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
NASA NTRS2017

The Iodine Satellite (iSAT) Propellant Feed System - Design and Development

The development, modeling, and testing of components and subsystems required to feed iodine propellant to a 200-W Hall thruster and cathode are described. This work aims to address design deficiencies and issues associated with the propellant feed system that were revealed by an integrated thruster-cathode-feed system test. The feed system design is modified to use materials that are more resistant to the highly-reactive nature of iodine propellant. Dynamic modeling indicates that the inclusion of additional constraints on feed system tubing will reduce the vibrationally-induced stresses that occur during launch. Full spacecraft thermal modeling show that the feed system heater power levels are sufficient to heat the tank and propellant lines to operating temperatures, where iodine in the tank is sublimed to supply propellant for operation and the tubing is elevated in temperature to keep propellant from redepositing to block the flow. Experiments are conducted to demonstrate that is it possible through the application of heating to clear an iodine deposit blocking the flow. Deposits in the low-pressure portion of the system near the exit to vacuum are shown to be relatively easy to remove in this manner while blockages forming upstream nearer to the higher-pressure propellant tank require significantly more effort to remove. Fluid flow modeling of the feed system is performed, exhibiting some qualitative agreement with experimental data. However, the highly viscous nature of the fluid flow and the dependence of the component flow coefficients on the Reynolds number are likely causes of the generally-poor quantitative agreement between the modeling results and experimentally-measured fluid flow properties.

Polzin, Kurt A.↗
NASA NTRSDate not supplied

Assessing the Feasibility of Using Various Earth Observations to Monitor Environmental Trends Associated with River Flow Impediments Near Energy Intake Structures

Grassing events, characterized by the release of river vegetation in large quantities, are often responsible for major flow impediments surrounding industrial water intake structures. These sudden flow impediments halt critical municipal functions, such as cooling at energy generating facilities, leading to widespread energy disruptions for surrounding communities. Because the origin and cause of grassing events in rivers are largely unknown, employees at these facilities can only reactively respond after they occur by pausing energy generation to manually remove the accumulation of aquatic vegetation. In 2023, the Dresden Generating Station, in collaboration with the United States Geological Survey, began looking into using remote sensing methodologies to locate floating aquatic vegetation surrounding their intake structure located along the Kankakee River in northern Illinois. This project contributed to this effort by conducting a case study that compared the performance, practicality, and feasibility of using various Earth observations (Landsat 9 OLI-2, Landsat 8 OLI, Sentinel-2 MSI, and DOVE PlanetScope) and vegetation indices (NDVI, EVI, SAVI, and GCI) to monitor aquatic vegetation in the Kankakee River. Additionally, this study incorporated several environmental metrics to identify potential triggers for grassing events. Results from this project found that while Landsat sensors provide a lower spatial resolution than commercial satellite Earth observations, areas of aquatic vegetation were similarly identified; therefore, it may not be necessary to acquire expensive, high-resolution datasets for continued monitoring. Finally, preliminary correlation analyses showed a potential negative relationship between the river discharge and the presence of aquatic vegetation (-0.875 correlation coefficient), suggesting that periods of low flow could lead to large releases in aquatic vegetation.

Marisa Smedsrud↗