Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Limit Branch”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

A tale of a tail: a tidally disrupting ultra-diffuse galaxy in the M81 group

ABSTRACT We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper Suprime-Cam survey of the M81 Group. It has an average surface brightness of μg ∼ 32 mag arcsec−2 and can be traced for over a degree on the sky (60 kpc at the distance of F8D1) with our current imagery. We revisit the main body properties of F8D1 using deep multiband imagery acquired with MegaCam on CFHT and measure effective radii of 1.7–1.9 kpc, central surface brightnesses of 24.7–25.7 mag, and a stellar mass of ∼7 × 107M⊙. Assuming a symmetric feature on the other side of the galaxy, we calculate that 30–36 per cent of F8D1’s present-day luminosity is contained in the tail. We argue that the most likely origin of F8D1’s disruption is a recent close passage to M81, which would have stripped its gas and quenched its star formation. As the only UDG that has so far been studied to such faint surface brightness depths, the unveiling of F8D1’s tidal disruption is important. It leaves open the possibility that many other UDGs could be the result of similar processes, with the most telling signatures of this lurking below current detection limits.

Žemaitis, Rokas (ORCID:0000000285660491)↗

Upgrades to the Thomson Scattering Diagnostic on the Lithium Tokamak eXperiment-Beta (LTX-$β$)

We describe the upgrades to the Thomson scattering (TS) diagnostic on the Lithium Tokamak eXperiment-beta (LTX- β ) to provide complete radial coverage of the plasma. LTX- β is a spherical tokamak with plasma limited at a major radius of R = 14 cm and R = 66 cm. Flat temperature profiles with peaked density were observed in LTX- β using an existing multipoint TS diagnostic with radial coverage of the plasma core. There is a need for more comprehensive core and edge TS measurements to thoroughly examine beam fueling, beam heating, and transport under low recycling conditions with liquid lithium walls and a lower collisionality scrape-off layer (SOL). The core TS set up provides plasma profiles from R = 40–60 cm. To extend the TS diagnostic coverage: 1) we added a TS viewing optics to cover the SOL region at five radial points from R = 57–64 cm, and 2) another TS viewing optics in the HFS region from R = 15–45 cm to obtain TS data from ten radial positions while making use of the original laser beam. The main focus of this article is on the new HFS TS setup that is equipped with better-matched collection optics located in the horizontal midplane where the laser polarization lies and views from a scattering angle of 130°–150° with respect to the laser path. A set of unequal-length Y-branched fiber bundles transports the TS light from two radial points to each polychromator. The new HFS view is better matched in spot size, scattering angle, and throughput to the laser and polychromators and has new V-grooved style viewing dumps. The new HFS TS view, along with the core and SOL TS views, will provide a full radial measurement coverage from R = 15–64 cm to explore the flat temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast GPU-Based Generation of Large Graph Networks From Degree Distributions

Synthetically generated, large graph networks serve as useful proxies to real-world networks for many graph-based applications. The ability to generate such networks helps overcome several limitations of real-world networks regarding their number, availability, and access. Here, we present the design, implementation, and performance study of a novel network generator that can produce very large graph networks conforming to any desired degree distribution. The generator is designed and implemented for efficient execution on modern graphics processing units (GPUs). Given an array of desired vertex degrees and number of vertices for each desired degree, our algorithm generates the edges of a random graph that satisfies the input degree distribution. Multiple runtime variants are implemented and tested: 1) a uniform static work assignment using a fixed thread launch scheme, 2) a load-balanced static work assignment also with fixed thread launch but with cost-aware task-to-thread mapping, and 3) a dynamic scheme with multiple GPU kernels asynchronously launched from the CPU. The generation is tested on a range of popular networks such as Twitter and Facebook, representing different scales and skews in degree distributions. Results show that, using our algorithm on a single modern GPU (NVIDIA Volta V100), it is possible to generate large-scale graph networks at rates exceeding 50 billion edges per second for a 69 billion-edge network. GPU profiling confirms high utilization and low branching divergence of our implementation from small to large network sizes. For networks with scattered distributions, we provide a coarsening method that further increases the GPU-based generation speed by up to a factor of 4 on tested input networks with over 45 billion edges.

97 MATHEMATICS AND COMPUTING↗

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementing belowground controls on nutrient uptake in ELMv2-SPRUCE improves representation of a boreal peatland ecosystem

Boreal peatlands store 13 %–32 % of the global soil carbon (C) stock, a service dependent on plant-mycorrhizal fungi associations. In these nutrient poor systems, ectomycorrhizal and ericoid mycorrhizal fungi supply up to >80 % of the nutrient requirements of their plant hosts, partly with mined nitrogen (N) and phosphorus (P) from soil organic matter that are otherwise inaccessible to plants. Despite the ecological significance, mycorrhizal associations are only represented in a few land surface or ecosystem models. We modify the peatland branch of version 2 of the Energy Exascale Earth System Land Model (ELMv2-SPRUCE) to replace the default photosynthesis-driven inorganic N and P (NP) uptake process with a more realistic representation of the process via three pathways: (1) direct inorganic NP uptake by uncolonized fine roots, (2) indirect inorganic NP acquisition and (3) indirect NP acquisition from organic sources by mycorrhizal roots. We systematically evaluated the performance of the default and modified models with field observations from a whole ecosystem warming and carbon dioxide fertilization experimental site: Spruce and Peatland Responses Under Changing Environment (SPRUCE), in northern Minnesota, USA. The modified model reduces the underestimation of the growth response of shrubs in the default model to warming from 40 %–80 % to 17 %–35 % and reduces the overall relative absolute error on C fluxes from 1.61 to 1.54 in calibration. Improvements on modeled shrub growths and shrub-moss community net ecosystem exchanges are also seen in validation. The improved growth response of shrubs to warming is accompanied by several-fold increase in direct inorganic NP uptake and decrease in fungal colonization rate. The modified model simulates a smaller magnitude of transition of the ecosystem from C sink to C source under warming due to alleviation of plant nutrient limitation. Equifinality analysis shows the newly added parameters in the modified model can be constrained by the observed C fluxes. Sensitivity analysis shows the newly added parameters have stronger statistical interactions than the preexisting parameters in the default model. Overall, the modified model is an improvement over the default ELMv2-SPRUCE and will be a useful tool for understanding boreal peatland change.

Wang, Yaoping [Oak Ridge National Laboratory (ORNL↗

Sequential polymer infusion into solid substrates (SPISS): Impact of processing on sorbent CO 2 adsorption properties

Solid sorbents made of small amine molecules and polyamines infused into mesoporous substrates are promising materials for CO 2 capture technologies. To date, their preparation is mainly based on wet infusion with focus on increasing amine content by varying the structure of the amine sorbent and designing solid substrates of various pore parameters. Less explored in the field are changes in processing to afford efficient CO 2 sorbents. In this study, branched poly(ethylenimine) (bPEI, M w = 800 Da) alone and blends with linear poly(propylenimine) (LPPI, M n = 6,700 Da) are infused into solid SBA-15 substrates by a method that varies the solution processing called sequential polymer infusion into solid substrates (SPISS). The reference 40 % bPEI-SBA-15 samples are split by two methods: split batch in dry suspension (SBD) and split batch in liquid suspension (SBL). Sequentially, alcoholic 10% of bPEI (non-blends) and 10% of LPPI (blends) solutions are introduced to afford the desired products. Under dry conditions, the resulting 50% bPEI-SBA-15 SBD & SBL sorbents display high CO 2 capacities up to 3.47 mmol CO 2 /g SiO 2 for simulated flue gas (10% CO 2 ) and 2.62 mmol CO 2 /g SiO 2 for direct air capture (DAC, 400 ppm CO 2 ). Under humid DAC conditions the CO 2 performance is further enhanced with an uptake of 4.62 mmol CO 2 /gSiO 2 and amine efficiency of 0.22 mmol CO 2 /mmol N. Subjected to extended temperature swing adsorption kinetic cycling (20 cycles), the SPISS samples display stable working capacities and retain over 70% (blends) and over 90% (non-blends) of their initial 12 h adsorption performance. LPPI is demonstrated to be an effective water sorption limiting agent using dynamic vapor sorption measurements. Solid state NMR techniques reveal important insights into the dynamics of the amine polymers confined into SBA-15 pores, as impacted by processing conditions. The results suggest that the conformation of the polymers is different depending on the processing method, displaying relatively tight (SBD) and loose (SBL) packing. Finally, the simple solution processing approaches presented here show that processing variations may guide the design of solid amine sorbents with desirable properties relevant for integration into CO 2 capture technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement placement in electric power transmission and distribution grids: Review of concepts, methods, and research needs

Sensing and measurement systems are quintessential to the safe and reliable operation of electric power grids. Their strategic placement is of ultimate importance because it is not economically viable to install measurement systems on every node and branch of a power grid, though they need to be monitored. An overwhelming number of strategies have been developed to meet oftentimes multiple conflicting objectives. The prime challenge in formulating the problem lies in developing a heuristic or an optimisation model that, though mathematically tractable and constrained in cost, leads to trustworthy technical solutions. Further, large-scale, long-term deployments pose additional challenges because the boundary conditions change as technologies evolve. For instance, the advent of new technologies in sensing and measurement, as well as in communications and networking, might impact the cost and performance of available solutions and shift initially set conditions. Also, the placement strategies developed for transmission grids might not be suitable for distribution grids, and vice versa, because of unique characteristics; therefore, the strategies need to be flexible, to a certain extent, because no two power grids are alike. Despite the extensive literature on the present topic, the focus of published works tends to be on a specific subject, such as the optimal placement of measurements to ensure observability in transmission grids. There is a dearth of work providing a comprehensive picture for developing optimal placement strategies. Because of the ongoing efforts on the modernisation of electric power grids, there is a need to consolidate the status quo while exposing its limitations to inform policymakers, industry stakeholders, and researchers on the research-and-development needs to push the boundaries for innovation. Accordingly, this paper first reviews the state-of-the-art considering both transmission and distribution grids. Then, it consolidates the key factors to be considered in the problem formulation. Finally, it provides a set of perspectives on the measurement placement problem, and it concludes with future research directions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Superallowed 0 + → 0 + nuclear β decays: 2020 critical survey, with implications for V ud and CKM unitarity

A new critical survey of all half-life, decay-energy and branching-ratio measurements related to 23 superallowed 0 + → 0 + β decays is presented. Included are 222 individual measurements of comparable precision obtained from 174 published references. Compared with our last survey in 2015, we have added results from 28 new publications and eliminated an approximately equal number whose results have been superseded by much more precise modern data. We obtain world-average ft values for each of the 21 transitions that have a complete set of data, then apply radiative and isospin-symmetry-breaking corrections to extract “corrected” Ft values. Fifteen of these Ft values now have a precision of 0.3% or better and all take the same value within statistics, as expected from conservation of the vector current. Their average, Ft¯, when combined with the muon lifetime, yields the up-down quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix, V ud = 0.97373 ± 0.00031. This is lower than our 2015 result by one standard deviation and its uncertainty is increased by 50%. This is a consequence, not of any shifts in the experimental data, but of new calculations for the radiative corrections. The lower V ud value now leads to greater tension in the top-row test of unitarity in the CKM matrix. Updates in experimental data have independently led to a factor-of-two tighter limit being set on the possible existence of a scalar interaction. In conclusion, the new limit on Fierz interference is b F ≤ 0.0033 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Wave energy converter buoy for Arctic observations

Ocean observation buoys are currently powered-constrained by battery storage capacity or available solar power. Power constraints limit the number of measurements that can be made and the lifetime of the buoys. Ocean surface wave energy could be used for power production, but wave energy converters are not yet commercially available for ocean observation buoys. Here we present the design and testing of a drifting wave energy converter buoy using a pendulum transmission system (PTS). The wave energy converter buoy was designed to operate in Arctic temperatures and wave conditions, but it could be used in areas with warmer temperatures and larger waves. Field tests measured a maximum power production of 5 watts in waves with a significant wave height of 0.42 meters and a ten minute average power production of 37 milliwatts. Power production increased with the energy period and significant wave height. The energy harvesting capabilities of the PTS showed the utility of adding a wave energy conversion device to a drifting ocean observation buoy.

Arctic observations↗

User Acceptance of Shared Autonomous Vehicles

A dissertation is proposed to explore user acceptance of shared autonomous vehicles (SAVs). SAVs are facing limited user acceptance. To systematically tackle the user acceptance barriers of SAVs, the main problem can be disintegrated into two sub-problems of user acceptance of autonomous vehicles (AVs) and ridesharing. The comfort of the ride experience in AVs is a determinant of user acceptance. Understanding the influential factors and developing methodologies to quantify human comfort in AVs are essential to facilitating future research to improve human comfort in AVs. The current pooled rideshare (PR) service closely resembles the anticipated future of SAVs. Understanding why users prefer or refuse to use PR at the current stage prepares SAVs for broader acceptance in the future. Until now, a series of peer-reviewed publications have been published to achieve the technical goals. Two simulator-based user studies were conducted to instrument the research on human comfort in AVs. Statistical analysis was performed to identify the crucial vehicular behavioral factors of human comfort in AVs. The influential factors of human comfort in AVs and methodologies to quantify and detect human comfort in AVs were investigated. Two survey-based studies were deployed to facilitate the investigation of user acceptance of rideshare services. The influential factors of users' willingness to consider PR were explored and identified, and the choice behaviors in ridesharing services were comprehensively modeled and analyzed. The proposed research answers a series of fundamental questions regarding the user acceptance of SAVs. For the branch of user acceptance of AVs, the research generated guidelines for improving passenger comfort in AVs by identifying a series of autonomous driving factors of passenger comfort. The research also provides fundamental tools to estimate human comfort levels for future research and in-AV applications. For the branch of user acceptance of PR, the research provided user acceptance-aware vehicle, service, and policy design insights that can promote the usage of PR.

Su, Haotian↗

From hydrogen to helium: the spectral evolution of white dwarfs as evidence for convective mixing

ABSTRACT We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the Gaia-DR2 catalogue, the Sloan Digital Sky Survey (SDSS), and the Galaxy Evolution Explorer (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20 000 to 9000 K, we find that 22 per cent of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5σ. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs – which include the full overshoot region – lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60 per cent of white dwarfs must have a hydrogen mass larger than MH/MWD = 10−10, another 25 per cent have masses in the range MH/MWD = 10−14–10−10, and 15 per cent have less hydrogen than MH/MWD = 10−14. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch and accretion of planetesimals on to white dwarfs.

Cunningham, Tim↗

Marine renewable energy for Arctic observations

Arctic observations are becoming increasingly valuable as researchers investigate climate change and its associated concerns, such as decreasing sea ice and increasing ship traffic. Networks of sensors with frequent sampling capabilities are needed to run forecast models, improve navigation, and inform climate research. Sampling frequency and deployment duration are currently constrained by battery power limitations. In-situ power generation using marine renewable energy sources such as waves and currents can be used to circumvent this constraint. Wave and current resources vary spatially and temporally in the Arctic, with some locations and seasons being better suited for marine renewable energy power generation. Locations and seasons with small resources may still be able to use marine renewable energy because of the low power requirements of the instruments. In this study, we describe the wave and current resources in the Arctic, outline the electricity generation developments that are needed to utilize the resources, and suggest use cases. Wave and current energy converters developed to power observations in the Arctic could also be used to power observations at lower latitudes. Marine renewable energy has the potential to decrease dependence on batteries and improve data collection capabilities in the Arctic; however, this would require the development of new low power technologies that can operate in extreme Arctic environments.

16 TIDAL AND WAVE POWER↗

An Experimental and Theoretical Study of the High Temperature Reactions of the Four Butyl Radical Isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl,sec-butyl, iso-butyl, andtert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH 3 or C 2 H 5 to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T< 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. Furthermore, this study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. Overall, the results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Free-Standing Boron-Doped Diamond Grid Electrode for Fundamental Spectroelectrochemistry

Spectroelectrochemistry (SEC) is a powerful technique that enables a variety of redox properties to be studied, including formal potential (E o ), thermodynamic values (ΔG, ΔH, ΔS), diffusion coefficient (D), electron transfer stoichiometry (n), and others. SEC requires an electrode which light can pass through while maintaining sufficient electrical conductivity. This has been traditionally composed of metal or metal oxide films atop transparent substrates like glass, quartz, or metallic mesh. Robust electrode materials like boron-doped diamond (BDD) could help expand the environments in which SEC can be performed, but most designs are limited to thin films (~100–200 nm) on transparent substrates less resilient than free-standing BDD. Here, this work presents a free-standing BDD grid electrode (G-BDD) for fundamental SEC measurements, using the well-characterized Fe(CN) 6 3–/4– redox couple as proof-of-concept. With a combination of cyclic voltammetry (CV), thin-layer SEC, and chronoabsorptometry, several of the redox properties mentioned above were calculated and compared. For E o' , n, and D, similar results were obtained when comparing the CV [E o' = +0.279 (±0.002) V vs Ag/AgCl; n = 0.97; D = 4.1 × 10 –6 cm 2 ·s –1 ] and SEC [E o' = +0.278 (±0.001) V vs Ag/AgCl; n = 0.91; D = 5.2 × 10 –6 cm 2 ·s –1 ] techniques. Both values align with what has been previously reported. To calculate D from the SEC data, modification of the classical equation used in chronoabsorptometry was required to accommodate the G-BDD electrode geometry. Overall, this work expands on the applicability of SEC techniques and BDD as a versatile electrode material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic response of additively manufactured Ti–5Al–5V–5Mo–3Cr as a function of heat treatment

Both shock and shockless compression experiments were performed on laser powder bed fusion (LPBF) Ti–5Al–5V–5Mo–3Cr (Ti-5553) to peak compressive stresses near 15 GPa. Experiments were performed on the as-built material, containing a purely β (body centered cubic) microstructure, and two differing heat treatments resulting in a dual phase α (hexagonal close packed) and β microstructure. The Hugoniot, Hugoniot elastic limit (HEL), and spallation strength were measured and compared to wrought Ti-6Al-4V (Ti-64). The results indicate the LPBF Ti-5553 Hugoniot response is similar between heat treatments and to Ti-64. The HEL stress observed in the LPBF Ti-5553 was considerably higher than Ti-64, with the as-built, fully β alloy exhibiting the largest values. The spallation strength of the LPBF Ti-5553 was also similar to Ti-64. Clear evidence of initial porosity serving as initiation sites for spallation damage was observed when comparing computed tomography measurements before and after loading. Post-mortem scanning electron microscopy images of the recovered spallation samples showed no evidence of retained phase changes near the spall plane. The spall plane was found to have kinks aligned with the loading direction near areas with large concentrations of twin-like, crystallographic defects in the as-built condition. For the heat-treated samples, the concentrations of twin-like, crystallographic defects were absent, and no preference for failure at the interface between the α and β phases was observed.

36 MATERIALS SCIENCE↗

DECOVALEX-2023: Task F1 Final Report

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 75 th percentile for all benchmarks. At the 90 th 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.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Mercury Speciation of Hanford 241-AP-107 Tank Waste Samples - 20384

The concentration of key mercury species in the Hanford tank waste has not been previously measured. Based on process knowledge, an estimated inventory of 2,100 kg of mercury is assumed and was distributed in varied amounts across 149 single-shell tanks (SSTs) and 28 double-shell tanks (DSTs). Limited analyses of individual SSTs and DSTs for total mercury were all below detection, yielding upper limits in concentration estimates based solely on instrument detection limits, resulting in an uncertain mercury inventory for Hanford tanks. Accurate and reliable chemical speciation and quantification is desired to anticipate, predict, and abate the mass movements of various mercury species through the Hanford waste processing flowsheet. Highly sensitive separation and quantification methods for total, elemental, and monomethyl mercury species, in tank waste, have been developed at the Pacific Northwest National Laboratory's (PNNL) Radiochemical Processing Laboratory (RPL). A non-radiological environmental method to quantify mercury species has been adapted for application to Hanford tank 241-AP-107 raw supernatant feed, filtered feed, and cesium-decontaminated samples. These samples represent feed inventory and pretreatment process effluents prior to immobilization by vitrification. This method separates and preconcentrates total, elemental, and monomethyl mercury onto selective solid substrates with little-to-no radioactive background. The solid samples can then be released from strict radiological control, allowing for analytical work to be conducted in low-level radiological facilities. Separation of mercury species was demonstrated at PNNL's RPL in Richland, Washington before transport and analytical quantification using cold vapor atomic fluorescence spectrometry (CVAFS) at PNNL's Marine Science Laboratory (MSL) in Sequim, Washington. Additional studies were conducted to assess the stability of monomethyl mercury in distilled matrices and preconcentrated onto the Carbotrap{sup R} solid substrate to assess sample degradation over time during transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Improving and Automating Building Model Data Exchange

There are many instances throughout a project’s lifecycle where there arises a need for quick and accurate risk assessment of building designs. For example, an unexpected design change during construction may necessitate structural engineers to perform a seismic risk assessment on analytical models of the updated building design using high fidelity structural analysis software, such as ANSYS or Abaqus. However, the efficiency of such workflows often depends upon the interoperability of architectural design software and structural analysis software. When the quality of this interoperability is lacking or even non-existent, the efficiency of virtual engineering workflows is hampered, which increases project costs. A McGraw Hill industry survey of professional users of Building Information Modeling (BIM) technologies found that there is high demand for BIM interoperability for structural analysis, but that the value/difficulty ratio is currently too low for practical use. There have been efforts by the academic community to facilitate model data exchange between the architectural design and structural analysis domains, but such solutions have not been widely adopted by industry, face technical challenges, and oftentimes are limited in applicability for users of various BIM software. Therefore, INL is developing capabilities to improve, automate, and generalize model data exchange between architectural BIM software (e.g., Revit) and structural analysis software (e.g., SAP2000, ANSYS). The goal is to help expedite and automate as much of the pre-processing step for creating analytical models in finite element analysis software as reasonably as possible. Such a "BIM-to-FEA" conversion tool should provide direct benefit to end-users through accuracy, automation, quick turn-around, and wide applicability. To generalize the application of this BIM-to-FEA conversion tool and increase its useability among the many different commercial BIM software currently used by industry, the program is being developed with the concept of openBIM. OpenBIM is the application of non-proprietary, open data standards that allow for BIM model data exchange in a format that is accessible, retainable, and useable for all users. The most widely used open, non-proprietary data exchange format for BIM is the Industry Foundation Classes (IFC) schema. IFC is developed by buildingSMART international and is ISO certified (ISO 16739-1:2018). The BIM-to-FEA conversion tool is being developed for compatibility with typical commercial building designs of steel framed structures. The tool is currently capable of importing architectural BIM data of framed building structures, recognizing and extracting the aspects of the model that are required for structural analysis, adjusting the connectivity of frame members, and finally exporting to an analytical model stored in the IFC format. The exported IFC analytical model can then be imported into various openBIM compliant software, such as SAP2000. Such capabilities have already been tested on commercial software, as shown above, and continue to be improved. Work is underway to test the conversion on various commercial BIM software, develop a user-friendly interface, incorporate the program into the broader DeepLynx data warehouse project being developed by INL, and to eventually open-source the tool for the benefit of the community. Future development of the tool envisions the ability for efficient iterative risk assessment of generative building designs, all within a workflow utilizing open-source tools. One such open-source tool will be MOOSE, an advanced finite element analysis tool developed at INL. The conversion tool will also branch out from typical commercial building designs and will aim to incorporate nuclear construction. The aim will be to convert both structural and non-structural components of nuclear facilities, such as curved concrete containment structures and piping systems, respectively.

97 MATHEMATICS AND COMPUTING↗