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.
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.
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.
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.
