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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1994

Diffusive modification of primary zoning in olivine cores in Angrite LEW 87051

Angrite LEW 87051 consists of large olivine crystals set in a fine-grained groundmass that clearly represents a crystallized melt. A few olivines contain Ca-poor, Cr-rich cores that crystallized from a very different melt than the outer part of the crystals constituting the majority of olivine in LEW 87051. We evaluate a model in which the cores formed through fractional crystallization of one melt, then were incorporated into a different melt as xenocrysts, whereupon the original zoning patterns were modified by diffusion. Using a similar approach, we calculate zoning patterns for the cores that would result from perfect fractional crystallization, compare them with the observed zoning, and determine whether the differences could result from diffusive modification consistent with known diffusion rates for Ca, Mn, and Cr. Using distribution coefficients from the 1400 C, IW + 1 experiments, we computed CaO, Cr2O3, and MnO abundances in the hypothetical parent melt by inverting the olivine at the centers of the cores. We further assumed that the primary zoning profile for CaO is essentially unmodified, because the diffusion rate of Ca in olivine is slow. We carried out the fractional crystallization calculation until the calculated Ca content was that observed at break in zoning profiles at the outer edge of the core. We then normalized the distance of this calculated profile to the length of the observed profile in the olivine core and calculated profiles for MnO and Cr2O3. The CaO zoning profile agrees well with the observed profile. The observed MnO profile is slightly higher than the calculated profile near the edge. Diffusion calculations indicate that reversal of the general trend of primary zoning through diffusion would require that diffusion of Cr is 5-10x faster than that of Mn.

Mikouchi, T.↗
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↗
NASA NTRS1980

Liquid-phase electroepitaxy - Dopant segregation

A theoretical model is presented which accounts for the dopant segregation in liquid-phase electroepitaxy in terms of dopant transport in the liquid phase (by electromigration and diffusion), the growth velocity, and the Peltier effect at the substrate-solution interface. The contribution of dopant electromigration to the magnitude of the effective segregation coefficient is dominant in the absence of convection; the contribution of the Peltier effect becomes significant only in the presence of pronounced convection. Quantitative expressions which relate the segregation coefficient to the growth parameters also permit the determination of the diffusion constant and electromigration mobility of the dopant in the liquid phase. The model was found to be in good agreement with the measured segregation characteristics of Sn in the electroepitaxial growth of GaAs from Ga-As solutions. For Sn in Ga-As solution at 900 C the diffusion constant was found to be 4 x 10 to the -5 sq cm/s and the electromigration velocity (toward the substrate with a positive polarity 2 x 10 to the -5 cm/s current density of 10 A/sq cm.

Lagowski, J.↗
NASA NTRS1986

Comparison of analytical and experimental performance of a wind-tunnel diffuser section

Wind tunnel diffuser performance is evaluated by comparing experimental data with analytical results predicted by an one-dimensional integration procedure with skin friction coefficient, a two-dimensional interactive boundary layer procedure for analyzing conical diffusers, and a two-dimensional, integral, compressible laminar and turbulent boundary layer code. Pressure, temperature, and velocity data for a 3.25 deg equivalent cone half-angle diffuser (37.3 in., 94.742 cm outlet diameter) was obtained from the one-tenth scale Altitude Wind Tunnel modeling program at the NASA Lewis Research Center. The comparison is performed at Mach numbers of 0.162 (Re = 3.097x19(6)), 0.326 (Re = 6.2737x19(6)), and 0.363 (Re = 7.0129x10(6)). The Reynolds numbers are all based on an inlet diffuser diameter of 32.4 in., 82.296 cm, and reasonable quantitative agreement was obtained between the experimental data and computational codes.

Shyne, R. J.↗
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↗
NASA NTRS1973

Experimental and analytical study of a conically diffused flow with a nearly separated boundary layer

Turbulence measurements were obtained in the nearly separated flow in a 13 deg total angle of divergence conical diffuser coupled to a constant area tailpipe. Air at 207 newtons per square centimeter and 308 K provided an inlet velocity of about 51 meters per second at an inlet unit Reynolds number of 63.7 million per meter. Very high longitudinal turbulence intensities accompanied the diffusion process with peak values approaching 40 percent when normalized by the local centerline velocity. Predictions of the pressure recovery coefficient using a mixing length concept were good in the early stages of diffusion. In the latter stages of diffusion satisfactory predictions of the pressure recovery were obtained with an empirical method.

Boldman, D. R.↗
NASA NTRS2002

Solutal Convection Around Growing Protein Crystal and Diffusional Purification in Space

This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities. Onset of convection was considered analytically and numerically. Crystal growth was characterized by slow surface incorporation kinetics, i.e. growth kinetic coefficient beta (cm/s) small as compared to the typical bulk diffusion rate, D(sub 1)/h, where D(sub 1) is diffusivity of major crystallizing protein and h is the crystal size. Scaling type analysis predicted two laws on how the convection rate, v, essentially the Peclet number, Pe exactly equal to vh/D(sub 1), depends on dimensionless kinetic coefficient a exactly equal to beta h/D(sub 1). Namely: Pe = C(sub 2/5)(aRa(sup 2/5)) and Pe = C(sub 1) aRa. Here, Reynolds number Ra = rho(sub 1)(sup 0)gh(sup 3)(rho(sub p) - rho(sub w))/rho(sup p)rho(sub 1)vD(sub 1), v being solution viscosity. The constants C(sub 2/5), exactly equal to 0.28 and C(sub 1) exactly equal to 10(exp -2) found from the full scale computer simulation for a cylindrical crystal inside big cylindrical vessel. The linear boundary conditions connecting protein and impurity concentration at the interface with the flux to/from the interface was applied. No-slip condition for Navier-Shocker equations was employed. With these conditions, flow and concentration distributions were calculated. Validity of the Pe(Ra) dependencies follows for wide range of parameters for which numerical calculations have been accomplished and presented by various points.

Chernov, A. A.↗
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↗
NASA NTRS1973

The application of far IR lasers to the investigation of the optical properties of materials

The effect of diffusion on the saturation and absorption properties of molecular gases is considered. The analysis starts with a simple derivation of the multi-level rate equations. Incorporating the diffusion term, calculations are carried out both for the homogeneous and inhomogeneous broadening cases. Both graphical and analytical solutions are obtained for the absorption coefficient and the saturation intensity. It is shown that the effect of diffusion is characterized by a diffusion function, which varies inversely with the radius of the beam at low pressures and small beam sizes. The analysis is done in terms of two parameters k1a and k2a, which are simple functions of the intensity of the beam I, the effective relaxation rate w, and the effective diffusion constant D the a is designated as the radius of the beam.

Stafsudd, O. M.↗
NASA NTRS1973

A thick-thin approximation for radiative transfer.

First- and second-order formulations of a thick-thin approximate technique for application to the solution of energy transport when nongray radiation is important are presented. The virtue of the technique is that it eliminates the cumbersome integral terms which usually appear. It is shown that the temperature distributions calculated with the technique agree well with the results of more rigorous calculations if a proper value of the critical optical depth, which separates the thick and thin regions of the spectrum, is used. The proper value of the critical optical depth is less than one and may be less than 0.5. It is also shown that the jump boundary conditions which are frequently employed with the Rosseland, or diffusion, approximation should incorporate the Planck mean absorption coefficient for a nongray gas. This is in contradiction to suggestions appearing in the literature.

Olstad, W. B.↗
NASA NTRS1979

The hydrogenation of metals upon interaction with water

Hydrogen evolution at 600 deg and 5 x 10 to the 7th power - 10 to the 6th power torr from AVOOO Al samples, which were pickled in 10 percent NaOH, is discussed. An H evolution kinetic equation is derived for samples of equal vol. and different surfaces (5 and 20 sq cm). The values of the H evolution coefficient K indicated an agreement with considered H diffusion mechanism through an oxide layer. The activation energy for the H evolution process, obtained from the K-temp. relation, was 13,000 2000 cal/g-atom.

Andreyev, L. A.↗
NASA NTRS1980

Effects of grain boundaries in polycrystalline solar cells

In polycrystalline solar cells, the short-circuit current is reduced (compared to single-crystal cells) due to recombination of minority carriers at the grain boundaries. The magnitude of this reduction is calculated for the case of a monochromatic beam of light, assuming that both the light beam and the grain boundaries are perpendicular to a p-n junction. The photoresponse of the junction is calculated in terms of the absorption coefficient of the light, the minority-carrier diffusion length, the surface recombination velocity at the grain boundary and the distance of the light beam from the grain boundary. The analysis gives an expression for an effective grain boundary width which depends on wavelength. This width determines the reduction in short-circuit current of solar cells made from material with columnar-type grains.

Zook, J. D.↗
NASA NTRS1982

Stagnation point heat transfer augmentation due to free stream turbulence

A model has been proposed for the momentum eddy diffusivity induced by free stream turbulence intensity and integral length scale. The eddy diffusivity model is applied to the stagnation point of a cylinder situated in a uniform crossflow in the presence of free stream turbulence. A numerical solution of the governing momentum and energy equations with the proposed eddy diffusivity model yielded results for the skin friction coefficient and the Nusselt number. The numerical predictions of the present work are compared with experimental data and the agreement between the two is seen to be very good.

Gorla, R. S. R.↗