Publications

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.

A well-field within a uranium (U) plume in the groundwater-surface water transition zone was monitored for a 3 year period for water table elevation and dissolved solutes. The plume discharges to the Columbia River, which displays a dramatic spring stage surge resulting from snowmelt. Groundwater exhibits a low hydrologic gradient and chemical differences with river water. River water intrudes the site in spring. Specific aims were to assess the impacts of river intrusion on dissolved uranium (Uaq), specific conductance (SpC), and other solutes, and to discriminate between transport, geochemical, and source term heterogeneity effects. Time series trends for Uaq and SpC were complex and displayed large temporal and well-to-well variability as a result of water table elevation fluctuations, river water intrusion, and changes in groundwater flow directions. The wells were clustered into subsets exhibiting common behaviors resulting from the intrusion dynamics of river water and the location of source terms. Hot-spots in Uaq varied in location with increasing water table elevation through the combined effects of advection and source term location. Heuristic reactive transport modeling with PFLOTRAN demonstrated that mobilized Uaq was transported between wells and source terms in complex trajectories, and was diluted as river water entered and exited the groundwater system. While Uaq time-series concentration trends varied significantly from year-to-year as a result of climate-caused differences in the spring hydrograph, common and partly predictable response patterns were observed that were driven by water table elevation, and the extent and duration of river water intrusion.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report describes specific GDSA activities in fiscal year 2015 (FY2015) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code (Hammond et al., 2011) and the Dakota uncertainty sampling and propagation code (Adams et al., 2013). Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through the engineered barriers and natural geologic barriers to a well location in an overlying or underlying aquifer. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

Abstract not provided.

Simulations of thermal-hydrology were carried out for the emplacement of spent nuclear fuel canisters and cesium and strontium capsules using the PFLOTRAN simulator. For the cesium and strontium capsules the analysis looked at disposal options such as different disposal configurations and surface aging of waste to reduce thermal effects. The simulations studied temperature and fluid flux in the vicinity of the borehole. Simulation results include temperature and vertical flux profiles around the borehole at selected depths. Of particular importance are peak temperature increases, and fluxes at the top of the disposal zone. Simulations of cesium and strontium capsule disposal predict that surface aging and/or emplacement of the waste at the top of the disposal zone reduces thermal effects and vertical fluid fluxes. Smaller waste canisters emplaced over a longer disposal zone create the smallest thermal effect and vertical fluid fluxes no matter the age of the waste or depth of emplacement.

Abstract not provided.

Deep geological disposal of nuclear waste in clay/shale/argillaceous rock formations has received much consideration given its desirable attributes such as isolation properties (low permeability), geochemically reduced conditions, slow diffusion, sorbtive mineralogy, and geologically widespread (Jové Colón et al., 2014). There is a wealth of gained scientific expertise on the behavior of clay/shale/ argillaceous rock given its focus in international nuclear waste repository programs that includes underground research laboratories (URLs) in Switzerland, France, Belgium, and Japan. Jové Colón et al. (2014) have described some of these investigative efforts in clay rock ranging from site characterization to research on the engineered barrier system (EBS). Evaluations of disposal options that include nuclear waste disposition in clay/shale/argillaceous rock have determined that this host media can accommodate a wide range of waste types. R&D work within the Used Fuel Disposition Campaign (UFDC) assessing thermal effects and fluid-mineral interactions for the disposition of heat-generating waste have so far demonstrated the feasibility for the EBS and clay host rock to withstand high thermal loads. This report represents the continuation of disposal R&D efforts on the advancement and refinement of coupled Thermal-Hydrological-Mechanical-Chemical (THMC), hydrothermal experiments on clay interactions, used fuel degradation (source term), and thermodynamic modeling and database development. The development and implementation of a clay/shale/argillite reference case described in Jové Colón et al. (2014) for FY15 will be documented in another report (Mariner et al. 2015) – only a brief description will be given here. This clay reference case implementation is the result of integration efforts between the GDSA PA and disposal in argillite work packages. The assessment of sacrificial zones in the EBS is being addressed through experimental work along with 1D reactive-transport and reaction path modeling. The focus of these investigations into the nature of sacrificial zones is to evaluate the chemical effects of heterogeneous chemical reactions at EBS interfaces. The difference in barrier material types and the extent of chemical reactions within these interfacial domains generates changes in mineral abundances. These mineralogical alterations also result in volume changes that, although small, could affect the interface bulk porosity. As in previous deliverables, this report is structured according to various national laboratory contributions describing R&D activities applicable to clay/shale/argillite media.

Abstract not provided.

Abstract not provided.

This paper describes and demonstrates two methods of providing a priori information to the surface-based time-lapse three-dimensional electrical resistivity tomography (ERT) problem for monitoring stage-driven or tide-driven surface water intrusion into aquifers. First, a mesh boundary is implemented that conforms to the known location of the water table through time, thereby enabling the inversion to place a sharp bulk conductivity contrast at that boundary without penalty. Second, a nonlinear inequality constraint is used to allow only positive or negative transient changes in EC to occur within the saturated zone, dependent on the relative contrast in fluid electrical conductivity between surface water and groundwater. A 3-D field experiment demonstrates that time-lapse imaging results using traditional smoothness constraints are unable to delineate river water intrusion. The water table and inequality constraints provide the inversion with the additional information necessary to resolve the spatial extent of river water intrusion through time.

Abstract not provided.

This study develops a new framework of facies-based data assimilation for characterizing spatial distribution of hydrofacies and estimating their associated hydraulic properties. This framework couples ensemble data assimilation with transition probability-based geostatistical model via a parameterization based on a level set function. The nature of ensemble data assimilation makes the framework efficient and flexible to be integrated with various types of observation data. The transition probability-based geostatistical model keeps the updated hydrofacies distributions under geological constrains. The framework is illustrated by using a two-dimensional synthetic study that estimates hydrofacies spatial distribution and permeability in each hydrofacies from transient head data. Our results show that the proposed framework can characterize hydrofacies distribution and associated permeability with adequate accuracy even with limited direct measurements of hydrofacies. Our study provides a promising starting point for hydrofacies delineation in complex real problems.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this study, we use PFLOTRAN, a highly scalable, parallel, flow, and reactive transport code to simulate the concentrations of 3H, 3He, CFC-11, CFC-12, CFC-113, SF6, 39Ar, and the mean groundwater age in heterogeneous fields on grids with an excess of 10 million nodes. We utilize this computational platform to simulate the concentration of multiple tracers in high-resolution, heterogeneous 2D and 3D domains, and calculate tracer-derived ages. Tracer-derived ages show systematic biases toward younger ages when the groundwater age distribution contains water older than the maximum tracer age. The deviation of the tracer-derived age distribution from the true groundwater age distribution increases with increasing heterogeneity of the system. However, the effect of heterogeneity is diminished as the mean travel time gets closer to the tracer age limit. Age distributions in 3D domains differ significantly from 2D domains. 3D simulations show decreased mean age, and less variance in age distribution for identical heterogeneity statistics. High-performance computing allows for investigation of tracer and groundwater age systematics in high-resolution domains, providing a platform for understanding and utilizing environmental tracer and groundwater age information in heterogeneous 3D systems.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

directly, rather than through simplified abstractions. It also a llows for complex representations of the source term, e.g., the explicit representation of many individual waste packages (i.e., meter - scale detail of an entire waste emplacement drift). This report fulfills the Generic Disposal System Analysis Work Packa ge Level 3 Milestone - Performance Assessment Modeling and Sensitivity Analyses of Generic Disposal System Concepts (M 3 FT - 1 4 SN08080 3 2 ).

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Threats to water distribution systems include release of contaminants and Denial of Service (DoS) attacks. A better understanding, and validated computational models, of the flow in water distribution systems would enable determination of sensor placement in real water distribution networks, allow source identification, and guide mitigation/minimization efforts. Validation data are needed to evaluate numerical models of network operations. Some data can be acquired in real-world tests, but these are limited by 1) unknown demand, 2) lack of repeatability, 3) too many sources of uncertainty (demand, friction factors, etc.), and 4) expense. In addition, real-world tests have limited numbers of network access points. A scale-model water distribution system was fabricated, and validation data were acquired over a range of flow (demand) conditions. Standard operating variables included system layout, demand at various nodes in the system, and pressure drop across various pipe sections. In addition, the location of contaminant (salt or dye) introduction was varied. Measurements of pressure, flowrate, and concentration at a large number of points, and overall visualization of dye transport through the flow network were completed. Scale-up issues that that were incorporated in the experiment design include Reynolds number, pressure drop across nodes, and pipe friction and roughness. The scale was chosen to be 20:1, so the 10 inch main was modeled with a 0.5 inch pipe in the physical model. Controlled validation tracer tests were run to provide validation to flow and transport models, especially of the degree of mixing at pipe junctions. Results of the pipe mixing experiments showed large deviations from predicted behavior and these have a large impact on standard network operations models.3