Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inference accelerators”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Indicator-directed Dynamic Power Management for Iterative Workloads on GPU-Accelerated Systems

Modern high-performance and warehouse computing centers show strong interest in minimizing system power consumption while satisfying customers’ quality of service (QoS). Dynamic voltage and frequency scaling (DVFS) is effective for achieving this goal. Nevertheless, automating the process online and making it transparent to users must address three major challenges: (1) Complexity — today’s hardware components (e.g., CPUs, GPUs, memory, network, etc.) can be configured in several or dozens of frequency/voltage states for satisfying divergent system demands. Given their combination and the emergence of heterogeneity, searching the optimal configuration in the design space online can be timing consuming. (2) QoS guarantee — user-defined objectives such as power constraint and performance target must be monitored, predicted and ensured at the best effort. (3) Adaptability — various known and unknown workloads run on systems. Workloads characteristics should be quickly determined and configurations dynamically adjusted in accord with workloads and QoS. In this work, we focus on applications exhibiting an interesting feature – iterative or periodic, which is common among conventional HPC and emerging machine learning workloads. We propose an online dynamic power-performance (ODPP) management framework to dynamically adjust GPU DVFS configurations to meet performance and power objectives and constraints, without any code annotation or intrusion. Particularly, ODPP extracts the performance and power indicators for applications from their resources utilization profiles in a short episode. It further automatically constructs an accurate model that infers from the indicators how the application's performance and power vary with GPU core and memory frequencies. Aided with the model, for both seen and unseen applications, ODPP can quickly determine the most appropriate DVFS configuration for their execution. We evaluate ODPP on an NVIDIA GPU using multiple exascale computing (ECP) and deep learning applications.

Zou, Pengfei↗
DOE OSTI2022

Overview of the LANL Multi-Probe Radiography Project [Slides]

Future radiographic facilities for the U.S. defense program will be required to provide more information as simulation codes improve in both physics’ fidelity and resolution. A possible approach is to use more types of probe beams in addition to, or instead of X-rays, generated by 20-MeV electron accelerators from a small number of directions. High power short-pulse laser systems can generate beams of protons, neutrons and electrons, as well as X-rays. The cost of these systems is falling rapidly. So, it can be imagined that deploying multiple short-pulse lasers along with other, more traditional probes, will become feasible. In this project, we are following three paths to determine if such an approach will succeed for cm-scale objects. The first is an experimental one to determine if the presence of multiple short-pulse probes cause interference with each other, especially while radiographing dynamic objects. The second leg of this project is to determine if having multiple types of probes really does give more information on composition. Finally, an overall assessment of the viability of this approach will be made. Examples from recent experiments at the Omega EP laser will be presented. The initial approach to evaluating radiographs with multiple probes using the Bayesian Inference Engine (BIE) also will be given.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Interplanetary Protons versus Interacting Protons in the 2017 September 10 Solar Eruptive Event

We analyze the relativistic proton emission from the Sun during the eruptive event on 2017 September 10, which caused a ground-level enhancement (GLE 72) registered by the worldwide network of neutron monitors. Using the neutron monitor data and interplanetary transport modeling both along and across interplanetary magnetic field (IMF) lines, we deduce parameters of the proton injection into the interplanetary medium. The inferred injection profile of the interplanetary protons is compared with the profile of the >100 MeV γ -ray emission observed by the Fermi Large Area Telescope, attributed to pion production from the interaction of >300 MeV protons at the Sun. GLE 72 started with a prompt component that arrived along the IMF lines. This was followed by a more prolonged enhancement caused by protons arriving at the Earth across the IMF lines from the southwest. The interplanetary proton event is modeled using two sources—one source at the root of the Earth-connected IMF line and another source situated near the solar western limb. The maximum phase of the second injection of interplanetary protons coincides with the maximum phase of the prolonged >100 MeV γ -ray emission that originated from a small area at the solar western limb, below the current sheet trailing the associated coronal mass ejection (CME). A possible common source of interacting protons and interplanetary protons is discussed in terms of proton acceleration at the CME bow shock versus coronal (re-)acceleration in the wake of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

An initial framework for the rapid qualification of long-term creep rupture strength via microstructural modeling

This report describes the development and testing of a new method for extrapolating short-term creep rupture test data to predict long-term rupture strength. The goal of this work is to reduce the time required to qualify new materials for nuclear service by reducing the lead time required for dedicated, long-term material testing to establish key long-term material properties. The new approach described here uses a physics-based model to predict the long-term creep rupture strength of 316H stainless steel using only short-term test data. The key idea is to use Bayesian inference to find the statistical distribution of the model parameters that best explain the short-term rupture data. Because the model is physics-based these parameters are all microstructural quantities that can be measured through detailed material characterization experiments. The Bayesian prior distributions provide a means for incorporating this characterization data into the final model to improve the accuracy of the long-term model predictions. However, where such data is not available the process still produces an accurate model based on an uniformed prior. Our hypothesis is that this approach more accurately extrapolates the short-term test data when compared to current, empirical methods. The report proves this hypothesis using actual long-term rupture data available for 316H, including tests with rupture times greater than 200,000 hours. The general approach developed here could be applied to other materials and other time-dependent material properties. Applying this new technique to develop long-term qualified material properties, potentially in conjunction with other accelerated qualification approaches like staggered qualification test programs, could greatly reduce the time required to qualify new materials for nuclear service.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Massive, Position-Resolving, High-Energy-Resolution Detector for Non-Accelerator Cosmic and Intensity Frontier Particle Physics

We proposed to take the first steps in the development of a detector that promises energy resolution of tens of eV FWHM combined with robust nuclear-recoil discrimination and high fidelity, mm-precision position reconstruction for applications in non-accelerator particle physics at the Cosmic and Intensity Frontiers. The detector would obtain these excellent resolutions by sensing athermal phonons produced by particle interactions in crystalline, dielectric targets using a sensitive, highly multiplexable superconducting phonon sensor, the kinetic inductance detector (KID). This detector would be applicable to: the search for low-mass particle dark matter candidates with masses below 5~GeV via direct detection of scattering of dark matter particles with terrestrial nuclei; detection of coherent elastic neutrino-nucleus scattering to test for new physics such as a non-standard value of the weak nuclear charge, non-standard neutrino interactions (perhaps driven by a neutrino magnetic moment), or the existence of sterile neutrinos; and, searches for neutrinoless double-beta decay. During the funding period, we demonstrated scaling up of the detector concept from a 22-mm by 22-mm by 1-mm, 1-g prototype with 0.9~keV FWHM energy resolution to a 75-mm diameter by 1-mm, 9-g prototype while improving the inferred energy resolution to 0.7~keV~FWHM. In the process, we solved many problems associated with scaling device fabrication to large wafers and vastly reduced the fraction of the surface occupied by inactive but energy-absorbing metal. We also demonstrated a new technique that substantially simplifies the process of characterizing a new detector. These results provide a good foundation for future work scaling up the design to 4-mm thickness substrates and improving the resolution to reach 0.035~keV~FWHM, yielding a detector with compelling potential for the above applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Neutrino Beam Monitoring

Accelerator facilities produce neutrino beams from meson decays in a decay volume. Experiments measure event rates that depend on flux, cross sections, and detector response, so the flux is predicted using hadron production and beamline modeling and constrained by beam instrumentation, since near detectors alone cannot separate flux from cross section. Proton, hadron, and muon monitors can track the parent particle distributions and beam conditions, providing the inputs needed for flux predictions in long-baseline experiments such as NOvA, T2K, and DUNE. This talk reviews how beam monitors are used in practice to understand neutrino flux. Proton beam monitors tell where the beam hits the target and how stable it is. Farther downstream, hadron and muon monitors sample particles produced in meson decays. Because those muons come from the same parents as the neutrinos, their profiles reveal focusing, alignment shifts, and other changes in the beam, and they are routinely used to detect problems and guide flux predictions. The muon information can be used more quantitatively; for example, to infer the parent meson phase space, and fast radiation-hard timing detectors can add sensitivity to the momentum dependence of the focusing. These developments show both how tightly beam measurements can constrain the flux and where the current limits still lie. These approaches complement monitored-beam concepts, in which the decay region is instrumented to detect charged leptons from meson decays and to measure the neutrino flux directly.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗
DOE OSTI2026

Neutrino Beam Monitoring

Accelerator facilities produce neutrino beams from meson decays in a decay volume. Experiments measure event rates that depend on flux, cross sections, and detector response, so the flux is predicted using hadron production and beamline modeling and constrained by beam instrumentation, since near detectors alone cannot separate flux from cross section. Proton, hadron, and muon monitors can track the parent particle distributions and beam conditions, providing the inputs needed for flux predictions in long-baseline experiments such as NOvA, T2K, and DUNE. This talk reviews how beam monitors are used in practice to understand neutrino flux. Proton beam monitors tell where the beam hits the target and how stable it is. Farther downstream, hadron and muon monitors sample particles produced in meson decays. Because those muons come from the same parents as the neutrinos, their profiles reveal focusing, alignment shifts, and other changes in the beam, and they are routinely used to detect problems and guide flux predictions. The muon information can be used more quantitatively; for example, to infer the parent meson phase space, and fast radiation-hard timing detectors can add sensitivity to the momentum dependence of the focusing. These developments show both how tightly beam measurements can constrain the flux and where the current limits still lie. These approaches complement monitored-beam concepts, in which the decay region is instrumented to detect charged leptons from meson decays and to measure the neutrino flux directly.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗
DOE OSTI2025

Ultra-Fast Non-Volatile Resistive Switching Devices with Over 512 Distinct and Stable Levels for Memory and Neuromorphic Computing

Low-current multilevel programmability with inherent non-volatility and high stability of resistance states is required for both multi-bit memory storage and deep learning accelerators but is difficult to achieve. Here, in a resistive switching system, this work realizes >512 (>9 bits) distinct non-volatile conductance levels with stable retention for each state with current levels down to the nanoampere range, highly promising for potential integration with small processing nodes with ultra-low power consumption requirements. This is achieved by demonstrating a new thin film design concept that encompasses three key features: an ultra-thin epitaxial oxygen ionic switching layer that provides a tunable energy barrier at the bottom electrode, an overcoat amorphous layer that acts as an ion migration barrier for stable state retention, and a partial conductive filament as a localized electronic transport channel to the epitaxial switching layer. A large dynamic resistance range of up to seven orders of magnitude is achieved with reset-free transitions among intermediate states, and programmability is demonstrated with ultra-fast (20 ns) pulses. Artificial neural network (ANN) simulations, based on the experimental performance and its non-idealities, demonstrate close-to-ideal inference accuracies for various Modified National Institute of Standards and Technology (MNIST) data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CSB-RNN: A Faster-Than-Realtime RNN Acceleration Framework with Compressed Structured Blocks

Recurrent Neural Networks (RNN) is widely applied to temporal sequence analysis, where real-time performance is usually in demand. However, RNN suffers a heavy computational workload as the model comes with a large weight matrix. To alleviate the pain, model compression (pruning) schemes have been proposed for RNN that pruning the redundant (near-zero) weight-values. On the one hand, the non-structured pruning methods achieve a considerable pruning rate while bringing the computational irregularity, which is un-friendly to parallel-hardware. On the other hand, the existing structured pruning methods consider the hardware parallelism; However, they suffer a poor pruning rate due to the restrict constraints on pruning structure. This paper presents CSB-RNN, an optimized full-stack RNN framework with the novel compressed structured block (CSB) technique. The CSB-pruned RNN model comes with both fine-granularity that benefits the pruning rate and regular structure that facilitates the hardware-parallelism. Further, we propose a novel hardware architecture for inferencing the CSB-pruned model. Different from conventional parallel hardware, this architecture solves the block-workload imbalance issue and achieves an over 95% hardware utilization. With the experiments on 10 RNN models in 5 application domains, the CSB-RNN realizes 7×-20× lossless compression and up to 50× acceptable lossy-compression, which is 2×-7× to the prior art. With the addition of the novel hardware, the compressed-RNN inference reaches a super real-time latency of 10-400µs with FPGA implementation.

Shi, Runbin↗
DOE OSTI2023

Semi-supervised graph neural networks for pileup noise removal

Abstract The high instantaneous luminosity of the CERN Large Hadron Collider leads to multiple proton–proton interactions in the same or nearby bunch crossings (pileup). Advanced pileup mitigation algorithms are designed to remove this noise from pileup particles and improve the performance of crucial physics observables. This study implements a semi-supervised graph neural network for particle-level pileup noise removal, by identifying individual particles produced from pileup. The graph neural network is firstly trained on charged particles with known labels, which can be obtained from detector measurements on data or simulation, and then inferred on neutral particles for which such labels are missing. This semi-supervised approach does not depend on the neutral particle pileup label information from simulation, and thus allows us to perform training directly on experimental data. The performance of this approach is found to be consistently better than widely-used domain algorithms and comparable to the fully-supervised training using simulation truth information. The study serves as the first attempt at applying semi-supervised learning techniques to pileup mitigation, and opens up a new direction of fully data-driven machine learning pileup mitigation studies.

43 PARTICLE ACCELERATORS↗
NASA NTRS2006

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗
DOE OSTI2024

Fast and Flexible Inference Framework for Continuum Reverberation Mapping Using Simulation-based Inference with Deep Learning

Continuum reverberation mapping (CRM) of active galactic nuclei (AGN) monitors multiwavelength variability signatures to constrain accretion disk structure and supermassive black hole (SMBH) properties. The upcoming Vera Rubin Observatory’s Legacy Survey of Space and Time will survey tens of millions of AGN over the next decade, with thousands of AGN monitored with almost daily cadence in the deep drilling fields. However, existing CRM methodologies often require long computation time and are not designed to handle such large amounts of data. In this paper, we present a fast and flexible inference framework for CRM using simulation-based inference (SBI) with deep learning to estimate SMBH properties from AGN light curves. We use a long short-term memory summary network to reduce the high dimensionality of the light curve data and then use a neural density estimator to estimate the posterior of SMBH parameters. Using simulated light curves, we find SBI can produce more accurate SMBH parameter estimation with 10 3 –10 5 times speed up in inference efficiency compared to traditional methods. The SBI framework is particularly suitable for wide-field CRM surveys as the light curves will have identical observing patterns, which can be incorporated into the SBI simulation. We explore the performance of our SBI model on light curves with irregular-sampled, realistic observing cadence and alternative variability characteristics to demonstrate the flexibility and limitation of the SBI framework.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Statistical survey of pitch angle distributions in core (0-50 eV) ions from Dynamics Explorer 1: Outflow in the auroral zone, polar cap, and cusp

Core (0-50 eV) ion pitch angle measurements from the retarding ion mass spectrometer on Dynamics Explorer 1 are examined with respect to magnetic disturbance, invariant latitude, magnetic local time, and altitude for ions H(+), He(+), O(+), M/Z = 2 (D(+) or He(++)), and O(++). Included are outflow events in the auroral zone, polar cap, and cusp, separated into altitude regions below and above 3 R(sub E). In addition to the customary division into beam, conic, and upwelling distributions, the high-latitude observations fall into three categories corresponding to ion bulk speeds that are (1) less than, (2) comparable to, or (3) faster than that of the spacecraft. This separation, along with the altitude partition, serves to identify conditions under which ionospheric source ions are gravita- tionally bound and when they are more energetic and able to escape to the outer magnetosphere. Features of the cleft ion fountain inferred from single event studies are clearly identifiable in the statistical results. In addition, it is found that the dayside pre-noon cleft is a dayside afternoon cleft, or auroral zone, becomes an additional source for increased activity. The auroral oval as a whole appears to be a steady source of escape velocity H(+), a steady source of escape velocity He(+) ions for the dusk sector, and a source of escape velocity heavy ions for dusk local times primarily during increased activity. The polar cap above the auroral zone is a consistent source of low-energy ions, although only the lighter mass particles appear to have sufficient velocity, on average, to escape to higher altitudes. The observations support two concepts for outflow: (1) The cleft ion fountain consists of ionospheric plasma of 1-20 eV energy streaming upward into the magnetosphere where high-latitude convection electric fields cause poleward dispersion. (2) The auroral ion fountain involves field-aligned beams which flow out along auroral latitude field lines; and, in addition, for late afternoon local times, they experience additional acceleration such that the ion energy distribution tends to exceed the detection range of the instrument (greater than 50-60 eV).

Giles, B. L.↗
DOE OSTI2026

RTN-124: Photometric Redshifts for the Vera C. Rubin Observatory Data Preview

We present the photometric redshifts (photo-z) inferred using algorithms implemented in the Redshift Assessment Infrastructure Layers (RAIL) for the NSF-DOE Vera C. Rubin Observatory Data Preview 2 (DP2). We produce a compilation of reference redshift catalog using spectroscopic, grism and many band photometric redshift dataset hosted on the LIneA Photo-z Server. We curate training and testing set for assessing the scientific and technical performance of Rubin photo-z. The algorithm applied to the object catalog are FlexZBoost, BPZ, kNN, GPz, DNF and TPz; with a combination of 6-band and 4-band photo-z depending on availability of u and y photometry. The redshift point estimates and uncertainty estimation in tabular format through the Large Survey DataBase (LSDB).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗
DOE OSTI2025

Inference of response functions with the help of machine-learning algorithms

Response functions are a key quantity to describe the near-equilibrium dynamics of strongly interacting many-body systems. Recent techniques that attempt to overcome the challenges of calculating these ab initio have employed expansions in terms of orthogonal polynomials. We employ a neural network prediction algorithm to reconstruct a response function 𝑆⁡(𝜔) defined over a range in frequencies 𝜔. Here, we represent the calculated response function as a truncated Chebyshev series whose coefficients can be optimized to reduce the representation error. We compare the quality of response functions obtained using coefficients calculated using a neural network (NN) algorithm with those computed using the Gaussian integral transform (GIT) method. In the regime where only a small number of terms in the Chebyshev series are retained, we find that the NN scheme outperforms the GIT method.

Kurkcuoglu, Doga Murat [Fermi National Accelerator↗
DOE OSTI2024

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗