This report describes an LDRD-supported experimental-theoretical collaboration on the enhanced low-dose-rate sensitivity (ELDRS) problem. The experimental work led to a method for elimination of ELDRS, and the theoretical work led to a suite of bimolecular mechanisms that explain ELDRS and is in good agreement with various ELDRS experiments. The model shows that the radiation effects are linear in the limit of very low dose rates. In this limit, the regime of most concern, the model provides a good estimate of the worst-case effects of low dose rate ionizing radiation.
This document describes the 2003 SNL ASCI Software Quality Engineering (SQE) assessment of twenty ASCI application code teams and the results of that assessment. The purpose of this assessment was to determine code team compliance with the Sandia National Laboratories ASCI Applications Software Quality Engineering Practices, Version 2.0 as part of an overall program assessment.
An increase in photocurrent has been observed at silicon electrodes coated with nanostructured porous silica films as compared to bare, unmodified silicon. Ultimately, to utilize this effect in devices such as sensors or microchip power supplies, the physical phenomena behind this observation need to be well characterized. To this end, Electrochemical Impedance Spectroscopy (EIS) was used to characterize the effect of surfactant-templated mesoporous silica films deposited onto silicon electrodes on the electrical properties of the electrode space-charge region in an aqueous electrolyte solution, as the electrical properties of this space-charge region are responsible for the photobehavior of semiconductor devices. A significant shift in apparent flat-band potential was observed for electrodes modified with the silica film when compared to bare electrodes; the reliability of this data is suspect, however, due to contributions from surface states to the overall capacitance of the system. To assist in the interpretation of this EIS data, a series of measurements at Pt electrodes was performed with the hope of decoupling electrode and film contributions from the EIS spectra. Surprisingly, the frequency-dependent impedance data for Pt electrodes coated with a surfactant-templated film was nearly identical to that observed for bare Pt electrodes, indicating that the mesoporous film had little effect on the transport of small electrolyte ions to the electrode surface. Pore-blocking agents (tetraalkylammonium salts) were not observed to inhibit this transport process. However, untemplated (non-porous) silica films dramatically increased film resistance, indicating that our EIS data for the Pt electrodes is reliable. Overall, our preliminary conclusion is that a shift in electrical properties in the space-charge region induced by the presence of a porous silica film is responsible for the increase in observed photocurrent.
The waters of the Pecos River in New Mexico must be delivered to three primary users: (1) The Pecos River Compact: each year a percentage of water from natural river flow must be delivered to Texas; (2) Agriculture: Carlsbad Irrigation District has a storage and diversion right and Fort Sumner Irrigation District has a direct flow diversion right; and, (3) Endangered Species Act: an as yet unspecified amount of water is to support Pecos Bluntnose Shiner Minnow habitat within and along the Pecos River. Currently, the United States Department of Interior Bureau of Reclamation, the New Mexico Interstate Stream Commission, and the United States Department of the Interior Fish and Wildlife Service are studying the Pecos Bluntnose Shiner Minnow habitat preference. Preliminary work by Fish and Wildlife personnel in the critical habitat suggest that water depth and water velocity are key parameters defining minnow habitat preference. However, river flows that provide adequate preferred habitat to support this species have yet to be determined. Because there is a limited amount of water in the Pecos River and its reservoirs, it is critical to allocate water efficiently such that habitat is maintained, while honoring commitments to agriculture and to the Pecos River Compact. This study identifies the relationship between Pecos River flow rates in cubic feet per second (cfs) and water depth and water velocity.
We have developed infrastructure, utilities and partitioning methods to improve data partitioning in linear solvers and preconditioners. Our efforts included incorporation of data repartitioning capabilities from the Zoltan toolkit into the Trilinos solver framework, (allowing dynamic repartitioning of Trilinos matrices); implementation of efficient distributed data directories and unstructured communication utilities in Zoltan and Trilinos; development of a new multi-constraint geometric partitioning algorithm (which can generate one decomposition that is good with respect to multiple criteria); and research into hypergraph partitioning algorithms (which provide up to 56% reduction of communication volume compared to graph partitioning for a number of emerging applications). This report includes descriptions of the infrastructure and algorithms developed, along with results demonstrating the effectiveness of our approaches.
A laser safety and hazard analysis was performed for the airborne AURA (Big Sky Laser Technology) lidar system based on the 2000 version of the American National Standard Institute's (ANSI) Standard Z136.1, for the Safe Use of Lasers and the 2000 version of the ANSI Standard Z136.6, for the Safe Use of Lasers Outdoors. The AURA lidar system is installed in the instrument pod of a Proteus airframe and is used to perform laser interaction experiments and tests at various national test sites. The targets are located at various distances or ranges from the airborne platform. In order to protect personnel, who may be in the target area and may be subjected to exposures, it was necessary to determine the Maximum Permissible Exposure (MPE) for each laser wavelength, calculate the Nominal Ocular Hazard Distance (NOHD), and determine the maximum 'eye-safe' dwell times for various operational altitudes and conditions. It was also necessary to calculate the appropriate minimum Optical Density (ODmin) of the laser safety eyewear used by authorized personnel who may receive hazardous exposures during ground base operations of the airborne AURA laser system (system alignment and calibration).
The Advanced Concepts Group (ACG) at Sandia National Laboratories is exploring the use of Red Teaming to help intelligence analysts with two key processes: determining what a piece or pieces of information might imply and deciding what other pieces of information need to be found to support or refute hypotheses about what actions a suspected terrorist organization might be pursuing. In support of this effort, the ACG hosted a terrorism red gaming event in Albuquerque on July 22-24, 2003. The game involved two 'red teams' playing the roles of two terrorist cells - one focused on implementing an RDD attack on the DC subway system and one focused on a bio attack against the same target - and two 'black teams' playing the role of the intelligence collection system and of intelligence analysts trying to decide what plans the red teams might be pursuing. This exercise successfully engaged human experts to seed a proposed compute engine with detailed operational plans for hypothetical terrorist scenarios.
The Accurate Time-Linked data Acquisition System (ATLAS II) is a small, lightweight, time-synchronized, robust data acquisition system that is capable of acquiring simultaneous long-term time-series data from both a wind turbine rotor and ground-based instrumentation. This document is a user's manual for the ATLAS II hardware and software. It describes the hardware and software components of ATLAS II, and explains how to install and execute the software.
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.
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.
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.
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).
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).
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.
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.