Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Injection testing”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Water Coolant Pool

Mechanistic assessments of radionuclide release during postulated accidents are expected to be included in advanced reactor license applications. The mechanistic source term (MST) provides an opportunity for vendors to realistically evaluate the radiological consequences of an incident, and may aid in justifying reduced emergency planning zones and plant sites. However, the development of MSTs for advanced nuclear reactors is challenging because there are numerous phenomena that can affect the transport and retention of radionuclides. As part of a trial MST assessment for a metal-fueled, pool-type sodium cooled fast reactor (SFR), led by Argonne National Laboratory, a simplified radionuclide transport code (SRT code) was developed, which includes models to estimate the quantity of fission product aerosols scrubbed in the sodium pool during postulated accident scenarios. In a pool-type SFR, when fission products are released into the coolant pool due to failure of fuel pins, most of the radionuclides are scrubbed by the coolant pool, but some have the potential to migrate to the cover gas region through entrainment within gas bubbles. The SRT code contains a model that evaluates this scrubbing behavior and calculates the fraction of fission product aerosols that reach the cover gas. Due to a lack of available validation data for sodium pool scrubbing, the U.S. Department of Energy funded an experiment at the University of Wisconsin-Madison to measure aerosol scrubbing by injecting air bubbles containing aerosol into a coolant pool. Prior to performing an experiment with liquid sodium, a water loop experiment was performed. Their experiment evaluated the effect of changing the aerosol size, aerosol density, aerosol concentration, bubble size, and pool depth on the aerosol scrubbing efficiency of the pool. In this benchmark experiment, the base tests were conducted by repeated tests of isolated bubbles. Afterwards, more prototypic tests with bubble swarms were performed to evaluate the interactions between the bubbles. The bubble swarm test was able to confirm that a larger amount of aerosol scrubbing occurred than the single bubble test. It was also confirmed that as the bubble size, aerosol density, and pool height increase, the extent of pool scrubbing also increases and does not change with the aerosol concentration. In addition, since the degree of scrubbing is the lowest at aerosol sizes between 0.01 and 1 μm, that is, the largest amount of aerosol is emitted, it was confirmed that the analysis of this size in MST is the most important. This benchmark experiment informs the direction of future sodium experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Alfven eigenmodes and fast ion transport in negative triangularity DIII-D plasmas

The first energetic particle experiments in negative triangularity tokamak plasmas have been carried out on DIII-D. Alfvén eigenmode (AE) activity and associated fast ion transport comparable to that in positive triangularity is observed during the current ramp portion of all plasmas with early beam heating indicating negative triangularity does not confer a special advantage with respect to AE induced transport. In these discharges, a range of mode activity driven by the sub-Alfvénic (Vbeam/VA < 0.5) 80 kV neutral beams is found including beta induced Alfvén acoustic eigenmodes, beta induced Alfvén eigenmodes, reversed shear Alfvén eigenmodes (RSAEs) and toroidicity induced Alfvén eigenmodes (TAEs). Mode intermittency and possibly chirping appears to be more common than comparable positive triangularity and/or oval discharges but overall, the unstable spectra and mode amplitudes observed on magnetics, CO2 interferometry, and electron cyclotron emission in a set of matched positive and negative triangularity cases at moderate beam power is similar. Large levels of Alfvén eigenmode induced fast ion transport are found in both positive and negative triangularity with up to 70% central fast ion pressure deficits relative to classical predictions early during the current ramp phase for discharges with 3 MW injected 80 kV neutral beam power. The deficit in both cases is reduced toward zero as the current penetrates and eventually reaches classical levels during current flattop with qmin ≈ 1. Fast ion transport in negative triangularity plasmas measured using a beam modulation technique show very similar levels to those measured in oval plasmas over the range of tested beam powers (Pbeam ≈ 2–7 MW). This similarity is found despite quite different unstable spectra in the highest beam power case where a mixture of coherent and quasi-coherent modes in the TAE frequency range are observed in negative triangularity and a spectrum of more typical narrowband TAEs and RSAEs is observed for the oval plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gain Stability of Hamamatsu R5912-MOD Photo-Multipliers at Low Temperature

The Photon Detection System is a crucial component of the ICARUS detector in the Short Baseline Neutrino (SBN) Program at Fermi National Accelerator Laboratory (FNAL). It consists of 360 Photo Multiplier Tubes (PMTs) 8” Hamamatsu R5912-MOD and, since June 2020, it has been operated at the liquid Argon cryogenic temperature. During these years, the PMTs have shown gain losses as an aging effect. Using a climatic chamber at the INFN Sezione di Catania, we confirmed a stable gain for a single PMT 8” Hamamatsu 5912-MOD when operated at room temperature, and a persistent gain loss at temperatures around -70°C, although far from the liquid Argon one. No gain recovery was obtained bringing the PMT back to room temperature. We suspect that dynode deterioration due to temperature gradient plays a significant role in this phenomenon. During our test, the PMT was illuminated by a 520 nm laser source and operated in current mode. The laser beam was split into two optical fibers, the first one transmitting light through Neutral Density filters and then diffusing it over the surface of the PMT photocathode. The other optical fiber was used for an independent reference measurement of the injected light via a bolometric photodetector. In this presentation, the main technical characteristics of the measurement system are shown together with preliminary results which need further investigation to elucidate the underlying mechanisms driving the gain loss and to propose new mitigation strategies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Scale-Up and Modeling Efforts Using an Omniphobic Surface Treatment for Mitigating Solids Deposition

Abstract Gas hydrates, waxes and asphaltenes represent some of the most significant flow assurance challenges, especially in subsea lines, where treatment options can be limited. Currently, complete avoidance is the primary strategy for hydrate management, while chemical or mechanical flushing/pigging may be utilized for other solids. Each system must also be treated with individualized solutions, as there have been no proved one-size-fits-all technologies demonstrated to date. As an alternative to constant chemical injection or thermodynamic controls such as insulation or heating, a robust omniphobic surface treatment material has been developed which has been shown in previous studies to significantly reduce the adhesion of flow assurance solids, resulting in lower risk for deposition and plugging of gas hydrates, waxes, and asphaltenes. As part of a Department of Energy study, laboratory scale tests were performed on a variety of apparatuses ranging from micromechanical force (micron scale) to laboratory flow loop testing (meter scale). While the results from these tests have been promising, there can be some disconnects between lab-scale observations and field-scale testing. Several factors can result in discrepancies between small-scale bench-top studies and deployment-scale efforts. These include flow conditions (Reynolds number, flow regime), chemistry (synthetic vs. crude oils, gas composition, and additives), physical differences (pipeline material, size, and geometry), etc. This work will describe the actions taken to better understand the deployment and performance of a material such as this omniphobic surface treatment in field scenarios, and to build confidence for deployment in a production scenario. These efforts include modeling field production scenarios using models developed from lab observations, application of the material to field-scale equipment to better understand application challenges, and expanded survivability and longevity testing. Understanding the bridge between laboratory testing and full-scale deployment allows for better technique and risk mitigation for field scale testing.

03 NATURAL GAS↗
DOE OSTI2022

ThunderSecure: deploying real-time intrusion detection for 100G research networks by leveraging stream-based features and one-class classification network

Nowadays, data generated by large-scale scientific experiments are on the scale of petabytes per month. These data are transferred through dedicated high-bandwidth networks (40/100G) across distributed sites for processing, storage, and analysis. Like general purpose networks, research networks experience intrusions. However, monitoring anomalies in such high-speed network traffics is challenging given current cyber-infrastructure. Moreover, traditional network intrusion detection systems (NIDS) are signature based. However, anomaly patterns are difficult to define and that rulesets are often not updated frequently enough to reflect the changes of attack behaviors. We present ThunderSecure, a high-throughput, unsupervised learning-based intrusions detection system for 100G research networks. ThunderSecure implements an efficient packet processing and detection pipeline using multi-cores and GPUs. It extracts statistical and temporal features from real-time network data streams and feeds them to a one-class anomaly detection network. A baseline of normal distribution will be created based on the training observation. Testing traffic deviated from the learned profile will be marked as anomalies. We trained ThunderSecure on hundreds of billions of science data packets mirrored from two 100G network connections at Fermi National Accelerator Laboratory. The detection performance was evaluated on traffic captured from the same research network days and weeks after the training with different types of attack flows injected. Results show that ThunderSecure can recognize science data traffic captured long after the training and made nearly certain detection on the segment of the streams where anomalous flows were injected.

100G research network↗
DOE OSTI2021

Model predictive combustion control of a Gasoline Compression Ignition engine

Gasoline Compression Ignition is a novel combustion concept that derives its superiority from the high compression ratio of a compression ignition engine as well as the properties of gasoline fuel, such as longer ignition delay and higher volatility compared to diesel. Here, this combustion concept was experimentally tested on a 12.4L Class 8 truck engine. Based on these experimental data, prior efforts by the authors focused on the development of an engine model for a heavy-duty engine operating on a low-reactivity fuel. This engine model was leveraged within this study to investigate a combustion control strategy at different engine conditions and injection methods and was augmented to incorporate cycle-to-cycle combustion variations. State estimation is performed by means of a Kalman filter which feeds into a model predictive controller. The model predictive controller chooses control actions based on a predefined cost function under consideration of bounds reflecting physical constraints. The engine model was utilized to establish a state-space model that serves the Kalman filter and model predictive controller for estimation and prediction. A comparative study investigating control actions and engine behavior was performed with and without limiting in-cylinder peak pressure as well as combustion noise, which is of particular interest for early pilot injection strategies. In addition, the proposed control architecture was investigated at two different levels of cycle-to-cycle variations and compared to the performance of a control structure with input disturbance rejection. For increased cycle-to-cycle variations, disturbance estimation reduces state fluctuations and control effort. In general, this investigation highlights control aspects specific to a compression-ignited combustion regime with low-reactivity fuel. The control algorithm is able to maintain the desired references for brake mean effective pressure and combustion phasing while controlling peak in-cylinder pressure and combustion noise.

42 ENGINEERING↗
DOE OSTI2026

Searching for synchrotron emission from the geminga TeV halo using the planck satellite

Pulsars convert a significant fraction of their total spin-down power into very high-energy electrons, leading to the formation of TeV halos. While these halos are well characterized at TeV energies, it remains unclear whether pulsars also accelerate electrons efficiently at lower energies and how these particles propagate through their surrounding environments. We aim to test whether synchrotron emission from $\sim 50$–$300 \, \textrm{GeV}$ electrons around the Geminga pulsar can be detected in the frequency range observed by the Planck satellite. This would help constrain low-energy particle acceleration and diffusion in the vicinity of pulsars. We model the expected synchrotron emission from Geminga’s TeV halo based on various diffusion and injection spectrum scenarios and compare these predictions to publicly available multi-frequency Planck data. We find no conclusive evidence of spatially extended synchrotron emission associated with Geminga in any of Planck’s frequency bands. Our calculations show that even under favorable diffusion and injection conditions, the predicted synchrotron flux lies well below Planck’s measured background levels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗
DOE OSTI2023

Performance, Combustion and Emissions Evaluation of Liquid Phase Port-Injected LPG on a Single Cylinder Heavy-Duty Spark Ignited Engine

Liquefied petroleum gas (LPG), like many other alternative fuels, has witnessed increased adoption in the last decade, and its use is projected to rise as stricter emissions regulations continue to be applied. However, much of its use is limited to dual fuel applications, gaseous phase injection, light-duty passenger vehicle applications, or scenarios that require conversion from gasoline engines. Therefore, to address these limitations and discover the most efficient means of harnessing its full potential, more research is required in the development of optimized fuel injection equipment for liquid port and direct injection, along with the implementation of advanced combustion strategies that will improve its thermal efficiency to the levels of conventional fuels. This paper focuses on the development of a liquid phase port-injection system for LPG, the design of a reference piston, and the baseline evaluation of the performance, combustion, and emissions characteristics of a single cylinder research engine to establish a benchmark comparable to existing LPG engines. A sweep of start of injection (SOI) timing is performed by injecting liquid LPG at several closed and open intake valve timings, which demonstrates no significant variation in engine performance, but accounts for a 10% reduction in bsCO with the optimal SOI timing. Spark timing sweep demonstrates the 50% burn crank angle location related to maximum brake torque (MBT) point with a brake thermal efficiency (BTE) of ~34% for the tested load case. The effect of equivalence ratio is also presented with optimal SOI timing at MBT condition. The engine starts exhibiting knocking combustion at 140kPa intake manifold air pressure (IMAP) with a peak torque of 253Nm and a 5% reduction in brake specific fuel consumption compared to the naturally aspirated scenario.

Fosudo, Toluwalase↗
DOE OSTI2020

Adaptive Neuro Fuzzy Inference System for Cyber-Intrusion Detection in a Smart Grid

The evolution of the power grid has brought increasing deployment of advance metering infrastructure, penetration of intelligent electronic devices, and integration of physical power system components with information and communications technologies. With the fast-expanding connectivity, cyber vulnerabilities arise due to the use of internet-based communication systems. These systems are targets of cyber-intrusions which attempt to disturb the normal power system functions. Traditional intrusion detection algorithms have been developed without an explicit model of the cyber components. In this paper, an algorithm to detect false data injections in the power system is proposed considering both cyber and physical models of the power system. The algorithm is based on an Adaptive Neuro Fuzzy Inference System (ANFIS) which collects information from state variables of the cyber-physical system to meet the performance requirements of the grid. Simulations of the proposed approach using the IEEE 13-bus test system validate the effectiveness of this artificial intelligence-based algorithm.

Adaptive Neuro Fuzzy Inference System (ANFIS), cyb↗
DOE OSTI2022

Voltage Stability Constrained Moving Target Defense Against Net Load Redistribution Attacks

Moving target defense (MTD) using distributed flexible AC transmission system (D-FACTS) devices is a promising defense strategy to detect stealthy false data injection (FDI) attacks against the power system state estimation. However, all existing studies myopically perturb the reactance of D-FACTS lines without considering the system voltage stability. In this paper, we first illustrate voltage instability induced by MTDs in a three-bus system. To address this issue, we further propose a novel MTD framework that explicitly considers system voltage stability by using continuation power flow and voltage stability indices. We mathematically derive the sensitivity matrix of voltage stability index to line impedance, on which an optimization problem for maximizing voltage stability index is formulated. This framework is tested on the IEEE 14-bus and the IEEE 118-bus transmission systems, in which net load redistribution attacks are launched by sophisticated attackers. Here, the simulation results show the effectiveness of the proposed framework in circumventing the voltage instability while maintaining the detection effectiveness of MTD. We conduct case studies with and without the proposed framework under different MTD planning and operational methods. The impacts of the proposed two methods on attack detection effectiveness and system economic metrics are also revealed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗
DOE OSTI2024

In Situ Plasma processing of SRF cuperconducting cavities at JLAB, 2024 Update

Jefferson Lab has an ongoing R&D program in plasma processing. The experimental program investi-gated processing using argon/oxygen and heli-um/oxygen gas mixtures. Plasma processing is a com-mon technique where the free oxygen produced by the plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the secondary emission coefficient. The initial focus of the effort was processing C100 cavities by injecting RF power into the high order mode (HOM) coupler ports. We also developed the methods for establishing a plasma in C75 cryomodules where the RF power is injected via the fundamental power-coupler. Four C100 cryomodules were in situ processed in the CE-BAF accelerator in May 2023 with the cryomodules returning to an operational status in Sept. 2023. The overall operational energy gain for the four cryomod-ules was 49 MeV. Methods, systems and results from processing cryomodules in the CEBAF accelerator and vertical test results are presented. Current status and future plans are discussed.

Powers, T.↗
DOE OSTI2026

Gigacycle Fatigue Strength Evaluation of Welded 316L Stainless Steels for Mercury Target Vessel

At the Materials and Life Science Experimental Facility (MLF) in J-PARC, liquid mercury target for the pulsed spallation neutron source is in operation. An enclosure vessel for the liquid mercury target made of type 316L stainless steel (SS316L) suffers two kinds of cyclic stress during operation. One is the thermal stress due to the internal heating and swings by proton beam trip. The other is the impulsive stress by the pressure waves generated by the proton beam injection. The total number of loading cycles for the former is ∼104, and the latter is ∼4 × 108 for a year operation in the J-PARC mercury target vessel. The target vessel is assembled by an electron beam welding (EBW) and a gas tungsten arc welding (GTAW). However, fatigue data of welded SS316L up to gigacycle is limited. Ultrasonic fatigue testing, applying load cycles by utilizing ultrasonic resonance, for the welded SS316L was performed to investigate the effect of welding on fatigue behavior up to a gigacycle. The result showed that the fatigue strength degradation by EBW and EBW with GTAW were not recognized up to 109 cycles. Crack initiation in welded specimens nucliated at off-center areas of the specimen whereas the cracks in base metal specimen originated at the specimen center.

Naoe, Takashi [Japan Atomic Energy Agency (JAEA)]↗
DOE OSTI2026

High‑performance cold‑climate heat pump using tandem compressors with and without vapor injection: Laboratory investigation and field demonstration

Cold‑climate air‑source heat pumps must maintain heating capacity and stable operation at low ambient temperatures for residential applications. This study adapts tandem single-speed compressors—proven in commercial systems —for residential CCHPs, integrating VI and validating performance across laboratory and extreme field conditions to achieve cold-climate targets. That is to achieve ≥75% of the rated heating capacity at −25.0 °C relative to 8.3 °C and heating coefficient of performance (COP) greater than 4.0 at 8.3 °C. Laboratory testing confirmed that both configurations met these targets while operating within acceptable discharge‑temperature limits. Relative to the non‑VI configuration, the VI system provided up to 15% higher heating capacity and up to 9% higher heating COP under identical conditions, yielding a 5.6–9.1% increase in heating seasonal performance factor (HSPF). Field tests in Ohio and Alaska, with minimum outdoor temperatures of −25.0 °C and −34.0 °C, respectively, showed stable operation, minimal supplemental heating (<3.2%), low defrost penalties, and consistent output. Tandem‑compressor architectures proves promising for cold-climate residential heating.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗
DOE OSTI2023

Self-consistent investigation of density fueling needs on ITER and CFETR utilizing the new Pellet Ablation Module

Abstract Self-consistent modeling using the stability, transport, equilibrium, and pedestal (STEP) workflow in the OMFIT integrated modeling framework (predicting pedestal with EPED, core profiles with TGYRO, current profile with ONETWO, and EFIT for equilibrium) suggests ITER and future devices such as China Fusion Engineering Test Reactor (CFETR) Zhuang et al (2019 Nucl. Fusion 59 112010) will benefit from high-density operation (Greenwald limit fraction f g w ≈ 0.7−1.3). Regimes with an operational density near the Greenwald limit will likely need peaked density profiles so that the pedestal density remains below the Greenwald limit. Peaked density profiles can be achieved with the help of pellet injection. A flexible Pellet Ablation Module (PAM), which predicts the density source based on a comprehensive analytical pellet ablation model, has been developed for predicting pellet fueling for transport studies, and has been incorporated into the STEP workflow for predictive modeling. This workflow is applied to DIII-D and finds good agreement with experiments. On ITER the effect of pellet fueling is examined in an advanced inductive scenario, where a fusion gain of up to Q = 9 is predicted with strong central pellet fueling. On CFETR, with a mid-radius density source, an average of 1.5 × 10 22 electrons s −1 are required to achieve the density and temperature profiles necessary for the 1000 MW advanced scenario with a tritium burn-up fraction of ∼ 3 % .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗