Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source range”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Temperature dependence of nitrate-reducing Fe(II) oxidation by Acidovorax strain BoFeN1 – evaluating the role of enzymatic vs. abiotic Fe(II) oxidation by nitrite

ABSTRACT Fe(II) oxidation coupled to nitrate reduction is a widely observed metabolism. However, to what extent the observed Fe(II) oxidation is driven enzymatically or abiotically by metabolically produced nitrite remains puzzling. To distinguish between biotic and abiotic reactions, we cultivated the mixotrophic nitrate-reducing Fe(II)-oxidizing Acidovorax strain BoFeN1 over a wide range of temperatures and compared it to abiotic Fe(II) oxidation by nitrite at temperatures up to 60°C. The collected experimental data were subsequently analyzed through biogeochemical modeling. At 5°C, BoFeN1 cultures consumed acetate and reduced nitrate but did not significantly oxidize Fe(II). Abiotic Fe(II) oxidation by nitrite at different temperatures showed an Arrhenius-type behavior with an activation energy of 80±7 kJ/mol. Above 40°C, the kinetics of Fe(II) oxidation were abiotically driven, whereas at 30°C, where BoFeN1 can actively metabolize, the model-based interpretation strongly suggested that an enzymatic pathway was responsible for a large fraction (ca. 62%) of the oxidation. This result was reproduced even when no additional carbon source was present. Our results show that at below 30°C, i.e. at temperatures representing most natural environments, biological Fe(II) oxidation was largely responsible for overall Fe(II) oxidation, while abiotic Fe(II) oxidation by nitrite played a less important role.

Microbiology↗
DOE OSTI2025

Design and Study of Inductively Coupled Plasma Chamber Components Using the SupRISE Test Device at DIII-D

The DIII-D National Fusion Facility aims to increase the auxiliary heating power for the tokamak by upgrading the Neutral Beam Injection (NBI) system. In collaboration with North Carolina State University, the conventional ‘arc and filament’ NBI ion sources will be converted to inductively coupled plasma (ICP) sources which utilize radio frequency (RF) coupling to maximize reliability for high power operation. In support of this initiative, a full scale test device (SupRISE, Superior Radiofrequency Ion Source Experiment) is currently under construction at the DIII-D Facility. In preparation for the construction of a full scale prototype that can be installed on the DIII-D NBI system, experiments on SupRISE will be conducted to determine the optimal RF frequency for high power coupling, Faraday shield slit configuration, and ICP chamber length. SupRISE is comprised of an approximately 30 x 70 cm quartz dielectric vessel with an internal Faraday shield enclosed in a secondary vacuum chamber to ensure structural stability of the dielectric. Actively cooled front and back plates are designed to reduce the thermal stresses on the plasma facing components and mate with the existing accelerator used by the NBI system at DIII-D. 50 kW of RF power will be coupled to the plasma through the quartz over a variable frequency range of 4-8 MHz to sustain a plasma density of ∼ 10 18 m −3 for a 10 s ON, 210 s OFF duty cycle. Various modelling efforts have been employed to simulate the thermal and stress profiles over the primary components of the SupRISE device as well as the inductance behavior of the RF antenna. These simulation results and the final design for SupRISE will be presented. An additional reduced-scale predecessor ICP source (called RISE) has been used to develop a predictive match model that will be applied to frequency optimization studies on SupRISE. Furthermore, the outcomes of this research and complementary efforts at North Carolina State University are essential for the incorporation of ICP NBI positive ion sources at the DIII-D facility.

DIII-D↗
DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗
DOE OSTI2025

Analysis of measurements from an array of radioxenon samplers near to Hartlepool Nuclear Power Station

As part of a scientific research and development project, the radionuclide fingerprint of an operating advanced gas-cooled reactor 25 (AGR) has been studied across several facets (Goodwin et al., 2024). One part of this project was to deploy an array of 26 radioxenon samplers to the region for a period of around 1 year, to measure any radioxenon emissions from the reactors of the 27 Hartlepool nuclear power station at a range of tens of kilometres away. The array of 3 sensors was operational for around 12 28 months from March 2022 and detected many occurrences of isotopes of radioxenon. Here we provide a detailed analysis and 29 interpretation of the data and where possible, attribution of detections to a source or region. Whilst a key part of this work is to 30 measure any emissions of radioxenon from Hartlepool, this work presents one of the most comprehensive efforts to determine the 31 source of a great number of (mostly) 133Xe detections. A combination of different types of atmospheric dispersion modelling 32 techniques, including the use of stack monitoring data from nearby civil radioxenon-emitting nuclear facilities has enabled the 33 majority of detections to be attributed to one or more possible sources. Whilst emissions from Hartlepool have been detected on 34 the systems, the majority of detections are associated with a medical isotope production facility in Fleurus, Belgium (IRE).

Advanced gas-cooled reactor↗
DOE OSTI2023

Determining the influence of temperature and pressure on the structural stability in a polyurea elastomer

The structural response of amorphous elastomers to applied pressure and temperature is of significant interest given their potential for blast and impact protection applications. While understanding disordered structures is often challenging, this study focuses on the ordered correlations that are present at the atomic level even in amorphous materials. This local order was characterized in polyurea elastomers with two molecular weights up to 6 GPa and ~280 °C using in situ energy dispersive X-ray diffraction at the Advanced Photon Source via total scattering techniques. The resulting structure factors and pair distribution functions indicate that the applied pressure reduced the free volume but had little effect on covalent bonds for both molecular weights. Temperature-induced order-disorder transitions that occur readily at ambient pressures within the temperature range studied here were not observed, their onset likely delayed by the limited chain mobility associated with the free volume reduction. In conclusion, this study demonstrates that total scattering offers insights into the high-pressure structure of disordered polymers that are otherwise inaccessible using more conventional X-ray scattering techniques.

42 ENGINEERING↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2025

Energy transfer between localized emitters in photonic cavities from first principles

Radiative and nonradiative resonant couplings between defects are ubiquitous phenomena in photonic devices used in classical and quantum information technology applications. In this work, we present a first-principles approach to enable quantitative predictions of the energy transfer between defects in photonic cavities, beyond the dipole-dipole approximation and including the many-body nature of the electronic states. As an example, we discuss the energy transfer from a dipolelike emitter to an 𝐹 center in MgO in a spherical cavity. We show that the cavity can be used to controllably enhance or suppress specific spin-flip and spin-conserving transitions. Specifically, we predict that an ∼10–100 enhancement in the resonant energy transfer rate can be gained in the case of the 𝐹 center in MgO at ∼10 nm distances from a dipolar source, using rather moderate cavity with quality factor 𝑄 ∼ 400. We also show that a similar suppression in the transfer rate can be achieved by off-tuning the cavity resonance relative to the emitter transition energy. The framework presented here is general and readily applicable to a wide range of devices where localized emitters are embedded in microspheres, core-shell nanoparticles, and dielectric Mie resonators. Hence, our approach paves the way to predict how to control energy transfer in quantum memories and in ultrahigh-density optical memories, and in a variety of quantum information platforms.

First-principles calculations↗
DOE OSTI2026

Formation of niobium hydride precipitates in superconducting qubits

We report evidence for the formation of niobium hydride phase within niobium films on silicon substrates in superconducting qubits fabricated at Rigetti Computing. For this study, we combined complementary techniques—including room-temperature and cryogenic atomic force microscopy (AFM), synchrotron x-ray diffraction, and time-of-flight secondary ion mass spectroscopy (ToF-SIMS)—to directly reveal the existence of niobium hydride precipitates on the surface of superconducting qubits. Upon cryogenic cooling, we observed variation in the size and morphology of the hydrides, ranging from small (∼5 nm) irregular shapes to large (∼10–100 nm) domains within the Nb grains, which were fully converted to niobium hydrides. Since niobium hydrides are nonsuperconducting and can easily change in size and location upon different cooldowns to cryogenic temperature, our finding highlights a previously unknown source of decoherence in superconducting qubits. This contributes to quasiparticle losses, offering a potential explanation for changes in qubit performance upon cooldowns. Finally, by leveraging the RF performance of a 3D bulk Nb resonator, we quantify RF dissipation in a superconducting qubit caused by hydrogen concentration variation, and propose a practical engineering pathway to mitigate the formation of Nb hydrides for superconducting qubit applications.

Sung, Zuhawn [Fermi National Accelerator Laborator↗
DOE OSTI2024

Reviving QFT in ( 2 + 1 )-dimensional de Sitter spacetime

We consider a conformally coupled scalar quantum field theory (QFT) on ( 2 + 1 )-dimensional static Einstein universe R × S 2 and write down the free theory Hilbert space. We explain that this theory is secretly a QFT in ( 2 + 1 )-dimensional de Sitter space because all the quantum observables experience “quantum revivals,” which naturally restricts the timelike R to the appropriate de Sitter time range. Our construction circumvents the causal obstruction to formulating QFT in de Sitter due to event horizons. There are not any in static Einstein. The “unitary gauge” description of the theory is realized by the zonal harmonics P ℓ ( n ^ · n ^ ′ ) . We verify that interactions with conformally invariant external sources are mediated only by these modes. Hence these modes comprise the complete basis of the “bulk” theory. When the theory is cut off in the UV, the basis dimension scales as the Bekenstein-Hawking formula. Published by the American Physical Society 2024

D’Amico, Guido (ORCID:0000000281831214)↗
DOE OSTI2020

Creating a Detector Response Matrix for the Micro-Calorimeter

To better understand the results from the micro calorimeter, a detector response matrix (DRM) was generated. This will allow to predict source responses, which in turn will improve the confidence in the instrument. To generate a DRM, the simulation package Rapid Adaptable Multi-threaded Particle and Radiation Transport Simulation (RAMPART) that utilizes underlying Geant4 libraries for particle transport was used to generate the simulation. A flat energy input spectrum ranging from 0 keV to 3000 keV was used to generate the response. 37 x 10 9 particles were simulated. In this report, the details and findings are presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Language Independent Static Analysis (LISA)

Software is becoming increasingly important in nearly every aspect of global society and therefore in nearly every aspect of national security as well. While there have been major advancements in recent years in formally proving properties of program source code during development, such approaches are still in the minority among development teams, and the vast majority of code in this software explosion is produced without such properties. In these cases, the source code must be analyzed in order to establish whether the properties of interest hold. Because of the volume of software being produced, automated approaches to software analysis are necessary to meet the need. However, this software boom is not occurring in just one language. There are a wide range of languages of interest in national security spaces, including well-known languages such as C, C++, Python, Java, Javascript, and many more. But recent years have produced a wide range of new languages, including Nim, (2008), Go (2009), Rust (2010), Dart (2011), Kotlin (2011), Elixir (2011), Red (2011), Julia (2012), Typescript (2012), Swift (2014), Hack (2014), Crystal (2014), Ballerina (2017) and more. Historically, automated software analyses are implemented as tools that intermingle both the analysis question at hand with target language dependencies throughout their code, making re-use of components for different analysis questions or different target languages impractical. This project seeks to explore how mission-relevant, static software analyses can be designed and constructed in a language-independent fashion, dramatically increasing the reusability of software analysis investments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Long Range Community Planning Process for Fusion Sciences and Plasma Physics (Final Technical Report)

This grant facilitated the APS-DPP Community Planning Process (DPP-CPP) which was a strategic planning process with the goals of identifying scientific and technological opportunities in the fields of Plasma Physics and Fusion Energy Science, and making consensus recommendations for a strategy to address these opportunities to the DOE Office of Fusion Energy Science Advisory Committee (FESAC). The DPP-CPP was initiated by the APS-DPP executive committee in order to help FESAC fulfill a charge from DOE to develop a strategic plan for the DOE Office of Fusion Energy Science (FES). In this charge, the DOE Deputy Director for Science, Dr. Stephen Binkley, requests “that the Fusion Energy Sciences Advisory Committee (FESAC) undertake a new long-range strategic planning activity for the Fusion Energy Sciences (FES) program. The strategic planning activity - to encompass the entire FES research portfolio (namely, burning plasma science and discovery plasma science) - should identify and prioritize the research required to advance both the scientific foundation needed to develop a fusion energy source, as well as the broader FES mission to steward plasma science.” The process as set out in the proposal (and below) was successful in meeting all of its stated goals culminating in the final report to the FESAC subcommittee given the charge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Gambusia holbrooki Survive Shear Stress, Pressurization and Avoid Blade Strike in a Simulated Pumped Hydroelectric Scheme

Pumped hydroelectric energy storage (PHES) projects are being considered worldwide as a means of achieving political renewable energy targets in a way that stabilises baseload energy supply from often intermittent renewable energy sources. Unlike a conventional hydroelectric system that only pass water in a downstream direction, a feature of PHES is that it relies on the bi-directional flow of water. In some cases this flow can be across different waterbodies or catchments, posing a risk of inadvertently expanding the range of aquatic animals like fish. The risk of this happening depends on the likelihood of survival of individuals, which remains poorly understood for turbines that are pumping rather than generating. This study quantified the survival of a globally widespread and invasive poeciliid fish, Eastern gambusia (Gambusia holbrooki) when exposed hydraulic stresses characteristic of what would be experienced through a PHES during the pumping phase. A shear flume and hyperbaric chamber were used to expose fish to different strain rates and rapid and sustained pressurisation. A blade strike model was also used to predict survival of fish passing through a Francis dual turbine / pump. The ranges simulated were based on design and operational conditions provided for a PHES scheme being proposed in south-eastern Australia. All gambusia tested survived extremely high (up to 7600 kPa gauge pressure) pressurisation, high levels of shear stress (up to 1853 s -1 ), and the majority (> 93 %) were unlikely to be struck by a turbine blade. Given their tolerance to these extreme simulated stresses, we conclude that gambusia will likely survive passage through a PHES scheme. Therefore, where a new PHES poses the risk of inadvertently expanding the range of gambusia or similar poeciliid species, measures to minimise their spread or mitigate their ecosystem impacts should be considered.

25 ENERGY STORAGE↗
DOE OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Reacceleration of Galactic Cosmic Rays beyond the Knee at the Termination Shock of a Cosmic-Ray-driven Galactic Wind

The origin of cosmic rays (CRs) above the knee in the spectrum is an unsolved problem. We present a wind model in which interstellar gas flows along a nonrotating, expanding flux tube with a changing speed and cross-sectional area. CRs from Galactic sources, such as supernova remnants, which are coupled to the plasma via Alfvén waves, provide the main pressure source for driving this outflow. These CRs are then subject to diffusive shock reacceleration at the Galactic wind termination shock, which is located at a distance ~200 kpc. Some of the highest-energy reaccelerated particles propagate upstream against the wind and can contribute to the petaelectronvolt to exaelectronvolt range of the spectrum. We analyze the conditions under which efficient reacceleration can occur and find that rigidities ~10–40 PV can be obtained and that the termination shock may account for between 10% and 50% of the proton spectrum measured in IceCube/IceTop experiment. In our model, the termination shock is unable to fully explain the CR spectrum in the petaelectronvolt to exaelectronvolt range. The highest-energy particles that escape downstream from our termination shock, and similar shocks surrounding most galaxies, can be further accelerated by intergalactic shock fronts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magellan/M2FS and MMT/Hectochelle Spectroscopy of Dwarf Galaxies and Faint Star Clusters within the Galactic Halo*

We present spectroscopic data for 16,369 stellar targets within and/or toward 38 dwarf spheroidal galaxies and faint star clusters within the Milky Way halo environment. All spectra come from observations with the multiobject, fiber-fed echelle spectrographs M2FS at the Magellan/Clay telescope or Hectochelle at the MMT, reaching a typical limiting magnitude G ≲ 21. Data products include processed spectra from all observations and catalogs listing estimates—derived from template model fitting—of line-of-sight velocity (median uncertainty 1.4 km s –1 ) effective temperature (255 K), (base-10 logarithm of) surface gravity (0.59 dex in cgs units), [Fe/H] (0.4 dex) and [Mg/Fe] (0.27 dex) abundance ratios. The sample contains multiepoch measurements for 3720 sources, with up to 15 epochs per source, enabling studies of intrinsic spectroscopic variability. The sample contains 6087 likely red giant stars (based on surface gravity), and 4492 likely members (based on line-of-sight velocity and Gaia-measured proper motion) of the target systems. The number of member stars per individual target system ranges from a few, for the faintest systems, to ~850 for the most luminous. For most systems, our new samples extend over wider fields than have previously been observed; of the likely members in our samples, 820 lie beyond 2 times the projected half-light radius of their host system, and 42 lie beyond 5 R half .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Measurement of Dielectric Loss in Silicon Nitride at Centimeter and Millimeter Wavelengths

This work presents a suite of measurement techniques for characterizing the dielectric loss tangent across a wide frequency range from similar to 1 GHz to 150 GHz using the same test chip. In the first method, we fit data from a microwave resonator at different temperatures to a model that captures the two-level system (TLS) response to extract and characterize both the real and imaginary components of the dielectric loss. The inverse of the internal quality factor is a second measure of the overall loss of the resonator, where TLS loss through the dielectric material is typically the dominant source. The third technique is a differential optical measurement at 150 GHz. The same antenna feeds two microstrip lines with different lengths that terminate in two microwave kinetic inductance detectors (MKIDs). The difference in the detector response is used to estimate the loss per unit length of the microstrip line. Our results suggest a larger loss for SiNx at 150 GHz of tan delta similar to 4 x 10(-3) compared to 2.0 x 10(-3) and greater than or similar to 1 x 10(-3) measured at similar to 1 GHz using the other two methods. These measurement techniques can be applied to other dielectrics by adjusting the microstrip lengths to provide enough optical efficiency contrast and other mm/sub-mm frequency ranges by tuning the antenna and feedhorn accordingly.

Dielectrics loss↗
DOE OSTI2023

Numerical relativity estimates of the remnant recoil velocity in binary neutron star mergers

We present, for the first time, recoil velocity estimates for binary neutron star mergers using data from numerical relativity simulations. We find that binary neutron star merger remnants can have recoil velocity of the order of a few tens of kilometers per second and as high as 150 kilometers per second in our dataset. These recoils are attained due to equivalent contributions from the anisotropic gravitational wave emission as well as the asymmetric ejection of dynamical matter during the merger. Here, we provide fits for net recoil velocity as well as its ejecta component as a function of the amount of ejected matter, which may be useful when constraints on the ejected matter are obtained through electromagnetic observations. We also estimate the mass and spin of the remnants and find them to be in the range [2.34,3.38]M ⊙ and [0.63, 0.82], respectively, for our dataset.

79 ASTRONOMY AND ASTROPHYSICS↗