Publications

Abstract not provided.

Abstract not provided.

This paper presents a nonlinear control design technique that capitalizes on an hour glass (HG) variable geometry wave energy converter (WEC). The HG buoy is assumed to operate in the heave motion of the wave. The unique interaction between the HG buoy and the wave creates a nonlinear cubic storage effect that produces actual energy storage or reactive power during operation. A multi-frequency Bretschneider spectrum wave excitation input is reviewed for the HG design both with constant and varying steepness angle profiles which demonstrates further increased power generation. Numerical simulations are performed to demonstrate the increase in power generation with changing sea states. The objective is to increase the power generation from multi-frequency nonlinear dynamic sources.

We evaluate the use of reference modules for monitoring effective irradiance in PV power plants, as compared with traditional plane-of-array (POA) irradiance sensors, for PV monitoring and capacity tests. Common POA sensors such as pyranometers and reference cells are unable to capture module-level irradiance nonuniformity and require several correction factors to accurately represent the conditions for fielded modules. These problems are compounded for bifacial systems, where the power loss due to rear side shading and rear-side plane-of-array (RPOA) irradiance gradients are greater and more difficult to quantify. The resulting inaccuracy can have costly real-world consequences, particularly when the data are used to perform power ratings and capacity tests. Here we analyze data from a bifacial single-axis tracking PV power plant, (175.6 MWdc) using 5 meteorological (MET) stations, located on corresponding inverter blocks with capacities over 4 MWdc. Each MET station consists of bifacial reference modules as well pyranometers mounted in traditional POA and RPOA installations across the PV power plant. Short circuit current measurements of the reference modules are converted to effective irradiance with temperature correction and scaling based on flash test or nameplate short circuit values. Our work shows that bifacial effective irradiance measured by pyranometers averages 3.6% higher than the effective irradiance measured by bifacial reference modules, even when accounting for spectral, angle of incidence, and irradiance nonuniformity. We also performed capacity tests using effective irradiance measured by pyranometers and reference modules for each of the 5 bifacial single-axis tracking inverter blocks mentioned above. These capacity tests evaluated bifacial plant performance at ∼3.9% lower when using bifacial effective irradiance from pyranometers as compared to the same calculation performed with reference modules.

Sedimentary-hosted geothermal energy systems are permeable structural, structural-stratigraphic, and/or stratigraphic horizons with sufficient temperature for direct use and/or electricity generation. Sedimentary-hosted (i.e., stratigraphic) geothermal reservoirs may be present in multiple locations across the central and eastern Great Basin of the USA, thereby constituting a potentially large base of untapped, economically accessible energy resources. Sandia National Laboratories has partnered with a multi-disciplinary group of collaborators to evaluate a stratigraphic system in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The goal of this study is to inform an optimized strategy for subsequent exploration and development of this and analogous resources. Building from prior Nevada Play Fairway Analysis (PFA), this team is primarily 1) collecting additional geophysical data, 2) employing novel joint geophysical inversion/modeling techniques to update existing 3D geologic models, and 3) integrating the geophysical results to produce a working, geologically constrained thermo-hydrological reservoir model. Prior PFA work highlights Steptoe Valley as a favorable resource basin that likely has both sedimentary and hydrothermal characteristics. However, there remains significant uncertainty on the nature and architecture of the resource(s) at depth, which increases the risk in exploratory drilling. Newly acquired gravity, magnetic, magnetotelluric, and controlled-source electromagnetic data, in conjunction with new and preceding geoscientific measurements and observations, are being integrated and evaluated in this study for efficacy in understanding stratigraphic geothermal resources and mitigating exploration risk. Furthermore, the influence of hydrothermal activity on sedimentary-hosted reservoirs in favorable structural settings (i.e., whether fault-controlled systems may locally enhance temperature and permeability in some deep stratigraphic reservoirs) will also be evaluated. This paper provides details and current updates on the course of this study in-progress.

Abstract not provided.

Density fluctuations in compressible turbulent boundary layers cause aero-optical distortions that affect the performance of optical systems such as sensors and lasers. The development of models for predicting the aero-optical distortions relies on theory and reference data that can be obtained from experiments and time-resolved simulations. This paper reports on wall-modeled large-eddy simulations of turbulent boundary layers over a flat plate at Mach 3.5, 7.87, and 13.64. The conditions for the Mach 3.5 case match those for the DNS presented by Miller et al.1 The Mach 7.87 simulation match those inside the Hypersonic Wind Tunnel at Sandia National Laboratories. For the Mach 13.64, the conditions inside the Arnold Engineering Development Complex Hypervelocity Tunnel 9 are matched. Overall, adequate agreement of the velocity and temperature as well as Reynolds stress profiles with reference data from direct numerical simulations is obtained for the different Mach numbers. For all three cases, the normalized root-mean-square optical path difference was computed and compared with data obtained from the reference direct numerical simulations and experiments, as well as predictions obtained with a semi-analytical relationship by Notre Dame University. Above Mach five, the normalized path difference obtained from the simulations is above the model prediction. This provides motivation for future work aimed at evaluating the assumptions behind the Notre Dame model for hypersonic boundary layer flows.

Decarbonizing natural gas networks is a challenging enterprise. Replacing natural gas with renewable hydrogen is one option under global consideration to decarbonize heating, power and residential uses of natural gas. Hydrogen is known to degrade fatigue and fracture properties of structural steels, including pipeline steels. In this study, we describe environmental testing strategies aimed at generating baseline fatigue and fracture trends with efficient use of testing resources. For example, by controlling the stress intensity factor (K) in both K-increasing and K-decreasing modes, fatigue crack growth can be measured for multiple load ratios with a single specimen. Additionally, tests can be designed such that fracture tests can be performed at the conclusion of the fatigue crack growth test, further reducing the resources needed to evaluate the fracture mechanics parameters utilized in design. These testing strategies are employed to establish the fatigue crack growth behavior and fracture resistance of API grade steels in gaseous hydrogen environments. In particular, we explore the effects of load ratio and hydrogen partial pressure on the baseline fatigue and fracture trends of line pipe steels in gaseous hydrogen. These data are then used to test the applicability of a simple, universal fatigue crack growth model that accounts for both load ratio and hydrogen partial pressure. The appropriateness of this model for use as an upper bound fatigue crack growth is discussed.

It is impossible in practice to comprehensively test even small software programs due to the vastness of the reachable state space; however, modern cyber-physical systems such as aircraft require a high degree of confidence in software safety and reliability. Here we explore methods of generating test sets to effectively and efficiently explore the state space for a module based on the Traffic Collision Avoidance System (TCAS) used on commercial aircraft. A formal model of TCAS in the model-checking language NuSMV provides an output oracle. We compare test sets generated using various methods, including covering arrays, random, and a low-complexity input paradigm applied to 28 versions of the TCAS C program containing seeded errors. Faults are triggered by tests for all 28 programs using a combination of covering arrays and random input generation. Complexity-based inputs perform more efficiently than covering arrays, and can be paired with random input generation to create efficient and effective test sets. A random forest classifier identifies variable values that can be targeted to generate tests even more efficiently in future work, by combining a machine-learned fuzzing algorithm with more complex model oracles developed in model-based systems engineering (MBSE) software.

Abstract not provided.

The postclosure criticality safety assessment for the direct disposal of dual-purpose canisters (DPCs) in a geologic repository includes considerations of transient criticality phenomena. The power pulse from a hypothetical transient criticality event in an unsaturated alluvial repository is evaluated for a DPC containing 37 spent pressurized water reactor (PWR) assemblies. The scenario assumes that the conditions for baseline criticality are achieved through flooding with groundwater and progressive failure of neutron absorbing media. A preliminary series of steady-state criticality calculations is conducted to characterize reactivity feedback due to absorber degradation, Doppler broadening, and thermal expansion. These feedback coefficients are used in an analysis with a reactor kinetics code to characterize the transient pulse given a positive reactivity insertion for a given length of time. The time-integrated behavior of the pulse can be used to model effects on the DPC and surrounding barriers in future studies and determine if transient criticality effects are consequential.

A survey of cadmium plated field return hardware showed ubiquitous cadmium whisker growth. The most worn and debris-covered hardware showed the densest whisker growth. Whiskers were often found growing in agglomerates of nodules and whiskers. The hardware was rinsed with alcohol to transfer whiskers and debris from the hardware to a flat stub. Fifty whiskers were studied individually by scanning electron microscopy (SEM), including energy dispersive spectroscopy (EDS) and electron backscatter diffraction (EBSD). Most of the whiskers were single crystal, though three were found to contain grain boundaries at kinks. The whiskers ranged from 5 to 600 μm in length and 80 pct showed a <1 ¯ 2 1 ¯ 0> type growth direction. This growth direction facilitates the development of low energy side faces of the whisker, (0001) and {1010}.

Abstract not provided.

A new small-scale pressure vessel with a 5×5 fuel assembly and axially truncated PWR hardware was created to simulate commercial vacuum drying processes. This test assembly, known as the Dashpot Drying Apparatus, was built to focus on the drying of a single PWR dashpot and surrounding fuel. Drying operations were simulated for three tests with the DDA based on the pressure and temperature histories observed in the HBDP. All three tests were conducted with an empty guide tube. One test was performed with deionized water as the fill fluid. The other two tests used 0.2 M boric acid as the fill fluid to accurately simulate spent fuel pool conditions. These tests proved the capability of the DDA to mimic commercial drying processes on a limited scale and detect the presence of bulk and residual water. Furthermore, for all tests, pressure remained below the 0.4 kPa (3 Torr) rebound threshold for the final evacuation step in the drying procedure. Results indicate that after bulk fluid is removed from the pressure vessel, residual water is verifiably measured through confirmatory measurements of pressure and water content using a mass spectrometer. The final pressure rebound behaviors for the three tests conducted were well below the established regulatory limit of less than 0.4 kPa (3 Torr) within 30 minutes of isolation. The water content measurements across all tests showed that despite observing high water content within the DDA vessel at the beginning of the vacuum isolations, the water content drastically drops to below 1,200 ppmv after the isolations were conducted. The data and operational experience from these tests will guide the next evolution of experiments on a prototypic-length scale with multiple surrogate rods in a full 17×17 PWR assembly. The insight gained through these investigations is expected to support the technical basis for the continued safe storage of spent nuclear fuel into long term operations.

Abstract not provided.

There is a need to perform offline anomaly detection in count data streams to simultaneously identify both systemic changes and outliers, simultaneously. We propose a new algorithmic method, called the Anomaly Detection Pipeline, which leverages common statistical process control procedures in a novel way to accomplish this. The method we propose does not require user-defined control or phase I training data, automatically identifying regions of stability for improved parameter estimation to support change point detection. The method does not require data to be normally distributed, and it detects outliers relative to the regimes in which they occur. Our proposed method performs comparably to state-of-the-art change point detection methods, provides additional capabilities, and is extendable to a larger set of possible data streams than known methods.

This paper focuses on the development and testing of spoofing detection and localization techniques that rely only on clock deviations to identify threat signals. Detection methods that rely on dynamic receiver geometries to triangulate threat locations or signal geometry to identify spoofing are not considered here. Instead this paper focuses on single antenna receivers and assumes the receiver tracks only the inauthentic signal. The quality of the receiver clock has a significant impact on the performance of the receiver tracking loops. Low quality clocks have frequency instabilities that inherently limit the sensitivity of the receiver to slow growing errors. Some clocks provide better frequency stabilities but have a higher white frequency noise that can induce false detections. Because of these trends, various detection methods are tested with four types of receiver and transmitter clocks of varying quality.

Abstract not provided.

This paper reports on Catalyst usability and initial adoption by SPARC analysts. The use case approach highlights the analysts' perspective. Impediments to adoption can be due to deficiencies in software capabilities, or analysts may identify mundane inconveniences and barriers that prevent them from fully leveraging Catalyst. With that said, for many analyst tasks Catalyst provides enough relative advantage that they have begun applying it in their production work, and they recognize the potential for it to solve problems they currently struggle with. The findings in this report include specific issues and minor bugs in ParaView Python scripting, which are viewed as having straightforward solutions, as well as a broader adoption analysis.

Abstract not provided.

The emerging use of the physics-based athermal recombination-corrected displacement per atom (arc-dpa) model for the displacement damage efficiency has motivated a re-evaluation of the historical empirically-derived GaAs damage response function with the purpose of highlighting needs for future analytical and experimental work. The 1-MeV neutron damage equivalence methodology used in the ASTM E-722 standard for GaAs has been re-evaluated using updated nuclear data. This yielded a higher fidelity representation of the GaAs displacement kerma and, through the use of the refined PKA recoil energy-dependent damage efficiency model, an updated 1-MeV(GaAs) displacement damage function. This re-evaluation included use of the Norgett-Robinson-Torrens (NRT) model for an updated threshold treatment, rather than the sharp-threshold Kinchin-Pease model used in the current ASTM standard. The underlying nuclear data evaluations have been updated to use the ENDF/VIII.0 ⁷⁵As and TENDL-2019 ⁷¹Ga/⁶⁹Ga evaluations. The displacement kerma and 1-MeV-equivalent damage responses were calculated using a modified NJOY-2016 code which allowed for refinements in some of the damage models. This paper shows that an updated displacement damage function, based upon the latest nuclear data, is consistent with the experimental data used to develop the current ASTM E-722 GaAs standard. Using a double ratio approach to compare the available experimental data with the calculated response, the average legacy double ratio was found to be 0.97±0.05 and the average updated double ratio was found to be 0.94 ±0.05.

A crucial component of field testing is the utilization of numerical models to better understand the system and the experimental data being collected. Meshing and modeling field tests is a complex and computationally demanding problem. Hexahedral elements cannot always reproduce experimental dimensions leading to grid orientation or geometric errors. Voronoi meshes can match complex geometries without sacrificing orthogonality. As a result, here we present a high-resolution 3D numerical study for the BATS heater test at the WIPP that compares both a standard non-deformed cartesian mesh along with a Voronoi mesh to match field data collected during a salt heater experiment.

Using neural networks to solve variational problems, and other scientific machine learning tasks, has been limited by a lack of consistency and an inability to exactly integrate expressions involving neural network architectures. We address these limitations by formulating a polynomial-spline network, a novel shallow multilinear perceptron (MLP) architecture incorporating free knot B-spline basis functions into a polynomial mixture-of-experts model. Effectively, our architecture performs piecewise polynomial approximation on each cell of a trainable partition of unity while ensuring the MLP and its derivatives can be integrated exactly, obviating a reliance on sampling or quadrature and enabling error-free computation of variational forms. We demonstrate hp-convergence for regression problems at convergence rates expected from approximation theory and solve elliptic problems in one and two dimensions, with a favorable comparison to adaptive finite elements.

Abstract not provided.

This paper applies sensitivity and uncertainty analysis to compare two model alternatives for fuel matrix degradation for performance assessment of a generic crystalline repository. The results show that this model choice has little effect on uncertainty in the peak 129I concentration. The small impact of this choice is likely due to the higher importance of uncertainty in the instantaneous release fraction and differences in epistemic uncertainty between the alternatives.

We present a distributed Brillouin fiber sensor that operates by exciting a series of discrete lasing modes. This approach provides inherently wide dynamic range (5m) while the narrow linewidth lasing modes enable low noise (8n/Hz)

Abstract not provided.

Abstract not provided.

Abstract not provided.

We computationally explore the optical and elastic modes necessary for acoustoelectrically enhanced Brillouin interactions. The large simulated piezoelectric (k2 « 6%) and optome-chanica (lgol ≈ 8000 (rad/s)√m) coupling theoretically predicts a performance enhancement several order of magniude in brrloum-based photoine tecnnologres.

Femtosecond laser electronic excitation tagging (FLEET) is a powerful unseeded velocimetry technique typically used to measure one component of velocity along a line, or two or three components from a dot. In this Letter, we demonstrate a dotted-line FLEET technique which combines the dense profile capability of a line with the ability to perform two-component velocimetry with a single camera on a dot. Our set-up uses a single beam path to create multiple simultaneous spots, more than previously achieved in other FLEET spot configurations. We perform dotted-line FLEET measurements downstream of a highly turbulent, supersonic nitrogen free jet. Dotted-line FLEET is created by focusing light transmitted by a periodic mask with rectangular slits of 1.6 × 40 mm2 and an edge-to-edge spacing of 0.5 mm, then focusing the imaged light at the measurement region. Up to seven symmetric dots spaced approximately 0.9 mm apart, with mean full-width at half maximum diameters between 150 and 350 µm, are simultaneously imaged. Both streamwise and radial velocities are computed and presented in this Letter.

Propagating thermal runaway events are a significant threat to utility-scale storage installations. A propagating thermal runaway event is a cascading series of failures in which energy released from a failed cell triggers subsequent failures in nearby cells. Without intervention, propagation can turn an otherwise manageable single cell failure into a full system conflagration. This study presents a method of mitigating the severity of propagating thermal runaway events in utility-scale storage systems by leveraging the capabilities of a module-interfaced power conversion architecture. The method involves strategic depletion of storage modules to delay or arrest propagation, reducing the total thermal energy released in the failure event. The feasibility of the method is assessed through simulations of propagating thermal runaway events in a 160 kW/80 kWh energy storage system.

Traditional deep learning (DL) models are powerful classifiers, but many approaches do not provide uncertainties for their estimates. Uncertainty quantification (UQ) methods for DL models have received increased attention in the literature due to their usefulness in decision making, particularly for high-consequence decisions. However, there has been little research done on how to evaluate the quality of such methods. We use statistical methods of frequentist interval coverage and interval width to evaluate the quality of credible intervals, and expected calibration error to evaluate classification predicted confidence. These metrics are evaluated on Bayesian neural networks (BNN) fit using Markov Chain Monte Carlo (MCMC) and variational inference (VI), bootstrapped neural networks (NN), Deep Ensembles (DE), and Monte Carlo (MC) dropout. We apply these different UQ for DL methods to a hyperspectral image target detection problem and show the inconsistency of the different methods' results and the necessity of a UQ quality metric. To reconcile these differences and choose a UQ method that appropriately quantifies the uncertainty, we create a simulated data set with fully parameterized probability distribution for a two-class classification problem. The gold standard MCMC performs the best overall, and the bootstrapped NN is a close second, requiring the same computational expense as DE. Through this comparison, we demonstrate that, for a given data set, different models can produce uncertainty estimates of markedly different quality. This in turn points to a great need for principled assessment methods of UQ quality in DL applications.

Long running applications often have varying behaviors, here called phases. While considerable work in computer architecture has been done in identifying application phases based on how the hardware is being exercised, comparatively less work has been focused on identifying application phases based on regions of source code being executed. In this paper we introduce a new methodology and an efficient tool framework, IncProf, for observing and capturing the time-varying source execution behavior of applications, and for then deducing application phases from the resulting data. Uses of this capability include simply better understanding the varying behavior of long running applications, and for efficiently tracking deployed application performance in the future by providing information to identify good instrumentation points.

The performance of APCs with relatively low compressive strength and poor stability under hydrothermal conditions make them less than desirable for the DPC use case. Meanwhile, grossite as a primary filler material or as a modifier has resulted in marked improvements in the properties of several DPC cement filler candidates. Grossite CAPCs have substantial mechanical strength even after irradiation. However, the significant decrease in strength observed post-irradiation requires further investigation before it is advanced as a material for the use case. As a modifier, grossite improves strength and set times of the APC and WAPC cements. Hibonite CAPCs also show considerable promise although their degradation under hydrothermal conditions is a potentially significant liability. Finally, with recent improvements in working time and compressive strength, the WAPCs remain in contention as viable candidates for the DPC use case.

Hardware design in high-performance computing (HPC) is often highly experimental. Exploring new designs is difficult and time-consuming, requiring lengthy vendor cooperation. RISC-V is an open-source processor ISA that improves the accessibility of chip design, including the ability to do hardware/software co-design using open-source hardware and tools. Co-design allows design decisions to easily flow across the hardware/software boundary and influence future design ideas. However, new hardware designs require corresponding software to drive and test them. Conventional operating systems like Linux are massively complex and modification is time-prohibitive. In this paper, we describe our port of the Kitten lightweight kernel operating system to RISC-V in order to provide an alternative to Linux for conducting co-design research. Kitten's small code base and simple resource management policies are well matched for quickly exploring new hardware ideas that may require radical operating system modifications and restructuring. Our evaluation shows that Kitten on RISC-V is functional and provides similar performance to Linux for single-core benchmarks. This provides a solid foundation for using Kitten in future co-design research involving RISC-V.

Fabric-attached memory (FAM) is proposed to enable the seamless integration of directly accessible memory modules attached to the shared system fabric, which will provide future systems with flexible memory integration options, mitigate underutilization, and facilitate data sharing. Recently proposed interconnects, such as Gen-Z and Compute Express Link (CXL), define security, correctness, and performance requirements of fabric-attached devices, including memory. These initiatives are supported by most major system and processor vendors, bringing widespread adoption of FAM-enabled systems one step closer to reality and security concerns to the forefront. This paper discusses the challenges for adapting secure memory implementations to FAM-enabled systems for the first time in literature. Specifically, we observe that handling the security metadata used to protect fabric-attached memories needs to be done deliberately to eliminate unintentional integrity check failures and/or security vulnerabilities, caused by an inconsistent view of the shared security metadata across nodes. Our scheme, Minerva, elegantly adapts secure memory implementations to support FAM-enabled systems with negligible performance over-heads (3.8% of an ideal scheme), compared to the performance overhead (99.5% of an ideal scheme) for a scheme that uses conventional invalidation-based cache coherence to ensure the consistency of security metadata across nodes.

Transient stability control of power systems is based on actions that are taken automatically following a disturbance to ensure that the system remains in synchronism. Examples of such measures include generator rejection and the insertion of dynamic braking resistors. Methods like these are designed to rapidly absorb excess energy or otherwise alter the generation-demand balance at key points in the system. While these methods are often effective, they lack the ability to inject real power to compensate for a deficit. Utility-scale inverter-based resources, particularly energy storage systems, enable bidirectional modulation of real power with the bandwidth necessary to provide synchronizing torque. These resources, and the control strategies they enable, have garnered substantial research interest. This paper provides a critical review of research on real power modulation strategies for transient stability control. The design of these control strategies is heavily informed by the methods used to assess changes in the transient stability margins. Rigorously assessing these changes is difficult because the dynamics of large-scale power systems are inherently nonlinear. The well-known equal-area criterion is physically intuitive, but conceptual extensions are necessary for multi-machine systems. So-called direct methods of transient stability analysis offer a more general alternative; however, these methods require many simplifying assumptions and have difficulty incorporating detailed system dynamics. In this paper, we discuss data-driven methods for offline stability assessment based on Koopman operator theory.

Metal additive manufacturing allows for the fabrication of parts at the point of use as well as the manufacture of parts with complex geometries that would be difficult to manufacture via conventional methods (milling, casting, etc.). Additively manufactured parts are likely to contain internal defects due to the melt pool, powder material, and laser velocity conditions when printing. Two different types of defects were present in the CT scans of printed AlSi10Mg dogbones: spherical porosity and irregular porosity. Identification of these pores via a machine learning approach (i.e., support vector machines, convolutional neural networks, k-nearest neighbors’ classifiers) could be helpful with part qualification and inspections. The machine learning approach will aim to label the regions of porosity and label the type of porosity present. The results showed that a combination approach of Canny edge detection and a classification-based machine learning model (k-nearest neighbors or support vector machine) outperformed the convolutional neural network in segmenting and labeling different types of porosity.

In this paper we introduce and study the STREAMINGCYCLES problem, a random order streaming version of the Boolean Hidden Hypermatching problem that has been instrumental in streaming lower bounds over the past decade. In this problem the edges of a graph G, comprising n/l disjoint length-l cycles on n vertices, are partitioned randomly among n players. Every edge is annotated with an independent uniformly random bit, and the players' task is to output, for some cycle in G, the sum (modulo 2) of the bits on its edges, after one round of sequential communication.Our main result is an lO(l) lower bound on the communication complexity of STREAMINGCYCLES, which is tight up to constant factors in the exponent. Applications of our lower bound for STREAMINGCYCLES include an essentially tight lower bound for component collection in (almost) random order graph streams, making progress towards a conjecture of Peng and Sohler [SODA'18] and the first exponential space lower bounds for random walk generation.

Transportation of sodium-bonded spent fuel appears to present no unique challenges. Storage systems for this fuel should be designed to keep water, both liquid and vapor, from contacting the spent fuel. This fuel is not suitable for geologic disposal; therefore, how the spent sodium bonded fuel will be processed and the characteristics of the final disposal waste form(s) need to be considered. TRISO spent fuel appears to present no unique challenges in terms of transportation, storage, or disposal. If the graphite block is disposed of with the TRISO spent fuel, the ¹⁴C and ³H generated would need to be considered in the postclosure performance assessment. Salt waste from the molten salt reactor has yet to be transported or stored and might be a challenge to dispose of in a non-salt repositories. Like sodium-bonded spent fuel, how the salt will be treated and the characteristics of the final disposal waste form(s) need to be considered. In addition, radiolysis in the frozen salt waste form continues to generate gas, which presents a hazard. Both HALEU and high-enriched uranium SNF are currently being stored and transported by the DOE. Disposal of fuels with enrichments greater than 5% was included in the disposal plan for Yucca Mountain. The increased potential for criticality associated with the higher enriched SNF is mitigated by additional criticality control measures. Fuels that are similar to some ATFs were part of the disposal plan for Yucca Mountain. Some of the properties of these fuels (swelling, generation of ¹⁴C) would have to be considered as part of a postclosure performance assessment.

OpenMP 5.0 added support for reductions over explicit tasks. This expands the previous reduction support that was limited primarily to worksharing and parallel constructs. While the scope of a reduction operation in a worksharing construct is the scope of the construct itself, the scope of a task reduction can vary. This difference requires syntactical means to define the scope of reductions, e.g., the task_reduction clause, and to associate participating tasks, e.g., the in_reduction clause. Furthermore, the disassociation of the number of threads and the number of tasks creates space for different implementations in the OpenMP runtime. In this work, we provide insights into the behavior and performance of task reduction implementations in GCC/g++ and LLVM/Clang. Our results indicate that task reductions are well supported by both compilers, but their performance differs in some cases and is often determined by the efficiency of the underlying task management.

This paper presents a run-to-run (R2R) controller for mechanical serial sectioning (MSS). MSS is a destructive material analysis process which repeatedly removes a thin layer of material and images the exposed surface. The images are then used to gain insight into the material properties and often to construct a 3-dimensional reconstruction of the material sample. Currently, an experience human operator selects the parameters of the MSS to achieve the desired thickness. The proposed R2R controller will automate this process while improving the precision of the material removal. The proposed R2R controller solves an optimization problem designed to minimize the variance of the material removal subject to achieving the expected target removal. This optimization problem was embedded in an R2R framework to provide iterative feedback for disturbance rejection and convergence to the target removal amount. Since an analytic model of the MSS system is unavailable, we adopted a data-driven approach to synthesize our R2R controller from historical data. The proposed R2R controller is demonstrated through simulations. Future work will empirically demonstrate the proposed R2R through experiments with a real MSS system.

This paper reports the experimental comparison of two silicon photomultipliers (SiPMs): the MicroFJ-30035 by ONSemi and the ASD-NUV3S-P by AdvanSiD, in terms of gain, dark count rate, and crosstalk probability. SiPMs are solid state photon detectors that enable high sensitivity light readout. They have low-voltage power requirements, small form factor, and are durable. For these reasons, they are being considered as replacements for vacuum photomultiplier tubes in some applications. However, their performance relies on several parameters, which need to be carefully characterized to enable their high-fidelity simulation and SiPM-based design of devices capable to operate in harsh environments. The parameters tend to vary between manufacturers and processing technologies. In this work, we have compared the MicroFJ and ASD SiPMs in terms of gain, dark count rate, and crosstalk probability. We found that the dark count rate of the MicroFJ was 16% higher than the ASD. Also, the gain of the MicroFJ is 3.5 times higher than the ASD. Finally, the crosstalk probability of the ASD 1.96 times higher than the MicroFJ. Our findings are in good agreement with manufacturer reported values.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations for arbitrary and unconventional tensor network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have a single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve eficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

We propose the use of balanced iterative reducing and clustering using hierarchies (BIRCH) combined with linear regression to predict the reduced Young's modulus and hardness of highly heterogeneous materials from a set of nanoindentation experiments. We first use BIRCH to cluster the dataset according to its mineral compositions, which are derived from the spectral matching of energy-dispersive spectroscopy data through the modular automated processing system (MAPS) platform. We observe that grouping our dataset into five clusters yields the best accuracy as well as a reasonable representation of mineralogy in each cluster. Subsequently, we test four types of regression models, namely linear regression, support vector regression, Gaussian process regression, and extreme gradient boosting regression. The linear regression and Gaussian process regression provide the most accurate prediction, and the proposed framework yields R2 = 0.93 for the test set. Although the study is needed more comprehensively, our results shows that machine learning methods such as linear regression or Gaussian process regression can be used to accurately estimate mechanical properties with a proper number of grouping based on compositional data.

Dynamical systems containing contact/impact between parts can be modeled as piecewise-smooth reduced-order models. The most common example is freeplay, which can manifest as a loose support, worn hinges, or backlash. Freeplay causes very complex, nonlinear responses in a system that range from isolated resonances to grazing bifurcations to chaos. This can be an issue because classical solution methods, such as direct time integration (e.g., Runge-Kutta) or harmonic balance methods, can fail to accurately detect some of the nonlinear behavior or fail to run altogether. To deal with this limitation, researchers often approximate piecewise freeplay terms in the equations of motion using continuous, fully smooth functions. While this strategy can be convenient, it may not always be appropriate for use. For example, past investigation on freeplay in an aeroelastic control surface showed that, compared to the exact piecewise representation, some approximations are not as effective at capturing freeplay behavior as other ones. Another potential issue is the effectiveness of continuous representations at capturing grazing contacts and grazing-type bifurcations. These can cause the system to transition to high-amplitude responses with frequent contact/impact and be particularly damaging. In this work, a bifurcation study is performed on a model of a forced Duffing oscillator with freeplay nonlinearity. Various representations are used to approximate the freeplay including polynomial, absolute value, and hyperbolic tangent representations. Bifurcation analysis results for each type are compared to results using the exact piecewise-smooth representation computed using MATLAB® Event Location. The effectiveness of each representation is compared and ranked in terms of numerical accuracy, ability to capture multiple response types, ability to predict chaos, and computation time.

Sandia National Laboratories has developed technology enabling novel downhole electrochemical assessment in extreme downhole environments. High-temperature high-pressure (HTHP) electrodes selectively sensitive to hydrogen (H+), chloride (Cl-), iodide (I-) and overall ionic strength (Reference Electrode+-) have been demonstrated in representative geothermal environments (225°C and 103 bar in surrogate geothermal brine). This 2-year program is a collaboration effort between Sandia and Thermochem, Inc. with the goal of taking the prototype sensors and developing them into a commercial product that is operable up to 300°C and 345 bar. The Sandia-developed prototype HTHP chemical sensor package creates a capability that has never been possible to date. This technology is desired by the geothermal industry to fill a gap in available downhole real-time measurements. Only limited sensors are available that operate at the extreme temperatures and pressures found in geothermal wells. For the purpose of this paper, high temperature is defined as temperatures exceeding 200°C and high pressure is defined as pressures exceeding 35 bar. Chemical sensors exceeding these parameters and sized appropriately for downhole applications do not exist. The current Thermochem two-phase downhole sampling tool (rated to 350 °C) will be re-configured to accept the sensors. A downhole tool with an integrated pH real-time sensor capable of operation at 300°C and 345 bar does not exist and as such, the developed technology will provide the geothermal industry with data that would otherwise not be possible such as vertical in-situ pH-profiling of geothermal wells. The pH measurement was chosen as the first chemical sensor focus since it is one of the fundamental measurements required to understand downhole chemistry, scaling and corrosion processes.

Abstract not provided.

Magnetic skyrmions are nanoscale whirls of magnetism that can be propagated with electrical currents. The repulsion between skyrmions inspires their use for reversible computing based on the elastic billiard ball collisions proposed for conservative logic in 1982. In this letter, we evaluate the logical and physical reversibility of this skyrmion logic paradigm, as well as the limitations that must be addressed before dissipation-free computation can be realized.