Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

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

The EGS Collab project: Outcomes and lessons learned from hydraulic fracture stimulations in crystalline rock at 1.25 and 1.5 km depth

With the goal of better understanding stimulation in crystalline rock for improving enhanced geothermal systems (EGS), the EGS Collab Project performed a series of stimulations and flow tests at 1.25 and 1.5 km depths. The tests were performed in two well-instrumented testbeds in the Sanford Underground Research Facility in Lead, South Dakota, United States. The testbed for Experiment 1 at 1.5 km depth contained two open wells for injection and production and six instrumented monitoring wells surrounding the targeted stimulation zone. Four multi-step stimulation tests targeting hydraulic fracturing and nearly year-long ambient temperature and chilled water flow tests were performed in Experiment 1. The testbed for Experiments 2 and 3 was at 1.25 km depth and contained five open wells in an outwardly fanning five-spot pattern and two fans of well-instrumented monitoring wells surrounding the targeted stimulation zone. Experiment 2 targeted shear stimulation, and Experiment 3 targeted low-flow, high-flow, and oscillating pressure stimulation strategies. Hydraulic fracturing was successful in Experiments 1 and 3 in generating a connected system wherein injected water could be collected. However, the resulting flow was distributed dynamically, and not entirely collected at the anticipated production well. Thermal breakthrough was not observed in the production well, but that could have been masked by the Joule-Thomson effect. Shear stimulation in Experiment 2 did not occur - despite attempting to pressurize the fractures most likely to shear - because of the inability to inject water into a mostly-healed fracture, and the low shear-to-normal stress ratio. The EGS Collab experiments are described to provide a background for lessons learned on topics including induced seismicity, the correlation between seismicity and permeability, distributed and dynamic flow systems, thermoelastic and pressure effects, shear stimulation, local geology, thermal breakthrough, monitoring stimulation, grouting boreholes, modeling, and system management.

15 - GEOTHERMAL ENERGY↗

DECOVALEX-2019: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX Project is an international research collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative, now running for almost 30 years, was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel, and also for a range of other sub-surface engineering activities. DECOVALEX emphasizes joint analysis and comparative modeling of state-of-the-art field and laboratory experiments, across a range of host rock options and repository designs. Participating research teams are from radioactive waste management organizations, national research institutes, regulatory agencies, universities, and consulting groups, providing a wide range of perspectives and solutions to these complex problems. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2019, started in 2016 and ended in 2019. Modeling teams from 13 international partner organizations participated in the comparative evaluation of seven modeling tasks involving complex field and/or laboratory experiments. Furthermore, this Virtual Special Issue on DECOVALEX-2019 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes.

58 GEOSCIENCES↗

DECOVALEX-2023: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX initiative is an international research collaboration (www.decovalex.org), initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of the complex perturbations and coupled processes in geologic repositories and how these impact long-term performance predictions. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2023, started in early 2020 and ended in late 2023. More than fifty research teams associated with 17 international DECOVALEX partner organizations participated in the comparative evaluation of eight modeling tasks covering a wide range of spatial and temporal scales, geological formations, and coupled processes. This Virtual Special Issue on DECOVALEX-2023 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes. While primarily focused on radioactive waste, much of the work included here has wider application to many geoengineering topics.

Coupled processes↗

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗

Brine Availability Test in Salt (BATS) FY24 Update

This report summarizes fiscal year 2024 (FY24) activities centered around a series of field tests in bedded salt at the Waste Isolation Pilot Plant (WIPP) funded by the Office of Spent Fuel and Waste Science and Technology in the Spent Fuel and Waste Disposition (SFWD) program of the US Department of Energy’s Office of Nuclear Energy (DOE-NE). High-level Purpose of Experiments: The Brine Availability Test in Salt (BATS) field tests are revealing both brine occurrence (i.e., where, and how much) and brine migration (i.e., how easily it moves) in the excavation damaged zone (EDZ). This understanding is foundational to develop a safety case for a future heat-generating waste repository in salt, and to starting up a generic repository program in salt to buy down risk. BATS seeks to predict how much brine can flow into both ambient and heated excavations (e.g., boreholes or rooms) in salt. This work is educating and empowering new repository scientists on two fronts: “design and execution of field tests” and “prediction and modeling of coupled processes.” DOE-NE capabilities in salt have grown and been tested through international modeling and benchmarking exercises (e.g., DECOVALEX, RANGERS, KOMPASS, and MEASURES; see Mills et al., 2024). The hands-on expertise we are building is a necessary step towards large-scale disposal demonstrations and eventual implementation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Roles of Tidal Marshes in the Estuarine Biochemical Processes: A Numerical Modeling Study

Observations suggest that the existence of tidal marsh can alter the oxygen and nutrient dynamics in adjacent water bodies, but assessing the impacts of large tidal marshes on an estuary is challenging. In this study, we use a modeling approach to investigate the roles of tidal marshes on the estuarine biochemical processes. The marsh model, which simulates the ecological functions of marshes at seasonal and annual time-scales, is embedded inside an unstructured-grid three-dimensional hydrodynamic and eutrophication model (SCHISM-ICM). This modeling system simulates the growth and metabolism of the tidal marshes and links biological processes to nutrient dynamics in the water column and sediment. This model dynamically simulates nutrient recycling and physical transport of the materials between marshes and open water through wetting-drying processes. This coupled model system is validated and successfully applied to the York River Estuary. Model results suggest that tidal marshes influence the local diurnal dissolved oxygen (DO) cycle by exporting dissolved organic carbon and high sediment oxygen demand in the marsh system through the tidal exchange. The high deposition rates of organics and diurnal DO cycle enhance the sediment release of phosphorus. On the other hand, marshes tend to decrease dissolved inorganic nitrogen in the water column by settling particulate nutrients and enhancing the denitrification process. The study demonstrates that tidal marshes exert substantial impacts on the estuarine biochemical processes. The developed tidal marsh model enhances eutrophication modeling and advances the understanding of the feedback effects between marsh biogeochemistry and estuarine eutrophication processes on a systemic scale.

54 ENVIRONMENTAL SCIENCES↗

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

Accurate assessment of land–atmosphere coupling in climate models requires high-frequency data output

Land–atmosphere (L–A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centered climate anomalies. Local L–A coupling (LoCo) metrics capture relevant L–A processes, highlighting the impact of soil and vegetation states on surface flux partitioning and the impact of surface fluxes on boundary layer (BL) growth and development and the entrainment of air above the BL. A primary goal of the Climate Process Team in the Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) project is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional Earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L–A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land–atmosphere interactions span timescales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability in behavioral regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here, we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyze, and better understand L–A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behavior. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant timescales, (ii) enable demonstration of inter-model spread in L–A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

54 ENVIRONMENTAL SCIENCES↗

Dynamical phase-field model of coupled electronic and structural processes

Abstract Many functional and quantum materials derive their functionality from the responses of both their electronic and lattice subsystems to thermal, electric, and mechanical stimuli or light. Here we propose a dynamical phase-field model for predicting and modeling the dynamics of simultaneous electronic and structural processes and the accompanying mesoscale pattern evolution under static or ultrafast external stimuli. As an illustrative example of application, we study the transient dynamic response of ferroelectric domain walls excited by an ultrafast above-bandgap light pulse. We discover a two-stage relaxational electronic carrier evolution and a structural evolution containing multiple oscillational and relaxational components across picosecond to nanosecond timescales. The phase-field model offers a general theoretical framework which can be applied to a wide range of functional and quantum materials with interactive electronic and lattice orders and phase transitions to understand, predict, and manipulate their ultrafast dynamics and rich mesoscale evolution dynamics of domains, domain walls, and charges.

36 MATERIALS SCIENCE↗

Hydro-mechanical behavior of heated bentonite buffer for geologic disposal of high-level radioactive waste: A bench-scale X-ray computed tomography investigation

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Further, compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Cloud Process Coupling and Time Integration in the E3SM Atmosphere Model

Abstract In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 s. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step and, in particular, on the frequency with which cloud‐related processes are coupled with each other and with the model dynamics.

54 ENVIRONMENTAL SCIENCES↗

Near Repository Unsaturated Alluvium Disposal Modeling with Improved Geological Realism

The Disposal Research and Development (R&D) Program of the US Department of Energy (DOE) office of Nuclear Energy (NE-8) Spent Fuel and Waste Science and Technology (SFWST) Campaign is to conduct R&D on disposal of spent nuclear fuel (SNF) and high-level waste (HLW). The goal of the Geologic Disposal Safety Assessment (GDSA) within this project is to develop a disposal system modeling and analysis capability that supports the integrated modeling of coupled processes controlling disposal system performance of deep geologic repositories, including uncertainty. This report describes a specific activity in the Fiscal Year 2024 (FY24) associated with the GDSA Repository Systems Analysis (RSA) work package in collaboration with the GDSA Geologic Modeling work package at Los Alamos National Laboratory (LANL). The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and repository system performance models in crystalline, argillite, salt, and unsaturated alluvium potential host-rock environments, and to simulate and analyze these generic repository concepts and models using GDSA Framework toolkit, and other tools as needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data for publication Wan et al. (2020, JAMES), Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling

This dataset contains simple condensation mdoel output and CAM4 model output that were used to produce figures in the following paper: Wan, H., Woodward, C. S., Zhang, S., Vogl, C. J., Stinis, P., Gardner, D. J., Rasch, P. J., Zeng, X., Larson, V. E., and Singh, B. (2020): Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001982

54 ENVIRONMENTAL SCIENCES↗

GDSA Repository Systems Analysis Investigations in FY 2024

The Disposal Research and Development (R&D) Program of the US Department of Energy (DOE) office of Nuclear Energy (NE-8) Spent Fuel and Waste Science and Technology (SFWST) Campaign is to conduct R&D on disposal of spent nuclear fuel (SNF) and high-level waste (HLW). The goal of the Geologic Disposal Safety Assessment (GDSA) within this project is to develop a disposal system modeling and analysis capability that supports the integrated modeling of coupled processes controlling disposal system performance of deep geologic repositories, including uncertainty. This report describes specific activities in the Fiscal Year (FY) 2024 associated with the GDSA Repository Systems Analysis (RSA) work package. The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and repository system performance models in crystalline, argillite, salt, and unsaturated alluvium potential host-rock environments, and to simulate and analyze these generic repository concepts and models using GDSA Framework toolkit, and other tools as needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analytical Modeling of Biomass Transport and Feeding Systems

The processing of biomass solids in a biorefinery consists of pretreatment, enzyme hydrolysis / concurrent fermentation of sugars to ethanol, product recovery, and drying. Sustainable operation requires a front end that transforms wet solids into a pumpable slurry. Otherwise the biorefinery will suffer unscheduled shut-downs and inefficient operation due to solids that obstruct pumps and other equipment and resist mixing in a bioreactor. Downtime in pioneer biorefineries due to interruptions from materials handling problems has been 50% or more, leading to unsustainable manufacturing processes. This work addresses new technology, predictive computational models, and definition of operational conditions that result in formation of slurries of corn stover at up to 300 g/L using low enzyme loadings (1 to 3 FPU cellulase/g) before the biomass (corn stover) enters the pretreatment step. A team of researchers from Purdue University, Idaho National Laboratory (INL), Forest Concepts, AdvanceBio, Argonne National Laboratory, and DOE BETO have combined their knowledge in agricultural and biological engineering, bioprocess engineering, mechanical engineering, chemical engineering, agricultural economics, materials engineering and enzyme and microbial technology to address the challenge of making lignocellulose flow. This team effort has resulted in the development and validation of conditions that employ low levels of commercial enzyme in an agitated bioreactor to which corn stover pellets are added resulting in formation of slurries at high solids loadings, before pretreatment. This approach overcomes challenges caused by handling of dry, particulate biomass materials at the front end of the biorefinery. The subsequent materials handling issues cause obstruction at pumps, pipes and valves. Formation of high loadings slurries with low yield stress, as reported here, significantly decreases the potential for process interruption and enhances plant operability. Key advances in the knowledge of how slurry formation occurs is reported here and in recently published journal papers. We found that pellets are needed to achieve high solids loading, and that commercial enzymes are effective in forming slurries of corn stover particles from pellets that have not been pretreated. Our work has resulted in models that predict solids behavior for formation of compressed solids and pellets that in turn facilitate slurries made of high concentrations of corn stover particles. A computational model was developed that gives mechanistic insights into properties of particles and mixing process that gives the slurry rheology needed to facilitate pumping. Hence, the corn stover may be pumped into a pretreatment reactor in place of auguring in solids against high pressure which is a root cause of interruptions at the front end of a biorefinery. Subsequent mixing in enzyme and microbial bioreactors results in conversion of lignocellulose to sugars in a biorefinery in agitated bioreactors, with flows in and out of the vessels being less likely to be interrupted due to plugging or materials handling problems. The obtained data coupled to process models, techno-economic assessment (TEA) and Life Cycle Analysis (LCA) were used to assess whether this approach is practical. These results are based on a foundation of laboratory characterization and pilot runs. The NREL biochemical sugar model was utilized to carry out techno-economic analysis of enzyme catalyzed liquefaction followed by enzyme hydrolysis. The minimum sugar selling price was between 17.5 and 18.3 ¢/pound or about the same as calculated by the NREL model for dilute acid pretreatment followed by enzyme hydrolysis. Life cycle analysis (LCA) based on Argonne’s Greet Model showed the enzyme catalyzed route had the lowest greenhouse gas emissions of the three combinations studied (i.e., enzyme, enzyme mimetic, and enzyme + mimetic combined). GHG emissions for enzyme-based corn stover liquefaction step, alone, were about 21 g CO 2 -equivalent/kg of liquefied slurry. We believe this approach will further enhance operability of a pioneer biorefinery, and bring large-scale conversion of lignocellulosic biomass to low carbon footprint biofuels closer to implementation.

09 BIOMASS FUELS↗

Reduced-order modeling of near-field THMC coupled processes for nuclear waste repositories in shale

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke’s law models) are implemented to describe mechanical characteristics of the system. Overall, our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗