Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

CHARACTERISTICS OF CYCLE-TO-CYCLE COMBUSTION VARIABILITY AT PARTIAL-BURN LIMITED AND MISFIRE LIMITED SPARK TIMING UNDER HIGHLY DILUTED CONDITIONS

Diluted combustion with exhaust gas recirculation (EGR) has been widely employed to improve the fuel economy of spark ignition engines. The combustion kinetics, however, are affected and the flame propagation speed is decreased. In order to compensate for this adverse effect, the spark timing needs to be recalibrated to achieve maximum brake torque (MBT). At high levels of EGR dilution, the spark timing is constrained by two ignition limits: 1) the partial-bum limit where the spark timing is retarded from MBT and 2) the misfire limit where the spark timing is too advanced. This work uses a probabilistic framework to capture the differences between both ignition limits. In particular, it introduces the concept of a nominal indicated mean effective pressure (IMEP) distribution based on the stochastic properties of the cycle-to-cycle variability (CCV) at nominal stable conditions. By defining a nominal band where fully burned cycles occur with high probability, we introduce a cycle classification method that can be used to 1) determine the level of randomness of misfire and partial-burn events, and 2) measure CCV. The new CCV metric based on the density of the nominal band is compared with the traditional coefficient of variation of IMEP (CoVIMEP). It is shown that the nominal band concept, together with the CoVIMEP, can help to discern between partial-burn limited and misfire limited conditions. Furthermore, the Kullback-Leibler divergence is used to demonstrate that the IMEP distribution is significantly different between nominal and partial-burn/misfire limited conditions. Experiments are carried at various EGR levels and spark timings while recoding in-cylinder pressure at steady state. Although the emphasis of this work is to characterize the differences of both ignition limits from a probabilistic point of view, similarities between partial-burn cycles at either limiting conditions are also discussed.

Maldonado, Bryan↗
DOE OSTI2020

Requirements and Conceptual Design of Off-gas Systems for the Reprocessing of Metallic Fuels

An assessment has been conducted to determine how key regulations regarding volatile radionuclide emissions to the atmosphere may apply to the off-gas streams associated with electrochemical reprocessing. The scope of this assessment was based upon a generic electrochemical reprocessing scheme with a throughput rate of 200 MTIHM/y applied to metallic fuel discharged from a sodium fast reactor (SFR), but the findings are able to be translated to other advanced nuclear scenarios as merited. Air dispersion modeling was performed using the EPA CAP-88 model and evaluated the uncontrolled decontamination factors (DFs) that would be required to achieve regulatory compliance with the dose-based limits set forth by EPA regulation 40 CFR 190.10(a). These DFs were compared to those required by fuel cycle–based limits set forth by EPA regulation 40 CFR 190.10(b). Two theoretical sites with disparate climatological conditions were selected for air dispersion modeling (Idaho and Tennessee). The radionuclides modeled included 3 H, 85 Kr, 129 I, and selected alpha-emitting transuranic isotopes (referred to here as 239 Pu-TRU <1y ). It was found that the fuel cycle-based limits in 40 CFR 190.10(b) are most restrictive for 85 Kr and 239 Pu-TRU <1y , with DFs of 3 and 6.1E+09, respectively. The dose-based limit as derived from 40 CFR 190.10(a) could require mitigation of tritium in some scenarios, with an estimated DF of about 3 for the reference scenarios. The fuel cycle-based limit for 129 I resulted in a DF of about 240 for the reference scenario. The need for iodine mitigation based on dose to the public depended upon the physical form of iodine as either particulate or vapor-phase species. Emission of iodine from the facility as a vapor necessitated DFs of about 2 but emission as a particulate would require DFs >6,000 to meet thyroid dose-based limits. Effects of physical form on needed iodine mitigation are significant, but the understanding of speciation of iodine both during electrochemical reprocessing and after release to the atmosphere is limited. The electrochemical processing unit operations were evaluated to identify potential release points for the volatile radionuclides and to assess the potential for retention of the radionuclides within the process (thus decreasing the need for mitigation). Mitigation strategies for 3 H, 85 Kr, 129 I, and 239 Pu-TRU <1y were identified. In all cases, there are reasonably achievable pathways to regulatory compliance, although in some cases additional R&D is merited to verify the chemical speciation of these isotopes and to develop and demonstrate potential treatment technologies for this application. Whether or not additional off-gas controls (beyond common operations such as HEPA filtration and oxygen and moisture control) are needed for any of these regulated or volatile radionuclides depends on the (a) type of facility (NRC-regulated or DOE), (b) used fuel process rate, (c) used fuel burnup and composition, (d) speciation and retention of volatile radionuclides in the process and in the cell gas cleanup system, (e) site-specific parameters such as location, meteorology, stack height, and site boundaries, and (f) levels of conservatism and safety factors used in assessing compliance to air emissions regulations. Performance of this assessment revealed several areas where information is lacking or additional research is required in order to better determine if or what kinds of off-gas control might be needed. First, and most significantly, the understanding of the chemical speciation and physical form and partitioning of iodine during electrochemical processing operations is lacking and prevents the ability to accurately assess the potential iodine mitigation requirements. Future research in this area should be multifaceted and include thermodynamic modeling of iodine speciation in different process steps, experiments to quantify the kinetics of vapor-phase and melt-phase transitions, bench-scale experiments to determine the potential chemical and physical form of iodine emissions from the electrorefining process, and verification of iodine behavior with experiments utilizing operational facilities. Similarly, an improved understanding of iodine behavior in the environment after release from the facility stack will be required to refine dose estimations, as particulate and vapor-phase emissions result in significantly different doses to the MEI.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗
NASA NTRS2020

The Flux of Impact Ejecta on the Lunar Surface from Scaling Considerations: Implications for Operational Hazards and Geomorphic Forcing

The impact cratering process has been critical to the evolution of the Moon’s surface over its geologic history and remains an important ongoing process today. Impact events have a major local effect, but also excavate ejecta particles that re-impact the lunar surface over a wide area. Quantifying the flux of ejecta to a given point on the Moon is the subject of this work. We also estimate how this flux is partitioned into different particle sizes and different ejecta velocities. Motivation: There are two main factors motivating this work. First, and most critically, is the assessment of the hazard posed by impact ejecta for future surface exploration (i.e., to infrastructure, spacesuits, etc.). LROC observations of new craters have led to the reemphasized need to consider this hazard. In fact, a hazard assessment of this type was made prior to Apollo, although some of the underlying assumptions of that work are now clearly obsolete (see [4]). We also now know much more about the impactor flux, scaling of impact events, and scaling of ejecta than was known in the 1960's, so revisiting this hazard assessment is appropriate.We note that also have recently revisited the earlier hazard estimates and independently revised them downward using an entirely different analytical approach. The second motivation is that several recent papers have argued that the flux of distal ejecta is the controlling factor in how fast the lunar surface evolves. For this reason, improving understanding of the ejecta mass flux and how the flux translates into geomorphic work is of interest. To be clear, it is obvious that the ejecta mass flux is much larger than the primary impactor mass flux – indeed, this is self-evident because the craters excavated by hypervelocity impacts are much larger than their impactors. On the other hand, the energy delivered by a given primary to the surface is larger than the sum of the energy delivered by all its associated ejecta, as required by conservation, aggravated by the fact that not all of an impactor’s kinetic energy is partitioned into ejecta excavation. If distal ejecta and secondaries control lunar geomorphic evolution, this suggests that re-impacting ejecta must more efficiently translate their energy into geomorphic work than primaries. It is also easy to imagine the relative efficiency of primary and secondary impacts to do geomorphic work varying with the size of the primary. Considering the details of this process is thus of significant interest for lunar geomorphology.

Fassett, C. I.↗
NASA NTRS2003

Radiation Tests of the Extravehicular Mobility Unit Space Suit for the International Space Station Using Energetic Protons

Measurements using silicon detectors to characterize the radiation transmitted through the EMU space suit and a human phantom have been performed using 155 and 250 MeV proton beams at LLUMC. The beams simulate radiation encountered in space, where trapped protons having kinetic energies on the order of 100 MeV are copious. Protons with 100 MeV kinetic energy and above can penetrate many centimeters of water or other light materials, so that astronauts exposed to such energetic particles will receive doses to their internal organs. This dose can be enhanced or reduced by shielding - either from the space suit or the self-shielding of the body - but minimization of the risk depends on details of the incident particle flux (in particular the energy spectrum) and on the dose responses of the various critical organs. Data were taken to characterize the beams and to calibrate the detectors using the beam in a treatment room at LLUPTF, in preparation for an experiment with the same beams incident on detectors placed in a human phantom within the EMU suit. Nuclear interactions of high-energy protons in various materials produce a small flux of highly ionizing, low-energy secondary radiation. Secondaries are of interest for their biological effects, since they cause doses and especially dose-equivalents to increase relative to the values expected simply from ionization energy loss along the Bragg curve. Because many secondaries have very short ranges, they are best measured in passive track detectors such as CR-39. The silicon detector data presented here are intended to supplement the CR-39 data in regions where silicon has greater sensitivity, in particular the portion of the LET spectrum below 5 keV/micron. The results obtained in this study suggest that optimizing the radiation shielding properties of space suits is a formidable task. The naive assumption that adding mass can reduce risk is not supported by the data, which show that reducing the dose delivered at or near the skin by low-energy particles may increase the dose delivered by energetic particles to points deeper in the body.

Zeitlin, C.↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
NASA NTRS2005

Early immune response and regulation of IL-2 receptor subunits

Affymetrix oligonucleotide arrays were used to monitor expression of 8796 genes and probe sets in activated T-cells; analysis revealed that 217 genes were significantly upregulated within 4 h. Induced genes included transcription factors, cytokines and their receptor genes. Analysis by semi-quantitative RT-PCR confirmed the significant induction of IL-2, IL-2R(gamma) and IL-2R(alpha). Forty-eight of the 217 induced genes are known to or predicted to be regulated by a CRE promoter/enhancer. We found that T-cell activation caused a significant increase in CREB phosphorylation furthermore, inhibition of the PKC pathway by GF109203 reduced CREB activation by 50% and inhibition of the PKA pathway caused a total block of CREB phosphorylation and significantly reduced IFN(gamma), IL-2 and IL-2R(alpha) gene expression by approximately 40% (p<0.001). PKC(theta) plays a major role in T-cell activation: inhibition of PKC significantly reduced the expression of IFN(gamma), IL-2 and IL-2R(alpha). Since PKC blocked activation of CREB, we studied potential cross-talk between the PKC and the PKA/MAPK pathways, PMA-stimulated Jurkat cells were studied with specific signal pathway inhibitors. Extracellular signal-regulated kinase-2 (ERK2) pathway was found to be significantly activated greater than seven-fold within 30 min; however, there was little activation of ERK-1 and no activation of JNK or p38 MAPK. Inhibition of the PKA pathway, but not the PKC pathway, resulted in inhibition of ERK1/2 activation at all time points, inhibition of MEK1 and 2 significantly blocked expression of IL-2 and IL-2R(alpha). Gene expression of IL-2R(alpha) and IFN(gamma) was dependent on PKA in S49 wt cells but not in kin- mutants. Using gel shift analysis, we found that forskolin activation of T-cells resulted in activation of AP1 sites; this increase in nuclear extract AP1 was significantly blocked by MEK1 inhibitor U0126. Taken together, these results suggest that the PKA in addition to PKC and MAPK pathways plays a role in early T-cell activation and induction of IL-2, IL-2R(alpha) and IFN(gamma) gene expression.

Non-NASA Center↗
DOE OSTI2026

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

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

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2024

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

The molecular responses of skeletal muscle satellite cells to continuous expression of IGF-1: implications for the rescue of induced muscular atrophy in aged rats

Approximately 50% of humans older than 85 years have physical frailty due to weak skeletal muscles. This indicates a need for determining mechanisms to combat this problem. A critical cellular factor for postnatal muscle growth is a population of myogenic precursor cells called satellite cells. Given the complex process of sarcopenia, it has been postulated that, at some point in this process, a limited satellite cell proliferation potential could become rate-limiting to the regrowth of old muscles. It is conceivable that if satellite cell proliferative capacity can be maintained or enhanced with advanced age, sarcopenia could potentially be delayed or prevented. Therefore, the purposes of this paper are to describe whether IGF-I can prevent muscular atrophy induced by repeated cycles of hindlimb immobilization, increase the in vitro proliferation in satellite cells from these muscles and, if so, the molecular mechanisms by which IGF-I mediates this increased proliferation. Our results provide evidence that IGF-I can enhance aged muscle regrowth possibly through increased satellite cell proliferation. The results also suggest that IGF-I enhances satellite cell proliferation by decreasing the cell cycle inhibitor, p27Kip1, through the PI3'-K/Akt pathway. These data provide molecular evidence for IGF-I's rescue effect upon aging-associated skeletal muscle atrophy.

NASA Discipline Musculoskeletal↗
DOE OSTI2021

Spot the difference: hydrogen adsorption and dissociation on unsupported platinum and platinum-coated transition metal carbides

Hydrogenation reactions are involved in several processes in heterogeneous catalysis. Platinum is the best-known catalyst; however, there are limitations to its practical use. Therefore, it is necessary to explore alternative materials and transition metal carbides (TMCs) have emerged as potential candidates. We explore the possibility of using cheap TMCs as supports for a Pt monolayer, aiming to reduce the amount of the noble metal in the catalyst without a significant loss of its activity towards H 2 dissociation. Hence, analyzing H 2 dissociation from a fundamental point of view is a necessary step towards a further practical catalyst. By means of periodic DFT calculations, we analyze H 2 adsorption and dissociation on Pt/β-Mo 2 C and Pt/α-WC surfaces, as a function of hydrogen surface coverage ( Θ H ), resembling a more realistic model of a catalyst. H 2 dissociation rates were analyzed as a function of the reaction temperature. Overall, the results show that Pt/C-WC and Pt/Mo-Mo 2 C have a Pt-like behavior for H 2 dissociation at Θ H > 1/2 ML. At a particular temperature of 298 K, Pt/C-WC and Pt/Mo-Mo2C have low energy barriers for H 2 * → 2H* (0.13 and 0.11 eV, respectively), close to the value of Pt (0.06 eV). For the highest coverage, i.e. Θ H = 1, Pt/C-WC has a lower activation energy and a higher reaction rate than Pt. Finally, the H 2 dissociation rate is higher in Pt/Mo-Mo 2 C than in Pt when increasing the temperature above 298 K. Our results put Pt/C-WC and Pt/Mo-Mo 2 C under the spotlight as potential catalysts for H 2 dissociation, with a similar performance to Pt, paving the way for further experimental and/or theoretical studies, addressing the capability of Pt/TMC as practical catalysts in hydrogenation reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Capturing the generation and structural transformations of molecular ions

Molecular ions are ubiquitous and play pivotal roles in many reactions, particularly in the context of atmospheric and interstellar chemistry. However, their structures and conformational transitions, particularly in the gas phase, are less explored than those of neutral molecules owing to experimental difficulties. A case in point is the halonium ions, whose highly reactive nature and ring strain make them short-lived intermediates that are readily attacked even by weak nucleophiles and thus challenging to isolate or capture before they undergo further reaction. Here we show that mega-electronvolt ultrafast electron diffraction (MeV-UED), used in conjunction with resonance-enhanced multiphoton ionization, can monitor the formation of 1,3-dibromopropane (DBP) cations and their subsequent structural dynamics forming a halonium ion. We find that the DBP + cation remains for a substantial duration of 3.6 ps in aptly named ‘dark states’ that are structurally indistinguishable from the DBP electronic ground state. The structural data, supported by surface-hopping simulations and ab initio calculations, reveal that the cation subsequently decays to iso-DBP + , an unusual intermediate with a four-membered ring containing a loosely bound bromine atom, and eventually loses the bromine atom and forms a bromonium ion with a three-membered-ring structure. We anticipate that the approach used here can also be applied to examine the structural dynamics of other molecular ions and thereby deepen our understanding of ion chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Orbital Processing of Eutectic Rod-Like Arrays

The eutectic is one of only three solidification classes that exist. The others are isostructural and peritectic-class reactions, respectively. Simplistically, in a binaryeutectic phase diagram, a single liquid phase isothermally decomposes to two solid phases in a cooperative manner. The melting point minimum at the eutectic composition, isothermal solidification temperature, near-isocompositional solidification and refined solidification microstructure lend themselves naturally to such applications as brazing and soldering; industries that eutectic alloys dominate. Interest in direct process control of microstructures has led, more recently, to in-situ eutectic directional solidification with applications in electro-magnetics and electro-optics. In these cases, controlled structural refinement and the high aspect ratio and regularity of the distributed eutectic phases is highly significant to the fabrication and application of these in-situ natural composites. The natural pattern formation and scaling of the dispersed phase on a sub-micron scale has enormous potential application, since fabricating bulk materials on this scale mechanically has proven to be particularly difficult. It is thus of obvious importance to understand the solidification of eutectic materials since they are of great commercial significance. The dominant theory that describes eutectic solidification was derived for diffusion-controlled growth of alloys where both solid eutectic phases solidify metallically, i.e. without faceting at the solidification interface. Both high volume fraction (lamellar) and low volume fraction (rod-like) regular metallic arrays are treated by this theory. Many of the useful solders and brazements, however, and most of the regular in-situ composites are characterized by solidification reactions that are faceted/non-faceted in nature, rather than doubly non-faceted (metallic). Further, diffusion-controlled growth conditions are atypical terrestrially since gravitationally-driven convection is pervasive. As a consequence, it is important to determine whether these faceted/non-faceted composites behave in the same manner as their doubly non-faceted counterparts, particularly in the presence of convection. Prior analytical convective sensitivity testing of this theory predicted insensitivity. Prior experimental testing of this theory offered broad-based agreement between theory and experiment, though most results were for high volume fraction lamellar eutectics that solidified without faceting at the solidification interface. Directional solidification experiments of low volume fraction rod eutectics under damped (microgravity or magnetic field) conditions, however, have demonstrated significant sensitivity, challenging this fundamental theory. More recent theories have been proposed which introduce kinetic undercooling, faceting, fluid shear of the solute redistribution zone and the possibility that the interface composition is not the same as the bulk liquid composition. This program tests the established and proposed analytical theories and addresses the origins of discrepancies between the experimental and analytical results.

Larson, David J., Jr.↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗