Engineering topics

Jordan, Amy

Publications and source records attributed to Jordan, Amy.

8 records

DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Preliminary Radiological Risk Assessment Model for Disposition of Remote-Handled Transuranic Wastes at Los Alamos National Laboratory Area G - 20116

The U.S. Department of Energy (DOE) operates a low-level radioactive waste (LLW) disposal site at Material Disposal Area G, in Los Alamos, New Mexico, USA. Area G has been the primary LLW disposal site for Los Alamos National Laboratory (LANL) since the 1960's. In addition to LLW, Area G is host to a variety of other wastes, the disposition of which must be determined before closure of the site. A probabilistic Radiological Risk Assessment (RRA) for Area G is used in order to support decision making regarding some wastes that are not addressed in the extant Area G Performance Assessment (PA) and Composite Analysis (CA). Between 1979 and 1987, 33 special shafts were augered into the Bandelier Tuff at Area G. This volcanic tuff is present across Pajarito Plateau on the eastern slopes of the Jemez Mountains, and varies widely in its consistency, from weakly indurated non-welded layers to welded layers that uphold the mesa cliffs of the Plateau. These mesas are home to LANL, Area G, and the townsites of Los Alamos and White Rock, with residences about 1400 m from Area G. The 33 Shafts were lined with steel casing, and contain remote-handled (RH) transuranic wastes (TRU) resulting from experiments and analysis performed in special glove boxes at the Chemistry and Metallurgy Research (CMR) facility at LANL. Some of these wastes originated as used nuclear fuel. The purpose of the Area G RRA is to evaluate the potential future risk to humans and the environment from the RH TRU in the 33 Shafts in the context of the risk associated with the surrounding wastes at Area G. The analysis is responsive to expectations outlined in DOE Order 458.1, Radiation Protection of the Public and the Environment, and is informed by the Manual and Guidance accompanying DOE O 435.1, Radioactive Waste Management. Because the waste meets the definition of TRU, the regulatory context necessarily takes into consideration the regulation governing the disposal of TRU from the U.S. Environmental Protection Agency (EPA): 40 CFR 191, Environmental Radiation Protection Standards for Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Radioactive Wastes. Given the broader regulatory context for the RRA, the analysis is subject to different assumptions from those made in the existing DOE O 435.1 PA and CA, such as allowing for future occupation of the site. The analysis begins with a comprehensive evaluation of features, events, processes, and exposure scenarios (FEPS) for Area G and the wastes it contains. These FEPSs are screened to eliminate from further consideration those of extremely low probability and/or consequence, and a conceptual site model (CSM) is subsequently developed. The scope and structure of the Area G RRA Model is informed by this CSM, and the Area G RRA Model is developed using the GoldSim systems analysis modeling platform. This paper presents the initial version of a defensible, transparent, and reasonably realistic model, which is based on the state of knowledge of the wastes, the site, and the FEPSs that govern contaminant transport from wastes into the environment and subsequent exposures to humans and other biota. Probabilistic model input distributions represent uncertainties inherent in the real and modeled systems. The results of the Area G RRA Model inform decisions regarding the disposition of the RH TRU in the 33 Shafts. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

A Preliminary Probabilistic Performance Assessment Model for Radiological and Chemical Contamination at the West Valley Site, New York - 19295

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States at the WNYNSC. NFS built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (640 metric tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations halted in 1972 and never restarted, leaving behind radioactive and chemical wastes. The reprocessing wastes include high-level waste (HLW), so U.S. Department of Energy (DOE) was compelled to complete certain waste management activities under the West Valley Demonstration Project (WVDP) Act of 1980, including decommissioning of WVDP facilities. As collaborating agencies, NYSERDA and the DOE are tasked with making decisions about decommissioning and risk reduction for the West Valley Site. Neptune and Company, Inc. (Neptune) was contracted to develop a probabilistic performance assessment (PPA) model to assist the agencies in their decision making process for decommissioning the WVDP and WNYNSC. One important tool that is needed in order to inform the decision-making process is a science-based model of the contamination at the West Valley Site and its potential implications for human health and the environment. The West Valley PPA Model (the Model), developed using the GoldSim modeling software, is a tool intended to provide support for decision making with uncertainty, in a manner that is transparent, defensible, and robust. The PPA Model includes contaminant transport and health effects components, and is organized around geographically-organized contaminated facilities. These include the waste disposal areas licensed by the U.S. Nuclear Regulatory Commission and the State of New York, a waste tank farm for storage of HLW resulting from reprocessing operations, and several areas contaminated with radioactive and chemical constituents. The Model evaluates contaminant transport from these sources to points of exposure across the site and into receiving surface waters and sediments. Hypothetical people and wildlife are exposed to contamination at these locations, and the effects of these exposures are evaluated. Contaminant transport processes to be evaluated in the Model include groundwater and surface water transport, contaminant translocation by plants and animals, diffusion in subsurface air and water, and erosion of the Site. The evaluation of exposures to people in this preliminary model is limited to a resident farmer scenario, and ecological assessment is performed at the level of a screening analysis. The results of these preliminary evaluations inform future model developments. Model results are subjected to sensitivity analysis in order to determine those pathways and parameters that are most significant in influencing the results. This information allows analysts and decision makers to focus on those aspects of Site behavior and processes that contribute most to overall uncertainty. With this information, the Supplemental Environmental Impact Statement (SEIS) can inform defensible decisions regarding remediation of the Site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Explicitly Accounting for Multiple Spatial Scales in Data Informing Kd Distribution Development - 18414

This work provides an example of incorporating information about spatial scale into input parameter distribution development for probabilistic performance assessment (PA) models and highlights the importance of having an interdisciplinary team working together in the distribution development process. The novel approach is demonstrated for contaminant partitioning coefficients for water/rock (Kds), providing a specific context for discussing the approach and challenges, but the strategy is general and applies to many other PA model inputs. In this example, a PA model of a radioactive waste disposal site is assumed to have large cells representing spatial scales of up to thousands of square meters in area and several to tens of meters in depth. Spatial variability in Kds is expected to occur at scales much smaller than the units in the model domain over which values are applied, and therefore values are intended to represent possible averages over the modeled volume (i.e., spatial scale represented by a PA model cell). Information about possible Kd values for each material category comes from a variety of literature sources, and the reported results used to inform the distribution are available at a wide range of scales (e.g., individual batch experiment results from small-scale samples, averages based on samples collected over site-scales, averages based on samples from large regions, etc.). Effectively addressing the mismatch in spatial scale between available information and the PA Model is a challenge within PA distribution development and it takes effort from database managers, subject matter experts, modelers, and statisticians. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Transient Unsaturated Zone Modeling of Radioactive Waste Pits at Los Alamos National Laboratory - 87544 - 18432

Waste disposals at Material Disposal Area G (MDA G) at Los Alamos National Laboratory (LANL) occur primarily in a series of long, rectangular pits into which layers of waste have been emplaced beginning in the late 1950's. MDA G is situated in north-central New Mexico on top of Mesita del Buey, a narrow 3 km long mesa in the Parajito Plateau. The mesa is composed primarily of volcanic tuff from two principle volcanic eruptions in the Jemez Mountains which deposited the Otowi and Tshirege Members of the Bandelier Tuff. The Cerros del Rio Basalt underlay the tuff units and extends down for approximately 200 meters. The water table lies within the basalts at around 300 m below the surface of the mesa. The local climate is temperate and semiarid, averaging around 34 cm of precipitation per year. Ambient subsurface moisture conditions are dry, with water content typically below 10%. During disposal operations, open pits are subject to enhanced infiltration compared to the final covered configuration due to the unvegetated disturbed surface and precipitation events. Additionally, water added to the pits during disposal operations for tasks like waste compaction and dust suppression introduces additional moisture to the system during disposal operations that could potentially provide a pathway for contaminants in the waste zone to migrate through the unsaturated zone. Modeling was conducted using the HYDRUS-2D/3D software to evaluate the spatial extent of saturation changes and fluxes across geologic layer interfaces in the unsaturated zone throughout the life cycle of a waste pit. A general 2-dimensional framework was developed that included a cross section of a waste pit and the underlying geologic units. Several geological layerings were developed to represent the variation across the mesa at Area G. The model was run stochastically for hundreds of realizations to understand how variations in hydraulic parameters (e.g., porosity, saturated hydraulic conductivity, Van Genuchten parameters, etc.) and operational parameters (e.g., pit open period duration, pit depth, initial saturation of the waste zone, etc.) affect the extent of saturation changes and fluxes in the deep unsaturated zone. Distributions for these parameters were developed based on available site data. Results of the unsaturated zone modeling indicate that a significant amount of water can reach the lower unsaturated zone under some circumstances, particularly when the amount of water introduced during pit operations is high. A statistical study of modeling results was performed to develop a regression that describes the variations of water flux in time at the geological layer interfaces. These statistical model results will be used to evaluate solute transport in a future effort. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

FEHM V3.0 - V3.X

The numerical background of the FEHM computer code can be traced to the early 1970s when it was used to simulate geothermal and hot dry rock reservoirs. The primary use over a number of years was to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the potential Yucca Mountain repository. Today FEHM is used to simulate groundwater and contaminant flow and transport in deep and shallow, fractured and un-fractured porous media throughout the US DOE complex. FEHM has proved to be a valuable asset on a variety of projects of national interest including Environmental Remediation of the Nevada Test Site, the LANL Groundwater Protection Program, geologic CO2 sequestration, Enhanced Geothermal Energy (EGS) programs, Oil and Gas production, Nuclear Waste Isolation, and Arctic Permafrost. Subsurface physics has ranged from single fluid/single phase fluid flow when simulating basin scale groundwater aquifers to complex multifluid/ multi-phase fluid flow that includes phase change with boiling and condensing in applications such as unsaturated zone surrounding nuclear waste storage facility or leakage of CO2/brine through faults or wellbores. The numerical method used in FEHM is the control volume method (CV) for fluid flow and heat transfer equations which allows FEHM to exactly enforce energy/mass conservation; while an option is available to use the finite element (FE) method for displacement equations to obtain more accurate stress calculations. In addition to these standard methods, an option to use FE for flow is available, as well as a simple Finite Difference scheme.

Pawar, Rajesh↗