Engineering Papers⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

MPACT Theory Manual (V.4.3)

This manual presents the theory underlying the three-dimensional (3D) whole-core, pin-resolved neutron transport calculation methodologies employed in the MPACT code. MPACT’s primary goal is to provide accurate sub-pin power distributions in a computationally efficient manner. To accomplish this, MPACT offers several different transport method options. The 2D/1D method, in which 3D problems are decomposed into an axial stack of radial “slices,” is currently the most commonly used. In this option, two-dimensional (2D) planar solutions are provided by the method of characteristics (MOC), and axial solutions are provided via one-dimensional (1D) approximate diffusion, or P3 solutions. The radial and axial solutions are coupled by (i) axial and radial transverse leakages, and (ii) a global 3D coarse mesh finite difference (CMFD) solve, which provides both acceleration and stability to the solution iteration scheme. The subsequent chapters of this manual present a range of topics, including MOC, CMFD, axial nodal transport solvers, 2D/1D, self-shielding, depletion, thermal-hydraulics, and transient methods. Development of the underlying theory of the 2D/1D method is presented, diagrams are included to highlight important algorithmic flow, and important concepts are discussed as appropriate. This manual is intended to be self-sufficient, but references to published articles and other materials are included for further reading. MPACT is a relatively new code, with new capabilities and many computational methods that did not exist until recently. This manual is an even newer document, with chapters written by several different code contributors working under time constraints. For these reasons, the manual is not yet complete and finalized. Future versions of this manual will address current deficiencies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

In Situ High-Temperature Ultrafast Electron Diffraction through Integrated Furnace and MEMS Platforms

Temperature fundamentally governs phase stability, defect evolution, and transport behavior in materials. Despite its central role, direct measurements of structural evolution at elevated temperatures on ultrafast timescales have remained limited. Here, we report the design, integration, and validation of 2 complementary in situ heating platforms that substantially extend the thermal operating range of ultrafast electron diffraction (UED). A compact furnace-type heating stage enables stable diffraction measurements from room temperature to 800 K with ±0.1 K stability under ultrahigh vacuum, achieved through multi-sensor feedback control, dual air-cooling channels, and a thermally isolated motion stage. In parallel, a microelectromechanical system (MEMS)-based heating platform provides rapid thermal response and access to extreme temperatures ≥1,373 K with ±0.1 K stability over hundreds-micrometer regions while supporting simultaneous electrical biasing for electrothermal coupling studies. Absolute temperature calibration is established using diffraction-based thermometry via aluminum lattice expansion and independently validated through in situ melting of bismuth thin films. UED measurements further reveal pronounced temperature-dependent nonequilibrium lattice dynamics in bismuth, including modifications to electron–phonon coupling and Debye–Waller behavior, as well as enhanced ultrafast diffuse scattering in aluminum at elevated temperatures. Together, these developments establish a practical framework for quantitative, time-resolved studies of temperature-driven kinetics and nonequilibrium structural dynamics under extreme thermal environments.

Bai, Qianqian [Chinese Academy of Sciences (CAS), ↗
DOE OSTI2022

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Superionic fluoride gate dielectrics with low diffusion barrier for two-dimensional electronics

Exploration of new dielectrics with a large capacitive coupling is an essential topic in modern electronics when conventional dielectrics suffer from the leakage issue near the breakdown limit. Here, to address this looming challenge, we demonstrate that rare-earth metal fluorides with extremely low ion migration barriers can generally exhibit an excellent capacitive coupling over 20 μF cm -2 (with an equivalent oxide thickness of ~0.15 nm and a large effective dielectric constant near 30) and great compatibility with scalable device manufacturing processes. Such a static dielectric capability of superionic fluorides is exemplified by MoS 2 transistors exhibiting high on/off current ratios over 10 8 , ultralow subthreshold swing of 65 mV dec -1 and ultralow leakage current density of ~10 -6 A cm -2 . Therefore, the fluoride-gated logic inverters can achieve notably higher static voltage gain values (surpassing ~167) compared with a conventional dielectric. Furthermore, the application of fluoride gating enables the demonstration of NAND, NOR, AND and OR logic circuits with low static energy consumption. In particular, the superconductor–insulator transition at the clean-limit Bi 2 Sr 2 CaCu 2 O 8+δ can also be realized through fluoride gating. In conclusion, our findings highlight fluoride dielectrics as a pioneering platform for advanced electronic applications and for tailoring emergent electronic states in condensed matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Se diffusion in CdTe thin films for photovoltaics

Manipulating CdSeTe bandgrading to enhance photocurrent and carrier lifetime is an essential part of high-performance CdTe photovoltaics (PVs). Here, we examine Se diffusion kinetics in single-crystal and polycrystalline CdTe during deposition, thermal annealing, and CdCl 2 treatments. Se distributions are determined by dynamic secondary-ion-mass spectroscopy and Auger electron spectroscopy depth profiling and coupled with electron backscatter diffraction images of the crystalline structure. Effective bulk and grain boundary diffusion coefficients are determined by analytical models and discussed in the context of processing and film morphology. Se is found to diffuse in CdTe at much higher rates during CdCl 2 treatments than with thermal processing alone. GB diffusion also occurs at a significantly faster rate than bulk diffusion. As a result of these two effects, the near interface bulk and GB Se diffusion during CdCl 2 treatments dominates the bandgrading profiles in CdTe PVs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials

Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.

14 SOLAR ENERGY↗
DOE OSTI2026

Structured electrolytes facilitate Grotthuss-type transport for enhanced proton-coupled electron transfer reactions

Concentrated hydrogen-bonded electrolytes (CoHBEs) are structured, electrochemically stable, less-volatile alternatives to aqueous and dilute nonaqueous electrolytes, however, with high viscosities that limit molecular diffusion. This work provides an understanding of the proton conduction mechanism in CoHBEs based on mixtures of acids and azoles and establishes a link between the structurally dictated transport properties and the proton-coupled electron transfer (PCET) reaction rates that can be leveraged for enhancing electrochemical reactions. Diffusion and relaxation NMR studies suggest a breaking of the viscosity–conductivity tradeoff, where at high azole concentrations (>45 mol%), Grotthuss transport is more likely with lowered proton transfer energy barriers between the azole and the acid according to the machine learning (ML) accelerated ab initio path integral MD (AI-PIMD) simulations. Proton conduction pathways are found to be switchable between the hydrogen bonding networks of the acid and the azole, with imidazole chain forming structures better facilitating Grotthuss hopping. Supported by small-angle neutron scattering studies, the chains are found to have six member molecules on average with maximum of 3 to 4 imidazole/imidazoliums at 50 to 60 mol%. Despite their high viscosities, the measured PCET rates for quinones and phenazines measured in the protic CoHBEs present relatively high electron transfer rate constants (k 0 ~ 10 −4 cm/s), validated by rotating disc electrode and scanning electrochemical microscopy measurements. The results demonstrate that strategic tuning of hydrogen-bond donor–acceptor interactions enables the decoupling of proton transport and viscosity, thereby impacting PCET reactions.

electrokinetics↗
DOE OSTI2021

Efficient diffusion of superdense lithium via atomic channels for dendrite-free lithium–metal batteries

The non-uniform aggregation of fast-diffused Li on an anode surface would aggravate its tip-effect-induced nucleation/growth, leading to the notorious dendrite growth in Li metal batteries (LMBs). Tuning the Li diffusion on the anode surface has been regarded previously as a mainstream method to induce its uniform deposition, while the diffusion of Li in the anode bulk is usually ignored. In this study, conceptually different from the classic surface modification, we propose a molecular tunnelling strategy to construct atomic channels in graphite bulk, which enables the fast diffusion of superdense Li. Density functional theory calculations and ab initio molecular dynamics simulations prove that the bulk diffusion through atomic channels could become a new and dominating path. Its reversible and efficient diffusion has been further visualized by in situ transmission electron microscopy. As a result, when coupled with high-loading LiFePO 4 cathodes (20 mg cm –2 ), a high areal capacity and 100% capacity retention are achieved over 370 cycles. Through this work a new strategy is developed based on the bulk-diffusion of superdense Li for dendrite-free LMBs, which can be pervasive in other high-performance energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2024

All-Perovskite Multicomponent Nanocrystal Superlattices

Nanocrystal superlattices (NC SLs) have long been sought as promising metamaterials, with nanoscale-engineered properties arising from collective and synergistic effects among the constituent building blocks. Lead halide perovskite (LHP) NCs come across as outstanding candidates for SL design, as they demonstrate collective light emission, known as superfluorescence, in single- and multicomponent SLs. Thus far, LHP NCs have only been assembled in single-component SLs or coassembled with dielectric NC building blocks acting solely as spacers between luminescent NCs. Here, we report the formation of multicomponent LHP NC-only SLs, i.e., using only CsPbBr 3 NCs of different sizes as building blocks. The structural diversity of the obtained SLs encompasses the ABO 6 , ABO 3 , and NaCl structure types, all of which contain orientationally and positionally locked NCs. For the selected model system, the ABO 6 -type SL, we observed efficient NC coupling and Förster-like energy transfer from strongly confined 5.3 nm CsPbBr 3 NCs to weakly confined 17.6 nm CsPbBr 3 NCs, along with characteristic superfluorescence features at cryogenic temperatures. Spatiotemporal exciton dynamics measurements reveal that binary SLs exhibit enhanced exciton diffusivity compared to single-component NC assemblies across the entire temperature range (from 5 to 298 K). The observed coherent and incoherent NC coupling and controllable excitonic transport within the solid NC SLs hold promise for applications in quantum optoelectronic devices.

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

Switching CO 2 Electroreduction Selectivity Between C 1 and C 2 Hydrocarbons on Cu Gas-Diffusion Electrodes

We report that regulating the selectivity toward a target hydrocarbon product is still the focus of CO 2 electroreduction. Here, we discover that the original surface Cu species in Cu gas-diffusion electrodes plays a more important role than the surface roughness, local pH, and facet in governing the selectivity toward C 1 or C 2 hydrocarbons. The selectivity toward C 2 H 4 progressively increases, while CH 4 decreases steadily upon lowering the Cu oxidation species fraction. At a relatively low electrodeposition voltage of 1.5 V, the Cu gas-diffusion electrode with the highest Cu δ+ /Cu 0 ratio favors the pathways of hydrogenation to form CH 4 with maximum Faradaic efficiency of 65.4% and partial current density of 228 mA cm -2 at -0.83 V vs RHE. At 2.0 V, the Cu gas-diffusion electrode with the lowest Cu δ+ /Cu 0 ratio prefers C–C coupling to form C 2+ products with Faradaic efficiency topping 80.1% at -0.75 V vs RHE, where the Faradaic efficiency of C 2 H 4 accounts for 46.4% and the partial current density of C 2 H 4 achieves 279 mA cm -2 . This work demonstrates that the selectivity from CH 4 to C 2 H 4 is switchable by tuning surface Cu species composition of Cu gas-diffusion electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Multiplicative noise and the diffusion of conserved densities

Stochastic fluid dynamics governs the long time tails of hydrodynamic correlation functions, and the critical slowing down of relaxation phenomena in the vicinity of a critical point in the phase diagram. In this work we study the role of multiplicative noise in stochastic fluid dynamics. Multiplicative noise arises from the dependence of transport coefficients, such as the diffusion constants for charge and momentum, on fluctuating hydrodynamic variables. We study long time tails and relaxation in the diffusion of a conserved density (model B), and a conserved density coupled to the transverse momentum density (model H). Careful attention is paid to fluctuation-dissipation relations. We observe that multiplicative noise contributes at the same order as non-linear interactions in model B, but is a higher order correction to the relaxation of a scalar density and the tail of the stress tensor correlation function in model H.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spin-caloric resistance of Dirac plasma in a graphene Corbino device

The thermal resistance of a spin-polarized hydrodynamic Dirac plasma in graphene is considered. Here, a mechanism for the coupling of heat and spin flows is discussed, demonstrating that spin diffusion and spin thermocurrent modify viscous dissipation, leading to a significant enhancement of thermal resistance. Practical calculations are then presented for graphene devices in the Corbino geometry.

Graphene↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗