Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Peak Refinement”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2022

Rapid water dynamics structures the OH-stretching spectra of solitary water in ionic liquids and dipolar solvents

In a recent study [J. Phys. Chem. B 126, 4584–4598 (2022)], we have used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands, common to all high-polarity solvents, were assigned to water in several solvation states. Here, in the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [J. Chem. Phys. 128, 224511–224512 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared with the inhomogeneous limit. This narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.

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

Quantifying the impacts of marine aerosols over the southeast Atlantic Ocean using a chemical transport model: implications for aerosol–cloud interactions

The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , a New Mineral Species from the Sisim Placer Zone, Eastern Sayans, Russia, and the Torryweiserite–Ferrotorryweiserite Series

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , occurs as small grains (≤20 µm) among droplet-like inclusions (up to 50 μm in diameter) of platinum-group minerals (PGM), in association with oberthürite or Rh-bearing pentlandite, laurite, and a Pt-Pd-Fe alloy (likely isoferroplatinum and Fe-Pd-enriched platinum), hosted by placer grains of Os-Ir alloy (≤0.5 mm) in the River Ko deposit. The latter is a part of the Sisim placer zone, which is likely derived from ultramafic units of the Lysanskiy layered complex, southern Krasnoyarskiy kray, Russia. The mineral is opaque, gray to brownish gray in reflected light, very weakly bireflectant, not pleochroic to weakly pleochroic (grayish to light brown tints), and weakly anisotropic. The calculated density is 5.93 g·cm –3 . Mean results (and ranges) of four WDS analyses are: Ir 18.68 (15.55–21.96), Rh 18.34 (16.32–20.32), Pt 0.64 (0.19–1.14), Ru 0.03 (0.00–0.13), Os 0.07 (0.02–0.17), Fe 14.14 (13.63–14.64), Ni 13.63 (12.58–14.66), Cu 4.97 (3.42–6.41), Co 0.09 (0.07–0.11), S 29.06 (28.48–29.44), and total 99.66 wt.%. They correspond to the following formula calculated for a total of 31 atoms per formula unit: (Rh 3.16 Ir 1.72 Pt 0.06 Ru 0.01 Os 0.01 ) Σ 4.95 (Fe 4.48 Ni 4.11 Cu 1.38 Co 0.03 ) Σ10.00 S 16.05 . The results of synchrotron micro-Laue diffraction studies indicate that ferrotorryweiserite is trigonal; its probable space group is $R\bar{3}m$ (#166) based on its Ni-analog, torryweiserite. The unit-cell parameters refined from 177 reflections are a = 7.069 (2) Å, c = 34.286 (11) Å, V = 1484 (1) Å 3 , and Z = 3. The c:a ratio is 4.8502. The strongest eight peaks in the X-ray diffraction pattern derived from results of micro-Laue diffraction study [d in Å(hkil)(I)] are 2.7950 ($20\bar{2}5$) (100); 5.7143 (0006) (60); 1.7671 ($22\bar{4}0$) (44.4); 3.0486 ($20\bar{2}1$) (39.4); 5.7650 ($10\bar{1}2$) (38.6); 2.5956 ($20\bar{2}7$) (37.8); 3.0058 ($11\bar{2}6$) (36.5); and 1.5029 ($4\bar{2}\bar{2}12$) (35.3). Ferrotorryweiserite and the associated PGM crystallized from microvolumes of residual melt at late stages of crystallization of grains of Os- and Ir-dominant alloys occurred in lode zones of chromitites of the Lysanskiy layered complex. In a particular case, the residual melt is disposed peripherally around a core containing a disequilibrium association of magnesian olivine (Fo 72.9–75.6 ) and albite (Ab 81.6–86.4 ), with the development of skeletal crystals of titaniferous augite: Wo 40.8–43.2 En 26.5–29.3 Fs 20.3–22.6 Aeg 6.9–9.5 (2.82–3.12 wt.% TiO 2 ). Ferrotorryweiserite represents the Fe-dominant analog of torryweiserite. We also report occurrences of ferrotorryweiserite in the Marathon deposit, Coldwell Complex, Ontario, Canada, and infer the existence of the torryweiserite–ferrotorryweiserite solid solution in other deposits and complexes.

58 GEOSCIENCES↗
DOE OSTI2026

Correlation between infrared spectra features and coverage of different adsorption sites for the NO/Pd(111) system

This study investigates the adsorption of nitric oxide (NO) on Pd(111) surfaces at 200 K using infrared reflection absorption spectroscopy (IRRAS). Peak positions and areas are used as proxies for tracking coverage of different species, serving as groundwork for in-situ experiments that need real-time tracking of surface species. Here, we derived a mathematical model correlating NO coverage with dose, consistent with the Molecular Langmuir Model. Distinct correlations are observed between coverage and spectral features (peak position and peak area for the NO vibration): at lower coverages (θ<0.6), the total integrated peak areas are linearly corelated to coverage; at higher coverages (θ>0.6), the peak center for a compressed-hollow site becomes the primary coverage indicator, showing a linear relationship between wavenumber and coverage. Note that this study addresses coverages below saturation. These findings refine our understanding of NO adsorption and establish a foundation for real-time in situ transient kinetics studies of chemical reactions involving NO on Pd(111).

25 ENERGY STORAGE↗
DOE OSTI2022

An Exact Algorithm for the Linear Tape Scheduling Problem

Magnetic tapes are often considered as an outdated storage technology, yet they are still used to store huge amounts of data. Their main interests are a large capacity and a low price per gigabyte, which come at the cost of a much larger file access time than on disks. With tapes, finding the right ordering of multiple file accesses is thus key to performance. Moving the reading head back and forth along a kilometer long tape has a non-negligible cost and unnecessary movements thus have to be avoided. However, the optimization of tape request ordering has rarely been studied in the scheduling literature, much less than I/O scheduling on disks. For instance, minimizing the average service time for several read requests on a linear tape remains an open question. Therefore, in this paper, we aim at improving the quality of service experienced by users of tape storage systems, and not only the peak performance of such systems. To this end, we propose a reasonable polynomial-time exact algorithm while this problem and simpler variants have been conjectured NP-hard. We also refine the proposed model by considering U-turn penalty costs accounting for inherent mechanical accelerations. Then, we propose a low-cost variant of our optimal algorithm by restricting the solution space, yet still yielding an accurate suboptimal solution. Finally, we compare our algorithms to existing solutions from the literature on logs of the mass storage management system of a major datacenter. This allows us to assess the quality of previous solutions and the improvement achieved by our low-cost algorithm. Aiming for reproducibility, we make available the complete implementation of the algorithms used in our evaluation, alongside the dataset of tape requests that is, to the best of our knowledge, the first of its kind to be publicly released.

Honoré, Valentin↗
DOE OSTI2021

Structure of the E. coli agmatinase, SPEB

Agmatine amidinohydrolase, or agmatinase, catalyzes the conversion of agmatine to putrescine and urea. This enzyme is found broadly across kingdoms of life and plays a critical role in polyamine biosynthesis and the regulation of agmatine concentrations. Here we describe the high-resolution X-ray crystal structure of the E . coli agmatinase, SPEB. The data showed a relatively high degree of pseudomerohedral twinning, was ultimately indexed in the P 3 1 space group and led to a final model with eighteen chains, corresponding to three full hexamers in the asymmetric unit. There was a solvent content of 38.5% and refined R/R free values of 0.166/0.216. The protein has the conserved fold characteristic of the agmatine ureohydrolase family and displayed a high degree of structural similarity among individual protomers. Two distinct peaks of electron density were observed in the active site of most of the eighteen chains of SPEB. As the activity of this protein is known to be dependent upon manganese and the fold is similar to other dinuclear metallohydrolases, these peaks were modeled as manganese ions. The orientation of the conserved active site residues, in particular those amino acids that participate in binding the metal ions and a pair of acidic residues (D153 and E274 in SPEB) that play a role in catalysis, are similar to other agmatinase and arginase enzymes and is consistent with a hydrolytic mechanism that proceeds via a metal-activated hydroxide ion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Desmearing small-angle scattering data by central moment expansions of instrument resolution

This study introduces an enhanced numerical technique tailored specifically for refining 1D small-angle scattering (SAS) intensity profiles affected by smearing. Our primary objective is to address the resolution blurring commonly encountered in SAS data, particularly in systems with clearly defined correlation peaks whose spread aligns with the width of the resolution function at corresponding Q positions. Unlike previous approaches that expanded the SAS intensity using central moment expansion, the new method focuses on expanding the resolution function itself, thus eliminating artificial oscillations observed in smeared spectra due to limitations inherent in our earlier algorithm. This method is straightforward to implement, computationally efficient and consistently performs well in numerical benchmarking. To illustrate its effectiveness, we present a case study of a lamellar phase characterized by distinct peaks in its small-angle neutron scattering intensities.

central moment expansions↗
DOE OSTI2023

SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0

Abstract. Water flowing under ice sheets and glaciers can have a strong influence on ice dynamics, particularly through pressure changes, suggesting that a comprehensive ice sheet model should include the effect of basal hydrology. Modeling subglacial hydrology remains a challenge, however, mainly due to the range of spatial and temporal scales involved – from subglacial channels to vast subglacial lakes. Additionally, networks of subglacial drainage channels dynamically evolve over time. To address some of these challenges, we have developed an adaptive mesh refinement (AMR) model based on the Chombo software framework. We extend the model proposed by Sommers et al. (2018) with a small but significant change to accommodate the transition from unresolved to resolved flow features. We handle the strong nonlinearities present in the equations by resorting to an efficient nonlinear full approximation scheme multigrid (FAS-MG) algorithm. We outline the details of the algorithm and present convergence analysis results demonstrating its good performance. Additionally, we present results validating our approach, using test cases from the Subglacial Hydrology Model Intercomparison Project (SHMIP) (de Fleurian et al., 2018). We finish by presenting a more complex, 100 km-by-100 km synthetic test case with peaks and valleys that we use to investigate the effective pressure distribution as the number of AMR levels increases. These preliminary results suggest that a minimum spatial resolution is needed to properly capture channel features, but additional work is required to precisely quantify this and its impact on accurately modeling the coupled ice sheet–hydrology system. The efficiency of our approach, relying on localized refinement, is also demonstrated. Future work will include coupling the SUbglacial Hydrology MOdel (SUHMO) with the BISICLES AMR ice sheet model (Cornford et al., 2013), both built on the same numerical framework.

58 GEOSCIENCES↗
DOE OSTI2024

ASb 3 Mn 9 O 19 (A = K or Rb): New Mn-Based 2D Magnetoplumbites with Geometric and Magnetic Frustration

Magnetoplumbites are one of the most broadly studied families of hexagonal ferrites, typically with high magnetic ordering temperatures, making them excellent candidates for permanent magnets. However, magnetic frustration is rarely observed in magnetoplumbites. Herein, the discovery, synthesis, and characterization of the first Mn-based magnetoplumbite, as well as the first magnetoplumbite involving pnictogens (Sb), ASb 3 Mn 9 O 19 (A = K or Rb) are reported. The Mn 3+ ( S = 2) cations, further confirmed by DC magnetic susceptibility and X-ray photoelectron spectroscopy, construct three geometrically frustrated sublattices, including Kagome, triangular, and puckered honeycomb lattices. Magnetic properties measurements revealed strong antiferromagnetic spin–spin coupling as well as multiple low-temperature magnetic features. Heat capacity data does not show any prominent λ-anomaly, suggesting minimal associated magnetic entropy. Moreover, neutron powder diffraction (NPD) implied the absence of long-range magnetic ordering in KSb 3 Mn 9 O 19 down to 3 K. However, several magnetic peaks are observed in RbSb 3 Mn 9 O 19 at 3 K, corresponding to an incommensurate magnetic structure. Interestingly, strong diffuse scattering is seen in the NPD patterns of both compounds at low angles and is analyzed by reverse Monte Carlo refinements, indicating short-range spin ordering related to frustrated magnetism as well as 2D magnetic correlations in ASb 3 Mn 9 O 19 (A = K or Rb).

2-D magnetic correlation↗
DOE OSTI2022

Magnetic excitations and interactions in the Kitaev hyperhoneycomb iridate β –Li 2 ⁢IrO 3

Here, we present a thorough experimental study of the three-dimensional hyperhoneycomb Kitaev magnet β –Li 2 ⁢IrO 3 , using a combination of inelastic neutron scattering (INS), time-domain terahertz spectroscopy (TDTS), and heat capacity measurements. The main results include a massive low-temperature reorganization of the INS spectral weight that evolves into a broad peak centered around 12 meV, and a distinctive peak in the terahertz data at 2.8(1) meV. A detailed comparison to powder-averaged spin-wave theory calculations reveals that the positions of these two features are controlled by the anisotropic Γ coupling and the Heisenberg exchange J, respectively. The refined microscopic spin model places β –Li 2 ⁢IrO 3 in close proximity to the Kitaev spin liquid phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Quenched disorder in the triangular lattice antiferromagnet YbZn 2⁢ GaO 5

We investigate the crystal electric field (CEF) excitations of Yb 3+ ions in powder samples of the triangular-lattice rare-earth-based antiferromagnet YbZn 2⁢ GaO 5 using inelastic neutron scattering (INS). Three CEF excitations from the ground-state Kramers doublet were observed, each exhibiting significant broadening beyond instrumental resolution. Combining temperature-dependent INS and neutron powder diffraction, we identify a significant static contribution to this broadening and attribute it to heterogeneous coordination of Yb 3+ ions due to Ga 3+ /Zn 2+ site mixing. Rietveld refinement of neutron powder diffraction indicates that 35% of Ga occupies the Zn site and 60% of Zn occupies the Ga site. We show with a point charge model for the CEF Hamiltonian that heterogeneous coordination of Yb3+ ions leads to broadened CEF peaks. First-principles calculations demonstrate that the random Ga 3+ /Zn 2+ distribution can produce the distortions of the YbO 6 octahedra observed from neutron diffraction. Because the documented heterogeneity will extend to exchange interactions, our results suggest that disorder is a significant factor in the unusual magnetism previously reported in YbZn 2 ⁢GaO 5 , including broad low-energy magnetic excitations and the absence of magnetic ordering down to 0.3 K.

Crystal field excitations↗
DOE OSTI2024

Four Years of Atmospheric Boundary Layer Height Retrievals Using COSMIC-2 Satellite Data

This work aimed to study the atmospheric boundary layer height (ABLH) from COSMIC-2 refractivity data, endeavoring to refine existing ABLH detection algorithms and scrutinize the resulting spatial and seasonal distributions. Through validation analyses involving different ground-based methodologies (involving data from lidar, ceilometer, microwave radiometers, and radiosondes), the optimal ABLH determination relied on identifying the lowest refractivity gradient negative peak with a magnitude at least $τ$% times the minimum refractivity gradient magnitude, where $τ$ is a fitting parameter representing the minimum peak strength relative to the absolute minimum refractivity gradient. Different $τ$ values were derived accounting for the moment of the day (daytime, nighttime, or sunrise/sunset) and the underlying surface (land or sea). Results show discernible relations between ABLH and various features, notably, the land cover and latitude. On average, ABLH is higher over oceans (≈1.5 km), but extreme values (maximums > 2.5 km, and minimums < 1 km) are reached over intertropical lands. Variability is generally subtle over oceans, whereas seasonality and daily evolution are pronounced over continents, with higher ABLHs during daytime and local wintertime (summertime) in intertropical (middle) latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Implementation of an Orificing Optimization Algorithm in the DASSH Subchannel Analysis Code

The Ducted Assembly Steady-State Heat transfer code (DASSH) performs full-core subchannel thermal hydraulics calculations in liquid metal fast reactors. One of the applications of subchannel codes is to optimize coolant flow orificing. As a design activity, the primary task is to determine the best way to divide assemblies into groups and distribute coolant flow rates among them. This report documents an algorithm implemented in DASSH to automatically optimize coolant orificing. Over the course of multiple iterations, DASSH determines the orifice grouping and flow distribution that minimizes peak coolant, clad, or fuel temperatures across all timesteps for a user-specified number of assembly groups. The total coolant flow rate in the reactor is constrained to achieve the specified core-average outlet temperature. The flow rate to each orifice group may also be constrained by the allowable pressure drop. The distribution of coolant flow among groups is accelerated using a predictor-corrector algorithm based on interpolated results from single-assembly parametric calculations. The assembly orificing grouping is initially predicted based on assembly power but can be refined if results demonstrate that an assembly would fit better in another group. The algorithm is demonstrated with two case studies. The first is a simple model for a reactor core consisting of just fuel assemblies; the pin power distributions are specified to create a situation where the initial assembly grouping prediction is suboptimal. This example is used to describe the initial grouping, demonstrate convergence over multiple iterations, and highlight the impact of regrouping. Then, the algorithm is applied to minimize peak clad and fuel temperatures in an example sodium-cooled fast reactor, the Versatile Test Reactor. The multicycle optimization confirms prior calculations for the reference core design. The example highlights how optimizing for different peak temperatures affects the results and demonstrates the use of the pressure drop constraint to limit the maximum flow rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

MILK : a Python scripting interface to MAUD for automation of Rietveld analysis

Modern diffraction experiments ( e.g. in situ parametric studies) present scientists with many diffraction patterns to analyze. Interactive analyses via graphical user interfaces tend to slow down obtaining quantitative results such as lattice parameters and phase fractions. Furthermore, Rietveld refinement strategies ( i.e. the parameter turn-on-off sequences) tend to be instrument specific or even specific to a given dataset, such that selection of strategies can become a bottleneck for efficient data analysis. Managing multi-histogram datasets such as from multi-bank neutron diffractometers or caked 2D synchrotron data presents additional challenges due to the large number of histogram-specific parameters. To overcome these challenges in the Rietveld software Material Analysis Using Diffraction ( MAUD ), the MAUD Interface Language Kit ( MILK ) is developed along with an updated text batch interface for MAUD . The open-source software MILK is computer-platform independent and is packaged as a Python library that interfaces with MAUD . Using MILK , model selection ( e.g. various texture or peak-broadening models), Rietveld parameter manipulation and distributed parallel batch computing can be performed through a high-level Python interface. A high-level interface enables analysis workflows to be easily programmed, shared and applied to large datasets, and external tools to be integrated with MAUD . Through modification to the MAUD batch interface, plot and data exports have been improved. The resulting hierarchical folders from Rietveld refinements with MILK are compatible with Cinema: Debye–Scherrer , a tool for visualizing and inspecting the results of multi-parameter analyses of large quantities of diffraction data. In this manuscript, the combined Python scripting and visualization capability of MILK is demonstrated with a quantitative texture and phase analysis of data collected at the HIPPO neutron diffractometer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2025

Harvesting Subsea Water Motion to Provide Clean Electrical Power

This TEAMER project had three main objectives related to characterizing the performance of a vortex-induced vibration (VIV) marine energy harvesting device called the Whatever Input To Torsion (WITT). The first was to characterize the motion of the WITT device as it vibrates on top of a pipe in a natural tidal flow. The second was to measure how much electricity the WITT device generates when it is deployed on top of a pipe that is secured to the seabed with a lander. The third objective was to measure how much electricity the WITT device generates in the same configuration but with a longer pipe. This TEAMER project succeeded in satisfying the first objective by measuring the motion of the WITT device with an inertial measurement unit (IMU) mounted in the WITT housing. The motion was described by examining changes in the acceleration, pitch, roll, and yaw, and displacement. The frequency of the motion was characterized using spectral analysis. The TEAMER project succeeded in reaching the second objective by measuring the power produced by the WITT. It is reported here as the average and maximum power in hourly intervals. WITT Energy decided to not have us pursue the third objective of testing a longer pipe because the first test was unexpectedly successful at having the system vibrate at the target frequency. A longer pipe would have resulted in a decreased frequency and therefore not have been a useful test. Based on these findings, the decision was made to keep the same pipe length for both tests and improve on other aspects of the test structure that had failed during the first test. The key finding of this TEAMER project is that the WITT is capable of generating electricity in moderate flow speeds when it is in a housing that is attached to a 2.5 m pipe and secured to the seabed by a lander. This report contains the first power values for a VIV instrument that generates electricity from a pendulum swinging at the top of a long pipe. The experiment demonstrated that the vibration frequency of the device varied with the tidal current speed, which changed throughout a given tidal cycle. This meant that the desired peak frequency was not sustained for more than one hour. The success of this TEAMER test provides evidence that continued research in the area of electricity generation from VIV holds promise as a source of marine energy for offshore applications. Future tests should focus on refining the pipe length equations and improving the power take off system of the WITT power electronics. In addition, experiments aimed at expanded testing and evaluation (e.g. longer than one month) are needed to measure the power production at higher tidal current speeds and test the system's durability.

Branch, Ruth A. (ORCID:0000000202356719)↗
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↗