Publications

This report describes work done in FY2003 under Advanced and Exploratory Studies funding for Advanced Weapons Controllers. The contemporary requirements and envisioned missions for nuclear weapons are changing from the class of missions originally envisioned during development of the current stockpile. Technology available today in electronics, computing, and software provides capabilities not practical or even possible 20 years ago. This exploratory work looks at how Weapon Electrical Systems can be improved to accommodate new missions and new technologies while maintaining or improving existing standards in nuclear safety and reliability.

A concurrent computational and experimental investigation of thermal transport is performed with the goal of improving understanding of, and predictive capability for, thermal transport in microdevices. The computational component involves Monte Carlo simulation of phonon transport. In these simulations, all acoustic modes are included and their properties are drawn from a realistic dispersion relation. Phonon-phonon and phonon-boundary scattering events are treated independently. A new set of phonon-phonon scattering coefficients are proposed that reflect the elimination of assumptions present in earlier analytical work from the simulation. The experimental component involves steady-state measurement of thermal conductivity on silicon films as thin as 340nm at a range of temperatures. Agreement between the experiment and simulation on single-crystal silicon thin films is excellent, Agreement for polycrystalline films is promising, but significant work remains to be done before predictions can be made confidently. Knowledge gained from these efforts was used to construct improved semiclassical models with the goal of representing microscale effects in existing macroscale codes in a computationally efficient manner.

The work discussed in this report was supported by a Campus Fellowship LDRD. The report contains three papers that were published by the fellowship recipient and these papers form the bulk of his dissertation. They are reproduced here to satisfy LDRD reporting requirements.

As MEMS transducers are scaled up in size, the threshold is quickly crossed to where magnetoquasistatic (MQS) transducers are superior for force production compared to electroquasistatic (EQS) transducers. Considerable progress has been made increasing the force output of MEMS EQS transducers, but progress with MEMS MQS transducers has been more modest. A key reason for this has been the difficulty implementing efficient lithographically-fabricated magnetic coil structures. The contribution of this study is a planar multilayer polyphase coil architecture which provides for the lithographic implementation of efficient stator windings suitable for linear magnetic machines. A millimeter-scale linear actuator with complex stator windings was fabricated using this architecture. The stators of the actuator were fabricated using a BCB/Cu process, which does not require replanarization of the wafer between layers. The prototype stator was limited to thin copper layers (3 {micro}m) due to the use of evaporated metal at the time of fabrication. Two layers of metal were implemented in the prototype, but the winding architecture naturally supports additional metal layer pairs. It was found in laboratory tests that the windings can support very high current densities of 4 x 10{sup 9}A/m{sup 2} without damage. Force production normal to the stator was calculated to be 0.54 N/A. For thin stators such as this one, force production increases approximately linearly with the thickness of the windings and a six-layer stator fabricated using a newly implemented electroplated BCB/Cu process (six layers of 15 {micro}m thick metal) is projected to produce approximately 8.8 N/A.

A laser safety hazard evaluation and pertinent output measurements were performed (June 2003 through August 2003) on several VITAL-2 Variable Intensity Tactical Aiming Light--infrared laser, associated with the Proforce M-4 system used in force-on-force exercises. The VITAL-2 contains two diode lasers presenting 'Extended Source' viewing out to a range on the order of 1.3 meters before reverting to a 'Small Source' viewing hazard. Laser hazard evaluation was performed in concert with the ANSI Std. Z136.1-2000 for the safe use of lasers and the ANSI Std. Z136.6-2000 for the safe use of lasers outdoors. The results of the laser hazard analysis for the VITAL-2, indicates that this Tactical Aiming IR laser presents a Class 1 laser hazard to personnel in the area of use. Field measurements performed on 71 units confirmed that the radiant outputs were at all times below the Allowable Emission Limit and that the irradiance of the laser spot was at all locations below the Maximum Exposure Limit. This system is eye-safe and it may be used under current SNL policy in force-on-force exercises. The VITAL-2 Variable Intensity Tactical Aiming Light does not present a laser hazard greater than Class 1, to aided viewing with binoculars.

Alloying element loss from the weld pool during laser spot welding of stainless steel was investigated experimentally and theoretically. The experimental work involved determination of work-piece weight loss and metal vapor composition for various welding conditions. The transient temperature and velocity fields in the weld pool were numerically simulated. The vaporization rates of the alloying elements were modeled using the computed temperature profiles. The fusion zone geometry could be predicted from the transient heat transfer and fluid flow model for various welding conditions. The laser power and the pulse duration were the most important variables in determining the transient temperature profiles. The velocity of the liquid metal in the weld pool increased with time during heating and convection played an increasingly important role in the heat transfer. The peak temperature and velocity increased significantly with laser power density and pulse duration. At very high power densities, the computed temperatures were higher than the boiling point of 304 stainless steel. As a result, evaporation of alloying elements was caused by both the total pressure and the concentration gradients. The calculations showed that the vaporization occurred mainly from a small region under the laser beam where the temperatures were very high. The computed vapor loss was found to be lower than the measured mass loss because of the ejection of tiny metal droplets owing to the recoil force exerted by the metal vapours. The ejection of metal droplets has been predicted by computations and verified by experiments.

Detailed experiments involving extensive high resolution transmission electron microscopy (TEM) revealed significant microstructural differences between Cu sulfides formed at low and high relative humidity (RH). It was known from prior experiments that the sulfide grows linearly with time at low RH up to a sulfide thickness approaching or exceeding one micron, while the sulfide initially grows linearly with time at high RH then becomes sub-linear at a sulfide thickness less than about 0.2 microns, with the sulfidation rate eventually approaching zero. TEM measurements of the Cu2S morphology revealed that the Cu2S formed at low RH has large sized grains (75 to greater than 150 nm) that are columnar in structure with sharp, abrupt grain boundaries. In contrast, the Cu2S formed at high RH has small equiaxed grains of 20 to 50 nm in size. Importantly, the small grains formed at high RH have highly disordered grain boundaries with a high concentration of nano-voids. Two-dimensional diffusion modeling was performed to determine whether the existence of localized source terms at the Cu/Cu2S interface could be responsible for the suppression of Cu sulfidation at long times at high RH. The models indicated that the existence of static localized source terms would not predict the complete suppression of growth that was observed. Instead, the models suggest that the diffusion of Cu through Cu2S becomes restricted during Cu2S formation at high RH. The leading speculation is that the extensive voiding that exists at grain boundaries in this material greatly reduces the flux of Cu between grains, leading to a reduction in the rate of sulfide film formation. These experiments provide an approach for adding microstructural information to Cu sulfidation rate computer models. In addition to the microstructural studies, new micro-patterned test structures were developed in this LDRD to offer insight into the point defect structure of Cu2S and to permit measurement of surface reaction rates during Cu sulfidation. The surface reaction rate was measured by creating micropatterned Cu lines of widths ranging from 5 microns to 100 microns. When sulfidized, the edges of the Cu lines show greater sulfidation than the center, an effect known as microloading. Measurement of the sulfidation profile enables an estimate of the ratio of the diffusivity of H2S in the gas phase to the surface reaction rate constant, k. Our measurements indicated that the gas phase diffusivity exceeds k by more than 10, but less than 100. This is consistent with computer simulations of the sulfidation process. Other electrical test structures were developed to measure the electrical conductivity of Cu2S that forms on Cu. This information can be used to determine relative vacancy concentrations in the Cu2S layer as a function of RH. The test structures involved micropatterned Cu disks and thin films, and the initial measurements showed that the electrical approach is feasible for point defect studies in Cu2S.

Abstract not provided.

Abstract not provided.

In the epitaxial lateral overgrowth of GaN, mass transport and the effects of crystal-growth kinetics lead to a wide range of observed feature growth rates depending on the dimensions of the masked and exposed regions. Based on a simple model, scaling relationships are derived that reveal the dynamic similarity of growth behavior across pattern designs. A time-like quantity is introduced that takes into account the varying transport effects, and provides a dimensionless time basis for analyzing crystal growth kinetics in this system. Illustrations of these scaling relationships are given through comparison with experiment. Published by Elsiver B.V.

Presented here are principles by which Sandia National Laboratories conducts its partnering activities with private industry.

The views of state of art in verification and validation (V & V) in computational physics are discussed. These views are described in the framework in which predictive capability relies on V & V, as well as other factors that affect predictive capability. Some of the research topics addressed are development of improved procedures for the use of the phenomena identification and ranking table (PIRT) for prioritizing V & V activities, and the method of manufactured solutions for code verification. It also addressed development and use of hierarchical validation diagrams, and the construction and use of validation metrics incorporating statistical measures.

Estimates of mass transfer timescales from 316 solute transport experiments reported in 35 publications are compared to the pore-water velocities and residence times, as well as the experimental durations. New tracer experiments were also conducted in columns of different lengths so that the velocity and the advective residence time could be varied independently. In both the experiments reported in the literature and the new experiments, the estimated mass transfer timescale (inverse of the mass-transfer rate coefficient) is better correlated to residence time and the experimental duration than to velocity. Of the measures considered, the experimental duration multiplied by 1 + β (where β is the capacity coefficient, defined as the ratio of masses in the immobile and mobile domains at equilibrium) best predicted the estimated mass transfer timescale. This relation is consistent with other work showing that aquifer and soil material commonly produce multiple timescales of mass transfer.

Given a finite set of points in Euclidean space, we can ask what is the minimum number of times a piecewise-linear path must change direction in order to pass through all of them. We prove some new upper and lower bounds for the rectilinear version of this problem in which all motion is orthogonal to the coordinate axes. We also consider the more general case of arbitrary directions.

Deep X-ray lithography based techniques such as LIGA (German acronym representing Lithographie, Galvanoformung, and Abformung) are being currently used to fabricate net-shape components for microelectromechanical systems (MEMS). Unlike other microfabrication techniques, LIGA lends itself to a broad range of materials, including metals, alloys, polymers, as well as ceramics and composites. Currently, Ni and Ni alloys are the materials of choice for LIGA microsystems. While Ni alloys may meet the structural requirements for MEMS, their tribological (friction and wear) behavior poses great challenges for the reliable operation of LIGA-fabricated MEMS. Typical sidewall morphologies of LIGA-fabricated parts are described, and their role in the tribological behavior of MEMS is discussed. The adaptation of commercial plasma-enhanced chemical vapor deposition to coat the sidewalls of LIGA-fabricated parts with diamond-like nanocomposite is described.

The propagation of a 30 kA, 3.5 Mev electron beam which was focused into gas and plasma-filled cells was discussed. Gas cells which were used for X-ray radiography were produced using pulsed-power accelerators, onto a high atomic number target to generate bremsstrahlung radiation. The effectiveness of beam focusing using neutral gas, partially ionized gas, and fully ionized (plasma-filled) cells was investigated using numerical simulation. It was observed in an optimized gas cell that an initial plasma density approaching 1016 cm-3 was sufficient to prevent significant net currents and the subsequent beam sweep.

Our national security, economic prosperity, and national well-being are dependent upon a set of highly interdependent critical infrastructures. Examples of these infrastructures include the national electrical grid, oil and natural gas systems, telecommunication and information networks, transportation networks, water systems, and banking and financial systems. Given the importance of their reliable and secure operations, understanding the behavior of these infrastructures - particularly when stressed or under attack - is crucial. Models and simulations can provide considerable insight into the complex nature of their behaviors and operational characteristics. These models and simulations must include interdependencies among infrastructures if they are to provide accurate representations of infrastructure characteristics and operations. A number of modeling and simulation approaches under development today directly address interdependencies and offer considerable insight into the operational and behavioral characteristics of critical infrastructures.

Acoustic testing using commercial sound system components is becoming more popular as a cost effective way of generating the required environment both in and out of a reverberant chamber. This paper will present the development of such a sound system that uses a state-of-the-art random vibration controller to perform closed-loop control in the reverberant chamber at Sandia National Laboratories. Test data will be presented that demonstrates narrow-band controlability, performance and some limitations of commercial sound generation equipment in a reverberant chamber.

We present two ways in which dynamic self-assembly can be used to perform computation, via stochastic protein networks and self-assembling software. We describe our protein-emulating agent-based simulation infrastructure, which is used for both types of computations, and the few agent properties sufficient for dynamic self-assembly. Examples of protein-network-based computation and self-assembling software are presented. We describe some novel capabilities that are enabled by the inherently dynamic nature of the self-assembling executable code. © Springer-Verlag 2004.

Finite difference equations are derived for the simulation of dielectric waveguides using an Hz -Ez formulation defined on a nonuniform triangular grid. The resulting equations may be solved as a banded eigenproblem for waveguide structures of arbitrary shape composed of regions of piecewise constant isotropic dielectric, and all transverse fields then computed from the solutions. Benchmark comparisons are presented for problems with analytic solutions, as well as a sample calculation of the propagation loss of a hollow Bragg fiber.

The structure of laminar inverse diffusion flames (IDF) of methane and ethylene in air was studied using a cylindrical co-flowing burner. IDF were similar to normal diffusion flames, except that the relative positions of the fuel and oxidizer were reversed. Radiation from soot surrounding the IDF masked the reaction zone in visible images. As a result, flame heights determined from visible images were overestimated. The height of the reaction zone as indicated by OH LIF was a more relevant measure of height. The concentration and position of PAH and soot were observed using LIF and laser-induced incandescence (LII). PAH LIF and soot LII indicated that PAH and soot are present on the fuel side of the flame, and that soot is located closer to the reaction zone than PAH. Ethylene flames produced significantly higher PAH LIF and soot LII signals than methane flames, which was consistent with the sooting propensity of ethylene. The soot and PAH were present on the fuel side of the reaction zone, but the soot was closer to the reaction zone than the PAH. This is an abstract of a paper presented at the 30th International Symposium on combustion (Chicago, IL 7/25-30/2004).

This paper describes an array of in-plane piezoelectric actuator segments laminated onto a comer-supported substrate to create a thin bimorph for reflector applications. An electric field distribution over the actuator segments causes the segments to expand or contract, thereby effecting plate deflection. To achieve a desired bimorph shape, the shape is first expressed as a two-dimensional series expansion. Then, using coefficients from the series expansion, an inverse problem is solved that determines the electric field distribution realizing the desired plate shape. A static example is presented where the desired deflection shape is a paraboloid. Copyright © 2004 by ASME.

An important component of ubiquitous computing is the ability to quickly sense the dynamic environment to learn context awareness in real-time. To pervasively capture detailed information of movements, we present a decentralized algorithm for feature extraction within a wireless sensor network. By approaching this problem in a distributed manner, we are able to work within the real constraint of wireless battery power and its effects on processing and network communications. We describe a hardware platform developed for low-power ubiquitous wireless sensing and a distributed feature extraction methodology which is capable of providing more information to the user of events while reducing power consumption. We demonstrate how the collaboration between sensor nodes can provide a means of organizing large networks into information-based clusters. © Springer-Verlag 2004.

Algorithms for multivariate image analysis and other large-scale applications of multivariate curve resolution (MCR) typically employ constrained alternating least squares (ALS) procedures in their solution. The solution to a least squares problem under general linear equality and inequality constraints can be reduced to the solution of a non-negativity-constrained least squares (NNLS) problem. Thus the efficiency of the solution to any constrained least square problem rests heavily on the underlying NNLS algorithm. We present a new NNLS solution algorithm that is appropriate to large-scale MCR and other ALS applications. Our new algorithm rearranges the calculations in the standard active set NNLS method on the basis of combinatorial reasoning. This rearrangement serves to reduce substantially the computational burden required for NNLS problems having large numbers of observation vectors. Copyright © 2005 John Wiley & Sons, Ltd.

We report the cooling of nitric oxide molecules in a single collision between an argon atom and an NO molecule at collision energies of 5.65 ± 0.36 kJ/mol and 14.7 ± 0.9 kJ/mol in a crossed molecular beam apparatus. We have produced in significant numbers (∼108 molecules cm -3 per quantum state) translationally cold NO(2Π 1/2, v′ = 0, j′ = 7.5) molecules in a specific quantum state with an upper-limit laboratory-frame rms velocity of 14.8 ± 1.1 m/s, corresponding to a temperature of 406 ± 28 mK. The translational cooling results from the kinematic collapse of the velocity distribution of the NO molecules after collision. Increasing the collision energy by increasing the velocity of the argon atoms, as we do here, does shift the scattering angle at which the cold molecules appear, but does not result in an experimentally measurable change in the velocity spread of the cold NO. This is entirely consistent with our analysis of the kinematics of the scattering which predicts that the velocity spread will actually decrease with increasing argon atom velocity. © EDP Sciences, Società, Italiana di Fisica, Springer-Verlag 2004.

Abstract not provided.

It seems well understood that supercomputer simulation is an enabler for scientific discoveries, weapons, and other activities of value to society. It also seems widely believed that Moores Law will make progressively more powerful supercomputers over time and thus enable more of these contributions. This paper seeks to add detail to these arguments, revealing them to be generally correct but not a smooth and effortless progression. This paper will review some key problems that can be solved with supercomputer simulation, showing that more powerful supercomputers will be useful up to a very high yet finite limit of around 1021 FLOPS (1 Zettaflops. The review will also show the basic nature of these extreme problems. This paper will review work by others showing that the theoretical maximum supercomputer power is very high indeed, but will explain how a straightforward extrapolation of Moores Law will lead to technological maturity in a few decades. The power of a supercomputer at the maturity of Moores Law will be very high by todays standards at 1016-1019 FLOPS (100 Petaflops to 10 Exaflops, depending on architecture , but distinctly below the level required for the most ambitious applications. Having established that Moores Law will not be that last word in supercomputing, this paper will explore the nearer term issue of what a supercomputer will look like at maturity of Moores Law. Our approach will quantify the maximum performance as permitted by the laws of physics for extension of current technology and then find a design that approaches this limit closely. We study a "multi-architecture" for supercomputers that combines a microprocessor with other "advanced" concepts and find it can reach the limits as well. This approach should be quite viable in the future because the microprocessor would provide compatibility with existing codes and programming styles while the "advanced" features would provide a boost to the limits of performance.

Many practical combustion devices and uncontrolled fires involve high Reynolds number nonpremixed turbulent flames that feature non-equilibrium finite-rate chemistry effects, e.g., local flame extinction and reignition, where enhanced transport of mass and heat away from the flame due to rapid turbulent mixing exceeds the local burning rate. Probability density function methods have shown promise in predicting piloted nonpremixed CH4-air flames over a range of Reynolds numbers and varying degrees of flame extinction and reignition. A study was carried out to quantify and characterize the kinetics of localized extinction and reignition in the Sandia flames D, E, and F, for which detailed velocity and scalar data exists. PDF methods in large eddy simulation to predict the filtered mass density function (FMDF) was used. A simple idealized mixing simulation was performed of a nonpremixed turbulent fuel jet in an air co-flow. Mixing statistics from the Monte Carlo-based FMDF solution of the chemical species scalar were compared to those from a more traditional Eulerian mixing simulation using gradient transport-based subgrid closure models. The FMDF solution will be performed with the Euclidian minimum spanning tree mixing model that uses the phenomenological connection between physical space and state space for mixing events. This is an abstract of a paper presented at the 30th International Symposium on Combustion (Chicago, IL 7/25-30/2004).

The synthesis, characterization, and separations capability of defect-free, thin-film zeolite membranes were presented. The one-micron thick sodium-aluminosilicate films of Silicalite-1 and ZSM-5 were synthesized by hydrothermal methods on either disk- or tube-supports. Techniques for growing membranes on both Al2O3 substrates as well as oxide-coated stainless steel substrates were presented. The resulting defect-free zeolite films had high flux rates at room temperature (∼ 10-7 mole/Pa-sec-sq m) and showed selective separations (3-7) between pure gases of H2 and CH4, O2, N2, CO2, CO, SF6. Results from mixed gas studies showed similar flux rates as pure gases with enhanced selectivity (15-50) for H2. The selectivity through both Silicalite-1 and ZSM-5 membranes was compared and contrasted for several gas mixtures. Data comparisons for defect-free and "defect-filled" membranes were also discussed. Under operation, the flow through these membranes quickly reached its maximum value and was stable over long periods of time. Results from experiments at high temperatures, ≤ 300°C, were compared with the data obtained at room temperature. This is an abstract of a paper presented at the 228th ACS National Meeting (Philadelphia, PA, 8/22-26/2004).

A new constitutive model relating proton conductivity to water content in a polymer electrolyte or membrane is presented. Our constitutive model is based on Faraday's law and the Nernst-Einstein equation; and it depends on the molar volumes of dry membrane and water but otherwise requires no adjustable parameters. We derive our constitutive model in two different ways. Predictions of proton conductivity as a function of membrane water content computed from our constitutive model are compared with that from a representative correlation and other models as well as experimental data from the literature and those obtained in our laboratory using a 4-point probe. Copyright © 2004 by ASME.

A coupled-physics analysis code has been developed to simulate the electrical, thermal, and mechanical responses of surface micromachined (SMM) actuators. Our objective is to optimize the design and performance of these micro actuators. Since many new designs of these electro-thermal actuators have shuttles or platforms between beams, calculating the local Joule heating requires a multi-dimensional electrostatics analysis. Moreover, the electrical solution is strongly coupled to the temperature distribution since the electrical resistivity is temperature dependent. Thus, it is essential to perform a more comprehensive simulation that solves the coupled electrostatics, thermal, and mechanical equations. Results of the coupled-physics analyses will be presented. Copyright © 2004 by ASME.

The process of removing liquid water droplets in polymer electrolyte fuel cells (PEFC) is examined using a simple analytical model and two-dimensional simulations. Specifically, the stability of a droplet adhering to the wall of the cathode flow channel is examined as a function of the geometry of the flow channel, the applied pressure gradient, and the wetting properties. The result is a prediction of the critical droplet size as a function of the difference between the advancing and receding contact angles, or contact angle hysteresis. The analytical model is shown to qualitatively predict this stability limit when compared to two-dimensional simulation results. The simulations are performed using both Arbitrary Lagrangian Eulerian (ALE) methods and level set methods. The ALE and level set predictions are shown to be in good agreement. Copyright © 2004 by ASME.

Surface micromachined structures with high aspect ratios are often utilized as sensor platforms in microelectromechanical systems (MEMS) devices. These structures generally fail by suction or adhesion to the underlying substrate during operation, or related initial processing. Such failures represent a major disadvantage in mass production of MEMS devices with highly compliant structures. Fortunately, most suction failures can be prevented or repaired in a number of ways. Passive approaches implemented during fabrication or release include: (1) utilizing special low adhesion coatings and (2) processing with low surface energy rinse agents. These methods, however, increase both the processing time and cost and are not entirely effective. Active approaches, such as illuminating stiction-failed microstructures with pulsed laser irradiation, have proven to be very effective for stiction repair [1-5]. A more recent and promising method, introduced by Gupta et al. [6], utilized laser-induced stress waves to repair stiction-failed microstructures. This approach represents a logical extension of the laser spallation technique for debonding thin films from substrates [7-9]. The method transmits stress waves into MEMS structures by laser-irradiating the back side of the substrate opposite the stiction-failed structures. This paper presents an experimental study that compares the stress wave repair method with the thermomechanical repair method on identical arrays of stiction-failed cantilevers. Copyright © 2004 by ASME.

X-ray radiography has long been recognized as a valuable tool for detecting internal features and flaws. Recent developments in microfabrication and composite materials have extended inspection requirements to the resolution limits of conventional radiography. Our work has been directed toward pushing both detection and measurement capabilities to a smaller scale. Until recently, we have used conventional contact radiography, optimized to resolve small features. With the recent purchase of a nano-focus (sub-micron) x-ray source, we are now investigating projection radiography, phase contrast imaging and micro-computed tomography (μ-CT). Projection radiography produces a magnified image that is limited in spatial resolution mainly by the source size, not by film grain size or detector pixel size. Under certain conditions phase contrast can increase the ability to resolve small features such as cracks, especially in materials with low absorption contrast. Micro-computed tomography can provide three-dimensional measurements on a micron scale and has been shown to provide better sensitivity than simple radiographs. We have included applications of these techniques to small-scale measurements not easily made by mechanical or optical means. Examples include void detection in meso-scale nickel MEMS parts, measurement of edge profiles in thick gold lithography masks, and characterization of the distribution of phases in composite materials. Our work, so far, has been limited to film. Copyright © 2004 by ASME.

Three-dimensional seismic wave propagation within a heterogeneous isotropic poroelastic medium is simulated with an explicit, time-domain, finite-difference algorithm. A system of thirteen, coupled, first-order partial differential equations is solved for the velocity vector components, stress tensor components, and pressure associated with solid and fluid constituents of the composite medium. A massively parallel computational implementation, utilizing the spatial domain decomposition strategy, allows investigation of large-scale earth models and/or broadband wave propagation within reasonable execution times.

In Combinatorial Optimization, one is frequently faced with linear programming (LP) problems with exponentially many constraints, which can be solved either using separation or what we call compact optimization. The former technique relies on a separation algorithm, which, given a fractional solution, tries to produce a violated valid inequality. Compact optimization relies on describing the feasible region of the LP by a polynomial number of constraints, in a higher dimensional space. A commonly held belief is that compact optimization does not perform as well as separation in practice. In this paper,we report on an application in which compact optimization does in fact largely outperform separation. The problem arises in structural proteomics, and concerns the comparison of 3-dimensional protein folds. Our computational results show that compact optimization achieves an improvement of up to two orders of magnitude over separation. We discuss some reasons why compact optimization works in this case but not, e.g., for the LP relaxation of the TSP. © Springer-Verlag 2004.

This report summarizes a series of structural calculations that examine effects of raising the Waste Isolation Pilot Plant repository horizon from the original design level upward 2.43 meters. These calculations allow evaluation of various features incorporated in conceptual models used for performance assessment. Material presented in this report supports the regulatory compliance re-certification, and therefore begins by replicating the calculations used in the initial compliance certification application. Calculations are then repeated for grid changes appropriate for the new horizon raised to Clay Seam G. Results are presented in three main areas: 1. Disposal room porosity, 2. Disturbed rock zone characteristics, and 3. Anhydrite marker bed failure. No change to the porosity surface for the compliance re-certification application is necessary to account for raising the repository horizon, because the new porosity surface is essentially identical. The disturbed rock zone evolution and devolution are charted in terms of a stress invariant criterion over the regulatory period. This model shows that the damage zone does not extend upward to MB 138, but does reach MB 139 below the repository. Damaged salt would be expected to heal in nominally 100 years. The anhydrite marker beds sustain states of stress that promote failure and substantial marker bed deformation into the room assures fractured anhydrite will sustain in the proximity of the disposal rooms.

This report documents the results obtained during a one-year Laboratory Directed Research and Development (LDRD) initiative aimed at investigating coupled structural acoustic interactions by means of algorithm development and experiment. Finite element acoustic formulations have been developed based on fluid velocity potential and fluid displacement. Domain decomposition and diagonal scaling preconditioners were investigated for parallel implementation. A formulation that includes fluid viscosity and that can simulate both pressure and shear waves in fluid was developed. An acoustic wave tube was built, tested, and shown to be an effective means of testing acoustic loading on simple test structures. The tube is capable of creating a semi-infinite acoustic field due to nonreflecting acoustic termination at one end. In addition, a micro-torsional disk was created and tested for the purposes of investigating acoustic shear wave damping in microstructures, and the slip boundary conditions that occur along the wet interface when the Knudsen number becomes sufficiently large.

Motivated by observations about job runtimes on the CPlant system, we use a trace-driven microsimulator to begin characterizing the performance of different classes of allocation algorithms on jobs with different communication patterns in space-shared parallel systems with mesh topology. We show that relative performance varies considerably with communication pattern. The Paging strategy using the Hilbert space-filling curve and the Best Fit heuristic performed best across several communication patterns.

We have made progress in developing a new statistical mechanics approach to designing self organizing systems that is unique to SNL. The primary application target for this ongoing research has been the development of new kinds of nanoscale components and hardware systems. However, this research also enables an out of the box connection to the field of software development. With appropriate modification, the collective behavior physics ideas for enabling simple hardware components to self organize may also provide design methods for a new class of software modules. Our current physics simulations suggest that populations of these special software components would be able to self assemble into a variety of much larger and more complex software systems. If successful, this would provide a radical (disruptive technology) path to developing complex, high reliability software unlike any known today. This high risk, high payoff opportunity does not fit well into existing SNL funding categories, as it is well outside of the mainstreams of both conventional software development practices and the nanoscience research area that spawned it. This LDRD effort was aimed at developing and extending the capabilities of self organizing/assembling software systems, and to demonstrate the unique capabilities and advantages of this radical new approach for software development.

Biological systems create proteins that perform tasks more efficiently and precisely than conventional chemicals. For example, many plants and animals produce proteins to control the freezing of water. Biological antifreeze proteins (AFPs) inhibit the solidification process, even below the freezing point. These molecules bond to specific sites at the ice/water interface and are theorized to suppress solidification chemically or geometrically. In this project, we investigated the theoretical and experimental data on AFPs and performed analyses to understand the unique physics of AFPs. The experimental literature was analyzed to determine chemical mechanisms and effects of protein bonding at ice surfaces, specifically thermodynamic freezing point depression, suppression of ice nucleation, decrease in dendrite growth kinetics, solute drag on the moving solid/liquid interface, and stearic pinning of the ice interface. Stearic pinning was found to be the most likely candidate to explain experimental results, including freezing point depression, growth morphologies, and thermal hysteresis. A new stearic pinning model was developed and applied to AFPs, with excellent quantitative results. Understanding biological antifreeze mechanisms could enable important medical and engineering applications, but considerable future work will be necessary.

An estimate of the distribution of fatigue ranges or extreme loads for wind turbines may be obtained by separating the problem into two uncoupled parts, (1) a turbine specific portion, independent of the site and (2) a site-specific description of environmental variables. We consider contextually appropriate probability models to describe the turbine specific response for extreme loads or fatigue. The site-specific portion is described by a joint probability distribution of a vector of environmental variables, which characterize the wind process at the hub-height of the wind turbine. Several approaches are considered for combining the two portions to obtain an estimate of the extreme load, e.g., 50-year loads or fatigue damage. We assess the efficacy of these models to obtain accurate estimates, including various levels of epistemic uncertainty, of the turbine response.

The quantitative analysis of ammonia binding sites in the Davison (Type 3A) zeolite desiccant using solid-state {sup 15}N MAS NMR spectroscopy is reported. By utilizing 15N enriched ammonia ({sup 15}NH{sub 3}) gas, the different adsorption/binding sites within the zeolite were investigated as a function of NH{sub 3} loading. Using {sup 15}N MAS NMR multiple sites were resolved that have distinct cross-polarization dynamics and chemical shift behavior. These differences in the {sup 15}N NMR were used to characterize the adsorption environments in both the pure 3A zeolite and the silicone-molded forms of the desiccant.

This SAND report provides the technical progress through October 2003 of the Sandia-led project, 'Carbon Sequestration in Synechococcus Sp.: From Molecular Machines to Hierarchical Modeling,' funded by the DOE Office of Science Genomes to Life Program. Understanding, predicting, and perhaps manipulating carbon fixation in the oceans has long been a major focus of biological oceanography and has more recently been of interest to a broader audience of scientists and policy makers. It is clear that the oceanic sinks and sources of CO2 are important terms in the global environmental response to anthropogenic atmospheric inputs of CO2 and that oceanic microorganisms play a key role in this response. However, the relationship between this global phenomenon and the biochemical mechanisms of carbon fixation in these microorganisms is poorly understood. In this project, we will investigate the carbon sequestration behavior of Synechococcus Sp., an abundant marine cyanobacteria known to be important to environmental responses to carbon dioxide levels, through experimental and computational methods. This project is a combined experimental and computational effort with emphasis on developing and applying new computational tools and methods. Our experimental effort will provide the biology and data to drive the computational efforts and include significant investment in developing new experimental methods for uncovering protein partners, characterizing protein complexes, identifying new binding domains. We will also develop and apply new data measurement and statistical methods for analyzing microarray experiments. Computational tools will be essential to our efforts to discover and characterize the function of the molecular machines of Synechococcus. To this end, molecular simulation methods will be coupled with knowledge discovery from diverse biological data sets for high-throughput discovery and characterization of protein-protein complexes. In addition, we will develop a set of novel capabilities for inference of regulatory pathways in microbial genomes across multiple sources of information through the integration of computational and experimental technologies. These capabilities will be applied to Synechococcus regulatory pathways to characterize their interaction map and identify component proteins in these pathways. We will also investigate methods for combining experimental and computational results with visualization and natural language tools to accelerate discovery of regulatory pathways. The ultimate goal of this effort is develop and apply new experimental and computational methods needed to generate a new level of understanding of how the Synechococcus genome affects carbon fixation at the global scale. Anticipated experimental and computational methods will provide ever-increasing insight about the individual elements and steps in the carbon fixation process, however relating an organism's genome to its cellular response in the presence of varying environments will require systems biology approaches. Thus a primary goal for this effort is to integrate the genomic data generated from experiments and lower level simulations with data from the existing body of literature into a whole cell model. We plan to accomplish this by developing and applying a set of tools for capturing the carbon fixation behavior of complex of Synechococcus at different levels of resolution. Finally, the explosion of data being produced by high-throughput experiments requires data analysis and models which are more computationally complex, more heterogeneous, and require coupling to ever increasing amounts of experimentally obtained data in varying formats. These challenges are unprecedented in high performance scientific computing and necessitate the development of a companion computational infrastructure to support this effort. More information about this project, including a copy of the original proposal, can be found at www.genomes-to-life.org

Military test and training ranges operate with live fire engagements to provide realism important to the maintenance of key tactical skills. Ordnance detonations during these operations typically produce minute residues of parent explosive chemical compounds. Occasional low order detonations also disperse solid phase energetic material onto the surface soil. These detonation remnants are implicated in chemical contamination impacts to groundwater on a limited set of ranges where environmental characterization projects have occurred. Key questions arise regarding how these residues and the environmental conditions (e.g., weather and geostratigraphy) contribute to groundwater pollution impacts. This report documents interim results of experimental work evaluating mass transfer processes from solid phase energetics to soil pore water. The experimental work is used as a basis to formulate a mass transfer numerical model, which has been incorporated into the porous media simulation code T2TNT. This report documents the results of the Phase III experimental effort, which evaluated the impacts of surface deposits versus buried deposits, energetic material particle size, and low order detonation debris. Next year, the energetic material mass transfer model will be refined and a 2-d screening model will be developed for initial site-specific applications. A technology development roadmap was created to show how specific R&D efforts are linked to technology and products for key customers.

A Micro Electro Mechanical System (MEMS) typically consists of micron-scale parts that move through a gas at atmospheric or reduced pressure. In this situation, the gas-molecule mean free path is comparable to the geometric features of the microsystem, so the gas flow is noncontinuum. When mean-free-path effects cannot be neglected, the Boltzmann equation must be used to describe the gas flow. Solution of the Boltzmann equation is difficult even for the simplest case because of its sevenfold dimensionality (one temporal dimension, three spatial dimensions, and three velocity dimensions) and because of the integral nature of the collision term. The Direct Simulation Monte Carlo (DSMC) method is the method of choice to simulate high-speed noncontinuum flows. However, since DSMC uses computational molecules to represent the gas, the inherent statistical noise must be minimized by sampling large numbers of molecules. Since typical microsystem velocities are low (< 1 m/s) compared to molecular velocities ({approx}400 m/s), the number of molecular samples required to achieve 1% precision can exceed 1010 per cell. The Discrete Velocity Gas (DVG) method, an approach motivated by radiation transport, provides another way to simulate noncontinuum gas flows. Unlike DSMC, the DVG method restricts molecular velocities to have only certain discrete values. The transport of the number density of a velocity state is governed by a discrete Boltzmann equation that has one temporal dimension and three spatial dimensions and a polynomial collision term. Specification and implementation of DVG models are discussed, and DVG models are applied to Couette flow and to Fourier flow. While the DVG results for these benchmark problems are qualitatively correct, the errors in the shear stress and the heat flux can be order-unity even for DVG models with 88 velocity states. It is concluded that the DVG method, as described herein, is not sufficiently accurate to simulate the low-speed gas flows that occur in microsystems.

CommAspen is a new agent-based model for simulating the interdependent effects of market decisions and disruptions in the telecommunications infrastructure on other critical infrastructures in the U.S. economy such as banking and finance, and electric power. CommAspen extends and modifies the capabilities of Aspen-EE, an agent-based model previously developed by Sandia National Laboratories to analyze the interdependencies between the electric power system and other critical infrastructures. CommAspen has been tested on a series of scenarios in which the communications network has been disrupted, due to congestion and outages. Analysis of the scenario results indicates that communications networks simulated by the model behave as their counterparts do in the real world. Results also show that the model could be used to analyze the economic impact of communications congestion and outages.

This report is the latest in a continuing series that highlights the recent technical accomplishments associated with the work being performed within the Materials and Process Sciences Center. Our research and development activities primarily address the materials-engineering needs of Sandia's Nuclear-Weapons (NW) program. In addition, we have significant efforts that support programs managed by the other laboratory business units. Our wide range of activities occurs within six thematic areas: Materials Aging and Reliability, Scientifically Engineered Materials, Materials Processing, Materials Characterization, Materials for Microsystems, and Materials Modeling and Simulation. We believe these highlights collectively demonstrate the importance that a strong materials-science base has on the ultimate success of the NW program and the overall DOE technology portfolio.

The catalytic combustion of natural gas has been the topic of much research over the past decade. Interest in this technology results from a desire to decrease or eliminate the emissions of harmful nitrogen oxides (NOX) from gas turbine power plants. A low-pressure drop catalyst support, such as a ceramic monolith, is ideal for this high-temperature, high-flow application. A drawback to the traditional honeycomb monoliths under these operating conditions is poor mass transfer to the catalyst surface in the straight-through channels. 'Robocasting' is a unique process developed at Sandia National Laboratories that can be used to manufacture ceramic monoliths with alternative 3-dimensional geometries, providing tortuous pathways to increase mass transfer while maintaining low pressure drops. This report details the mass transfer effects for novel 3-dimensional robocast monoliths, traditional honeycomb-type monoliths, and ceramic foams. The mass transfer limit is experimentally determined using the probe reaction of CO oxidation over a Pt / {gamma}-Al{sub 2}O{sub 3} catalyst, and the pressure drop is measured for each monolith sample. Conversion versus temperature data is analyzed quantitatively using well-known dimensionless mass transfer parameters. The results show that, relative to the honeycomb monolith support, considerable improvement in mass transfer efficiency is observed for robocast samples synthesized using an FCC-like geometry of alternating rods. Also, there is clearly a trade-off between enhanced mass transfer and increased pressure drop, which can be optimized depending on the particular demands of a given application.

This document introduces the use of Trilinos, version 3.1. Trilinos has been written to support, in a rigorous manner, the solver needs of the engineering and scientific applications at Sandia National Laboratories. Aim of this manuscript is to present the basic features of some of the Trilinos packages. The presented material includes the definition of distributed matrices and vectors with Epetra, the iterative solution of linear system with AztecOO, incomplete factorizations with IFPACK, multilevel methods with ML, direct solution of linear system with Amesos, and iterative solution of nonlinear systems with NOX. With the help of several examples, some of the most important classes and methods are detailed to the inexperienced user. For the most majority, each example is largely commented throughout the text. Other comments can be found in the source of each example. This document is a companion to the Trilinos User's Guide and Trilinos Development Guides. Also, the documentation included in each of the Trilinos' packages is of fundamental importance.