SEARCH · Engineering Papers

Results for “Agreements”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗

DOE OSTI2024

Impact of LWR assembly structural features on cladding burst behavior under LOCA conditions

This provides an initial scoping study on clad balloon and burst behavior for burnup extension of reactor fuel. The associated issues with burnup extension are fuel fragmentation, relocation, and dispersal in the event of cladding failure. The general finding of this work is that the structural features, spacer grids and mixing vanes, locally suppress cladding deformation but have little impact on the overall clad performance during loss-of-coolant accidents. The work detailed in a previous report by Capps et al. focused on core optimization via neutronics, thermal hydraulics and thermomechanical analysis for burnups beyond 62 GWD/tU and enrichments above 5%. Uncertainty of fuel fragmentation relocation and dispersal in high-burnup rods during accident conditions was also investigated. The dispersal aspect of fuel fragmentation depends on cladding rupture. Thus, assessing uncertainties in the rupture behavior is helpful in estimating the dispersal of high-burnup fuel. This study builds on the previous work by assessing the impact of assembly structural features on cladding balloon and burst behavior in a full-length fuel rod. In this work, the BISON fuel performance code was used to generate 2D radial and height meshes containing structural features commonly used in nuclear fuel assemblies. First, meshes were generated with spacer grids. Results were then compared to the cladding burst temperature and balloon strain results from the previous work. A mesh sensitivity study was performed to ensure that mixing vanes and spacer grid effects were appropriately considered, resulting in a more refined mesh than the previous study. The balloon deformation and burst times of the cladding were compared to the original case. Consideration was also given to the effect of rod initial pressure. In conclusion, 3D quarter rod simulations were also performed and found good agreement with the 2D simulations in clad deformation and reasonable agreement in burst times.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

DOE OSTI2022

Measurement of cold nuclear matter effects for inclusive J / ψ in p +Au collisions at s NN = 200 GeV

Measurement by the STAR experiment at RHIC of the cold nuclear matter (CNM) effects experienced by inclusive J / ψ at mid-rapidity in 0-100% p +Au collisions at s NN = 200 GeV is presented. Such effects are quantified utilizing the nuclear modification factor, R p Au , obtained by taking a ratio of J / ψ yield in p +Au collisions to that in p + p collisions scaled by the number of binary nucleon-nucleon collisions. The differential J / ψ yield in both p + p and p +Au collisions is measured through the dimuon decay channel, taking advantage of the trigger capability provided by the Muon Telescope Detector in the RHIC 2015 run. Consequently, the J / ψ R p Au is derived within the transverse momentum ( p T ) range of 0 to 10 GeV/ c . A suppression of approximately 30% is observed for p T < 2 GeV/ c , while J / ψ R p Au becomes compatible with unity for p T greater than 3 GeV/ c , indicating the J / ψ yield is minimally affected by the CNM effects at high p T . Comparison to a similar measurement from 0-20% central Au+Au collisions reveals that the observed strong J / ψ suppression above 3 GeV/ c is mostly due to the hot medium effects, providing strong evidence for the formation of the quark-gluon plasma in these collisions. Several model calculations show qualitative agreement with the measured J / ψ R p Au , while their agreement with the J / ψ yields in p + p and p +Au collisions is worse.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2022

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Field-scale dynamics of planting dates in the US Corn Belt from 2000 to 2020

Crop planting dates are a dynamic feature of agricultural systems that respond to short- and long-term climate signals, crop and cultivar selection, and technology changes. Planting date records are essential for yield gap analyses, accurate crop modeling, and tracking farmer adaptations to weather and climate change. Although planting dates have high variation at local scales due to heterogeneity in farm resources and decision-making, available long-term data on planting dates is largely restricted to aggregated regional statistics or, at best, satellite-derived datasets with limited spatiotemporal extent and at resolutions unable to distinguish individual fields (> 250 m). Here, we generated retrospective annual field-scale (30 m) planting date maps for both maize and soybeans spanning 2000-2020 across a 12 state region in the United States Corn Belt based on Landsat satellite data and a large ground sample of over 28,000 maize and soybean fields. Using training data from 2015-2020 for model selection, we found that planting date predictions improved with harmonic regression of Landsat data and additional annual weather covariates. The preferred random forests model approximately doubled performance compared to a null model based on state median planting dates, capturing 47% of field-level variation for maize (mean absolute error, MAE = 7.4 days) and 44% for soybeans (MAE = 7.5 days) against held-out ground truth test data for 2008-2014. We also evaluated the full 2000-2020 dataset with state agricultural statistics, finding strong agreement with median planting dates for maize (R 2 = 0.76, MAE = 4.4 days) and slightly lower agreement for soybeans (R 2 = 0.65, MAE = 5.4 days) when aggregated to the state level. We then used this new dataset to analyze environmental determinants of planting dates at a finer-scale than previously possible, controlling for unobserved variation at the sub-state district level. We found that during 2000-2020, each standard deviation increase in rainfall delayed planting by ~ 2.5 days, and fields with higher soil productivity ratings tended to be planted earlier. We did not find meaningful trends over the last two decades in planting dates for maize or soybeans, in contrast to trends towards earlier planting dates late last century and predicted for this period in climate adaptation studies. We hypothesize increases in early season rainfall may have inhibited these shifts towards earlier planting. Remotely sensed planting dates will be a useful tool for yield gap analyses, crop simulation modeling, and ongoing assessment of climate adaptation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Experimental comparison of pyranometer, reflectometer, and spectrophotometer methods for the measurement of roofing product albedo

Albedo (solar reflectance) can be measured outdoors with a pyranometer or indoors with a hemispherical reflectometer or a spectrophotometer. The current study evaluates these methods and their applicability to roofing materials by measuring and comparing the ASTM E1918 (pyranometer), non-ASTM E1918A (alternative pyranometer), ASTM C1549 (reflectometer), and E903 (spectrophotometer) albedos of 10 roofing products, including three single-ply membranes, one asphalt shingle, three roofing aggregates, and three high-profile tiles. It uses full-size (4 m × 4 m) assemblies in the E1918 and E1918A trials; corrects the E1918 albedos to remove shadow and background errors; and evaluates C1549 and E903 albedos with two different irradiance spectra, one global horizontal and the other beam normal. E1918A albedos matched E1918 corrected (E1918_cor) albedos to within 0.036. Furthermore, agreement between C1549 air mass 1 global horizontal (C1549_G1) albedo and E1918_cor albedo was within 0.015 for membrane and shingle coupons, and within 0.031 for aggregates, though the uncertainty in the latter agreement was up to 0.043. After minor corrections, C1549_G1 albedos of flat or slightly convex tile chips were 0.032–0.052 higher than the E1918_cor albedos of their corresponding high-profile tile assemblies because high-profile surfaces have concavities. Switching to the beam-normal albedo C1549_1.5E increased the C1549 albedos of the 10 tested products by 0.004–0.054, with the largest increases (C1549_1.5E – C1549_G1) accruing to spectrally selective cool colors. Using C1549_1.5E albedo to characterize a chip of a cool-colored tile compounded the error induced by representing a high-profile surface by a flat or slightly convex specimen.

14 SOLAR ENERGY↗

DOE OSTI2024

Heat transfer coefficients of moving particle beds from flow-dependent thermal conductivity and near-wall resistance

Accurate determination of heat transfer coefficients for flowing packed particle beds is essential to the design of particle heat exchangers and other thermal and thermochemical equipment. While such dense granular flows mostly fall into the well-known plug-flow regime, the discrete nature of granular materials alters the thermal transport processes in both the near-wall and bulk regions of flowing particle beds from their stationary counterparts. As a result, heat transfer correlations based on the stationary particle bed thermal conductivity could be inadequate for flowing particles in a heat exchanger. Most earlier works have achieved a reasonable agreement with experiments by treating granular heat transfer media as a plug-flow continuum with a near-wall thermal resistance in series. However, the thermal conductivity values of the continuum were often obtained from measurements on stationary beds owing to the difficulty of flowing bed measurements. In this work, it was found that the properties of a stationary bed are highly sensitive to the method of particle packing and there is a decrease in the particle bed thermal conductivity and increase in the near-wall thermal resistance, measured as an effective air gap thickness, on the onset of particle flow. These variations in thermal conductivity of stationary and flowing particle beds can lead to errors in heat transfer coefficient calculations. Therefore, the heat transfer coefficients for granular flows were calculated using experimentally determined flowing particle bed thermal conductivity and near-wall air gap for ceramic particles – CARBO CP 40/100 (mean diameter = 275 µm), HSP 40/70 (404 µm) and HSP 16/30 (956 µm); at velocities of 5–15 mm·s –1 ; and temperatures of 300–650 °C. The thermal conductivity and air gap values for CP 40/100 and HSP 40/70 were further used to calculate heat transfer coefficients across different particle bed temperatures and velocities for different parallel-plate heat exchanger dimensions. Furthermore, these calculations, which show good agreement with measured HTC values reported in literature, can be used as a guide for heat exchanger designs. Graphical abstract

14 SOLAR ENERGY↗

DOE OSTI2023

Optimal white-noise stochastic forcing for linear models of turbulent channel flow

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.

Mechanics↗

DOE OSTI2021

Modelling of parallel dynamics of a pellet-produced plasmoid

The parallel expansion of a dense, pellet-produced plasmoid is modelled with parameters relevant to pellet fuelling experiments in the Wendelstein7-X stellarator. Good agreement is found between the analytical theory and more detailed modelling. In particular, much of the energy deposited in the pellet by the ambient plasma is transferred to the pellet ions by the ambipolar electric field during the expansion. The validity of the hydrodynamic treatment of the plasmoid and the ambient plasma is discussed.

Physics↗

DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗

DOE OSTI2023

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Radiation-Enhanced Anion Transport in Hematite

The influence of radiation-induced (1 MeV energy H+ to ~0.1 dpa at 450°C), non-equilibrium point defect populations on mass transport is studied with an integrated campaign of experimental and theoretical methods. Using epitaxial thin films of hematite with embedded 18O tracer layers and nanoscale atom probe tomography measurements, it is shown that anion self-diffusion is enhanced by at least 2 orders of magnitude under irradiation compared to thermal diffusion alone. Complementary scanning transmission electron microscopy of vacuum annealed specimens, reveals associated microstructural changes in the oxide films, including local phase transformation to Fe3O4 and the development of nanoscale voids from vacancy coalescence. Point defect formation and migration energies were computed from density functional theory and applied within the context of chemical rate theory to analyze contributions from both interstitial and vacancy mechanisms to self-diffusion in thermal and irradiation conditions. Comparisons are made between calculated, literature and newly measured self-diffusion values, revealing good agreement on the magnitude of radiation-enhanced anion diffusion. Further, the model suggests a transition from vacancy to interstitialcy mechanisms at low temperatures and high oxygen activity, consistent with the varied activation energies reported from prior studies.

iron oxides, isotopic tracers, atom probe tomograp↗

DOE OSTI2024

A Detailed Reaction Mechanism for Thiosulfate Oxidation by Ozone in Aqueous Environments

The ozone oxidation, or ozonation, of thiosulfate is an important reaction for wastewater processing, where it is used for remediation of mining effluents, and for studying aerosol chemistry, where its fast reaction rate makes it an excellent model reaction. Although thiosulfate ozonation has been studied since the 1950s, challenges remain in developing a realistic reaction mechanism that can satisfactorily account for all observed products with a sequence of elementary reaction steps. Here, we present novel measurements using trapped microdroplets to study the pH-dependent thiosulfate ozonation kinetics. We detect known products and intermediates, including SO32-, SO42-, S3O62-, and S4O62-, establishing agreement with the literature. However, we identify S2O42- as a new reaction intermediate and find that the currently accepted mechanism does not directly explain observed pH effects. Thus, we develop a new mechanism, which incorporates S2O42- as an intermediate and uses elementary steps to explain the pH dependence of thiosulfate ozonation. The proposed mechanism is tested using a kinetic model benchmarked to the experiments presented here, then compared to literature data. We demonstrate good agreement between the proposed thiosulfate ozonation mechanism and experiments, suggesting that the insights in this paper can be leveraged in wastewater treatment and in understanding potential climate impacts.

Deal, Alexandra M↗

DOE OSTI2020

Probing the Electronic Structure of a Thorium Nitride Complex by Solid-State 15N NMR Spectroscopy

The solid-state 15N NMR powder spectra of the thorium nitride complex, [K(18-crown-6)(THF)2][(R2N)3Th(m-15N)(Th(NR2)3] ([K][1-15N], R = SiMe3), and the thorium amide complex, [Th(NR2)3(15NH2)] (2-15N) were recorded. The spectrum for [K][1-15N] represents the first reported solid-state 15N NMR data for an actinide nitride complex. The experimentally measured tensor spans are 847 ppm for [K][1-15N] and 237 ppm for 2-15N. Both shielding tensors exhibit a near-zero asymmetry parameter, which for [K][1-15N] is consistent with a local rotational symmetry of its 15N-labelled nitride ligand. For 2-15N, the lack of asymmetry can be rationalized by a quasi-free Th-NH2 bond rotation in the solid-state. DFT calculations overestimate the tensor span somewhat for [K][1-15N], but provide isotropic shifts in good agreement with both the solid-state and solution values for both complexes. Natural localized molecular orbital (NLMO) analyses of the nuclear shielding reveal that the larger tensor span in [K][1-15N] vs. 2-15N is primarily a consequence of more pronounced covalency of the N-Th bonds, and large spin-orbit coupling due to significant Th 5f orbital contribution to those bonds, impacting the principal components of the shielding tensor perpendicular to the Th-N-Th axis. Overall, our analysis confirms the involvement of the 5f orbitals in Th-N multiple bonds, and further demonstrates the value of solid-state NMR spectroscopy for interrogating actinide-ligand bonding.

actinides, covalent bonding, Nuclear Magnetic Reso↗