Fu, Pengcheng; Schoenball, Martin; Ajo-Franklin, Jonathan B.; Chai, Chengping; Maceira, Monica; Morris, Joseph P.; Wu, Hui; Knox, Hunter; Schwering, Paul C.; White, Mark D.; Burghardt, Jeffrey A.; Strickland, Christopher E.; Johnson, Timothy C.; Vermeul, Vince R.; Sprinkle, Parker; Roberts, Benjamin; Ulrich, Craig; Guglielmi, Yves; Cook, Paul J.; Dobson, Patrick F.; Wood, Todd; Frash, Luke P.; Ingraham, Mathew D.; Pope, Joseph S.; Smith, Megan M.; Neupane, Ghanashyam; Doe, Thomas W.; Roggenthen, William M.; Horne, Roland; Singh, Ankush; Zoback, Mark D.; Wang, Herb; Condon, Kate; Ghassemi, Ahmad; Chen, Hao; Mcclure, Mark W.; Vandine, George; Blankenship, Douglas A.; Kneafsey, Timothy J.
The final version of the above article was posted prematurely on 16 July 2021, owing to a technical error. The final, corrected version of record will be made fully available at a later date.
Direct Numerical Simulations (DNS) are performed to investigate the process of spontaneous ignition of hydrogen flames at laminar, turbulent, adiabatic and non-adiabatic conditions. Mixtures of hydrogen and vitiated air at temperatures representing gas-turbine reheat combustion are considered. Adiabatic spontaneous ignition processes are investigated first, providing a quantitative characterization of stable and unstable flames. Results indicate that, in hydrogen reheat combustion, compressibility effects play a key role in flame stability and that unstable ignition and combustion are consistently encountered for reactant temperatures close to the mixture's characteristic crossover temperature. Furthermore, it is also found that the characterization of the adiabatic processes is also valid in the presence of non-adiabaticity due to wall heat-loss. Finally, a quantitative characterization of the instantaneous fuel consumption rate within the reaction front is obtained and of its ability, at auto-ignitive conditions, to advance against the approaching turbulent flow of the reactants, for a range of different turbulence intensities, temperatures and pressure levels.
Laplace, T.A.; Goldblum, B.L.; Manfredi, J.J.; Brown, J.A.; Bleuel, D.L.; Brand, C.A.; Gabella, G.; Gordon, J.; Brubaker, Erik B.
Background: Organic scintillators are widely used for neutron detection in both basic nuclear physics and applications. While the proton light yield of organic scintillators has been extensively studied, measurements of the light yield from neutron interactions with carbon nuclei are scarce. Purpose: Demonstrate a new approach for the simultaneous measurement of the proton and carbon light yield of organic scintillators. Provide new carbon light yield data for the EJ-309 liquid and EJ-204 plastic organic scintillators. Method: A 33-MeV H+2 beam from the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory was impinged upon a 3-mm-thick Be target to produce a high-flux, broad-spectrum neutron beam. The double time-of-flight technique was extended to simultaneously measure the proton and carbon light yields of the organic scintillators, wherein the light output associated with the recoil particle was determined using np and nC elastic scattering kinematics. Results: The proton and carbon light yield relations of the EJ-309 liquid and EJ-204 plastic organic scintillators were measured over a recoil energy range of approximately 0.3 to 1 MeV and 2 to 5 MeV, respectively, for EJ-309, and 0.2 to 0.5 MeV and 1 to 4 MeV, respectively, for EJ-204. Conclusions: These data provide new insight into the ionization quenching effect in organic scintillators and key input for simulation of the response of organic scintillators for both basic science and a broad range of applications.
International Journal of High Performance Computing Applications
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
Progress in numerical weather and climate prediction accuracy greatly depends on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increased 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 hardware, application-level and algorithm-level resilience approaches of particular relevance to time-critical numerical weather and climate prediction systems. A selection of applicable existing strategies is analysed, featuring interpolation-restart and compressed checkpointing for the numerical schemes, in-memory checkpointing, user-level failure mitigation 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 analysed and some recommendations outlined for future developments.
Two events of magnitude (mb) 3.6-3.8 occurred in southern North Korea (NK) on 27 June 2019 and 11 May 2020. Although these events were located ~330-400 km from the known nuclear test site, the fact that they occurred within the territory of NK, a country with a recent history of underground nuclear tests, made them events of interest for the monitoring community. Weused P/Lg ratios from regional stations to categorize seismic events that occurred in NK from 2006 to May 2020, including these two recent events, the six declared NK nuclear tests, and the cavity collapse and triggered earthquakes that followed the 3 September 2017 nuclear explosion. We were able to separate the cavity collapse from the population of nuclear explosions. However, based on P/Lg ratios, the distinction between the earthquakes and the cavity collapse is ambiguous. The performed discriminant analyses suggest that combining Pg/Lg and Pn/Lg ratios results in improved discriminant power compared with any of the ratio types alone. We used the two ratio types jointly in a quadratic discriminant function and successfully classified the six declared nuclear tests and the triggered earthquakes that followed the September 2017 explosion. Our analyses also confirm that the recent southern events of June 2019 and May 2020 are both tectonic earthquakes that occurred naturally.
This report summarizes the international collaboration work conducted by Sandia and funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-20SN010303062. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), and model comparison (DECOVALEX). Lastly, the report summarizes a newly developed working group on the development of scenarios as part of the performance assessment development process, and the activities related to the Nuclear Energy Agency (NEA) Salt club and the US/German Workshop on Repository Research, Design and Operations.
Natural and anthropogenic infrasound may travel vast distances, making it an invaluable resource for monitoring phenomena such as nuclear explosions, volcanic eruptions, severe storms, and many others. Typically, these waves are captured using pressure sensors, which cannot encode the direction of arrival—critical information when the source location is not known beforehand. Obtaining this information therefore requires arrays of sensors with apertures ranging from tens of meters to kilometers depending on the wavelengths of interest. This is often impractical in locations that lack the necessary real estate (urban areas, rugged regions, or remote islands); in any case, it requires multiple power, digitizer, and telemetry deployments. Here, the theoretical basis behind a compact infrasound direction of arrival sensor based on the acoustic metamaterials is presented. This sensor occupies a footprint that is orders of magnitude smaller than the span of a typical infrasound array. The diminutive size of the unit greatly expands the locations where it can be deployed. The sensor design is described, its ability to determine the direction of arrival is evaluated, and further avenues of study are suggested.
The software package developed by Sandia National Laboratories is intended to allow the integration of Simulink models into emulations of control networks. To accomplish this, three programs are included: Simulink S-Function, Data Broker, and End Point
We present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current-voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia's Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn's Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16 times 10 -{10} and 1.88 times 10 -{10} rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV's non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.
Despite their wide use in terahertz (THz) research and technology, the application spectra of photoconductive antenna (PCA) THz detectors are severely limited due to the relatively high optical gating power requirement. This originates from poor conversion efficiency of optical gate beam photons to photocurrent in materials with subpicosecond carrier lifetimes. Here we show that using an ultra-thin (160 nm), perfectly absorbing low-temperature grown GaAs metasurface as the photoconductive channel drastically improves the efficiency of THz PCA detectors. This is achieved through perfect absorption of the gate beam in a significantly reduced photoconductive volume, enabled by the metasurface. This Letter demonstrates that sensitive THz PCA detection is possible using optical gate powers as low as 5 μW-three orders of magnitude lower than gating powers used for conventionalPCAdetectors.We show that significantly higher optical gate powers are not necessary for optimal operation, as they do not improve the sensitivity to the THz field. This class of efficient PCA THz detectors opens doors for THz applications with low gate power requirements.
Accelerated growth of the additive manufacturing (AM) industry in recent years is accompanied by a rising need for methods to quickly assess quality at-scale. Current practices for quality inspection include nondestructive test methods and destructive testing of witness coupons, which are artifacts built alongside the actual part. However, these methods can be costly and time-consuming. Recognizing this need, the Additive Manufacturing Center of Excellence (AM CoE) initiated a project led by its partner, Auburn University, to develop rapid testing procedure using asbuilt samples tested in torsion to quantitatively assess build quality. The presented work developed a rapid testing procedure using as-built samples tested in torsion to quantify small variances for assessing build quality.
Nagasaka, Cocoro A.; Kozma, Karoly; Russo, Chris J.; Alam, Todd M.; Nyman, May
Removal of radioactive Cs from sodium-rich solutions is a technical challenge that goes back to post World War II nuclear waste storage and treatment; and interest in this topic was reinvigorated by the Fukushima-Daiichi nuclear power plant disaster, 10 years ago. Since the 1960′s there has been considerable focus on layered Zr phosphates as robust inorganic sorbents for separation of radionuclides such as Cs. Here we present synthesis and characterization, and direct comparison of Cs sorption capacity and selectivity of four related materials: 1) crystalline α-Zr phosphate and α-Hf phosphate, and 2) amorphous analogues of these. Powder X-ray diffraction, thermogravimetry, solid-state 31P magic angle spinning nuclear magnetic resonance (MAS-NMR) spectroscopy, and compositional analysis (inductively coupled plasma optical emission spectroscopy and mass spectroscopy, ICP OES and ICP MS) provided formulae; respectively M(HPO4)2⋅1H2O and M(HPO4)2⋅4H2O (M = Hf, Zr) for crystalline and amorphous analogues. Maximum Cs loading, competitive Cs-Na selectivity and maximum Cs-Na loading followed by the above characterizations plus 133Cs MAS-NMR spectroscopy revealed that amorphous analogues are considerably better Cs-sorbents (based on maximum Cs-loading and selectivity over Na) than the well-studied crystalline Zr-analogue. Additionally, crystalline α-Hf phosphate is better Cs-sorbent than crystalline α-Zr phosphate. All these studies consistently show that Hf phosphate is less crystallize than Zr phosphate, when obtained under similar or identical synthesis conditions. We attribute this to lower solubility of Hf phosphate compared to Zr phosphate, preventing ‘defect healing’ during the synthesis process.
This project was a follow-on to the Sandia National Laboratories (SNL) and the Laboratory for Laser Energetics (LLE) ARPA-E ALPHA project entitled “Demonstrating Fuel Magnetization and Laser Heating Tools for Low-Cost Fusion Energy”. The primary purpose of this follow-on project was to obtain additional data at the OMEGA facility to help better understand how MagLIF, a platform that has already demonstrated the scientific viability of magneto-inertial fusion, scales across a factor of 1000 in driver energy. A secondary aspect of this project was to extend simulations and analysis at SNL to cover a wider magneto-inertial fusion (MIF) parameter space and test scaling of those models across this wide range of input energies and conditions of the target. This work was successful in improving understanding of how key physics elements of MIF scales and improves confidence in setting requirements for fusion gain with larger drivers. The OMEGA experiments at the smaller scale verified the hypothesis that preheating the fuel plays a significant role in introducing wall contaminants that mix into the fuel and significantly degrade fusion performance. This contamination not only impacts target performance but the optimal input conditions for the target. However, analysis at the Z-scale showed that target performance at high preheat levels is limited by the Nernst effect, which advects magnetic flux from the hot spot, reducing magnetic insulation and consequently reduces the temperature of the fuel. The combination of MagLIF experiments at the disparate scales of OMEGA and Z along with a multiscale 3D simulation analysis has led to new insight into the physical mechanisms responsible for limiting target performance and provides important benchmarks to assess target scaling more generally for MIF schemes. Finally, in addition to the MagLIF related work, a semi-analytic model of liner driven Field Reversed Configuration (FRC) was developed that predicts the fusion gain for such systems. This model was also validated with 2D radiation magneto-hydrodynamic simulations and predicts that fusion gains of near unity could be driven by the Z machine.