SEARCH · Engineering Papers

Results for “subsurface diffusion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗

DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

DOE OSTI2026

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗

DOE OSTI2021

Rapid Alloy Surface Engineering through Closed-Vessel Reagent Pyrolysis

For rapid surface engineering of Cr-containing alloys by low-temperature nitrocarburization, we introduce a process based on pyrolysis of solid reagents, e.g., urea, performed in an evacuated closed vessel. Upon heating to temperatures high enough for rapid diffusion of interstitial solute, but low enough to avoid second-phase precipitation, the reagent is pyrolyzed to a gas atmosphere containing molecules that (i) activate the alloy surface by stripping away the passivating Cr 2 O 3 -rich surface film (diffusion barrier) and (ii) rapidly infuse carbon and nitrogen into the alloy. We demonstrate quantitatively that this method can generate a subsurface zone with concentrated carbon and nitrogen comparable to what can be accomplished by established (e.g., gas-phase- or plasma-based) methods, but with significantly reduced processing time. As another important difference to established gas-phase processing, the interaction of gas molecules with the alloy surface can have auto-catalytic effects by altering the gas composition in a way that accelerates solute infusion by providing a high activity of HNCO. The new method lends itself to rapid experimentation with a minimum of laboratory equipment.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Drought-induced changes in groundwater-surface water exchange at Lake Mead area

This study focuses on the Lake Mead region in the southwestern United States, a key water reservoir serving over 25 million people and agricultural lands across several states. The area has experienced recurring anthropogenic droughts since the early 2000s. We investigate the hydrological response of the coupled surface-groundwater system to the 2020–2022 drought, one of the most severe on record. To this end, we use Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data to quantify vertical land motion caused by the elastic response of the crust to water-mass loss in the lake vicinity. Next, we apply an inverse elastic load modeling framework to quantify water loss. We further assess possible hydraulic connectivity between Lake Mead and adjacent groundwater reservoirs. We detect ground uplift of up to 8 mm/yr near the lake center, likely due to crustal rebound from reduced water-mass loading. We estimated the total water storage loss at 3.03 ± 0.25 km 3 /yr across a 3150 km 2 area surrounding Lake Mead. Groundwater accounts for approximately a third of that, being 0.94 ± 0.32 km 3 /yr. In addition, observed time lags of 6–98 days between lake and groundwater level responses, corresponding to a lateral diffusivity of 3.2–86 m 2 /s, suggest spatially variable connectivity between the lake and aquifers. These findings highlight that drought impacts propagate through the subsurface within interconnected systems, resulting in reduced buffering capacity of groundwater resources following droughts, and emphasizing the need for more integrated surface and groundwater management strategies to enhance resilience under climate and anthropogenic stressors.

58 GEOSCIENCES↗

DOE OSTI2020

Deciphering atomistic mechanisms of the gas-solid interfacial reaction during alloy oxidation

Gas-solid interfacial reaction is critical to many technological applications from heterogeneous catalysis to stress corrosion cracking. A prominent question that remains unclear is how gas and solid interact beyond chemisorption to form a stable interphase for bridging subsequent gas-solid reactions. Here, we report real-time atomic-scale observations of Ni-Al alloy oxidation reaction from initial surface adsorption to interfacial reaction into the bulk. We found distinct atomistic mechanisms for oxide growth in O 2 and H 2 O vapor, featuring a “step-edge” mechanism with severe interfacial strain in O 2 , and a “subsurface” one in H 2 O. Ab initio density functional theory simulations rationalize the H 2 O dissociation to favor the formation of a disordered oxide, which promotes ion diffusion to the oxide-metal interface and leads to an eased interfacial strain, therefore enhancing inward oxidation. Our findings depict a complete pathway for the Ni-Al surface oxidation reaction and delineate the delicate coupling of chemomechanical effect on gas-solid interactions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗

DOE OSTI2020

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗

DOE OSTI2022

Evaluation of subsurface transport processes of delayed gas signatures applicable to underground nuclear explosions

Radioactive gas signatures from underground nuclear explosions (UNEs) result from gas-migration processes occurring in the subsurface. The processes considered in this study either drive or retard upward migration of gases from the detonation cavity. The relative importance of these processes is evaluated by simulating subsurface transport in a dual-permeability medium for the multi-tracer Noble Gas Migration Experiment (NGME) originally intended to study some aspects of transport from a UNE. For this experiment, relevant driving processes include weak two-phase convection driven by the geothermal gradient, over pressuring of the detonation cavity, and barometric pumping while gas sorption, dissolution, radioactive decay, and usually diffusion represent retarding processes. From deterministic simulations we found that over-pressuring of the post-detonation chimney coupled with barometric pumping produced a synergistic effect amplifying the tracer-gas reaching the surface. Bounding simulations indicated that the sorption and dissolution of gases, tending to retard transport, were much smaller than anticipated by earlier laboratory studies. The NGME observations themselves show that differences in gas diffusivity have a larger effect on influencing upward transport than do the combined effects of tracer-gas sorption and dissolution, which is consistent with a Sobol’ sensitivity analysis. Both deterministic simulations and those considering parametric uncertainties of transport-related properties predict that the excess in concentration of SF 6 compared to 127 Xe as might be captured in small volumetric samples should be much smaller than the order-of-magnitude contrast found in the large-volume gas samples taken at the site. While extraction of large-volume subsurface gas samples is shown to be capable of distorting in situ gas compositions, the highly variable injection rate of SF 6 into the detonation cavity relative to that of 127 Xe at the start of the field experiment is the most likely explanation for the large difference in observed concentrations.

58 GEOSCIENCES↗

DOE OSTI2025

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗

DOE OSTI2023

Degeneracy in Intercalated Pb Phases under Buffer-Layer Graphene on SiC(0001) and Diffuse Moiré Spots in Surface Diffraction

First-principles density functional theory (DFT) is used to analyze the stability of Pb intercalated phases under buffer layer graphene on SiC(0001) as a function of the supercell size, Pb coverage, and degree of Pb ordering. By comparing the chemical potentials of such two-dimensional Pb structures, here we find that there is a family of structurally distinct thermodynamically preferred Pb subsurface configurations with minute stability differences. These differences are comparable to the thermal energies at about 450 °C, where the Pb intercalated phases are grown. High-resolution surface-diffraction experiments using Spot Profile Analysis Low-Energy Electron Diffraction (SPA-LEED) confirm this high degree of degeneracy of the Pb intercalated phases from broad, low-intensity moiré spots observed exclusively from intercalated Pb. The low intensity of the moiré spots implies the coexistence of structurally different subsurface Pb phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Using Radiogenic Noble Gas Nuclides to Identify and Characterize Rock Fracturing

Abstract Fracture‐released radiogenic noble gas nuclides are used to identify locations and constrain the volume of new fracture creation during subsurface detonations. Real‐time, in situ noble gases and reactive gases were monitored using a field‐deployed mass spectrometer and automated sampling system in a multilevel borehole array. Released gases were measured after two different detonations having distinct energy, pressure, and gas volume characteristics. Explosive‐derived gases (N 2 O, CO 2 ) and excess radiogenic 4 He and 40 Ar above atmospheric background are used to identify locations of gas transport and new fracture creation after each detonation. Fracture‐released radiogenic 4 He is used to constrain the volume of newly created fractures with a model of helium release from fracturing. Explosive by‐product gas was observed in multiple locations both near and distal to the shot locations for both detonations. Radiogenic 4 He and 40 Ar release from rock damage was observed in locations near the detonation after the second, more powerful detonation. Observed 4 He response is consistent with a model of diffusive release from newly created fractures. Volume of new fractures estimated from the 4 He release ranges from 1 to 5 m 2 with apertures ranging from 0.1 to 1 m. Our results provide evidence that radiogenic noble gases released during fracture creation can be identified at the field scale in real time and used to identify timing and location of fracture creation during deformation events. This technique could be useful in subsurface science and engineering problems where the location and amount of newly created rock fracturing is of interest including fault rupture, mine safety, subsurface detonation monitoring and reservoir stimulation.

58 GEOSCIENCES↗

DOE OSTI2023

Effect of the Molecular Structure of Surface Vanadia on Activity and Regenerability of VO $x$ /In 2 O 3 Catalysts for CO 2 -Assisted Oxidative Dehydrogenation of Propane

Our recent work has reported that higher propylene selectivity and improved stability can be achieved by combining redox-active VO $x$ and basic In 2 O 3 for CO 2 -assisted oxidative dehydrogenation of propane (CO 2 -ODHP). In the present work, we continued to explore the stability and regenerability of V/In catalysts. In particular, our interest lies in identifying the effect of mono- and polyvanadate on catalytic performance and regenerability. A V/In catalyst with an increased proportion of monovanadate was prepared using the Schlenk line under moisture-free conditions (V/In–S), while the fully polymerized vanadate catalyst was prepared through a regular impregnation (V/In) for comparison. The Schlenk-line-prepared catalyst, namely, V/In–S, not only exhibits a 17–30% enhanced propylene yield at high temperatures (500–540 °C) over V/In but also presents improved stability and regenerability with nearly 88% activity recovered after regeneration in O 2 . Detailed characterizations have been performed to reveal the catalyst structure–performance relationship, including chemisorption (NH 3 /CO 2 -temperature-programmed desorption, NH 3 /CO 2 -TPD), H 2 -temperature-programmed reduction (H 2 -TPR), and spectroscopic studies [Raman spectroscopy, UV–vis diffuse reflectance spectroscopy (UV–vis DRS), near-ambient-pressure X-ray photoelectron spectroscopy (NAP-XPS), and high-sensitivity low-energy ion scattering (HS-LEIS)]. Characterization results demonstrate that compared with polyvanadates, monovanadates lead to strengthened interaction with In 2 O 3 and a more stabilized V/In surface and subsurface, as well as improved redox properties of VO $x$ . These advantages give rise to the observed enhancement in activity, stability, and regenerability. In conclusion, these findings advance the understanding of the relationship between the activity/stability and the molecular structure of surface oxide species (vanadia) and the interplay between acid–base interactions and redox properties of mixed metal-oxide catalysts for efficient CO 2 -ODHP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

RAIS Paducah Background Values

The Paducah Background Values (https://rais.ornl.gov/tools/pgdp_background.html) are from Methods for Conducting Risk Assessments and Risk Evaluations at the Paducah Gaseous Diffusion Plant, Volume 1, Human Health, which sought to establish provisional background values for groundwater and soil. These values were developed by considering agreements reached between the DOE and the regulatory agencies during comment resolution meetings, in the Federal Facility Agreement, and at technical meetings. The Paducah background concentrations are available for soil (surface and subsurface) and groundwater (filtered and unfiltered).

Dolislager, Fred [Oak Ridge National Laboratory (O↗

DOE OSTI2021

Strengthening of Calcite Assemblages Through Chemical Complexation Reactions

Abstract Utilization of subsurface reservoirs, including fluid extraction/injection, can induce stress changes and modify in‐situ chemical equilibrium causing subcritical fracturing and deformation. Here, we show that chemical complexation reactions at the newly created fracture surfaces within the crack tip limit fracturing during consolidation of granular calcite assemblages. Previously, we showed that analogous chemical complexation reactions can increase in‐situ fracture toughness in pre‐fractured single crystal calcite. To test the chemical complexation effect under subsurface conditions, we consolidated samples under increasing hydrostatic pressure to induce widespread intergranular fracturing in the presence of fluid containing common anions. Measured consolidation, acoustic emissions, and microfracturing demonstrate that deformation correlates with complexation affinity for calcium in calcite and aqueous anions. We observed that deformation is lowest in the presence of strongly complexing anions (e.g., sulphate), provided sufficient anion diffusion rates. This understanding is important for predicting reservoir effective strength in the presence of complex aqueous fluids.

Choens, R. C.↗

DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗