Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Nuclear power plant (NPP) risk assessment is broadly separated into disciplines of nuclear safety, security, and safeguards. Different analysis methods and computer models have been constructed to analyze each of these as separate disciplines. However, due to the complexity of NPP systems, there are risks that can span all these disciplines and require consideration of safety-security (2S) interactions which allows a more complete understanding of the relationship among these risks. A novel leading simulator/trailing simulator (LS/TS) method is introduced to integrate multiple generic safety and security computer models into a single, holistic 2S analysis. A case study is performed using this novel method to determine its effectiveness. The case study shows that the LS/TS method avoided introducing errors in simulation, compared to the same scenario performed without the LS/TS method. A second case study is then used to illustrate an integrated 2S analysis which shows that different levels of damage to vital equipment from sabotage at a NPP can affect accident evolution by several hours.

Abstract not provided.

Abstract not provided.

The additive manufacture of compositionally graded Al/Cu parts by laser engineered net shaping (LENS) is demonstrated. The use of a blue light build laser enabled deposition on a Cu substrate. The thermal gradient and rapid solidification inherent to selective laser melting enabled mass transport of Cu up to 4 mm from a Cu substrate through a pure Al deposition, providing a means of producing gradients with finer step sizes than the printed layer thicknesses. Divorcing gradient continuity from layer or particle size makes LENS a potentially enabling technology for the manufacture of graded density impactors for ramp compression experiments. Printing graded structures with pure Al, however, was prevented by the growth of Al2Cu3 dendrites and acicular grains amid a matrix of Al2Cu. A combination of adding TiB2 grain refining powder and actively varying print layer composition suppressed the dendritic growth mode and produced an equiaxed microstructure in a compositionally graded part. Material phase was characterized for crystal structure and nanoindentation hardness to enable a discussion of phase evolution in the rapidly solidifying melt pool of a LENS print.

The formation of magnesium chloride-hydroxide salts (magnesium hydroxychlorides) has implications for many geochemical processes and technical applications. For this reason, a thermodynamic database for evaluating the Mg(OH)2–MgCl2–H2O ternary system from 0 °C–120 °C has been developed based on extensive experimental solubility data. Internally consistent sets of standard thermodynamic parameters (ΔGf°, ΔHf°, S°, and CP) were derived for several solid phases: 3 Mg(OH)2:MgCl2:8H2O, 9 Mg(OH)2:MgCl2:4H2O, 2 Mg(OH)2:MgCl2:4H2O, 2 Mg(OH)2:MgCl2: 2H2O(s), brucite (Mg(OH)2), bischofite (MgCl2:6H2O), and MgCl2:4H2O. First, estimated values for the thermodynamic parameters were derived using a component addition method. These parameters were combined with standard thermodynamic data for Mg2+(aq) consistent with CODATA (Cox et al., 1989) to generate temperature-dependent Gibbs energies for the dissolution reactions of the solid phases. These data, in combination with values for MgOH+(aq) updated to be consistent with Mg2+-CODATA, were used to compute equilibrium constants and incorporated into a Pitzer thermodynamic database for concentrated electrolyte solutions. Phase solubility diagrams were constructed as a function of temperature and magnesium chloride concentration for comparisons with available experimental data. To improve the fits to the experimental data, reaction equilibrium constants for the Mg-bearing mineral phases, the binary Pitzer parameters for the MgOH+ — Cl− interaction, and the temperature-dependent coefficients for those Pitzer parameters were constrained by experimental phase boundaries and to match phase solubilities. These parameter adjustments resulted in an updated set of standard thermodynamic data and associated temperature-dependent functions. The resulting database has direct applications to investigations of magnesia cement formation and leaching, chemical barrier interactions related to disposition of heat-generating nuclear waste, and evaluation of magnesium-rich salt and brine stabilities at elevated temperatures.

Poster for CCN

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The 2022 National Defense Strategy of the United States listed climate change as a serious threat to national security. Climate intervention methods, such as stratospheric aerosol injection, have been proposed as mitigation strategies, but the downstream effects of such actions on a complex climate system are not well understood. The development of algorithmic techniques for quantifying relationships between source and impact variables related to a climate event (i.e., a climate pathway) would help inform policy decisions. Data-driven deep learning models have become powerful tools for modeling highly nonlinear relationships and may provide a route to characterize climate variable relationships. In this paper, we explore the use of an echo state network (ESN) for characterizing climate pathways. ESNs are a computationally efficient neural network variation designed for temporal data, and recent work proposes ESNs as a useful tool for forecasting spatiotemporal climate data. However, ESNs are noninterpretable black-box models along with other neural networks. The lack of model transparency poses a hurdle for understanding variable relationships. We address this issue by developing feature importance methods for ESNs in the context of spatiotemporal data to quantify variable relationships captured by the model. We conduct a simulation study to assess and compare the feature importance techniques, and we demonstrate the approach on reanalysis climate data. In the climate application, we consider a time period that includes the 1991 volcanic eruption of Mount Pinatubo. This event was a significant stratospheric aerosol injection, which acts as a proxy for an anthropogenic stratospheric aerosol injection. We are able to use the proposed approach to characterize relationships between pathway variables associated with this event that agree with relationships previously identified by climate scientists.

Abstract not provided.

Interpenetrating lattices consist of two or more interwoven but physically separate sub-lattices with unique behaviors derived from their multi-body construction. If the sublattices are constructed or coated with an electrically conducting material, the close proximity and high surface area of the electrically isolated conductors allow the two lattices to interact electromagnetically either across the initial dielectric filled gap or through physical contact. Changes in the size of the dielectric gap between the sub-lattices induced by deformation can be measured via capacitance or resistance, allowing a structurally competent lattice to operate as a force or deformation sensor. In addition to resistive and capacitive deformation sensing, this work explores capacitance as a fundamental metamaterial property and the environmental sensing behaviors of interpenetrating lattices.

Study of subcooled pool boiling experiments performed using a dielectric coolant to test effects of variations in heater surface configuration on pool boiling characteristics.

The rise of grid modernization has been prompted by the escalating demand for power, the deteriorating state of infrastructure, and the growing concern regarding the reliability of electric utilities. The smart grid encompasses recent advancements in electronics, technology, telecommunications, and computer capabilities. Smart grid telecommunication frameworks provide bidirectional communication to facilitate grid operations. Software-defined networking (SDN) is a proposed approach for monitoring and regulating telecommunication networks, which allows for enhanced visibility, control, and security in smart grid systems. Nevertheless, the integration of telecommunications infrastructure exposes smart grid networks to potential cyberattacks. Unauthorized individuals may exploit unauthorized access to intercept communications, introduce fabricated data into system measurements, overwhelm communication channels with false data packets, or attack centralized controllers to disable network control. An ongoing, thorough examination of cyber attacks and protection strategies for smart grid networks is essential due to the ever-changing nature of these threats. Previous surveys on smart grid security lack modern methodologies and, to the best of our knowledge, most, if not all, focus on only one sort of attack or protection. This survey examines the most recent security techniques, simultaneous multi-pronged cyber attacks, and defense utilities in order to address the challenges of future SDN smart grid research. The objective is to identify future research requirements, describe the existing security challenges, and highlight emerging threats and their potential impact on the deployment of software-defined smart grid (SD-SG).

Seismic waveform data recorded at stations can be thought of as a superposition of the signal from a source of interest and noise from other sources. Frequency-based filtering methods for waveform denoising do not result in desired outcomes when the targeted signal and noise occupy similar frequency bands. Recently, denoising techniques based on deep-learning convolutional neural networks (CNNs), in which a recorded waveform is decomposed into signal and noise components, have led to improved results. These CNN methods, which use short-time Fourier transform representations of the time series, provide signal and noise masks for the input waveform. These masks are used to create denoised signal and designaled noise waveforms, respectively. However, advancements in the field of image denoising have shown the benefits of incorporating discrete wavelet transforms (DWTs) into CNN architectures to create multilevel wavelet CNN (MWCNN) models. The MWCNN model preserves the details of the input due to the good time–frequency localization of the DWT. Here, we use a data set of over 382,000 constructed seismograms recorded by the University of Utah Seismograph Stations network to compare the performance of CNN and MWCNN-based denoising models. Evaluation of both models on constructed test data shows that the MWCNN model outperforms the CNN model in the ability to recover the ground-truth signal component in terms of both waveform similarity and preservation of amplitude information. Model evaluation of real-world data shows that both the CNN and MWCNN models outperform standard band-pass filtering (BPF; average improvement in signal-to-noise ratio of 9.6 and 19.7 dB, respectively, with respect to BPF). Evaluation of continuous data suggests the MWCNN denoiser can improve both signal detection capabilities and phase arrival time estimates.

Explosion sources have been observed to generate significant shear-wave energy despite their isotropic nature. To investigate this phenomenon, we conduct an analysis of the seismic data collected as part of the Source Physics Experiment (SPE): Dry Alluvium Geology (DAG) and investigate the generation of shear-wave energy via scattering. The data were produced by three underground chemical explosions and consist of three-component seis-mograms, which were recorded by the DAG Large-N array. Synthetic tests suggest that for the DAG experiments, small-scale stochastic heterogeneities, defined as features with correlation lengths of 10–100s of meters, are more effective than large-scale geologic structure (scales >1–10 km) at reproducing the scattering of explosion generated wavefields observed at DAG. We analyze the seismic data for spatially variable ratios between transversely and radially polarized seismic energy, and then estimate the mean free path of P and S waves. All analyses are conducted within a frequency band of 5–50 Hz. The ratio of transversely to radially polarized energy is the highest in the east and west portion of the Large-N array. In addition, the magnitude of the estimated S-wave mean free path is shorter in the eastern portion of the Large-N array. This variation indicates that the eastern area of the DAG array is where more scattering is occurring, suggesting azimuthal dependence of P-to-P and P-to-S scattering. This azimuthal dependence of P-to-S scattering can have implications for explosion discrimination based on spectral ratios of seismic wave types, because the general assumption is that explosions do not generate shear-wave energy. Synthetic tests modeling only larger-scale geologic structure had lower transversely polarized energy (only four stations showing a transversely to radially polarized energy ratio greater than 1) and fewer stations (<10) displaying shorter (<300 m) mean free paths than what was observed in the DAG data results.

Abstract not provided.

This study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3 % to 25 %. Differences in the excitation forces and torques, however, are contained within 12 %. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The rotating bending fatigue (RBF) behavior (fully reversed, R = −1) of additively manufactured (AM) Ti–6Al–4V alloy produced via laser powder bed fusion (PBF-L) was investigated with respect to different microstructures achieved through novel heat treatments. The investigation herein seeks to elucidate the effect of microstructure by controlling variables that can affect fatigue behavior in Ti–6Al–4V, such as chemistry, porosity, and surface roughness. In order to control these variables, different hot isostatic pressing (HIP) treatments at 800 °C, 920 °C, and 1050 °C with a 920 °C temper were applied to three sets of Ti–6Al–4V cylinders that originated from the same PBF-L build, such that there were 30 tests per condition. After HIP treatment, the specimens were machined and tested. The highest runout stress was achieved after sub-β transus HIP at 800 °C for 2 h at 200 MPa of pressure. A significant drop in fatigue strength was attributed to large prior-β grains and grain boundary α resulting from super-β transus HIP treated specimens. For the sub-β transus HIP specimens, differences in fatigue strength were attributed to α lath thickness, relative dislocation density, and dislocation boundary strengthening.

Naturally occurring uranium complicates monitoring for occupational exposures. There are several retroactive methods that can be used to monitor for occupational exposures, with benefits and drawbacks to each. Analysis of uranium in urine by mass spectrometry and alpha spectrometry is compared, and methods of determining an occupational exposure are presented. Furthermore, the minimum detectable concentrations from each analysis and a method for intake determination based on the analytical results are compared for various solubility types and mixtures. Mass spectrometry with radiochemical separation was found to be the most sensitive analysis for detecting occupational exposures to anthropogenic mixtures based on minimum detectable doses calculated from the proposed method for intake determination.

The critical stress for cutting of a void and He bubble (generically referred to as a cavity) by edge and screw dislocations has been determined for FCC Fe0.70Cr0.20Ni0.10—close to 300-series stainless steel—over a range of cavity spacings, diameters, pressures, and glide plane positions. The results exhibit anomalous trends with spacing, diameter, and pressure when compared with classical theories for obstacle hardening. These anomalies are attributed to elastic anisotropy and the wide extended dislocation core in low stacking fault energy metals, indicating that caution must be exercised when using perfect dislocations in isotropic solids to study void and bubble hardening. In many simulations with screw dislocations, cross-slip was observed at the void/bubble surface, leading to an additional contribution to strengthening. We refer to this phenomenon as cavity cross-slip locking, and argue that it may be an important contributor to void and bubble hardening.

Abstract not provided.

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. We present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

Abstract not provided.

Abstract not provided.

Radiation and radioactive substances result in the production of radioactive wastes which require safe management and disposal to avoid risks to human health and the environment. To ensure permanent safe disposal, the performance of a deep geological repository for radioactive waste is assessed against internationally agreed risk-based standards. Assessing postclosure safety of the future system's evolution includes screening of features, events, and processes (FEPs) relevant to the situation, their subsequent development into scenarios, and finally the development and execution of safety assessment (SA) models. Global FEP catalogs describe important natural and man-made repository system features and identify events and processes that may affect these features into the future. By combining FEPs, many of which are uncertain, different possible future system evolution scenarios are derived. Repository licensing should consider both the reference or “base” evolution as well as alternative futures that may lead to radiation release, pollution, or exposures. Scenarios are used to derive and consider both base and alternative evolutions, often through production of scenario-specific SA models and the recombination of their results into an assessment of the risk of harm. While the FEP-based scenario development process outlined here has evolved somewhat since its development in the 1980s, the fundamental ideas remain unchanged. A spectrum of common approaches is given here (e.g., bottom–up vs. top–down scenario development, probabilistic vs. bounding handling of uncertainty), related to how individual numerical models for possible futures are converted into a determination as to whether the system is safe (i.e., how aleatoric uncertainty and scenarios are integrated through bounding or Monte Carlo approaches).

Abstract not provided.

Abstract not provided.

Hydrogen powered locomotives are being explored to reduce emissions in rail applications. The risks of operations like refueling should be understood to ensure safe environments for workers and members of the public. Sensitivity analyses were conducted using HyRAM+ to identify major drivers of risk and compare effects of system parameters on individual risk. The consequences of jet fires from full-bore leaks dominated the risk, compared to explosions or smaller leaks. Pipe size, leak detection capability, and leak frequencies of system components greatly affected risk while overpressure modeling parameters and ambient conditions had little effect. The effects of personal protective equipment (PPE) materials on individual risk were quantified by reducing the individual’s exposure time or absorbed thermal dose. PPE only showed a risk reduction in low-risk cases. This study highlighted target areas for risk mitigation, including leak detection equipment and component maintenance, and indicated that the minimal effects of other parameters on risk may not justify prescriptive requirements for refueling operators.

Typical QRAs provide deterministic estimates and understanding of risks posed but are constructed using significant assumptions and uncertainties due to limited data availability and historical momentum of using nominal estimates. This report presents a hydrogen QRA analysis using HyRAM+ that incorporates uncertainty with Latin hypercube sampling and sensitivity analysis using linear regression.

Abstract not provided.

Abstract not provided.

This study conducts a comparative analysis, using non-equilibrium Green’s functions (NEGF), of two state-of-the-art two-well (TW) Terahertz Quantum Cascade Lasers (THz QCLs) supporting clean 3-level systems. The devices have nearly identical parameters and the NEGF calculations with an abrupt-interface roughness height of 0.12 nm predict a maximum operating temperature (T_max) of ~ 250 K for both devices. However, experimentally, one device reaches a T_max of ~ 250 K and the other a T_max of only ~ 134 K. Both devices were fabricated and measured under identical conditions in the same laboratory, with high quality processes as verified by reference devices. The main difference between the two devices is that they were grown in different MBE reactors. Our NEGF-based analysis considered all parameters related to MBE growth, including the maximum estimated variation in aluminum content, growth rate, doping density, background doping, and abrupt-interface roughness height. From our NEGF calculations it is evident that the sole parameter to which a drastic drop in T_max could be attributed is the abrupt-interface roughness height. We can also learn from the simulations that both devices exhibit high-quality interfaces, with one having an abrupt-interface roughness height of approximately an atomic layer and the other approximately a monolayer. However, these small differences in interface sharpness are the cause of the large performance discrepancy. This underscores the sensitivity of device performance to interface roughness and emphasizes its strategic role in achieving higher operating temperatures for THz QCLs. We suggest Atom Probe Tomography (APT) as a path to analyze and measure the (graded)-interfaces roughness (IFR) parameters for THz QCLs, and subsequently as a design tool for higher performance THz QCLs, as was done for mid-IR QCLs. Our study not only addresses challenges faced by other groups in reproducing the record T_max of ~ 250 K and ~ 261 K but also proposes a systematic pathway for further improving the temperature performance of THz QCLs beyond the state-of-the-art.

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

Abstract not provided.

Abstract not provided.

We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

This study investigates high performance electrochromic windows used on a passive house and residential dwelling to IECC 2021 (i.e., IECC dwelling). In the lab, the electrochromic film switches transmitted solar heat gain coefficient (SHGC) from 0.09 to 0.7 and visible transmittance from 0.15 to 0.82 with power consumption of 1.23 W/m2 during switching times less than 3 minutes. We extrapolate these results to a window assembly. Building energy models of the houses were evaluated in Santa Fe, New Mexico. A Monte Carlo analysis for 2020, 2040, 2060, and 2080 was conducted for Shared Socioeconomic Pathways 2-4.5, 3-7.0, and 5-8.5. Cases with and without the electrochromic windows and with and without electricity were used to determine energy use intensity and hours beyond thermal safety thresholds. The passive house showed 1.3-3.1% mean energy savings and the IECC dwelling 4.4-5.1% with electrochromic efficiency benefits growing into the future for both cases. Even so, overall savings decrease into the future for the passive house, due to growth in cooling load being dominant, conversely overall energy savings increase into the future for the IECC dwelling due to heating loads being dominant. For thermal resilience, the passive house exhibited a mean percent decrease of 0.02-0.31% hours in the extreme caution (i.e., > 32.2 ∘C, ≤ 39.4 ∘C) range while the IECC dwelling exhibited 0.38-4.38%. The study therefore shows that electrochromic windows will have smaller benefits for the passive house in comparison to the IECC dwelling. The relationship between electrochromic windows is shown to have a complex relationship between house efficiency and climate change by these results.