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

A comprehensive diffusion mobility database comprising 23 elements for magnesium alloys

We report that reliable experimental diffusion coefficients of 10 key alloying elements in Mg obtained by the present authors together with experimental data in the literature enabled us to perform a systematic test of the reliability of diffusion coefficients obtained from DFT calculations. The computed activation energy values were found to be quite accurate (mostly within 0.2 eV) but the computed pre-factors were less reliable. Such insights allowed us to develop a practical and yet robust strategy to perform diffusion mobility assessments by adopting the computed activation energy while fitting only the pre-factor when available experimental data are limited to a narrow temperature range. The overall good agreement between the DFT data and experimental data also gave us the confidence to employ the computed data for those that were still missing or inaccessible from experimental measurements. A systematic assessment of both the measured and computed diffusion data in hcp Mg was performed using the above holistic approach to yield the most comprehensive open Mg mobility database to date, comprising 23 elements (Mg, Ag, Al, Be, Ca, Cd, Ce, Cu, Fe, Ga, Gd, In, La, Li, Mn, Nd, Ni, Pu, Sb, Sn, U, Y, Zn). This more reliable mobility database will contribute to future development of advanced Mg alloys. The holistic approach developed in this study will be very beneficial to the future establishment of reliable mobility databases for other alloy systems as well.

36 MATERIALS SCIENCE↗
NASA NTRS1979

The dynamical halo and the variation of cosmic-ray path length with energy

It is shown that the dynamical halo model offers a natural explanation for the form of the variation of the cosmic-ray path length with energy. The variation above about 1 GeV per nucleon can be understood as due to the variation of the diffusion coefficient, and hence the resident time, with energy. The flattening of the curve below 1 GeV per nucleon is seen to mark a transition to a convection-dominated regime where the diffusion coefficient is no longer the determining parameter. A fit to the observations yields a halo outflow velocity of 8 km/s. An attempt is made to determine the overall scale of the halo and the diffusion coefficient using recent Be-10 flux measurements, but the data do not agree well enough to pin down these variables to within less than four or five orders of magnitude.

Jones, F. C.↗
NASA NTRS1981

Solute partitioning under continuous cooling conditions as a cooling rate indicator

A model of solute partitioning in a finite body under conditions of continuous cooling is developed for the determination of cooling rates from concentration profile data, and applied to the partitioning of zirconium between ilmenite and ulvospinel in the Apollo 15 Elbow Crater rocks. Partitioning in a layered composite solid is described numerically in terms of concentration profiles and diffusion coefficients which are functions of time and temperature, respectively; a program based on the model can be used to calculate concentration profiles for various assumed cooling rates given the diffusion coefficients in the two phases and the equilibrium partitioning ratio over a range of temperatures. In the case of the Elbow Rock gabbros, the cooling rates are calculated from measured concentration ratios 10 microns from the interphase boundaries under the assumptions of uniform and equilibrium initial conditions at various starting temperatures. It is shown that the specimens could not have had uniform concentrations profiles at the previously suggested initial temperature of 1350 K. It is concluded that even under conditions where the initial temperature, grain sizes and solute diffusion coefficients are not well characterized, the model can be used to estimate the cooling rate of a grain assemblage to within an order of magnitude.

Onorato, P. I. K.↗
NASA NTRS1974

M558 radioactive tracer diffusion

This experiment was performed in Skylab 3 with two objectives in mind. First, the experimental self-diffusion coefficients for liquid zinc were to be determined in a convection-free environment. Secondly the reduction in convective mixing in earth gravity by going into the zero-gravity environment of space was to be estimated. The experiment was designed to utilize high temperatures and linear thermal gradients provided by the M518 Multipurpose Electric Furnace, and the radioactivity of zinc-65 of 245-day half-life to investigate self-diffusion in liquid zinc. The distribution of zinc-65 tracer, after melting, maintaining at soak temperature for 1 hour of soak time and then resolidifying, was obtained by sample sectioning. The concentration of activity of each section (microcurie-gram) was plotted against positions along the sample axial and radial position. Experimental data and theoretical results from solution of Fick's law of diffusion in one dimensional were compared. Samples tested on earth showed very rapid diffusion. Diffusion coefficient in unit gravity was 50 times the zero-gravity diffusion coefficient of Skylab.

Ukanwa, A. O.↗
NASA NTRS1976

Absorption of radiation belt particles by the inner satellites of Jupiter

The study of trapped particle absorption by the inner Jovian satellites is reviewed from the viewpoint of radiation belt physics. Both pre- and post-Pioneer work is discussed but the emphasis is on methods used to deduce radial diffusion coefficients of particle transport from particle data. The phenomenon of particle absorption as observed by experiments on Pioneers 10 and 11 is considered; absorption effects are found to depend on the satellite, and on particle energy and species. Approximate diffusion coefficients derived from the data are found to follow a steeper spatial dependence than previously expected. The assumptions and limitations of absorption analysis and diffusion coefficient estimation are pointed out.

Mogro-Campero, A.↗
NASA NTRS2014

Improved Algorithms for Accurate Retrieval of UV - Visible Diffuse Attenuation Coefficients in Optically Complex, Inshore Waters

Photochemical processes driven by high-energy ultraviolet radiation (UVR) in inshore, estuarine, and coastal waters play an important role in global bio geochemical cycles and biological systems. A key to modeling photochemical processes in these optically complex waters is an accurate description of the vertical distribution of UVR in the water column which can be obtained using the diffuse attenuation coefficients of down welling irradiance (Kd()). The Sea UV Sea UVc algorithms (Fichot et al., 2008) can accurately retrieve Kd ( 320, 340, 380,412, 443 and 490 nm) in oceanic and coastal waters using multispectral remote sensing reflectances (Rrs(), Sea WiFS bands). However, SeaUVSeaUVc algorithms are currently not optimized for use in optically complex, inshore waters, where they tend to severely underestimate Kd(). Here, a new training data set of optical properties collected in optically complex, inshore waters was used to re-parameterize the published SeaUVSeaUVc algorithms, resulting in improved Kd() retrievals for turbid, estuarine waters. Although the updated SeaUVSeaUVc algorithms perform best in optically complex waters, the published SeaUVSeaUVc models still perform well in most coastal and oceanic waters. Therefore, we propose a composite set of SeaUVSeaUVc algorithms, optimized for Kd() retrieval in almost all marine systems, ranging from oceanic to inshore waters. The composite algorithm set can retrieve Kd from ocean color with good accuracy across this wide range of water types (e.g., within 13 mean relative error for Kd(340)). A validation step using three independent, in situ data sets indicates that the composite SeaUVSeaUVc can generate accurate Kd values from 320 490 nm using satellite imagery on a global scale. Taking advantage of the inherent benefits of our statistical methods, we pooled the validation data with the training set, obtaining an optimized composite model for estimating Kd() in UV wavelengths for almost all marine waters. This optimized composite set of SeaUVSeaUVc algorithms will provide the optical community with improved ability to quantify the role of solar UV radiation in photochemical and photobiological processes in the ocean.

UV↗
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
DOE OSTI2023

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamics of the O ( 4 ) critical point in QCD: Critical pions and diffusion in model G

We present a detailed study of the finite momentum dynamics of the O ( 4 ) critical point of QCD, which lies in the dynamic universality class of “model G.” The critical scaling of the model is analyzed in multiple dynamical channels. For instance, the finite momentum analysis allows us to precisely extract the pion dispersion curve below the critical point. The pion velocity is in striking agreement with the predictions relation and static universality. The pion damping rate and velocity are both consistent with the dynamical critical exponent ζ = 3 / 2 of model G. Similarly, although the critical amplitude for the diffusion coefficient of the conserved O ( 4 ) charges is small, it is clearly visible both in the restored phase and with finite explicit symmetry breaking, and its dynamical scaling is again consistent with ζ = 3 / 2 . We determine a new set of universal dynamical critical amplitude ratios relating the diffusion coefficient to a suitably defined order parameter relaxation time. We also show that in a finite volume simulation, the chiral condensate diffuses on the coset manifold in a manner consistent with dynamical scaling, and with a diffusion coefficient that is determined by the transport coefficients of hydrodynamic pions. Finally, the amplitude ratios (together with other nonuniversal amplitudes also reported here) compile all relevant information for further studies of model G both in and out of equilibrium. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Dynamics of the middle atmosphere in winter (Dynamics). Interrelation between the different variations of turbulent diffusion and ionospheric absorption originating in the middle atmosphere

The turbulent diffusion coefficient was computed from the parameters of sporadic E layers using the wind shear theory of midlatitude sporadic E and models of the ionosphere as well as that of the neutral upper atmosphere. The turbulent diffusion coefficient obtained for the period of circulation disturbances associated with stratospheric warmings and for the intervals of the winter anomaly indicate changes similar to the ionospheric absorption of radio waves, in the former case decreased, in the latter case increased values. This may hit at the role of turbulent transport in the formation of these anomalies. On the basis of these findings, a seasonal variation of the turbulent diffusion coefficient having a minimum in summer and an increase of this parameter with increasing geomagnetic activity are anticipated.

Bencze, P.↗
NASA NTRS1980

The effect of convection on the propagation of relativistic galactic electrons

The paper presents exact analytical solutions which describe steady-state transport of relativistic electrons subject to diffusion, convection, and radiation losses. The diffusion coefficient is spatially and energy dependent; the bulk convection velocity is considered to be constant and directed outward from the galactic plane. The electron spectral behavior with respect to galactic height and energy indicates that convection provides a 'break' in the steady-state energy spectrum without recourse to distribution of sources. Radio observations of NGC 891 and NGC 4631 galaxies can be used to estimate their diffusion coefficient and its energy dependence. Application of analytical solutions to our Galaxy indicates that the diffusion coefficient and outflow velocity deduced from relativistic electron measurements and radio observations of their synchrotron radiation agree with values based on the transport of the cosmic ray nucleonic component.

Lerche, I.↗
DOE OSTI2021

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

The Martian climate and energy balance models with CO2/H2O atmospheres

The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.

Hoffert, M. I.↗
NASA NTRS1977

Variations in eddy diffusion and associated transport processes

Variation in the eddy diffusion coefficient constitutes a momentum source. This variation near the turbopause affects both the composition and the wind and temperature fields of the thermosphere. Two types of calculations are carried out, for horizontal winds considered self-consistently in a three-dimensional dynamic model and horizontal winds artificially suppressed to simulate a one-dimensional model. Momentum transfer associated with variations in eddy diffusion coefficients represents the source function in the analysis presented. It is shown that for long-term variations in the thermosphere, the compositional dependence on variations in eddy diffusion coefficients can be entirely different in one and three dimensional models with horizontal winds included. Horizontal winds can greatly reduce the effects of eddy diffusion in general and the ratio between He and O amplitudes in particular.

Mayr, H. G.↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
DOE OSTI2021

Electrochemical measurement and analysis of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 in molten eutectic LiCl-KCl

For this work, electrochemical measurements of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 were obtained in molten LiCl-KCl eutectic at 773 K using a three-electrode configuration. Tungsten rods were used as the working and counter electrodes. Two quasi- and one reference electrode(s) (glassy carbon, nickel wire, and Ni/NiO) were used to collect electrochemical measurements of the analyte salts. Cyclic voltammetry data were analyzed to calculate values of diffusion coefficient, exchange current density, and charge transfer resistance and were determined to be on the order of 10 –5 cm 2 s –1 , 10 –2 A cm –2 , and 10 0 Ω; respectively. In general, the values of diffusion coefficients were found to be consistent with those reported of high-temperature molten salts. Relatively large values of exchange current density corresponded with smaller values of charge transfer resistance. These values were found to be reasonable in comparison to results available in the literature. The measured cyclic voltammograms were normalized with respect to both lithium reduction and chloride ion oxidation potentials. Such a normalization technique is effective for comparing experimentally obtained cyclic voltammetry data to those that have been published. An analysis of experimentally obtained results indicates the relation of electrode reactions (e.g., differences of reduction potentials) are independent of the choice of reference electrode. Additionally, the choice of reference electrode did not affect the electrochemical window, exchange current density, or diffusion coefficient values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗