Publications

Abstract not provided.

This document describes the implementation level changes in the source code and input files of Sandia National Laboratories Environmental Fluid Dynamics Code (SNL-EFDC) that are necessary for including pH effects into algae-growth dynamics. The document also gives a brief introduction to how pH effects are modeled into the algae-growth model. The document assumes that the reader is aware of the existing algae-growth model in SNL-EFDC. The existing model is described by James, Jarardhanam and more theoretical considerations behind modeling pH effects are presented therein. This document should be used in conjunction with the original EFDC manual and the original water-quality manual.

Increasing interest in marine hydrokinetic (MHK) energy has spurred to significant research on optimal placement of emerging technologies to maximize energy conversion and minimize potential effects on the environment. However, these devices will be deployed as an array in order to reduce the cost of energy and little work has been done to understand the impact these arrays will have on the flow dynamics, sediment-bed transport and benthic habitats and how best to optimize these arrays for both performance and environmental considerations. An "MHK-friendly" routine has been developed and implemented by Sandia National Laboratories (SNL) into the flow, sediment dynamics and water-quality code, SNL-EFDC. This routine has been verified and validated against three separate sets of experimental data. With SNL-EFDC, water quality and array optimization studies can be carried out to optimize an MHK array in a resource and study its effects on the environment. The present study examines the effect streamwise and spanwise spacing has on the array performance. Various hypothetical MHK array configurations are simulated within a trapezoidal river channel. Results show a non-linear increase in array-power efficiency as turbine spacing is increased in each direction, which matches the trends seen experimentally. While the sediment transport routines were not used in these simulations, the flow acceleration seen around the MHK arrays has the potential to significantly affect the sediment transport characteristics and benthic habitat of a resource. Evaluation Only. Created with Aspose.Pdf.Kit. Copyright 2002-2011 Aspose Pty Ltd Evaluation Only. Created with Aspose.Pdf.Kit. Copyright 2002-2011 Aspose Pty Ltd

Abstract not provided.

Oxygen gas injection has been studied as one method for rapidly generating hydrogen gas from a uranium hydride storage system. Small scale reactors, 2.9 g UH{sub 3}, were used to study the process experimentally. Complimentary numerical simulations were used to better characterize and understand the strongly coupled chemical and thermal transport processes controlling hydrogen gas liberation. The results indicate that UH{sub 3} and O{sub 2} are sufficiently reactive to enable a well designed system to release gram quantities of hydrogen in {approx} 2 seconds over a broad temperature range. The major system-design challenge appears to be heat management. In addition to the oxidation tests, H/D isotope exchange experiments were performed. The rate limiting step in the overall gas-to-particle exchange process was found to be hydrogen diffusion in the {approx}0.5 {mu}m hydride particles. The experiments generated a set of high quality experimental data; from which effective intra-particle diffusion coefficients can be inferred.

A one-dimensional isobaric reactor model is used to simulate hydrogen isotope exchange processes taking place during flow through a powdered palladium bed. This simple model is designed to serve primarily as a platform for the initial development of detailed chemical mechanisms that can then be refined with the aid of more complex reactor descriptions. The one-dimensional model is based on the Sandia in-house code TPLUG, which solves a transient set of governing equations including an overall mass balance for the gas phase, material balances for all of the gas-phase and surface species, and an ideal gas equation of state. An energy equation can also be solved if thermodynamic properties for all of the species involved are known. The code is coupled with the Chemkin package to facilitate the incorporation of arbitrary multistep reaction mechanisms into the simulations. This capability is used here to test and optimize a basic mechanism describing the surface chemistry at or near the interface between the gas phase and a palladium particle. The mechanism includes reversible dissociative adsorptions of the three gas-phase species on the particle surface as well as atomic migrations between the surface and the bulk. The migration steps are more general than those used previously in that they do not require simultaneous movement of two atoms in opposite directions; this makes possible the creation and destruction of bulk vacancies and thus allows the model to account for variations in the bulk stoichiometry with isotopic composition. The optimization code APPSPACK is used to adjust the mass-action rate constants so as to achieve the best possible fit to a given set of experimental data, subject to a set of rigorous thermodynamic constraints. When data for nearly isothermal and isobaric deuterium-to-hydrogen (D {yields} H) and hydrogen-to-deuterium (H {yields} D) exchanges are fitted simultaneously, results for the former are excellent, while those for the latter show pronounced deviations at long times. These discrepancies can be overcome by postulating the presence of a surface poison such as carbon monoxide, but this explanation is highly speculative. When the method is applied to D {yields} H exchanges intentionally poisoned by known amounts of CO, the fitting results are noticeably degraded from those for the nominally CO-free system but are still tolerable. When TPLUG is used to simulate a blowdown-type experiment, which is characterized by large and rapid changes in both pressure and temperature, discrepancies are even more apparent. Thus, it can be concluded that the best use of TPLUG is not in simulating realistic exchange scenarios, but in extracting preliminary estimates for the kinetic parameters from experiments in which variations in temperature and pressure are intentionally minimized.

Abstract not provided.

Abstract not provided.

Cost-effective implementation of microalgae as a solar-to-chemical energy conversion platform requires extensive system optimization; computer modeling can bring this to bear. This work uses modified versions of the U.S. Environmental Protection Agency's (EPA's) Environmental Fluid Dynamics Code (EFDC) in conjunction with the U.S. Army Corp of Engineers water-quality code (CE-QUAL) to simulate hydrodynamics coupled to growth kinetics of algae (Phaeodactylum tricornutum) in open-channel raceways. The model allows the flexibility to manipulate a host of variables associated with raceway-design, algal-growth, water-quality, hydrodynamic, and atmospheric conditions. The model provides realistic results wherein growth rates follow the diurnal fluctuation of solar irradiation and temperature. The greatest benefit that numerical simulation of the flow system offers is the ability to design the raceway before construction, saving considerable cost and time. Moreover, experiment operators can evaluate the impacts of various changes to system conditions (e.g., depth, temperature, flow speeds) without risking the algal biomass under study.

This document describes the sediment transport subroutines and input files for the Sandia National Laboratories Environmental Fluid Dynamics Code (SNL-EFDC). Detailed descriptions of the input files containing data from Sediment Erosion at Depth flume (SEDflume) measurements are provided along with the description of the source code implementing sediment transport. Both the theoretical description of sediment transport employed in SNL-EFDC and the source code are described. This user manual is meant to be used in conjunction with the EFDC manual (Hamrick 1996) because there will be no reference to the hydrodynamics in EFDC. Through this document, the authors aim to provide the necessary information for new users who wish to implement sediment transport in EFDC and obtain a clear understanding of the source code.

Uncertainty in site characterization arises from a lack of data and knowledge about a site and includes uncertainty in the boundary conditions, uncertainty in the characteristics, location, and behavior of major features within an investigation area (e.g., major faults as barriers or conduits), uncertainty in the geologic structure, as well as differences in numerical implementation (e.g., 2-D versus 3-D, finite difference versus finite element, grid resolution, deterministic versus stochastic, etc.). Since the true condition at a site can never be known, selection of the best conceptual model is very difficult. In addition, limiting the understanding to a single conceptualization too early in the process, or before data can support that conceptualization, may lead to confidence in a characterization that is unwarranted as well as to data collection efforts and field investigations that are misdirected and/or redundant. Using a series of numerical modeling experiments, this project examined the application and use of information criteria within the site characterization process. The numerical experiments are based on models of varying complexity that were developed to represent one of two synthetically developed groundwater sites; (1) a fully hypothetical site that represented a complex, multi-layer, multi-faulted site, and (2) a site that was based on the Horonobe site in northern Japan. Each of the synthetic sites were modeled in detail to provide increasingly informative 'field' data over successive iterations to the representing numerical models. The representing numerical models were calibrated to the synthetic site data and then ranked and compared using several different information criteria approaches. Results show, that for the early phases of site characterization, low-parameterized models ranked highest while more complex models generally ranked lowest. In addition, predictive capabilities were also better with the low-parameterized models. For the latter iterations, when more data were available, the information criteria rankings tended to converge on the higher parameterized models. Analysis of the numerical experiments suggest that information criteria rankings can be extremely useful for site characterization, but only when the rankings are placed in context and when the contribution of each bias term is understood.

Abstract not provided.

The fate of contaminants after a dispersal event is a major concern, and waterways may be particularly sensitive to such an incident. Contaminants could be introduced directly into a water system (municipal or general) or indirectly (Radiological Dispersal Device) from aerial dispersion, precipitation, or improper clean-up techniques that may wash contamination into storm water drains, sewer systems, rivers, lakes, and reservoirs. Most radiological, chemical, and biological contaminants have an affinity for sediments and organic matter in the water system. If contaminated soils enter waterways, a plume of contaminated sediments could be left behind, subject to remobilization during the next storm event. Or, contaminants could remain in place, thus damaging local ecosystems. Suitable planning and deployment of resources to manage such a scenario could considerably mitigate the severity of the event. First responses must be prearranged so that clean-up efforts do not increase dispersal and exacerbate the problem. Interactions between the sediment, contaminant, and water cycle are exceedingly complex and poorly understood. This research focused on the development of a risk-based model that predicts the fate of introduced contaminants in surface water systems. Achieving this goal requires integrating sediment transport with contaminant chemical reactions (sorption and desorption) and surface water hydrodynamics. Sandia leveraged its existing state-of-the-art capabilities in sediment transport measurement techniques, hydrochemistry, high performance computing, and performance assessment modeling in an effort to accomplish this task. In addition, the basis for the physical hydrodynamics is calculated with the EPA sponsored, public domain model, Environmental Fluid Dynamics Code (EFDC). The results of this effort will enable systems analysis and numerical simulation that allow the user to determine both short term and long-term consequences of contamination of waterways as well as to help formulate preventative and remedial strategies.

Given pre-existing Groundwater Modeling System (GMS) models of the Horonobe Underground Research Laboratory (URL) at both the regional and site scales, this work performs an example uncertainty analysis for performance assessment (PA) applications. After a general overview of uncertainty and sensitivity analysis techniques, the existing GMS site-scale model is converted to a PA model of the steady-state conditions expected after URL closure. This is done to examine the impact of uncertainty in site-specific data in conjunction with conceptual model uncertainty regarding the location of the Oomagari Fault. A heterogeneous stochastic model is developed and corresponding flow fields and particle tracks are calculated. In addition, a quantitative analysis of the ratio of dispersive to advective forces, the F-ratio, is performed for stochastic realizations of each conceptual model. Finally, a one-dimensional transport abstraction is modeled based on the particle path lengths and the materials through which each particle passes to yield breakthrough curves at the model boundary. All analyses indicate that accurate characterization of the Oomagari Fault with respect to both location and hydraulic conductivity is critical to PA calculations. This work defines and outlines typical uncertainty and sensitivity analysis procedures and demonstrates them with example PA calculations relevant to the Horonobe URL. Acknowledgement: This project was funded by Japan Nuclear Cycle Development Institute (JNC). This work was conducted jointly between Sandia National Laboratories (SNL) and JNC under a joint JNC/U.S. Department of Energy (DOE) work agreement. Performance assessment calculations were conducted and analyzed at SNL based on a preliminary model by Kashima, Quintessa, and JNC and include significant input from JNC to make sure the results are relevant for the Japanese nuclear waste program.

Given pre-existing Groundwater Modeling System (GMS) models of the Horonobe Underground Research Laboratory (URL) at both the regional and site scales, this work performs an example uncertainty analysis for performance assessment (PA) applications. After a general overview of uncertainty and sensitivity analysis techniques, the existing GMS sitescale model is converted to a PA model of the steady-state conditions expected after URL closure. This is done to examine the impact of uncertainty in site-specific data in conjunction with conceptual model uncertainty regarding the location of the Oomagari Fault. In addition, a quantitative analysis of the ratio of dispersive to advective forces, the F-ratio, is performed for stochastic realizations of each conceptual model. All analyses indicate that accurate characterization of the Oomagari Fault with respect to both location and hydraulic conductivity is critical to PA calculations. This work defines and outlines typical uncertainty and sensitivity analysis procedures and demonstrates them with example PA calculations relevant to the Horonobe URL.

Incorporating results from a previously developed finite element model, an uncertainty and parameter sensitivity analysis was conducted using preliminary site-specific data from Horonobe, Japan (data available from five boreholes as of 2003). Latin Hypercube Sampling was used to draw random parameter values from the site-specific measured, or approximated, physicochemical uncertainty distributions. Using pathlengths and groundwater velocities extracted from the three-dimensional, finite element flow and particle tracking model, breakthrough curves for multiple realizations were calculated with the semi-analytical, one-dimensional, multirate transport code, STAMMT-L. A stepwise linear regression analysis using the 5, 50, and 95% breakthrough times as the dependent variables and LHS sampled site physicochemical parameters as the independent variables was used to perform a sensitivity analysis. Results indicate that the distribution coefficients and hydraulic conductivities are the parameters responsible for most of the variation among simulated breakthrough times. This suggests that researchers and data collectors at the Horonobe site should focus on accurately assessing these parameters and quantifying their uncertainty. Because the Horonobe Underground Research Laboratory is in an early phase of its development, this work should be considered as a first step toward an integration of uncertainty and sensitivity analyses with decision analysis.

Abstract not provided.