Publications

Recently, significant progress has been made in using finite-difference analysis cod es to simulate the responses of complex structures due to direct lightning. Advances have been made in interfacing a finite-difference code with commercial computer aided design tools, in suppressing a weak instability associated with the thin-wire algorithm for modeling conductors much smaller than a cell size, and in visualizing the results with color movies. Preliminary comparisons between the results of the finite-difference code and the results obtained during a recent rocket-triggered lightning test are also presented.

Sandia National Laboratories and ICI Explosives USA have worked together since 1987 to develop computer modeling techniques for Rock Blasting. A result of this effort is the computer program DMC (Distinct Motion Code) which was developed for two-dimensional simulation of rock motion following a blast (Taylor and Preece, 1989 1992). This program has been used to study blasting-induced rock motion resulting from oil shale mining and has been coupled with a gas flow computation capability for better treatment of the explosive behavior. This past year it has been customized for simulations of bench blasting in coat mines and rock quarries (Preece and Knudsen, 1992b). The explicit descretized nature of DMC gives it an advantage over previous blast modeling programs because subtle differences, such as row delay timing, have an influence on the results. This paper will present a DMC study of the influence on percent cast of row delay timing in a typical coal mine bench blast.

The Yucca Mountain Site Characterization Project is studying Yucca Mountain in southwestern Nevada as a potential site for a high-level nuclear waste repository. Site characterization includes surface-based and underground testing. Analyses have been performed to design site characterization activities with minimal impact on the ability of the site to isolate waste, and on tests performed as part of the characterization process. One activity of site characterization is the construction of an Exploratory Studies Facility, consisting of underground shafts, drifts, and ramps, and the accompanying surface pad facility and roads. The information in this report addresses the following topics: (1) a discussion of the potential effects of surface construction water on repository-performance, and on surface and underground experiments; (2) one-dimensional numerical calculations predicting the maximum allowable amount of water that may infiltrate the surface of the mountain without affecting repository performance; and (3) two-dimensional numerical calculations of the movement of that amount of surface water and how the water may affect repository performance and experiments. The results contained herein should be used with other site data and scientific/engineering judgement in determining controls on water usage at Yucca Mountain. This document contains information that has been used in preparing Appendix I of the Exploratory Studies Facility Design Requirements document for the Yucca Mountain Site Characterization Project.

The focus of this paper is on changes in perceptions of the risks associated with nuclear waste management over time. In particular, we are interested in the kinds of change that take place when the management programs, and those who are charged with implementing them, are subject to intensive public debate over an extended period of time. We are undertaken an over-time study of perceived risks in Colorado and New Mexico by implementing sequential random household surveys in each state, timed at six month intervals. This study employs three of these surveys, spanning the period from summer, 1990 to summer, 1991. Using these data, we examine the dynamics that may underlie variations in perceived risks over time. In particular, our analysis is focused on changes in the roles played by (1) basic political orientations (i.e. political ideology) and (2) trust in those who advocate conflicting policy positions.

Midway Valley, located at the eastern base of the Yucca Mountain in southwestern Nevada, is the preferred location of the surface facilities for the potential high-level nuclear waste repository at Yucca Mountain. One goal in siting these surface facilities is to avoid faults that could produce relative displacements in excess of 5 cm in the foundations of the waste-handling buildings. This study reviews existing geologic and geophysical data that can be used to assess the potential for surface fault rupture within Midway Valley. Dominant tectonic features in Midway Valley are north-trending, westward-dipping normal faults along the margins of the valley: the Bow Ridge fault to the west and the Paintbrush Canyon fault to the east. Published estimates of average Quaternary slip rates for these faults are very low but the age of most recent displacement and the amount of displacement per event are largely unknown. Surface mapping and interpretive cross sections, based on limited drillhole and geophysical data, suggest that additional normal faults, including the postulated Midway Valley fault, may exist beneath the Quaternary/Tertiary fill within the valley. Existing data, however, are inadequate to determine the location, recency, and geometry of this faulting. To confidently assess the potential for significant Quaternary faulting in Midway Valley, additional data are needed that define the stratigraphy and structure of the strata beneath the valley, characterize the Quaternary soils and surfaces, and establish the age of faulting. The use of new and improved geophysical techniques, combined with a drilling program, offers the greatest potential for resolving subsurface structure in the valley. Mapping of surficial geologic units and logging of soil pits and trenches within these units must be completed, using accepted state-of-the-art practices supported by multiple quantitative numerical and relative age-dating techniques.

Sandia National Laboratories (SNL) is a Department of Energy multiprogram engineering and scientific facility with unique design, development, and test capabilities arising from their work in nuclear weapons, energy resources, defense systems, nuclear safeguards, and specialized scientific endeavors. To support these programs, they have developed instrumentation and telemetry expertise not available elsewhere. This technology is applicable to projects in government and industry. Since the 1950s, they have applied our technical competence to meet difficult challenges with innovative solutions to data acquisition and telemetry problems. Sandia - with experience in fields as diverse as parachute design and plasma physics, geology and rocket guidance, human factors and high-speed aerodynamics, non-destructive testing and satellite communications - can use the power of synergism among our many disciplines to solve your complex problems of data and acquisition and analysis. SNL solves difficult data acquisition problems for extreme environments with expertise in advanced telemetry techniques, high data rate telemetry design, specialized electronics packaging, MIL-STD-1553 communications, instrumentation development, real-time data analysis, project management, specialized testers and data encryption.

Pretest analysis of a heated block test, proposed for the Exploratory Studies Facility at Yucca Mountain, Nevada, was conducted in this investigation. Specifically, the study focuses on the evaluation of the various designs to drill holes and cut slots for the block. The thermal/mechanical analysis was based on the finite element method and a compliant-joint rock-mass constitutive model. Based on the calculated results, relative merits of the various test designs are discussed.

Photovoltaic energy systems have historically been treated as a bulk power generation source for the future. However, utilities and other agencies involved with electrification throughout the world are beginning to find photovoltaics a least-cost option to meet specific loads both for themselves and their customers, in both off-grid and grid-connected applications. These expanding markets offer the potential of hundreds of megawatts of sales in the coming decade, but a strategy addressing both industrial growth and user acceptance is necessary to capitalize on this opportunity. 11 refs.

Phase mixing of transverse oscillations changes the nature of the ion hose instability from an absolute to a convective instability. The stronger the phase mixing, the faster an electron beam reaches equilibrium with the guiding ion channel. This is important for long distance propagation of relativistic electron beams where it is desired that transverse oscillations phase mix within a few betatron wavelengths of injection and subsequently an equilibrium is reached with no further beam emittance growth. In the linear regime phase mixing is well understood and results in asymptotic decay of transverse oscillations as 1/Z{sup 2} for a Gaussian beam and channel system, Z being the axial distance measured in betatron wavelengths. In the nonlinear regime (which is likely mode of propagation for long pulse beams) results of the spread mass model indicate that phase mixing is considerably weaker than in the regime. In this paper we consider this problem of phase mixing in the nonlinear regime. Results of the spread mass model will be shown along with a simple analysis of phase mixing for multiple oscillator models. Particle simulations also indicate that phase mixing is weaker in nonlinear regime than in the linear regime. These results will also be shown. 3 refs., 4 figs.

Photovoltaic (PV) systems are increasing in popularity in the northern latitudes and in the arctic regions in the state of Alaska. This increased interest and the high cost of providing electric power in these remote areas have prompted the Alaska Energy Authority (AEA) to request assistance from the Photovoltaic Design Assistance Center at Sandia National Laboratories. A project to investigate the feasibility of using PV-Diesel hybrid power systems in small villages in Alaska was started in 1989. Data acquisition systems (DAS) were designed and installed in selected villages to obtain resource and load information. The DAS is described and village electrical and resource data are presented. Simulations were run using the collected village data and actual cost data provided by the AEA. Results of the simulations and the economic analysis are presented. 5 refs., 8 figs.

Fabrication of high-efficiency silicon solar cells in an industrial environment requires a different optimization than in a laboratory environment. Strategies are presented for process development of high-efficiency silicon solar cells, with a goal of simplifying technology transfer into an industrial setting. The strategies emphasize the use of statistical experimental design for process optimization, and the use of baseline processes and cells for process monitoring and quality control. 8 refs.

The purpose of this paper is to develop an analytical model to convert ballistic limit curves obtained from flat projectile experiments to ballistic limit curves based on equivalent diameter spheres. Results from a test program involving flat plat projectiles conducted at Sandia National Laboratories are compared against the predicted performance of equivalent spherical projectiles as determined from the Wilkinson and Cour-Palais penetration equations. The developed method demonstrates good correlation of the ballistic limit of the shield concept for the flat plate projectiles to the theoretical ballistic limit for equivalent spheres as predicted by the penetration equations. 3 refs.

PRA studies are being extended to include a wider spectrum of reactor plants than was considered in NUREG-1150. There is a need for computationally simple models for Direct Containment Heating (DCH) that could be used for screening studies aimed at identifying potentially significant contributors to risk. This paper discusses two adiabatic equilibrium models that are candidates for the task. The first, a 1-cell model, places a true upper bound on DCH loads. This upper bound, however, often far exceeds reasonable expectations of containment loads based on best estimate CONTAIN calculations or experiment observations. In this paper, a 2-cell model is developed that largely captures the major mitigating features of containment compartmentalization, thus providing more reasonable estimates of the containment load. Predictions of the equilibrium models are compared with experiment data from the Limited Flight Path (LFP) test series conducted at Sandia National Laboratories.

The CONTAIN code is currently being used to predict containment thermal hydraulic conditions during design basis and severe accidents for advanced light water reactor (ALWR) designs such as the Westinghouse AP600. In the AP600 design, a passive containment cooling system (PCCS) is used for reducing long-term overpressure during accidents. CONTAIN models for heat and mass transfer within the AP600 containment and outer air cooling channel are verified by comparing recent CONTAIN calculations to integral test data obtained by Westinghouse in their PCCS Integral Test Facility. The comparison includes test in which the outer containment wall is both dry and wet, that is, the wet tests involve an evaporative water film that enhances heat transfer as will be the case for AP600. The appropriateness of the heat and mass transfer analogy methodology used in the CONTAIN code is demonstrated. Code model limitations are discussed along with model development plans and applications for AP600.

An additive three step process has been developed for patterned deposition of Cu onto poly(tetrafluoroethylene) (PTFE). The first step involves patterned irradiation with X-rays or electrons which is thought to cross link the PTFE surface; step two involves chemical etching with the result that only the non-irradiated areas are etched; and step three involves selective chemical vapor deposition (CVD) of Cu onto the etched surface at 200 C using (hexafluoroacetylacetonato)Cu(I) trimethylphosphine ((hfac)Cu(PMe{sub 3})). The non-irradiated areas of the surface are covered by a continuous, dense Cu film with X-ray photoelectron spectra show to contain only surface impurities that are easily removed by a short Ar ion sputter. The irradiated areas show the presence of only C and F, characteristic of PTFE.

The first experiment of the Integral Effects Test (IET-1) series was conducted to investigate the effects of high pressure melt ejection (HPME) on direct containment heating (DCH). A 1:10 linear scale model of the Zion reactor pressure vessel (RPV), cavity, instrument tunnel, and subcompartment structures were constructed in the Surtsey Test Facility at Sandia National Laboratories (SNL). The RPV was modelled with a melt generator that consisted of a steel pressure barrier, a cast MgO crucible, and a thin steel inner liner. The melt generator/crucible had a semi-hemispherical bottom head containing a graphite limitor plate with a 3.5 cm exit hole to simulate the ablated hole in the RPV bottom head that would be formed by tube ejection in a severe nuclear power plant (NPP) accident. The reactor cavity model contained 3.48 kg of water with a depth of 0.9 cm that corresponded to condensate levels in the Zion plant. A steam driven iron oxide/aluminum/chromium thermite was used to simulate HPME. A relatively small steam explosion occurred in the cavity during IET-1. Steam blowthrough entrained debris into the Surtsey vessel resulting in a peak pressure increase in Surtsey of 98 kPa. The Surtsey vessel had been previously inerted with N{sub 2}. The total debris mass ejected into the Surtsey vessel was 43 kg. The hydrogen concentration was 3.1 mol.% in the vessel at equilibrium. The concentration measured inside the subcompartment structures immediately following HPME transient was 20.7 mol.% H{sub 2}. 4 refs., 17 figs., 5 tabs.

During the RADLAC II open-air beam propagation experiments this last year three separate optical diagnostics were used. (1) Streak cameras were deployed to measure matched beam radius and centroid versus time. (2) Three gated, intensified TV cameras were used to image the beam from the end of the propagation range. They gave beam radius and centroid for three slices of the beam over a five meter propagation length. (3) Open shutter cameras were developed to give the time-averaged beam position over the entire propagation range. Data from all three diagnostics were digitized, stored in files on a computer, and post-processed to give temporally and spatially resolved beam size and position. These diagnostics used beam induced air-fluorescence as the mechanism to provide a prompt signal representative of the beam current density. Previous experiments and analysis have shown that the radiation is prompt with the intensity proportional to the beam current density for high energy, high current electron beams propagating in full density air.

The technical issues brought about by recent federal mandates are reviewed and discussed. Progress made in the elimination of CFCs is briefly reviewed. The problems, implications, and status of pending anti-lead legislation and taxation are discussed at length. Recommendations are made for the enactment of rational, fair, and orderly legislation and taxation.

This paper describes current research and development on a miniaturized sensing system for use during in situ characterization of nuclear waste storage tanks. Sandia is designing this sensing system as a tool for a large robotic arm that is deployed through an access port in the top of a storage tank. While the robot arm scans the sensing package over the waste, a distributed computing system acquires sensor data, correlates the data with the position of the robot, and produces maps of the chemical and radiological contents of the tanks in real time. We have built and demonstrated a first prototype system containing eight sensors. 53 refs.

Two revisions of the CONTAIN code, CONTAIN 1.11 and 1.12 , have recently been released. The purpose of this paper is to highlight the new features of these revisions and to discuss other new code features currently under development. The features of CONTAIN 1.11 discussed here include a quasi-mechanistic concrete outgassing model, the connected structure option for heat conduction between compartments, and a new approach for modeling forced convective heat transfer. The direct containment heating (DCH) models released as part of CONTAIN 1.12 are also discussed. New code features currently under development include a revised gas combustion model and a new multifield DCH model. New features of the revised combustion model include the treatment of spontaneous recombination and diffusion flames. CONTAIN plant calculations comparing the old and the revised combustion models are presented. The new features of the multifield DCH model are discussed, and demonstration calculations using this model to analyze a small scale experiment are presented.

The Department of Energy`s (DOE`s) Office of Civilian Radioactive Waste Management (OCRWM) is in the process of developing a new generation of casks to transport spent fuel from commercial nuclear reactor facilities to federal waste facilities. In evaluating the needs of the cask development program a number of unresolved technical issues with potential impacts on the transportation system were identified. This paper provides three samples of issues being addressed by the Cask Systems Development Program for technical resolution: (1) burn-up credit, (2) containment source term evaluation, and (3) weeping.

SMILE is a coaxial Self Magnetically Insulated Transmission Line voltage adder. It replaces the original beam line of the RADLAC II accelerator by a 12.5 m long cathode electrode. The anode electrode remains practically the same, consisting of the original eight insulating stacks or feeds which are connected with equal diameter stainless steel cylinders. The beam is produced at the end of the accelerator and is free of all the possible instabilities associated with accelerating gaps and magnetic vacuum transport. Annular beams with {beta}{perpendicular} {le} 0.1 and radius r{sub b} {le} 1 cm were routinely obtained and extracted from a small magnetically-immersed foilless electron diode. Results of the experimental evaluation are presented and compared with design parameters and numerical simulation predictions. 4 refs.

Abstract not provided.

In this paper, measurements on the quasi-isentropic compression of tungsten to stress levels of 250 GPa are reported. Results of these experiments have been compared to those obtained under shock loading conditions to comparable stresses. These experiments have allowed the determination of temperature, pressure, and loading rate effects on the dynamic yield strength of tungsten up to 250 GPa. These results show that the dynamic yield strength of tungsten is dependent on the loading rate with the strength being higher for the relatively slower rates of loading along the quasi-isentropic. The pressure dependence of the yield strength of tungsten is determined nearly independent of temperature effects from quasi-isentropic loading experiments to 250 GPa, because the temperature rise in an quasi-loading experiment is much lower than those associated with shock loading experiments.

Continuum dynamics codes are categorized as Lagrangian or Eulerian according to the motion of the mesh. A Lagrangian code`s mesh moves with the material, so no mass flows between cells. An Eulerian code`s mesh is stationary, so mass flows between the cells. Eulerian codes have improved to the point where they are routinely used to solve a broad variety of large deformation solid and fluid dynamics problems ranging from air flow over an airplane wing to meteor impact on space structures. This presentation will concentrate on multi-fluid Eulerian codes capable of modeling transient were propagation in solids. These codes use a two-step process to integrate the physics across a time step. The first step, referred to as the Lagrangian step, integrates the physics on a Lagrangian mesh across the time step. The field values are then at the new time, but they are on the distorted Lagrangian mesh. The second step, referred to as the remap step, remaps the data on the distorted Lagrangian mesh back to the original Eulerian mesh thus completing one time step. The algorithms used in the first step are similar to those used in modern Lagrangian codes but they must be extended to handle multi-material cells. The algorithms used in the second step are complex and must be very carefully chosen to minimize errors. These algorithms include second-order, monotone advection equations to calculate the quantities flowing between cells. They also require algorithms that construct material interfaces inside multi-material cells. The strength and limitations of currently used numerical techniques will be discussed. New code development activities that combine the best features on both Lagrangian and Elueian codes will also be discussed. These new codes will employ the strengths of both technologies to address problems that cannot be adequately solved at this time.

The Hybrid Thin-Slot Algorithm (HTSA) integrates a transient integral-equation solution for an aperture in an infinite plane into a finite-difference time-domain (FDTD) code. The technique was introduced for linear apertures and was extended to include wall loss and lossy internal gaskets. A general implementation for arbitrary thin slots is briefly described here. The 3-D FDTD-code TSAR was selected for the implementation. The HTSA does not provide universal solutions to the narrow slot problem, but has merits appropriate for particular applications. The HTSA is restricted to planar slots, but can solve the important case that both the width and depth of the slot are narrow compared to the FDTD spatial cell. IN addition, the HTSA is not bound to the FDTD discrete spatial and time increments, and therefore, high-resolution solutions for the slot physics are possible. The implementation of the HTSA into TSAR is based upon a ``slot data file`` that includes the cell indices where the desired slots are exist within the FDTD mesh. For an HTSA-defined slot, the wall region local to the slot is shorted, and therefore, to change the slot`s topology simply requires altering the file to include the desired cells. 7 refs.

The one-electron energy levels of icosahedral boron clusters have been calculated as a function of intericosahedral spacing maintaining the intraicosahedral spacing of {alpha}-boron. For crystalline lattice constants greater than 1.25 times the equilibrium one, band overlap occurs with concomitant metallic behavior. At smaller lattice constants, orbitals(bands) associated with bonds to adjacent icosahedra are lowered and orbitals(bands) associated with ``antibonds`` are raised. Four bands which were three quarters full become empty, while three bands which were empty become filled. This leads to an energy gap between the filled states and the empty states which accounts for the experimentally observed insulating behavior of this elemental material with three valence electrons per atom.

Verifying the velocity accuracy of a GPS receiver or an integrated GPS/INS system in a dynamic environment is a difficult proposition when many of the commonly used reference systems have velocity uncertainities of the same order of magnitude or greater than the GPS system. The results of flight tests aboard an aircraft in which multiple reference systems simultaneously collected data to evaluate the accuracy of an integrated GPS/INS system are reported. Emphasis is placed on obtaining high accuracy estimates of the velocity error of the integrated system in order to verify that velocity accuracy is maintained during both linear and circular trajectories. Three different reference systems operating in parallel during flight tests are used to independently determine the position and velocity of an aircraft in flight. They are a transponder/interrogator ranging system, a laser tracker, and GPS carrier phase processing. Results obtained from these reference systems are compared against each other and against an integrated real time differential based GPS/INS system to arrive at a set of conclusions about the accuracy of the integrated system.

Computer-aided molecular design methods were used to tailor binding sites for small substrate molecules, including CO{sub 2} and methane. The goal is to design a cavity, adjacent to a catalytic metal center, into which the substrate will selectively bind through only non-bonding interactions with the groups lining the binding pocket. Porphyrins are used as a basic molecular structure, with various substituents added to construct the binding pocket. The conformations of these highly-substituted porphyrins are predicted using molecular mechanics calculations with a force field that gives accurate predictions for metalloporhyrins. Dynamics and energy-minimization calculations of substrate molecules bound to the cavity indicate high substrate binding affinity. The size, shape and charge-distribution of groups surrounding the cavity provide molecular selectivity. Specifically, calculated binding energies of methane, benzene, dichloromethane, CO{sub 2} and chloroform vary by about 10 kcal/mol for metal octaethyl-tetraphenylporphyrins (OETPPs) with chloroform, dichloromethane, and CO{sub 2} having the lowest. Significantly, a solvent molecule is found in the cavity in the X-ray structures of Co- and CuOETPP crystals obtained from dichloromethane. 5 refs., 3 figs., 3 tabs.

Single Event Upset Imaging utilizes the scanning of a micro-focused MeV ion beams across an integrated circuit to test the upset response of the circuit to energetic heavy ions. Using this technique, the position dependence of logic state upsets, as well as the charge collection efficiency of an integrated circuit, can be directly measured with micron resolution. We present in this paper a review of a series of measurements carried out on the TA670 16K static random access memory chip which display this technique`s capabilities.

Short communication.

Hot cracking, or solidification cracking, is one of the most extensively studied phenomenon in welding metallurgy. The efforts made to identify the causes of this type of cracking have been driven by the negative commercial and engineering consequences resulting from the formation of these defects. Through judicious weld joint design, the mechanical restraint can be minimized, but it can never be entirely eliminated simply because metals expand and contract when heated and cooled, respectively. The localized nature of heat input in fusion welding insures a non-homogeneous thermal field being applied to the parts being welded, resulting in the development of strains in the as-solidifying weld metal. With the inherent limitations on the mechanical restraint factor, much research has gone into identifying those alloy compositions which minimize the microstructural factor required for hot cracking to occur. Examples from the author`s own research are presented as a tutorial to show how differential thermal analysis techniques have been used to study the chemical/microstructural factors associated with solidification and fusion zone hot cracking in nickel based engineering alloys. References to other uses of these techniques in related welding metallurgy studies are also given.

Dislocation formation in InAs{sub 1-x}Sb{sub x} buffer layers grown by metal-organic chemical vapor deposition is shown to be reproducibly enhanced by p-type doping at levels greater than or equal to the intrinsic carrier concentration at the growth temperature. To achieve a carrier concentration greater than 2 {times} 10{sup 18} cm{sup {minus}3}, the intrinsic carrier concentration of InSb at 475 C, p-type doping with diethylzinc was used. Carrier concentrations up to 6 {times} 10{sup 18} cm{sup {minus}3} were obtained. The zinc doped buffer layers have proven to be reproducibly crack free for InAs{sub 1-x}Sb{sub x} step graded buffer layers with a final composition of x = 0.12 and a strained layer superlattice with an average composition of x = 0.09. These buffer layers have been used to prepare SLS infrared photodiodes. The details of the buffer layer growth, an explanation for the observed Fermi level effect and the growth and characterization of an infrared photodiode are discussed.

The use of coatings on carbon-carbon materials to reduce the oxidation of carbon is of interest for the production of non-ablative aerospace structures. The arc-jet ground test facility can produce the high energy oxidizing environment necessary to simulate hypersonic flight in which to test candidate coatings. The test conditions usually required are characterized by material temperature and length of time the material remains at that temperature. Material specimens were exposed to high energy supersonic air exhausting from the NASA-Ames Research Center`s 20-MW arc-jet facility. The carbon-carbon materials were heated to required temperatures with arc-heated air for specified lengths of time. This report describes the test methodology and observations of those tests.

We described a new family of versatile, cascadable, optical switches with different functional characteristics -- latching, non-latching, and bistable -- using a single epitaxial structure base don the monolithic integration of photothyristors and surface-emitting layers. High performance optical switching characteristics have been achieved for all three switch archetypes. We also demonstrate the AND, OR, NAND, NOR, and INVERT optical logic functions using monolithic switch structures. 7 refs.

A boundary integral equation method for steady unsaturated flow in nonhomogeneous porous media is presented. Steady unsaturated flow in porous media is described by the steady form of the so-called Richards equation, a highly nonlinear Fokker-Planck equation. By applying a Kirchhoff transformation and employing an exponential model for the relation between capillary pressure and hydraulic conductivity, the flow equation is rendered linear in each subdomain of a piece-wise homogeneous material. Unfortunately, the transformation results in nonlinear conditions along material interfaces, giving rise to a jump in the potential along these boundaries. An algorithm developed to solve the nonhomogeneous flow problem is described and verified by comparison to analytical and numerical solutions. The code is applied to examine the moisture distribution in a layered porous medium due to infiltration from a strip source, a model for infiltration from shallow ponds and washes in arid regions.

Structural system identification is undergoing a period of renewed interest. Probabilistic approaches to physical parameter identification in analysis finite element models make uncertainty in test results an important issue. In this paper, we investigate this issue with a simple, though in many ways representative, structural system. The results of two modal parameter identification techniques are compared and uncertainty estimates, both through bias and random errors, are quantified. The importance of the interaction between test and analysis is also highlighted. 25 refs.

Pre-exposure induced stress corrosion cracking (SCC) of an Al-Li-Cu, AA 2090, was studied using a variety of test techniques. Results from SCC testing in a simulated isolated pit solution are correlated with electrochemical corrosion rate data obtained for individual phases in the subgrain boundary region. These experimental data, combined with existing data on the crevice chemistry of isolated pits in Al-Li alloys and X-ray diffraction studies of solid corrosion products formed in crevice environments are used to propose a model for pre-exposure induced cracking based on anodic dissolution along subgrain boundaries. Key features of the model are selective dissolution of the subgrain boundary T{sub 1} phase (Al{sub 2}CuLi) at the crack tip and passivation of crack walls by the formation of an Li{sub 2}[Al{sub 2}(OH){sub 6}]{sub 2}{center_dot}CO{sub 3}{center_dot}nH{sub 2}O barrier film.

The transportation risk evaluation code RADTRAN 4 is designed to evaluate doses and risks associated with the transportation of radioactive materials (Ne92). RADTRAN 4 may be used to calculate dose consequences for incident-free transportation and dose risks for accidents. Consequences of normal (or incident-free) transportation include doses to crew members, persons at stops, and members of the public sharing a route segment (on-link) and residing near the segment (off-link) during normal transportation. These dose estimates are not multiplied by a probability factor and, hence, are referred to as dose consequences. Calculated doses that might be incurred during accidents are multiplied by the probabilities of those accidents, and hence are referred to as dose risks. RADTRAN 4 includes a LINK option that allows the user to characterize each link or segment of a transportation route in greater detail than that provided by average or default values for route-related parameters.

This paper explains how an induction coilgun works and presents the factors which go into its design. Our purpose is to obtain algebraic expressions which, although crude, provide useful predictors of behavior, illustrate the dependence on various parameters, and suggest ways to optimize the design. Detailed prediction of the gun`s behavior can be obtained from simulation codes, such as SLINGSHOT.

Very high driving pressures (tens or hundreds of GPa), are required to accelerate flier plats to hypervelocities. This loading pressure pulse on the fiber plates must be nearly shockless to prevent the plate from melting or vaporizing. This is accomplished by using graded-density impactors referred to as ``pillows.`` When this graded-density material is used to impact a flier-plate in a modified two-stage light gas gun, nearly shockless megabar pressures are introduced into the flier plate. The pressure pulses must also be tailored to prevent spallation of the flier-plate. This technique has been used to launch nominally 1-mm-thick aluminum, magnesium and titanium (gram-size) intact plates to 10.4 km/s, and 0.5-mm-thick aluminum and titanium (half-gram size) intact plates to 12.2 km/s. This is the highest mass-velocity capability attained with laboratory launchers to data, and should open up new regimes of impact physics and lethality studies related to space sciences for laboratory investigations. 14 refs.

Relativistic high current electron beams can be transported long distances across the geomagnetic field using the IFR (Ion focused Regime) technique. IFR is a method of providing strong electrostatic focusing and guiding of the beam. The guiding is sufficiently strong to allow the beam to transport any angle with respect to geomagnetic field. In the IFR method, first an ionizing laser (or any ionizing method) is used to create a preionized cylindrical channel.

Very high pressure and acceleration is necessary to launch flier plates to hypervelocities. In addition, the high pressure loading must be uniform, structured, and shockless, i.e., time-dependent to prevent the flier plate from either fracturing or melting. In this paper, a novel technique is described which allows the use of 100 GPa megabar loading pressures and 10{sup 9}-g acceleration to launch intact flier plates to velocities of 12.2 km/s. The technique has been used to launch nominally 1-mm thick aluminum, magnesium and titanium alloy plates to velocities over 10 km/s, and 0.5-mm thick aluminum and titanium alloy plates to velocities of 12.2 km/s.

We present a learning algorithm designed to improve robot path planning. The algorithm relies on an existing path planner to provide solutions to difficult tasks. From these solutions, it learns a sparse network of useful robot subgoals which guide and support fast planning. We analyze the algorithm theoretically by developing some general techniques useful in characterizing behaviors of probabilistic learning. We also demonstrate the effectiveness of the algorithm empirically with an existing path planner in practical environments. The learning algorithm not only reduces the time cost of existing planners, but also increases their capability in solving difficult tasks. 7 refs.

An understanding of the state of stress on faults is important for pre- and postclosure performance considerations for the potential high-level radioactive waste repository at Yucca Mountain. This paper presents the results of three-dimensional numerical analyses that provide estimates of the state of stress through time (10,000 years) along three major faults in the vicinity of the potential repository due to thermal stresses resulting from waste emplacement. It was found, that the safety factor for slip close to the potential repository increases with time after waste emplacement. Possible fault slip is predicted above and below the potential repository for certain loading conditions and times. In general, thermal loading reduces the potential for slip in the vicinity of the potential repository.

This paper presents a summary of the conduct and findings of the Exploratory Studies Facility Alternatives Study (ESF-AS). The Exploratory Studies Facility (ESF) is being planned for use in the characterization of a site for a potential high-level nuclear waste repository at Yucca Mountain, NV. The purpose of the ESF-AS were to identify and rank order ESF-repository options and to improve understanding of the favorable or unfavorable features of the ESF design. The analysis resulted in the ranking of 34 options, in accordance with the extent to which each option could achieve the objectives. Additional findings regarding design features that were identified as key elements in an option`s ability to provide good overall performance are also discussed.

The eutectic Sn-Pb solder alloy is discussed with respect to alloy development options to improve the thermomechanical fatigue behavior of solder joints. Eutectic Sn-Pb solder joints fail through the development of a heterogeneous coarsened band of recrystallized and coarsened Pb- and Sn-rich phases. All imposed deformation concentrates solely into this thin region, accelerating fatigue failures. The development of solder alloys is currently being undertaken to improve the fatigue characteristics of eutectic Sn-Pb solder. New alloys must retain wetting and manufacturability characteristics similar to eutectic Sn-Pb. The options discussed to improve fatigue life include: creating a fine superplastic microstructure, small alloy additions to homogenize the microstructure, carbon reinforced composite solder, dispersed second phase precipitates that break up the solder microstructure, and using different solder alloys to replace eutectic Sn-Pb.

The large deformation elastic response of a plane woven Kevlar fabric is investigated analytically and experimentally. The analysis assumes the undeformed geometry to be a sequence of interlaced arcs of circles that reverse at each yarn midpoint, and each yarn is modeled as an extensible elastica subject to certain compatibility conditions. Deflection-force relations for the fabric are determined in terms of the initial weave geometry and the elastic properties of the individual yarns. The theoretical results agree well with the results of experiments performed on a fabric woven from 400 denier Kevlar yarns under conditions of uniaxial loading in both warp and fill directions.

The Faceted Stretched-Membrane Dish Program is part of a DOE-sponsored effort to develop a commercial 25 kWe dish/Stirling system employing a twelve-facet dish concentrator. The facets will utilize the stretched-membrane technology originated in the heliostat development program. Each facet is constructed with a thin metal membrane stretched over both sides of a steel ring. When a small vacuum is induced between the membranes they assume a parabolic contour capable of concentrating sunlight at a predetermined focal length. A reflective polymer film is attached to the face of the facet of the facet to enhance the optical performance. During Phase II of the Faceted Stretched-Membrane Dish Program, Science Applications International Corp. and Solar Kinetics, Inc., constructed prototype 3.5-meter facets utilizing different design approaches to demonstrate their manufacturability and optical performance. Sandia engaged in a program to determine the on-sun performance of the facets (for f/Ds of 2.7 to 3.0). A uniformly distributed slope error was used as the basis for comparison. Flux arrays based on slope error from a computer model were compared to a measured flux array for each facet. The slope error for the facet was determined by the value that would produce a modeled array with the minimum mean square difference to the measured array. The facet produced by SAIC demonstrated uniform slope errors of 2.2 to 3.0 milliradians with peak flux intesities of 334 to 416 kW/m{sup 2}. The SKI facet had slope errors of 1.6 to 1.9 milliradians with peak flux intesities of 543 to 1186 kW/m{sup 2}.

A 1-MeV neutron damage equivalence methodology and damage function have been developed for GaAs based on a recoil-energy dependent damage efficiency and the displacement kerma. This method, developed using life-time degradation in GaAs LEDs in a variety of neutron spectra, is also shown to be applicable to carrier removal. A validated methodology, such as this, is required to ensure and evaluate simulation fidelity in the neutron testing of GaAs semiconductors.