The goal of this report is to provide a comprehensive status report of the research & development conducted in the context of the DARMA project by the end of the first quarter of fiscal year 2020. It follows in particular [LBS+19] and [PL19].
We summarize the narrow slot algorithms, including the thick electrically small depth case, conductive gaskets, the deep general depth case, multiple fasteners along the length, and finally varying slot width.
Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.
Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process ows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newtons method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h-adaptivity and dynamic load balancing are some of Arias more advanced capabilities.
The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.
SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.
SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported 'classical' phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.
Sensors continue to decrease in size and power. This report presents results of a market survey conducted in February 2020 for commercial off-the-self sensors with optimal size, weight, and power to be carried onboard a small unmanned aircraft system. For this report, Sandia National Laboratories considered sensors that can detect an object in three dimensions. The sensors that were researched are broken into three categories: radio detection and ranging sensors, stereo camera sensors, and light detection and ranging sensors.
What it is: A roughly spherical balloon constructed from light duty painter's plastic (0.31 mil high density polyethylene) and darkened with air float charcoal powder. Balloons typically range from 12-40 ft across depending on mission needs. How it works: Sunlight shines on the balloon, heating the air inside. The density difference due to the hot air in the balloon is sufficient to lift it up to 80,000 ft in the air
Sandia National Laboratories (SNL) has developed and produced a simple add-on kit (Pathogen Management Kit, or PMK) that can quickly add on to ventilators of all types to disinfect exhaled air to keep healthcare workers safe.
Sandia National Laboratories has hired Itasca Consulting Group, Inc., the authors of the FLAC3D geomechanics software, to couple FLAC3D with TOUGH3, the porous media flow solver. The work is being done to enable a coupled mechanical-thermal-hydraulic analysis of a potential criticality event in a dual purpose cannister (DPC). The U.S. Department of Energy Office of Spent Fuel and Waste Science & Technology is investigating the performance of DPCs for direct geological disposal of spent nuclear fuel. Post closure criticality control is an important aspect of this investigation. Over geological timescales, it is envisioned that the canister and canister overpack will develop fractures due to stress corrosion processes. A breach in the canister could allow groundwater to fill the canister. Fresh water is a neutron moderator; thus, if the canister internals and fuel assemblies have been sufficiently degraded, a criticality event could occur. Such an event would release enough energy to boil the water between the fuel rods and pressurize the cannister. This internal pressurization may cause the initial fractures in the canister and overpack to grow. It is important to understand the change in hydraulic transmissivity between the canister and surroundings for two reasons: first, because it may control the potential for and frequency of subsequent criticality events; second, because it will control the release of radionuclides from the canister. The motivation for this work is to better understand the potential for periodic criticality events, cannister damage, and release of radionuclides during a criticality event in a DPC.
This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.
Sandia's GEMINI-Scout Mine Rescue Robot is an unmanned ground vehicle designed to enter potentially hazardous environments to explore, assess, and evaluate dangerous situations first responders may face when conducting a rescue mission. GEMINI is approximately four feet long and two feet tall, which enables the robot to maneuver through small locations on rough terrains caused by earthquakes, fires, or radiological incidents. GEMINI uses track propulsion to climb stairs, travel through gravel and sand pits, pivot in place, and traverse 45-degree climbs with few problems. Furthermore, the vehicle's dual tracked-chassis design allows it to operate in hostile, dark, muddy, high-temperature, and explosive debris-strewn environments, while maintaining efficient ground mobility. The mobility and modularity of the vehicle allow for easy integration of sensors to conduct gas and temperature sensing and offers pan/tilt, zoom color, and thermal camera video streaming capabilities. The vehicle is also able to carry a payload of about 50 pounds of batteries and can handle an additional 200 pounds of payload, whether for additional diagnostics, supplies, or clothing for those trapped in an effected area. GEMINI is remotely operated through a wireless connection and an onboard computer running a customized embedded control application, which directly communicates to all onboard components except for the audio and video systems. When line of sight is not possible, operators use a shockresistant fiber optic cable to ensure continuous functionality of the vehicle. This allows for direct local control of the vehicle, which streams collected data back to the operator for enhanced situational awareness. In addition, the vehicle incorporates safety features such as explosion proof housing to ensure safe electronic operations in hazardous gas or flooded environments; a four-channel video link and two-way audio to ensure located survivors can communicate with operators; and an MSHA-approved multi-gas sensor to monitor air quality.
The Center for Disease Control has recommended that to reduce potential exposure to COVID-19 the public should wear cloth face coverings in public settings where other social distancing measures are difficult to maintain. These face coverings and other Emulated-Personal Protective Equipment (E-PPE) can be made by using Commonly Available Materials (CAMs). As E-PPE recommendations continue to flood the media, a Sandia COVID-19 LDRD effort, the Sandia E-PiPEline Team, systematically evaluated E-PPE design options considering their effectiveness, durability, build difficulty, build cost, and comfort. Using qualitative and semi-quantitative evaluation tools, results of the investigation are presented here to provide guidelines for home and office construction of E-PPE.
The Center for Disease Control has recommended that to reduce potential exposure to COVID-19 the public should wear cloth face coverings in public settings where other social distancing measures are difficult to maintain. These face coverings and other Emulated-Personal Protective Equipment (E-PPE) can be made by using Commonly Available Materials (CAMs). As E-PPE recommendations continue to flood the media, a Sandia COVID-19 LDRD effort, the Sandia E-PiPEline Team, systematically evaluated E-PPE design options considering their effectiveness, durability, build difficulty, build cost, and comfort. Using qualitative and semi-quantitative evaluation tools, results of the investigation are presented here to provide guidelines for home and office construction of E-PPE.
The Center for Disease Control has recommended the public to wear cloth face coverings in public settings that reduce potential exposure to COVID-19 where other social distancing measures are difficult to maintain (e.g., grocery stores and pharmacies) especially in areas of significant community based transmission. These face coverings and other Emulated-Personal Protective Equipment (EPPE) can be made by using Commonly Available Materials (CAMs). As part of the Sandia COVID-19 LDRD effort (funded under the Materials Science Investment Area), the Sandia E-PiPEline task evaluated E-PPE design options for face coverings and face shields considering their effectiveness, durability, build difficulty, build cost, and comfort. Observations from this investigation are presented here to provide guidelines for home construction of E-PPE. This executive summary includes a brief roadmap of the analysis methodology, two one-page handouts geared to be distributed to the public at large (one for E-PPE face coverings and one for E-PPE face shields), and additional observations regarding the potential solutions for E-PPE face coverings and face shields included to further support the one-page handouts.
Traditional dual-channel phase-monopulse and amplitude-monopulse antenna systems might electrically steer their difference-channel nulls by suitably adjusting characteristics of their constituent beams or lobes. A phase-monopulse systems' null might be steered by applying suitable relative phase shifts. An amplitude-monopulse systems' null might be steered by applying a suitable relative beam amplitude scaling. The steering of the null might be employed by a continuously mechanically-scanning antenna to stabilize the null direction over a series of radar pulses.
This report describes the results of a seven day effort to assist subject matter experts address a problem related to COVID-19. In the course of this effort, we analyzed the 29K documents provided as part of the White House's call to action. This involved applying a variety of natural language processing techniques and compression-based analytics in combination with visualization techniques and assessment with subject matter experts to pursue answers to a specific question. In this paper, we will describe the algorithms, the software, the study performed, and availability of the software developed during the effort.
Purpose: To remind personnel of potential power strip electrical hazards at the office and home Before you reach for that power strip, read this! This Lessons Learned Snapshot highlights the very real potential for electrical contact with power strip usage. It happens more often than you might think! But there are ways to reduce this potential.
For my internship project, I chose to evaluate the training program for the Sandia National Laboratories (SNL) Recorded Information Management (RIM) Program. The purpose of the evaluation was to identify: Where the current training meets federal and corporate requirements and if not, what gaps exist; If the training is accessible in a consistent location; If the training contains a consistent message about the responsibility of managing documents; If the training is organized by records management subject areas.
The research presented in this paper is to be a contribution to a larger research project. The project was designed to test the safety standards of Sandia National Laboratories' ACRR. The MCNP model of the ACRR core is being used to simulate the environment. To test the safety standards, the dimensions and density of the fuel in the core were varied randomly. Each rod had its own dimension and density assigned to it randomly. Ten different cases of the ACRR core were randomly created. The different cases were studied using 640-group neutron fluxes. The fluxes were compared to see how the variations affected the system.
In response to anticipated resource shortfalls related to the treatment and testing of COVID-19, many communities are planning to build additional facilities to increase capacity. These facilities include field hospitals, testing centers, mobile manufacturing units, and distribution centers. In many cases, these facilities are intended to be temporary and are designed to meet an immediate need. When deciding where to place new facilities many factors need to be considered, including the feasibility of potential locations, existing resource availability, anticipated demand, and accessibility between patients and the new facility. In this project, a facility location optimization model was developed to integrate these key pieces of information to help decision makers identify the best place, or places, to build a facility to meet anticipated resource demands. The facility location optimization model uses the location of existing resources and the anticipated resource demand at each location to minimize the distance a patient must travel to get to the resource they need. The optimization formulation is presented below. The model was designed to operate at the county scale, where patients are grouped per county. This assumption can be modified to integrate other scales or include individual patients.
Genome editing technologies, particularly those based on zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR (clustered regularly interspaced short palindromic repeat DNA sequences)/Cas9 are rapidly progressing into clinical trials. Most clinical use of CRISPR to date has focused on ex vivo gene editing of cells followed by their re-introduction back into the patient. The ex vivo editing approach is highly effective for many disease states, including cancers and sickle cell disease, but ideally genome editing would also be applied to diseases which require cell modification in vivo. However, in vivo use of CRISPR technologies can be confounded by problems such as off-target editing, inefficient or off-target delivery, and stimulation of counterproductive immune responses. Current research addressing these issues may provide new opportunities for use of CRISPR in the clinical space. In this review, we examine the current status and scientific basis of clinical trials featuring ZFNs, TALENs, and CRISPR-based genome editing, the known limitations of CRISPR use in humans, and the rapidly developing CRISPR engineering space that should lay the groundwork for further translation to clinical application.
Hansen, Nils H.; Liao, Handong; Kang, Shiqing; Zhang, Feng; Yang, B.
The influence of ozone addition on the low-temperature oxidation of dimethyl ether (DME) was investigated experimentally in an atmospheric-pressure jet-stirred reactor, over the temperature range of 400–800 K. Detailed speciation information was obtained by employing synchrotron vacuum ultraviolet photoionization mass spectrometry. Experimental results revealed that the ozone addition had a positive influence on the production of the highly reactive intermediates. Moreover, the low-temperature reactivity of DME was significantly enhanced, which resulted in the broadening of the temperature window of fuel consumption and intermediates formation at lower temperatures. Therefore, novel experimental data of the low temperature regime (400–500 K) could be obtained. The data set of this special temperature regime yielded insights into the DME low-temperature kinetics, which were further supported with modeling analysis based on two existing DME models (Metcalfe et al., 2013; Wang et al., 2015) combined with an ozone sub-mechanism (Zhao et al., 2016). The analysis showed that temperature-sensitive reactions such as the second oxygen channel could be nearly “frozen” at this low temperature (T < 440 K). Furthermore, the production of some intermediates was found to be strongly governed by reaction pairs, such as CH3OCH2 + O2 = CH3OCH2O2 and CH3OCH2 + O2 = 2CH2O + OH for the CH2O formation. This finding could be useful for examining branching ratios in both models, and the analysis suggested the further modification of the branching ratios for the oxygen addition to CH3OCH2O2 pathways and the CH3OCH2O2 self-reactions were required. Finally, the influences of the O3 addition in the sensitive reactions of the fuel initial low-temperature oxidation were investigated in this work. It was interesting to note that O3 addition could change the dominating reactions in the initial low-temperature oxidation, by the addition of some O3-related pathways with relatively high sensitivity.
Holland, Rayne; Khan, M.A.H.; Chhantyal-Pun, Rabi; Orr-Ewing, Andrew J.; Percival, Carl J.; Taatjes, Craig A.; Shallcross, Dudley E.
Perfluorooctanoic acid, PFOA, is one of the many concerning pollutants in our atmosphere; it is highly resistant to environmental degradation processes, which enables it to accumulate biologically. With direct routes of this chemical to the environment decreasing, as a consequence of the industrial phase out of PFOA, it has become more important to accurately model the effects of indirect production routes, such as environmental degradation of precursors; e.g., fluorotelomer alcohols (FTOHs). The study reported here investigates the chemistry, physical loss and transport of PFOA and its precursors, FTOHs, throughout the troposphere using a 3D global chemical transport model, STOCHEM-CRI. Moreover, this investigation includes an important loss process of PFOA in the atmosphere via the addition of the stabilised Criegee intermediates, hereby referred to as the "Criegee Field. " Whilst reaction with Criegee intermediates is a significant atmospheric loss process of PFOA, it does not result in its permanent removal from the atmosphere. The atmospheric fate of the resultant hydroperoxide product from the reaction of PFOA and Criegee intermediates resulted in a ≈0.04 Gg year-1 increase in the production flux of PFOA. Furthermore, the physical loss of the hydroperoxide product from the atmosphere (i.e., deposition), whilst decreasing the atmospheric concentration, is also likely to result in the reformation of PFOA in environmental aqueous phases, such as clouds, precipitation, oceans and lakes. As such, removal facilitated by the "Criegee Field" is likely to simply result in the acceleration of PFOA transfer to the surface (with an expected decrease in PFOA atmospheric lifetime of ≈10 h, on average from ca. ≈80 h without Criegee loss to 70 h with Criegee loss).
High-temperature optical analysis of three different InGaN/GaN multiple quantum well (MQW) light-emitting diode (LED) structures (peak wavelength λp = 448, 467, and 515 nm) is conducted for possible integration as an optocoupler emitter in high-density power electronic modules. The commercially available LEDs, primarily used in the display (λp = 467 and 515 nm) and lighting (λp = 448 nm) applications, are studied and compared to evaluate if they can satisfy the light output requirements in the optocouplers at high temperatures. The temperature- and intensity-dependent electroluminescence (T-IDEL) measurement technique is used to study the internal quantum efficiency (IQE) of the LEDs. All three LEDs exhibit above 70% IQE at 500 K and stable operation at 800 K without flickering or failure. At 800 K, a promising IQE of above 40% is observed for blue for display (BD) (λp = 467 nm) and green for display (GD) (λp = 515 nm) samples. The blue for light (BL) (λp = 448 nm) sample shows 24% IQE at 800 K.
Complex metal hydrides provide high-density hydrogen storage, which is essential for vehicular applications. However, the practical application of these materials is limited by thermodynamic and kinetic barriers present during the dehydrogenation and rehydrogenation processes as new phases form inside parent phases. An improved understanding of the mixed-phase mesostructures and their interfaces will assist in improving cyclability. In this work, the phase evolution during hydrogenation of lithium nitride and dehydrogenation of lithium amide with lithium hydride is probed with scanning transmission X-ray microscopy at the nitrogen K edge. With this technique, core–shell structures are observed in particles of both partially hydrogenated Li3N and partially dehydrogenated LiNH2 + 2LiH. To generate these structures, the rate-limiting step must shift from internal hydrogen diffusion during hydrogenation to the formation of hydrogen gas at the surface during desorption.
The Microgrid Design Toolkit (MDT) supports decision analysis for new ("greenfield") microgrid designs as well as microgrids with existing infrastructure. The current version of MDT includes two main capabilities. The first capability, the Microgrid Sizing Capability (MSC), is used to determine the size and composition of a new, grid connected microgrid in the early stages of the design process. MSC is focused on developing a microgrid that is economically viable when connected to the grid. The second capability is focused on designing a microgrid for operation in islanded mode. This second capability relies on two models: the Technology Management Optimization (TMO) model and Performance Reliability Model (PRM).
Sandia National Laboratories is part of the government test and evaluation team for the Defense Advanced Research Projects Agency Collection and Monitoring via Planning for Active Situational Scenarios program. The program is designed to better understand competition in the area between peace and conventional conflict when adversary actions are subtle and difficult to detect. For the purposes of test and evaluation, Sandia conducted a range of activities for the program: creation of the Grey Zone Test Range; design of the data stream for a user experiment conducted with U.S. Indo-Pacific Command; design, implementation, and execution of the formal evaluation; and analysis and summary of the evaluation results. This report details Sandia's activities and provides additional information on the Grey Zone Test Range urban simulation environment developed to evaluate the performer technologies.
Telecommuting at Sandia and within the Federal workforce has been optional for years. With the COVID-19 crisis, social distancing makes this a necessity that may be with us for some weeks and even months to come. All of us - whether new to telecommuting or an experienced remote worker - are challenged with how to efficiently and effectively adapt to this new work environment, while maintaining work-life balance to adequately care for ourselves and loved ones.
This Sandia National Laboratories, New Mexico Environmental Restoration Operations (ER) Consolidated Quarterly Report (ER Quarterly Report) fulfills all quarterly reporting requirements set forth in the Compliance Order on Consent. Table I-1 lists the six sites remaining in the corrective action process.
Ronevich, Joseph A.; Song, Eun J.; Feng, Zhili; Wang, Yanli; D'Elia, Christopher; Hill, Michael R.
Fatigue crack growth rates (FCGR) of multiple X100 pipeline steel welds and heat affected zones were measured in high-pressure hydrogen gas to investigate their behavior compared to lower strength pipeline welds. A total of five high strength welds and two heat affected zones (HAZ) were examined all of which were fabricated using the same X100 base material. Different welding wires and techniques were used to fabricate the welds to provide a variety of end products to evaluate susceptibility to fatigue in high pressure hydrogen gas. Residual stresses were measured for each weld and HAZ using the slitting method and the effect of residual stress on the stress intensity factor, Kres, was determined. Using Kres, the fatigue crack growth rate curves were corrected to remove the effects of residual stress by examining the influence of Kres on stress ratio, R. Comparisons were then made between the high strength welds, which were corrected for residual stress, and lower strength welds from the literature. It was found that the higher strength welds and heat affected zones exhibited comparable fatigue crack growth rates to lower strength welds, as the FCGR data of the high strength welds overlaid the lower strength welds. This suggests that despite distinct differences in strength and microstructure between the different welds, hydrogen-assisted fatigue crack growth susceptibility is similar. A comparison was made between the Kres measured in extracted coupons and residual stress estimates provided in relevant welded pipe assessment standards such as API 579-1/ASME FFS-1. It was found the residual stress values in the test coupons extracted from welded pipe were significantly lower than those expected in the intact welded pipes and highlights the importance in quantifying and removing coupon residual stresses when fatigue crack growth rates are measured and including expected weld joint residual stress when making structural assessments.
Li, Tao; Zhou, Bo; Frank, Jonathan H.; Dreizler, Andreas; Bohm, Benjamin
Abstract: The development of high-speed volumetric laser-induced fluorescence measurements of formaldehyde (CH 2O -LIF) using a pulse-burst laser operated at a repetition rate of 100kHz is presented. A novel laser scanning system employing an acousto-optic deflector (AOD) enables quasi-4D CH 2O -LIF imaging at a scan frequency of 10kHz. The diagnostic capability of time-resolved volumetric imaging is demonstrated in a partially premixed DME/air lifted turbulent jet flame near the flame base. Simultaneous imaging of laser beam profiles is performed to account for the laser pulse energy fluctuation and laser sheet inhomogeneity. With the accurate registration of laser sheet positions, the volumetric reconstruction of CH 2O -LIF signals is performed within a detection volume of 17.3×11.9×2.3mm3 with an average out-of-plane spatial resolution of 250μm. A surface detection algorithm with adaptive thresholding is used to determine the global maximum intensity gradient by calculating gradient percentiles. The flame topology characteristics are investigated by evaluating the 3D curvatures of CH 2O surfaces. Curvatures calculated using 2D data systematically underestimate the full 3D curvature due to the lack of out-of-plane information. The inner surfaces near the turbulent fuel jet exhibit higher probabilities of large mean curvature than the outer surfaces. The saddle and cylindrical structures are dominant on both the inner and outer surfaces and the elliptic structures occur with lower probability. The results suggest that the damping of turbulent fluctuations by the temperature increase through the CH 2O region reduces the curvature, but the local structure topology remains self-similar. Graphic abstract: [Figure not available: see fulltext.].
Data assimilation techniques are investigated for integrating high-speed high-resolution experimental data into large-eddy simulations. To this end, an ensemble Kalman filter is employed to assimilate velocity measurements of a turbulent jet at a Reynolds number of 13,500 into simulations. The goal of the current work is to examine the behavior of the assimilation algorithm for state estimation of turbulent flows that are of relevance to engineering applications. This is accomplished by investigating the impact that localization, measurement uncertainties, assimilation frequency, data sparsity and ensemble size have on the estimated state vector. For the flow configuration and computational setup considered in this study an optimal value of the localization radius is identified, which minimizes the error between experimental data and state vector. The impact of experimental uncertainties on the state estimation is demonstrated to provide solution bounds on the assimilation algorithm. It is found that increasing the number of ensembles has a positive impact on the state estimation. In comparison, decreasing the assimilation frequency or reducing the experimental data available for assimilation is found to have a negative impact on the state estimation. These findings demonstrate the viability of assimilating measurements into numerical simulations to improve state estimates, to support parameter evaluations and to guide model assessments.
Ultrafast x-ray imagers developed at Sandia National Laboratories are a transformative diagnostic tool in inertial confinement fusion and high energy density physics experiments. The nanosecond time scales on which these devices operate are a regime with little precedent, and applicable characterization procedures are still developing. This paper presents pulsed x-ray characterization of the Icarus imager under a variety of illumination levels and timing modes. Results are presented for linearity of response, absolute sensitivity, variation of response with gate width, and image quality.
Sauer, Kirsten; Rock, Marlena; Caporuscio, Florie; Hardin, Ernest H.
The stability of neutron absorber composite materials at hydrothermal conditions was tested in a series of two-week experiments to mimic spent nuclear fuel disposal. Coupons, composed of boron carbide (B4C) sintered with and encased in aluminum, were increasingly altered in experiments at 150, 230, and 300 °C and pressures of 150 bar. Alteration of aluminum to boehmite (γ-AlO(OH)) and hydrogen gas generation occurred over the range of investigated temperatures, but is most significant at 300 °C. The formation of boron-bearing mineral phases was not detected; however, aqueous boron was present in the reaction fluids.
By engineering atomic geometries composed of nearly 1000 atomic segments embedded in micro-resonators, we observe Bragg resonances induced by the atomic lattice at the telecommunication wavelength. The geometrical arrangement of erbium atoms into a lattice inside a silicon nitride (SiN) microring resonator reduces the scattering loss at a wavelength commensurate with the lattice. We confirm dependency of light emission to the atomic positions and lattice spacing and also observe Fano interference between resonant modes in the system.
Quiroz-Arita, Carlos; Blaylock, Myra L.; Gharagozloo, Patricia E.; Bark, David; Prasad Dasi, Lakshmi; Bradley, Thomas H.
Turbulent mixing in pilot-scale cultivation systems influences the productivity of photoautotrophic cultures. We studied turbulent mixing by applying particle image velocimetry and acoustic doppler velocimetry to pilot-scale, flat-panel photobioreactor, and open-channel raceway. Mixing energy inputs were varied from 0.1 to 2.1 W·m−3. The experimental results were used to quantify turbulence and to validate computational fluid dynamics models, from which Lagrangian representations of the fluid motion in these reactors were derived. The results of this investigation demonstrated that differences in mixing energy input do not significantly impact the structure of turbulence and the light/dark cycling frequencies experienced by photoautotrophic cells within the reactors. The experimental and computational results of our research demonstrated that well-mixed conditions exist in pilot-scale, flat-panel photobioreactors and open-channel raceways, even for relatively low mixing energy inputs.
Two different cell constructions are currently being investigated for the LFP/Graphite system. These cells have undergone continuous cycling until the target of 1000 cycles was reached.
Goldberg, Benjamin M.; Chng, Tat L.; Naphade, Maya; Adamovich, Igor V.; Starikovskaia, Svetlana M.
Electric-field-induced second-harmonic generation, or E-FISH, has received renewed interest as a nonintrusive tool for probing electric fields in gas discharges and plasmas using ultrashort laser pulses. An important contribution of this work lies in establishing that the E-FISH method works effectively in the nanosecond regime, yielding field sensitivities of about a kV/cm at atmospheric pressure from a 16 ns pulse. This is expected to broaden its applicability within the plasma community, given the wider access to conventional nanosecond laser sources. A Pockels-cell-based pulse-slicing scheme, which may be readily integrated with such nanosecond laser systems, is shown to be a complementary and cost-effective option for improving the time resolution of the electric field measurement. Using this scheme, a time resolution of ∼3 ns is achieved, without any detriment to the signal sensitivity. This could prove invaluable for nonequilibrium plasma applications, where time resolution of a few nanoseconds or less is often critical. Finally, we take advantage of the field vector sensitivity of the E-FISH signal to demonstrate simultaneous measurements of both the horizontal and vertical components of the electric field.
This report details the test setup, process, and results for radiated susceptibility testing of multicrystalline silicon photovoltaic (PV) modules as part of the EMP-Resilient Electric Grid Grand Challenge Laboratory Directed Research and Development (LDRD) project at Sandia National Laboratories. Testing was conducted over October 10-17, 2019, where 8 photovoltaic modules were exposed to E1 transient pulses with peak field levels up to 100 kV/m. Modules were terminated in a resistive load representing connected components. State of health testing conducted via I-V curve tracing of the photovoltaic modules showed no observable loss of device function due to large electric field transients. Differential mode currents were measured on the order of 10's of amps for up to a microsecond following the radiated field pulse. Common mode currents took the form of a damped sinusoid with a maximum peak of 10's to 100's of amps with a resonance near 60 MHz.
Climate change and its impacts on average temperature, water supply, agriculture, coastal flooding, biodiversity, social and economic stability, human migration, and overall global stability is one of the leading threats emerging for humanity, and is likely to increase as a threat for decades to come. This report represents a snapshot of climate change data, information, and understanding, as of 2018, and can serve as a kind of benchmark for changes going forward in time. This report covers temperature change and heat effects, atmospheric and ocean circulation, freshwater supply changes, sea level rise and flooding, extreme climate events, oceanic deoxygenation, land degredation, social and economic changes, migration, and food. Data and information on all these issues represent a compelling body of knowledge supporting the scientific hypothesis describing climate change, and an important mile marker in our effort to track the unfolding nature of the problem. A slide presentation that summarizes the content of this report is in Appendix A and can serve as an executive summary for the report.
The IEC 61853 PV module energy rating standard requires measuring module power (and hence efficiency) over a matrix of irradiance and temperature conditions. These matrix points represent nearly the full range of operating conditions encountered in the field in all but the most extreme locations, and create an opportunity to develop alternative approaches to existing modelssuch as the single-diode models and the Sandia Array Performance Modelfor calculating system performance. This report begins by discussing the bilinear interpolation and extrapolation method from IEC 61853-3, and then describes four existing model-based methods that could be used with matrix measurements. Then a new model is developed and all options are compared according to seven objective criteria using the matrix measurements of four PV modules of different technologies. The results show that the new model is an excellent candidate for launching power matrix-based PV system simulations.
Space rendezvous and proximity operations are increasing in numbers, enabling inspections, diagnostics, and maintenance of on-orbit systems. Because collision, loss of control, and unintended damage can impact the system under examination -- and at the extreme, cause system break-up and space debris -- the safety practices for rendezvous and proximity operations can have significant implications for national security. This study examines the applicability of the Always/Never surety framework, which was developed for United States nuclear weapons, as a model safety basis for unmanned space proximity operations. This unclassified framework has understandable safety approaches and principles and focuses on a system being always safenever unsafe. The authors consider that the adapting the framework might present a means for standardization across government and commerce, encouraging a consistent approach and a set of clarifying safety principles and applications for rendezvous and proximity operations. The framework also offers a consistent taxonomy, presents safety and reliability requirements organized by four environment categories, defines accident or abnormal conditions, contributes a strategy for identifying hostile and tactical environments, and enables decision-making for determining if conditions are safe for proximity space operations.
This report examines the sensitivity of salt cavern Mechanical Integrity Tests (MIT) to uncertainties in key test parameters. MIT's are used by cavern operators to detect and quantify leak rates in access wells in underground salt storage caverns and involve a suite of measured and assumed parameters that have a direct impact on the sensitivity of the testing to detect actual leaks from the cavern storage system. Determining the sensitivity of the testing to these different parameters provides a basis for understanding the results from, and informing the design criteria for, this type of testing. Without fully understanding the sensitivity of the test to the testing parameters, it is possible that the test results may not accurately reflect the integrity of the cavern system; an actual leak may be missed, or an intact system may be interpreted as leaking. This report reviews the main parameters included in MITs and examines how selected changes in their values can impact test results. The deviations used in the sensitivity analyses were designed to be within the ranges believed to be similar to those which may be encountered during testing. The results show that small, plausible fluctuations in some of the parameters measured values can have a significant impact on the testing results. Of the parameters studied here, the sensitivity analyses showed the order of importance to be (from highest to lowest): nitrogen-oil interface depth measurement, well bore temperature, well head pressure, and finally the internal geometry of the testing interval.
Milestone Description: Enhance Nalu-Wind's actuator disc model through hardening, documenting, stress-testing, verifying, and validating. Existing workflows will be improved by reducing the data output stream, and by making the analysis capabilities more modular and generally better. These model capabilities are needed by other A2e areas, namely Wake Dynamics, AWAKEN, and VV&UQ.
The international export of handheld spectroscopy detectors by the National Nuclear Security Administration (NNSA) to partner states will provide state regulatory authorities and nuclear material owners a way to improve accountancy for accidental gains of nuclear material, including the provision of reports to the IAEA, in order to meet their safeguards agreements. International Atomic Energy Agency (IAEA) safeguards agreements for non-nuclear weapons states requires accountancy for all nuclear material. As defined in Article XX of the IAEA statute, nuclear material includes source materials: "uranium containing the mixture of isotopes occurring in nature," and special fissionable material: "plutonium-239; uranium-233; uranium enriched in the isotopes 235 or 233". For IAEA Member States to meet their requirements under comprehensive safeguards agreements (CSA), safeguards are to be applied on "all source or special fissionable material," which "includes all nuclear material subject to IAEA safeguards". Therefore, accidental gains and losses of nuclear material must be reported to the IAEA. An accidental gain occurs when a state unexpectedly adds nuclear material to their inventory by various means such as seizing smuggled material or the discovery of legacy items previously unaccounted for. The material type and quantity must be added to the State's inventory by updating domestic records and then communicated to the IAEA.
The availability of repair garage infrastructure for hydrogen fuel cell vehicles is becoming increasingly important for future industry growth. Ventilation requirements for hydrogen fuel cell vehicles can affect both retrofitted and purpose-built repair garages and the costs associated with these requirements can be significant. A hazard and operability study (HAZOP) was performed to identify key risk-significant scenarios related to hydrogen vehicles in a repair garage. Detailed simulations and modeling were performed using appropriate computational tools to estimate the location, behavior, and severity of hydrogen release based on key HAZOP scenarios. This work compares current fire code requirements to an alternate ventilation strategy to further reduce potentially hazardous conditions. Overall, the amount of flammable mass of hydrogen at any one time in the simulation is low compared to the total mass of hydrogen released, due to the low flow rate of a low pressure release. It is shown that position, direction, and velocity of ventilation have a significant impact on the amount of instantaneous flammable mass in the domain.
Continuing previous efforts to investigate and develop the Unclassified Radioisotope Algorithm, the goal of the FY19-FY20 effort was to develop a prototype detector system which uses the algorithm to confirm warhead attributes related to the presence of either weapons grade plutonium (WGPu) or highly enriched uranium (HEU). The final deliverable is a prototype attribute measurement system built with common, commercially available gamma radiation detector components, capable of confirming the presence of specific, complex radioactive sources of interest, without the collection and storage of gamma energy spectra. This is accomplished by processing each gamma pulse as it is received, applying weight values based on the energy and incrementing or decrementing scalar counters which can be compared with expected values to determine if the measured source is consistent with WGPu or HEU. This report documents the design of the prototype system as well as the development of the algorithm and performance testing results. While the previously conceptualized, simple algorithm resulted in a prohibitive amount of false positives, the goal for a simple attribute measurement system capable of verifying Ba-133 and Ra-226 (weapons grade plutonium and highly enriched uranium surrogate testing sources) at over 95% accuracy with sub 5% false positive rate was demonstrated.
There are numerous vehicles which utilize alternative fuels, or fuels that differ from typical hydrocarbons such as gasoline and diesel, throughout the world. Alternative vehicles include those running on the combustion of natural gas and propane as well as electrical drive vehicles utilizing batteries or hydrogen as energy storage. Because the number of alternative fuels vehicles is expected to increase significantly, it is important to analyze the hazards and risks involved with these new technologies with respect to the regulations related to specific transport infrastructure, such as bridges and tunnels. This report focuses on hazards presented by hydrogen fuel cell electric vehicles that are different from traditional fuels. There are numerous scientific research and analysis publications on hydrogen hazards in tunnel scenarios; however, compiling the data to make conclusions can be a difficult process for tunnel owners and authorities having jurisdiction over tunnels. This report provides a summary of the available literature characterizing hazards presented by hydrogen fuel cell electric vehicles, including light-duty, medium and heavy-duty, as well as buses. Research characterizing both worst-case and credible scenarios, as well as risk-based analysis, is summarized. Gaps in the research are identified to guide future research efforts to provide a complete analysis of the hazards and recommendations for the safe use of hydrogen fuel cell electric vehicles in tunnels.
The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr(WO4)2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr(WO4)2. In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y2O3 (charge compensation).
The U.S. Nuclear Regulatory Commission (NRC) has interacted with vendors pursuing the commercialization of micro-reactors (i.e., reactors capable of producing about 1 MW(th) to 20 MW(th) of energy from nuclear fission). It is envisioned that micro-reactors could be assembled and fueled in a factory and shipped to a site. Many of the sites are expected to be remote locations requiring off-grid power or in some cases military bases. The objective of this effort is to explore the technical issues and the approach required to reach a finding of "reasonable assurance of public health and safety" for this new and different class of reactors. The analysis performed here leverages available micro-reactor design and testing data available from national laboratory experience as well as commercial design information to explore technical issues. Some factors considered include source term, accidents that would need to be analyzed, and the extent of the probabilistic risk assessment (PRA). The technical evaluation was performed within the framework of the Licensing Modernization Project (LMP) to identify licensing basis events, classification of structures, systems and components, and defense-in-depth needed to provide regulatory certainty. With this framework and technical evaluation in mind, the scope and content of a micro-reactor licensing application is discussed.
Under its Grid Modernization Initiative, the U.S. Department of Energy (DOE), in collaboration with energy industry stakeholders developed a multi-year research plan to support modernizing the electric grid. One of the foundational projects for accelerating modernization efforts is information and communications technology interoperability. A key element of this project has been the development of a methodology for engaging ecosystems related to grid integration to create roadmaps that advance the ease of integration of related smart technology. This document is the product of activities undertaken in 2017 through 2019. It provides a Cybersecurity Plan describing the technology to be adopted in the project with details as per the GMLC Call document.
A new commercial design of BAM friction tester that utilizes a programmable servo motor was tested in comparison to the traditional cam-driven model. Displacement and velocity profiles were analyzed for both designs; significant differences were found between the two designs, most notable of which is that the traditional cam-driven unit has plate velocities that can be 50–75 % greater than those seen on the servo-driven model. Five energetic materials were also tested on each machine, including PETN, RDX, HMX, CL-20, and HNAB. Results from the servo-driven model generally showed slightly less sensitivity when compared to the cam driven model, though the magnitude of the difference is not significant enough to require modification of safe handling procedures for the materials tested.
IEEE Journal of Emerging and Selected Topics in Power Electronics
Augustine, Sijo; Reno, Matthew J.; Brahma, Sukumar M.; Lavrova, Olga
This report presents a novel fault detection, characterization, and fault current control algorithm for a standalone solar-photovoltaic (PV) based dc microgrids. The protection scheme is based on the current derivative algorithm. The overcurrent and current directional/differential comparison based protection schemes are incorporated for the dc microgrid fault characterization. For a low impedance fault, the fault current is controlled based on the current/voltage thresholds and current direction. Generally, the droop method is used to control the power-sharing between the converters by controlling the reference voltage. In this article, an adaptive droop scheme is also proposed to control the fault current by calculating a virtual resistance R droop , and to control the converter output reference voltage. For a high impedance fault, differential comparison method is used to characterize the fault. These algorithms effectively control the converter pulsewidth and reduce the flow of source current from a particular converter, which helps to increase the fault clearing time. Additionally, a trip signal is sent to the corresponding dc circuit breaker (DCCB), to isolate the faulted converter, feeder or a dc bus. The dc microgrid protection design procedure is detailed, and the performance of the proposed method is verified by simulation analysis.
In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.
In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution modeling and simulation (M&S) tools for modeling of light water reactors (LWRs). This hub, the Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called the Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including CTF (the CASL version of Coolant Boiling in Rod Arrays— Three Field (COBRA-TF)), Michigan Parallel Characteristics Transport (MPACT), BISON, and Materials Performance and Optimization (MPO) Advanced Model for Boron Analysis (MAMBA). As part of its M&S efforts, CASL has identified various challenge problems, including Crud Induced Power Shift (CIPS), Crud-Induced Localized Corrosion (CILC), Pellet-Cladding Interaction (PCI), and Departure from Nucleate Boiling (DNB). This work addresses CASL milestone L2:VVI.P19.03, which focuses on uncertainty quantification of crud, which is relevant to both CIPS and CILC. This is achieved through an analysis and separate effects validation of the thermal hydraulic phenomenon known as subcooled boiling. As part of this work, various sources of experimental data are examined and compared to different options for empirical modeling of subcooled boiling. Through this analysis, a complete understanding of the underlying models and their implementation details are understood. A subset of these data are incorporated into a separate effects validation study of CTF. The Westinghouse Advanced Loop Tester (WALT) and Rohsenow experiments are modeled, and it is shown that the newly-implemented Gorenflo correlation is more accurate than the existing Chen and Thom correlations.
The selection of austenitic stainless steels for hydrogen service is challenging since there are few intrinsic metrics that relate alloy composition to hydrogen degradation. One such metric, explored here, is intrinsic stacking fault energy. Stacking fault energy has an influence on the character and structure of dislocations and on the formation of secondary crystalline phases created during mechanical deformation in austenitic alloys. In this work, a data-driven model for the intrinsic stacking fault energy of common austenitic stainless steel alloys is applied to compare the relative degradation of tensile performance in the presence of hydrogen. A transition in the tensile reduction of area of both 300-series and manganese stabilized stainless steels is observed at a calculated stacking fault energy of approximately 43 mJ m-2, below which pronounced hydrogen degradation on tensile ductility is observed. The model is also applied to suggest alloying strategies for low nickel austenitic stainless steels for hydrogen service. Lastly, through this investigation, we find that calculated intrinsic stacking fault energy is a high-throughput screening metric that enables the ranking of the performance of a diverse range of austenitic stainless steel compositions, as well as the identification of new alloys, with regard to hydrogen compatibility.
This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. The convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.
Terahertz semiconductor quantum-cascade lasers (QCLs) are widely implemented with metallic cavities that support low-loss plasmonic optical modes at long wavelengths. However, resonant optical modes in such cavities suffer from poor radiative characteristics due to their subwavelength transverse dimensions. Consequently, single-mode terahertz QCLs with metallic cavities and large (> 100 mW) output power have only been realized in the surface-emitting configuration that affords a large radiating surface. Here, we demonstrate a method to enhance radiative outcoupling from such plasmonic lasers for high-power emission in the edge-emitting (end-fire or longitudinal) direction. Single-sided plasmon waves propagating in vacuum are resonantly excited in surrounding medium of metallic cavities with the QCL semiconductor medium. The vacuum guided plasmon waves with a large wavefront phase-lock multiple metallic cavities longitudinally, which leads to intense radiation in multiple directions, including that in the longitudinal direction in a narrow single-lobed beam. The multicavity array radiates predominantly in a single spectral mode. A peak-power output of 260 mW and a slope efficiency of 303 mW/A are measured for the end-fire beam from a 3.3 THz QCL operating at 54 K in a Stirling cooler. Single-mode operation and lithographic tuning across a bandwidth of ∼ 150 GHz are demonstrated.
This is an addendum to the Sierra/SolidMechanics 4.56 User's Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State's International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 4.56 User's Guide should be referenced for most general descriptions of code capability and use.
This user's guide documents capabilities in Sierra/SolidMechanics which remain "in-development" and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.56 User's Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.
This document describes how to build a 5-gore, 5.8 m diameter heliotrope solar hot air balloon. This is a fairly straightforward process, but it is painstaking. When making the balloons, make sure not to wear anything that can snag the material (badges, etc). Sharp objects or corners should not be present. When laying out, folding, and cutting gores, it is best to wear socks instead of shoes. Tape should never be pulled off of a balloon. If it accidentally touches the balloon material, it should be left in place or cut free. Also, when adding tape (either intentionally or not), no sticky parts should be left. Sticky parts should either be cut free or taped over. Otherwise, the sticky part will grab the balloon envelope and tear it. You are building a 20 ft sphere out of material thinner than a grocery bag — the best guidance is just to use common sense.
The study of infrasound (acoustic) and gravity waves sources and propagation in the atmosphere of a planet gives us precious insight on atmosphere dynamics, climate, and even internal structure. The implementation of modern pressure sensors with high rate sampling on stratospheric balloons is improving their study. We analyzed the data from the National Aeronautics and Space Administration Ultra Long Duration Balloon mission (16 May to 30 June 2016). Here, we focus on the balloon's transit of the Andes Mountains. We detected gravity waves that are associated to troposphere convective activity and mountain waves. An increase of the horizontal wavelengths from 50 to 70 km with increasing distance to the mountains is favoring the presence of mountain waves. We also report on the detection of infrasounds generated by the mountains in the 0.01–0.1 Hz range with a pressure amplitude increase by a factor 2 relative background signal. Besides, we characterized the decrease of microbaroms power when the balloon was flying away from the ocean coast. These observations suggest, in a way similar to microseisms for seismometers, that microbaroms are the main background noise sources recorded in the stratosphere even far from the ocean sources. Finally, we observed a broadband signal above the Andes, between 0.45 and 2 Hz, probably associated with a thunderstorm. The diversity of geophysical phenomena captured in less than a day of observation stresses the interest of high rate pressure sensors on board long-duration balloon missions.
We present a theoretical model that predicts the peak strength of polycrystalline metals based on the activation energy (or stress) required to cause deformation via amorphization. Building on extensive earlier work, this model is based purely on materials properties, requires no adjustable parameters, and is shown to accurately predict the strength of four exemplar metals (fcc, bcc, and hcp, and an alloy). This framework reveals new routes for design of more complex high-strength materials systems, such as compositionally complex alloys, multiphase systems, nonmetals, and composite structures.
Terahertz (THz) technology has shown promise for several applications, but limitations in sources and detectors have prevented broader adoption. Existing THz detectors are rigid, planar, and fabricated using complex technology, making it difficult to integrate into systems. Here we demonstrate THz detectors fabricated by inkjet printing on submicrometer thick, ultraflexible substrates. By developing p- and n-type carbon nanotube inks, we achieve optically thick p–n junction and p-type devices, enabling antenna-free pixels for THz imaging. By further designing and folding the printed devices, we realize origami-inspired architectures with improved performance over single devices, achieving a noise-equivalent power of 12 nW/Hz1/2 at room temperature with no voltage bias. Our approach opens avenues for nonplanar, foldable, deployable, insertable, and retractable THz detectors for applications in nondestructive inspection.
Presented in this document are the theoretical aspects of capabilities contained in the Sierra/SM code. This manuscript serves as an ideal starting point for understanding the theoretical foundations of the code. For a comprehensive study of these capabilities, the reader is encouraged to explore the many references to scientific articles and textbooks contained in this manual. It is important to point out that some capabilities are still in development and may not be presented in this document. Further updates to this manuscript will be made as these capabilities come closer to production level.
Barnes, Philip M.; Wallace, Laura M.; Saffer, Demian M.; Bell, Rebecca E.; Underwood, Michael B.; Fagereng, Ake; Meneghini, Francesca; Savage, Heather M.; Rabinowitz, Hannah S.; Morgan, Julia K.; Kutterolf, Steffen; Hashimoto, Yoshitaka; Engelmann De Oliveira, Christie H.; Noda, Atsushi; Crundwell, Martin P.; Shepherd, Claire L.; Woodhouse, Adam D.; Harris, Robert N.; Wang, Maomao; Henrys, Stuart; Barker, Daniel H.N.; Petronotis, Katerina E.; Bourlange, Sylvain M.; Clennell, Michael B.; Cook, Ann E.; Dugan, Brandon E.; Elger, Judith; Fulton, Patrick M.; Gamboa, Davide; Greve, Annika; Han, Shuoshuo; Hupers, Andre; Ikari, Matt J.; Ito, Yoshihiro; Kim, Gil Y.; Koge, Hiroaki; Lee, Hikweon; Li, Xuesen; Luo, Min; Malie, Pierre R.; Moore, Gregory F.; Mountjoy, Joshu J.; Mcnamara, David D.; Paganoni, Matteo; Screaton, Elizabeth J.; Shankar, Uma; Shreedharan, Srisharan; Solomon, Evan A.; Wang, Xiujuan; Wu, Hung-Yu; Pecher, Ingo A.; Levay, Leah J.; Nole, Michael A.
Slow slip events (SSEs) accommodate a significant proportion of tectonic plate motion at subduction zones, yet little is known about the faults that actually host them. The shallow depth (<2 km) of well-documented SSEs at the Hikurangi subduction zone offshore New Zealand offers a unique opportunity to link geophysical imaging of the subduction zone with direct access to incoming material that represents the megathrust fault rocks hosting slow slip. Two recent International Ocean Discovery Program Expeditions sampled this incoming material before it is entrained immediately down-dip along the shallow plate interface. Drilling results, tied to regional seismic reflection images, reveal heterogeneous lithologies with highly variable physical properties entering the SSE source region. These observations suggest that SSEs and associated slow earthquake phenomena are promoted by lithological, mechanical, and frictional heterogeneity within the fault zone, enhanced by geometric complexity associated with subduction of rough crust.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user's guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.
Benacchio, Tommaso; Bonaventura, Luca; Altenbernd, Mirco; Cantwell, Chris D.; Duben, Peter D.; Gillard, Mike; Giraud, Luc; Goddeke, Dominik; Raffin, Erwan; Teranishi, Keita T.; Wedi, Nils
Numerical weather and climate prediction rates as one of the scientific applications whose accuracy improvements greatly depend on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increasing average frequency of hardware and software failures forces users to review their algorithms and systems in order to protect simulations from breakdown. This report surveys approaches for fault-tolerance in numerical algorithms and system resilience in parallel simulations from the perspective of numerical weather and climate prediction systems. A selection of existing strategies is analyzed, featuring interpolation-restart and compressed checkpointing for the numerics, in-memory checkpointing, user-level failure mitigation-based and backup-based methods for the systems. Numerical examples showcase the performance of the techniques in addressing faults, with particular emphasis on iterative solvers for linear systems, a staple of atmospheric fluid flow solvers. The potential impact of these strategies is discussed in relation to current development of numerical weather prediction algorithms and systems towards the exascale. Trade-offs between performance, efficiency and effectiveness of resiliency strategies are analyzed and some recommendations outlined for future developments.
Siera/SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra/SM.
Planck’s law predicts the distribution of radiation energy, color and intensity, emitted from a hot object at thermal equilibrium. The Law also sets the upper limit of radiation intensity, the blackbody limit. Recent experiments reveal that micro-structured tungsten can exhibit significant deviation from the blackbody spectrum. However, whether thermal radiation with weak non-equilibrium pumping can exceed the blackbody limit in the far field remains un-answered experimentally. Here, we compare thermal radiation from a micro-cavity/tungsten photonic crystal (W-PC) and a blackbody, which are both measured from the same sample and also in-situ. We show that thermal radiation can exceed the blackbody limit by >8 times at λ=1.7 μm resonant wavelength in the far-field. Our observation is consistent with a recent calculation by Wang and John performed for a 2D W-PC filament. This finding is attributed to non-equilibrium excitation of localized surface plasmon resonances coupled to nonlinear oscillators and the propagation of the electromagnetic waves through non-linear Bloch waves of the W-PC structure. This discovery could help create super-intense narrow band thermal light sources and even an infrared emitter with a laser-like input-output characteristic.
The attachment of dopant precursor molecules to depassivated areas of hydrogen-terminated silicon templated with a scanning tunneling microscope (STM) has been used to create electronic devices with sub-nanometer precision, typically for quantum physics demonstrations, and to dope silicon past the solid-solubility limit, with potential applications in microelectronics and plasmonics. However, this process, which we call atomic precision advanced manufacturing (APAM), currently lacks the throughput required to develop sophisticated applications because there is no proven scalable hydrogen lithography pathway. Here, we demonstrate and characterize an APAM device workflow where STM lithography has been replaced with photolithography. An ultraviolet laser is shown to locally heat silicon controllably above the temperature required for hydrogen depassivation. STM images indicate a narrow range of laser energy density where hydrogen has been depassivated, and the surface remains well-ordered. A model for photothermal heating of silicon predicts a local temperature which is consistent with atomic-scale STM images of the photo-patterned regions. Finally, a simple device made by exposing photo-depassivated silicon to phosphine is found to have a carrier density and mobility similar to that produced by similar devices patterned by STM.
Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.56 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to Sierra/SD, User's Notes. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer notes manual, the user's notes and of course the material in the open literature.
This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.
We report on the experimental observation of a decreased self-injection threshold by using laser pulses with circular polarization in laser wakefield acceleration experiments in a nonpreformed plasma, compared to the usually employed linear polarization. A significantly higher electron beam charge was also observed for circular polarization compared to linear polarization over a wide range of parameters. Theoretical analysis and quasi-3D particle-in-cell simulations reveal that the self-injection and hence the laser wakefield acceleration is polarization dependent and indicate a different injection mechanism for circularly polarized laser pulses, originating from larger momentum gain by electrons during above threshold ionization. This enables electrons to meet the trapping condition more easily, and the resulting higher plasma temperature was confirmed via spectroscopy of the XUV plasma emission.
Moderate-temperature thermal sources (100° to 400°C) that radiate waste heat are often the by-product of mechanical work, chemical or nuclear reactions, or information processing. We demonstrate conversion of thermal radiation into electrical power using a bipolar grating-coupled complementary metal-oxide-silicon (CMOS) tunnel diode. A two-step photon-assisted tunneling charge pumping mechanism results in separation of charge carriers in pn-junction wells leading to a large open-circuit voltage developed across a load. Electrical power generation from a broadband blackbody thermal source has been experimentally demonstrated with converted power densities of 27 to 61 microwatts per square centimeter for thermal sources between 250° and 400°C. Scalable, efficient conversion of radiated waste heat into electrical power can be used to reduce energy consumption or to power electronics and sensors.
Low-salinity waterflooding (LSWF) has proven to improve oil recovery in carbonate formations through rock wettability alteration, although the underlying mechanism remains elusive. Multivalent ionic exchange and calcite dissolution have usually been investigated using geochemical analysis in secondary coreflooding. In this work, coreflooding, in tertiary mode, coupled with a surface reactivity analysis approach was employed to investigate the interplay of wettability alteration mechanisms such as mineral dissolution, electrostatic bond attraction, and the effect of pH at in situ conditions. Improved oil recovery (IOR) in tertiary mode observed by coreflooding in Indiana limestone rocks showed an ionic strength dependence, that is, reducing brine ionic strength resulted in an increase in oil recovery. Coreflooding results showed that the seawater and low-salinity brines deprived of Mg2+ ions resulted in the lowest IOR in tertiary mode, indicating the significance of Mg2+ on IOR in limestone rocks. Similar results were observed through the contact angle measurement showing the limestone rock wettability state dependence on ionic strength and the effect of Mg2+ ions. Surface reactivity analysis showed an increase in solution pH, Ca2+ and Mg2+ ions concentration in the effluent solution from the coreflooding in tertiary mode using low salinity brines (about 40 and 20% increase in the effluent composition for Ca2+ and Mg2+, respectively). These changes in solution composition were used to calculate the in situ oil-brine and rock-brine zeta potential using a validated surface complexation model, showing the changes of zeta potential as brine is injected into limestone rocks. The results show that using seawater-like brine in tertiary mode resulted in no mineral dissolution or ionic exchange. However, improved oil recovery (IOR) using such seawater-like brine was due to wettability alteration caused by reduced electrostatic bond attraction associated with Mg2+ ions [from 2.6 × 10-13 (mol/m2)2 for formation water salinity to 1.5 × 10-13 (mol/m2)2 for seawater salinity]. Using low-salinity brines in tertiary mode improved oil recovery by mineral dissolution, resulting in oil desorption and an increase in solution pH. The increase in solution pH also resulted in reduced electrostatic bond attraction which lead to rock wettability alteration using low-salinity brines.
Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.
Molten Salt Reactors (MSRs) have seen a resurgence of interest in the past decade around the world. Support for these activities is provided from both national and private sources. The largest difference from the 2011 GIF MSR PR&PP evaluation consequently is the transition from evaluating academic systems focused on exploring the technical potential of MSRs to those of companies and countries focusing on near-term deployment. A wide variety of designs currently exist ranging from solid to liquid-fueled designs, with salt processing on-site or off-site, and a variety of fuel choices. As such, the proliferation resistance and physical protection aspects will have significantly more variation depending on reactor design than the other advanced reactors. The rapid introduction and evolution of innovative MSR designs inevitably means that technology specific details of overview reports, such as this one, become rapidly outdated. Consequently, this report focuses on essential features required for any MSR rather than specific design aspects.