Publications

Abstract not provided.

In this report, we investigate how manufacturing conditions result in the warpage of moderate density PMDI polyurethane foam (12-50 lb/ft³) when they are released from a mold. We have developed a multiphysics modeling framework to simulate the manufacturing process including resin injection, foaming and mold filling, gelation of the matrix, elevated cure, vitrification, cool down, and demolding. We have implemented this framework within the Sierra Mechanics Finite Element Code Suite. We couple Aria for flow, energy conservation, and foaming/curing kinetics with Adagio for the nonlinear viscoelastic solid response in a multi-staged simulation process flow. We calibrate a model for the PMDI-10S (10 lb/ft³ free rise foam) through a suite of characterization data presented here to calibrate the solid cure behavior of the foam. The model is then used and compared to a benchmark experiment, the manufacturing and warpage over 1 year of a 10 cm by 10 cm by 2.5 cm foam "staple". This component features both slender and thick regions that warp considerably differently over time. Qualitative agreement between the model and the experiment is achieved but quantitative accuracy is not.

Abstract not provided.

Abstract not provided.

Sandia National Laboratories has tested and evaluated a new digitizer, the Affinity, manufactured by Guralp Systems Ltd. These digitizers are used to record sensor output for seismic and infrasound monitoring applications. The purpose of the digitizer evaluation was to measure the performance characteristics in such areas as power consumption, input impedance, sensitivity, full scale, self-noise, dynamic range, system noise, response, passband, and timing. The Affinity digitizer is Guralp's latest release in their digitizer product line. The Affinity is available with either 4 or 8 channels at 24 bit resolution. In addition to the 24 bit channels, 16 multiplexed low resolution channels are provided. Other features include the means to accept multiple types of timing sources (e.g. GPS, NTP and PTP) and a web page interface for command and control of the unit.

Sandia National Laboratories has tested and evaluated a new digitizer, the Q330HR, manufactured by Quanterra. These digitizers are used to record sensor output for seismic and infrasound monitoring applications. The purpose of the digitizer evaluation was to measure the performance characteristics in such areas as power consumption, input impedance, sensitivity, full scale, self-noise, dynamic range, system noise, response, passband, and timing. The Q330HR is Quanterra’s improved Q330 datalogger with a 26 bits of resolution on channels 1-3 and a 24 bits of resolution on channels 4-6 (26 bit is optional). The Quanterra Q330HR is being evaluated for potential use U.S. Air Force seismic monitoring systems as part of their Next Generation Qualification effort.

Sandia National Laboratories has tested and evaluated a new digitizer, the Centaur, manufactured by Nanometrics. These digitizers are used to record sensor output for seismic and infrasound monitoring applications. The purpose of the digitizer evaluation was to measure the performance characteristics in such areas as power consumption, input impedance, sensitivity, full scale, self-noise, dynamic range, system noise, response, passband, and timing. The Centaur digitizer is Nanometrics’ replacement for their Taurus digitizer and marks Nanometrics first 6 channel, 24 bit resolution system. Other improvements include LED status indicators on the top of the unit, providing basic status of the core systems of a seismic station (e.g. timing, sensor SOH, storage, etc), an optional wifi system allowing password protected access to the unit and a web interface for monitoring and configuration of the unit. The Nanometrics Centaur is being evaluated for potential use U.S. Air Force seismic monitoring systems as part of their Next Generation Qualification effort.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The levelized cost of energy for an offshore wind plant consisting of floating vertical-axis wind turbines is studied in this report. A 5 MW Darrieus vertical-axis wind turbine rotor is used as the study turbine as this architecture was determined to have the greatest ability to reduce the system cost. The rotor structural design was used with blade manufacturing cost model studies to estimate its cost. A two-bladed, carbon fiber rotor was selected in this analysis since the lower topside mass resulted in a reduction of the platform costs which exceeded the increased rotor cost. A direct-drive, medium efficiency drivetrain was designed which represents 25% of the costs and 45% of the mass of the combined rotor/drivetrain system. A direct-drive, permanent magnet generator drivetrain was selected due to the improved reliability of this type of system, while the cost was not significantly higher than for geared drivetrains. A platform was designed by first identifying the optimal architecture for the vertical-axis wind turbine at a water depth of 150 m. A survey was performed of floating platform types, and six characteristic designs were analyzed which span the range of stability mechanisms available to floating systems. A multi-cellular tension-leg platform was identified as the lowest cost platform which additionally provided some interesting performance benefits. The small motions of the tension-leg platform benefit the system energy capture while limiting inertial loads placed on the rotor’s tower and blades. A final design was produced for the multi-cellular tension-leg platform considering operational fatigue, storm wind and wave conditions, and tow-out design cases. The driving design load was stability during tow-out while ballasting the platform. System levelized cost of energy was calculated, including operational expenses and balance of system costs estimated for the wind plant. Opportunities for reduction in the component costs are predicted and used to make projections of the system levelized cost of energy for future developments. The opportunities and challenges for floating vertical-axis wind turbines are identified by the system design and levelized cost of energy analysis.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The moisture absorption behavior of two fiber reinforced composite materials was evaluated in a unidirectional manner The flat materials were exposed to varying humidity and temperature conditions inside of an environmental chamber in order to determine their effective moisture equilibrium (M m ) and moisture absorption rate (D z ). Two-ply (thin) and four-ply (thick) materials were utilized to obtain M,,, and Dz, respectively. The results obtained from laboratory work were then compared to modeling data to better understand the material properties. Predictions capabilities were built to forecast the maximum moisture content, time required for saturation, and the moisture content at any given humidity and temperature. A case study was included to demonstrate this capability. Also of interest were cubed samples to investigate directionality preferences in water immersion studies. Several coatings were evaluated for their water permeation properties. Further dissemination authorized to the Department of Energy and DOE contractors only; other requests shall be approved by the originating facility or higher DOE programmatic authority.

This report compares ATLOG modeling results for the response of a finite-length dissipative buried conductor interacting with a conducting ground to a measurement taken November 2016 at the High-Energy Radiation Megavolt Electron Source (HERMES) facility. We use the ATLOG frequency-domain method based on transmission line theory. Estimates of the impedance per unit length and admittance per unit length for a cable laying in a PVC pipe embedded in a concrete block are reported. Current wave shapes from both a single conductor and composite differential mode and antenna mode arrangements are close to those observed in the experiments.

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.

Narratives about water resources have evolved, transitioning from a sole focus on physical and biological dimensions to incorporate social dynamics Recently, the importance of understanding the visibility of water resources through media coverage has gained attention. This study leverages recent advancements in natural language processing (NLP) methods to characterize and understand patterns in water narratives, specifically in 4 local newspapers in Utah and Georgia. Analysis of the corpus identified coherent topics on a variety of water resources issues, including weather and pollution. Closer inspection of the topics revealed temporal and spatial variations in coverage, with a topic on hurricanes exhibiting cyclical patterns whereas a topic on tribal issues showed coverage predominantly in the western newspapers. We also analyzed the dataset for sentiments, identifying similar categories of words on trust and fear emerging in the narratives across newspaper sources. An analysis of novelty, transience, and resonance using Kullback-Leibler Divergence techniques revealed that topics with high novelty generally contained high transience and marginally high resonance over time. Although additional analysis needs to be conducted, the methods explored in this analysis demonstrate the potential of NLP methods to characterize water narratives in media coverage.

In this work, various material models were studied for their ability to simulate puncture in a thin aluminum 7075-T651 plate due to low-velocity probe impact. Material models were generated by mixing and matching various work hardening laws with different failure criteria, and several hybrid material models were investigated. Finite element simulations of aluminum impact-response, based on each material model, were employed to predict the energy required for puncture and final plate tear-out geometry. Probes of different size and shape were used to impose various loading regimes, and numerical predictions were compared to experimental results from a previous study. It was found that no single combination of hardening and failure laws yielded universally accurate data, but that several material models could be used more reliably than others. Further, the importance of obtaining unique parameter-sets for work-hardening and failure criteria was illustrated.

Abstract not provided.

Abstract not provided.

Abstract not provided.

CSP Dish systems are parabolic mirror structures that track the sun in two axes, focusing the Direct Normal insolation (DNI) to a point or spot on a boom or tripod mounted to the tracking dish structure. This focused light is typically utilized in-situ to operate a heat engine, such as a Stirling cycle, Brayton cycle, or Rankine cycle engine to make electricity. Other dish systems have been used to generate steam for a centralized engine fed by multiple dishes, or to operate thermochemical processes for industrial use, storage, or creation of fuels. Because the dish is always pointing at the sun, a well-designed dish system has a very high concentration ratio, allowing the generation of high temperatures, leading to high thermodynamic efficiencies.

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.

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.

Schedule Management Optimization (SMO) is a tool for automatically generating a schedule of project tasks. Project scheduling is traditionally achieved with the use of commercial project management software or case-specific optimization formulations. Commercial software packages are useful tools for managing and visualizing copious amounts of project task data. However, their ability to automatically generate optimized schedules is limited. Furthermore, there are many real-world constraints and decision variables that commercial packages ignore. Case-specific optimization formulations effectively identify schedules that optimize one or more objectives for a specific problem, but they are unable to handle a diverse selection of scheduling problems. SMO enables practitioners to generate optimal project schedules automatically while considering a broad range of real-world problem characteristics. SMO has been designed to handle some of the most difficult scheduling problems – those with resource constraints, multiple objectives, multiple inventories, and diverse ways of performing tasks. This report contains descriptions of the SMO modeling concepts and explains how they map to real-world scheduling considerations.

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.

The atmospheric dispersion of contaminants in the wake of a large urban structure is a challenging fluid mechanics problem of interest to the scientific and engineering communities. Magnetic Resonance Velocimetry (MRV) is a relatively new technique that leverages diagnostic equipment used primarily by the medical field to make 3D engineering measurements of flow and contaminant dispersal. SIERRA/Fuego, a computational fluid dynamics (CFD) code at Sandia National Labs is employed to make detailed comparisons to the dataset to evaluate the quantitative and qualitative accuracy of the model. The comparison exercise shows good comparison between model and experimental results, with the wake region downstream of the tall building presenting the most significant challenge to the quantitative accuracy of the model. Model uncertainties are assessed through parametric variations. Some observations are made in relation to the future utility of MDV and CFD, and some productive follow-on activities are suggested that can help mature the science of flow modeling and experimental testing.

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.

Full waveform inversion allows the seismologist to utilize an entire waveform and all the information it contains to help image the 3-D structure of the interior of the earth. This report summarizes the basic theory that has been developed in full waveform seismic inversion, primarily related to computation of sensitivity kernels. It then describes the implementation of this theory using Sandia Geophysics Department's Parelasti code, a 3-D full waveform elastic simulation algorithm. Finally, the code is validated using synthetics from simple homogeneous elastic earth models.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Many earth materials and minerals are seismically anisotropic; however, due to the weakness of anisotropy and for simplicity, the earth is often approximated as an isotropic medium. Specific circumstances, such as in shales, tectonic fabrics, or oriented fractures, for example, require the use of anisotropic simulations in order to accurately model the earth. This report details the development of a new massively parallel 3-D full seismic waveform simulation algorithm within the principle coordinate system of an orthorhombic material, which is a specific form of anisotropy common in layered, fractured media. The theory and implementation of Pararhombi is described along with verification of the code against other solutions.

We invert far field infrasound data for the equivalent seismo-acoustic time domain moment tensor to assess the effects of variable atmospheric models as well as to quantify the relative contributions of two presumed source phenomena. The infrasound data was produced by a series of underground chemical explosions that were conducted during the Source Physics Experiment, (SPE) which was originally designed to study explosion-generated seismo-acoustic signal phenomena. The goal of the work presented herein is two-fold: the first goal is to investigate the sensitivity of the estimated time domain moment tensors to variability of the estimated atmospheric model. The second goal is to determine the relative contribution of two possible source mechanisms to the observed infrasonic wave field. Rather than using actual atmospheric observations to estimate the necessary atmospheric Green's functions, we build a series of atmospheric models that rely on publicly available, regional atmospheric observations and the assumption that the acoustic energy results from a linear combination of an underground isotropic explosion and surface spall. The atmospheric observations are summarized and interpolated onto a 3D grid to produce a model of sound speed at the time of the experiment. For each of four SPE acoustic datasets that we invert, we produced a suite of three atmospheric models, based on ten years of regional meteorological observations: an average model, which averages the atmospheric conditions for ten years prior to each SPE event, as well as two extrema models. We find that the inversion yields relatively repeatable results for the estimated spall source. Conversely, the estimated isotropic explosion source is highly variable. This suggests that the majority of the observed acoustic energy is produced by the spall source and/or our modeling of the elastic energy propagation, and it's subsequent conversion to acoustic energy via linear elastic-to-acoustic coupling at the free surface, is too simplistic.

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.

Due to the weight of overburden and tectonic forces, the solid earth is subject to an ambient stress state. This stress state is quasi-static in that it is generally in a state of equilibrium. Typically, seismology assumes this ambient stress field has a negligible effect on wave propagation. However, two basic theories have been put forward to describe the effects of ambient stress on wave propagation. Dahlen and Tromp (2002) expound a theory based on perturbation analysis that largely supports the traditional seismological view that ambient stress is negligible for wave propagation. The second theory, espoused by Korneev and Glubokovskikh (2013) and supported by some experimental work, states that perturbation analysis is inappropriate since the elastic modulus is very sensitive to the ambient stress states. This brief report reformulates the equations given by Korneev and Glubokovskikh (2013) into a more compact form that makes it amenable to statement in terms of a pre-stress form of Hooke's Law. Furthermore, this report demonstrates the symmetries of the pre-stress modulus tensor and discusses the reciprocity relationship implied by the symmetry conditions.

Due to the weight of overburden and tectonic forces, the solid earth is subject to an ambient stress state. This stress state is quasi-static in that it is generally in a state of equilibrium. Typically, seismology assumes this ambient stress field has a negligible effect on wave propagation. However, two basic theories have been put forward to describe the effects of ambient stress on wave propagation. Dahlen and Tromp (2002) expound a theory based on perturbation analysis that largely supports the traditional seismological view that ambient stress is negligible for wave propagation. The second theory, espoused by Korneev and Glubokovskikh (2013) and supported by some experimental work, states that perturbation analysis is inappropriate since the elastic modulus is very sensitive to the ambient stress states. This brief report reformulates the equations given by Korneev and Glubokovskikh (2013) into a more compact form that makes it amenable to statement in terms of a pre-stress form of Hooke's Law. Furthermore, this report demonstrates the symmetries of the pre-stress modulus tensor and discusses the reciprocity relationship implied by the symmetry conditions.

Waves propagating through natural materials such as ocean water encounter spatial variations in material properties that cannot easily be predicted or known in advance. Deterministic wave simulation algorithms must assume that all properties throughout the model space are precisely known. However, a stochastic wave simulation tool can parameterize the material as a stochastic medium with a certain probability distribution and correlation length. This report documents the addition of spatial stochastic variability into Paracousti-UQ, Sandia Geophysics Department's 3-D full waveform acoustic algorithm within stochastic media. The ability of the code to replicate Monte Carlo solutions in 1-D spatially variable media is also evaluated.

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.

This report is an outcome of the ASC CSSE Level 2 Milestone 6362: Analysis of Re- silient Asynchronous Many-Task (AMT) Programming Model. It comprises a summary and in-depth analysis of resilience schemes adapted to the AMT programming model. Herein, performance trade-offs of a resilient-AMT prograrnming model are assessed through two ap- proaches: (1) an analytical model realized by discrete event simulations and (2) empirical evaluation of benchmark programs representing regular and irregular workloads of explicit partial differential equation solvers. As part of this effort, an AMT execution simulator and a prototype resilient-AMT programming framework have been developed. The former permits us to hypothesize the performance behavior of a resilient-AMT model, and has undergone a verification and validation (V&V) process. The latter allows empirical evaluation of the perfor- mance of resilience schemes under emulated program failures and enabled the aforementioned V&V process. The outcome indicates that (1) resilience techniques implemented within an AMT framework allow efficient and scalable recovery under frequent failures, that (2) the abstraction of task and data instances in the AMT programming model enables readily us- able Application Program Interfaces (APIs) for resilience, and that (3) this abstraction enables predicting the performance of resilient-AMT applications with a simple simulation infrastruc- ture. This outcome will provide guidance for the design of the AMT programming model and runtime systems, user-level resilience support, and application development for ASC's next generation platforms (NGPs).

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document consolidates the work performed by Sandia National Laboratories and the US Nuclear Regulatory Commission in participation of Program IRIS: “Improving the Robustness of the Assessment Methodologies for Structures Impacted by Missiles”. Three round-robin benchmark exercises on improving the robustness of the assessment of structures impacted by large missiles at medium to high velocities were organized by either the IAGE Subgroup on Ageing of Concrete Structures of the Organization for Economic Co-operation and Development Nuclear Energy Agency (NEA) or Électricité de France (EDF). The objectives of the exercises were to develop guidance for conducting impact analyses including issues related to computer codes, modeling approaches, and analysis techniques. The full project was comprised of three phases: Phase I, impact of walls; Phase II, impact of larger structures; and Phase III, transmission of shock and vibration to internal components.

Abstract not provided.

Abstract not provided.

This report provides details of the algorithms in the Bloodhound package for infrasound data analysis. The report provides a detailed description of the algorithms, general instructions on tuning Bloodhound for different signal types, and a complete listing of all input parameters and the complete output schema. Several Jupyter notebooks are provided with the distribution for illustrating how to use Bloodhound for different workflows.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report presents computational analyses that simulate the structural response of caverns at the Strategic Petroleum Reserve Bryan Mound site. The cavern field comprises 20 caverns. Five caverns (1, 2, 4, and 5; 3 was later plugged and abandoned) were acquired from industry and have unusual shapes and a history dating back to 1946. The other 16 caverns (101-116) were leached according to SPR standards in the mid-1980s and have tall cylindrical shapes. The history of the caverns and their shapes are simulated in a 3-D geomechanics model of the site that predicts deformations, strains, and stresses. Historical wellhead pressures are used to calculate cavern pressures up through July 2016. Because of the extent of heterogeneous creep behavior observed throughout the Bryan Mound site, a set of cavern-specific creep coefficients was developed to produce better matches with measured cavern closure and surface subsidence. For this new implementation of the model, there are two significant advances: the use of the multimechanism deformation (M-D) salt creep model to evaluate both steady-state and transient salt creep; and the creation of finite element mesh geometries for the caverns that nearly exactly match the geometries obtained through sonar measurements. The results of the finite element model are interpreted to provide information on the current and future status of subsidence, well integrity, cavern stability, and drawdown availability.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Analysis of quartz sandstones shows that grain-scale crushing (fracture and rearrangement) and associated sealing of fractures contribute significantly to consolidation. The crushing strength (P*) for granular material is defined by laboratory experiments conducted at strain rates of 10−4 to 10−5 s−1 and room temperature. Based on experiments, many sandstones would require burial depths in excess of the actual maximum burial depth to create observed microstructure and density. We use experiments and soil mechanics principles to determine rate laws for brittle consolidation of fine-grained quartz sand to better estimate in situ failure conditions of porous geomaterials. Experiments were conducted on St. Peter sand utilizing different isostatic consolidation and creep load paths at temperatures to 200 °C and at strain rates of 10−4 to 10−10 s−1. Experiment results are consistent with observed rate dependence of consolidation in soils, and P* for sand can be identified by the change in the dependence of consolidation rate with stress, allowing the extrapolation of P* determined in the laboratory to geologic rates and temperatures. Additionally, normalized P* values can be described by a polynomial function to quantify temperature, stress, and strain-rate relationships for the consolidation of porous geomaterials by subcritical cracking. At geologic loading rates, P* for fine-grained quartz sand is achieved within ~3-km burial depth, and thus, shear-enhanced compaction under nonisostatic stress can occur at even shallower depths. These results demonstrate that time and temperature effects must be considered for predicting the brittle consolidation of sediments in depositional basins, petroleum reservoirs, and engineering applications.

Will quantum computation become an important milestone in human progress? Passionate advocates and equally passionate skeptics abound. IEEE already provides useful, neutral forums for state-of-the-art science and engineering knowledge as well as practical benchmarks for quantum computation evaluation. But could the organization do more.

Recent advances in polymer synthesis have allowed remarkable control over chain microstructure and conformation. Capitalizing on such developments, here we create well-controlled chain folding in sulfonated polyethylene, leading to highly uniform hydrated acid layers of subnanometre thickness with high proton conductivity. The linear polyethylene contains sulfonic acid groups pendant to precisely every twenty-first carbon atom that induce tight chain folds to form the hydrated layers, while the methylene segments crystallize. The proton conductivity is on par with Nafion 117, the benchmark for fuel cell membranes. We demonstrate that well-controlled hairpin chain folding can be utilized for proton conductivity within a crystalline polymer structure, and we project that this structure could be adapted for ion transport. This layered polyethylene-based structure is an innovative and versatile design paradigm for functional polymer membranes, opening doors to efficient and selective transport of other ions and small molecules on appropriate selection of functional groups.

This monthly report is intended to communicate the status of North Slope ARM facilities managed by Sandia National Labs.

Measurements of energy balance components (energy intake, energy expenditure, changes in energy stores) are often plagued with measurement error. Doubly-labeled water can measure energy intake (EI) with negligible error, but is expensive and cumbersome. An alternative approach that is gaining popularity is to use the energy balance principle, by measuring energy expenditure (EE) and change in energy stores (ES) and then back-calculate EI. Gold standard methods for EE and ES exist and are known to give accurate measurements, albeit at a high cost. We propose a joint statistical model to assess the measurement error in cheaper, non-intrusive measures of EE and ES. We let the unknown true EE and ES for individuals be latent variables, and model them using a bivariate distribution. We try both a bivariate Normal as well as a Dirichlet Process Mixture Model, and compare the results via simulation. Our approach, is the first to account for the dependencies that exist in individuals’ daily EE and ES. We employ semiparametric regression with free knot splines for measurements with error, and linear components for error free covariates. We adopt a Bayesian approach to estimation and inference and use Reversible Jump Markov Chain Monte Carlo to generate draws from the posterior distribution. Based on the semipar-ameteric regression, we develop a calibration equation that adjusts a cheaper, less reliable estimate, closer to the true value. Along with this calibrated value, our method also gives credible intervals to assess uncertainty. A simulation study shows our calibration helps produce a more accurate estimate. Our approach compares favorably in terms of prediction to other commonly used models.

We have demonstrated a laboratory based, high bandwidth atom interferometer instrument and have performed an incipient gravity measurement with a fractional statistical uncertainty of ag/g = 4.4 x 10^-6 where g is the acceleration due to gravity. We have designed, constructed, and optimised numerous laser systems for this purpose, most notably a powerful Raman laser (12 W at 780 nm) which will allow large momentum transfer and large detuning interferometry to be carried out at high bandwidth for the first time. This experiment is a general purpose test bed for exploring the fundamental limitations of atom interferometer techniques.

There are currently 2,462 dual-purpose canisters (DPCs) containing spent nuclear fuel (SNF) across the United States. Repackaging DPCs into specialized disposal canisters could be financially and operationally costly with additional radiological, operational safety, and management risks. There are several approaches to facilitate direct disposal of DPCs and demonstrate acceptable repository performance. A promising approach is to fill the void space within the DPCs with a material that would significantly limit the potential for criticality through limiting moderation and/or the addition of neutron absorbers in the interstitial spaces within the fuel assemblies and baskets. An acceptable filler would demonstrably show that the probability of criticality in DPCs during the disposal period of interest to be below the probability threshold for inclusion in repository performance assessment. Based on previous work conducted by domestic and international organizations, two approaches were identified as potentially viable for introduction of fillers into DPCs as liquids that would eventually solidify: (1) molten metal fillers introduced at higher temperatures, and (2) resins or cement slurries that solidify at lower temperatures.

This report documents the completion of milestone STPM12-4 Kokkos Training Bootcamp. The goal of this milestone was to hold a combined tutorial and hackathon bootcamp event for the Kokkos community and prospective users. The Kokkos Bootcamp event was held on-site at Oak Ridge National Lab from July 24 — July 27, 2018. There were over 40 registered participants from 12 institutions, including 7 Kokkos project staff from SNL, LANL, and ORNL. The event consisted of a roughly a two-day tutorial session including hands exercises, followed by 1.5 days of intensive porting work on codes that the participants brought explore, port, and optimize the use of Kokkos with the help of Kokkos project experts.

Acoustic waves with a wide range of frequencies are generated by lightning strokes during thunderstorms, including infrasonic waves (0.1 to 20 Hz). The source mechanism for these low-frequency acoustic waves is still debated, and studies have so far been limited to ground-based instruments. Here we report the first confirmed detection of lightning-generated infrasound with acoustic instruments suspended at stratospheric altitudes using a free-flying balloon. We observe high-amplitude signals generated by lightning strokes located within 100 km of the balloon as it flew over the Tasman Sea on 17 May 2016. The signals share many characteristics with waveforms recorded previously by ground-based instruments near thunderstorms. The ability to measure lightning activity with high-altitude infrasound instruments has demonstrated the potential for using these platforms to image the full acoustic wavefield in the atmosphere. Furthermore, it validates the use of these platforms for recording and characterizing infrasonic sources located beyond the detection range of ground-based instruments.

We demonstrate the ultrahigh extinction operation of a silicon photonic (SiP) amplitude modulator (AM) employing a cascaded Mach-Zehnder interferometer. By carrying out optimization sweeps without significantly degrading the extinction, the SiP AM is robust to environment changes and maintained >52 dB extinction for >6 hrs.

The purpose of this study was to first assess the sensitivity of the parameters of elasticplastic material models with anisotropic yield to choices in the calibration procedure. Two models were considered: Hill's 1948 and Barlat's Y1d2004-18p. Subsequently, it was shown that calibration choices can have an effect on the values of the stress and strain at the ultimate point in uniaxial specimen responses. Finally, the calibrated Barlat model was able to reasonably reproduce the load-deflection and strain fields of a validation specimen that experienced multiaxial states of stress. Overall, it was found that the Barlat model resulted in a closer fit to the measurements and that the parameters of the calibration procedure should be varied to assess the sensitivity of the results.

Here, we provide a demonstration that gas-kinetic methods incorporating molecular chaos can simulate the sustained turbulence that occurs in wall-bounded turbulent shear flows. The direct simulation Monte Carlo method, a gas-kinetic molecular method that enforces molecular chaos for gas-molecule collisions, is used to simulate the minimal Couette flow at Re = 500 . The resulting law of the wall, the average wall shear stress, the average kinetic energy, and the continually regenerating coherent structures all agree closely with corresponding results from direct numerical simulation of the Navier-Stokes equations. Finally, these results indicate that molecular chaos for collisions in gas-kinetic methods does not prevent development of molecular-scale long-range correlations required to form hydrodynamic-scale turbulent coherent structures.