Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Opposed jet burner studies of effects of CO, CO2, and N2 air-contaminants on hydrogen-air diffusion flames

The blowoff/restore characteristics for jets of various H2/N2 mixtures opposed to jets of air contaminated by N2, CO, and CO2 have been determined using a counterflow diffusion flame formed by a tubular opposed jet burner. Both blowoff and restore limits are found to be sensitive to fuel and air composition. Empirically derived variations in the limits of the average mass flux of incoming H2 with percent contaminant, at fixed incoming fuel and H2/O2 inputs, are used to quantify the effects of oxygen dilution, flame augmentation, and flame retardation by N2, CO, and CO2 contaminants. The implications of the results are discussed.

Guerra, Rosemary↗
DOE OSTI2023

Diffusion, atomic transport, and ordering in Al-Zr alloys: FCC and liquid phases

Additive manufacturing of materials with controlled microstructure demands knowledge of atomic scale properties near the solid-liquid transition state. Many of these properties are not affordable by experimental techniques and computer modeling is the possible solution to the problem. In this paper, we present the results of an extended atomistic study of intrinsic atomic transport due to vacancy diffusion in FCC and L12 solid phases and diffusion in the liquid phase of Al-Zr alloys. A deceleration of the overall self-diffusion was observed when Zr was added to Al. The effect was stronger in the solid and weaker in the liquid. Additionally, the effect was strongly temperature dependent in the solid phases, but not in the liquid. Atomic transport was chemically biased: transport of Zr atoms was significantly slower than that of Al atoms, and this bias effect was stronger in the solid phases. The overall diffusion and chemical ordering processes in the liquid state were five to six orders in magnitude faster than in the solid. Chemical short-range order parameters in the liquid saturated at values close to those in the ordered L12 structure of Al3Zr. Chemical and structural ordering in the solid phases was negligible over the modeled microsecond time scale. Here, the results are discussed in view of optimizing additive manufacturing parameters for the controlled formation of metastable L1 2 precipitates.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Nonneutral and quasi-neutral diffusion of weakly ionized multiconstituent plasma

A theory is developed for quasi-neutral and nonneutral ambipolar diffusion in a weakly ionized multiconstituent plasma. Magnetic field effects are disregarded. The effects of negative ions and various positive and negative ion diffusion coefficients are identified. It is found that the presence of negative ions leads to enhanced electron diffusion and a retarded diffusion of positive ions. In some cases electrons may be forcibly propelled from the area of perturbations in ionic concentrations. The ambipolar electric field is observed to cause negative ions to flow in a countergradient fashion where the gradient is small. Electrons tend toward a more homogeneous spatial distribution in the nonneutral case than in the quasi-neutral case.

Hill, R. J.↗
DOE OSTI2020

Hybridizing pseudo-Hamiltonians and non-local pseudopotentials in diffusion Monte Carlo

An accurate treatment of effective core potentials (ECPs) requires care in continuum quantum Monte Carlo (QMC) methods. While most QMC studies have settled on the use of familiar non-local (NL) pseudopotentials with additional localization approximations, these approaches have been shown to result in moderate residual errors for some classes of molecular and solid state applications. Here, we revisit an idea proposed early in the history of QMC ECPs that does not require localization approximations, namely, a differential class of potentials referred to as pseudo-Hamiltonians. We propose to hybridize NL potentials and pseudo-Hamiltonians to reduce residual non-locality of existing potentials. We derive an approach to recast pseudopotentials for 3d elements as hybrid pseudo-Hamiltonians with optimally reduced NL energy. We demonstrate the fidelity of the hybrid potentials by studying atomic ionization potentials of Ti and Fe and the binding properties of TiO and FeO molecules with diffusion Monte Carlo (DMC). We show that localization errors have been reduced relative to potentials with the same NL channels for Sc–Zn by considering the DMC energy change with respect to the choice of approximate localization. While localization error decreases proportionate to the reduced NL energy without a Jastrow, with a Jastrow, the degree of reduction decreases at higher filling of the d-shell. Our results suggest that a subset of existing ECPs may be recast in this hybrid form to reduce the DMC localization error. They also point to the prospect of further reducing this error by generating ECPs within this hybrid form from the start.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Wave Driven Exothermic Heating in the Mesopause Region

A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Hickey, Michael P.↗
NASA NTRS1982

Gas-phase hydrogen permeation through alpha-titanium - Surface film and dimensional effects

The process of hydrogen transport through alpha-Ti involves simultaneous diffusion and phase boundary reactions at both surfaces, with the relative effect on each surface depending on the extent of surface contamination as well as the physical dimensions of the titanium membrane used. It is shown by the present study that hydrogen permeation in commercially pure alpha-Ti increases exponentially with temperature and is dependent on the first power of the input pressure, whether the surface is as-polished, preoxidized or prenitrided. Permeation decreases in the case of the as-polished condition if nitride or oxide films are formed at the surface in contact with source hydrogen, while increasing slightly for the same condition if such films are formed at the hydrogen exit surface.

Shah, K. K.↗
NASA NTRS1993

Effect of Long-Term 1093 K Exposure to Air or Vacuum on the Structure of Several Wrought Superalloys

Long-term 1093 K heat treatments of three commercial superalloy sheet materials were undertaken in air -IAI and vacuum. With either exposure, significant precipitation of second phases occurred in the Co-base Haynes(R) Alloy 188 (HA 188) and the Ni-base Haynes(R) Alloy 230 (HA 230); however, much less precipitation was found in the exposed Ni-base alloy Inconel(R) 617 (IN 617). Although some grain growth occurred in HA 198, no changes in the grain size of either HA 230 or IN 617 were observed after 22,500 h at temperature. Oxidation during air heat treatments led to weight gain due to the formation of chromia + spinel scales and surface-connected grain boundary pits/oxides in all three superalloys. Both the weight gain and depth of intergranular attack were dependent on the square root of time, which is indicative of diffusion-controlled phenomena. Because many alloy samples had neighbors in close proximity, most vacuum heat treated specimens did not suffer significant loss of volatile elements. However, some exposed samples were subjected to unrestricted vacuum heat treatments, allowing estimates of volatilization to be made. Based on the data for HA 188, the weight loss during 1093 K vacuum exposure was diffusion controlled once the inhibiting effects of surface films on the as-received alloys were broken down.

J.D. Whittenberger↗
NASA NTRS1981

Double-beam optical method and apparatus for measuring thermal diffusivity and other molecular dynamic processes in utilizing the transient thermal lens effect

A sample material was irradiated by relatively high power, short pulses from a dye laser. Energy from the pulses was absorbed by the sample material, thereby forming a thermal lens in the area of absorption. The pulse repetition rate was chosen so that the thermal lens is substantially dissipated by the time the next pulse reaches the sample material. A probe light beam, which in a specific embodiment is a relatively low power, continuous wave (cw) laser beam, irradiated the thermal lens formed in the sample material. The intensity characteristics of the probe light beam subsequent to irradiation of the thermal lens is related to changes in the refractive index of the sample material as the thermal lens is formed and dissipated.

Moacanin, J.↗
NASA NTRS2018

Effect of the Reconnection Electric Field on Electron Distribution Functions in the Diffusion Region of Magnetotail Reconnection

Electron distribution functions in the electron diffusion region during symmetric magnetic reconnection are investigated by means of theory and fully kinetic simulations. Crescent‐like striations are formed in distribution functions in the velocity plane perpendicular to the magnetic field. Using an analytical current sheet, we solve the equation of motion for electrons, and derive the shape of a crescent distribution, as a function of the distance from the neutral line, field gradients, and the reconnection electric field. Each crescent is tilted in the velocity plane because of the acceleration by the reconnection electric field, and multiple stripes appear due to multiple meandering bounces. Applying the theory to distribution functions observed in Earth's magnetotail, we deduce the amplitude of the reconnection electric field.

Bessho, N.↗
DOE OSTI2023

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Short-range thermal magnon diffusion in magnetic garnet

Here, using the spin Seebeck effect (SSE), we study the propagation distance of thermally induced spin currents inside a magnetic insulator thin film in the short-range regime. We disambiguate spin currents driven by temperature and chemical potential gradients by comparing the SSE signal before and after adding a heat-sinking capping layer on the same device. We report that the measured spin decay behavior near the heat source is well accounted for by a diffusion model where the magnon diffusion length is in submicron range, in other words, two orders of magnitude smaller than previous estimates inferred from the long-range behavior. Our results highlight the caveat in applying a diffusive theory to describe thermally generated magnon transport, where a single decay length may not capture the behavior on all length scales.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Diffusive-convective physical vapor transport of PbTe from a Te-rich solid source

Crystal growth of PbTe by physical vapor transport (sublimation) in a closed ampoule is governed by the vapor species in thermal equilibrium with the solid compound. Deviations from stoichiometry in the source material cause diffusion limitation of the transport rate, which can be modified by natural (gravity-driven) convection. Mass-transport experiments have been performed using Te-rich material wherein sublimation rates have been measured in order to study the effects of natural convection in diffusion-limited vapor transport. Linear velocities for both crystal growth and evaporation (back sublimation) have been measured for transport in the direction of gravity, horizontally, and opposite to gravity. The experimental results are discussed in terms of both the one-dimensional diffusive-advective model and current, more sophisticated theory which includes natural convection. There is some evidence that convection effects from radial temperature gradients and solutal density gradients have been observed.

Zoutendyk, J.↗
NASA NTRS1992

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗
NASA NTRS1993

Diffusion of Na and K in the uppermost regolith of Mercury

We examine the conditions under which ions impacting the Hermean surface can act as a regional source of enhanced atmospheric column through ion implantation and subsequent release at the surface. At most latitudes, energetic (a few keV) Na or K ions which impact the nightside surface, are released quickly (well before noon) upon warming of the surface. We show that the relative sunrise/sunset difference produced by ion implantation is ((n)(sub SR) - (n)(sub SS)/(n)(sub Ave) = f, where f is the fraction of the photo-ions recycled, (n)(sub SR) is the average zenith column above the sunrise portion of the illuminated hemisphere, (n)(sub SS) is the average zenith column in the sunset portion, and (n)(sub Ave) is the average column over the sunlit disk. Thus, to produce a large sunrise/sunset difference via ion implantation and subsequent release requires efficient (close to total) recycling. We show that the most extensive set of available data reduced to Na column abundance does not show any sunrise enhancement. We argue that the K data do not permit an unambiguous interpretation in favor of sunrise/sunset differences. We further find that if an efficient surface loss process for the alkali is not operating after sunrise, the initial Na and K distribution will relax into the bulk of the solid. Preserved abundance gradients of Na and K in lunar glasses that the Arrhenius coefficients for impact glasses are likely more modest than those for laboratory glasses derived from rock by nonimpact processes, but the effect of the more modest diffusion rates is only to delay the efficient loss of Na by a few Earth days. We argue that implantation can lead to observable regional increases in the observed Na or K column densities only if it occurs at very high latitudes, where diffusion is slow. It is typically lower energy ions which impact at high latitude and these are both more numerous than the high energy ions and possess smaller average penetration depths; thus there are additional reasons to favor a high-latitude locus for any possible prompt return related increases in zenith column. We find that we cannot rule out sputtering as a source process, as the low-energy ions are efficient sputterers. Finally, we argue that the observed Na/K ratio in the atmopshere may be the expression of their different asymptotic rates of loss from the interiors of the regolith grains. If so, this allows us to fix the importance of sputtering and photon stimulated desorption relative to impact vaporization.

Killen, R. M.↗
NASA NTRS2004

Downscattering due to Wind Outflows in Compact X-ray Sources: Theory and Interpretation

A number of recent lines of evidence point towards the presence of hot, outflowing plasma from the central regions of compact Galactic and extragalactic X-ray sources. Additionally, it has long been noted that many of these sources exhibit an "excess" continuum component, above approx. 10 keV, usually attributed to Compton Reflection from a static medium. Motivated by these facts, as well as by recent observational constraints on the Compton reflection models - specifically apparently discrepant variability timescales for line and continuum components in some cases - we consider possible of effects of out-flowing plasma on the high-energy continuum spectra of accretion powered compact objects. We present a general formulation for photon downscattering diffusion which includes recoil and Comptonization effects due to divergence of the flow. We then develop an analytical theory for the spectral formation in such systems that allows us to derive formulae for the emergent spectrum. Finally we perform the analytical model fitting on several Galactic X-ray binaries. Objects which have been modeled with high-covering-fraction Compton reflectors, such as GS1353-64 are included in our analysis. In addition, Cyg X-3, is which is widely believed to be characterized by dense circumstellar winds with temperature of order 10(exp 6) K, provides an interesting test case. Data from INTEGRAL and RXTE covering the approx. 3 - 300 keV range are used in our analysis. We further consider the possibility that the widely noted distortion of the power-law continuum above 10 keV may in some cases be explained by these spectral softening effects.

Titarchuk, Lev↗