This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.
DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75th percentile for all benchmarks. At the 90th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.
The association of ionizable polymers strongly affects their motion in solutions, where the constraints arising from clustering of the ionizable groups alter the macroscopic dynamics. The interrelation between the motion on multiple length and time scales is fundamental to a broad range of complex fluids including physical networks, gels, and polymer-nanoparticle complexes where long-lived associations control their structure and dynamics. Using neutron spin echo and fully atomistic, multimillion atom molecular dynamics (MD) simulations carried out to times comparable to that of chain segmental motion, the current study resolves the dynamics of networks formed by suflonated polystryene solutions for sulfonation fractions 0 ≤ f ≤ 0.09 across time and length scales. The experimental dynamic structure factors were measured and compared with computational ones, calculated from MD simulations, and analyzed in terms of a sum of two exponential functions, providing two distinctive time scales. These time constants capture confined motion of the network and fast dynamics of the highly solvated segments. A unique relationship between the polymer dynamics and the size and distribution of the ionic clusters was established and correlated with the number of polymer chains that participate in each cluster. The correlation of dynamics in associative complex fluids across time and length scales, enabled by combining the understanding attained from reciprocal space through neutron spin echo and real space, through large scale MD studies, addresses a fundamental long-standing challenge that underline the behavior of soft materials and affect their potential uses.
Methanol emerges as a compelling renewable fuel for decarbonizing engine applications due to a mature industry with high production capacity, existing distribution infrastructure, low carbon intensity and favorable cost. Methanol's high flame speed and high autoignition resistance render it particularly well-suited for spark-ignition (SI) engines. Previous research showed a distinct phenomenon, known deflagration-based knock in methanol combustion, whereby knocking combustion was observed albeit without end-gas autoignition. This work studies the implications of deflagration-based knock on noise emissions by investigating the knock intensity and combustion noise at knock-limited operation of methanol in a single-cylinder direct-injection SI engine operated at both stoichiometric and lean (λ = 2.0) conditions. Results are compared against observations from a premium-grade gasoline. Experiments show that methanol's end-gas autoignition occurs at lean conditions, leading to the typical autoignition-based knock as that occurring with premium-grade gasoline. However, at stoichiometric conditions, knock-limited operation is achieved with deflagration-based knock. Noise of deflagration-based knock has lower variability than that of autoignition-based knock and it does not seem to be an issue at the engine speed tested experimentally in this paper (1400 rpm). However, computational fluid dynamic large eddy simulations show that deflagration-based knock may lead to high noise levels at 2000 rpm. Deflagration-based knock is insensitive to changing spark timings, so new knock mitigation strategies are required, such as adjusting the spark energy and/or adding dilution. Finally, this study shows that deflagration-based-knock may be directly impacted by the flame speed, occurring more frequently with faster-burning fuels or under conditions that elevate flame speeds, like rich-stoichiometric operation. The finding bears implications on renewable e-fuels, such as ethanol, methanol and hydrogen.
The ability to visualize x-ray and neutron emission from fusion plasmas in 3D is critical to understand the origin of the complex shapes of the plasmas in experiments. Unfortunately, this remains challenging in experiments that study a fusion concept known as Magnetized Liner Inertial Fusion (MagLIF) due to a small number of available diagnostic views. Here, we present a basis function-expansion approach to reconstruct MagLIF stagnation plasmas from a sparse set of x-ray emission images. A set of natural basis functions is “learned” from training volumes containing quasi-helical structures whose projections are qualitatively similar to those observed in experimental images. Tests on several known volumes demonstrate that the learned basis outperforms both a cylindrical harmonic basis and a simple voxel basis with additional regularization, according to several metrics. Two-view reconstructions with the learned basis can estimate emission volumes to within 11% and those with three views recover morphology to a high degree of accuracy. The technique is applied to experimental data, producing the first 3D reconstruction of a MagLIF stagnation column from multiple views, providing additional indications of liner instabilities imprinting onto the emitting plasma.
In this article we present a quantitative analysis of the second positive system of molecular nitrogen and the first negative system of the molecular nitrogen cation excited in the presence of ionizing radiation. Optical emission spectra of atmospheric air and nitrogen surrounding 210Po sources were measured from 250 to 400 nm. Multi-Boltzmann and non-Boltzmann vibrational distribution spectral models were used to determine the vibrational temperature and vibrational distribution function of the emitting N2(C3Πu) and N2+(B2Σ+u) states. A zero-dimensional kinetic model, based on the electron energy distribution function (EEDF) and steady-state excitation and de-excitation of N2(X1Σ+g), N2+(B2Σ+u), N2+(X2Σ+g), N4+, O2+, and N2(C3Πu, v), was developed for the prediction of the relative spectral intensity of both the N2+(B2Σ+u → X2Σ+g) emission band and the vibrational bands of N2(C3Πu → B3Πg) for comparison with the experimental data.
Shahili, Mohammad; Addamane, Sadhvikas J.; Kim, Anthony D.; Curwen, Christopher A.; Kawamura, Jonathan H.; Williams, Benjamin S.
Design strategies for improving terahertz (THz) quantum cascade lasers (QCLs) in the 5-6THz range are investigated numerically and experimentally, with the goal of overcoming the degradation in performance that occurs as the laser frequency approaches the Reststrahlen band. Two designs aimed at 5.4THz were selected: one optimized for lower power dissipation and one optimized for better temperature performance. The active regions exhibited broadband gain, with the strongest modes lasing in the 5.3-5.6THz range, but with other various modes observed ranging from 4.76 to 6.03THz. Pulsed and continuous-wave (cw) operation is observed up to temperatures of 117K and 68K, respectively. In cw mode, the ridge laser has modes up to 5.71THz - the highest reported frequency for a THz QCL in cw mode. The waveguide loss associated with the doped contact layers and metallization is identified as a critical limitation to performance above 5THz.
ORIGEN is one of the main transmutation software packages used in nuclear engineering Modeling Tungsten Boride Neutronics in ORIGEN for Z-Facilityproblems. For the case of this study, tungsten borides are studied using a coupled framework between MCNP and the ORIGEN package of scale. The input used four compositions of tungsten boride: WB with natural boron- 10 abundance, WB with 80wt% B-10 per isotope of boron, WB4 with natural boron-10 abundance, and WB4 with 80wt% B-10 per isotope of boron. Isotopic inventories were produced for WB which show the time dependent change up to 2 years after a 6-Month irradiation. This will allow for further studies of the materials to assess things material composition changes, dose contribution, and waste management requirements.
Understanding temperature-dependent material decomposition and structural deformation induced by combined thermal-mechanical environments is critical for safety qualification of hardware under accident scenarios. Seeing in with X-rays elucidated the physics necessary to develop X-ray strain and thermometry diagnostics for use in optically opaque environments. Two parallel thermometry schemes were explored: X-ray fluorescence and X-ray diffraction of inorganic doped ceramics– colloquially known as thermographic phosphors. Two parallel surface strain techniques–Path-Integrated Digital Image Correlation and Frequency Multiplexed Digital Image Correlation–were demonstrated. Finally, preliminary demonstration of time-resolved digital volume correlation was performed by taking advantage of limited view reconstruction techniques. Additionally, research into blended ceramic-metal coatings was critical to generating intrinsic thermographic patterns for the future combination of X-ray strain and thermometry measurements.
The properties of defects in n-p-n Si bipolar junction transistors (BJTs) caused by 17-MeV Si ions are investigated via current-voltage, low-frequency (LF) noise, and deep level transient spectroscopy (DLTS) measurements. Four prominent radiation-induced defects in the base-collector junction of these transistors are identified via DLTS. At least two defect levels are observed in temperature-dependent LF 1/f noise measurements, one that is similar to a prominent defect in DLTS and another that is not. Defect microstructures are discussed. Our results show that DLTS and 1/f noise measurements can provide complementary information about defects in linear bipolar devices.
Long-duration energy storage (LDES) is critical to a stable, resilient, and decarbonized electric grid. While batteries are emerging as important LDES devices, extended, high-power discharges necessary for cost-competitive LDES present new materials challenges. Focusing on a new generation of low-temperature molten sodium batteries, we explore here unique phenomena related to long-duration discharge through a well-known solid electrolyte, NaSICON. Specifically, molten sodium symmetric cells at 110 °C were cycled at 0.1 A cm−2 for 1-23 h discharges. Longer discharges led to unstable overpotentials, reduced resistances, and decreased electrolyte strength, caused by massive sodium penetration not observed in shorter duration discharges. Scanning electron microscopy informed mechanisms of sodium penetration and even “healing” during shorter-duration cycling. Importantly, these findings show that traditional, low-capacity, shorter-duration tests may not sufficiently inform fundamental materials phenomena that will impact LDES battery performance. This case highlights the importance that candidate LDES batteries be tested under pertinent long-duration conditions.
Sandia National Laboratories (SNL) requested a measure of possible human error for each state verification method for a safety mechanism performed at partnering production agencies. A team of three human factors individuals were tasked with conducting observations during site visits of both production agencies in order to complete a Human Reliability Analysis (HRA). A HRA will be used because it provides both qualitative and quantitative reports of human error. This report is the first phase of that effort, which will describe the methods which occur at one of the production agencies.
Abstract: The tetrahedron, as the simplest platonic shape, is a profound building block with the potential to create intricate superstructures. Noteworthy designs utilizing tetrahedral building blocks include the Sierpiński tetrahedron (the most fundamental three-dimensional fractal), a one-dimensional helical structure known as the tetrahelix, and various crystalline and quasicrystalline packings. Historically, the practicality of tetrahedral superstructures has been evident, providing stable, well-defined frameworks for various constructions, including truss bridges, tower cranes, and electricity transmission line pylons. In the field of self-assembled nanocrystal superlattices, tetrahedral nanocrystals, as building blocks, occupy a unique place among all the possible nanoscale particles. Mathematical models, simulation work, and experimental studies using nanocrystals in the laboratory have suggested that self-assembled structures derived from nanoscale tetrahedral building blocks are notably intricate, giving rise to new horizons of high-entropy nanocrystal superlattices. An important implication from previous works is that such tetrahedral nanocrystal superlattices form through highly delicate interparticle interactions, emphasizing the importance of the fine features of these nanocrystals. In this article, we summarize the advances in superlattices assembled from tetrahedral nanocrystals. We first define the tetrahedron and tetrahedron analogues based on Conway’s transformation and graph theory, underscoring their relevance to the crystallization process producing tetrahedral nanocrystals. Then, we showcase previous reports on the synthesis of tetrahedral nanocrystals and the resulting nanocrystal superstructures. Finally, we conclude by offering insights and perspective into the chemical and architectural intricacy that could emerge from tetrahedral nanocrystals. Graphical abstract: (Figure presented.).
Johnson, Dylan M.; Khakhum, Nittaya; Wang, Min; Warner, Nikole L.; Jokinen, Jenny D.; Comer, Jason E.; Lukashevich, Igor S.
Lymphocytic choriomeningitis virus (LCMV) and Lassa virus (LASV) share many genetic and biological features including subtle differences between pathogenic and apathogenic strains. Despite remarkable genetic similarity, the viscerotropic WE strain of LCMV causes a fatal LASV fever-like hepatitis in non-human primates (NHPs) while the mouse-adapted Armstrong (ARM) strain of LCMV is deeply attenuated in NHPs and can vaccinate against LCMV-WE challenge. Here, we demonstrate that internalization of WE is more sensitive to the depletion of membrane cholesterol than ARM infection while ARM infection is more reliant on endosomal acidification. LCMV-ARM induces robust NF-κB and interferon response factor (IRF) activation while LCMV-WE seems to avoid early innate sensing and failed to induce strong NF-κB and IRF responses in dual-reporter monocyte and epithelial cells. Toll-like receptor 2 (TLR-2) signaling appears to play a critical role in NF-κB activation and the silencing of TLR-2 shuts down IL-6 production in ARM but not in WE-infected cells. Pathogenic LCMV-WE infection is poorly recognized in early endosomes and failed to induce TLR-2/Mal-dependent pro-inflammatory cytokines. Following infection, Interleukin-1 receptor-associated kinase 1 (IRAK-1) expression is diminished in LCMV-ARM- but not LCMV-WE-infected cells, which indicates it is likely involved in the LCMV-ARM NF-κB activation. By confocal microscopy, ARM and WE strains have similar intracellular trafficking although LCMV-ARM infection appears to coincide with greater co-localization of early endosome marker EEA1 with TLR-2. Both strains co-localize with Rab-7, a late endosome marker, but the interaction with LCMV-WE seems to be more prolonged. These findings suggest that LCMV-ARM’s intracellular trafficking pathway may facilitate interaction with innate immune sensors, which promotes the induction of effective innate and adaptive immune responses.
