SEARCH · Search NASA

Results for “diffusion processes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Two reference time scales for studying the dynamic cavitation of liquid films

Two formulas, one for the characteristic time of filling a void with the vapor of the surrounding liquid, and one of filling the void by diffusion of the dissolved gas in the liquid, are derived. By comparing these time scales with that of the dynamic operation of oil film bearings, it is concluded that the evaporation process is usually fast enough to fill the cavitation bubble with oil vapor; whereas the diffusion process is much too slow for the dissolved air to liberate itself and enter the cavitation bubble. These results imply that the formation of a two phase fluid in dynamically loaded bearings, as often reported in the literature, is caused by air entrainment. They further indicate a way to simplify the treatment of the dynamic problem of bubble evolution.

Sun, D. C.↗

NASA NTRS1975

The variation of solar proton energy spectra and size distribution with heliolongitude

A statistical study of the initial phases of a large number of solar particle events provides an approach for an investigation of the variability in the production process and coronal diffusion. A description is given of the procedures which are suitable for identifying the parent flare of an event, taking into account the heliocentric distribution of the flare-associated events. The determination of the proton energy spectrum at the sun is discussed along with the variation of the spectral index with the associated flare heliocentric longitude, questions concerning the source spectra and the acceleration process, and the size distribution of flare associated particle events.

Van Hollebeke, M. A. I.↗

DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗

DOE OSTI2023

Rubidium and potassium isotopic variations in chondrites and Mars: Accretion signatures and planetary overprints

As moderately volatile elements, isotopes of Rb and K can trace volatilization processes in planetary bodies. Rubidium isotopic data are however very scarce, especially for non-carbonaceous meteorites. Here, in this study, we report combined Rb and K isotopic data (δ 87/85 Rb and δ 41/39 Κ) for 7 ordinary, 6 enstatite, and 4 Martian meteorite falls to understand the causes for the variations in volatile abundances and isotopic compositions. Bulk Rb and K isotopic compositions of planetary bodies are estimated to be (Table 1): Mars +0.10 ± 0.03 ‰ for Rb and -0.26 ± 0.05 ‰ for K, bulk OCs $-0.12^{+0.15}_{-0.24}$ ‰ for Rb and $-0.72^{+0.28}_{-0.41}$ ‰ for K, bulk ECs $-0.02^{+0.29}_{-0.26}$ ‰ for Rb and $-0.33^{+0.37}_{-0.23}$ ‰ for K. The bulk K isotopic compositions of subgroup OCs are estimated to be $-0.72^{+0.26}_{-0.55}$ ‰ for H chondrites, $-0.71^{+0.23}_{-0.39}$ ‰ for L chondrites, and $-0.77^{+0.63}_{-0.30}$ ‰ for LL chondrites. A broad correlation between the Rb and K isotopic compositions of planetary bodies is observed. The correlation follows a slope that is consistent with kinetic evaporation and condensation processes, suggesting volatility-controlled mass-dependent isotope fractionation (as opposed to nucleosynthetic anomalies). Individual ordinary and enstatite chondrites show large Rb and K isotopic variations (-1.02 to +0.29 ‰ for Rb and -0.91 to -0.15 ‰ for K). Samples of lower metamorphic grades display correlated elemental and isotopic fractionations between Rb and K, while samples of higher metamorphic grades show great scatter, suggesting that chondrite parent-body processes have decoupled the two elements and their isotopes at the sample scale. Several processes could have contributed to the observed isotopic variations of Rb and K, including (i) chondrule “nugget effect”, (ii) volatilization during parent-body thermal metamorphism (heat-induced vaporization and gas transport within parent bodies), (iii) thermal diffusion during parent-body metamorphism, and (iv) impact/shock heating. Quantitative modeling of the first two processes suggests that neither of them could produce isotopic variations large enough to explain the observed isotopic variations. Volatilization during parent-body thermal metamorphism [the scenario (ii)], which has been commonly invoked to explain the isotopic variations of volatile elements, is gas transport-limited and its effect on isotopic fractionations of moderately volatile elements should be negligible. Modeling of diffusion processes suggests that (iii) could produce K isotopic variation comparable to the observed variation. The large isotopic variations in non-carbonaceous meteorites are thus most likely due to diffusive redistribution of K and Rb during metamorphism and/or shock-induced heating and vaporization.

58 GEOSCIENCES↗

DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗

DOE OSTI2024

Lamination of >21% Efficient Perovskite Solar Cells with Independent Process Control of Transport Layers and Interfaces

Transport layer and interface optimization is critical for improving the performance and stability of perovskite solar cells (PSCs) but is restricted by the conventional fabrication approach of sequential layer deposition. While the bottom transport layer is processed with minimum constraints, the narrow thermal and chemical stability window of the halide perovskite (HP) layer severely restricts the choice of top transport layer and its processing conditions. To overcome these limitations, we demonstrate lamination of HPs—where two transport layer-perovskite half-stacks are independently processed and diffusion-bonded at the HP-HP interface—as an alternative fabrication strategy that enables self-encapsulated solar cells. Power conversion efficiencies (PCE) of >21% are realized using cells that incorporate a novel transport layer combination along with dual-interface passivation via self-assembled monolayers, both of which are uniquely enabled by the lamination approach. This is the highest reported PCE for any laminated PSC encapsulated between glass substrates. We further show that this approach expands the processing window beyond traditional fabrication processes and is adaptable for different transport layer compositions. The laminated PSCs retained >75% of their initial PCE after 1000 h of 1-sun illumination at 40 °C in air using an all-inorganic transport layer configuration without additional encapsulation. Furthermore, a laminated 1 cm 2 device maintained a V oc of 1.16 V. Finally, the scalable lamination strategy in this study enables the implementation of new transport layers and interfacial engineering approaches for improving performance and stability.

14 SOLAR ENERGY↗

DOE OSTI2018

Reactivity of heterogeneous surfaces: Modeling precursor–surface interaction using absorbing Markov chains

Heterogeneous reactive surfaces involving a spatial distribution of reactive sites are ubiquitous in thin film growth, from defect nucleation to step flow growth. In this work, the interaction of gas phase species with heterogeneous reactive surfaces was explored using a model combining absorbing Markov chains and local balance equations. Through this model, single particle statistics can be extracted of the adsorption-diffusion-desorption/reaction process even beyond the transport-limited regime, where the surface coverage of adsorbates, and hence adsorbate-adsorbate interaction, is not negligible, and local rates can be connected with macroscopic quantities, such as the sticking probability and desorption rates. In the transport-limited regime, a similarity law was identified for the sticking probability for both step flow growth and nucleation on randomly distributed reactive sites in which the sticking probability is a function of the product of the local diffusion and reaction probabilities. Lastly, the model can be applied to extract local reactivity maps providing the probability that a particle adsorbing on a given site reacts with the surface. This allows one to explore the development of capture zones and the impact of site-blocking and desorption on the local sticking probabilities.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1992

Chemical vapor deposition fluid flow simulation modelling tool

Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.

Bullister, Edward T.↗

DOE OSTI2020

Understanding chemical and physical mechanisms in atomic layer deposition

Atomic layer deposition (ALD) is a powerful tool for achieving atomic level control in the deposition of thin films. However, several physical and chemical phenomena can occur which cause deviation from “ideal” film growth during ALD. Understanding the underlying mechanisms that cause these deviations is important to achieving even better control over the growth of the deposited material. Herein, we review several precursor chemisorption mechanisms and the effect of chemisorption on ALD growth. We then follow with a discussion on diffusion and its impact on film growth during ALD. Together, these two fundamental processes of chemisorption and diffusion underlie the majority of mechanisms which contribute to material growth during a given ALD process, and the recognition of their role allows for more rational design of ALD parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Collective Summary of sCO2 Materials Development Part II: High-Temperature Alloy Corrosion Behavior within Compact Heat Exchangers.

Materials compatibility is a major concern for the design and operation of supercritical carbon dioxide (sCO2) power cycles. Two areas were recently identified for which very little prior knowledge was available. These were the behavior of polymer materials in low-temperature sCO2 environments as well as the high-temperature corrosion behavior for alloys that have been fabricated into compact heat exchangers. A critical need exists to increase understanding in both areas. As such, the Sandia sCO2 Materials Development program initiated a series of experiments for both areas in FY19. The progress that has been made in understanding the behavior of polymers in sCO2 was summarized in the Part I report, while this Part II report describes progress in the area of compact heat exchanger corrosion. For the compact heat exchangers being developed for use in supercritical CO2 power cycles, alloy corrosion is important to understand as it may lead to reduced flow path area, reduced heat transfer, and blocked flow channels. The fabrication of these heat exchangers involve thermal processes such as brazing or diffusion bonding, where the influence of these processes on the corrosion behavior of the alloys is unknown. Candidate alloys for these heat exchangers are typically evaluated for corrosion behavior as witness coupons, without being subjected to these thermal processes. To close this gap in understanding, a series of experiments were completed that utilized small heat exchanger samples subjected to flowing CO2 for 500 hours at 750°C. Small sections of the heat exchanger samples were characterized before and after the 500 hour exposure in order to characterize the oxide growth in the channels as well as the hardness of the alloys. Results are compared to witness coupons of the same alloys.

36 MATERIALS SCIENCE↗

NASA NTRS1992

Convective effects during the physical vapor transport process. I - Thermal convection

The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra(T) ranging from 1.80 to 1.92 x 10 exp 6. Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is 10 increase the threshold of instability. For typical ground based conditions time dependent oscillatory convection can occur. This results in nonuniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows regions of both periodic and quasi-periodic states. Low gravity conditions can effectively reduce convective effects, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Duval, Walter M. B.↗

DOE OSTI2023

Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation

The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? For this work, we investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and timescales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and timescales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Explosive nucleation and growth of Pb islands on Ge(111) below room temperature via collective diffusion

Nucleation is a fundamental process in nature controlling phase transitions and pattern formation. Classical nucleation is based on the gradual aggregation of diffusing atoms after a critical cluster size is reached. Contrary to these expectations, we have observed in real time the formation of perfect Pb(111) islands with more than 10 5 atoms emerging out of the compressed wetting layer within a few seconds during growth of Pb on Ge(111) below room temperature. Pb deposited on Ge(111) in a temperature range between −30 °C and 10 °C exhibited explosive nucleation of height-selected islands upon reaching a critical coverage of 1.33 ± 0.07 ML with respect to Ge(111). Island nucleation was fueled by the 2% compression of the Pb wetting layer at the critical coverage. The island areas grew linearly with time, exhibiting collective diffusion, where thousands of atoms follow correlated, non-random walk diffusion. The total growth rate was higher with increasing temperature. Density functional theory simulations of the chemical potential and binding sites of the Pb/Ge(111) system give insight into the role of the compression of the wetting layer and the strain within the growing Pb layers, which explain this nonclassical behavior of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Morphogenic Growth 3D Printing

Inspired by nature's morphogenesis, a new 3D printing process –growth printing (GP)– takes advantage of a self‐propagating curing front to produce 3D polymeric parts following a growth‐like development plan. The propagation of the curing front is driven by the exothermic polymerization of dicyclopentadiene (DCPD), which transforms the liquid resin into a stiff polymer as it propagates at 1 mm s −1 . GP is triggered when a heated initiator contacts the uncured liquid resin in an open container. The initiator nucleates the frontal polymerization reaction and the isotropic radial propagation of the growth front. Simultaneously, the initiator is moved up across the free surface of the resin, pulling the cured object out of the uncured resin. The motion trajectory of the initiator with respect to the free resin surface controls the growth morphology of the 3D part. An inverse design algorithm is developed to produce 3D parts by modeling the reaction‐diffusion‐driven solidification process. This process has substantial energy savings and high printing speeds.

3D printing↗

DOE OSTI2024

Tunable Growth of Layered Double Hydroxide Nanosheets through Hydrothermal Conversion of Atomic Layer Deposition Seed Layers

To enable the design and manufacturing of hierarchical nanomaterial architectures, there is a need for synthesis and processing methods that can enable tunable geometric control at the nanoscale while maintaining conformality on complex 3-D templates. Here, in this study, we explore the programmable control of vertically oriented Zn-Al layered double hydroxide (LDH) nanosheet arrays using atomic layer deposition (ALD) to deposit a seed layer of Al 2 O 3 , which is subsequently consumed and converted into the LDH phase under hydrothermal growth conditions. We demonstrate tunable control over the spacing and length of the nanosheets by varying the thickness of the initial ALD seed layer with subnanometer precision. This can be viewed as a nanoscale titration reaction, where Al acts as the limiting reagent during the hydrothermal synthesis of the nanosheets. Elemental mapping demonstrates the dynamic evolution of the resulting morphology, which is driven by surface diffusion and nucleation processes. The conformal nature of ALD allows for hierarchical growth of nanosheets on the surface of a variety of nonplanar substrate geometries, including microposts, paper fibers, and porous ceramic supports. This illustrates the power of ALD to enable bottom-up growth of 3-D nanoarchitectures with tunable geometries by controlling nucleation and growth in subsequent solution reactions.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Fractography of the high temperature hydrogen attack of a medium carbon steel

Microscopic fracture processes were studied which are associated with hydrogen attack of a medium carbon steel in a well-controlled, high-temperature, high-purity hydrogen environment. Exposure to a hydrogen pressure and temperature of 3.5 MN/m2 and 575 C was found to degrade room temperature tensile properties with increasing exposure time. After 408 hr, yield and ultimate strengths were reduced by more than 40 percent and elongation was reduced to less than 2 percent. Initial fissure formation was found to be associated with manganese rich particles, most probably manganese oxide, aligned in the microstructure during the rolling operation. Fissure growth was found to be associated with a reduction in carbide content of the microstructure and was inhibited by the depletion of carbon. The interior surfaces of sectioned fissures or bubbles exhibit both primary and secondary cracking by intergranular separation. The grain surfaces were rough and rounded, suggesting a diffusion-associated separation process. Specimens that failed at room temperature after exposure to hydrogen were found to exhibit mixed mode fracture having varying amounts of intergranular separation, dimple formation, and cleavage, depending on exposure time.

Melson, H. G.↗

DOE OSTI2019

Chemical explosive mode analysis of lean blowout in a gas turbine combustor

In this study, stabilization and lean blowout (LBO) mechanisms in a gas turbine combustor are numerically studied based on large eddy simulation (LES) data and the chemical explosive mode analysis (CEMA). Three CEMA-based local combustion modes, that is, auto-ignition, diffusion-assisted ignition, and extinction, are employed to quantify interactions between local diffusion and chemical processes in complex turbulent fields. The capability of CEMA-based local modes in extinction detection is first demonstrated in a one-dimensional non-premixed flamelet undergoing transient extinction. Flame structure and dynamics in the gas turbine combustor are then examined in detail using CEMA and the local combustion modes for both the stably burning and lean blowout conditions. Local structure comprised of all the three modes are identified, and each mode is found to play an important role in determining the overall burning rate. The heat release rate attributed to each local mode remains nearly constant at the stably burning condition. In contrast, the flame during LBO shows a two-stage extinction behavior: at the initial stage, the extinction mode is largely enhanced and becomes dominant over the assisted-ignition model; at the second stage, all three modes are suppressed, signifying global extinction. It is therefore suggested that the excessive local extinction mode, due to an increased amount of heat and radical losses through diffusion, finally leads to lean flame blowout.

chemical explosive mode↗