Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Pyrolyzed carbon as a mechanical material is promising for applications in harsh environments. In this work, we characterized the material and developed novel processes for fabricating carbon composite micro-electromechanical systems (CMEMS) structures. A novel method of increasing Young's modulus and the conductivity of pyrolyzed AZ 4330 was demonstrated by loading the films with graphene oxide prior to pyrolysis. By incorporating 2 wt.% graphene stiffeners into the film, a 65% increase in Young's modulus and 11% increase in conductivity were achieved. By reactive ion etching pyrolyzed blanket AZ 50XT thick film photoresist, a high aspect ratio process was demonstrated with films >7.5um thick. Two novel multi-level, volume-scalable CMEMS processes were developed on 6" diameter wafers. Young's modulus of 23 GPa was extracted from nanoindentation measurements of pyrolyzed AZ 50XT films. The temperature-dependent resistance was characterized from room temperature to 500C and found to be nearly linear over this range. By fitting the results of self-heated bridges in an inert ambient, we calculated that the bridges survived to 1000C without failure. Transmission electron microscopy (TEM) results showed the film to be largely amorphous, containing some sub-micrometer sized graphite crystallites. This was consistent with our Raman analysis, which also showed the film to be largely sp² bonded. The calculated average density of pyrolyzed AZ 4330 films was 1.32 g/cm². Thin level of disorder and the conductivity of thin film resistors were found to unchanged by 2Mrad gamma irradiation from a Co₆₀ source. Thin film pyrolyzed carbon resistors were hermetically sealed in a nitrogen ambient in 24-pin dual in-line packages (DIP's). The resistance was measured periodically and remained constant over 6 months' time.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Thermodynamic quantities, such as pressure and internal energy, and their derivatives, are used in many applications. Depending on application, a natural set of quantities related to one of four thermodynamic potentials are typically used. For example, hydro-codes use internal energy derived quantities and Equation of State work often uses Helmholtz free energy quantities. When performing work spanning over several fields, transformations between one set of quantities and another set of quantities are often needed. A short, but comprehensive, review of such transformations are given in this report.

Abstract not provided.

Here, we describe a new method to measure the activation energy required to remove a strongly-bound membrane-associated protein from a lipid membrane (anchoring energy). It is based on measuring the rate of release of a liposome-bound protein during centrifugation on a sucrose gradient as a function of time and temperature. The method was used to determine anchoring energy for the soluble dengue virus envelope protein (sE) strongly bound to 80:20 POPC:POPG liposomes at pH 5.5. We also measured the binding energy of sE at the same pH for the initial, predominantly reversible, phase of binding to a 70:30 PC:PG lipid bilayer. The anchoring energy (37 +/- 1.7 kcal/mol, 20% PG) was found to be much larger than the binding energy (7.8 +/- 0.3 kcal/mol for 30% PG, or est. 7.0 kcal/mol for 20% PG). This is consistent with data showing that free sE is a monomer at pH 5.5, but assembles into trimers after associating with membranes. But, trimerization alone is insufficient to account for the observed difference in energies, and we conclude that some energy dissipation occurs during the release process. This new method to determine anchoring energy should be useful to understand the complex interactions of integral monotopic proteins and strongly-bound peripheral membrane proteins with lipid membranes.

The goal of the Domestic Nuclear Detection Office (DNDO) Algorithm Improvement Program (AIP) is to facilitate gamma-radiation detector nuclide identification algorithm development, improvement, and validation. Accordingly, scoring criteria have been developed to objectively assess the performance of nuclide identification algorithms. In addition, a Microsoft Excel spreadsheet application for automated nuclide identification scoring has been developed. This report provides an overview of the equations, nuclide weighting factors, nuclide equivalencies, and configuration weighting factors used by the application for scoring nuclide identification algorithm performance. Furthermore, this report presents a general overview of the nuclide identification algorithm scoring application including illustrative examples.

Performance portability on heterogeneous high-performance computing (HPC) systems is a major challenge faced today by code developers: parallel code needs to execute correctly as well as with high performance on machines with different architectures, operating systems, and software libraries. The Finite Element Method (FEM) is a popular and flexible method for discretizing partial differential equations arising in a wide variety of scientific, engineering, and industry applications that require HPC. This paper presents some preliminary results pertaining to our development of a performance portable implementation of the FEM-based Albany code. Performance portability is achieved using the Kokkos library of Trilinos. We present performance results for two different physics simulations modules in Albany: the Aeras global atmosphere dynamical code and the FELIX land-ice solver. As a result, numerical experiments show that our single code implementation gives reasonable performance across two multi-core/many-core architectures: NVIDIA GPUs and multi-core CPUs.

We introduce a discretization for a nonlocal diffusion problem using a localized basis of radial basis functions. The stiffness matrix entries are assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, sparse, symmetric positive definite stiffness matrix. We demonstrate that both the continuum and discrete problems are well-posed and present numerical results for the convergence behavior of the radial basis function method. As a result, we explore approximating the solution to anisotropic differential equations by solving anisotropic nonlocal integral equations using the radial basis function method.

Abstract not provided.

Geologic carbon storage in deep saline aquifers is a promising technology for reducing anthropogenic emissions into the atmosphere. Dissolution of injected CO₂ into resident brines is one of the primary trapping mechanisms generally considered necessary to provide long-term storage security. Given that diffusion of CO₂ in brine is woefully slow, convective dissolution, driven by a small increase in brine density with CO₂ saturation, is considered to be the primary mechanism of dissolution trapping. Previous studies of convective dissolution have typically only considered the convective process in the single-phase region below the capillary transition zone and have either ignored the overlying two-phase region where dissolution actually takes place or replaced it with a virtual region with reduced or enhanced constant permeability. Our objective is to improve estimates of the long-term dissolution flux of CO₂ into brine by including the capillary transition zone in two-phase model simulations. In the fully two-phase model, there is a capillary transition zone above the brine-saturated region over which the brine saturation decreases with increasing elevation. Our two-phase simulations show that the dissolution flux obtained by assuming a brine-saturated, single-phase porous region with a closed upper boundary is recovered in the limit of vanishing entry pressure and capillary transition zone. For typical finite entry pressures and capillary transition zone, however, convection currents penetrate into the two-phase region. As a result, this removes the mass transfer limitation of the diffusive boundary layer and enhances the convective dissolution flux of CO₂ more than 3 times above the rate assuming single-phase conditions.

We present a verification and validation analysis of a coordinate-transformation-based numerical solution method for the two-dimensional axisymmetric magnetic diffusion equation, implemented in the finite-element simulation code ALEGRA. The transformation, suggested by Melissen and Simkin, yields an equation set perfectly suited for linear finite elements and for problems with large jumps in material conductivity near the axis. The verification analysis examines transient magnetic diffusion in a rod or wire in a very low conductivity background by first deriving an approximate analytic solution using perturbation theory. This approach for generating a reference solution is shown to be not fully satisfactory. A specialized approach for manufacturing an exact solution is then used to demonstrate second-order convergence under spatial refinement and tem- poral refinement. For this new implementation, a significant improvement relative to previously available formulations is observed. Benefits in accuracy for computed current density and Joule heating are also demonstrated. The validation analysis examines the circuit-driven explosion of a copper wire using resistive magnetohydrodynamics modeling, in comparison to experimental tests. The new implementation matches the accuracy of the existing formulation, with both formulations capturing the experimental burst time and action to within approximately 2%.

Herein we develop a quantitative dye dequenching technique for the measurement of polymersome fusion, using it to characterize the salt mediated, mechanically-induced fusion of polymersomes with polymer, lipid, and so-called stealth lipid vesicles. While dye dequenching has been used to quantitatively explore liposome fusion in the past, this is the first use of dye dequenching to measure polymersome fusion of which we are aware. In addition to providing quantitative results, dye dequenching is ideal for detecting fusion in instances where DLS results would be ambiguous, such as low yield levels and size ranges outside the capabilities of DLS. The dye chosen for this study was a cyanine derivative, 1,1′-dioctadecyl-3,3,3′,3′-tetramethylindotricarbocyanine iodide (DiR), which proved to provide excellent data on the extent of polymersome fusion. Using this technique, we have shown the limited fusion capabilities of polymersome/liposome heterofusion, notably DOPC vesicles fusing with polymersomes at half the efficiency of polymersome homofusion and DPPC vesicles showing virtually no fusion. In addition to these key heterofusion experiments, we determined the broad applicability of dye dequenching in measuring kinetic rates of polymersome fusion; and showed that even at elevated temperatures or over multiple weeks' time, no polymersome fusion occurred without agitation. Stealth liposomes formed from DPPC and PEO-functionalized lipid, however, fused with polymersomes and stealth liposomes with relatively high efficiency, lending support to our hypothesis that the response of the PEO corona to salt is a key factor in the fusion process. This last finding suggests that although the conjugation of PEO to lipids increases vesicle biocompatibility and enables their longer circulation times, it also renders the vesicles subject to destabilization under high salt and shear (e.g. in the circulatory system) that may lead to, in this case, fusion.

With vibrant colours and simple, roomerature processing methods, electrochromic polymers have attracted attention as active materials for flexible, low-power-consuming devices. However, slow switching speeds in devices realized to date, as well as the complexity of having to combine several distinct polymers to achieve a full-colour gamut, have limited electrochromic materials to niche applications. Here we achieve fast, high-contrast electrochromic switching by significantly enhancing the interaction of light-propagating as deep-subwavelength-confined surface plasmon polaritons through arrays of metallic nanoslits, with an electrochromic polymer-present as an ultra-thin coating on the slit sidewalls. The switchable configuration retains the short temporal charge-diffusion characteristics of thin electrochromic films, while maintaining the high optical contrast associated with thicker electrochromic coatings. We further demonstrate that by controlling the pitch of the nanoslit arrays, it is possible to achieve a full-colour response with high contrast and fast switching speeds, while relying on just one electrochromic polymer.

We have investigated the atomic structure of graphene/Ir(111) supported platinum clusters with on average fewer than 40 atoms by means of surface x-ray diffraction (SXRD), grazing incidence small angle x-ray scattering (GISAXS), and normal incidence x-ray standing waves (NIXSW) measurements, in comparison with density functional theory calculations (DFT). GISAXS revealed that the clusters with 1.3 nm diameter form a regular array with domain sizes of 90 nm. SXRD shows that the 1-2 monolayer high, (111) oriented Pt nanoparticles grow epitaxially on the graphene support. From the combined analysis of the SXRD and NIXSW data, a three-dimensional (3D) structural model of the clusters and the graphene support can be deduced which is in line with the DFT results. For the clusters grown in ultrahigh vacuum the lattice parameter is reduced by (4.6±0.1)% compared to bulk platinum. The graphene layer undergoes a strong Pt adsorption induced buckling, caused by a rehybridization of the carbon atoms below the cluster. In situ observation of the Pt clusters in CO and O2 environments revealed a reversible change of the clusters' strain state while successively dosing CO at room temperature and O2 at 575 K, pointing to a CO oxidation activity of the Pt clusters.

The High Optical Access (HOA) trap was designed in collaboration with the Modular Universal Scalable Ion-trap Quantum Computer (MUSIQC) team, funded along with Sandia National Laboratories through IARPA's Multi Qubit Coherent Operations (MQCO) program. The design of version 1 of the HOA trap was completed in September 2012 and initial devices were completed and packaged in February 2013. The second version of the High Optical Access Trap (HOA-2) was completed in September 2014 and is available at IARPA's disposal.

Many ideas for liquid surface PFCs are for divertors. First walls are likely to be more challenging technologically because long flow paths are necessary for fast flowing systems and the first wall must be an integral structure with the blanket. Maximum tolerable heat loads are a critical concern. This paper describes several processes at work in walls with fast-flowing or slow-flowing liquid plasma-facing surfaces, and the considerations imposed by heat transfer and the power balance for the PFC as well as the structure needed for an integrated first wall and blanket, and uses thermal modeling of a generic PFC structure to illustrate the issues and support the conclusions.

In this study, collection of mosquitoes and testing for vector-borne viruses is a key surveillance activity that directly influences the vector control efforts of public health agencies, including determining when and where to apply insecticides. Vector control districts in California routinely monitor for three human pathogenic viruses including West Nile virus (WNV), Western equine encephalitis virus (WEEV), and St. Louis encephalitis virus (SLEV). Reverse transcription quantitative polymerase chain reaction (RT-qPCR) offers highly sensitive and specific detection of these three viruses in a single multiplex reaction, but this technique requires costly, specialized equipment that is generally only available in centralized public health laboratories. We report the use of reverse transcription loop-mediated isothermal amplification (RT-LAMP) to detect WNV, WEEV, and SLEV RNA extracted from pooled mosquito samples collected in California, including novel primer sets for specific detection of WEEV and SLEV, targeting the nonstructural protein 4 (nsP4) gene of WEEV and the 3’ untranslated region (3’-UTR) of SLEV. Our WEEV and SLEV RT-LAMP primers allowed detection of <0.1 PFU/reaction of their respective targets in <30 minutes, and exhibited high specificity without cross reactivity when tested against a panel of alphaviruses and flaviviruses. Furthermore, the SLEV primers do not cross-react with WNV, despite both viruses being closely related members of the Japanese encephalitis virus complex. The SLEV and WEEV primers can also be combined in a single RT-LAMP reaction, with discrimination between amplicons by melt curve analysis. Although RT-qPCR is approximately one order of magnitude more sensitive than RT-LAMP for all three targets, the RT-LAMP technique is less instrumentally intensive than RT-qPCR and provides a more cost-effective method of vector-borne virus surveillance.

In this study, the diffusion of water and ions in the interlayer region of smectite clay minerals represents a direct probe of the type and strength of clay–fluid interactions. Interlayer diffusion also represents an important link between molecular simulation and macroscopic experiments. Here we use molecular dynamics simulation to investigate trends in cation and water diffusion in montmorillonite interlayers, looking specifically at the effects of layer charge, interlayer cation and cation charge (sodium or calcium), water content, and temperature. For Na-montmorillonite, the largest increase in ion and water diffusion coefficients occurs between the one-layer and two-layer hydrates, corresponding to the transition from inner-sphere to outer-sphere surface complexes. Calculated activation energies for ion and water diffusion in Na-montmorillonite are similar to each other and to the water hydrogen bond energy, suggesting the breaking of water–water and water–clay hydrogen bonds as a likely mechanism for interlayer diffusion. A comparison of interlayer diffusion with that of bulk electrolyte solutions reveals a clear trend of decreasing diffusion coefficient with increasing electrolyte concentration, and in most cases the interlayer diffusion results are nearly coincident with the corresponding bulk solutions. Trends in electrical conductivities computed from the ion diffusion coefficients are also compared.

Ternary polymer brushes consisting of polystyrene, poly(methyl methacrylate), and poly(4-vinylpyridine) have been synthesized. These brushes laterally phase separate into several distinct phases and can be tailored by altering the relative polymer composition. Self-consistent field theory has been used to predict the phase diagram and model both the horizontal and vertical phase behavior of the polymer brushes. All phase behaviors observed experimentally correlate well with the theoretical model.

In this study, a mechanical model is introduced for predicting the initiation and evolution of complex fracture patterns without the need for damage variable or law. The model, a continuum variant of Newton's second law, uses integral rather than partial differential operators where the region of integration is over finite domain. The force interaction is derived from a novel nonconvex strain energy density function, resulting in a nonmonotomic material model. The resulting equation of motion is proved to be mathematically well-posed. The model has the capacity to simulate nucleation and growth of multiple, mutually interacting dynamic fractures. In the limit of zero region of integration, the model reproduces the classic Griffith model of brittle fracture. The simplicity of the formulation avoids the need for supplemental kinetic relations that dictate crack growth or the need for an explicit damage evolution law.

Environmental contours describing extreme sea states are generated as the input for numerical or physical model simulations as a part of the standard current practice for designing marine structures to survive extreme sea states. These environmental contours are characterized by combinations of significant wave height (Hs) and either energy period (Te) or peak period (Tp) values calculated for a given recurrence interval using a set of data based on hindcast simulations or buoy observations over a sufficient period of record. The use of the inverse first-order reliability method (I-FORM) is a standard design practice for generating environmental contours. This paper develops enhanced methodologies for data analysis prior to the application of the I-FORM, including the use of principal component analysis (PCA) to create an uncorrelated representation of the variables under consideration as well as new distribution and parameter fitting techniques. These modifications better represent the measured data and, therefore, should contribute to the development of more realistic representations of environmental contours of extreme sea states for determining design loads for marine structures.

High-temperature Arrhenius ignition delay time correlations are useful for revealing the underlying parameter dependencies of combustion models, for simplifying and optimizing combustion mechanisms for use in engine simulations, for scaling experimental data to new conditions for comparison purposes, and for guiding in experimental design. We have developed a scaling relationship for Fatty Acid Methyl Ester (FAME) ignition time data taken at high temperatures in 4%O2/Ar mixtures behind reflected shocks using an aerosol shock tube:τign[ms]=2.24×10-6[ms](P[atm])-0.41(φ)0.30(Cn)-0.61exp37.1[kcal/mol]Ru[kcal/mol K]T[K]Additionally, we have combined our ignition delay time data for methyl decanoate, methyl palmitate, methyl oleate, and methyl linoleate with other experimental results in the literature in order to derive fuel-specific oxygen-mole-fraction scaling parameters for these surrogates. In this article, we discuss the significance of the parameter values, compare our correlation to others found in the literature for different classes of fuels, and contrast the above expression's performance with correlations obtained using leading FAME kinetic models in 4%O2/Ar mixtures.

Nanosecond pulsed laser irradiation was used to fabricate colored, mechanically robust oxide "tags" on 304L stainless steel. Immersion in simulated seawater solution, salt fog exposure, and anodic polarization in a 3.5% NaCl solution were employed to evaluate the environmental resistance of these oxide tags. Single layer oxides outside a narrow thickness range (~100-150 nm) are susceptible to dissolution in chloride containing environments. The 304L substrates immediately beneath the oxides corrode severely-attributed to Cr-depletion in the melt zone during laser processing. For the first time, multilayered oxides were fabricated with pulsed laser irradiation in an effort to expand the protective thickness range while also increasing the variety of film colors attainable in this range. Layered films grown using a laser scan rate of 475 mm/s are more resistant to both localized and general corrosion than oxides fabricated at 550 mm/s. In the absence of pre-processing to mitigate Cr-depletion, layered films can enhance environmental stability of the system.

High temperature operation (250-340 K) of short-wavelength interband cascade infrared photodetectors (ICIPs) with InAs/GaSb/Al0.2In0.8Sb/GaSb superlattice absorbers has been demonstrated with a 50% cutoff wavelength of 2.9 μm at 300 K. Two ICIP structures, one with two and the other with three stages, were designed and grown to explore this multiple-stage architecture. At λ = 2.1 μm, the two- and three-stage ICIPs had Johnson-noise-limited detectivities of 5.1 × 109 and 5.8 × 109cm Hz1/2/W, respectively, at 300 K. The better device performance of the three-stage ICIP over the two-stage ICIP confirmed the advantage of more stages for this cascade architecture. An Arrhenius activation energy of 450 meV is extracted for the bulk resistance-area product, which indicates the dominance of the diffusion current at these high temperatures.

Laser-based failure analysis techniques demonstrate the ability to quickly and non-intrusively screen deep ultraviolet light-emitting diodes (LEDs) for electrically-active defects. In particular, two laser-based techniques, light-induced voltage alteration and thermally-induced voltage alteration, generate applied voltage maps (AVMs) that provide information on electrically-active defect behavior including turn-on bias, density, and spatial location. Here, multiple commercial LEDs were examined and found to have dark defect signals in the AVM indicating a site of reduced resistance or leakage through the diode. The existence of the dark defect signals in the AVM correlates strongly with an increased forward-bias leakage current. This increased leakage is not present in devices without AVM signals. Transmission electron microscopy analysis of a dark defect signal site revealed a dislocation cluster through the pn junction. The cluster included an open core dislocation. Even though LEDs with few dark AVM defect signals did not correlate strongly with power loss, direct association between increased open core dislocation densities and reduced LED device performance has been presented elsewhere [M. W. Moseley et al., J. Appl. Phys. 117, 095301 (2015)].

During the course of their careers, welding engineers and welding metallurgists are often confronted with questions regarding welding process and properties that on the surface appear to be simple and direct, but are in fact quite challenging. These questions generally mask an underlying complexity whose underpinnings in scientific and applied research predate even the founding of the American Welding Society, and previous Comfort A. Adams lectures provide ample and fascinating evidence of the breadth and depth of this complexity. Using these studies or their own experiences and investigations as a basis, most welding and materials engineers have developed engineering tools to provide working approaches to these day-to-day questions and problems. In this article several examples of research into developing working approaches to welding problems are presented.

In this study, the modeling of mixing by molecular diffusion is a central aspect for transported probability density function (tPDF) methods. In this paper, the newly-proposed shadow position mixing model (SPMM) is examined, using a DNS database for a temporally evolving di-methyl ether slot jet flame. Two methods that invoke different levels of approximation are proposed to extract the shadow displacement (equivalent to shadow position) from the DNS database. An approach for a priori analysis of the mixing-model performance is developed. The shadow displacement is highly correlated with both mixture fraction and velocity, and the peak correlation coefficient of the shadow displacement and mixture fraction is higher than that of the shadow displacement and velocity. This suggests that the composition-space localness is reasonably well enforced by the model, with appropriate choices of model constants. The conditional diffusion of mixture fraction and major species from DNS and from SPMM are then compared, using mixing rates that are derived by matching the mixture fraction scalar dissipation rates. Good qualitative agreement is found, for the prediction of the locations of zero and maximum/minimum conditional diffusion locations for mixture fraction and individual species. Similar comparisons are performed for DNS and the IECM (interaction by exchange with the conditional mean) model. The agreement between SPMM and DNS is better than that between IECM and DNS, in terms of conditional diffusion iso-contour similarities and global normalized residual levels. It is found that a suitable value for the model constant c that controls the mixing frequency can be derived using the local normalized scalar variance, and that the model constant a controls the localness of the model. A higher-Reynolds-number test case is anticipated to be more appropriate to evaluate the mixing models, and stand-alone transported PDF simulations are required to more fully enforce localness and to assess model performance.

Much effort has been made to model hydrogen releases from leaks during potential failures of hydrogen storage systems. A reduced-order jet model can be used to quickly characterize these flows, with low computational cost. Notional nozzle models are often used to avoid modeling the complex shock structures produced by the underexpanded jets by determining an "effective" source to produce the observed downstream trends. In this work, the mean hydrogen concentration fields were measured in a series of subsonic and underexpanded jets using a planar laser Rayleigh scattering system. The experimental data was compared to a reduced order jet model for subsonic flows and a notional nozzle model coupled to the jet model for underexpanded jets. The values of some key model parameters were determined by comparisons with the experimental data. The coupled model was also validated against hydrogen concentrations measurements for 100 and 200 bar hydrogen jets with the predictions agreeing well with data in the literature.

Much effort has been made to model hydrogen releases from leaks during potential failures of hydrogen storage systems. A reduced-order jet model can be used to quickly characterize these flows, with low computational cost. Notional nozzle models are often used to avoid modeling the complex shock structures produced by the underexpanded jets by determining an "effective" source to produce the observed downstream trends. In this work, the mean hydrogen concentration fields were measured in a series of subsonic and underexpanded jets using a planar laser Rayleigh scattering system. The experimental data was compared to a reduced order jet model for subsonic flows and a notional nozzle model coupled to the jet model for underexpanded jets. The values of some key model parameters were determined by comparisons with the experimental data. The coupled model was also validated against hydrogen concentrations measurements for 100 and 200 bar hydrogen jets with the predictions agreeing well with data in the literature.

An associated particle neutron generator is described that employs a negative ion source to produce high neutron flux from a small source size. Negative ions produced in an rf-driven plasma source are extracted through a small aperture to form a beam which bombards a positively biased, high voltage target electrode. Electrons co-extracted with the negative ions are removed by a permanent magnet electron filter. The use of negative ions enables high neutron output (100% atomic ion beam), high quality imaging (small neutron source size), and reliable operation (no high voltage breakdowns). The neutron generator can operate in either pulsed or continuous-wave (cw) mode and has been demonstrated to produce 106 D-D n/s (equivalent to ~108 D-T n/s) from a 1 mm-diameter neutron source size to facilitate high fidelity associated particle imaging.

In previous research, two-pass repeat-geometry synthetic aperture radar (SAR) coherent change detection (CCD) predominantly utilized the sample degree of coherence as a measure of the temporal change occurring between two complex-valued image collects. Previous coherence-based CCD approaches tend to show temporal change when there is none in areas of the image that have a low clutter-to-noise power ratio. Instead of employing the sample coherence magnitude as a change metric, in this paper, we derive a new maximum-likelihood (ML) temporal change estimate—the complex reflectance change detection (CRCD) metric to be used for SAR coherent temporal change detection. The new CRCD estimator is a surprisingly simple expression, easy to implement, and optimal in the ML sense. As a result, this new estimate produces improved results in the coherent pair collects that we have tested.

We develop and demonstrate a method to efficiently use density functional calculations to drive classical dynamics of complex atomic and molecular systems. The method has the potential to scale to systems and time scales unreachable with current ab initio molecular dynamics schemes. It relies on an adapting dataset of independently computed Hellmann–Feynman forces for atomic configurations endowed with a distance metric. The metric on configurations enables fast database lookup and robust interpolation of the stored forces. Here, we discuss mechanisms for the database to adapt to the needs of the evolving dynamics, while maintaining accuracy, and other extensions of the basic algorithm.

The 6th US/German Workshop on Salt Repository Research, Design, and Operation was held in Dresden. Germany on September 7-9, 2015. Over seventy participants helped advance the technical basis for salt disposal of radioactive waste. The number of collaborative efforts continues to grow and to produce useful documentation, as well as to define the state of the art for research areas. These Proceedings are divided into Chapters, and a list of authors is included in the Acknowledgement Section. Also in this document are the Technical Agenda, List of Participants, Biographical Information, Abstracts, and Presentations. Proceedings of all workshops and other pertinent information are posted on websites hosted by Sandia National Laboratories and the Nuclear Energy Agency Salt Club. The US/German workshops provide continuity for long-term research, summarize and publish status of mature areas, and develop appropriate research by consensus in a workshop environment. As before, major areas and findings are highlighted, which constitute topical Chapters in these Proceedings. In total, the scientific breadth is substantial and while not all subject matter is elaborated into chapter format, all presentations and abstracts are published in this document. In the following Proceedings, six selected topics are developed in detail.

We report on a new technique for obtaining off-Hugoniot pressure vs. density data for solid metals compressed to extreme pressure by a magnetically driven liner implosion on the Z-machine (Z) at Sandia National Laboratories. In our experiments, the liner comprises inner and outer metal tubes. The inner tube is composed of a sample material (e.g., Ta and Cu) whose compressed state is to be inferred. The outer tube is composed of Al and serves as the current carrying cathode. Another aluminum liner at much larger radius serves as the anode. A shaped current pulse quasi-isentropically compresses the sample as it implodes. The iterative method used to infer pressure vs. density requires two velocity measurements. Photonic Doppler velocimetry probes measure the implosion velocity of the free (inner) surface of the sample material and the explosion velocity of the anode free (outer) surface. These two velocities are used in conjunction with magnetohydrodynamic simulation and mathematical optimization to obtain the current driving the liner implosion, and to infer pressure and density in the sample through maximum compression. This new equation of state calibration technique is illustrated using a simulated experiment with a Cu sample. Monte Carlo uncertainty quantification of synthetic data establishes convergence criteria for experiments. Results are presented from experiments with Al/Ta, Al/Cu, and Al liners. Symmetric liner implosion with quasi-isentropic compression to peak pressure ∼1000 GPa is achieved in all cases. These experiments exhibit unexpectedly softer behavior above 200 GPa, which we conjecture is related to differences in the actual and modeled properties of aluminum.

The ACRR pulse is pneumatically driven by nitrogen in a system of pipes, valves and hoses up to the connection of the pneumatic system and mechanical linkages of the transient rod (TR). The main components of the TR pneumatic system are the regulator, accumulator, solenoid valve and piston-cylinder assembly. The purpose of this analysis is to analyze the flow of nitrogen through the TR pneumatic system in order to develop a motion profile of the piston during the pulse and be able to predict the pressure distributions inside both the cylinder and accumulators. The predicted pressure distributions will be validated against pressure transducer data, while the motion profile will be compared to proximity switch data. By predicting the motion of the piston, pulse timing will be determined and provided to the engineers/operators for verification. The motion profile will provide an acceleration distribution to be used in Razorback to more accurately predict reactivity insertion into the system.

The article provides information about an upcoming conference from the program chair. The Microscopy Society of America (MSA), the Microanalysis Society (MAS), and the International Metallographic Society (IMS) invite participation in Microscopy & Microanalysis 2016 in Columbus, Ohio, July 24 through July 28, 2016.

Controlling the quantum entanglement between parts of a many-body system is key to unlocking the power of quantum technologies such as quantum computation, high-precision sensing, and the simulation of many-body physics. The spin degrees of freedom of ultracold neutral atoms in their ground electronic state provide a natural platform for such applications thanks to their long coherence times and the ability to control them with magneto-optical fields. However, the creation of strong coherent coupling between spins has been challenging. Here we demonstrate a strong and tunable Rydberg-dressed interaction between spins of individually trapped caesium atoms with energy shifts of order 1 MHz in units of Planck's constant. This interaction leads to a ground-state spin-flip blockade, whereby simultaneous hyperfine spin flips of two atoms are inhibited owing to their mutual interaction. We employ this spin-flip blockade to rapidly produce single-step Bell-state entanglement between two atoms with a fidelity 81(2)%.

Localized stress variation in aluminum nitride (AlN) sputtered on patterned metallization has been monitored through the use of UV micro-Raman spectroscopy. This technique utilizing 325 nm laser excitation allows detection of the AlN E2(high) phonon mode in the presence of metal electrodes beneath the AlN layer with a high spatial resolution of less than 400 nm. The AlN film stress shifted 400 MPa from regions where AlN was deposited over a bottom metal electrode versus silicon dioxide. Across wafer stress variations were also investigated showing that wafer level stress metrology, for example using wafer curvature measurements, introduces large uncertainties for predicting the impact of AlN residual stress on the device performance.

Rebooting Computing (RC) is an effort in the IEEE to rethink future computers. RC started in 2012 by the co-chairs, Elie Track (IEEE Council on Superconductivity) and Tom Conte (Computer Society). RC takes a holistic approach, considering revolutionary as well as evolutionary solutions needed to advance computer technologies. Three summits have been held in 2013 and 2014, discussing different technologies, from emerging devices to user interface, from security to energy efficiency, from neuromorphic to reversible computing. The first part of this paper introduces RC to the design automation community and solicits revolutionary ideas from the community for the directions of future computer research. Energy efficiency is identified as one of the most important challenges in future computer technologies. The importance of energy efficiency spans from miniature embedded sensors to wearable computers, from individual desktops to data centers. To gauge the state of the art, the RC Committee organized the first Low Power Image Recognition Challenge (LPIRC). Each image contains one or multiple objects, among 200 categories. A contestant has to provide a working system that can recognize the objects and report the bounding boxes of the objects. The second part of this paper explains LPIRC and the solutions from the top two winners.

We demonstrate that metal-organic frameworks (MOFs) can catalyze hydrogenolysis of aryl ether bonds under mild conditions. Mg-IRMOF-74(I) and Mg-IRMOF-74(II) are stable under reducing conditions and can cleave phenyl ethers containing β-O-4, α-O-4, and 4-O-5 linkages to the corresponding hydrocarbons and phenols. Reaction occurs at 10 bar H2 and 120 °C without added base. DFT-optimized structures and charge transfer analysis suggest that the MOF orients the substrate near Mg2+ ions on the pore walls. Ti and Ni doping further increase conversions to as high as 82% with 96% selectivity for hydrogenolysis versus ring hydrogenation. Repeated cycling induces no loss of activity, making this a promising route for mild aryl-ether bond scission.

Particle image velocimetry measurements have been conducted for a Mach 0.8 flow over a wall-mounted hemisphere with a strongly separated wake. The shock foot was found to typically sit just forward of the apex of the hemisphere and move within a range of about ±10 deg. Conditional averages based upon the shock foot location show that the separation shock is positioned upstream along the hemisphere surface when reverse velocities in the recirculation region are strong and is located downstream when they are weaker. The recirculation region appears smaller when the shock is located farther downstream. No correlation was detected of the incoming boundary layer with the shock position nor with the wake recirculation velocities. These observations are consistent with recent studies concluding that, for large, strong separation regions, the dominant mechanism is the instability of the separated flow rather than a direct influence of the incoming boundary layer.

A thermal rectifier that utilizes thermal expansion to directionally control interfacial conductance between two contacting surfaces is presented. The device consists of two thermal reservoirs contacting a beam with one rough and one smooth end. When the temperature of reservoir in contact with the smooth surface is raised, a similar temperature rise will occur in the beam, causing it to expand, thus increasing the contact pressure at the rough interface and reducing the interfacial contact resistance. However, if the temperature of the reservoir in contact with the rough interface is raised, the large contact resistance will prevent a similar temperature rise in the beam. As a result, the contact pressure will be marginally affected and the contact resistance will not change appreciably. Owing to the decreased contact resistance of the first scenario compared to the second, thermal rectification occurs. A parametric analysis is used to determine optimal device parameters including surface roughness, contact pressure, and device length. Modeling predicts that rectification factors greater than 2 are possible at thermal biases as small as 3 K. Additionally, thin surface coatings are discussed as a method to control the temperature bias at which maximum rectification occurs.

Abstract not provided.

Polymer foam encapsulants provide mechanical, electrical, and thermal isolation in engineered systems. It can be advantageous to surround objects of interest, such as electronics, with foams in a hermetically sealed container in order to protect them from hostile environments or from accidents such as fire. In fire environments, gas pressure from thermal decomposition of foams can cause mechanical failure of sealed systems. In this work, a detailed uncertainty quantification study of polymeric methylene diisocyanate (PMDI)-polyether-polyol based polyurethane foam is presented and compared to experimental results to assess the validity of a 3-D finite element model of the heat transfer and degradation processes. In this series of experiments, 320 kg/m3 PMDI foam in a 0.2 L sealed steel container is heated to 1,073 K at a rate of 150 K/min. The experiment ends when the can breaches due to the buildup of pressure. The temperature at key location is monitored as well as the internal pressure of the can. Both experimental uncertainty and computational uncertainty are examined and compared. The mean value method (MV) and Latin hypercube sampling (LHS) approach are used to propagate the uncertainty through the model. The results of the both the MV method and the LHS approach show that while the model generally can predict the temperature at given locations in the system, it is less successful at predicting the pressure response. Also, these two approaches for propagating uncertainty agree with each other, the importance of each input parameter on the simulation results is also investigated, showing that for the temperature response the conductivity of the steel container and the effective conductivity of the foam, are the most important parameters. For the pressure response, the activation energy, effective conductivity, and specific heat are most important. The comparison to experiments and the identification of the drivers of uncertainty allow for targeted development of the computational model and for definition of the experiments necessary to improve accuracy.

Implications of conducting hardware-in-the-loop testing of a specific hydrokinetic turbine on controllable motor-generator sets or electromechanical emulation machines (EEMs) are explored. The emulator control dynamic equations are presented, methods for scaling turbine parameters are developed and evaluated, and experimental results are presented from three EEMs programmed to emulate the same vertical-axis fixed-pitch turbine. Although hardware platforms and control implementations varied, results show that each EEM is successful in emulating the turbine model at different power levels, thus demonstrating the general feasibility of the approach. However, performance of motor control under torque command, current command, or speed command differed. In a demonstration of the intended use of an EEM for evaluating a hydrokinetic turbine implementation, a power takeoff controller tracks the maximum power-point of the turbine in response to turbulence. Utilizing realistic inflow conditions and control laws, the emulator dynamic speed response is shown to agree well at low frequencies with numerical simulation but to deviate at high frequencies.

Verification of tightly coupled multiphysics computational codes is generally significantly more difficult than verification of single-physics codes. The case of coupled heat conduction and thermal radiation in an enclosure is considered, and it is extended to a manufactured solution verification test for enclosure radiation to a fully two-dimensional coupled problem with conduction and thermal radiation. Convergence results are shown using a production thermal analysis code. Convergence rates are optimal with a pairwise view-factor calculation algorithm.

We report on the thermodynamic properties of binary compound mixtures of model groups II-VI semiconductors. We use the recently introduced Stillinger-Weber Hamiltonian to model binary mixtures of CdTe and CdSe. We use molecular dynamics simulations to calculate the volume and enthalpy of mixing as a function of mole fraction. The lattice parameter of the mixture closely follows Vegard's law: a linear relation. This implies that the excess volume is a cubic function of mole fraction. A connection is made with hard sphere models of mixed fcc and zincblende structures. The potential energy exhibits a positive deviation from ideal soluton behaviour; the excess enthalpy is nearly independent of temperatures studied (300 and 533 K) and is well described by a simple cubic function of the mole fraction. Using a regular solution approach (combining non-ideal behaviour for the enthalpy with ideal solution behaviour for the entropy of mixing), we arrive at the Gibbs free energy of the mixture. The Gibbs free energy results indicate that the CdTe and CdSe mixtures exhibit phase separation. The upper consolute temperature is found to be 335 K. Finally, we provide the surface energy as a function of composition. It roughly follows ideal solution theory, but with a negative deviation (negative excess surface energy). This indicates that alloying increases the stability, even for nano-particles.

The breakup of liquids due to aerodynamic forces has been widely studied. However, the literature contains limited quantified data on secondary droplet sizes, particularly as a function of time. Here, a column of liquid water is subjected to a step change in relative gas velocity using a shock tube. A unique digital in-line holography (DIH) configuration is proposed which quantifies the secondary droplets sizes, three-dimensional position, and three-component velocities at 100 kHz. Results quantify the detailed evolution of the characteristic mean diameters and droplet size-velocity correlations as a function of distance downstream from the initial location of the water column. Accuracy of the measurements is confirmed through mass balance. These data give unprecedented detail on the breakup process which will be useful for improved model development and validation.

In this discussion paper, we explore different ways to assess the value of verification and validation (V&V) of engineering models. We first present a literature review on the value of V&V and then use value chains and decision trees to show how value can be assessed from a decision maker’s perspective. In this context, the value is what the decision maker is willing to pay for V&V analysis with the understanding that the V&V results are uncertain. The 2014 Sandia V&V Challenge Workshop is used to illustrate these ideas.

High-speed, time-resolved particle image velocimetry with a pulse-burst laser was used to measure the gas-phase velocity upstream and downstream of a shock wave-particle curtain interaction at three shock Mach numbers (1.19, 1.40, and 1.45), at a sampling rate of 37.5 kHz. The particle curtain, formed from free-falling soda-lime particles with diameters ranging from 300 - 355 μm, had a streamwise thickness of 3.5 mm and volume fraction of 9% at mid-height. Following impingement by a shock wave, a pressure difference was created between the upstream/downstream sides of the curtain, which accelerated flow through the curtain. Jetting of flow through the curtain was observed downstream once deformation of the curtain began, demonstrating a long-term unsteady effect. Using a control volume approach, the unsteady drag on the curtain was determined from velocity and pressure data. Initially, the pressure difference between the upstream and downstream sides of the curtain was the largest contributor to the total drag. The data suggests, however, that as time increases, the change in momentum flux could become the dominant component as the pressure difference decreases.

Pulse-burst particle image velocimetry (PIV) has been used to acquire time-resolved data at 37.5 kHz of the flow over a finite-width rectangular cavity at Mach 0.6, 0.8, and 0.94. Power spectra of the PIV data reveal four resonance modes that match the frequencies detected simultaneously using high-frequency wall pressure sensors. Velocity resonances exhibit spatial dependence in which the lowest-frequency acoustic mode is active within the recirculation region whereas the three higher modes are concentrated within the shear layer. Spatio-temporal cross-correlations were calculated from velocity data first bandpass filtered for specific resonance frequencies. The low-frequency acoustic mode shows properties of a standing wave without spatial correlation. Higher resonance modes are associated with alternating coherent structures whose size and spacing decrease for higher resonance modes and increase as structures convect downstream. The convection velocity appears identical for the high-frequency resonance modes, but it too increases with downstream distance. This is in contrast to the well-known Rossiter equation, which assumes a convection velocity constant in space.

Time-resolved particle image velocimetry (PIV) using a pulse-burst laser has been acquired of a supersonic jet issuing into a Mach 0.8 crossflow. Simultaneously, the final pulse pair in each burst has been imaged using conventional PIV cameras to produce an independent two-component measurement and two stereoscopic measurements. Each measurement depicts generally similar flowfield features with vorticity contours marking turbulent eddies at corresponding locations. Probability density functions of the velocity fluctuations are essentially indistinguishable but the precision uncertainty estimated using correlation statistics shows that the pulse-burst PIV data have notably greater uncertainty than the three conventional measurements. This occurs due to greater noise in the cameras and a smaller size for the final iteration of the interrogation window. A small degree of peak locking is observed in the aggregate of the pulse-burst PIV data set. However, some of the individual vector fields show peak locking to non-integer pixel values as a result of real physical effects in the flow. Even if peak locking results entirely from measurement bias, the effect occurs at too low a level to anticipate a significant effect on data analysis.

Stereoscopic particle image velocimetry was used to experimentally measure the recirculating flow within finite-span cavities of varying complex geometry at a freestream Mach number of 0.8. Volumetric measurements were made to investigate the side wall influences by scanning a laser sheet across the cavity. Each of the geometries could be classied as an open-cavity, based on L/D. The addition of ramps altered the recirculation zone within the cavity, causing it to move along the streamwise direction. Within the simple rectangular cavity, a system of counter-rotating streamwise vortices formed due to spillage from along the side wall, which caused the mixing layer to develop a steady spanwise waviness. The ramped complex geometry, due to the presence of leading edge and side ramps, appeared to suppress the formation of streamwise vorticity associated with side wall spillage, resulting in a much more two-dimensional mixing layer.

A previous study in the UK demonstrated that vibration response on a scaled-down model of a missile structure in a wind tunnel could be replicated in a laboratory setting with multiple shakers using an approach dubbed as impedance matching. Here we demonstrate on a full scale industrial structure that the random vibration induced from a laboratory acoustic environment can be nearly replicated at 37 internal accelerometers using six shakers. The voltage input to the shaker amplifiers is calculated using a regularized inverse of the square of the amplitude of the frequency response function matrix and the power spectral density responses of the 37 internal accelerometers. No cross power spectral density responses are utilized. The structure has hundreds of modes and the simulation is performed out to 4000 Hz.