SEARCH · Engineering Papers
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.
Here, the Chebyshev Rational Approximation Method (CRAM) algorithm has become the method of choice for modern nuclear depletion calculations. A simple modification of the CRAM yields not only depleted nuclide amounts but also time-averaged nuclide amounts, with little additional calculation cost.
Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.
Time-Averaged "Slow" Carbon Dioxide Flux Data from wind and gas files, with additional diagnostic data: 22m samples
Time-Averaged "Slow" Carbon Dioxide Flux Data from wind and gas files, with additional diagnostic data: 40m samples
Time-Averaged "Slow" Carbon Dioxide Flux Data from wind and gas files, with additional diagnostic data: 45m samples
This dataset contains microseismic data acquired during the Frisco pad flow test project led by Fervo Energy, conducted between July 17th - Aug 12th 2024, near the Utah FORGE geothermal site. The microseismic data was collected from various Utah FORGE wells: via Distributed Acoustic Sensing (DAS) fiber in 16B, and two 3-component geophones located in wells 56-32 and 78B. The dataset is structured in SEGY format, where the first six traces represent data from the geophones, and the remaining traces capture DAS data from well 16B. Each SEGY file in this dataset contains triggered microseismic events, with event initiation based on Short-Time Average over Long-Time Average (STA/LTA) detection criteria during the stimulation process. Files are grouped by time intervals and named following the structure "[WellPad][WellName][Month]_[Year]Divine_Trigger[EventNumber].sgy," indicating well pad, well name, date, and event number. Sampling parameters include a spatial sampling of approximately 2 meters for DAS channels and a temporal sampling rate of 2000 Hz, with each data record spanning 1.2 seconds. The files' coordinates are referenced to the location of the FORGE 16A-32 wellhead, positioned at UTM coordinates: Easting 334641.1891 m and Northing 4263443.693 m. The geographic coordinates for this origin are 38.50402147 latitude and -112.8963897 longitude, with an elevation of 1650.0249 meters above sea level.
Many models of the atmosphere simulate vertical profiles of observable fields, such as temperature, moisture, and vertical velocity. To assess the realism of such model output, it is important to evaluate the modeled profiles against observations. However, many extant evaluations consist of comparing time averaged profiles from the model and the observations. Such time averages can smooth out unrealistic features that appear in instantaneous profiles. In this way, the time averages may obscure certain model errors from detection. For instance, many models misrepresent turbulent processes, leading to profiles that are too smoothed out or too noisy. E.g., if a model has coarse vertical resolution and overly strong diffusion, then the model may be unable to produce observed sharp peaks. On the other hand, a model may have artificial approximations that produce spurious jumps in the vertical. Much can be learned by browsing through instantaneous snapshots of model output. But the volume of output produced by a global model is far too great to rely purely on visual inspection. In addition, we need a way to statistically characterize the profiles, so that anomalies can be detected more conveniently.
Atmospheric processes over the Southern Ocean have a profound influence on regional and global climate. This is a part of the world where global climate models perform particularly poorly, with models persistently overpredicting the amount of sunlight reaching the Earth's surface (Regayre et al. 2020). A major challenge when trying to improve the representation of atmospheric processes in this part of the globe is the scarcity of observations capable of constraining our understanding. During 2024/25 a U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility field campaign (CAPE-k) deployed a suite of instrumentation at the kennaook/Cape Grim atmospheric monitoring station in Tasmania, Australia to provide detailed cloud-aerosol observations in this region in order to address the greatest source of uncertainty in climate models – atmospheric aerosols and how they influence cloud formation (Carslaw et al. 2013). In support of the ARM observations, an Australian Research Council (ARC)-funded project “Southern Ocean aerosols: sources, sinks and impact on cloud properties,” led by the Queensland University of Technology, complemented the ARM measurements with a suite of chemical measurements at the kennaook/Cape Grim site in northwestern Tasmania to provide more information on the processes controlling aerosol formation and growth. The deployment of the University of York LIF-SO2 instrument (Temple et al. 2025) in the ARM mobile facility container was part of this chemistry-focused deployment. Increased observational efforts are critical for understanding sources and sinks of Southern Ocean aerosols, in particular the role of marine micro-organisms (phytoplankton, algae) on aerosols formation, their properties and growth, cloud droplet formation, and the conversion of cloud droplets to ice crystals and precipitation. The observational focus of the ARC-funded project is aerosol chemical composition and their gaseous precursors. The remote nature of the Southern Ocean poses a significant analytical challenge when studying gaseous precursors, as low concentrations are often below the limits of detection of commercial instrumentation. Historically this has meant that observations have had to be time-averaged over days to weeks to achieve the limits of detection required. Although useful when considering overall average levels of aerosol precursors, this time-averaging obscures temporal variability that can provide insight into the controlling processes. Over recent years advances in analytical techniques have made high-time-resolution measurements of key gas-phase aerosol precursors achievable with sufficiently low limits of detection.
TiSe 2 is a narrow-gap insulator with a rich array of unique properties. In addition to being a superconductor under certain modifications, it is commonly thought to be a rare realisation of an excitonic insulator. Below 200 K, TiSe 2 undergoes a transition from a high-symmetry ($P\bar{3}m1$) phase to a low-symmetry ($P\bar{3}c1$) charge density wave (CDW). Here we establish that it is indeed an insulator in both $P\bar{3}m1$ and $P\bar{3}c1$ phases. However, the insulating state is driven not by excitonic effects but by symmetry-breaking. In the CDW phase it is static. At high temperature, thermally driven instantaneous deviations from $P\bar{3}m1$ break the symmetry on the characteristic time scale of a phonon. Even though the time-averaged lattice structure assumes $P\bar{3}m1$ symmetry, the time-averaged energy band structure is closer to the CDW phase - a rare instance of a metal-insulator transition induced by dynamical symmetry breaking. We establish these conclusions from quasiparticle self-consistent GW (QSGW) and many-body calculations (QS$G\widehat{W}$), in combination with molecular dynamics simulations to capture the effects of thermal disorder. The many-body theory includes explicitly ladder diagrams in the polarizability, which incorporates excitonic effects in an ab initio manner. We find that the excitonic modification to the potential is weak, ruling out the possibility that TiSe 2 is an excitonic insulator.
The Contractor will build a liquid hydrogen (LH2) fueling station model that is capable of simulating the amount of hydrogen transferred from the liquid storage tank to the vehicle, depending on the station’s design and control scheme. This modeling work will allow the Parties to seek the design and control scheme of an LH2 station that can fill a targeted number of vehicles while meeting a time-average flow rate of 10 kg/min for each fill. While evaluating a high dispensing performance station design and control scheme, the Contractor will analyze the design and control scheme are economically viable.
Per- and polyfluoroalkyl substances (PFAS) are a group of industrial and commercial chemicals widely used throughout the world due to their beneficial chemical properties. Because of their widespread use, their chemical stability, and their ability to be transported over long distances through atmospheric deposition and movement through waterways, PFAS are found throughout most aquatic ecosystems; yet large sampling gaps exist among reservoir and river ecosystems in the desert southwest of the United States. In this study, we examine PFAS concentrations in the tissue of fish (catfish [channel and blue], common carp, smallmouth bass, northern pike, walleye, white crappie and white sucker) collected in northern New Mexico, including examining PFAS composition and concentration relative to trophic level distribution. We collected fish from two man-made reservoirs and from the Rio Grande. We then collected muscle and liver tissues from fish specimens, which were screened for 39 PFAS compounds. We detected PFAS compounds in most fish tissue sampled, including the biomagnification of PFAS compounds within liver samples, with PFOS concentrations ranged from 1.13 to 350.1 (64.4 average) times higher in the liver samples compared to muscle samples. Most PFAS concentrations within muscle samples were within the range of atmospheric transportation previously reported and average tissue concentrations of PFAS were calculated to be 2.02 ± 1.81 ng g -1 . Using stable isotopes as a predictor of trophic-foraging exposure and PFAS concentrations, we noted a correlation between enriched δ 15 N values, which had higher perfluorodecanoic acid concentrations.
A privacy-preserving average consensus algorithm is designed based on the Beaver triple technique against passive adversaries. The Beaver triple technique is integrated into a restructure of the discrete-time average consensus algorithm to preserve the privacy of initial values of agents in a multiagent system. The performance of the algorithm is theoretically analyzed.
X-ray reverberation, which exploits the time delays between variability in different energy bands as a function of Fourier frequency, probes the structure of the inner accretion disks and X-ray coronae of active galactic nuclei. We present a systematic X-ray spectroscopic and reverberation study of the high-Eddington-ratio narrow-line Seyfert 1 galaxy Ark 564, using over 900 ks of XMM-Newton and NuSTAR observations spanning 13 yr. The time-averaged spectra can be well described by a model consisting of a coronal continuum, relativistic disk reflection, warm Comptonization, and three warm absorbers. Leveraging the high X-ray brightness of Ark 564, we are able to resolve the time evolution of the spectra and contemporaneous reverberation lags. The soft-band lag relative to the continuum increases with the X-ray flux, while Fe Kα lags are detected in only a subset of epochs and do not correlate with soft lags. Models based on a lamppost corona and reflection from a standard thin disk can broadly reproduce the observed lag-energy spectra of low-flux epochs; however, additional reverberation from the warm Comptonized atmosphere is required to explain the soft lags observed in high-flux epochs. A vertically puffed-up inner disk and a variable, vertically extended corona can better explain the observed evolution of the lags and covariance spectra. Our study underscores the importance of multiepoch, multiband analyses for a comprehensive understanding the inner accretion disk and corona.
The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.
Materials improvement is a powerful approach to reducing loss and decoherence in superconducting qubits, because such improvements can be readily translated to large-scale processors. Recent work improved transmon coherence by using tantalum as a base layer and sapphire as a substrate. The losses in these devices are dominated by two-level systems with comparable contributions from both the surface and bulk dielectrics, indicating that both must be tackled to achieve substantial improvements in the state of the art. Here we show that replacing the substrate with high-resistivity silicon markedly decreases the bulk substrate loss, enabling 2D transmons with time-averaged quality factors (Q avg ) of 9.7 × 10 6 across 45 qubits. For our best qubit, we achieve a Q avg of 1.5 × 10 7 , reaching a maximum Q of 2.5 × 10 7 , corresponding to a lifetime (T 1 ) up to 1.68 ms. This low loss also allows us to observe decoherence effects related to the Josephson junction, and we use an improved, low-contamination junction deposition to achieve Hahn echo coherence times (T 2E ) exceeding T 1 . We achieve these materials improvements without any modifications to the qubit architecture, allowing us to readily incorporate standard quantum control gates. Here, we demonstrate single-qubit gates with 99.994% fidelity. The tantalum-on-silicon platform comprises a simple material stack that can potentially be fabricated at the wafer scale and therefore can be readily translated to large-scale quantum processors.
Since the turn of the millennium, the United States (U.S.) oil and natural gas (ONG) industry has nearly doubled its natural gas production rate. As a result, the ONG industry has recently come under increasing scrutiny for its contributions to greenhouse gas (GHG) emissions. Consequently, various solutions to this problem have been proposed and formulated to reduce the impacts of GHG emissions on the environment. West Virginia University (WVU) have found it important to research the impacts of recovering vented gas streams into prime-mover engines. The U.S. Department of Energy (DOE) and National Energy Technology Laboratory (NETL) have granted WVU funding to research and develop a “Methane Mitigator” (M2) - a “Scalable Vent Mitigation Strategy to Simultaneously Reduce Methane Emissions and Fuel Consumption from the Compression Industry.” One of the main areas of interest for this research was the collection of emissions from natural gas equipment into a Caterpillar G3508J natural gas compression engine. The parameters being analyzed from the engine were brake-specific emissions and power output. The emissions sources considered for this research were pneumatic controllers (PCs), reciprocating compressor vents, and the engine’s open crankcase breather. The compressor vent and PC emissions were simulated using a mass flow controller (MFC) and flowed into the engine using two separate methods: (1) directly into the air intake, and (2) through a retrofitted closed crankcase ventilation system (CCV), serving as a buffer volume. The crankcase emissions were quantified without the CCV, and the impact on exhaust emissions from circulating the crankcase gases into the intake was measured. The simulated compressor vent and PC flows from the MFC had limited effect on the steady state operation of the engine and resulting performance. When the simulated flows were fed directly into the engine’s air intake, the changes within the engine’s continuous performance and emission parameters were larger but lasted for shorter durations. Conversely, when the simulated flows were fed into the CCV before entering the air intake, the changes in the engine’s performance and emission parameters were less pronounced for continuous analysis but lasted for longer durations. In either case, the continuous emission changes in both emissions and performance varied in size depending on the test scenario being run, but the cycle average changes in emissions and performance showed little impact overall compared to the engine’s baseline operation. As a result, the inclusion of a CCV shows a decrease in baseline carbon dioxide equivalent (CO2-eq.) engine emissions (from combined exhaust and open crankcase) of almost 4%. Likewise, the CCV inclusion reduced baseline total methane (CH4) from combined exhaust and open crankcase by upwards of 16%. These atmospheric emissions only decreased further with the inclusions of collected PC and compressor vent flows. The resulting changes in time-averaged rated exhaust behavior (or lack thereof) prove that the proposed M2 system could likely be deployed at sites with modern lean-burn natural gas engines as a viable option for reducing and eliminating potential GHG sources that would have otherwise been unutilized as energy sources.
This document presents the as-run analysis of the Advanced Gas Reactor (AGR)-5/6/7 irradiation experiment. AGR-5/6/7 is the last of a series of experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory in support of the development and qualification of tri-structural isotropic low-enriched fuel for use in high-temperature gas-cooled reactors. The test train contained five separate capsules that were independently controlled and monitored. Each capsule contained multiple 24.91-mm-long and 12.25-mm-dimeter compacts filled with low-enriched uranium carbide/oxide tri-structural isotropic fuel particles. The objectives of the AGR-5/6/7 experiment were to: • Irradiate reference-design fuel particles to support fuel qualification. • Establish operating margins for the fuel, beyond normal operating conditions. • Provide irradiated-fuel performance data and irradiated-fuel samples for post-irradiation examination and safety testing. The primary objective of the AGR-5/6 test (Capsules 1, 2, 4, and 5) was to verify the successful performance of the reference-design fuel under normal operating conditions. The AGR-7 test (Capsule 3) was designed to explore fuel performance at higher temperatures. Its primary objective was to demonstrate the capability of the fuel to withstand conditions beyond normal operating conditions, in support of plant design and licensing. AGR-5/6/7 will also provide irradiated-fuel performance data based on the fission gas release from particles during irradiation. To achieve the test objectives, the AGR-5/6/7 experiment was irradiated in the northeast flux trap of the ATR with a planned duration of 500 effective full-power days. The northeast flux trap was selected because its larger diameter provided greater flexibility for test-train design compared to the Large B positions used for the AGR-1 and AGR-2 irradiations, significantly enhancing test capabilities for the combined irradiation campaigns. Due to delays in the ATR schedule, the AGR-5/6/7 irradiation was significantly shorter than the originally planned 13-cycle schedule. Irradiation began on February 16, 2018 and ended on July 22, 2020, spanning nine ATR cycles (162B–168A) over two and a half years. Thus, the AGR-5/6/7 fuel compacts were irradiated for a total of approximately 360.9 effective full-power days. Final burnup values, on a per-compact basis, ranged from 5.66 to 15.26% fissions per initial heavy metal atom, while fast fluence values ranged from 1.62 to 5.55 × 1025 n/m2 (E >0.18 MeV). Time-averaged volume-averaged fuel temperatures on a capsule basis at the end of irradiation ranged from 756°C in Capsule 5 to 1313°C in Capsule 3 excluding days with significantly lower temperature during the two short powered axial locator mechanism cycles, 163A and 167A. By the end of irradiation, 48 out of 54 installed thermocouples had failed (the bottom three capsules lost all thermocouples). During the first five cycles (162B – 165A), the fission-gas isotope release-rate-to-birthrate (R/B) ratios were stable in the 10-8–10-6 range, and no in-pile particle failures were observed based on the gross gamma counts. During this time, the high exposed kernel fraction and high fuel particle temperatures in Capsule 1 led to a maximum R/B value of around 2 ? 10-6 for Kr-85m. The fission gas release in all capsules started to increase from the second half of Cycle 166A, when a large number of in-pile particle failures occurred in Capsule 1 and a gas line problem in this capsule caused fission gas leakage at various degrees into the other four capsules. This gas line problem also prevented a fission gas release measurement for Capsule 1 during the last three cycles due to gas flow isolation. By the end of irradiation, it is estimated that approximately 15 particles failed in Capsule 3, which was considered possible because the experiment was designed to operate beyond the high-temperature gas-cooled reactor normal operating temperature range. A few hundred in-pile particle failures were estimated for Capsule 1 by the end of Cycle 166A, but the total number of failures is unknown due to the lack of fission gas release data in the later cycles. Additionally, four potential in-pile failures were identified for Capsule 2 during the last cycle, Cycle 168A. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on the absence of the typical spikes in gross gamma counts and low failure estimates using the AGR-3/4 R/B per exposed kernel model.
The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.