Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗

Sulphate resistance of low‐clinker engineered cementitious composites examined by MicroXRF imaging

Abstract Engineered cementitious composites (ECC) are a class of high‐performing fibre‐reinforced cementitious materials recognised for their increased ductility and durability compared to conventional cement‐based materials, owing to their autogenously controlled tight crack widths, even when subjected to high strains. To reduce ECC's environmental impact, this research examines the use of a low‐clinker binder − limestone‐calcined clay cement (LC3) − as an alternative to portland cement (PC), along with fly ash to further reduce the clinker proportion and the embodied CO 2 of the formulations. In conventional concrete, LC3 hydrates to a denser microstructure resulting from the synergistic reaction between limestone and calcined clay. At the lower water contents typical of ECC and with the presence of fly ash, the influence of the binder composition on the microstructure is difficult to anticipate. To examine the influence of these compositional variables on microstructure, permeability and durability, the sulphate resistance of LC3‐based ECC is explored. Specifically, the ECC‐LC3 blends are designed with high clinker replacement rate of 75% by mass of binder and contain either conventional fly ash or reclaimed fly ash at 50% by mass of binder. Expansion of ECC‐LC3 samples subjected to standard sodium sulphate test conditions was measured up to 12 months and the depth of penetration of sulphates into the ECC‐LC3 of varying compositions was quantified using micro‐X‐Ray Fluorescence (microXRF) imaging and modelling. The expansion results show that the ECC‐LC3 formulations performed better than the PC samples and can provide adequate resistance to external sulphate attack, even when reclaimed fly ashes are used in place of the conventional ash. In addition, the shallow penetration of sulphate into these cementitious composites demonstrates the low diffusion coefficients values that were determined using the quantitative data from MicroXRF imaging.

Microscopy↗

Physically Based Modeling of PEMFC Cathode Catalyst Layers: Effective Microstructure and Ionomer Structure–Property Relationship Impacts

Abstract This work presents a pseudo-two-dimensional proton exchange membrane fuel cell (PEMFC) model incorporating Nafion ionomer structure–property relationships in the cathode catalyst layer (CL) to capture and explain losses at low Pt loading. Structural data from neutron reflectometry and thin film Nafion conductivity measurements predict variations in the oxygen diffusion coefficient and ionic conductivity with changing CL ionomer thickness and Pt loading. By including these structure–property relationships, predicted polarization curves agree closely with previously published experimental data from cells with Pt loadings between 0.025 and 0.2 mg/cm2. Results demonstrate that structure–property relationships based on physically measurable ionomer and CL properties provide a feasible interpretation of PEMFC CL phenomena for a range of Pt loadings and help explain previously unaccounted-for losses at low Pt. Results also show that simulations must account for surface species coverage variations in order to properly capture the kinetic losses. Finally, results suggest that an increase in ionomer thickness surrounding the C/Pt surfaces may lead to improved cell performance due to improved ionic conductivity.

Electrochemistry↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗

In-Situ and Ex-Situ Studies on the Morphology Changes of Polymer Pipeline Materials for Use in Hydrogen Gas Environments

The US natural gas infrastructure is a national asset that could be used to deliver hydrogen and hydrogen blends of natural gas as a pathway to reduce carbon emissions. The distribution system comprises nearly 50% plastic pipe composed of medium- and high-density polyethylene materials (MDPE and HDPE). While these materials perform adequately for natural gas, research on their hydrogen compatibility is essential to understand if any immediate and long-term risks are associated with hydrogen addition. The Blended Gas CRADA, a HyBlend project, has established a comprehensive test method for evaluating MDPE and HDPE of various plastic resin compositions of pipeline material in pure hydrogen and 20% hydrogen/80% methane blends. Both in-situ and ex-situ measurements were performed to capture hydrogen-induced changes in the polyethylene material's crystalline, amorphous and their interphase regions. We investigated MDPE and HDPE pipeline materials made from different polymer resin systems to evaluate the effects of hydrogen gas. The materials were characterized by their density, diffusion coefficient, free volume ratio, and degree of crystallinity. Various advanced characterization methods, including in situ NMR, ex situ XRD, ex situ DSC, and ex situ TDA, were used to analyze the effects of changes in crystalline, amorphous, and interphase regions due to gas exposure. Time-dependent post-decompression quasi-static tensile tests were conducted to explore the effects of gas exposure time on the mechanical behavior of the pipe materials. This work will highlight the time sensitivities during and after gas exposure. The correlation between gas-induced polyethylene morphology changes and the associated material performance will be addressed for the intended applications. These studies will show that polyethylene resin composition and material exposure are important factors when considering whether hydrogen gas affects pipeline materials positively or negatively.

Simmons, Kevin L.↗

Visualizing energy transfer between redox-active colloids

Redox-active colloids (RACs) represent a novel class of energy carriers that exchange electrical energy upon contact. Understanding contact-mediated electron transfer dynamics in RACs offers insights into physical contact events in colloidal suspensions and enables quantification of electrical energy transport in nonconjugated polymers. Redox-based electron transport was directly observed in monolayers of micron-sized RACs containing ethyl-viologen side groups via fluorescence microscopy through an unexpected nonlinear electrofluorochromism that is quantitatively coupled to the redox state of the colloid. Via imaging studies, using this electrofluorochromism, the apparent charge transfer diffusion coefficient D CT of the RAC was easily determined. The visualization of energy transport within suspensions of redox-active colloids was also demonstrated. Our work elucidates fundamental mechanisms of energy transport in colloidal systems, informs the development of next-generation redox flow batteries, and may inspire new designs of smart active soft matter including conductive polymers for applications ranging from electrochemical sensors and organic electronics to colloidal robotics.

Qu, Alan Subing [Univ. of Illinois at Urbana-Champ↗

Learning Optimal Multigrid Smoothers via Neural Networks

Multigrid methods are one of the most efficient techniques for solving large sparse linear systems arising from partial differential equations (PDEs) and graph Laplacians from machine learning applications. One of the key components of multigrid is smoothing, which aims at reducing high-frequency errors on each grid level. However, finding optimal smoothing algorithms is problem-dependent and can impose challenges for many problems. In this paper, we propose an efficient adaptive framework for learning optimized smoothers from operator stencils in the form of convolutional neural networks (CNNs). Here, the CNNs are trained on small-scale problems from a given type of PDEs based on a supervised loss function derived from multigrid convergence theories and can be applied to large-scale problems of the same class of PDEs. Numerical results on anisotropic rotated Laplacian problems and variable coefficient diffusion problems demonstrate improved convergence rates and solution time compared with classical hand-crafted relaxation methods.

97 MATHEMATICS AND COMPUTING↗

T -matrix analysis of static Wilson line correlators from lattice QCD at finite temperature

Here, we utilize a previously constructed thermodynamic T -matrix approach to the quark-gluon plasma (QGP) to derive constaints on its input by using results on Wilson line correlators (WLCs) of a static quark-antiquark pair from 2 + 1-flavor lattice-QCD (lQCD) computations with realistic pion mass. The self-consistent T -matrix results, which include previous constraints from the lQCD equation of state in the light-parton sector, can describe the lQCD data for WLCs fairly well once refinements of its driving kernel are applied. In particular, the input potential requires less screening than what has been inferred from previous T-matrix analyses. Pertinent predictions for the spectral and transport properties of the QGP are discussed, including the spatial diffusion coefficient for heavy quarks; the latter turns out to have a rather weak temperature dependence, in approximate agreement with recent unquenched lQCD results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

QLBT: a linear Boltzmann transport model for heavy quarks in a quark-gluon plasma of quasi-particles

Abstract We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor $$R_\mathrm {AA}$$ R AA and the elliptic flow $$v_2$$ v 2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson $$R_\mathrm {AA}$$ R AA and $$v_2$$ v 2 , we extract the heavy quark transport parameter $$\hat{q}$$ q ^ and diffusion coefficient $$D_\mathrm {s}$$ D s in the temperature range of $$1-4~T_\mathrm {c}$$ 1 - 4 T c , and compare them with the lattice QCD results and other phenomenological studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unraveling TeV halos with the Cherenkov Telescope Array

Pulsars are observed to emit bright and spatially extended gamma-ray emission at multi-TeV energies. These so-called "TeV halos" are now understood to be a nearly universal feature of middle-aged pulsars. However, many of the key physical processes that govern these systems, particularly those affecting particle diffusion, remain poorly constrained. We aim to evaluate the ability of the Cherenkov Telescope Array (CTA) to probe the physical properties of TeV halos, with a focus on the nearby and well-studied case of the Geminga pulsar. We simulate gamma-ray emission from various TeV halo models, incorporating different assumptions for the injected electron spectrum, spin-down evolution, and energy-dependent diffusion. These models are then used to forecast CTA's sensitivity to spectral and spatial differences, based on realistic mock observations and instrument response simulations. We find that CTA will be able to distinguish between a wide range of TeV halo models that are currently consistent with existing data. In particular, CTA observations can constrain the normalization, energy dependence, and spatial extent of the diffusion coefficient surrounding Geminga, as well as the spectral shape of the injected electron population.

79 ASTRONOMY AND ASTROPHYSICS↗

Coarse-grained explicit-solvent molecular dynamics simulations of semidilute unentangled polyelectrolyte solutions

In this study, we present results from explicit-solvent coarse-grained molecular dynamics (MD) simulations of fully charged, salt-free, and unentangled polyelectrolytes in semidilute solutions. The inclusion of a polar solvent in the model allows for a more physical representation of these solutions at concentrations, where the assumptions of a continuum dielectric medium and screened hydrodynamics break down. The collective dynamic structure factor of polyelectrolytes, S(q, t), showed that at q > q*, where q* = 2π/ξ is the polyelectrolyte peak in the structure factor S(q) and ξ is the correlation length, the relaxation time obtained from fits to stretched exponential was $\tau$ KWW ~ q -3 , which describes unscreened Zimm-like dynamics. This is in contrast to implicit-solvent simulations using a Langevin thermostat where $\tau$ KWW ~ q -2 . At q < q*, a crossover region was observed that eventually transitions to another inflection point $\tau$ KWW ~ q -2 at length scales larger than ξ for both implicit- and explicit-solvent simulations. The simulation results were also compared to scaling predictions for correlation length, ξ ~ c$-½\atop{p}$, specific viscosity, η sp ~ c$½\atop{p}$, and diffusion coefficient, D ~ c$0\atop{p}$, where c p is the polyelectrolyte concentration. The scaling prediction for ξ holds; however, deviations from the predictions for η sp and D were observed for systems at higher c p , which are in qualitative agreements with recent experimental results. This study highlights the importance of explicit-solvent effects in molecular dynamics simulations, particularly in semidilute solutions, for a better understanding of polyelectrolyte solution behavior.

36 MATERIALS SCIENCE↗

Electrochemical Characterization of Bromine Reduction to Tribromide in Individual Nitrobenzene-in-Water Emulsion Droplets

We report an electrochemical study of bromine reduction to tribromide in a single nitrobenzene (NB) droplet during a single collision event. The feasibility of this study is based on the favorable distribution of Br 2 in NB rather than in water, which was demonstrated by voltammetry at the NB/water interface. The NB-in-water emulsions containing the ionic liquids trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)amide (ILPA) as the supporting electrolyte and sodium dodecyl sulfate (SDS) as the surfactant were prepared by a high-power ultrasonication. Its single droplet as an attoliter reactor, colliding on a 25 μm Pt UME, was monitored by chronoamperometry. The spike-type current transients of Br 2 reduction can elucidate the size and size distribution of NB, and measure the diffusion in the droplet as analyzed by a bulk electrolysis model. The frequency of collisions can provide diffusion coefficient of the droplets to the UME. Moreover, the i-t decay behavior can be simulated, which fits well with the experimental one.

25 ENERGY STORAGE↗

The Effect of Proton Conductivity of Fe–N–C–Based Cathode on PEM Fuel cell Performance

A model–based impedance spectroscopy is used to determine proton conductivity, oxygen transport parameter, double layer capacitance and oxygen reduction reaction (ORR) Tafel slope in the Fe–N–C cathode catalyst layer (CCL) of a PEM fuel cell. Experimental spectra of two cells differing by the membrane thickness only are processed using a physics–based model for PEMFC impedance. The spectra have been measured in the range of current densities from 25 to 800 mA cm -2 . The ORR Tafel slope of both the cells shows almost linear growth with the current density. In one of the cells, the CCL proton conductivity σp strongly decays at the current density of 100 mA cm -2 ; this decay is accompanied by the step growth of the double layer capacitance. Other minor variations of proton conductivity and double layer capacitance with the cell current occur also in a counterphase; presumed origin of this effect is discussed. The oxygen diffusion coefficient in the cathode exhibits explosive growth with the cell current. We attribute this effect to formation of temperature and pressure gradients in the CCL due to strongly non–uniform distribution of ORR rate in the electrode.

25 ENERGY STORAGE↗

Electrochemistry of the NaI-AlBr 3 Molten Salt System: A Redox-Active, Low-Temperature Molten Salt Electrolyte

NaI-AlBr 3 is a very appealing low melting temperature (<100 °C), salt system for use as an electrochemically-active electrolyte. This system was investigated for its electrochemical and physical properties with focus to energy storage considerations. A simple phase diagram was generated; at >100 °C, lower NaI concentrations had two partially miscible liquid phases, while higher NaI concentrations had solid particles. Considering the fully molten regime, electrical conductivities were evaluated over 5–25 mol% NaI and 110 °C–140 °C. Conductivities of 6.8–38.9 mS cm −1 were observed, increasing with temperature and NaI concentration. Effective diffusion coefficients of the I − /I 3 − redox species were found to decrease with both increasing NaI concentration and increasing applied potential. Regardless, oxidation current density at 3.6 V vs Na/Na + was observed to increase with increasing NaI concentration over 5–25 mol%. Finally, the critical interface between the molten salt electrolyte and electrode materials was found to significantly affect reaction kinetics. When carbon was used instead of tungsten, an adsorbed species, most likely I 2 , blocked surface sites and significantly decreased current densities at high potentials. This study shows the NaI-AlBr 3 system offers an attractive, low-temperature molten salt electrolyte that could be useful to many applied systems, though composition and electrode material must be considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potassium-Containing α-MnO 2 Nanotubes: The Impact of Hollow Regions on Electrochemistry

α-MnO 2 type materials have been studied as electrode materials in rechargeable batteries and electrocatalysts due to their 2 × 2 tunneled crystal structures capable of accommodating cations and their tunable physiochemical properties. In this study, we deliberately synthesized K + containing α-MnO 2 (K 0.9 Mn 8 O 16 ) hollow nanotubes varying the dimensions of the hollow regions and level of surface defects. The K 0.9 Mn 8 O 16 nanotube material samples have similar crystallinity, thermal stability, and average Mn oxidation state. Oxygen surface defects in the hollow regions were revealed through detailed studies using electron energy loss spectroscopy. The impact of the hollow regions and associated surface defects on the electrochemistry of K x Mn 8 O 16 were investigated using cyclic voltammetry, galvanostatic intermittent titration technique, and galvanostatic cycling. The K 0.9 Mn 8 O 16 nanotubes with a large hollow region (~30 nm) and higher level of surface defects show higher apparent lithium ion diffusion coefficients and lower polarization compared to the nanotubes with a small hollow region (~10 nm). In-situ lithiation demonstrated that the dimensions of the nanotube walls expanded, but the hollow region did not change in size as result of lithiation. Furthermore, this research demonstrates that tuning particle architecture and surface defects can positively impact functional behavior of electrochemical storage materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transport of Ligand Coordinated Iron and Chromium through Cation-Exchange Membranes

Fluxes of negatively charged ligand-coordinated iron, Fe(CN) 6 3/4- , and chromium, CrPDTA 1/2− , through two promising commercial cation-exchange membranes, Aquivion E87–05S and Fumasep E-620(K), were measured as functions of current density. The magnitude of the partial current density reached a maximum of − 43 μ A cm −2 at the maximum applied current density magnitude of − 43 mA cm −2 for Fe(CN) 6 3/4− transport through Aquivion, or 0.1% of the total current density. Fumasep E-620(K) blocks practically all crossover of both compounds. Both membranes sorb more Fe(CN) 6 3/4− and CrPDTA 1/2− than predicted by Donnan equilibrium, and low crossover rates can be attributed mainly to slow diffusion, not charge-based rejection of co-ions. The magnitude of the diffusion coefficient appears to correlate with hydraulic permeability. Although Aquivion E87–05S and Fumasep E-620(K) have significant and observable differences in membrane crossover rates, cells built with the DI-soaked membranes offer similarly high coulombic efficiency, indicating the relatively small contribution that crossover makes to inefficiency over a single cycle.

25 ENERGY STORAGE↗