SEARCH · Engineering Papers

Results for “diffusion couples”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗

DOE OSTI2023

High methane concentrations in tidal salt marsh soils: Where does the methane go?

Abstract Tidal salt marshes produce and emit CH 4 . Therefore, it is critical to understand the biogeochemical controls that regulate CH 4 spatial and temporal dynamics in wetlands. The prevailing paradigm assumes that acetoclastic methanogenesis is the dominant pathway for CH 4 production, and higher salinity concentrations inhibit CH 4 production in salt marshes. Recent evidence shows that CH 4 is produced within salt marshes via methylotrophic methanogenesis, a process not inhibited by sulfate reduction. To further explore this conundrum, we performed measurements of soil–atmosphere CH 4 and CO 2 fluxes coupled with depth profiles of soil CH 4 and CO 2 pore water gas concentrations, stable and radioisotopes, pore water chemistry, and microbial community composition to assess CH 4 production and fate within a temperate tidal salt marsh. We found unexpectedly high CH 4 concentrations up to 145,000 μmol mol −1 positively correlated with S 2− (salinity range: 6.6–14.5 ppt). Despite large CH 4 production within the soil, soil–atmosphere CH 4 fluxes were low but with higher emissions and extreme variability during plant senescence (84.3 ± 684.4 nmol m −2 s −1 ). CH 4 and CO 2 within the soil pore water were produced from young carbon, with most Δ 14 C‐CH 4 and Δ 14 C‐CO 2 values at or above modern. We found evidence that CH 4 within soils was produced by methylotrophic and hydrogenotrophic methanogenesis. Several pathways exist after CH 4 is produced, including diffusion into the atmosphere, CH 4 oxidation, and lateral export to adjacent tidal creeks; the latter being the most likely dominant flux. Our findings demonstrate that CH 4 production and fluxes are biogeochemically heterogeneous, with multiple processes and pathways that can co‐occur and vary in importance over the year. This study highlights the potential for high CH 4 production, the need to understand the underlying biogeochemical controls, and the challenges of evaluating CH 4 budgets and blue carbon in salt marshes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Defects of perovskite semiconductor CsPbBr3 investigated via photoluminescence and thermally stimulated current spectroscopies

Halide perovskites are essential materials for hard radiation detectors at ambient temperature. To improve detector performance, charge transport must be investigated and optimized. Using photoluminescence (PL) and thermally stimulated current (TSC) spectroscopies, we investigate photogenerated charge carriers in Bridgman-grown CsPbBr3 single crystals to understand the nature of charge transport. PL spectroscopy of these halide perovskites revealed the presence of strong emission bands at the band edge, which were attributed to free or bound excitons. It is shown that a wide broadening of the excitonic linewidth in these halide perovskites arises from strong exciton–phonon coupling, which is substantially dominated by longitudinal optical phonons via Fröhlich interaction. An additional contribution due to the presence of ionized impurities was also observed. Crystals with a detectable sensitivity to high-energy gamma radiation are characterized by a higher intensity and a narrower linewidth of the principal PL peak at 2.326 eV. Defect states beyond 2.214 eV have a negative impact on detector sensitivity to high-energy gamma radiation. TSC spectroscopy reveals an array of trap levels spanning 0.15–0.70 eV, attributed to intrinsic point defects and multiple extrinsic defects involving dopants or impurities. Defects identified included Cs and Br vacancies, as well as Pb interstitials with concentrations in the 1011–1016 cm−3 range. Understanding how the synthesis process impacts the types and concentrations of the defects present is currently under investigation. Elimination or suppression of the defect/trap states should result in halide perovskite materials with longer carrier diffusion lengths and improved detector characteristics.

Physics↗

DOE OSTI2026

Modeling glass degradation and release of radionuclides from vitrified waste for performance assessment simulations

The release of radionuclides initially encapsulated in a slowly degrading solid waste form and contained in an eventually corroding canister defines the source term for numerical simulations for the assessment of a geologic repository for high-level radioactive waste. While the details of waste degradation, canister corrosion, and dissolution and mobilization of the radionuclides in pore water include complex chemical reaction and transport processes that are coupled to the thermal, hydrological, microbiological, and mechanical conditions in the repository, the source-term model suitable for use in a numerical performance assessment model should be a defensible abstraction of these mechanisms. We developed a radiological source-term model and implemented it into a non-isothermal flow and transport simulator. While the proposed source-term model is applicable to various waste forms, canister systems, and disposal concepts, we specifically considered radionuclide releases from vitrified high-level waste placed in a cylindrical canister disposed in a deep vertical borehole repository. In this model, waste degradation is a function of temperature, and it can be adjusted to evaluate the influence of and propagate uncertainties in pH, passivation reactions, and chemical conditions as well as geometrical factors. The time-dependent, congruent release of safety-relevant radionuclides present in the decaying inventory is then calculated. Finally, the radionuclides are mobilized by diffusive and advective transport according to the thermo-hydraulic conditions prevailing in the near field of the repository, from where they migrate through the geosphere to the accessible environment. We examine the influence of the source-term model’s parameters on performance assessment calculations through sensitivity and uncertainty propagation analyses, identifying influential factors and confirming the upper bound of their impact. These considerations align with the overarching goal of repository design, which is to demonstrate that engineered and natural barriers can collectively delay radionuclide migration for timescales far exceeding human planning, thereby providing multiple, redundant barriers against environmental contamination.

iTOUGH2↗

DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Drought-induced changes in groundwater-surface water exchange at Lake Mead area

This study focuses on the Lake Mead region in the southwestern United States, a key water reservoir serving over 25 million people and agricultural lands across several states. The area has experienced recurring anthropogenic droughts since the early 2000s. We investigate the hydrological response of the coupled surface-groundwater system to the 2020–2022 drought, one of the most severe on record. To this end, we use Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data to quantify vertical land motion caused by the elastic response of the crust to water-mass loss in the lake vicinity. Next, we apply an inverse elastic load modeling framework to quantify water loss. We further assess possible hydraulic connectivity between Lake Mead and adjacent groundwater reservoirs. We detect ground uplift of up to 8 mm/yr near the lake center, likely due to crustal rebound from reduced water-mass loading. We estimated the total water storage loss at 3.03 ± 0.25 km 3 /yr across a 3150 km 2 area surrounding Lake Mead. Groundwater accounts for approximately a third of that, being 0.94 ± 0.32 km 3 /yr. In addition, observed time lags of 6–98 days between lake and groundwater level responses, corresponding to a lateral diffusivity of 3.2–86 m 2 /s, suggest spatially variable connectivity between the lake and aquifers. These findings highlight that drought impacts propagate through the subsurface within interconnected systems, resulting in reduced buffering capacity of groundwater resources following droughts, and emphasizing the need for more integrated surface and groundwater management strategies to enhance resilience under climate and anthropogenic stressors.

58 GEOSCIENCES↗

DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Dynamic entity formed by protein and its hydration water

The interaction between protein and water plays a pivotal role in shaping the structure, dynamics, and function of biomacromolecules. A comprehensive understanding of this intricate interplay necessitates a systematic evaluation of interaction strength and its consequential impact on the dynamics of proteins and water across diverse protein systems. Despite numerous works on understanding the dynamics of water and proteins and the coupling between them, there are still unanswered questions. Here, we combine neutron scattering and isotope labeling to probe the dynamics of proteins and their hydration water in a variety of protein systems. We consider proteins of different structures and varying thermostability as well as proteins within living cells with distinct growth temperatures. Simultaneous characterization of protein and hydration water dynamics across diverse systems was achieved. Moreover, we performed water sorption isothermal measurements on three representative proteins to correlate the observed dynamics with the strength of the interaction energies governing each system. The experimental results underscore that proteins manifesting stronger attractive interactions with water display diffusionlike dynamics with higher flexibility upon hydration, concomitant with a reduced mobility in hydration water. Significantly, our findings suggest that, in fact, it is the interaction between protein and its hydration water that facilitates the transfer of mobility from water to protein, with stronger interactions correlating to greater protein flexibility and slower hydration water diffusion. Published by the American Physical Society 2024

Ye, Yongfeng (ORCID:0009000792374198)↗

DOE OSTI2020

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗

DOE OSTI2023

Field reversal in low pressure, unmagnetized radio frequency capacitively coupled argon plasma discharges

In general, the radio frequency (rf) electric field within a sheath points toward the metal electrode in low pressure, unmagnetized rf electropositive capacitively coupled plasma (CCP) glow discharges. This is due to the large ratio of electron to ion mobility and the formation of an ion sheath. In this work, we studied, using fully kinetic particle-in-cell simulations, a reversed electric field induced by the strong secondary electron emission during the phase of sheath collapse in a high-voltage rf-driven low pressure CCP glow discharge. We explored the transition behavior of the formation of field reversal as a function of driving voltage amplitude and found that field reversal starts to form at around 750 V, for a discharge with an electrode spacing of 4 cm at 10 mTorr argon pressure driven at 13.56 MHz. Accordingly, the energy distribution function of electrons incident on the electrode shows peaks from around 3 to 10 eV while varying the driving voltage from 150 to 2000 V, showing potentially beneficial effects in plasma material processing where relatively directional electrons are preferred to solely thermal diffusion electrons.

Physics↗

DOE OSTI2022

Measurements of atoms and metastable species in N 2 and H 2 –N 2 nanosecond pulse plasmas

We report time-resolved, absolute number densities of metastable N 2 (A 3 Σ u + , v = 0, 1) molecules, ground state N 2 and H atoms, and rotational–translational temperature have been measured by tunable diode laser absorption spectroscopy and two-photon absorption laser-induced fluorescence in diffuse N 2 and N 2 –H 2 plasmas during and after a nanosecond pulse discharge burst. Comparison of the measurement results with the kinetic modeling predictions, specifically the significant reduction of the N 2 (A 3 Σ u + ) populations and the rate of N atom generation during the burst, suggests that these two trends are related. The slow N atom decay in the afterglow, on a time scale longer than the discharge burst, demonstrates that the latter trend is not affected by N atom recombination, diffusion to the walls, or convection with the flow. This leads to the conclusion that the energy pooling in collisions of N 2 (A 3 Σ u + ) molecules is a major channel of N 2 dissociation in electric discharges where a significant fraction of the input energy goes to electronic excitation of N 2 . Additional measurements in a 1% H 2 –N 2 mixture demonstrate a further significant reduction of N 2 (A 3 Σ u + , v = 0, 1) populations, due to the rapid quenching by H atoms accumulating in the plasma. Comparison with the modeling predictions suggests that the N 2 (A 3 Σ u + ) molecules may be initially formed in the highly vibrationally excited states. The reduction of the N 2 (A 3 Σ u + ) number density also diminishes the contribution of the energy pooling process into N2 dissociation, thus reducing the N atom number density. The rate of N atom generation during the burst also decreases, due to its strong coupling to N 2 (A 3 Σ u + , v) populations. On the other hand, the rate of H atom generation, produced predominantly by the dissociative quenching of the excited electronic states of N 2 by H 2 , remains about the same during the burst, resulting in a nearly linear rise in the H atom number density. Comparison of the kinetic model predictions with the experimental results suggests that the yield of H atoms during the quenching of the excited electronic state of N 2 by molecular H 2 is significantly less than 100%. The present results quantify the yield of N and H atoms in high-pressure H 2 –N 2 plasmas, which have significant potential for ammonia generation using plasma-assisted catalysis

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Pervasive symmetry-lowering nanoscale structural fluctuations in the cuprate La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4

The cuprate superconductors are among the most widely studied quantum materials, yet there remain fundamental open questions regarding their electronic properties and the role of structural degrees of freedom. Recent neutron and x-ray scattering measurements uncovered exponential scaling with temperature of the strength of orthorhombic fluctuations in the tetragonal phase of La 2−𝑥⁢ Sr 𝑥 ⁢Cu⁢O 4 and Tl 2 ⁢Ba 2 ⁢Cu⁢O 6+𝛿 , unusual behavior that closely resembles prior results for the emergence of superconducting fluctuations, and that points to a common origin rooted in inherent correlated structural inhomogeneity. Here, in this study, we use neutron and x-ray diffuse scattering to provide further insight into structural fluctuations in La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4 . We perform measurements up to temperatures that approach the melting point of the undoped parent compound, and we investigate the effects of in situ in-plane uniaxial elastic stress and plastic deformation. Our neutron scattering results for the parent compound La 2 ⁢Cu⁢O 4 reveal that short-range orthorhombic fluctuations persist to the maximum experimental temperature of nearly 1000 K, i.e., to a significant fraction of the crystallization temperature. At this temperature, the spatial characteristic length extracted from the momentum-space data is still about three lattice constants. The neutron scattering experiment enables quasistatic discrimination and reveals that the response is increasingly dynamic at higher temperatures. We also find that neither compressive elastic stress along the tetragonal [110] direction nor irreversible plastic deformation activating the ⟨100⟩⁢[001] slip system significantly alter the robust orthorhombic fluctuations. Overall, these results support the notion that these fluctuations couple to subtle, underlying inhomogeneity that underpins the cuprate phase diagram. Finally, in La 1.8⁢ Sr 0.2⁢ Cu⁢O 4 , we uncover low-energy structural fluctuations at a nominally forbidden reflection that are distinct from the orthorhombic distortions. While the origin of these fluctuations is not clear, they might be related to the presence of extended defects such as dislocations or stacking faults.

Spieker, R. J. [Univ. of Minnesota, Minneapolis, M↗

DOE OSTI2021

An unstructured C-grid type variational formulation for the sea ice dynamics

Historically, B-grid formulations of sea ice dynamics have been dominant because they have matched the grid type used by ocean models. The reason for the grid match is simple – it facilitates penetration of the curl of ice-ocean stress into the deep ocean with minimal numerical diffusivity because sea ice and ocean velocity are co-located. In recent years, as ocean models have increasingly progressed to C-grids, sea ice models have followed suit on quadrilateral meshes, but the implementation of an unstructured C-grid sea ice models is new. We present an unstructured C-grid discretization of the Elastic Viscous Plastic (EVP) rheology, where the velocity unknowns are discretized at the edges of the mesh cells with n-sides, where typically n is greater than or equal to four, rather than at the vertices, as in the B-grid. Our framework of choice is the Model for Prediction Across Scales (MPAS) within E3SM, the climate model of the U.S. Department of Energy, although our approach is general and could be applied to other models as well. While MPAS-Seaice is currently defined on a B-grid, MPAS-Ocean runs on a C-grid, hence interpolation operators are heavily used when coupled simulations are performed. In this work, we describe a mathematical formulation to transition the dynamics of MPAS-Seaice to a C-grid, in order to ultimately facilitate the coupling with MPAS-Ocean and reduce numerical errors associated with this communication. Numerical results are reported to highlight the features of the method.

58 GEOSCIENCES↗

DOE OSTI2025

Coherent Control over Nuclear Hyperpolarization Using an Optically Initializable Chromophore-Radical System

Chromophore radicals (CR) are emerging as important components for molecular quantum information science (QIS), especially in the context of quantum sensing. Here, we demonstrate that the optically hyperpolarized electrons in a 1,6,7,12-tetrakis(4-tert-butylphenoxy)-perylene-3,4,9,10-bis(dicarboximide) (tpPDI) covalently linked to a partially deuterated 1,3-bis(diphenylene)-d 16 -2-phenylallyl radical (BDPA-d 16 ) can be coherently manipulated via pulsed dynamic nuclear polarization (DNP) methods to transfer polarization to nuclear spins and back. Under light illumination at 85 K, electron hyperpolarization in BDPA is enhanced 2.1- to 2.4-fold over thermal polarization and lasts for more than 100 ms. By applying nuclear orientation via electron spin-locking (NOVEL) DNP, this optically amplified electron hyperpolarization was successfully transferred to a 1 H nuclear spin within the CR system and efficiently returned to the electron spin for readout via reverse-NOVEL. The NOVEL transfer efficiency of 65% amounts to a 688-fold nuclear spin hyperpolarization of the target nuclear spin, considering the 2.1-fold electron spin hyperpolarization. This reversible coherent manipulation of hyperpolarization transfer highlights the utility of CR systems to initialize and read out nuclear spin states in a disordered matrix at moderate cryogenic temperatures. Coupled with CRs’ environmental compatibility, tunability, and precise state initialization, these results highlight the promising role of nuclear spins in CRs for QIS applications, including quantum sensing and memory.

charge transfer↗

DOE OSTI2022

Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE

Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗