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 523 records · Page 29

Coupled Process/Device Modeling and Point Defect Engineering of Cu(In,Ga)Se 2 Solar Cells

Point defects directly impact solar cell device performance by limiting the carrier lifetime. In this work, density functional theory calculations are first used to determine the formation energy and diffusion energy barriers of dominant defects in Cu(In,Ga)Se 2 . Next, continuum reaction–diffusion models are developed to analyze the redistribution of defects during manufacturing processes. We estimate defect capture cross sections using a first-principles-based approach. Here, these cross sections are combined with our calculated defect profiles and trap energy levels to parameterize a Shockley–Read–Hall recombination model, which we implement into a device simulator to predict carrier lifetimes and device performance. In that way, a predictive technology computer aided design model is built to predict and optimize the performance of Cu(In,Ga)Se 2 solar cells.

14 SOLAR ENERGY↗

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗

In Situ Identification of Reaction Intermediates and Mechanistic Understandings of Methane Oxidation over Hematite: A Combined Experimental and Theoretical Study

Effective methane utilization for either clean power generation or value-added chemical production has been a subject of growing attention worldwide for decades, yet challenges persist mostly in relation to methane activation under mild conditions. Here, we report hematite, an earth-abundant material, to be highly effective and thermally stable to catalyze methane combustion at low temperatures (<500 °C) with a low light-off temperature of 230 °C and 100% selectivity to CO 2 . The reported performance is impressive and comparable to those of precious-metal-based catalysts, with a low apparent activation energy of 17.60 kcal·mol –1 . Our theoretical analysis shows that the excellent performance stems from a tetra-iron center with an antiferromagnetically coupled iron dimer on the hematite (110) surface, analogous to that of the methanotroph enzyme methane monooxygenase that activates methane at ambient conditions in nature. Isotopic oxygen tracer experiments support a Mars van Krevelen redox mechanism where CH 4 is activated by reaction with a hematite surface oxygen first, followed by a catalytic cycle through a molecular-dioxygen-assisted pathway. Surface studies with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and density functional theory (DFT) calculations reveal the evolution of reaction intermediates from a methoxy CH 3 –O–Fe, to a bridging bidentate formate b-HCOO–Fe, to a monodentate formate m-HCOO–Fe, before CO 2 is eventually formed via a combination of thermal hydrogen-atom transfer (HAT) and proton-coupled electron transfer (PCET) processes. Finally, the elucidation of the reaction mechanism and the intermediate evolutionary profile may allow future development of catalytic syntheses of oxygenated products from CH 4 in gas-phase heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proximate spin liquid and fractionalization in the triangular antiferromagnet KYbSe 2

The Heisenberg triangular-lattice quantum spin liquid and its phase transitions to nearby magnetic orders have received much theoretical attention, but clear experimental manifestations of these states are rare. Here we demonstrate that a spin-half delafossite material, namely, KYbSe 2 , shows close proximity to the triangular-lattice Heisenberg quantum spin liquid. Using neutron scattering, we identify a diffuse continuum with a sharp lower bound within the measured spectra. Applying entanglement witnesses to the data indicates multipartite entanglement spread between its neighbours, and an analysis of its magnetic-exchange couplings reveals close proximity to the theoretical quantum spin-liquid phase. The key features of the data are reproduced by Schwinger boson theory and tensor network calculations with a substantial next-nearest-neighbour coupling. The strength of the dynamical structure factor at the Brillouin-zone K point shows a scaling collapse down to 0.3 K, indicating the existence of a second-order quantum phase transition. Finally, comparing this with previous theoretical work suggests that the proximate phase at a larger next-nearest-neighbour coupling is a gapped $\mathbb{Z}_{2}$ spin liquid, resolving a long-debated issue. Similar content being viewed by others

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cavity molecular dynamics simulations of liquid water under vibrational ultrastrong coupling

We simulate vibrational strong coupling (VSC) and vibrational ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O–H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Furthermore, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H 2 O molecules especially under V-USC, which could play a role in ground-state chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics and transport in the boundary-driven dissipative Klein-Gordon chain

Motivated by experiments on chains of superconducting qubits, we consider the dynamics of a classical Klein-Gordon chain coupled to coherent driving and subject to dissipation solely at its boundaries. As the strength of the boundary driving is increased, this minimal classical model recovers the main features of the "dissipative phase transition" seen experimentally. Between the transmitting and non-transmitting regimes on either side of this transition (which support ballistic and diffusive energy transport respectively), we observe additional dynamical regimes of interest. These include a regime of superdiffusive energy transport at weaker driving strengths, together with a "dissipative breather crystal" regime at stronger driving strengths, which is characterized by emergent translation symmetry, ballistic energy transport, and coherent oscillations of a nonlinear normal mode. Here, we propose a non-local Lyapunov exponent as an experimentally measurable diagnostic of many-body chaos in this system, and more generally in open systems that are only coupled to an environment at their boundaries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗

First operation of a multi-channel Q-Pix prototype: measuring transverse electron diffusion in a gas time projection chamber

We report measurements of the transverse diffusion of electrons in P-10 gas (90% Ar, 10% CH 4 ) in a laboratory-scale time projection chamber (TPC) utilizing a novel pixelated signal capture and digitization technique known as Q-Pix. The Q-Pix method incorporates a precision switched integrating transimpedance amplifier whose output is compared to a threshold voltage. Upon reaching the threshold, a comparator sends a 'reset' signal, initiating a discharge of the integrating capacitor. The time difference between successive resets is inversely proportional to the average current at the pixel in that time interval, and the number of resets is directly proportional to the total collected charge. We developed a 16-channel Q-Pix prototype fabricated from commercial off-the-shelf components and coupled them to 16 concentric annular anode electrodes to measure the spatial extent of the electron swarm that reaches the anode after drifting through the uniform field of the TPC. The swarm is produced at a gold photocathode using pulsed UV light. The measured transverse diffusion agrees with simulations in PyBoltz across a range of operating pressures (200–1500 Torr). These results demonstrate that a Q-Pix readout can successfully reconstruct the ionization topology in a TPC.

47 OTHER INSTRUMENTATION↗

Printed circuit primary heat exchanger for power conversion in heat pipe microreactor applications

Here, in the present study, a model was developed for a recuperated air Brayton cycle for heat pipe microreactors, and the effect of the primary heat exchanger type assessed. An annular flow heat exchanger and a printed circuit heat exchanger (PCHE) were evaluated. A sub-sized diffusion bonded, PCHE test specimen was manufactured and tested with nitrogen as the working fluid to validate and improve the PCHE model. A comparison of the thermal cycle efficiency for the Brayton cycle coupled to a heat pipe microreactor using each of the two primary heat exchangers showed competitive performance for both options, with achievable cycle efficiencies of 34 % and 35 % for the annular and printed circuit heat exchangers, respectively. The impact of heat pipe length and additional reheating stages on cycle efficiency was studied. Increasing the heat pipe condenser length from 0.8 to 2.0 m resulted in an increase in cycle efficiency of 3.5 % and 3.2 % for the annular and printed circuit heat exchangers, respectively. The use of additional reheating stages showed a reduction in cycle efficiency since the heat pipe surface area limits the primary heat exchanger size and effectiveness. These results highlight the importance of the primary heat exchanger on the performance of a heat pipe microreactor. Overall, the cycle model, testing, and heat exchanger model showed the potential for the PCHE to improve the economics of a heat pipe microreactor relative to the annular flow heat exchanger, which also opens the possibility to use different heat transfer fluids such as supercritical CO 2 .

PCHE↗

Modelling and analysis of salt-convection effect on oxide reduction process for uranium oxides using smoothed particle hydrodynamics

Oxygen ion transport via molten salt continuum is a key kinetic feature of electrochemical reduction process for uranium oxides in molten salt. During the electrolysis, the evolution of oxygen ion concentration field in molten salt continuum appears from porous metallic layer to open bulk salt. Despite the important role affecting process throughput, its effect coupled with controlled salt movement was not adequately investigated. To capture the reactive oxygen ion transport phenomena, a multidimensional reduction model based on smoothed particle hydrodynamics (SPH) was combined with the advection-diffusion model. A pelletized oxide fuel bed capturing interconnected salt continuum is considered as a reactant material geometry, which allows forced convective delivery of molten salt through the reactant materials with a relatively low pressure head. Here, in this study, reduction time of a pelletized-fuel bed was quantified with respect to key design parameters of the reactant material configuration. The simulation results demonstrated the advantage of salt convection which effectively removes oxygen ion products from the reactant materials and consequently accelerates the reduction process. Moreover, the reduction-time correlation of a single pelletized fuel was developed as a function of relevant dimensionless parameters and was extended to the application for a pelletized-fuel bed. Finally, several insights and implications in the use of a pelletized oxide fuel were presented. The developed SPH-based framework enabled multi-physics analysis of governing dynamics in appreciating and screening candidate reactant material configurations in a cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stabilizing Aqueous Three-Dimensional Printed Constructs Using Chitosan-Cellulose Nanocrystal Assemblies

The assembly and binding of nanoparticles at the interfaces of aqueous two-phase systems enable the three-dimensional (3D) printing of all-aqueous naturally occurring materials. When a dispersion of cellulose nanocrystals (CNCs) in an aqueous solution of polyethylene glycol (PEG) is brought into contact with chitosan dissolved in an aqueous solution of dextran, the CNCs and chitosan diffuse to the interface between the two immiscible aqueous solutions, electrostatically interact, and form a solid, membranous layer sufficiently rapidly to 3D print tubules of one liquid in the other. Additionally, the diameter, length, spatial arrangement, and stability of the printed tubules can be broadly controlled. Adsorption and directional diffusion of ionic species across the membranous layer make heavy metal ion removal possible. The results present a platform for fabricating and developing all-aqueous compartmentalized systems where function can be independently coupled to the inherent functionality of the nanoparticles or ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Static and dynamic properties of multi-ionic plasma mixtures

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "Pseudo-Ions in Jellium" (PIJ), based on an iso-electronic assumption (iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Further, applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

Stabilizing and Improving Qubit Coherence by Engineering the Noise Spectrum of Two-Level Systems

Superconducting circuits are a leading platform for quantum computing. However, their coherence times are still limited and exhibit temporal fluctuations. Those phenomena are often attributed to the coupling between qubits and material defects that can be well described as an ensemble of two-level systems (TLSs). Among them, charge fluctuators inside amorphous oxide layers contribute to both low-frequency 1 / f charge noise and high-frequency dielectric loss, causing fast qubit dephasing and relaxation. Moreover, spectral diffusion from mutual TLS interactions varies the noise amplitude over time, fluctuating the qubit lifetime. Here, we propose to mitigate those harmful effects by engineering the relevant TLS noise spectral densities. Specifically, our protocols smooth the high-frequency noise spectrum and suppress the low-frequency noise amplitude via depolarizing and dephasing the TLSs, respectively. As a result, we predict a drastic stabilization in qubit lifetime and an increase in qubit pure dephasing time. Our detailed analysis of feasible experimental implementations shows that the improvement is not compromised by spurious coupling from the applied noise to the qubit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗