Publications

Abstract not provided.

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

Non-uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium (3H), chlorofluorocarbon-12 (CFC-12), and sulfur hexafluoride (SF6) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady-state 2-D synthetic aquifer and a transient 3-D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post-calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration-constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. This model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.

Abstract not provided.

Abstract not provided.

This document is a high-level test plan for small-scale experimental activities in boreholes in salt at WIPP.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report is the Task E specification (Revision 0) for DECOVALEX-2023. Task E is focused on understanding thermal, hydrological, mechanical and chemical (THMC) processes, especially related to predicting brine migration in heated salt. The main test case being used is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) in Carlsbad, New Mexico. This report provides short motivational background, a summary of relevant experiments and data, and a step-by-step plan for the analysis by the teams participating in Task E (Rev. 0 includes detailed description of steps 0 and 1). This document will be revised, and more detail will be added to later steps during DECOVALEX-2023.

Abstract not provided.

Abstract not provided.