Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Waste Repository”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Modeling the Ionic Strength Effect on Diffusion in Clay. The DR-A Experiment at Mont Terri

Solute diffusion in compacted clays depends on ionic strength through its control on the thickness of the electrical double layer (EDL) on the charged clay surfaces. In the DR-A field experiment (Mont Terri, Switzerland), synthetic porewater was circulated through a borehole for 189 days, leading to the out-diffusion of a variety of tracers into the Opalinus Clay. The borehole solution was then replaced with a higher-salinity solution for an additional 540 days, leading to the diffusion of Cs + , Ca 2+ , Mg 2+ , and Sr 2+ back into the borehole and to an increase in the out-diffusion of anions (I - , Br - ) and 3 H. The experimental results were interpreted using the CrunchClay code, which includes a mean electrostatic potential model for the EDL. The EDL corresponds to a second continuum in addition to bulk electrically neutral porewater. Species-specific diffusion (Nernst-Planck equation) occurs through both domains. A 1D radial model considered a single pore diffusion coefficient (Dp = 10-9 m2/s) for cations and 3H in the bulk porosity, and a smaller D p (3 × 10 -10 m 2 /s) for anions. Dp values in the EDL were smaller (10 -11 m 2 /s), except for Cs + and K + (5 × 10 -10 and 2 × 10 -10 m 2 /s, respectively). The model reproduced well the experimental results and showed the capability to consider temporal changes in geochemical conditions affecting the transport and retention of potentially important radionuclide contaminants (e.g., 137 Cs + , 90 Sr 2+ , 129 I - ) in underground geological nuclear waste repositories. Finally, coupled multicomponent diffusion together with the electrostatic properties of the charged surfaces are essential in the development of predictive models for ion transport in clays.

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

Predicting Effective Diffusion Coefficients in Mudrocks Using a Fractal Model and Small-Angle Neutron Scattering Measurements

The determination of effective diffusion coefficients of gases or solutes in the water–saturated pore space of mudrocks is time consuming and technically challenging. Yet, reliable values of effective diffusion coefficients are important to predict migration of hydrocarbon gases in unconventional reservoirs, dissipation of (explosive) gases through clay barriers in radioactive waste repositories, mineral alteration of seals to geological CO 2 storage reservoirs and contaminant migration through aquitards. In this study, small angle and very small angle neutron scattering techniques have been utilized to determine a range of transport properties in mudrocks, including porosity, pore size distributions and surface and volume fractal dimensions of pores and grains, from which diffusive transport parameters can be estimated. Using a fractal model derived from Archie's Law, we calculate effective diffusion coefficients from these parameters and compare them to laboratory–derived effective diffusion coefficients for CO 2 , H 2 , CH 4 and HTO on either the same or related mudrock samples. The samples include Opalinus Shale from the underground laboratory in Mont Terri, Switzerland; Boom Clay from a core drilled in Mol, Belgium and a marine claystone cored in Utah, USA. The predicted values were compared to laboratory diffusion measurements. The measured and modeled diffusion coefficients show good agreement, differing generally by less than factor 5. Furthermore, neutron or X–ray scattering analysis is therefore proposed as a novel method for fast, accurate estimation of effective diffusion coefficients in mudrocks, together with simultaneous measurement of multiple transport parameters including porosity, pore size distributions and surface areas, important for (reactive) transport modelling.

58 GEOSCIENCES↗
DOE OSTI2019

Evolution of the Percolation Threshold in Muds and Mudrocks During Burial

Abstract Muds and mudrocks are important barriers to flow and perform this role over their entire depositional history, from mediating sediment‐ocean methane exchange to preventing leakage from hydrocarbon reservoirs, CO 2 sequestration targets, or nuclear waste repositories. It is well established that the percolation threshold, which is the fluid saturation at which percolation occurs, is directly related to the pressure sealing capacity of a mud or mudrock. However, the evolution of the percolation threshold during burial and diagenesis is not well understood. Using a data set of natural muds and mudrocks and resedimented laboratory mixtures, we show that the percolation threshold is strongly controlled by porosity, achieving a minimum of 0.10–0.30 at a porosity of 0.40, and increasing as porosity increases or decreases from that value, to as high as 0.60 at deposition and >0.80 at porosity <0.10. This implies a highly disconnected pore network at high and low porosities.

Daigle, Hugh↗
DOE OSTI2021

Thermodynamic non-ideality and disorder heterogeneity in actinide silicate solid solutions

Non-ideal thermodynamics of solid solutions can greatly impact materials degradation behavior. We have investigated an actinide silicate solid solution system (USiO 4 –ThSiO 4 ), demonstrating that thermodynamic non-ideality follows a distinctive, atomic-scale disordering process, which is usually considered as a random distribution. Neutron total scattering implemented by pair distribution function analysis confirmed a random distribution model for U and Th in first three coordination shells; however, a machine-learning algorithm suggested heterogeneous U and Th clusters at nanoscale (~2 nm). The local disorder and nanosized heterogeneous is an example of the non-ideality of mixing that has an electronic origin. Partial covalency from the U/Th 5f–O 2p hybridization promotes electron transfer during mixing and leads to local polyhedral distortions. The electronic origin accounts for the strong non-ideality in thermodynamic parameters that extends the stability field of the actinide silicates in nature and under typical nuclear waste repository conditions.

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

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

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

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗
DOE OSTI2022

Frictional properties of Opalinus Clay: influence of humidity, normal stress and grain size on frictional stability

SUMMARY The Opalinus Clay (OPA) is a clay-rich formation considered as a potential host rock for radioactive waste repositories and as a caprock for carbon storage in Switzerland. Its very low permeability (10−19 to 10−21 m2) makes it a potential sealing horizon, however the presence of faults that may be activated during the lifetime of a repository project can compromise the long-term hydrological confinement, and lead to mechanical instability. Here, we have performed laboratory experiments to test the effect of relative humidity (RH), grain size (g.s.) and normal stress on rate-and-state frictional properties and stability of fault laboratory analogues corresponding to powders of OPA shaly facies. The sifted host rock powders at different grain size fractions (<63 μm and 63 < g.s. < 125 μm), at room (∼25 per cent) and 100 per cent humidity, were slid in double-direct shear configuration, under different normal stresses (5–70 MPa). We observe that peak friction, μpeak and steady-state friction, μss, depend on water vapour content and applied normal stress. Increasing relative humidity from ∼25 per cent RH (room humidity) to 100 per cent RH causes a decrease of frictional coefficient from 0.41 to 0.35. The analysis of velocity-steps in the light of rate-and-state friction framework shows that the stability parameter (a–b) is always positive (velocity-strengthening), and it increases with increasing sliding velocity and humidity. The dependence of (a–b) on slip rate is lost as normal stress increases, for each humidity condition. By monitoring the variations of the layer thickness during the velocity steps, we observe that dilation (Δh) is directly proportional to the sliding velocity, decreases with normal stress and is unaffected by humidity. Microstructural analysis shows that most of the deformation is accommodated within B-shear zones, and the increase of normal stress (σn) promotes the transition from strain localization and grain size reduction to distributed deformation on a well-developed phyllosilicate network. These results suggest that: (1) the progressive loss of velocity dependence of frictional stability parameter (a–b) at σn > 35 MPa is dictated by a transition from localized to distributed deformation and (2) water vapour content does not affect the deformation mechanisms and dilation, whereas it decreases steady-state friction (μss), and enhances fault stability.

58 GEOSCIENCES↗
DOE OSTI2015

Linking preferred orientations to elastic anisotropy in Muderong Shale, Australia

The significance of shales as unconventional hydrocarbon reservoirs, nuclear waste repositories, and geological carbon storage has opened new research frontiers in geophysics. Among many of its unique physical properties, elastic anisotropy has long been investigated by both experimental and computational approaches. Here we calculate elastic properties of Cretaceous Muderong Shale from Australia with a self-consistent averaging method based on microstructural information. The volume fraction and crystallographic preferred orientation distributions of constituent minerals are based on synchrotron X-ray diffraction experiments. Aspect ratios of minerals and pores, determined from scanning electron microscopy (SEM), are introduced in the self-consistent averaging. Our analysis suggests that phyllosilicates (i.e., illite-mica, illite-smectite, kaolinite, and chlorite) are dominant with ~70 vol.%. The shape of clay platelets displays an average aspect ratio of 0.05. These platelets are aligned parallel to the bedding plane with a rather high degree of preferred orientation. The estimated porosity at ambient pressure is ~17 vol.% and is divided into equiaxial pores and flat pores with an average aspect ratio of 0.01. Our model shows results (e.g. at pressure of ~50 MPa with C 11 = 26.3; C 13 =12.5; C 33 = 18.2; C 44 = 2.8; C 66 = 6.8 [GPa]) that compare satisfactorily with values derived from ultrasonic velocity measurements (C 11 = 26.6; C 13 = 16.2; C 33 = 18.3; C 44 = 4.5; C 66 = 8.8 [GPa]), confirming the validity and reliability of our approximations and averaging approach.

58 GEOSCIENCES↗
DOE OSTI2016

Molecular Dynamics Simulations of Anion Exclusion in Clay Interlayer Nanopores

© 2016, Clays and Clay Minerals. All rights reserved. The aqueous chemistry of water films confined between clay mineral surfaces remains an important unknown in predictions of radioelement migration from radioactive waste repositories. This issue is particularly important in the case of long-lived anionic radioisotopes ( 129 I - , 99 TcO - 4 , 36 Cl - ) which interact with clay minerals primarily by anion exclusion. For example, models of ion migration in clayey media do not agree as to whether anions are completely or partially excluded from clay interlayer nanopores. In the present study, this key issue was addressed for Cl - using MD simulations for a range of nanopore widths (6 to 15 Å) overlapping the range of average pore widths that exists in engineered clay barriers. The MD simulation results were compared with the predictions of a thermodynamic model (Donnan Equilibrium model) and two pore-scale models based on the Poisson-Boltzmann equation under the assumption that interlayer water behaves as bulk liquid water. The simulations confirmed that anion exclusion from clay interlayers is greater than predicted by the pore-scale models, particularly at the smallest pore size examined. This greater anion exclusion stems from Cl - being more weakly solvated in nano-confined water than it is in bulk liquid water. Anion exclusion predictions based on the Poisson-Boltzmann equation were consistent with the MD simulation results, however, if the predictions included an ion closest approach distance to the clay mineral surface on the order of 2.0 ± 0.8 Å. These findings suggest that clay interlayers approach a state of complete anion exclusion (hence, ideal semi-permeable membrane properties) at a pore width of 4.2 ± 1.5 Å.

Tournassat, Christophe↗
DOE OSTI2023

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗
DOE OSTI2016

A Control Variate Method for Probabilistic Performance Assessment. Improved Estimates for Mean Performance Quantities of Interest

We present a method of control variates for calculating improved estimates for mean performance quantities of interest, E(PQI) , computed from Monte Carlo probabilistic simulations. An example of a PQI is the concentration of a contaminant at a particular location in a problem domain computed from simulations of transport in porous media. To simplify the presentation, the method is described in the setting of a one- dimensional elliptical model problem involving a single uncertain parameter represented by a probability distribution. The approach can be easily implemented for more complex problems involving multiple uncertain parameters and in particular for application to probabilistic performance assessment of deep geologic nuclear waste repository systems. Numerical results indicate the method can produce estimates of E(PQI)having superior accuracy on coarser meshes and reduce the required number of simulations needed to achieve an acceptable estimate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Risk Assessment in Event Evaluation P-302

Incidents at locations such as nuclear power plants, waste repositories, and medical treatment facilities occur at many different times and under a variety of situations. Even the most modern, well run organizations may experience events that are precursors to more serious situations. Currently, the U.S. Nuclear Regulatory Commission evaluates a variety of such events at commercial nuclear plants via the Significance Determination Process (SDP), Management Directive 8.3, and the Accident Sequence Precursor (ASP) program. In these programs, the analyst calculates a conditional risk measure (e.g., conditional core damage probability) for an initiating event situation or an unplanned equipment outage. This calculation is called an “event evaluation.” An “event evaluation” uses a probabilistic risk assessment model to obtain a risk measure that is conditional on the situation which existed during an incident. The Risk Assessment in Event Evaluation material is intended to address how this evaluation is performed using modern risk analysis tools.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2018

M4SF-18LL010302082-Thermodynamic Database Development and Identification of Actinide Sequestration in Corrosion Products

This progress report (Level 4 Milestone Number M4SF-18LL010302082) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Activity Number SF-18LL010302081. The activity is focused on our long-term commitment of engaging our partners in international nuclear waste repository research. This includes participation in the following activities: actinide sequestration in steel corrosion products (Enrica Balboni), Nuclear Energy Agency Thermochemical Database (NEA TDB) Project (Cindy Atkins-Duffin), thermodynamic database collaborations (Tom Wolery), and surface complexation model international collaborations (Mavrik Zavarin).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

M4SF-19LL010302082-International Thermodynamics Activities

This progress report (Level 4 Milestone Number M4SF-19LL010302082) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Activity Number SF-19LL010302081. The activity is focused on our long-term commitment of engaging our partners in international nuclear waste repository research. This includes participation in the Nuclear Energy Agency Thermochemical Database (NEA TDB) Project (Cindy Atkins-Duffin) and development of methodologies for integrating US and international thermodynamic databases for use in SFWST Generic Disposal System Assessment (GDSA) efforts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Baseline Postirradiation Examination of the FUTURIX-FTA Experiments

Postirradiation Examination (PIE) results from the baseline examination of the FUTURIX-FTA irradiation experiment are presented in this report. The FUTURIX-FTA irradiation was designed to study the fuel performance of candidate transmutation fuels. The FUTURIX-FTA experiment was conducted in the true fast neutron spectrum conditions of the Phénix sodium fast test reactor in France. FUTURIX-FTA contained several compositions that were also irradiated as part of the AFC-1 test in cadmium filtered positions in the Idaho National Laboratory (INL) Advanced Test Reactor (ATR) that approximate the neutron spectrum of a fast reactor. This experiment was a joint collaboration between the Department of Energy (DOE) in the U.S. and the Commissariat á l'Energie Atomique et aux Energies Alternatives (CEA) in France. A joint agreement between DOE and CEA was established in 2004 to conduct an irradiation experiment of metallic and nitride fuels in the Phénix reactor at CEA. Both the FUTURIX-FTA tests and the AFC-1 tests were conducted for the DOE Advanced Fuels Campaign (AFC) that is part of Fuel Cycle Research & Development (FCRD) program. These programs seek to develop and demonstrate the technologies needed to transmute long-lived transuranic actinide isotopes contained in spent nuclear fuel via fast reactor technology. The long term goal of this work is to reduce the radiotoxicity of future high level waste repositories. Post-irradiation examination of irradiated fuels experiments provides data related to in reactor fuel performance and input into future fuel design choices. The compositions of the FUTURIX-FTA pins were designed to test if minor actinides could be incorporated into the fuel of a fast reactor without significantly changing fuel performance observed for previously well studied fuels. The fuel performance of metallic fuel in fast reactors has been well documented. Likewise the irradiation performance of nitride fuels in fast reactors has also been documented. However the existing database for nitride fuels is much smaller than the database for metal fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

High Fidelity Hybrid Method for In Situ Borehole Stress Determination Final Report

The state of stress in the earth is complicated and it is difficult to determine all three components and directions of the stress. However, the state of stress affects all activities which take place in the earth, from causing earthquakes on critically stressed faults, to affecting production from hydraulically fractured shale reservoirs, to determining closure rates around a subterranean nuclear waste repository. Current state of the art methods commonly have errors in magnitude and direction of up to 40%. This is especially true for the intermediate principal stress. This project seeks to better understand the means which are used to determine the state of stress in the earth and improve upon current methods to decrease the uncertainty in the measurement. This is achieved by a multipronged experimental investigation which is closely coupled with advanced constitutive and numeric modeling.

58 GEOSCIENCES↗
DOE OSTI2020

Water transport pathway in clay interlayer upon dehydration [Slides]

Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Our research question is, "Microscopically, how does the hydration/dehydration process occur?"

36 MATERIALS SCIENCE↗
DOE OSTI2020

Water transport pathway in clay interlayer upon dehydration [Slides]

Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Microscopically, how does the hydration/dehydration process occur?

36 MATERIALS SCIENCE↗