V-TOUGH: An enhanced version of the TOUGH code for the thermal and hydrologic simulation of large-scale problems in nuclear waste isolation

The TOUGH code developed at Lawrence Berkeley Laboratory (LBL) is being extensively used to numerically simulate the thermal and hydrologic environment around nuclear waste packages in the unsaturated zone for the Yucca Mountain Project. At the Lawrence Livermore National Laboratory (LLNL) we have rewritten approximately 80 percent of the TOUGH code to increase its speed and incorporate new options. The geometry of many problems requires large numbers of computational elements elements in order to realistically model detailed physical phenomena, and, as a result, large amounts of computer time are needed. In order to increase the speed of the code we have incorporated fast linear equation solvers, vectorization of substantial portions of code, improved automatic time stepping, and implementation of table look-up for the steam table properties. These enhancements have increased the speed of the code for typical problems by a factor of 20 on the Cray 2 computer. In addition to the increase in computational efficiency we have added several options: vapor pressure lowering; equivalent continuum treatments of fractures; energy and material volumetric, mass and flux accounting; and Stefan-Boltzmann radiative heat transfer. 5 refs

Toxic Phenolic Glycosides From \u3ci\u3ePopulus:\u3c/i\u3e Physiological Adaptations of the Western North American Tiger Swallowtail Butterfly, \u3ci\u3ePapilio Rutulus\u3c/i\u3e (Lepidoptera: Papilionidae)

The phenolic glycosides tremulacin and salicortin found in quaking aspen, Populus tremuloides, and other members of the Salicaceae, are known to be toxic to larvae of the Eastern tiger swallowtail butterfly, Papilio glaucus, but not to the Canadian tiger swallowtail, P. canadensis. Larvae of the western tiger swallowtail, P. rutulus, were not killed nor were their growth rates suppressed when fed a mixture of tremulacin and salicortin on black cherry leaves. When the Salicaceae adapted P. rutulus penultimate instar larvae were fed a combination of the two phenolic glycosides and the esterase inhibitor (DEF S,S,S-tributylphosphorotrithioate), growth was reduced more than 50OJo compared to controls, and half of the larvae died before completing the instaL Our results indicate that esterase detoxification mechanisms are involved in the western tiger swallowtail, P. rutulus, as is also known to be the case for the northern tiger swallowtail, P. canadensis. It is not known whether the same esterase isozyme is involved in both species. From an evolutionary perspective such information could help resolve whether the Salicaceae-adapted swallowtails species are a monophyletic group (perhaps due to isolation in the Beringial Pleistocene glacial refuge of Alaska)

On the movement of a liquid front in an unsaturated, fractured porous medium, Part 2, Mathematical theory

A simplified equation of motion is derived for the flow of liquid through an idealized one-dimensional fracture situated in an unsaturated imbibing porous medium. The equation is valid for the case where the matrix material has a much lower saturated conductivity than that of the fracture and the capillary tension in the matrix is sufficiently stronger than gravity. Asymptotic solutions are given for the motion of the liquid front in a parallel fracture system. With the introduction of natural time constants and dimensionless parameters, the flow behavior can be shown to possess various temporal flow regimes. This work is part of the Nevada Nuclear Waste Storage Project and is applicable to understanding some of the various physical parameters affecting liquid flow through a fracture in an unsaturated porous medium, and is particularly useful as a step in understanding the hydrological processes around a nuclear waste repository in an unsaturated environment as well as in other applications where unsaturated fracture flow conditions exist. The solutions are also relevant to numerical model verification. 10 refs., 2 tabs

Documentation of the Thermal Energy Balance Equation used in the USNT Model of the NUFT Flow and Transport Code

The main purpose of this document is to describe the thermal energy balance used in the USNT module of the NUFT code (Nitao, 1998) and its derivation. We also address several specific issues that have arisen within the Yucca Mountain Project

Numerical modeling of the thermal and hydrological environment around a nuclear waste package using the equivalent continuum approximation: Horizontal emplacement

In support of the investigations for an underground high-level nuclear waste repository at Yucca Mountain, Nevada, we have performed computer simulations of the immediate thermal and hydrological environment around a nuclear waste package. Calculations of this type will be needed for waste package design, performance assessment, and radionuclide transport analyses. Two dimensional computer simulations using a modified version of the TOUGH code were run for an idealized configuration derived from the COVE3 benchmarking effort consisting of a single spent fuel waste package with laterally periodic boundary conditions. The model domain extended downward to the water table and upward to the ground level. Fluid behavior in the rock was modeled using the equivalent continuum approximation. Runs were made with surface water influx rates at the surface set to 0.1, 0.5, and 1.0 mm/yr. A significant amount of code modification and development was needed in order to develop the capability to run these types of problems out to the long time spans required. 26 refs., 59 figs

Toxic Phenolic Glycosides From \u3ci\u3ePopulus:\u3c/i\u3e Physiological Adaptations of the Western North American Tiger Swallowtail Butterfly, \u3ci\u3ePapilio Rutulus\u3c/i\u3e (Lepidoptera: Papilionidae)

The phenolic glycosides tremulacin and salicortin found in quaking aspen, Populus tremuloides, and other members of the Salicaceae, are known to be toxic to larvae of the Eastern tiger swallowtail butterfly, Papilio glaucus, but not to the Canadian tiger swallowtail, P. canadensis. Larvae of the western tiger swallowtail, P. rutulus, were not killed nor were their growth rates suppressed when fed a mixture of tremulacin and salicortin on black cherry leaves. When the Salicaceae adapted P. rutulus penultimate instar larvae were fed a combination of the two phenolic glycosides and the esterase inhibitor (DEF S,S,S-tributylphosphorotrithioate), growth was reduced more than 50OJo compared to controls, and half of the larvae died before completing the instaL Our results indicate that esterase detoxification mechanisms are involved in the western tiger swallowtail, P. rutulus, as is also known to be the case for the northern tiger swallowtail, P. canadensis. It is not known whether the same esterase isozyme is involved in both species. From an evolutionary perspective such information could help resolve whether the Salicaceae-adapted swallowtails species are a monophyletic group (perhaps due to isolation in the Beringial Pleistocene glacial refuge of Alaska)

On the movement of a liquid front in an unsaturated, fractured porous medium, Part 1

The primary aim of this paper is to present approximate analytical solutions of the fracture flow which gives the position of the liquid fracture front as a function of time. These solutions demonstrate that the liquid movement in the fracture can be classified into distinctive time periods, or flow regimes. It is also shown that when plotted versus time using a log-log scale, the liquid fracture front position asymptotically approaches a series of line segments. Two-dimensional numerical simulations were run utilizing input data applicable to the densely welded, fractured tuff found at Yucca Mountain in order to confirm these observations. 19 refs., 15 figs., 8 tabs

Evaluation of thermo-hydrological performance in support of the thermal loading systems study

Heat generated as a result of emplacing spent nuclear fuel will significantly affect the pre- and post-closure performance of the Mined Geological Disposal System (MGDS) at the potential repository site in Yucca Mountain. Understanding thermo-hydrological behavior under repository thermal loads is essential in (a) planning and conducting the site characterization and testing program, (b) designing the repository and engineered barrier system, and (c) assessing performance. The greatest concern for hydrological performance is source of water that would contact a waste package, accelerate its failure rate, and eventually transport radionuclides to the water table. The primary sources of liquid water are: (1) natural infiltration, (2) condensate generated under boiling conditions, and (3) condensate generated under sub-boiling conditions. Buoyant vapor flow, occurring either on a sub-repository scale or on a mountain scale, any affect the generation of the second and third sources of liquid water. A system of connected fractures facilitates repository-heat-driven gas and liquid flow as well as natural infiltration. With the use of repository-scale and sub-repository-scale models, the authors analyze thermo-hydrological behavior for Areal Mass Loadings (AMLs) of 24.2, 35.9, 55.3, 83.4, and 110.5 MTU/acre for a wide range of bulk permeability. They examine the temporal and spatial extent of the temperature and saturation changes during the first 100,000 yr. They also examine the sensitivity of mountain scale moisture redistribution to a range of AMLs and bulk permeabilities. In addition, they investigate how boiling and buoyant, gas-phase convection influence thermo-hydrological behavior in the vicinity of emplacement drifts containing spent nuclear fuel

Double Diffusive Natural Convection in a Nuclear Waste Repository

In this study, we conduct a two-dimensional numerical analysis of double diffusive natural convection in an emplacement drift for a nuclear waste repository. In-drift heat and moisture transport is driven by combined thermal- and compositional-induced buoyancy forces. Numerical results demonstrate buoyancy-driven convective flow patterns and configurations during both repository heat-up and cool-down phases. It is also shown that boundary conditions, particularly on the drip-shield surface, have strong impacts on the in-drift convective flow and transport

A Markov-Chain Monte-Carlo Based Method for Flaw Detection in Beams

A Bayesian inference methodology using a Markov Chain Monte Carlo (MCMC) sampling procedure is presented for estimating the parameters of computational structural models. This methodology combines prior information, measured data, and forward models to produce a posterior distribution for the system parameters of structural models that is most consistent with all available data. The MCMC procedure is based upon a Metropolis-Hastings algorithm that is shown to function effectively with noisy data, incomplete data sets, and mismatched computational nodes/measurement points. A series of numerical test cases based upon a cantilever beam is presented. The results demonstrate that the algorithm is able to estimate model parameters utilizing experimental data for the nodal displacements resulting from specified forces