Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic separation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Radiation-Induced Effects on the Extraction Properties of Hexa-n-octylnitrilo-triacetamide (HONTA) Complexes of Americium and Europium

The candidate An(III)/Ln(III) separation ligand hexa-n-octylnitrilo-triacetamide (HONTA) was irradiated under envisioned SELECT (Solvent Extraction from Liquid waste using Extractants of CHON-type for Transmutation) process conditions (n-dodecane/0.1 M HNO 3 ) using a solvent test loop in conjunction with cobalt-60 gamma irradiation. The extent of HONTA radiolysis and complimentary degradation product formation was quantified by HPLC-ESI-MS/MS. Further, the impact of HONTA radiolysis on process performance was evaluated by measuring the change in 243 Am and 154 Eu distribution ratios as a function of absorbed gamma dose. HONTA was found to decay exponentially with increasing dose, affording a dose coefficient of d = (4.48 ± 0.19) × 10 –3 kGy –1 . Multiple degradation products were detected by HPLC-ESI-MS/MS with dioctylamine being the dominant quantifiable species. Both 243 Am and 154 Eu distribution ratios exhibited an induction period of ~70 kGy for extraction (0.1 M HNO 3 ) and back-extraction (4.0 M HNO 3 ) conditions, after which both values decreased with absorbed dose. The decrease in distribution ratios was attributed to a combination of the destruction of HONTA and ingrowth of dioctylamine, which is capable of interfering in metal ion complexation. Furthermore, the loss of HONTA with absorbed gamma dose was predominantly attributed to its reaction with the n-dodecane radical cation (R˙ + ). These R˙ + reaction kinetics were measured for HONTA and its 241 Am and 154 Eu complexes using picosecond pulsed electron radiolysis techniques. All three second-order rate coefficients (k) were essentially diffusion limited in n-dodecane indicating a significant reaction pathway: k(HONTA + R˙ + ) = (7.6 ± 0.8) × 10 9 M –1 s –1 , k(Am(HONTA) 2 + R˙ + ) = (7.1 ± 0.7) × 10 10 M –1 s –1 , and k(Eu(HONTA) 2 + R˙ + ) = (9.5 ± 0.5) × 10 10 M –1 s –1 . HONTA-metal ion complexation afforded an order-of-magnitude increase in rate coefficient. Additional nanosecond time-resolved measurements showed that both direct and indirect HONTA radiolysis yielded the short-lived (< 100 ns) HONTA radical cation and a second long-lived (µs) species identified as the HONTA triplet excited state. The latter was confirmed by a series of oxygen quenching picosecond pulsed electron measurements, affording a quenching rate coefficient of k( 3 [HONTA]* + O 2 ) = 2.2 × 10 8 M –1 s –1 . Overall, both the HONTA radical cation and triplet excited state are important precursors to the suite of measured HONTA degradation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Recent Advances in Nonclassical Crystallization: Fundamentals, Applications, and Challenges

We report crystals are ubiquitous and play an important role in our daily life. Recent studies on properties and formation of crystals have proven the attachment of particle ranging from ion pairs to well-crystallized nanocrystal is one of the most important ways for crystallization. Compared with the monomer-by-monomer (i.e., atoms, ions, or molecules) addition described in longstanding classical theories, so-called nonclassical crystallization pathways are more complex as the non-specific interaction of variation in system free energy and reaction kinetics. In view of the detailed insights into nonclassical crystallization pathways has broadened the scope of research across various disciplines such as geological events, biomineralization mechanisms, environmental remediation, and the development of functional materials, this paper reviews the current state of the art in nonclassical nucleation (the prenucleation clusters pathway, aggregation pathway, etc.) and growth pathways (oriented attachment and random attachment) based on several spotlights. In addition, the observed specific phenomena on crystallization in the last 10 years as well as recent advancements and achievements made on the synthesis of mesocrystals are also elucidated. With implications and perspectives of nonclassical crystallization researches proposed in this paper, this review aims to summarize research advancement regarding crystallizing systems. Emphasis will be placed on novel insights into nonclassical crystallization of materials, which are expected to provide completely new perspectives on the underlying crystallization mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Overcoming Fundamental Limitations in Adsorbent Design: Alkene Adsorption by Non-porous Copper(I) Complexes

Purifying alkenes from alkanes requires cryogenic distillation. This consumes energy equivalent to countries of ca. 5 million people. Replacing distillation with adsorption processes would significantly increase energy efficiency. Trade-offs between kinetics, selectivity, capacity, and heat of adsorption have prevented production of an optimal adsorbent. Here, we report adsorbents that overcome these trade-offs. [Cu-Br] 3 and [Cu-H] 3 are air-stable trinuclear complexes that undergo reversible solid-state inter-molecular rearrangements to produce dinuclear [Cu-Br∙(alkene)] 2 and [Cu-H∙(alkene)] 2 . The reversible solid-state rearrangement, confirmed in situ using powder X-ray diffraction, allows adsorbent design trade-offs to be overcome, coupling low heat of adsorption (-10 to -17 kJ mol -1 alkene ), high alkene:alkane selectivity (47; 29), and uptake capacity (>2.5 molalkene mol -1 Cu3 ). Most remarkably, [Cu-H] 3 displays fast uptake and regenerates capacity within 10 minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Selectivity of Electrochemical Ion Insertion into Manganese Dioxide Polymorphs

The ion insertion redox chemistry of manganese dioxide has diverse applications in energy storage, catalysis, and chemical separations. Unique properties derive from the assembly of Mn–O octahedra into polymorphic structures that can host protons and nonprotonic cations in interstitial sites. Despite many reports on individual ion-polymorph couples, much less is known about the selectivity of electrochemical ion insertion in MnO 2 . In this work, we use density functional theory to holistically compare the electrochemistry of A x MnO 2 (where A = H + , Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ , Al 3+ ) in aqueous and nonaqueous electrolytes. We develop an efficient computational scheme demonstrating that Hubbard-U correction has a greater impact on calculating accurate redox energetics than choice of exchange-correlation functional. Using PBE+U, we find that for nonprotonic cations, ion selectivity depends on the oxygen coordination environments inside a polymorph. When H + is present, however, the driving force to form hydroxyl bonds is usually stronger. In aqueous electrolytes, only three ion-polymorph pairs are thermodynamically stable within water’s voltage stability window (Na + and K + in α-MnO 2 , and Li + in λ-MnO 2 ), with all other ion insertion being metastable. We find Al 3+ may insert into the δ, R, and λ polymorphs across the full 2-electron redox of MnO 2 at high voltage; however, electrolytes for multivalent ions must be designed to impede the formation of insoluble precipitates and facilitate cation desolvation. We also show that small ions coinsert with water in α-MnO 2 to achieve greater coordination by oxygen, while solvation energies and kinetic effects dictate water coinsertion in δ-MnO 2 . Finally, taken together, these findings explain reports of mixed ion insertion mechanisms in aqueous electrolytes and highlight promising design strategies for safe, high energy density electrochemical energy storage, desalination batteries, and electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2021

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Illuminating the Photophysics of Oxygen Atom Transfer in Rare Earth Metal-Organic Complexes

Industrial separations of rare earth metals are reliant on selectivity based on ionic radii, which are time, energy, and resource consuming.1 The unique and discontinuous manifold of energy states associated with the f-electron configurations of the rare earth ions2 offer an opportunity to overcome some of these limitations by incorporating light into the separation process. Here, we use transient absorption spectroscopy to understand the differences in mechanism of light-initiated oxygen atom transfer (OAT) in yttrium (Y) and dysprosium (Dy) rare earth complexes. These metal-organic complexes are functionalized with hexfluoroacetyl-acetonate (hfac) and 4-methylmorpholine-N-oxide (NMMO). After light absorption in the presence of triphenyl phosphine (TPP), the NMMO in the Y-NMMO complex is replaced by the TPPO OAT reaction product. Because the Dy f-electron manifold is accessible to the ligand excited states, while there are no accessible states associated with f-electrons for Y, there are significant differences in reactivity. Transient absorption spectroscopy (TAS) was used to monitor changes after 340 nm photoexcitation of Dy/Y-NMMO or the product complex Dy/Y-TPPO. TAS is a pump-probe technique where a white light probe pulse is delayed relative to a visible pump pulse, which photoexcites a sample. Differences between the ground and excited states are monitored over time. With TAS, we access time scales from 100 fs-400 us to determine the photophysics of the complexes of interest. Photoexcitation of Y-NMMO (Fig. 1A) creates a broad photoinduced absorption (PIA) that evolves into a sharp positive feature over 4.2 ps centered at 450 nm. A secondary feature centered at 550 nm grows in with a 200 ps time constant and remains through the 5 ns time window of our experiment. Control experiments for Y-TPPO reveal that the feature at 450 nm is due to internal dynamics of the hfac ligand, while the 550 nm shoulder is due to the presence of the NMMO. Comparing kinetics at 550 nm for Y-NMMO and Dy-NMMO (Fig. 1B) shows that while Y exhibits a distinct growth in the NMMO feature, the Dy decays almost completely within 5 ns. Control experiments with Dy-TPPO prove that the ligand excited states centered on the hfac ligand are deactivated via energy transfer to the Dy manifold within 200 ps. This deactivation leads to significant differences in OAT reactivity between Y and Dy.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

We report quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyI)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI 2 (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI 2 ∙- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI 2 ∙- . Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the 3* ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Influence of americium complexation on the radiation-induced chemical reactivity of sulfophenyl bistriazinyl pyridine (SO 3 -Ph-BTP) towards the nitrate radical

Sulfophenyl bistriazinyl pyridine (SO 3 -Ph-BTP) is a hydrophilic organic ligand used to separate actinides and lanthanides. Dose accumulation and time-resolved studies have previously provided insight into the radiolytic stability and degradation pathways of SO 3 -Ph-BTP in reprocessing environments, but no study has yet addressed the impact of minor actinide complexation on the radiation chemistry of this ligand. To begin to fill this knowledge gap, a systematic, time-resolved study exploring the reactivity of the nitrate radical (NO 3 • ) with SO 3 -Ph-BTP in the presence of trivalent americium, Am(III), has been conducted, which demonstrates enhanced reactivity (an order of magnitude faster) upon metal complexation.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2022

The role of surface forces in environment-enhanced cracking of brittle solids

Fracture initiation and propagation in brittle materials is promoted in surface-reactive (sorptive) environments, a phenomenon known as subcritical crack growth (SCG). Laboratory measured crack-propagation velocity vs. stress intensity factor relationships typically exhibit highly nonlinear, multi-stage characteristics that are sensitive to environmental factors such as adsorbate concentration and temperature. For practical purposes, empirical relationships (e.g., a power law) have been used to describe this complex phenomenon. However, how the overall SCG behavior emerges from the underlying fundamental processes near the crack tip, such as the interaction of the crack surfaces separated by only a few nanometers and mass transport within the nano-confined space, is still not well understood. This paper develops a mechanistic, surface-force-based fracture theory (SFFT) which integrates surface force models, fluid transport models, and linear elastic fracture mechanics to quantitatively explain the multi-stage characteristics of SCG in brittle solids. A numerical model is developed based on SFFT and solved through an implicit partitioned scheme for efficiency and modularity. The results are validated by Wiederhorn's data on crack propagation in soda-lime glasses at a wide range of relative humidity levels. We show that, for the first time, the entire range of an SCG curve can be captured by a single physics-based model. The predicted SCG curves reveal that the development of repulsive disjoining pressure behind the crack tip can be responsible for the reduced apparent fracture toughness in a sorptive environment. In conclusion, the shape of the SCG curve, and its changes with respect to the environment, is found to critically depend on the assumed transport models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Savannah River Site High Level Radioactive Waste Tank Sample Thermolytic Hydrogen Generation - 20346

To support the revision of flammability calculations and controls for the Savannah River Site (SRS) High-Level Radioactive Waste (HLW) tanks, Savannah River National Laboratory (SRNL) conducted laboratory measurements with the goal of quantifying hydrogen produced in HLW by non-radiolytic chemical reactions (i.e., thermolysis). Testing was performed with HLW tank supernate samples that represented a cross section of types of waste at SRS, including tanks that store typical evaporator concentrate, dissolved salt-cake, dilute recycle stream from the Defense Waste Processing Facility (DWPF), evaporator concentrate of the recycle stream from DWPF, fresh waste from the canyon separations facility, and waste that has been through salt processing cesium removal. Non-radioactive simulant testing included the classes of organic compounds historically introduced into the HLW tanks that were recently shown to be the most active toward the thermolytic generation of hydrogen. The reaction rate equations developed for thermolytic hydrogen generation from each class of organic compound showed direct proportionality to organic compound concentration and hydroxide concentration. Simulant testing identified that highly concentrated waste with high hydroxide concentration had the highest rates of thermolytic hydrogen generation. Using the radioactive tank sample thermolysis measurements at high temperatures and the mechanistic salt dependence of the simulant models, a global model was developed for tank waste thermolytic hydrogen generation as a function of total organic carbon content of SRS HLW. Through the primary functionality of the global model (organic carbon and free hydroxide concentration), the relative reactivity of the organic carbon in the waste was quantified. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Imaging Phase Segregation in Nanoscale Li x CoO 2 Single Particles

Li x CoO 2 (LCO) is a common battery cathode material that has recently emerged as a promising material for other applications including electrocatalysis and as electrochemical random access memory (ECRAM). During charge– discharge cycling LCO exhibits phase transformations that are significantly complicated by electron correlation. While the bulk phase diagram for an ensemble of battery particles has been studied extensively, it remains unclear how these phases scale to nanometer dimensions and the effects of strain and diffusional anisotropy at the single-particle scale. Understanding these effects is critical to modeling battery performance and for predicting the scalability and performance of electrocatalysts and ECRAM. Here we investigate isolated, epitaxial LiCoO 2 islands grown by pulsed laser deposition. After electrochemical cycling of the islands, conductive atomic force microscopy (c-AFM) is used to image the spatial distribution of conductive and insulating phases. Above 20 nm island thicknesses, we observe a kinetically arrested state in which the phase boundary is perpendicular to the Li-planes; we propose a model and present image analysis results that show smaller LCO islands have a higher conductive fraction than larger area islands, and the overall conductive fraction is consistent with the lithiation state. Thinner islands (14 nm), with a larger surface to volume ratio, are found to exhibit a striping pattern, which suggests surface energy can dominate below a critical dimension. When increasing force is applied through the AFM tip to strain the LCO islands, significant shifts in current flow are observed, and underlying mechanisms for this behavior are discussed. The c-AFM images are compared with photoemission electron microscopy images, which are used to acquire statistics across hundreds of particles. Finally, the results indicate that strain and morphology become more critical to electrochemical performance as particles approach nanometer dimensions.

25 ENERGY STORAGE↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Enabling Microscale Processing: Combined Raman and Absorbance Spectroscopy for Microfluidic On-Line Monitoring

Microfluidics have many potential applications including characterization of chemical processes on a reduced scale, spanning the study of reaction kinetics using on-chip liquid–liquid extractions, sample pretreatment to simplify off-chip analysis, and for portable spectroscopic analyses. The use of in situ characterization of process streams from laboratory-scale and microscale experiments on the same chemical system can provide comprehensive understanding and in-depth analysis of any similarities or differences between process conditions at different scales. A well-characterized extraction of Nd(NO 3 ) 3 from an aqueous phase of varying NO 3– (aq) concentration with tributyl phosphate (TBP) in dodecane was the focus of this microscale study and was compared to an earlier laboratory-scale study utilizing counter current extraction equipment. Here, we verify that this same extraction process can be followed on the microscale using spectroscopic methods adapted for microfluidic measurement. Concentration of Nd (based on UV–vis) and nitrate (based on Raman) was chemometrically measured during the flow experiment, and resulting data were used to determine the distribution ratio for Nd. Extraction distributions measured on the microscale were compared favorably with those determined on the laboratory scale in the earlier study. Both micro-Raman and micro-UV–vis spectroscopy can be used to determine fundamental parameters with significantly reduced sample size as compared to traditional laboratory-scale approaches. This leads naturally to time, cost, and waste reductions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗