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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-Abelian electric field correlator at NLO for dark matter relic abundance and quarkonium transport

We perform a complete next-to-leading order calculation of the non-Abelian electric field correlator in a SU(N c ) plasma, which encodes properties of the plasma relevant for heavy particle bound state formation and dissociation, and is different from the correlator for the heavy quark diffusion coefficient. The calculation is carried out in the real-time formalism of thermal field theory and includes both vacuum and finite temperature contributions. By working in the R ξ gauge, we explicitly show the results are gauge independent, infrared and collinear safe. The renormalization group equation of this electric field correlator is determined by that of the strong coupling constant. Our next-to-leading order calculation can be directly applied to any dipole singlet-adjoint transition of heavy particle pairs. For example, it can be used to describe dissociation and (re)generation of heavy quarkonia inside the quark-gluon plasma well below the melting temperature, as well as heavy dark matter pairs (or charged co-annihilating partners) in the early universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2020

Electrical control of spin relaxation anisotropy during drift transport in a two-dimensional electron gas

Spin relaxation was studied in a two-dimensional electron gas confined in a wide GaAs quantum well. Recently, the control of the spin relaxation anisotropy by diffusive motion was first shown in D. Iizasa et al., arXiv:2006.08253. Here, we demonstrate electrical control by drift transport in a system with two subbands occupied. The combined effect of in-plane and gate voltages was investigated using time-resolved Kerr rotation. Furthermore, the measured relaxation time presents strong anisotropy with respect to the transport direction. For an in-plane accelerating electric field along [110], the lifetime was strongly suppressed irrespective of the applied gate voltage. Remarkably, for transport along [$1\bara{1}0$], the data shows spin lifetime that was gate dependent and longer than in the [110] direction regardless of the in-plane voltage. In agreement, independent results of anisotropic spin precession frequencies are also presented. Nevertheless, the long spin lifetime, strong anisotropy and drift response seen in the data are beyond the existing models for spin drift and diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
NASA NTRS1992

Implementation/validation of a low Reynolds number two-equation turbulence model in the Proteus Navier-Stokes code: Two-dimensional/axisymmetric

The implementation and validation of the Chien low Reynolds number k-epsilon turbulence model in the two dimensional axisymmetric version Proteus, a compressible Navier-Stokes computer code, are presented. The set of k-epsilon equations are solved by marching in time using a coupled alternating direction implicit (ADI) solution procedure with generalized first or second order time differencing. To validate Proteus and the k-epsilon turbulence model, laminar and turbulent computations were done for several benchmark test cases: incompressible fully developed 2-D channel flow; fully developed axisymmetric pipe flow; boundary layer flow over a flat plate; and turbulent Sajben subsonic transonic diffuser flows. Proteus results from these test cases showed good agreement with analytical results and experimental data. Detailed comparisons of both mean flow and turbulent quantities showed that the Chien k-epsilon turbulence model given good results over a wider range of turbulent flow than the Baldwin-Lomax turbulence model in the Proteus code with no significant CPU time penalty for more complicated flow cases.

Bui, Trong T.↗
NASA NTRS1994

Impact winter and the Cretaceous/Tertiary extinctions: Results of a Chicxulub asteroid impact model

The Chicxulub impact crater in Mexico is the site of the impact purported to have caused mass extinctions at the Cretaceous/Tertiary (K/T) boundary. 2-D hydrocode modeling of the impact, coupled with studies of the impact site geology, indiate that between 0.4 and 7.0 x 10(exp 17) g of sulfur were vaporized by the impact into anhydrite target rocks. A small portion of the sulfur was released as SO3 or SO4, which converted rapidly into H2SO4 aerosol and fell as acid rain. A radiative transfer model, coupled with a model of coagulation indicates that the aerosol prolonged the initial blackout period caused by impact dust only if the aerosol contained impurities. A larger portion of sulfur was released as SO2, which converted to aerosol slowly, due to the rate-limiting oxidation of SO2. Our radiative transfer calculations, combined with rates of acid production, coagulation, and diffusion indicate that solar transmission was reduced to 10-20% of normal for a period of 8-13 yr. This reduction produced a climate forcing (cooling) of -300 W/sq.m, which far exceeded the +8 W/sq.m greenhouse warming, caused by the CO2 released through the vaporization of carbonates, and therefore produced a decade of freezing and near-freezing temperatures. Several decades of moderate warming followed the decade of severe cooling due to the long residence time of CO2. The prolonged impact winter may have been a major cause of the K/T extinctions.

Pope, Kevin O.↗
DOE OSTI2020

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Testing of a 40-kWth Counterflow Particle-Supercritical Carbon Dioxide Narrow-Channel, Fluidized Bed Heat Exchanger

Particle-based primary heat exchangers (HXs) must deliver sCO2 fluid temperatures above 700°C to couple particle-based concentrating solar receivers and thermal energy storage (TES) sub-systems with efficient sCO2 power cycles. Particle-sCO2 HX designs have struggled to meet DOE cost targets (≤ $150/kWth) due to the amount of expensive nickel alloys necessary for manufacturing full-scale, particle-sCO2 HXs. Our team has demonstrated that mild bubbling fluidization of falling particles in a counterflow narrow-channel fluidized bed can reduce required HX surface area and thus, costs by increasing particle-wall heat transfer coefficients hT,w > 800 W m-2 K-1. This paper reports on the fabrication and testing of a stainless steel, particle-sCO2 HX with 12 fluidized-bed channels approximately 10.5 mm deep spaced between diffusion-bonded, micro-channel sCO2 plates. The HX with a core length of ≈0.56 m is fed with CARBOBEAD HSP particles through a short, fluidized freeboard zone just above the core. Testing to date in the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories has shown that parallel bed fluidization maintains uniform particle inventory across the instrumented channels. Heat transfer thermal duty between the particle and sCO2 flows exceeds 30 kWth with sCO2 inlet temperatures of 200ºC and particle inlet temperatures up to 440ºC and mass flow rates of 0.2 kg s-1 fluidized by counterflowing gas flow rates of 0.005 kg s-1. Tests at higher particle and sCO2 inlet temperatures (600ºC and 400ºC respectively) are targeted to achieve > 40 kWth with model-predicted overall heat transfer coefficients U > 400 W m-2 K-1.

14 SOLAR ENERGY↗
DOE OSTI2026

Electrochemical behavior of SnCl 2 and influence of Cu and Ni ions in molten LiCl−KCl−CaCl 2 eutectic

Reliable transport and thermodynamic data for multivalent ions in complex molten salts are scarce, limiting model fidelity for electrorefining and impurity control. Here, we report a comprehensive electrochemical characterization of SnCl₂ in LiCl–KCl–CaCl₂ (50.5–44.2–5.3 mol%) at 685 K, including the effects of Ni 2+ and Cu + impurities. Using cyclic voltammetry (CV), chronoamperometry (CA), and chronopotentiometry (CP), we quantified Sn 2+ and Ni 2+ diffusion with exceptional agreement across methods: Sn 2+ averaged (1.03 ± 0.10) × 10 −5 cm 2 s −1 , and Ni 2+ averaged (0.75 ± 0.19) × 10 −5 cm 2 s −1 . The tight confidence-interval overlap across CV, CA, and CP strengthens confidence in these values and is uncommon in molten chloride studies. Open-circuit-potential measurements provided standard apparent reduction potentials that closely match LiCl–KCl literature, indicating minimal shift with CaCl₂ present. The Sn 2+ /Sn couple behaves as a reversible two-electron soluble–insoluble process at 685 K; the Sn 4+ /Sn 2+ couple transitions to soluble–soluble behavior near 788 K, which may correlate with the decomposition of surface bound chlorostannates, though direct characterization remains to be established. In mixed systems, Cu+/Cu overlaps Sn 2+ /Sn, limiting Cusingle bondSn electroseparation, whereas the larger potential gap between Ni 2+ /Ni and Sn 2+ /Sn supports selective Ni removal. These internally consistent transport and thermodynamic data establish a validated basis for process modeling and optimization of Sn electrorefining and impurity management in LiCl–KCl–CaCl₂.

Berzins-Delahay↗
DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

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

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2020

Heavy flavor quenching and flow: the roles of initial condition, pre-equilibrium evolution, and in-medium interaction

Within an advanced Langevin-hydrodynamics framework coupled to a hybrid fragmentation-coalescence hadronization model, we study heavy flavor quenching and flow in relativistic heavy-ion collisions. We investigate how the initial heavy quark spectrum, the in-medium energy loss and hadronization mechanisms of heavy quarks, the evolution profile of the pre-equilibrium stage, the medium flow, and the temperature dependence of heavy quark diffusion coefficients influence the suppression and elliptic flow of heavy mesons at the RHIC and the LHC. Our results show that the different modeling of initial conditions, pre-equilibrium evolution, and in-medium interactions can individually yield uncertainties of approximately 10-40% in D meson suppression and flow at a low transverse momentum. Furthermore, we also find that proper combinations of collisional versus radiative energy loss, coalescence versus fragmentation in hadronization, and the inclusion of medium flow are the most important factors for describing the suppression and elliptic flow of heavy mesons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Disordered Weyl semimetal as an array of coupled Hubbard chains

We demonstrate that a disordered magnetic Weyl semimetal may be mapped onto a two-dimensional array of coupled replicated Hubbard chains, where the Hubbard 𝑈 is directly related to the variance of the disorder potential. This is a three-dimensional generalization of a similar mapping of the two-dimensional quantum Hall plateau transition to a one-dimensional Hubbard chain. Here, we demonstrate that this mapping leads to the conclusion that the Weyl semimetal becomes a diffusive metal with a nonzero density of states at arbitrarily weak disorder, in agreement with recent work. We also discuss the absence of localization in strongly disordered Weyl semimetals from the viewpoint of this mapping.

Mott insulators↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Observation of a Knotted Electron Diffusion Region in Earth's Magnetotail Reconnection

Magnetic reconnection is a fundamental plasma process that alters the magnetic field topology and releases magnetic energy. Most numerical simulations and spacecraft observations assume a two‐dimensional diffusion region, with the electron diffusion region (EDR) embedded in the same plane as the ion diffusion region (IDR) and a uniform guide field throughout. Using observations from Magnetospheric Multiscale mission, we report a non‐coplanar, knotted EDR in Earth's magnetotail current sheet. The reconnection plane of the knotted EDR deviates by approximately 38° from that of the IDR, with the guide field exhibiting both a 38° directional shift and a twofold increase in amplitude. Moreover, the Hall magnetic field is bipolar in the EDR but quadrupolar in the IDR, indicating different Hall current structures at electron and ion scales. These observations highlight the importance of three‐dimensional effects and illustrate the complexity of multiscale coupling between the EDR and IDR during reconnection studies.

diffusion region↗
DOE OSTI2020

Inverted Region in Bimolecular Electron Transfer in Solution Enabled by Delocalization

Rate constants for bimolecular electron transfer (ET) increased with driving force, -ΔG°, reached a plateau, and then decreased in an inverted region. This rate data was described well by electron transfer theory subject to a diffusion-controlled limit. These were for ET from radical anions of polydecylthiophene (P3DT) to a series of acceptors in THF solution. When the donor was the smaller anion of quaterthiophene (T4 •- ) the inverted region was much less prominent, and still less so for when the donor was the anion of bithiophene (T2 •- ). Description of the data using ET theory identifies smaller electronic couplings for the highly-delocalized P3DT anions as enabling the inverted behavior: The presence of a Marcus inverted region is a consequence of delocalized electronic states. The results of this study further imply that electronic couplings smaller than usually found for molecules in contact could boost efficiency of energy storage by electron transfer, and identifies size-mismatch as an important concept in control of electronic couplings.

14 SOLAR ENERGY↗