Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Overlay”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A subzone reconstruction algorithm for efficient staggered compatible remapping

Staggered-grid Lagrangian hydrodynamics algorithms frequently make use of subzonal discretization of state variables for the purposes of improved numerical accuracy, generality to unstructured meshes, and exact conservation of mass, momentum, and energy. For Arbitrary Lagrangian–Eulerian (ALE) methods using a geometric overlay, it is difficult to remap subzonal variables in an accurate and efficient manner due to the number of subzone–subzone intersections that must be computed. This becomes prohibitive in the case of 3D, unstructured, polyhedral meshes. A new procedure is outlined in this paper to avoid direct subzonal remapping. The new algorithm reconstructs the spatial profile of a subzonal variable using remapped zonal and nodal representations of the data. The reconstruction procedure is cast as an under-constrained optimization problem. Enforcing conservation at each zone and node on the remapped mesh provides the set of equality constraints; the objective function corresponds to a quadratic variation per subzone between the values to be reconstructed and a set of target reference values. Numerical results for various pure-remapping and hydrodynamics tests are provided. Ideas for extending the algorithm to staggered-grid radiation-hydrodynamics are discussed as well as ideas for generalizing the algorithm to include inequality constraints.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Center for Frontiers of Subsurface Energy Security (Final Report)

For the past two decades a slowly growing consensus has been that rising temperature is being caused by increases in carbon dioxide (CO 2 ) in the Earth’s atmosphere. What is less consensual is what to do about it. Solutions range from reducing emissions to various ways to capture CO 2 (from effluents of power plants, as plants, etc.) and store it away from the atmosphere. The problem is of such magnitude that it is likely that we will need many methods. Also some methods would cause draconian harm to global economies. This work is about the storage part of the reducing of CO 2 in the atmosphere. There have been many methods proposed for storing CO 2 : in vegetation, at the bottom of oceans, in land farms, and in subsurface formations. The latter is particularly prominent given the large estimated volume in porous formation in the Earth, and the 50-year old experience in injecting CO 2 for enhanced oil recovery. This report is about the science behind subsurface storage, called here geological carbon storage or GCS. Effective GCS requires that billions of metric tons of CO 2 be injected into storage reservoirs annually. GCS feasibility requires balancing the opposing impacts of large emplacement (injection) rates without compromising caprock seal integrity during injection. Once injected, the challenge is to assure the injected CO 2 remains permanently stored. Under most conditions, CO 2 is in a supercritical state, or scCO 2 . scCO 2 is immiscible with the subsurface brine and has a factor of ten smaller viscosity, leading to the development of two-phase fluid flow and scCO 2 bypassing. After the injection of large volumes of scCO 2 , the pore pressure, state of stress, chemical, thermal, and biological steady-state or equilibrium conditions in the subsurface are disrupted, which may lead to unpredictable behavior because the responses to these changes are often non-linear. Current challenges for efficient and reliable scCO 2 storage are: (1) sustaining large injection rates; (2) using available pore space efficiently, and (3) controlling undesired or unexpected behavior, such as scCO 2 leakage. Various aspects of these challenges are addressed by collaborative projects within Center for Frontiers in Subsurface Energy Security (CFSES). Challenge 1 – Sustaining large storage rates. To offset the annual CO 2 emissions of in the United States, ~7 billion metric tons of CO 2 must be injected into the storage reservoirs annually. Injected CO 2 forms a plume that starts spreading, displacing originally resident brine, though inefficiently so. Because of its low viscosity, and buoyant forces, a CO 2 plume rises to the reservoir rock-caprock interface. Injecting large volumes of scCO 2 at sufficiently high rates into subsurface storage formations leads to increases in pore pressures, potential expansion of the reservoir rock, and may result in fracturing of the reservoir rock and/or caprock. CFSES research uses existing and new experimental and modeling approaches to identify chemical-physical controls on the permeability and pore pressure dynamics, anticipating geomechanical (fracturing) events, predicting multi-phase flow patterns and trapping of CO 2 . Challenge 2 – Using pore space with unprecedented efficiency. Since large volumes of scCO 2 are to be stored, the volume available in storage reservoirs must be used efficiently. Storage efficiency is defined as a fraction of pore space occupied by carbon dioxide. Current estimates are that less than 5 percent of available pore volume is available for storage because scCO 2 , having a smaller density and viscosity than reservoir brine, tends to form viscous fingers, and bypass available pore space. This flow regime results in a CO 2 plume spreading over large areas in subsurface. The research goal in Challenge 2 is to advance the fundamental understanding of two-phase flow through porous media, including molecular-to meter- scale experimental and modeling approaches. This work leads to novel strategies to ensure that injected CO 2 occupies more than 5% of the pore space, ideally up to 50%. Challenge 3 – Controlling undesired or unexpected behavior. Following the injection of scCO 2 , the subsurface reservoir re-equilibrates to different physical and chemical states. This re-equilibration in space and time could lead to unexpected and/or undesired behavior: fracturing of reservoir rock and/or caprock, development of preferential flow paths, CO 2 leakage (to the surface or overlaying fresh water reservoirs), induced seismicity, and activation of faults. One of the main research goals of CFSES was to advance the fundamental understanding of how individual chemical and physical processes are coupled. We particularly focused on the chemical-mechanical coupling, and flow-mechanical coupling. We collected experimental and modeling data to determine the temporal and spatial scales of the coupled processes, which may lead to pore collapse caused by creep, activation of fractures and faults, chemically-induced fracture propagation, and self-focusing of scCO 2 flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗
DOE OSTI2021

Visualization, understanding, and mitigation of process-induced-membrane irregularities in gas diffusion electrode-based polymer electrolyte membrane fuel cells

Polymer electrolyte membrane fuel cells (PEMFC) show substantial promise for their application in electric vehicles. For large-scale manufacturing of PEMFCs, roll-to-roll coated gas-diffusion-electrodes (GDE) offer certain advantages over other production pathways. Procedures including hot pressing and coating an ionomer overlayer may be necessary for this manufacturing pathway to enable a suitable catalyst layer/membrane interface. The same procedures may potentially introduce membrane irregularities, especially when thin membranes are used. Limited understanding exists regarding if and to what extent such irregularities impact PEMFC performance and lifetime, and therefore be considered defects. In this study, NREL's customized fuel cell hardware that enables quasi in-situ infrared (IR) thermography studies was utilized to visualize spatial hydrogen crossover and identify membrane irregularities that originated from the GDE-based MEA fabrication process. The structure of these membrane irregularities was investigated by scanning electron microscopy (SEM) and its impact on initial H2/air performance was determined. Accelerated stress testing (AST) revealed that these irregularities develop into failure point locations. These results were validated across many MEAs with identified process-induced membrane irregularities. By selecting specific gas diffusion media properties and by fine tuning the MEA hot pressing parameters, the formation of such membrane irregularities was mitigated.

08 HYDROGEN↗
DOE OSTI2026

Atomic-Scale Characterization and Electronic Properties of Gold-Capped Niobium Films for Superconducting Qubits

Niobium thin films are central to superconducting qubits, but their complex native oxides contribute significantly to microwave losses. A promising mitigation strategy is to suppress oxide formation using engineered thin-film encapsulation layers. Recent experiments have shown that ∼10 nm metallic overlayers on Nb capacitor films can improve the energy relaxation time T1 of transmon qubits. Here, we present a comparative study of Au-capped Nb films fabricated in situ by molecular beam epitaxy and ex situ by sequential deposition benchmarked against bare Nb films. Using complementary structural, chemical, and spectroscopic techniques, we correlate the interface quality with superconducting electronic properties of the Au surface. Low-temperature scanning tunneling spectroscopy (STS) provides spatially resolved quasiparticle density-of-states maps. While both Au-capped Nb samples exhibit large areas with uniform, fully gapped density of states, clear differences emerge between the two interfaces. Compared to in situ Nb-Au, ex situ Nb-Au exhibits a reduced induced superconducting gap, broadened coherence peak, and localized in-gap states, consistent with a residual NbxOy layer at the interface. Supported by complementary structural and chemical analyses, these findings demonstrate that STS directly links nanoscale superconducting properties to interface preparation, highlighting the importance of controlled in situ encapsulation for minimizing dissipation and improving quantum coherence.

Makita, Junki [MIT; Temple U.] (ORCID:000000026889↗
DOE OSTI2025

Spatially Enhanced Electrostatic Doping in Graphene Realized via Heterointerfacial Precipitated Metals

Forming heavily-doped regions in 2D materials, like graphene, is a steppingstone to the design of emergent devices and heterostructures. Here, a selective-area approach is presented to tune the work-function and carrier density in monolayer graphene by spatially synthesizing sub-monolayer gallium beneath the 2D-solid. The localized metallic gallium is formed via precipitation from an underlying diamond-like carbon (DLC) film that is spatially implanted with gallium-ions. By controlling the interfacial precipitation process with annealing temperature, spatially precise ambipolar tuning of the graphene work-function is achieved, and the tunning effect preserved upon cooling to ambient conditions. Consequently, charge carrier densities from ≈1.8 × 10 10 cm -2 (hole-doped) to ≈7 × 10 13 cm -2 (electron-doped) are realized, confirmed by in situ and ex situ measurements. The theoretical studies corroborated the role of gallium at the heterointerface on charge transfer and electrostatic doping of the graphene overlayer. Specifically, sub-monolayer gallium facilitates heavy n-doping in graphene. Extending this doping strategy to other implantable elements in DLC provides a new means of exploring the physics and chemistry of highly-doped 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin Frustration and Unconventional Surface Spin Canting State in Van der Waals Ferromagnet/Antiferromagnet Heterostructures

Atomically flat surfaces of van der Waals (vdW) materials pave an avenue for addressing a long-standing fundamental issue of how a compensated antiferromagnet (AFM) surface frustrates a ferromagnetic (FM) overlayer in FM/AFM heterostructures. We investigate Fe5GeTe2/NiPS3 vdW heterostructures by characterizing AFM and FM spins separately. We find that in-plane zig-zag AFM NiPS3 develops three equivalent AFM domains, which are robust against external magnetic field and magnetic coupling with Fe5GeTe2. Moreover, evidence is provided of in-plane-AFM-induced perpendicular magnetic anisotropy (PMA) in adjacent Fe5GeTe2, and an unconventional out-of-plane surface spin canting state with the Fe5GeTe2 spins spatially turn from out-of-plane direction near the interface to in-plane direction away from the interface in Fe5GeTe2/NiPS3. The out-of-plane surface spin canting is a unique property of spin frustration in vdW magnetic heterostructures.

Wang, Tianye↗
DOE OSTI2019

Biogeographic classification of streams using fish community– and trait–environment relationships

Abstract Aim Present a hybrid biogeographic and environmentally constrained clustering approach to classify ~853,000 stream reaches in the eastern United States. Examine the frequency of stream typologies in the landscape relative to anthropogenic stressors to identify potential stream conservation needs. Location Eastern United States. Methods Fish communities at 956 least‐disturbed sampling reaches were characterized using taxonomic and functional composition of native species. Environmental variables summarized within stream reaches included stream discharge, channel gradient, hydrologic regime, summer water temperature and boundaries of freshwater ecoregions. Multivariate regression trees were used to relate environmental variables to fish communities while simultaneously developing taxa‐ and trait‐specific biogeographic classifications. We then overlaid stream classifications with indices of anthropogenic disturbances to evaluate rarity and risk of complete loss of natural representation of stream types. Results Taxonomically based classes represented a combination of ecoregional boundaries and environmental gradients, whereas temperature gradients were most important in differentiating functional composition. An optimal taxonomy‐based classification contained 13 classes and explained 26.4% of community variation, and the optimal trait‐based classification contained 6 classes, explaining 25.2% of community variation. Overlaying class maps with an anthropogenic disturbance index revealed substantial variation in the severity of degradation among classes relative to their abundance in the landscape, especially rare classes (i.e., large rivers and cold headwater systems). Main conclusions We demonstrate the value of integrating local, regional and historical biogeographic factors with taxonomic and functional characteristics of freshwater biota to classify ecosystems based on mechanistic biota–environment relationships. This approach is important to addressing the challenge of developing stream classifications that reconcile the environmental drivers of community trait composition with the constraints of biogeographic regionality at large spatial scales. Furthermore, the spatially explicit and high‐resolution maps disseminated here can uniquely inform management decisions regarding anthropogenic threats to the biodiverse freshwaters of the eastern United States.

Troia, Matthew J.↗
DOE OSTI2026

Towards Spatial Control of Reaction Selectivity on Photocatalysts using Area Selective Atomic Layer Deposition on Model Dual Site Electrocatalyst Platform

Photocatalytic water splitting is a promising route to low-cost H2; however, this approach is currently limited by solar-to-hydrogen (STH) conversion efficiencies in the sub-10% range. Z-scheme water splitting, in which H2 and O2 evolving particles are operated in separate compartments and electronically coupled by a soluble redox mediator, offers improvements in STH efficiency but introduces high rates of undesired side and back reactions. Nanoscopic oxide coatings (e.g., CrOx, and SiOx) have previously been used to selectively block undesired reactants from reaching active sites; however, a coating encapsulating the entire photocatalyst particle limits activity as the particle can no longer facilitate both half reactions. In contrast to photodeposition, which may produce non-uniform overlayers due to the difficulty of controlling local electrochemical reactivity on the surface, area-selective atomic layer deposition (AS-ALD) can be used to deposit conformal ultrathin oxide coatings while preserving access to non-growth areas of a substrate surface. To develop this technique for Z-scheme photocatalysts, we performed AS-ALD of TiO2 on a dual site planar electrocatalyst based on interdigitated arrays of Pt and Au. Self-assembled monolayers of 1-octadecanethiol (ODT) were used to block growth on the Au array, resulting in a patterned surface in which only the Pt sites were encapsulated with TiO2. These electrocatalysts demonstrated localized reaction selectivity for hydrogen evolution reaction (HER) over the TiO2/Pt sites, while the uncoated Au sites retained activity towards Fe(II)/Fe(III) redox (FeRR). Under independent potential control, these microelectrodes showed that selectively deposited TiO2 coatings can suppress the rate of back reactions on neighboring active sites by an order of magnitude compared to uncoated control samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Solution-Phase Synthesis of PdH 0.706 Nanocubes with Enhanced Stability and Activity toward Formic Acid Oxidation

We report palladium is one of the few metals capable of forming hydrides, with the catalytic properties being dependent on the elemental composition and spatial distribution of H atoms in the lattice. Herein, we report a facile method for the complete transformation of Pd nanocubes into a stable phase made of PdH 0.706 by treating them with aqueous hydrazine at a concentration as low as 9.2 mM. Using formic acid oxidation (FAO) as a model reaction, we systematically investigated the structure–catalytic property relationship of the resultant nanocubes with different degrees of hydride formation. The current density at 0.4 V was enhanced by four times when the nanocubes were completely converted from Pd to PdH 0.706 . On the basis of a set of slab models with PdH(100) overlayers on Pd(100), we conducted density functional theory calculations to demonstrate that the degree of hybrid formation could influence both the activity and selectivity toward FAO by modulating the relative stability of formate (HCOO) and carboxyl (COOH) intermediates. Furthermore, this work provides a viable strategy for augmenting the performance of Pd-based catalysts toward various reactions without altering the loading of this scarce metal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Quantum dots with split enhancement gate tunnel barrier control

We introduce a silicon metal-oxide-semiconductor quantum dot architecture based on a single polysilicon gate stack. The elementary structure consists of two enhancement gates separated spatially by a gap, one gate forming a reservoir and the other a quantum dot. We demonstrate, in three devices based on two different versions of this elementary structure, that a wide range of tunnel rates is attainable while maintaining single-electron occupation. A characteristic change in slope of the charge transitions as a function of the reservoir gate voltage, attributed to screening from charges in the reservoir, is observed in all devices, and is expected to play a role in the sizable tuning orthogonality of the split enhancement gate structure. The all-silicon process is expected to minimize strain gradients from electrode thermal mismatch, while the single gate layer should avoid issues related to overlayers (e.g., additional dielectric charge noise) and help improve yield. Finally, reservoir gate control of the tunnel barrier has implications for initialization, manipulation and readout schemes in multi-quantum dot architectures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A genomics-informed study of oxalate and cellulase regulation by brown rot wood-degrading fungi

Wood-degrading fungi that selectively remove carbohydrates (brown rot) combine Fenton-based oxidation and enzymatic hydrolysis to degrade wood. These two steps are incompatible in close proximity. To explain this, brown rot fungi may stagger oxidative reactions ahead of hydrolysis, but the scale and environmental controls for such a mechanism have not been resolved in solid wood. Here, we focused on one reaction control parameter, oxalate. In coordination with Fe 3+ -reducing compounds (e.g., 2,5-dimethoxyhydroquinone), oxalate can either promote Fenton chemistry by mobilizing Fe 3+ as mono-oxalates (facilitative) or inhibit Fenton chemistry (protective) by restricting reducibility and the formation of Fenton’s reagent as Fe 3+ /Fe 2 -(oxalate) 2,3 . We sectioned wood wafers colonized directionally by Postia placenta and Gloeophyllum trabeum to map end-to-end the expression of oxalate synthesis genes and to overlay enzyme activities, metabolites, and wood modifications. Near advancing hyphal fronts, oxaloacetase expression was up upregulated for both fungi, while regulation patterns of paralogous of isocitrate lyases and glyoxylate dehydrogenases varied, suggesting different physiological roles. Oxalate decarboxylase (ODC) expression in G. trabeum was induced in more decayed wood behind the hyphal front, but was constitutively expressed in all P. placenta sections. Relative ODC activities increased and oxalate levels stabilized in more decayed wood behind the hyphal front. Endoglucanase (EG) activity, on the other hand, peaked for both fungi in later decay stages. These oxalate optimization patterns are in line with previous whole-block ‘spiking’ experiments tracking oxalate, but we provide here information on its genetic controls across a spatial gradient. As a complement, we also demonstrate in vitro the plausibility of a protective role for oxalate, to emphasize that these fungi might be optimizing oxalate at a given level to maximize Fenton reactions but to minimize oxidative damage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Detection of magnetic circular dichroism in amorphous materials utilizing a single-crystalline overlayer

Physicists are fascinated with topological defects in solid-state materials, because by breaking the translational symmetry they offer emerging properties that are not present in their parental phases. For example, edge dislocations—the 2π phase-winding topological defects—in antiferromagnetic NiO crystals can exhibit ferromagnetic behaviors. Herein, we study how these defects could give rise to exotic topological orders when they interact with a high energy electron beam. To probe this interaction, we formed a coherent electron nanobeam in a scanning transmission electron microscope and recorded the far-field transmitted patterns as the beam steps through the edge dislocation core in [001] NiO. Surprisingly, we found the amplitude patterns of the <020> Bragg disks evolve in a similar manner to the evolution of an annular solar eclipse. Using the ptychographic technique, we recovered the missing phase information in the diffraction plane and revealed the topological phase vortices in the diffracted beams. Through atomic topological defects, the wave function of electrons can be converted from plane wave to electron vortex. This approach provides a new perspective for boosting the collection efficiency of magnetic circular dichroism spectra with high spatial resolution and understanding the relationship between symmetry breaking and exotic property of individual topological defect at atomic level.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Selective-Area Remote Epitaxy of ZnO Microrods Using Multilayer–Monolayer-Patterned Graphene for Transferable and Flexible Device Fabrications

Here, selective-area remote epitaxy (SA-REpi) is demonstrated for fabricating mechanically releasable position-controlled ZnO microrod (MR) arrays from donor wafers in an arrayed form. Intaglio-patterned graphene, consisting of basal single-layer graphene (SLG) overlayered with multilayer graphene (MLG) patterned with perforated holes, is transferred onto a GaN/Al 2 O 3 wafer on which the hydrothermal synthesis is performed for growing ZnO MRs. The basal SLG area exposed through the MLG pattern yields ZnO MRs, whereas the MLG plateau inhibits the growth. The noncovalent remote epitaxial heterointerface enables the release of the MR overlayer in the arrayed form, and the original source wafer is refurbished for reproducibly repeating the SA-REpi. Density-functional theory calculations suggest that localized surface charge density is induced on the surface of SLG by the underlying GaN across ultrathin SLG, which possibly provides a driving force for precursor adatoms and the following remote epitaxy of ZnO. In contrast, the induction of the charge density redistribution does not clearly occur through MLG; so, that keeps the surface of MLG nearly charge-neutral. The diameter and spacing of ZnO MRs are controlled in a designed way by changing the pattern geometries. High-resolution scanning transmission electron microscopy reveals the remote heteroepitaxial relationship at an atomic level. The remote epitaxy is expected to provide an ideal platform to transfer the addressable spatial arrays of nano- or micro-architecture semiconductor components to arbitrary target surfaces directly after the growth without the assembly procedures.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Innovative Use of Cr(VI) Plume Depictions and Pump-and-Treat Capture Analysis to Estimate Risks of Contaminant Discharge to Surface Water at Hanford Reactor Areas

The Hanford Site nuclear reactor operations required large quantities of high-quality cooling water, which was treated with chemicals including sodium dichromate dihydrate for corrosion control. Cooling water leakage, as well as intentional discharge of cooling water to ground during upset conditions, produced extensive groundwater recharge mounds consisting largely of contaminated cooling water and resulted in wide distribution of hexavalent chromium (Cr[VI]) contamination in the unconfined aquifer. The 2013 Cr(VI) groundwater plumes in the 100 Areas cover approximately 6 km2 (1500 acres), primarily in the 100-HR-3 and 100-KR-4 groundwater operable units (OUs). The Columbia River is a groundwater discharge boundary; where the plumes are adjacent to the Columbia River there remains a potential to discharge Cr(VI) to the river at concentrations above water quality criteria. The pump-and-treat systems along the River Corridor are operating with two main goals: 1) protection of the Columbia River, and 2) recovery of contaminant mass. An evaluation of the effectiveness of the pump-and-treat systems was needed to determine if the Columbia River was protected from contamination, and also to determine where additional system modifications may be needed. In response to this need, a technique for assessing the river protection was developed which takes into consideration seasonal migration of the plume and hydraulic performance of the operating well fields. Groundwater contaminant plume maps are generated across the Hanford Site on an annual basis. The assessment technique overlays the annual plume and the capture efficiency maps for the various pump and treat systems. The river protection analysis technique was prepared for use at the Hanford site and is described in detail in M.J. Tonkin, 2013. Interpolated capture frequency maps, based on mapping dynamic water level observed in observation wells and derived water levels in the vicinity of extraction and injection wells, are developed initially. Second, simulated capture frequency maps are developed, based on transport modelling results. Both interpolated and simulated capture frequency maps are based on operation of the systems over a full year. These two capture maps are then overlaid on the plume distribution maps for inspection of the relative orientation of the contaminant plumes with the capture frequency. To quantify the relative degree of protection of the river from discharges of Cr(VI) (and conversely, the degree of threat) at any particular location, a systematic method of evaluating and mapping the plume/capture relationship was developed. By comparing the spatial relationship between contaminant plumes and hydraulic capture frequency, an index of relative protectiveness is developed and the results posted on the combined plume/capture plan view map. Areas exhibiting lesser degrees of river protection are identified for remedial process optimization actions to control plumes and prevent continuing discharge of Cr(VI) to the river.

Miller, Chuck W.↗
DOE OSTI2015

Innovative use of Cr(VI) plume depictions and pump-and-treat capture analysis to estimate risks of contaminant discharge to surface water at Hanford reactor areas - 15315

The Hanford Site nuclear reactor operations required large quantities of high-quality cooling water, which was treated with chemicals including sodium dichromate dihydrate for corrosion control. Cooling water leakage, as well as intentional discharge of cooling water to ground during upset conditions, produced extensive groundwater recharge mounds consisting largely of contaminated cooling water and resulted in wide distribution of hexavalent chromium (Cr[VI]) contamination in the unconfined aquifer. The 2013 Cr(VI) groundwater plumes in the 100 Areas cover approximately 6 km{sup 2} (1500 acres), primarily in the 100-HR-3 and 100-KR-4 groundwater operable units (OUs). The Columbia River is a groundwater discharge boundary; where the plumes are adjacent to the Columbia River there remains a potential to discharge Cr(VI) to the river at concentrations above water quality criteria. The pump-and-treat systems along the River Corridor are operating with two main goals: 1) protection of the Columbia River, and 2) recovery of contaminant mass. An evaluation of the effectiveness of the pump-and-treat systems was needed to determine if the Columbia River was protected from contamination, and also to determine where additional system modifications may be needed. In response to this need, a technique for assessing the river protection was developed which takes into consideration seasonal migration of the plume and hydraulic performance of the operating well fields. Groundwater contaminant plume maps are generated across the Hanford Site on an annual basis. The assessment technique overlays the annual plume and the capture efficiency maps for the various pump and treat systems. The river protection analysis technique was prepared for use at the Hanford site and is described in detail in M.J. Tonkin, 2013. Interpolated capture frequency maps, based on mapping dynamic water level observed in observation wells and derived water levels in the vicinity of extraction and injection wells, are developed initially. Second, simulated capture frequency maps are developed, based on transport modelling results. Both interpolated and simulated capture frequency maps are based on operation of the systems over a full year. These two capture maps are then overlaid on the plume distribution maps for inspection of the relative orientation of the contaminant plumes with the capture frequency. To quantify the relative degree of protection of the river from discharges of Cr(VI) (and conversely, the degree of threat) at any particular location, a systematic method of evaluating and mapping the plume/capture relationship was developed. By comparing the spatial relationship between contaminant plumes and hydraulic capture frequency, an index of relative protectiveness is developed and the results posted on the combined plume/capture plan view map. Areas exhibiting lesser degrees of river protection are identified for remedial process optimization actions to control plumes and prevent continuing discharge of Cr(VI) to the river. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Smoothed particle hydrodynamics model for Landau-Lifshitz-Navier-Stokes and advection-diffusion equations

We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗