Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plume flow”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Orphaned oil and gas well methane emission rates quantified using Gaussian plume inversions of ambient observations

Abstract. Annually, ∼ 3.6 million abandoned oil and gas wells in the US emit a combined ∼ 2.6 Tg methane (CH4), adversely affecting climate and regional air quality. However, these estimates depend on emission factors derived from measuring subpopulations of wells that vary by orders of magnitude due to very limited field sampling and poorly characterized distributions. Currently, US protocols to remediate orphaned wells lacks standardized quantification methods needed to both prioritize plugging and account for emission reductions. Therefore, sensitive, reliable, affordable, and scalable CH4 flux quantification methods are needed. We report the use of a simple Gaussian plume method where the dispersion parameters are constrained by in situ ground measurements of CH4 concentration at four locations 7.5–49 m downwind of the orphan well as well as local winds to estimate the leak rate from an orphan well. We derive a flux of 10.53 ± 1.16 kg CH4 h−1 during a venting procedure in April 2023 that agrees with the directly measured volumetric flow rate of 9.00 ± 0.25 kg CH4 h−1. This is 71 % greater than the 5.3 kg CH4 h−1 flux measured 7 months prior. Additionally, we discovered a secondary leak through the surface casing inferred as 0.43–0.67 kg CH4 h−1 both by our ground Gaussian analysis and by transecting the plume with an uncrewed aerial system (UAS). We show that in situ determination of the dispersion parameters used in our Gaussian inversions allows us to measure methane emissions to 15 % accuracy, significantly reducing errors when compared to the standard practice of assuming stability class. Our results help develop simpler methods and protocols for robust orphan well emission quantification that can be used for reporting.

Follansbee, Emily↗
DOE OSTI2020

River Plume Liftoff Dynamics and Surface Expressions

The water surface expression of liftoff and its dependence on discharge are examined using numerical simulations with the Regional Ocean Modeling System (ROMS). Liftoff is the process by which buoyant river water separates from the bed and flows over denser saltwater. During low-discharge conditions liftoff occurs in the river and is accompanied by a change in the surface slope. During high-discharge conditions liftoff occurs outside the mouth and generates a ridge on the water surface. The location and height of the ridge can be described by analytical equations in terms of discharge, shelf slope, and river mouth aspect ratio. The offshore distance and height of the ridge are proportional to the river discharge and vary inversely with river mouth aspect ratio. For steep shelf slopes liftoff occurs close to the river mouth and generates a large ridge. The ridge is modified, but not eliminated, by the presence of tides. The water surface slope change at the ridge peak is large enough to be detected by the upcoming Surface Water and Ocean Topography (SWOT) altimeter and can be used to identify the liftoff location during high discharge. However, during low discharge the water surface slope change at the liftoff location is too small to be detected by SWOT. These results indicate that remote measurements of the presence or absence of the ridge may be useful to distinguish between low and high flows, and remote measurements of the ridge location or height could be used to estimate freshwater discharge.

42 ENGINEERING↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electrical monitoring of saline tracers to reveal subsurface flow pathways in a flat ditch-drained field

Time-lapse electrical resistivity imaging (ERI) combined with salt tracers can enhance our understanding of subsurface flow and transport in ditch-drained agroecosystems. In a nearly level field adjacent to a drainage ditch, we installed 192 electrodes in a 72-m2 plot along with five hydrometric stations to enable time-lapse ERI monitoring during storms. Immediately prior to a wintertime rainstorm, we applied 5 kg of KBr to a shallow trench located roughly 11 m upgradient from the ditch. Antecedent conditions prior to the storm were wet, assuring a hydrologically connected state that favored tracer dissolution and movement in the subsurface. We monitored the evolution of the tracer plume with ERI over a seven-day period. Four-dimensional (4D) delineation of precipitation and tracer transport using time-lapse ERI showed that the tracer plume moved preferentially downward through a discontinuity in the argillic horizon within the first 57 hrs. Over the following 79 hrs, the tracer plume exhibited steady, lateral movement towards the drainage ditch within a more permeable sandy horizon situated just below the argillic. A spatial moment analysis of logarithmic changes in conductivity in the sand layer revealed seepage velocities ranging from 0.0027 to 0.0482 m d-1; these values were positively correlated with hydraulic gradients, suggesting laminar flow processes typical of sandy aquifers. Increasing groundwater velocities with hydraulic gradients suggest that solute transfers from groundwater to ditch are likely to be faster in storms, especially during wet periods. Even so, predominantly slow lateral transport within the sandy horizon, even during storms, indicates that critical source areas of pollutant loss are likely within a few meters of the ditch.

Robinson, Judith L.↗
DOE OSTI2022

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗
DOE OSTI2020

Temporal flow variations interact with spatial physical heterogeneity to impact solute transport in managed river corridors

The interactions between surface water and groundwater in river corridors lead to temporal fluctuations in subsurface water fluxes which have a critical role on solute transport dynamics. Here, we develop a framework to analyze the relative impacts of different temporal frequencies of the flow field in a spatially heterogeneous aquifer on solute transport. Our analysis indicates that the advection-dispersion equation behaves as a low-pass filter by wiping out the effect of high-frequency velocity fluctuations on the first two spatial moments of the solute plume, namely its center of mass and spreading. The concepts discussed in the theoretical analysis are then applied to understand solute transport dynamics at the 300 Area of the Hanford site (USA) adjacent to the Columbia River. We examine the temporal behavior of the solute plume's spatial moments for different temporal frequencies utilizing geostatistical parameters estimated in the 300 Area. Due to the proximity to the Columbia river, groundwater fluxes at the Hanford site are highly dynamic resulting in a large range of characteristic temporal frequencies. Nonetheless, similar to the theoretical analysis, our results show that the effect of high-frequency fluctuations is filtered, with most of the solute transport dynamics being controlled by fluctuations characterized by a large characteristic period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computational Fluid Dynamics Simulations to Assess Spatial Variability and Optimal Ventilation Scenarios for Biological Laboratory Exposures

A significant amount of uncertainty exists regarding potential human exposure to laboratory biomaterials and organisms in Biosafety Level 2 (BSL-2) research laboratories. Computational fluid dynamics (CFD) modeling is proposed as a way to better understand potential impacts of different combinations of biomaterials, laboratory manipulations, and exposure routes on risks to laboratory workers. Here, in this study, we use CFD models to simulate airborne concentrations of contaminants in an actual BSL-2 laboratory under different configurations. Results show that ventilation configuration, sampling location, and contaminant source location can significantly impact airborne concentrations and exposures. Depending on the source location and airflow patterns, the transient and time-integrated concentrations varied by several orders of magnitude. Contaminant plumes from sources located near a return vent (or exhaust like a fume hood or ventilated biosafety cabinet) are likely to be more contained than sources that are further from the exhaust. Having a direct flow between the source and the exhaust (through-flow condition) may reduce potential exposures to individuals outside the air flow path. Designing a BSL-2 room with ventilation and airflow patterns that maximize through-flow conditions to the return/exhaust vents and minimize dispersion and mixing throughout the room is, therefore, recommended. CFD simulations can also be used to assist in characterizing the impacts of supply and return vent locations, room layout, and source locations on spatial and temporal contaminant concentrations. In addition, proper placement of particle sensors can also be informed by CFD simulations to provide additional characterization and monitoring of potential exposures in BSL-2 facilities.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Validation and Comparison of HyRAM Physics Models

The Hydrogen Risk Assessment Models (HyRAM) software version 3 uses a real gas equation of state rather than the Abel-Noble equation of state that is used in 2.0 and previous versions. This change enables the use of HyRAM 3 for cryogenic hydrogen flows, whereas the Abel-Noble equation of state is not accurate at low temperatures. HyRAM 3.1 results were compared to experimental data from the literature in order to demonstrate the accuracy of the physics models. HyRAM 3.1 results were also compared to HyRAM 2.0 for high-pressure, non-cryogenic flows to highlight the differences in predictions between the two major versions of HyRAM. Validation data sets are from multiple groups and span the range of HyRAM physics models, including tank blowdown, unignited dispersion jet plume, ignited jet flame, and accumulation and overpressure inside an enclosure. Both versions 2.0 and 3.1 of HyRAM are accurate for predictions of blowdowns, diffusion jets, and diffusion flames of hydrogen at pressures up to 900 bar, and HyRAM 3.1 also shows good agreement with cryogenic hydrogen data. Overall, HyRAM 3.1 improves on the accuracy of the physical models relative to HyRAM 2.0. In most cases, this reduces the conservatism in risk calculations using HyRAM.

08 HYDROGEN↗
DOE OSTI2021

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

CO2 Plume Imaging with Accelerated Deep Learning-based Data Assimilation Considering Multiple Realizations: Application to the Illinois Basin-Decatur Carbon Sequestration Project

We propose a fast and efficient deep learning workflow for near real-time data assimilation, forecasting and visualization of CO2 plume evolution in saline aquifer and demonstrate its application at a field site. Unlike the previous work, this study incorporates the impact of spatial heterogeneity using multiple realizations. In the proposed workflow, a neural network model utilizes available monitoring data such as downhole pressure measurements as input and predicts the propagating pressure ‘front’ using the diffusive time of flight (DTOF) map which is considered as representative reservoir image of the flow field. The DTOF is the arrival time of pressure front propagation, which can be computed by the Fast Marching Method rapidly without flow simulations. Reservoir model calibration can be implemented by selecting the training data samples that describe the predicted DTOF map based on observed data. The power and efficacy of our workflow is demonstrated by application to the Illinois Basin-Decatur Project.

CO2 plume imaging↗
DOE OSTI2024

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗
DOE OSTI2022

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hero Carbonsafe Phase 2 Project in the Columbia River Basalt Group: Technical Program Overview

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Hermiston is located near the center of the Columbia River Basalt Group (CRBG), which is one of the largest basalt flows in the US. Basalt CO2 storage has potential advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. associated decreases in pressure and CO2 migration risks, 3. reduced long-term monitoring requirements with respect to plume tracking, 4. widespread geographic distribution and, 5. large storage potential due to thickness, porosity, and CO2 interactions with basalt. For locations such as the Pacific Northwest (PNW), Hawaii, Iceland, India and Japan, whose localities are isolated from large sedimentary basins offering conventional saline storage options, basalt may offer the only feasible option for local CO2 storage. However, mineralization/basalt storage still has many uncertainties, as there are limited field-scale assessments of CO2 storage in basalt. There are significant uncertainties hindering the effective implementation of carbon capture utilization and storage (CCUS) in basalt. These include the lack of proven storage capacities, challenges in methodologies for modeling the area of review in igneous formations, limited understanding of mineralization kinetics and timing, and uncertainties in injectivity. Additionally, the domestic availability of specialized services and drilling expertise is constrained, and existing CCUS permitting and regulatory frameworks, originally developed for conventional saline reservoirs, may not adequately address the unique requirements of basalt systems. HERO CarbonSAFE is designed to address major research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts of the CRBG, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts.

58 GEOSCIENCES↗
DOE OSTI2023

Multiwell Fiber Optic Sensing Reveals Effects of CO2 Flow on Triggered Seismicity

Induced seismicity is one of the main risks for gigaton-scale geological storage of carbon dioxide (⁠CO 2 ⁠). Thus, passive seismic monitoring is often recommended as a necessary component of the monitoring systems for CO 2 storage projects, with a particular forcus on risk mitigation. We present the first field study, CO2CRC Otway Project Stage 3 (Victoria, Australia), where distributed acoustic sensing (DAS) enabled high-precision tracking of the induced seismicity triggered by a small CO 2 injection and also informed the reservoir models. In 610 days of passive seismic monitoring of the Stage 3 injection, we detected 17 microseismic events (maximum moment magnitude Mw 0.1) using five deep boreholes equipped with enhanced-sensitivity optical fiber. The DAS array has sensitivity sufficient for detection and location of induced events with M w ~ ₋2 in a monitoring borehole located up to 1500 m away. Thanks to the dense spatial sampling by the DAS, we were able to estimate the focal mechanisms for events with M w > ₋1.5; although the monitoring boreholes provided very limited angular coverage. The main cluster of the events has the same location and source mechanism as the one triggered by the previous CO 2 injection at the Otway Project site, Stage 2C. Surprizingly, the Stage 2C and Stage 3 events closely followed the actual movement of the CO 2 saturation plume front (not the pressure front), as observed using controlled-source reflection seismic images. The nature of the plume-fault interaction remains unclear, but some alteration of the fault gouge by CO 2 might be responsible for the faults’ reactivation by the pressure perturbation. Importantly, the seismogenic fault could not be identified in the seismic images and was only revealed by DAS observations, which also demonstrated the signature of fluid–rock interaction, that may control the CO 2 flow.

58 GEOSCIENCES↗
DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling CO 2 plume migration using an invasion-percolation approach that includes dissolution

Geological carbon dioxide (CO 2 ) sequestration has been proposed as a viable technique to decrease effective emissions of CO 2 into the atmosphere. However, the security of this sequestration is tied to our understanding of the long-term migration of CO 2 in subsurface. The dissolution of CO 2 in the reservoir brine is one of the main long-term trapping mechanisms. However, the assumptions used in large-scale reservoir simulations usually lead to an overestimation of the dissolution volume. We propose a modified approach based on the macroscopic invasion-percolation (MIP) theory that allows the dissolution of CO 2 into brine. We used a high-resolution geological model to compare the Darcy-, modified MIP-, and classic MIP-based simulation results. We observed a significant shrinkage in the nonaqueous plume volume when dissolution is considered during the MIP simulation. In the case of Darcy-based simulation, the plume was completely trapped inside the reservoir with limited migration even after a thousand-year simulation. On the other hand, the majority of the plume migrated out of the simulated reservoir in the case of MIP. Our approach provides more realistic estimation of the dissolution volume and nonaqueous plume extent while leveraging the computational efficiency enjoyed by MIP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Independent Review of Groundwater Remediation Strategy for Hexavalent Chromium and RDX Groundwater Plumes at Los Alamos National Laboratory (Rev. 1)

Site operations at the Los Alamos National Laboratory (LANL) resulted in the release of oxidized chromium, Cr(VI), into Sandia Canyon from cooling tower effluent from 1956 until 1972. The chromium traveled with the surface water approximately 3 miles downstream before migrating below ground surface. Chromium concentrations exceed 50 μg/L in the upper portion of the aquifer. Another LANL groundwater plume of concern is associated with RDX (Royal Demolition Explosives, 1,3,5-trinitro-1,3,5-triazine). Between 1951 and 1996, RDX was released to the mesa-top facilities' process water outfall, adjacent and underlying soils, and alluvial sediments, along with surface water in Cañon de Valle. Between 2000 and 2010, two remedial actions were deployed, removing much of the near-surface RDX, however, recharge due to precipitation has transported RDX into the perched-intermediate zone and into the regional aquifer. The report documents an independent technical review by scientists from the Department of Energy (DOE) Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) to provide recommendations for potential near term actions to address and optimize remediation for both the Cr(VI) and RDX plumes. The proposed near-term remedial actions include design of pump and treat systems for Cr(VI) and monitoring and study for natural attenuation for RDX. The review assesses existing data, conceptual and numerical modeling, and it recommends a technical integration process to support identifying and implementing strategic, effective and efficient remedies. The DOE Environmental Management Los Alamos Field Office (EM-LA) and their cleanup contractor Newport News Nuclear-BWTX, LLC Los Alamos (N3B) provided the information required for the review. Interviews were also conducted with regulators to obtain the full spectrum of technical, regulatory and scientific perspectives. The independent review team was impressed by the capabilities, experiences, innovativeness, and insightfulness of the technical representatives from both the regulator, the New Mexico Environment Department (NMED) and N3B. Incorporation of vadose zone flow pathways in the conceptual site model (CSM) and configuring the numerical modeling for the site was generally state-of-the-practice (or better). This could be considered state-of-the-art by addressing uncertainties related to spatial extent of hydraulic windows. The reviews from the regulators were thorough and often provided useful concepts for consideration and future/study resolutions. The overarching consensus recommendation of independent review team is that the LANL groundwater plumes should be addressed in context of the emerging "management of complex sites" paradigm.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗