Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Injection of CO2 into formations containing brine is proposed as a long-term sequestration solution. A significant obstacle to sequestration performance is the presence of existing wells providing a transport pathway out of the sequestration formation. To understand how heterogeneity impacts the leakage rate, we employ two dimensional models of the CO2 injection process into a sandstone aquifer with shale inclusions to examine the parameters controlling release through an existing well. This scenario is modeled as a constant-rate injection of super-critical CO2 into the existing formation where buoyancy effects, relative permeabilities, and capillary pressures are employed. Three geologic controls are considered: stratigraphic dip angle, shale inclusion size and shale fraction. In this study, we examine the impact of heterogeneity on the amount and timing of CO2 released through a leaky well. Sensitivity analysis is performed to classify how various geologic controls influence CO2 loss. A 'Design of Experiments' approach is used to identify the most important parameters and combinations of parameters to control CO2 migration while making efficient use of a limited number of computations. Results are used to construct a low-dimensional description of the transport scenario. The goal of this exploration is to develop a small set of parametric descriptors that can be generalized to similar scenarios. Results of this work will allow for estimation of the amount of CO2 that will be lost for a given scenario prior to commencing injection. Additionally, two-dimensional and three-dimensional simulations are compared to quantify the influence that surrounding geologic media has on the CO2 leakage rate.

Abstract not provided.

Abstract not provided.

This report evaluates the feasibility of high-level radioactive waste disposal in shale within the United States. The U.S. has many possible clay/shale/argillite basins with positive attributes for permanent disposal. Similar geologic formations have been extensively studied by international programs with largely positive results, over significant ranges of the most important material characteristics including permeability, rheology, and sorptive potential. This report is enabled by the advanced work of the international community to establish functional and operational requirements for disposal of a range of waste forms in shale media. We develop scoping performance analyses, based on the applicable features, events, and processes identified by international investigators, to support a generic conclusion regarding post-closure safety. Requisite assumptions for these analyses include waste characteristics, disposal concepts, and important properties of the geologic formation. We then apply lessons learned from Sandia experience on the Waste Isolation Pilot Project and the Yucca Mountain Project to develop a disposal strategy should a shale repository be considered as an alternative disposal pathway in the U.S. Disposal of high-level radioactive waste in suitable shale formations is attractive because the material is essentially impermeable and self-sealing, conditions are chemically reducing, and sorption tends to prevent radionuclide transport. Vertically and laterally extensive shale and clay formations exist in multiple locations in the contiguous 48 states. Thermal-hydrologic-mechanical calculations indicate that temperatures near emplaced waste packages can be maintained below boiling and will decay to within a few degrees of the ambient temperature within a few decades (or longer depending on the waste form). Construction effects, ventilation, and the thermal pulse will lead to clay dehydration and deformation, confined to an excavation disturbed zone within a few meters of the repository, that can be reasonably characterized. Within a few centuries after waste emplacement, overburden pressures will seal fractures, resaturate the dehydrated zones, and provide a repository setting that strongly limits radionuclide movement to diffusive transport. Coupled hydrogeochemical transport calculations indicate maximum extents of radionuclide transport on the order of tens to hundreds of meters, or less, in a million years. Under the conditions modeled, a shale repository could achieve total containment, with no releases to the environment in undisturbed scenarios. The performance analyses described here are based on the assumption that long-term standards for disposal in clay/shale would be identical in the key aspects, to those prescribed for existing repository programs such as Yucca Mountain. This generic repository evaluation for shale is the first developed in the United States. Previous repository considerations have emphasized salt formations and volcanic rock formations. Much of the experience gained from U.S. repository development, such as seal system design, coupled process simulation, and application of performance assessment methodology, is applied here to scoping analyses for a shale repository. A contemporary understanding of clay mineralogy and attendant chemical environments has allowed identification of the appropriate features, events, and processes to be incorporated into the analysis. Advanced multi-physics modeling provides key support for understanding the effects from coupled processes. The results of the assessment show that shale formations provide a technically advanced, scientifically sound disposal option for the U.S.

This report summarizes the results of a computer model that describes the behavior of pulsating heat pipes (PHP). The purpose of the project was to develop a highly efficient (as compared to the heat transfer capability of solid copper) thermal groundplane (TGP) using silicon carbide (SiC) as the substrate material and water as the working fluid. The objective of this project is to develop a multi-physics model for this complex phenomenon to assist with an understanding of how PHPs operate and to be able to understand how various parameters (geometry, fill ratio, materials, working fluid, etc.) affect its performance. The physical processes describing a PHP are highly coupled. Understanding its operation is further complicated by the non-equilibrium nature of the interplay between evaporation/condensation, bubble growth and collapse or coalescence, and the coupled response of the multiphase fluid dynamics among the different channels. A comprehensive theory of operation and design tools for PHPs is still an unrealized task. In the following we first analyze, in some detail, a simple model that has been proposed to describe PHP behavior. Although it includes fundamental features of a PHP, it also makes some assumptions to keep the model tractable. In an effort to improve on current modeling practice, we constructed a model for a PHP using some unique features available in FLOW-3D, version 9.2-3 (Flow Science, 2007). We believe that this flow modeling software retains more of the salient features of a PHP and thus, provides a closer representation of its behavior.

This document summarizes research and planning for the development of a numerical simulation capability for nonisothermal multiphase, multicomponent transport in heterogeneous deformable porous materials. Particular attention is given to describing a mathematical formulation for flow in deformable media and for numerical techniques for dealing with phase transitions. A development plan is formulated to provide a computational capability motivated by current and future needs in geosystems management for energy security.

Abstract not provided.

This report documents the results for the FY07 ASC Integrated Codes Level 2 Milestone number 2354. The description for this milestone is, 'Demonstrate level set free surface tracking capabilities in ARIA to simulate the dynamics of the formation and time evolution of a weld pool in laser welding applications for neutron generator production'. The specialized boundary conditions and material properties for the laser welding application were implemented and verified by comparison with existing, two-dimensional applications. Analyses of stationary spot welds and traveling line welds were performed and the accuracy of the three-dimensional (3D) level set algorithm is assessed by comparison with 3D moving mesh calculations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Analytical and numerical solutions are presented for steady evaporating flow in open microchannels having a rectangular cross section and a uniform depth. The flow, driven by the axial gradient of capillary pressure, generally consists of an entry region where the meniscus is attached to the top corners of the channel followed by a jump-like transition to a corner-flow region in which the meniscus progressively recedes into the bottom corners of the channel. Illustrative numerical solutions are used to guide the derivation of an easily applied analytical approximation for the maximum sustainable heat flux or capillary limit. © 2005 Elsevier Ltd. All rights reserved.

The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass-lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov-Galerkin streamline upwind schemes in the context of advection-diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection-dominated problems.

Abstract not provided.

The derivation and justification for various low-speed approximations to the fully compressible, Navier-Stokes equations are presented. A numerical formulation based on the finite element method is developed and implemented as an extension to the standard Boussinesq equations. Example steady and transient flow problems are simulated to examine the performance of the numerical algorithm and the solution differences with the more commonly studied Boussinesq approximation. © 2004 Elsevier B.V. All rights reserved.

Abstract not provided.

Abstract not provided.

This report presents a detailed multi-methods comparison of the spatial errors associated with finite difference, finite element and finite volume semi-discretizations of the scalar advection-diffusion equation. The errors are reported in terms of non-dimensional phase and group speeds, discrete diffusivity, artificial diffusivity, and grid-induced anisotropy. It is demonstrated that Fourier analysis (aka von Neumann analysis) provides an automatic process for separating the spectral behavior of the discrete advective operator into its symmetric dissipative and skew-symmetric advective components. Further it is demonstrated that streamline upwind Petrov-Galerkin and its control-volume finite element analogue, streamline upwind control-volume, produce both an artificial diffusivity and an artificial phase speed in addition to the usual semi-discrete artifacts observed in the discrete phase speed, group speed and diffusivity. For each of the numerical methods considered, asymptotic truncation error and resolution estimates are presented for the limiting cases of pure advection and pure diffusion. The Galerkin finite element method and its streamline upwind derivatives are shown to exhibit super-convergent behavior in terms of phase and group speed when a consistent mass matrix is used in the formulation. In contrast, the CVFEM method and its streamline upwind derivatives yield strictly second-order behavior. While this work can only be considered a first step in a comprehensive multi-methods analysis and comparison, it serves to identify some of the relative strengths and weaknesses of multiple numerical methods in a common mathematical framework.