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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Insights into Na + Diffusion in Silicon Modules under Operating Conditions: Measuring Low Concentrations by D-SIMS

Sodium induced shunting under an electric field is a challenging reliability issue in crystalline Si solar modules. The source of this Potential-Induced Degradation of the Shunting type (PID-s) is well understood and its influence on power loss has been intensively studied based on phenomenological models on cell or module level relating the experimental power-loss to stressing parameters (time, temperature, voltage). However, little is known about the Na ion migration kinetics, responsible for PID on a microscopic level, and its quantitative relation to the efficiency degradation. In this paper we present our investigations of sodium ion migration in Ethylene-Vinyl Acetate (EVA) and silicon through Dynamic Secondary Ion Mass Spectroscopy (D-SIMS). Each sample was annealed at field relevant temperatures from 60-90 °C to address typical migration mechanisms of common PV installations. Analysis of the SIMS migration profiles revealed a diffusivity constant D 0,EVA = 0.09 ± 0.14 cm 2 /s and an activation energy E A,EVA = 0.85 ± .04 eV for Na in EVA and diffusivities higher than extrapolated literature values in silicon (D 0,Si = (3.03 ± 2.42)x10 -5 cm 2 /s, and E A,Si = 0.98 ± 0.02 eV). The new insight will be included in a drift-diffusion based degradation model accounting for the partition coefficient across all relevant interfaces. This model can assist in predicting PID-failure in the field based on the given model stack and the diffusion of Na + through each material. Furthermore, this tool can be used for process optimization as well as material selection significantly reducing the cost and time to validate a technology.

14 SOLAR ENERGY↗
DOE OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparison of greenhouse gas emission estimates from six hydropower reservoirs using modeling versus field surveys

As with most aquatic ecosystems, reservoirs play an important role in the global carbon (C) cycle and emit greenhouse gases (GHG) as carbon dioxide (CO 2 ) and methane (CH 4 ). However, GHG emissions from reservoirs are poorly quantified, especially in temperate systems, resulting in high uncertainty. We compared reservoir C emission estimates and uncertainty of diffusive, ebullitive, and degassing pathways in six hydropower reservoirs in the southeastern United States among four data sources: two field-based surveys and two models (including the GHG Reservoir “G-res” Tool). We found that CH 4 diffusion was most similar across data sources (modeled minus observed, bias = - 21 g CO 2-eq m -2 y -1 ) and had low relative uncertainty (coefficient of variation, CV = 0.98). On the other hand, CO 2 diffusion was least consistent across data sources (bias = - 518 g CO 2-eq m -2 y -1 ). Both field surveys indicated strong negative CO 2 diffusion (i.e., CO 2 uptake) at all reservoirs, while G-res estimated positive CO 2 diffusion. By extension, total C emissions showed similar discrepancies, leading to high uncertainty in upscaling and interpreting reservoir source-sink dynamics. Finally, CH 4 ebullition had the highest relative uncertainty (CV = 2.77) due to high variability across sites. We discuss limitations of field surveys and these models, including temperature-based annualization methods, varying definitions of ebullition zones, low sampling resolution, and lack of dynamism. Future field efforts focused on capturing variability in CO 2 diffusion and CH 4 ebullition will be especially valuable in reducing uncertainty and improving models to advance our understanding reservoir GHG emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Temperature dependence of resistivity of carbon micro/nanostructures: Microscale spatial distribution with mixed metallic and semiconductive behaviors

The temperature coefficient of resistivity (θ T ) of carbon-based materials is a critical property that directly determines their electrical response upon thermal impulses. It could have metal- (positive) or semiconductor-like (negative) behavior, depending on the combined temperature dependence of electron density and electron scattering. Its distribution in space is very difficult to measure and is rarely studied. Here, for the first time, we report that carbon-based micro/nanoscale structures have a strong non-uniform spatial distribution of θ T . This distribution is probed by measuring the transient electro-thermal response of the material under extremely localized step laser heating and scanning, which magnifies the local θ T effect in the measured transient voltage evolution. For carbon microfibers (CMFs), after electrical current annealing, θ T varies from negative to positive from the sample end to the center with a magnitude change of >130% over <1 mm. This θ T sign change is confirmed by directly testing smaller segments from different regions of an annealed CMF. For micro-thick carbon nanotube bundles, θ T is found to have a relative change of >125% within a length of ~2 mm, uncovering strong metallic to semiconductive behavior change in space. In conclusion, our θ T scanning technique can be readily extended to nm-thick samples with μm scanning resolution to explore the distribution of θ T and provide a deep insight into the local electron conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development and Prospects of Halide Perovskite Single Crystal Films

Halide perovskite single crystal (HPSC) films have demonstrated extraordinary performance in solar cells, photodetectors, and lighting applications, owing to the high carrier mobility, long carrier diffusion length, tunable bandgap, and large light absorption coefficients of these materials. Here, the recent development of the major HPSC thin films is reviewed systematically, including MAPbI 3 , α-FAPbI 3 , MAPbBr 3 , CsPbBr 3 , with special emphasis on different classes of fabrication methods and thin-film characteristics. So finally, the conclusion and prospects in this field are discussed in detail. It is anticipated that the HPSC thin films will lead to a compelling future for the rising electronic applications.

MAPbI3↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2023

Intrinsic Berry curvature driven anomalous Nernst thermopower in the semimetallic Heusler alloy CoFeVSb

Understanding of spin-heat coupling mechanisms and magnetothermoelectric phenomena, including the anomalous Nernst effect (ANE), in emergent quaternary Heusler alloys is of practical importance for applications in thermal management and energy harvesting. Here, we demonstrate an intrinsic Berry curvature mediated anomalous Nernst thermopower in CoFeVSb, which orders magnetically at high temperature ( T C ≈ 850 K ) with a large saturation magnetization of ≈ 2.2 μ B / f . u . at room temperature. We show that the electron-electron elastic and electron-magnon inelastic scattering dominate longitudinal electrical transport at low temperatures ( T ≤ 50 K ), whereas the electron-phonon and electron-magnon scatterings govern it at higher T . The longitudinal thermopower is resulted mainly from the diffusive contribution with a very large longitudinal Seebeck coefficient ( 42 μ V K - 1 at 395 K). The value of the anomalous Nernst coefficient ( S ANE ) for CoFeVSb at room temperature is 0.039 μ V K - 1 which is higher than the compressively strained SrRu O 3 film ( 0.03 μ V K - 1 ) as well as the spin gapless semiconductor CoFeCrGa ( 0.018 μ V K - 1 ). On lowering T , both the ordinary Nernst coefficient and carrier mobility increase but an opposite trend is found for S ANE . Our ab initio simulations reveal the topological semimetallic nature of CoFeVSb with a pair of Weyl points. These Weyl crossings result in a significant contribution to the Berry curvature, leading to an intrinsic anomalous Hall conductivity ( σ x y AHE ) of ≈ 85 S/cm, which matches well with experiment (77 S/cm at 2 K). Our experimental findings and ab initio calculations support the dominance of the intrinsic Berry curvature in the observed ANE. The ratio of σ x y AHE to the transverse anomalous thermoelectric conductivity ( α x y ANE ) shows an increasing trend with T attaining a sizable fraction of k B e ( ≈ 0.35 k B e ) at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Novel Relaxation Time Approximation to the Relativistic Boltzmann Equation

In this work, we show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with micro- and macroscopic conservation laws if the relaxation time depends on energy or general matching conditions are applied. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Review of the first charged-particle transport coefficient comparison workshop

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4–6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron–ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. In this paper, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG

Plasma–wall interaction is one of the key research topics on the way to controlled fusion. To study the best operational designs with reduced heat and particle fluxes onto tokamak plasma facing components (PFCs) comprehensive plasma simulations are required. A recent implementation of a hybridized discontinuous Galerkin scheme into a new version of SolEdge code has the advantage of using magnetic equilibrium-free mesh. This allows us to conduct pioneering 2-D transport simulations of a full discharge in the WEST tokamak. In this work, we implemented plasma transport coefficients as functions of coordinate in the poloidal plane and neutral diffusion as a function of neutral mean free path. Moreover, the perpendicular convection flux terms were added to the code. Using the new features, a few test cases were investigated. Finally, the influence of nonconstant transport coefficients on the simulated particle and heat fluxes onto the WEST tokamak PFCs are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modularization of Ceramic Hollow Fiber Membrane Technology for Air Separation

This proposed project is aimed at studying high performance and economically competitive ceramic membrane technology for air separation and high-purity oxygen production using hollow fiber ceramic membrane stack and module technology. The design of the single permeate membrane is a hollow fiber substrate-supported thin tri-layer structure. Radially well-aligned micro-channels are embedded in the thick substrate and open at the inner surface of the substrate, enabling facile air/gas diffusion. The tri-layer structure of thin dense membrane layer (~ 10 µm) sandwiched by porous surface layer on either side is then built on the substrate using advanced fabrication process. The tri-layer structure design allows different materials to be used in different layers, where the materials of surface layers have high surface exchange coefficients while the material of dense layer has high bulk diffusivity. Such a synergetic combination leads to high permeation performance. The focus of the proposed project will be on membrane stack/module development using the developed single hollow fiber membranes, including: 1) fabrication and characterization of novel single hollow fiber membranes; 2) membrane stack design and assembly using fabricated single membranes; 3) stack modeling and analysis to guide membrane stack designs; 4) permeation performance testing and characterization of membrane stacks. The hollow fiber feature and simple sealing requirement enable very compact design of membrane stacks, which have excellent flexibility for further modularizations at different scales. The operations of such membrane stack and module may employ the exhaust heat from other components of Integrated Gasification Combined Cycle and oxy-combustion systems. Therefore, modularization of such an air separation membrane technology can be incorporated into the DOE’s REMS (radically engineered modular systems)-gasification skid and support the oxidant feed of an oxygen-blown REMS gasifier scaled to different ranges.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Wear Performance of incoloy 800HT and inconel 617 in various surface conditions for high-temperature gas-cooled reactor components

Incoloy™ 800HT and Inconel™ 617 are primary candidate structural materials for high-temperature gas-cooled reactors (HTGRs). Impurities (e.g., H 2 O, CH 4 ) in the coolant induce corrosion reactions that affect tribological behavior. We investigate the high-temperature tribological performance of the alloys under surface conditions including: conditioned in H 2 O-containing helium at elevated temperatures, hardened via the carbon-diffusion Kolsterising® treatment, and conditioned after hardening. Initial friction coefficients and volume removed serve as metrics for wear performance. In this work, we find that wear mechanisms including scuffing, glaze-oxide formation, and oxide-layer wear and break-through are governed by surface hardness, chemistry, and applied load. The formation of a glaze-oxide layer is achieved with an oxide of sufficient mechanical integrity if the oxide-substrate interface is strong and the load is low.

42 ENGINEERING↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗