Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field enhancement factor”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2023

Tracking Aerosol Convection Interactions Experiment (TRACER) Field Campaign Report

Convective clouds serve a critical role in the Earth’s energy and water cycles through their transport of heat, moisture, momentum, and chemical species through the troposphere driving the global circulation (e.g., Hartmann et al. 1984, Del Genio et al. 2012, Su et al. 2014). On more local scales, convective clouds impact the atmospheric heating profile through diabatic heating effects, removal of water from the atmospheric column through precipitation, and conditioning of the local environment impacting further development of clouds (e.g., Sullivan and Voigt 2021). These critical roles underscore the importance of realistic representation of convective processes across scales of models from large-eddy simulation (LES), to convection-permitting models (CPM; e.g., Kendon et al. 2020, Marinescu et al. 2021), to numerical weather prediction (NWP) models used for operational weather forecasting, to Earth system models used to predict climate sensitivity (Sanderson et al. 2011, Sherwood et al. 2014, Tomassini et al. 2014, Zhao et al. 2016, Cronin et al. 2017). A key component of improving model representation of convective clouds is better quantification and parameterization of updraft microphysics and dynamics, including their interactions with the surrounding environment and storm organization (Bony et al. 2015, Hagos and Houze 2016, Donner et al. 2016, Morrison et al. 2020). Aerosol is an important environmental factor that could affect convective clouds and precipitation since cloud droplet and ice formation processes are initiated by it. Andrae et al. (2004) hypothesized that aerosols associated with increased biomass burning particles acting as cloud condensation nuclei (CCN) result in smaller and more monodisperse cloud droplets leading to suppression of warm rain formation, ultimately leading to more cloud water being lofted above the freezing level based on observations in the Amazon region. The subsequent increase in latent heat release increases the buoyancy of rising convective parcels invigorating the deep convection. This work was followed by a description of the theoretical basis for this “cold-phase invigoration” by Rosenfeld et al. (2008), who argued that it could have a significant effect for deep convective clouds with warm cloud-bases. Several modeling studies (e.g., Khain et al. 2005, 2009, van den Heever et al. 2006, Fan et al. 2007, 2009, 2012, Lee et al. 2008, Storer et al. 2010, Lebo et al. 2012, Storer and van den Heever 2013, Chen et al. 2020, Dagan et al. 2022) have investigated these aerosol-convection interactions and the environmental factors that influence their relative importance and magnitude. More recently, several studies have indicated that “warm-phase invigoration”, the enhancement of convection through condensational heating, also appears to play a role in enhancing both shallow cumuli (Seiki and Nakajima 2014, Saleeby et al 2015) and deeper tropical convection (Lebo and Seinfeld 2011, Khain et al. 2012, Sheffield et al 2015, Fan et al. 2018, Igel and van den Heever 2021), as well as Houston thunderstorms (Fan et al. 2007, 2020). However, still other studies have provided additional evidence of systematic biases in simulated convective outflow ice size distribution properties, which are consistent with a lack of poorly understood secondary ice production within convective updrafts (e.g., Fridlind et al. 2017). To help address these critical gaps in our understanding of cloud processes, aerosol processes and aerosol-cloud interactions, the Tracking Aerosol Convection Interactions Experiment was designed building upon efforts by the Aerosol, Cloud, Precipitation and Climate (ACPC) Initiative (http://acpcintiative.org/), a joint effort of the International Geosphere-Biosphere Programme (IGBP) and the World Climate Research Program (WCRP) that focused on resolving uncertainties in the interactions between aerosol and clouds towards better understanding the role that these interactions play in the climate system. The TRACER campaign was motivated by recommendations from a number of pilot studies undertaken by ACPC (van den Heever et al. 2017, Fridlind et al. 2019, Hu et al. 2019, Fan et al. 2020, Marinescu et al. 2021, Hernandez-Deckers et al. 2022) that pointed towards the southeastern Texas region as a locale where aerosol-convection interactions could be studied owing to the copious occurrence of isolated convection during the summer months accompanied by diverse and significant sources of aerosols from both anthropogenic and natural sources. The TRACER campaign began on 01 October 2021 and extended through 30 September 2022 with an intensive operational period (IOP) during June-September 2022. Three main sites (Table 1) were managed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

Improving the Visible and Infrared Contrast Ratio of Microshutter Arrays

Three device improvements have been developed that dramatically enhance the contrast ratio of microshutters. The goal of a microshutter is to allow as much light through as possible when the shutters are in the open configuration, and preventing any light from passing through when they are in the closed position. The ratio of the transmitted light that is blocked is defined here as the contrast ratio. Three major components contribute to the improved performance of these microshutters: 1. The precise implementation of light shields, which protect the gap around the shutters so no light can leak through. It has been ascertained that without the light shield there would be a gap on the order of 1 percent of the shutter area, limiting the contrast to a maximum of 100. 2. The precise coating of the interior wall of each microshutter was improved with an insulator and metal using an angle deposition technique. The coating prevents any infrared light that finds an entrance on the surface of the microshutter cell from being emitted from a sidewall. Since silicon is in effect transparent to any light with a wavelength longer than .1 micrometer, these coatings are essential to blocking any stray signals when the shutters are closed. 3. A thin film of molybdenum nitride (MoN) was integrated onto the surface of the microshutter blade. This film provides the majority of light blockage over the microshutter and also ensures that the shutter can be operated over a wide temperature range by maintaining its flatness. These improvements were motivated by the requirements dictated by the James Webb Space Telescope NIRSpec instrument. The science goals of the NIRSpec require observing some of the very faintest objects in a given field of view that also may contain some very bright objects. To observe the faint objects, the light from the bright objects - which could be thousands of times brighter - must be completely blocked. If a closed microshutter is even slightly transmissive, a very bright object will still transmit a small signal, which can be larger than a signal from a very faint object transmitted through an open shutter. Since this situation can completely corrupt the results, it was necessary that the closed shutters be able to attenuate light by at least a factor of 2,000. There currently exist four flight-quality microshutter arrays that have been fully or are currently undergoing testing and the results support that the three improvements described above have successfully led to contrast levels greater than 50,000 in over 99 percent of the microshutters at an operating temperature of 35 K. Applications for these high-contrast microshutters are in the photomask generation and stepper equipment used to make integrated circuits and microelectromechanical (MEMS) devices. Since microshutters are a reconfigurable optical element, their versatility in these industries provides an improvement over printed masks and fixed projection alignment systems.

Jhabvala, Murzy↗
DOE OSTI2022

φ meson production in ρ + Al, ρ + Au, d + Au , and 3 He + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

We report small nuclear collisions are mainly sensitive to cold-nuclear-matter effects; however, the collective behavior observed in these collisions shows a hint of hot-nuclear-matter effects. The identified-particle spectra, especially the φ mesons which contain strange and antistrange quarks and have a relatively small hadronic-interaction cross section, are a good tool to study these effects. The PHENIX experiment has measured φ mesons in a specific set of small collision systems p+Al, p+Au, and 3 He+Au, as well as d+Au at $\sqrt{s_{NN}}$=200 GeV. The transverse-momentum spectra and nuclear-modification factors are presented and compared to theoretical-model predictions. The comparisons with different calculations suggest that quark-gluon plasma may be formed in these small collision systems at $\sqrt{s_{NN}}$=200 GeV. However, the volume and the lifetime of the produced medium may be insufficient for observing strangeness-enhancement and jet-quenching effects. The comparison with calculations suggests that the main production mechanisms of φ mesons at midrapidity may be different in p+Al versus p/d/ 3 He+Au collisions at $\sqrt{s_{NN}}$=200 GeV. While thermal quark recombination seems to dominate in p/d/ 3 He+Au collisions, fragmentation seems to be the main production mechanism in p+Al collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Properties of Ground Level Enhancement Events and the Associated Solar Eruptions During Solar Cycle 23

Solar cycle 23 witnessed the most complete set of observations of coronal mass ejections (CMEs) associated with the Ground Level Enhancement (GLE) events. We present an overview of the observed properties of the GLEs and those of the two associated phenomena, viz., flares and CMEs, both being potential sources of particle acceleration. Although we do not find a striking correlation between the GLE intensity and the parameters of flares and CMEs, the solar eruptions are very intense involving X-class flares and extreme CME speeds (average approx. 2000 km/s). An M7.1 flare and a 1200 km/s CME are the weakest events in the list of 16 GLE events. Most (80 %) of the CMEs are full halos with the three non-halos having widths in the range 167 to 212 degrees. The active regions in which the GLE events originate are generally large: 1290 msh (median 1010 msh) compared to 934 msh (median: 790 msh) for SEP-producing active regions. For accurate estimation of the CME height at the time of metric type II onset and GLE particle release, we estimated the initial acceleration of the CMEs using flare and CME observations. The initial acceleration of GLE-associated CMEs is much larger (by a factor of 2) than that of ordinary CMEs (2.3 km/sq s vs. 1 km/sq s). We confirmed the initial acceleration for two events for which CME measurements are available in the inner corona. The GLE particle release is delayed with respect to the onset of all electromagnetic signatures of the eruptions: type II bursts, low frequency type III bursts, soft X-ray flares and CMEs. The presence of metric type II radio bursts some 17 min (median: 16 min; range: 3 to 48 min) before the GLE onset indicates shock formation well before the particle release. The release of GLE particles occurs when the CMEs reach an average height of approx 3.09 R(sub s) (median: 3.18 R (sub s) ; range: 1.71 to 4.01 R (sub s) ) for well-connected events (source longitude in the range W20–W90). For poorly connected events, the average CME height at GLE particle release is ∼66 % larger (mean: 5.18 R (sub s) ; median: 4.61 R (sub s) ; range: 2.75–8.49 R (sub s) ). The longitudinal dependence is consistent with shock accelerations because the shocks from poorly connected events need to expand more to cross the field lines connecting to an Earth observer. On the other hand, the CME height at metric type II burst onset has no longitudinal dependence because electromagnetic signals do not require magnetic connectivity to the observer. For several events, the GLE particle release is very close to the time of first appearance of the CME in the coronagraphic field of view, so we independently confirmed the CME height at particle release. The CME height at metric type II burst onset is in the narrow range 1.29 to 1.8 R(sub s), with mean and median values of 1.53 and 1.47 R(sub s). The CME heights at metric type II burst onset and GLE particle release correspond to the minimum and maximum in the Alfven speed profile. The increase in CME speed between these two heights suggests an increase in Alfvenic Mach number from 2 to 3. The CME heights at GLE particle release are in good agreement with those obtained from the velocity dispersion analysis, including the source longitude dependence. We also discuss the implications of the delay of GLE particle release with respect to complex type III bursts by approx 18 min (median: 16 in; range: 2 to 44 min) for the flare acceleration mechanism. A similar analysis is also performed on the delay of particle release relative to the hard X-ray emission.

Geomagnetic Storms↗
DOE OSTI2024

Electrode Elastic Modulus as the Dominant Factor in the Capping Effect in Ferroelectric Hafnium Zirconium Oxide Thin Films

The discovery of ferroelectricity in hafnia based thin films has catalyzed significant research focused on understanding the ferroelectric property origins and means to increase stability of the ferroelectric phase. Prior studies have revealed that biaxial tensile stress via an electrode “capping effect” is a suspected ferroelectric phase stabilization mechanism. This effect is commonly reported to stem from a coefficient of thermal expansion (CTE) incongruency between the hafnia and top electrode. Despite reported correlations between ferroelectric phase fraction and electrode CTE, the thick silicon substrate dominates the mechanics and CTE-related stresses, negating any dominant contribution from an electrode CTE mismatch toward the capping effect. In this work, these discrepancies are reconciled, and the origin of these differences deriving from electrode elastic modulus, not CTE, is demonstrated. Pt/M/TaN/Hf 0.5 Zr 0.5 O 2 /TaN/Si devices, where M is platinum, TaN, iridium, tungsten, and ruthenium, were fabricated. Sin 2 (ψ)-based X-ray diffraction measurements of biaxial stress in the HZO layer reveal a strong correlation between biaxial stress, remanent polarization, and electrode elastic modulus. Conversely, a low correlation exists between the electrode CTE, HZO biaxial stress, and remanent polarization. A higher elastic modulus enhances the resistance to electrode elastic deformation, which intensifies the capping effect during crystallization, and culminates in the tandem restriction of out-of-plane hafnia volume expansion and preferential orientation of the polar c-axis normal to the plane. These behaviors concomitantly increase the ferroelectric phase stability and polarization magnitude. This work provides electrode material selection guidelines toward the development of high-performing ferroelectric hafnia into microelectronic devices, such as nonvolatile memories.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Raman-shifting an ArF excimer laser to generate new lines for obtaining optical diagnostic based information in flow fields

The application of tunable excimer lasers in combustion and flow diagnostics is almost routine nowadays. The properties of this laser system that enable density and temperature measurements in supersonic and hypersonic flow fields to be conducted are its high power, high repetition rate, and high spectral brightness. The limitation imposed by this system on these measurements is the paucity of lines in the wavelength region, the vacuum-ultraviolet, where species of interest, such as OH, N2, O2, H2, H2O, CO, NO, etc., are susceptible to electronic excitation to high-lying states. To circumvent this problem one normally resorts to nonlinear optical techniques such as frequency conversion via stimulated Raman scattering (SRS), more commonly known as Raman shifting or Raman mixing, to extend these nonintrusive and nonperturbing techniques to the shorter wavelengths in the VUV region and, for that matter, to longer wavelengths in the infrared region, if the need arises. The theoretical basis of SRS and its application are well documented in the literature. In essence, the Raman shift is a consequence of the inelastic scattering of the incident radiation by the sample. Most of the scattered radiation from the molecules of the sample is unchanged in frequency. However, a small fraction of the incident radiation is changed in frequency. This shift is a result of the fact that some of the incident photons on colliding with the molecules of the sample give up some of their energy and emerge with a lower energy resulting in the lower-frequency Stokes radiation. Other incident photons may increase their energy by colliding with the vibrationally excited molecules of the medium and emerge as higher-frequency antistokes radiation. The generation of the latter is the main objective of this project. The process, however, depends on several factors, including the beam quality of the pump laser, the cross-section of the gaseous medium, the gas pressure, and the ambient temperature of the gas near the focal region. Furthermore, since the Raman shifting process is polarization sensitive, it is necessary to have all of the laser energy in a single polarization. These factors were taken into consideration in the execution of the project. The implementation of the Raman shift was accomplished by focusing the 193 nm output of an ArF excimer laser (Lamda-Physik LPX 150) into a 1-meter long high pressure recirculating Raman cell filled with H2 gas. The laser system was modified in order to improve the mode quality of the pump beam to enhance the Raman shifting. To accomplish this feat, a prism beam expander and grating on the oscillator discharge provided wavelength tuning over the excimer gain profile. Furthermore, a triple-pass configuration, as opposed to unstable resonator optics, was employed in the operation of the amplifier cavity so that when the oscillator output radiation, focused by a 51-cm focal length fused silica lens through a 50 micron pinhole (serving as a spatial filter) and recollimated with a 25-cm focal length lens, was fed into the amplifier, it was injection locked, thereby providing tunable radiation with relatively low divergence. The forward scattered radiation emanating from the impingement of the modified pump beam on the Raman cell was detected using an energy meter after the latter had been separated from it using a dispersing prism.

Koker, Edmond B.↗
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
NASA NTRSDate not supplied

Enhancing Air Traffic Control Planning with Automatic Speech Recognition

The decisions made during the Federal Aviation Administration Air Traffic Control System Command Center's planning teleconferences hold significant sway over the National Airspace System. Held every two hours, these teleconferences convene air traffic managers and stakeholders from across the nation to discuss airspace conditions, weather, and constraints, leading to the formulation and adjustment of traffic management initiatives. Given the critical nature of these decisions, the need for accurate and efficient record-keeping is paramount. In recent years, the application of automatic speech recognition has gained popularity across diverse industries, including aviation. While traditional applications focus on transcribing air traffic control communication, this paper explores a unique application of automatic speech recognition by converting the audio from planning teleconferences into text transcriptions. This innovative approach addresses key challenges in the field, presenting potential benefits for quality assurance, real-time participation, and downstream natural language processing tasks. A notable breakthrough in the machine learning community, namely the transformer neural network architecture, forms the backbone of the proposed solution in this paper. The transformer architecture's role in this research represents a paradigm shift in the efficiency of automatic speech recognition models. By reducing the amount of in-domain training data required, this architecture allows for the fine-tuning of such models like Whisper, originally pretrained on vast English speech datasets. The adaptability of the transformer architecture proves invaluable in capturing the nuances of aviation terminology and specific language used in planning teleconferences. Leveraging the Whisper model as a baseline, our research details the fine-tuning and validation using a dataset comprising 20 hours of meticulously transcribed planning teleconferences. Notably, the baseline pretrained Whisper model exhibited a word error rate of 18.77%. Through the fine-tuning process, the model achieved a substantial improvement, demonstrating an impressive performance with a reduced word error rate of 6.82%. This substantial decrease in WER not only highlights the effectiveness of the transformer architecture but also emphasizes the practical advancements achieved through the application of automatic speech recognition in this specific domain. The utilization of automatic speech recognition in planning teleconferences in this work introduces several novelties. Firstly, the creation of text transcriptions offers a valuable tool for quality assurance and facilitates the efficient review of teleconferences. This is an important aspect of the proposed solution, given the time-sensitive and high-stakes nature of decisions made during these meetings. Furthermore, text-searchable transcriptions provide a streamlined approach for locating and validating critical information, potentially saving hours of manual effort in searching through audio recordings. Moreover, our research identifies a key use case for external facilities and stakeholders. In situations where attendance at the planning teleconference is not feasible, having access to text transcriptions in real-time or shortly after the teleconference ends, proves to be a time-saving and informative resource. This feature enhances collaboration and ensures that stakeholders can stay abreast of important discussions and decisions even in their absence. Despite the efficiency gains facilitated by the transformer architecture in automatic speech recognition technology, it is essential to acknowledge the human factors in data creation. Subject matter experts play a crucial role in accurately transcribing planning teleconferences due to the specificity and complexity of the information discussed. The research dataset, consisting of 20 hours of transcribed planning teleconferences, forms the foundation for fine-tuning and validating the Whisper model. The achieved word error rate of 6.82% demonstrates promising advancements, particularly in recognizing essential aviation terminology within the teleconferences. In conclusion, this paper presents a comprehensive exploration of the application of automatic speech recognition in Air Traffic Control System Command Center planning teleconferences, leveraging the transformer architecture for enhanced efficiency. The novel contributions lie in the improved accessibility of decision-making records, real-time participation opportunities for external stakeholders, and the potential for downstream natural language processing advancements. As the aviation industry continues to evolve, the integration of automatic speech recognition technologies holds the promise of revolutionizing decision-making processes and contributing to the overall safety and efficiency of air traffic management.

ATM↗
DOE OSTI2026

Advancing Artificial Intelligence with Liquid Argon Neutrino Experiments (Technical Report)

The grant allowed two main contributions: 1) The development of a first successful demonstration of the employment of Optimal Transport in liquid argon time projection chamber neutrino detectors. Optimal Transport, used in other contexts and specifically with LHC calorimetric data, was adapted to address a key particle identification challenge in LArTPCs: the separation of pi0 backgrounds from single-electrons produced in charged-current electron neutrino interactions. The work, leveraging ML methods such as k-nearest-neighbor (kNN) and support-vector-machine (SVM), showed an increase in background rejection of a factor of two or more. Work is now ongoing to incorporate this development in physics analyses for LArTPC experiments and more broadly expand the use of OT in LArTPC detectors including DUNE. This work was done in collaboration with the phenomenology group led by Nathaniel Craig at UCSB. 2) The deployment of NuGraph2, a graph neural network developed for LArTPC reconstruction, in the MicroBooNE experiment. NuGraph2 uses novel graph-neural-network methods on the rather simple LArTPC inputs of reconstructed hits, greatly simplifying the workflow compared to the use of waveform or signal-deconvolved wire ROIs. The network performed particle classification and was shown to address many challenging problems in LArTPC imaging including track-shower separation and the identification of protons and charged pions from primary muons. Our group collaborated with Giuseppe Cerati (FNAL scientist) who is one of the core developers of NuGraph2 to integrate this tool in MicroBooNE’s analysis framework. This consisted in tow key contributions: a) Studying performance on real data, which came with several months of iterations because the MC-trained version of the network was found to show significant bias that our group investigated and addressed. b) Integrating the output hit labeling of NuGraph2 into the existing particle tracking and shower reconstruction code. As a result of this work led by our team NuGraph2 is now enabling a suite of new analyses which benefit from enhanced capabilities and thus broader physics reach. The grant supported primarily the salary of UCSB graduate student Chuyue “Michaelia” Fang as well as partial summer salary support for PI Caratelli. Some funds were used for travel by Michaelia to ML related schools and conferences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Techno-economic life cycle assessment of CO 2 -EOR operations towards net negative emissions at farnsworth field unit

Optimizations of CO 2 Water Alternating Gas(WAG)- systems with multi-objectives of incremental recovery and maximization of CO 2 storage are challenging. Here, the incorporation of a total Greenhouse gas (GHG) life cycle assessment is mostly ignored leading to inaccurate estimation of overall net carbon emissions of their operations. In this study, the effect of a total GHG life cycle assessment on a multi-objective CO 2 -WAG optimization with integrated techno-economic assessment (TEA) which factors carbon tax credit is conducted. A life cycle assessment (LCA) was conducted utilizing a 20 -year optimized post history matched data from a high fidelity reservoir simulation model. Using data generated from the optimum result, a techno-economic life cycle analysis was further conducted. The first scenario classified as the base model had an estimated 81% of purchased CO 2 sequestered. The results through a comprehensive techno-economic LCA model yielded a net estimate of 73% of purchased CO 2 . The optimized forecasted model which considered key operational and reservoir factors such as WAG ratio, injection rates and periods, and well specification resulted in an improved sequestration of 92% of purchased CO 2 . However, this also dropped to 84% after taking it through LCA. These results clearly indicate a significant amount of net CO 2 is not accounted for when operations are not analyzed through LCA. From the LCA, direct flaring volumes of CO 2 , energy consumption and efficiency of unit equipment were noticed to be the major causes of these reductions. Considering ten main sources of energy as source of energy generation, a comparative techno-eco LCA was conducted. The results confirmed a lower net volume and NPV for energy sources with higher carbon footprints and vice versa. Thus, a total LCA of CO 2 -WAG greatly influences net storage factor of purchased CO 2 and hence project NPV where tax credit/incentives per ton of CO 2 sequestered is considered. Although operational conditions are optimized for best results, there are significant factors that leads to minimization of net storage factor. This study therefore provides an insightful information for optimizing CO 2 -WAG multi-objectives to achieve minimum GHG emission.

02 PETROLEUM↗
DOE OSTI2024

Identified charged-hadron production in 𝑝+Al, 3 He + Au, and Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and in U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV

The PHENIX experiment has performed a systematic study of identified charged-hadron (𝜋 ± , 𝐾 ± , 𝑝, $\bar{𝑝}$) production at midrapidity in 𝑝 + Al, 3 He + Au, and Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Identified charged-hadron invariant transverse-momentum (𝑝 𝑇 ) and transverse-mass (𝑚 𝑇 ) spectra are presented and interpreted in terms of radially expanding thermalized systems. The particle ratios of 𝐾/𝜋 and 𝑝/𝜋 have been measured in different centrality ranges of large (Cu + Au and U + U) and small (𝑝 + Al and 3 He + Au) collision systems. The values of 𝐾/𝜋 ratios measured in all considered collision systems were found to be consistent with those measured in 𝑝+𝑝 collisions. However, the values of 𝑝/𝜋 ratios measured in large collision systems reach the values of ≈0.6, which is a factor of ≈2 larger than in 𝑝 + 𝑝 collisions. These results can be qualitatively understood in terms of the baryon enhancement expected from hadronization by recombination. Identified charged-hadron nuclear-modification factors (𝑅 𝐴⁢𝐵 ) are also presented. Enhancement of proton 𝑅 𝐴⁢𝐵 values over meson 𝑅 𝐴⁢𝐵 values was observed in central 3 He + Au, Cu + Au, and U + U collisions. Finally, the proton 𝑅 𝐴⁢𝐵 values measured in the 𝑝 + Al collision system were found to be consistent with 𝑅 𝐴⁢𝐵 values of 𝜙, 𝜋 ± , 𝐾 ± , and 𝜋 0 mesons, which may indicate that the size of the system produced in 𝑝 + Al collisions is too small for recombination to cause a noticeable increase in proton production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Jet radius dependence of dijet momentum balance and suppression in Pb + Pb collisions at 5.02 TeV with the ATLAS detector

This paper describes a measurement of the jet radius dependence of the dijet momentum balance between leading back-to-back jets in 1.72 nb -1 of Pb + Pb collisions collected in 2018 and 255 pb -1 of pp collisions collected in 2017 by the ATLAS detector at the LHC. Both datasets were collected at $\sqrt{s_{NN}}$ = 5.02 TeV. Jets are reconstructed using the anti-k t algorithm with jet radius parameters R = 0.2, 0.3, 0.4, 0.5, and 0.6. The dijet momentum balance distributions are constructed for leading jets with transverse momentum p T from 100 to 562 GeV for R = 0.2, 0.3, and 0.4 jets, and from 158 to 562 GeV for R = 0.5 and 0.6 jets. The absolutely normalized dijet momentum balance distributions are constructed to compare measurements of the dijet yields in Pb + Pb collisions directly to the dijet cross sections in pp collisions. For all jet radii considered here, there is a suppression of more balanced dijets in Pb + Pb collisions compared with pp collisions, while for more imbalanced dijets there is an enhancement. There is a jet radius dependence to the dijet yields, being stronger for more imbalanced dijets than for more balanced dijets. Additionally, jet pair nuclear modification factors are measured. The subleading jet yields are found to be more suppressed than leading jet yields in dijets. A jet radius dependence of the pair nuclear modification factors is observed, with the suppression decreasing with increasing jet radius. These measurements provide new constraints on jet quenching scenarios in the quark-gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2011

Hydrofocusing Bioreactor Produces Anti-Cancer Alkaloids

A methodology for growing three-dimensional plant tissue models in a hydrodynamic focusing bioreactor (HFB) has been developed. The methodology is expected to be widely applicable, both on Earth and in outer space, as a means of growing plant cells and aggregates thereof under controlled conditions for diverse purposes, including research on effects of gravitation and other environmental factors upon plant growth and utilization of plant tissue cultures to produce drugs in quantities greater and at costs lower than those of conventional methodologies. The HFB was described in Hydro focus - ing Bioreactor for Three-Dimensional Cell Culture (MSC-22358), NASA Tech Briefs, Vol. 27, No. 3 (March 2003), page 66. To recapitulate: The HFB offers a unique hydrofocusing capability that enables the creation of a low-shear liquid culture environment simultaneously with the herding of suspended cells and tissue assemblies and removal of unwanted air bubbles. The HFB includes a rotating cell-culture vessel with a centrally located sampling port and an internal rotating viscous spinner attached to a rotating base. The vessel and viscous spinner can be made to rotate at the same speed and direction or different speeds and directions to tailor the flow field and the associated hydrodynamic forces in the vessel in order to obtain low-shear suspension of cells and control of the locations of cells and air bubbles. For research and pharmaceutical-production applications, the HFB offers two major benefits: low shear stress, which promotes the assembly of cells into tissue-like three-dimensional constructs; and randomization of gravitational vectors relative to cells, which affects production of medicinal compounds. Presumably, apposition of plant cells in the absence of shear forces promotes cell-cell contacts, cell aggregation, and cell differentiation. Only gentle mixing is necessary for distributing nutrients and oxygen. It has been postulated that inasmuch as cells in the simulated microgravitation of an HFB do not need to maintain the same surface forces as in normal Earth gravitation, they can divert more energy sources to growth and differentiation and, perhaps, to biosynthesis of greater quantities of desired medicinal compounds. Because one can adjust the HFB to vary effective gravitation, one can also test the effects of intermediate levels of gravitation on biosynthesis of various products. The potential utility of this methodology for producing drugs was demonstrated in experiments in which sandalwood and Madagascar periwinkle cells were grown in an HFB. The conditions in the HFB were chosen to induce the cells to form into aggregate cultures that produced anti-cancer indole alkaloids in amounts greater than do comparable numbers of cells of the same species cultured according to previously known methodologies. The observations made in these experiments were interpreted as suggesting that the aggregation of the cells might be responsible for the enhancement of production of alkaloids.

Gonda, Steve R.↗
DOE OSTI2021

Exploring and enhancing the accuracy of interior-scaled Perdew–Zunger self-interaction correction

The Perdew–Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange–correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z σ , which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional’s energy density. Here, the resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.

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

CubeSub

This presentation introduces and discusses the development of the CubeSub submersible concept, an Autonomous Underwater Vehicle (AUV) designed around the CubeSat satellite form factor. The presented work is part of the author's MSc thesis in Aerospace Engineering at the Royal Institute of Technology, Stockholm, Sweden, and was performed during an internship at the Mission Design Division of the NASA Ames Research Center, Moffett Field, CA. Still in the early stages of its development, the CubeSub is to become a submersible test-bed for technology qualified for underwater and space environments. With the long-term goal of exploring the underwater environments in outer space, such as the alleged subsurface ocean of Jupiter's moon Europa, a number of technology and operational procedures must be developed and matured. To assist in this, the CubeSub platform is introduced as a tool to allow engineers and scientists to easily test qualified technology underwater. A CubeSat is a class of miniaturized satellite built to a standardized size. The base size is 1U (U for unit), corresponding to a 100 x 100 x 113.5 cu mm cube. A 1U CubeSat can in other words easily be held in one hand. Stacking units give larger satellite sizes such as the also commonly used 1.5U, 2U and 3U. The CubeSat standard is in itself already well established and hundreds of CubeSats have to date been launched into space. Compatible technology is readily available and the know-how exists in the space industry, all of which makes it a firm ground to stand on for the CubeSub. The rationale behind using the CubeSat form factor is to make use of this pre-existing foundation, making the CubeSub easy to develop, modular and readily available. It will thereby aid in the process of maturing the concept of a fully space qualified submersible headed for outer space. As a further clarification, the CubeSub is itself not meant for outer space, but to facilitate development of such a vessel. Along with its uses as a testbed, the CubeSub also holds the potential to become a useful tool for exploration and experimentation here on Earth. A highly standardized system utilizing well-known hardware can reduce the cost and required work load for researchers wishing to perform experiments and exploration. Users could design sensors and experiments to comply with the already well established CubeSat standard, which are then carried by the CubeSub to the region of interest. This in turn means that the end users can focus more on formulating the experiment itself and less about how to get it where they want it. The CubeSub is designed to be built up by modules, which can be assembled in different configurations to fulfill different needs. Each module will be powered individually and intermodular communication will be wireless, removing the need for wiring. The inside of the cylindrical hull will be flooded with ambient water to enhance the interaction between payloads and surrounding environment. The overall torpedo-like shape is similar to that of a conventional AUV, slender and smooth. This is to make for a low drag, reduce the risk of snagging on surrounding objects and make it possible to deploy through an ice sheet via a narrow borehole or navigate in tight areas. To keep costs low and further accelerate development, rapid prototyping is utilized wherever possible. Full-scale prototypes are being constructed through 3D-printing and using COTS (Commercial Off-The-Shelf) components. 3D-printing is used both for the largest hull components and the relatively small and delicate propellers. Arduino boards are used for control and internal communication

AUV↗
DOE OSTI2024

A semi–automatic analytical methodology for characterizing the energy consumption of MRI systems using load duration curves

Background and purpose: Magnetic resonance imaging (MRI) scanners are a major contributor to greenhouse gas emissions from the healthcare sector, and efforts to improve energy efficiency and reduce energy consumption rely on quantification of the characteristics of energy consumption. The purpose of this work was to develop a semi-automatic analytical methodology for the characterization of the energy consumption of MRI systems using only the load duration curve (LDC). LDCs are a fundamental tool used across various fields to analyze and understand the behavior of loads over time. Methods: An electric current transformer sensor and data logger were installed on two 3T MRI scanners from two vendors, termed M1 (outpatient scanner) and M2 (inpatient/emergency scanner). Data was collected for 1 month (7/11/2023 to 8/11/2023). Active power was calculated, assuming a balanced three-phase system, using the average current measured across all three phases, a 480 V reference voltage for both machines, and vendor-provided power factors. An LDC was constructed for each system by sorting the active power values in descending order and computing the cumulative time (in units of percentage) for each data point. The first derivative of the LDC was then computed (LDC’), smoothed by convolution with a window function (sLDC’), and used to detect transitions between different system modes including (in descending power levels): scan, prepared-to-scan, idle, low-power, and off. The final, segmented LDC was used to measure time (% total time), total energy (kWh), and mean power (kW) for each system mode on both scanners. The method was validated by comparing mean power values, computed using the segmented 1-month LDC, for each nonproductive system mode (i.e., prepared-to-scan, idle, lower-power, and off) against power levels measured after a deliberate system shutdown was performed for each scanner (1 day worth of data). Results: The validation revealed differences in mean power values <1.4% for all nonproductive modes and both scanners. In the scan system mode, the mean power values ranged from 29.8 to 37.2 kW and the total energy consumed for 1 month ranged from 11 106 to 14 466 kWh depending on the scanner. Over the course of 1 month, the portion of time the scanners were in nonproductive modes ranged from 76% to 80% across scanners and the nonproductive energy consumption ranged from 8010 to 6722 kWh depending on the scanner. The M1 (outpatient) scanner consumed 99.9 and 183.9 kWh/day in idle mode for weekdays and weekends, respectively, because the scanner spent 23% more time proportionally in idle mode on the weekends. Conclusions: A semi-automatic method for quantifying energy consumption characteristics of MRI scanners was introduced and validated. This method is relatively simple to implement as it requires only power data from the scanners and avoids the technical challenges associated with extracting and processing scanner log files. Finally, the methodology enables quantitative evaluation of the power, time, and energy characteristics of MRI scanners in scan and nonproductive system modes, providing baseline data and the capability of identifying potential opportunities for enhancing the energy efficiency of MRI scanners.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗