Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vertical phase 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling of Sagging for 3D Printed Layers During the Curing Process

Here, a printed layer of silica-enforced poly(dimethylsiloxane)-co-(diphenylsiloxane) is modeled as a two-phase system consisting of air and polymer with an interface set up during the printing process. The structural geometry changes mostly due to the action of surface tension, while all material properties are strongly temperature dependent. Polymer flow is described using equations of the extended Herschel–Bulkley model, with parameters strongly dependent on temperature and degree of curing. Parameters of the model are determined using flow sweep measurements and separate experiments with vertical structure sagging at different temperatures. The curing process is modeled using dependencies between the curing rate, degree of curing, and temperature obtained in studies by differential scanning calorimetry. The developed model is used for simulations of printed structure deformations with different initial and boundary conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗
DOE OSTI2022

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗
DOE OSTI2022

Single-Step Fabrication of Au-Fe-BaTiO3 Nanocomposite Thin Films Embedded with Non-Equilibrium Au-Fe Alloyed Nanostructures

Nanocomposite thin film materials present great opportunities in coupling materials and functionalities in unique nanostructures including nanoparticles-in-matrix, vertically aligned nanocomposites (VANs), and nanolayers. Interestingly the nanocomposites processed through a non-equilibrium processing method, e.g., pulsed laser deposition (PLD), often possess unique metastable phases and microstructures that could not achieve using equilibrium techniques, and thus lead to novel physical properties. In this work, a unique three-phase system composed of BaTiO3 (BTO), with two immiscible metals, Au and Fe, is demonstrated. By adjusting the deposition laser frequency from 2 Hz to 10 Hz, the phase and morphology of Au and Fe nanoparticles in BTO matrix vary from separated Au and Fe nanoparticles to well-mixed Au-Fe alloy pillars. This is attributed to the non-equilibrium process of PLD and the limited diffusion under high laser frequency (e.g., 10 Hz). The magnetic and optical properties are effectively tuned based on the morphology variation. This work demonstrates the stabilization of non-equilibrium alloy structures in the VAN form and allows for the exploration of new non-equilibrium materials systems and their properties that could not be easily achieved through traditional equilibrium methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structure-Related Electronic and Magnetic Properties in Ultrathin Epitaxial Ni x Fe 3–x O 4 Films on MgO(001)

Off-stoichiometric Ni x Fe 3–x O 4 ultrathin films (x < 2.1) with varying Ni content x and thickness 16 (±2) nm were grown on MgO(001) by reactive molecular beam epitaxy. Synchrotron-based high-resolution X-ray diffraction measurements reveal vertical compressive strain for all films, resulting from a lateral pseudomorphic adaption of the film to the substrate lattice without any strain relaxation. Complete crystallinity with smooth interfaces and surfaces is obtained independent of the Ni content x. For x < 1 an expected successive conversion from Fe 3 O 4 to NiFe 2 O 4 is observed, whereas local transformation into NiO structures is observed for films with Ni content x > 1. However, angle-resolved hard X-ray photoelectron spectroscopy measurements indicate homogeneous cationic distributions without strictly separated phases independent of the Ni content, while X-ray absorption spectroscopy shows that also for x > 1, not all Fe 2+ cations are substituted by Ni 2+ cations. The ferrimagnetic behavior, as observed by superconducting quantum interference device magnetometry, is characterized by decreasing saturation magnetization due to the formation of antiferromagnetic NiO parts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2026

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Direct recovery of interfacial topography from coherent X-ray reflectivity: model calculations for a 1D interface

The use of coherent X-ray reflectivity to recover interfacial topography is described using model calculations for a one-dimensional interface. The results reveal that the illuminated topography can be recovered directly from the measured reflected intensities. This is achieved through an analysis of the Patterson function, the Fourier transform of the scattering intensity (as a function of lateral momentum transfer, Q // , at fixed vertical momentum transfer, Q z ). Specifically, a second-order Patterson function is defined that reveals the discrete set of separations and contrast factors (i.e., the product of changes in the effective scattering factor) associated with discontinuities in the effective interfacial topography. Finally, it is shown that the topography is significantly over-determined by the measurements, and an algorithm is described that recovers the actual topography through a deterministic sorting of this information.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Top-Down Strategy to Engineer ActiveLayer Morphology for Highly Efficient and Stable All-Polymer Solar Cells

A major challenge hindering the further development of all-polymer solar cells (all-PSCs) employing polymerized small-molecule acceptors is the relatively low fill factor (FF) due to the difficulty in controlling the active-layer morphology. The issues typically arise from oversized phase separation resulting from the thermodynamically unfavorable mixing between two macromolecular species, and disordered molecular orientation/packing of highly anisotropic polymer chains. Herein, a facile top-down controlling strategy to engineer the morphology of all-polymer blends is developed by leveraging the layer-by-layer (LBL) deposition. Optimal intermixing of polymer components can be achieved in the two-step process by tuning the bottom-layer polymer swelling during top-layer deposition. Consequently, both the molecular orientation/packing of the bottom layer and the molecular ordering of the top layer can be optimized with a suitable top-layer processing solvent. A favorable morphology with gradient vertical composition distribution for efficient charge transport and extraction is therefore realized, affording a high all-PSC efficiency of 17.0% with a FF of 76.1%. The derived devices also possess excellent long-term thermal stability and can retain >90% of their initial efficiencies after being annealed at 65 °C for 1300 h. In conclusion, these results validate the distinct advantages of employing an LBL processing protocol to fabricate high-performance all-PSCs.

14 SOLAR ENERGY↗
DOE OSTI2025

Unlocking enhanced gas capture via core scrambling of porous-organic cages

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Here, therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO 2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO 2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO 2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Gas capture↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

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

Twist Angle-Dependent Interlayer Exciton Lifetimes in van der Waals Heterostructures

In van der Waals (vdW) heterostructures formed by stacking two monolayers of transition metal dichalcogenides, multiple exciton resonances with highly tunable properties are formed and subject to both vertical and lateral confinement. We investigate how a unique control knob, the twist angle between the two monolayers, can be used to control the exciton dynamics. We observe that the interlayer exciton lifetimes in MoSe 2 /WSe 2 twisted bilayers (TBLs) change by one order of magnitude when the twist angle is varied from 1° to 3.5°. Additionally, using a low-energy continuum model, we theoretically separate two leading mechanisms that influence interlayer exciton radiative lifetimes. The shift to indirect transitions in the momentum space with an increasing twist angle and the energy modulation from the moiré potential both have a significant impact on interlayer exciton lifetimes. We further predict distinct temperature dependence of interlayer exciton lifetimes in TBLs with different twist angles, which is partially validated by experiments. While many recent studies have highlighted how the twist angle in a vdW TBL can be used to engineer the ground states and quantum phases due to many-body interaction, our studies explore its role in controlling the dynamics of optically excited states, thus, expanding the conceptual applications of “twistronics”.

36 MATERIALS SCIENCE↗
DOE OSTI2021

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Updraft dynamics and microphysics: on the added value of the cumulus thermal reference frame in simulations of aerosol–deep convection interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kinetic versus Thermodynamic Stability of LLZO in Contact with Lithium Metal

Li 7 La 3 Zr 2 O 12 (LLZO) garnet-based oxides are a promising class of solid electrolytes used as the separator in all-solid-state batteries (ASSBs). While LLZO is considered to have a wide electrochemical stability window, its intrinsic stability in contact with lithium metal is not sufficiently well understood, and there is still a debate on the key question of whether LLZO does or does not form passivation layers before and during cycling. Utilizing both in situ and operando X-ray photoelectron spectroscopy techniques, in this work we reveal the presence of a kinetic barrier to the reduction of LLZO by Li metal, with the extent of oxygen-deficient interphase (ODI) formation depending sensitively on the energetics of Li metal arriving at the Li vertical bar LLZO interface. Despite the clear presence of a kinetic barrier to reduction, the electrochemical response of the Li vertical bar LLZO interface is unchanged by the presence of the ODI, indicating that ODI formation during electrochemical cycling does not hinder charge transfer across the Li vertical bar LLZO interface. Overall, these results reveal that the reactivity of LLZO with Li metal depends not only on the material properties of the adjoining phases (i.e., surface purity and active contact) and their resulting thermodynamic stability but also on the energy input at the interface and the resulting reaction kinetics. Furthermore, the presence of a kinetic barrier to reduction highlights the additional complexities governing the reactivity of solid-state interfaces in ASSBs and underscores the importance of operando characterization of interfacial stability to design more robust, high-performance protection strategies for solid electrolytes in contact with reactive electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Electronic Structure of Excitons in Hematite Fe 2 O 3

We present a computational study of the electronic structure of excitons in hematite Fe 2 O 3 using density functional theory (DFT) and time-dependent DFT theory. Excitons resulting from solar light absorption are precursors to photogenerated charge carriers, which are the active species in solar-to-fuel redox conversions. Upon vertical excitation, the lowest energy exciton is found to have a multideterminantal wave function, describing an electron excited from O 2p states to Fe 3d states of the nearest basal plane, possibly supporting a recent XUV-derived exciton structure with a radius of a single Fe–O bond. Subsequently, the exciton self-traps into stable electron–hole pair structures whose wave functions are well-approximated by single excitations conveniently labeled “HOMO–LUMO”. Self-trapping yields electron–hole pair structures whereby the electron state (electron polaron-like) is separated from the hole state (hole polaron-like) by 3, 5, 7, 9, ... basal planes in structures referred to as Exc-3, Exc-5, Exc-7, Exc-9, ... The natural transition orbitals (electron–hole states) exhibit a strong “localized” character. The hole state is best described as a (FeO 6 ) + octahedral moiety carrying ~0.65 h + charge, and overall ~70% of the hole is assigned to O 2p atomic states, mostly on the four (4) equatorial O atoms of the moiety. The electron state is best described as a (FeFe) – moiety carrying ~0.80 e – charge in Fe 3d atomic states. The lattice distortions around the hole site exhibit Fe–O bond shortening due to the removal of electron density from an orbital state with O–Fe–O antibonding character. The electron site shows Fe–Fe distance shortening due to excess electron density added to an orbital state with in-phase Fe–Fe interactions. Furthermore, these excitonic structures with increasing electron–hole separation can be viewed as the onset of exciton separation into isolated electrons and hole polarons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Revisions of ORNL 188 W Process Based on Nonradiological Experiments

Tungsten-188 is in widespread use in the 188 W(t 1/2 = 69 d)/ 188 Re(t 1/2 =16.9 h) biomedical generator. Oak Ridge National Laboratory (ORNL) has been providing this product to the world since 1999. At ORNL, 188 W is produced via irradiation in ORNL’s High Flux Isotope Reactor (HFIR). Enriched 186 W targets in the form of sintered metallic pellets or rings achieve a compact loading in the irradiation vessel, providing a high yield per unit target. The enrichment of the target is >90% 186 W and this isotope undergoes double neutron capture to produce the desired 188W product. While 188 W is produced by neutron bombardment, 191 Os(t 1/2 = 15.4 d) is simultaneously produced as a byproduct requiring separation from 188 W by post irradiation treatment. In the current processing pathway, the irradiated W metal pellets or rings are first converted into an oxide form of WO 3 by heating the irradiated W metal target at 750°C in a quartz reaction vessel inside a vertical furnace under a constant flow of air. During heating, W metal reacts with oxygen in the air to produce WO 3 , which is soluble in 6 M NaOH for later purification process. This oxidation process also converts 188 Os (the decay daughter of 188 W) and 191 Os (the irradiation produced byproduct) into OsO 4 , a highly volatile and toxic gas. The gaseous effluents driven from the quartz reaction vessel are passed through a scrubbing array to remove OsO 4 before the air is discharged from the process. This heterogeneous oxidation method simultaneously achieves two goals: 1). converting metal target to a soluble oxide form and 2). separating volatile OsO 4 away from the solid WO 3 product, although harmful 191 OsO 4 is unfavorable but being taken care of by the down road scrubbing array. In the past twenty years the existing OsO 4 scrubbing array served well in preventing OsO 4 from being released into the environment, until July of 2020 when a minor amount of 191 Os was found to have been released into the environment which resulted in a standby of building 4501 from July 2020 to February 2021.Since October 20 of 2020 a team of researchers from groups of RSTD and NEFD were organized to work for the following tasks: (1) Fully understand the root causes of the July incident of 191 Os release and create a redesigned 191 Os scrubbing array to eliminate the possibility of 191 Os release in future 188 W processes; (2) In Phase-1 tests, select a correct air (O 2 ) flow rate in 188 W process, that provides sufficient oxygen for the oxidation of tungsten metal, while still allowing sufficient resonance time for 191 Os removal by the scrubbers; (3) In Phase-2 tests, confirm the efficiency of the redesigned scrubbing array to absorb excessive amount of non-rad OsO 4 (generated from Os metal powder) at the selected air flow; (4) In Phase-3 tests, simulate the hot cell 188 W process by heating pressed non-rad W-Os pellets in a quartz vessel of new design and the new array under operation parameters selected in Phase-1 and -2; (5) Complete a new operation procedure for 188 W process in hot cell and a TM report as a summary of the Os Mitigation Project. This TM report summarizes the investigation on root causes from technical aspects in July 191 Os release and the correspondent improvement towards the redesigned 191 Os scrubbing system. The report presents the efforts made on selection of optimal operational parameter with the new scrubbing array and the quartz vessel of new design, based on experiment data obtained in Phase-1, -2 and -3 tests. Further renovations for 188 W process and additional improvement of W target treatment will also be discussed.

07 ISOTOPE AND RADIATION SOURCES↗