Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion process”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Pore-scale modelling of Ostwald ripening

In a saturated solution with dispersed clusters of a second phase, the mechanism by which the larger clusters grow at the expense of the smaller ones is called Ostwald ripening. Although the mechanism is well understood in situations where multiple clusters of gas exist in a liquid solution, evolution is much more complicated to predict when the two phases interact with a solid matrix, since the solid plays a significant role in determining the shape of the interface. In this paper, we study capillary dominated regimes in porous media where the driving force is inter-cluster diffusion. By decomposing the Ostwald ripening mechanism into two processes that operate on different time scales – the diffusion of solute gas in the liquid and the readjustment of the shape of the gas–liquid interface to accommodate a change in mass – we develop a sequential algorithm to solve for the evolution of systems with multiple gas ganglia. In the absence of a solid matrix, thermodynamic equilibrium is reached when all of the gas phase aggregates to form one large bubble. In porous media on the other hand, we find that ripening can lead to equilibrium situations with multiple disconnected ganglia, and that evolution is highly dependent on initial conditions and the structure of the solid matrix. The fundamental difference between the two cases is in the relationship between mass and capillary pressure.

bio-inspired, mechanical behavior, carbon sequestr↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2021

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗
DOE OSTI2020

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗
DOE OSTI2022

Ultrasonic nondestructive evaluation of additively manufactured wear coatings

Here, this work pertains to the nondestructive evaluation of additively deposited coatings using ultrasound measurements. The specific objective was to evaluate using ultrasound surface wave measurements the quality of Stellite 21 coatings deposited on Inconel 718 substrates by the directed energy deposition (DED) additive manufacturing (AM) process. The surface wave speed and diffuse backscatter amplitude of the ultrasound waves were correlated with both cracking and warpage of the coating. This research is important, because, the integrity of DED-processed coatings currently requires a combination of destructive metallographic studies and X-ray computed tomography analysis, which are both expensive and time consuming. Instead of using a normal incidence configuration for the ultrasound measurements, the surface-wave approach allows the inspection of the near surface coating integrity. Three different frequencies were used to excite surface waves because each frequency has a different penetration depth. The signals obtained were used to quantify and compare the surface wave speed and surface wave-diffuse ultrasonic backscatter amplitude from five different DED-processed Stellite coated samples. Considering destructive metallurgical characterization as a reference, surface wave measurements were found to be effective for examining coating integrity. Further, the backscatter amplitude of the surface waves was correlated with the crack density and warpage. These flaws also change the stress state of the Stellite coating and consequently the surface wave propagation speed and its scattering behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Diffusion coefficients in the envelopes of white dwarfs

The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs and their spectral evolution. The diffusion coefficients of Paquette et al. have been widely used to model diffusion in white dwarfs. In this work, we perform new calculations of the coefficients of interdiffusion and ionic thermal diffusion with (1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate, and moderately coupled plasma and (2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of interdiffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars, as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (${T}_{\mathrm{eff}}$ < 15,000 K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion timescales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of ≳3 differences between calculations based on Paquette et al. and our model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Stochastic Acceleration of ~0.1–5 keV Pickup Ions in the Heliotail

We seek to understand the quantitative role of the dominant physical processes (charge-exchange, adiabatic heating, stochastic acceleration) governing the proton distribution in the heliotail using observations of hydrogen energetic neutral atoms (ENAs) from the Interstellar Boundary Explorer ( IBEX ). We solve the Parker transport equation for solar wind protons and pickup ions (PUIs) as they propagate from the termination shock (TS) down the heliotail, including charge-exchange between protons and neutral hydrogen atoms as source terms derived from an MHD-fluid and kinetic-neutral simulation of the heliosphere. We compute ENA fluxes at 1 au from the results of the proton transport model and compare them with IBEX observations. We find that, under the assumptions of our model, a stochastic acceleration process is needed to counteract the energy-dependent losses of ~0.1–5 keV PUIs from charge-exchange to reproduce IBEX data. The power-law velocity dependence of the diffusion coefficient (spectral index γ ) is limited to the range 0.67 < γ < 2, and the best fit to IBEX data appears close to γ ~ 1.25. The diffusion rate ~1.1 &#x00d7; 10 -8 km 2 s -3 ( v / v 0 ) 1.25 nearly balances the loss of ~0.1–5 keV PUIs by charge-exchange. Our analysis suggests that cyclotron resonance with two widely known incompressible MHD turbulence: namely, isotropic Kolmogorov and anisotropic Goldreich–Sridhar turbulence, as well as stochastic particle interactions with compressive waves are not by themselves the dominant diffusion mechanisms. However, some intermediate processes may be occurring due to the presence of PUIs.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Influence of Tabula Rasa on Process- and Light-Induced Degradation of Solar Cells Fabricated From Czochralski Silicon

Monocrystalline Si solar cells are fabricated from Czochralski (Cz) Si, which contains 10 17 –10 18 cm -3 oxygen atoms. Cz Si undergoes degradation during high-temperature thermal processing steps, such as dopant diffusion to form the p - n junction. This degradation in the bulk minority carrier lifetime can be related to the formation of oxygen precipitates. We found that a high-temperature annealing process known as tabula rasa (TR) not only mitigates process-induced degradation via oxygen precipitate nuclei dissolution, but also modifies subsequent light-induced degradation. We report on the bulk carrier lifetime of n - and p -type Cz Si after TR, which homogenizes the interstitial oxygen in the bulk Si to its monoatomic form in either an N 2 or O 2 environment. A control sample, which was not subjected to a TR processing step, experienced severe process-induced degradation during a boron emitter thermal budget as oxygen precipitates were formed in the Si bulk. These precipitates could be imaged using photoluminescence. Additionally, samples that underwent a TR processing step prior to the boron emitter thermal budget show efficient gettering of metallic impurities compared to the control sample, which showed a decline in the implied open-circuit voltage after the gettering step. Furthermore, modification of the interstitial oxygen bonding by TR had a strong effect on the light-induced degradation kinetics. Instead of a typically observed monotonic decrease, minority carrier lifetime increases first, followed by a nonmonotonic decrease over a ~20 h period. Finally, we conclude that by modifying the interstitial oxygen bonding via TR pretreatment, p -type Cz Si wafers become substantially resistant to harsh solar cell processes and strongly modified light-induced degradation, which would open alternative ways to mitigate this degradation mechanism.

14 SOLAR ENERGY↗
DOE OSTI2021

Formation of three-dimensional bicontinuous structures via molten salt dealloying studied in real-time by in situ synchrotron X-ray nano-tomography

Abstract Three-dimensional bicontinuous porous materials formed by dealloying contribute significantly to various applications including catalysis, sensor development and energy storage. This work studies a method of molten salt dealloying via real-time in situ synchrotron three-dimensional X-ray nano-tomography. Quantification of morphological parameters determined that long-range diffusion is the rate-determining step for the dealloying process. The subsequent coarsening rate was primarily surface diffusion controlled, with Rayleigh instability leading to ligament pinch-off and creating isolated bubbles in ligaments, while bulk diffusion leads to a slight densification. Chemical environments characterized by X-ray absorption near edge structure spectroscopic imaging show that molten salt dealloying prevents surface oxidation of the metal. In this work, gaining a fundamental mechanistic understanding of the molten salt dealloying process in forming porous structures provides a nontoxic, tunable dealloying technique and has important implications for molten salt corrosion processes, which is one of the major challenges in molten salt reactors and concentrated solar power plants.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Utilization of Tabula Rasa to Stabilize Bulk Lifetimes in n-Cz Silicon for High-Performance Solar Cell Processing

We investigate a high temperature, high cooling-rate anneal Tabula Rasa (TR) and report its implications on n-type Czochralski-grown silicon (n-Cz Si) for photovoltaic fabrication. Tabula Rasa aims at dissolving and homogenizing oxygen precipitate nuclei that can grow during the cell process steps and degrade the cell performance due to their high internal gettering and recombination activity. The Tabula Rasa thermal treatment is performed in a clean tube furnace with cooling rates >100 degrees C/s. We characterize the bulk lifetime by Sinton lifetime and photoluminescence mapping just after Tabula Rasa, and after the subsequent cell processing. After TR, the bulk lifetime surprisingly degrades to <; 0.1ms, only to recover to values equal or higher than the initial non-treated wafer (several ms), after typical high temperature cell process steps. Those include boron diffusion and oxidation; phosphorus diffusion/oxidation; ambient annealing at 850 degrees C; and crystallization annealing of tunneling-passivating contacts (doped polycrystalline silicon on 1.5 nm thermal oxide). The drastic lifetime improvement during high temperature cell processing is attributed to improved external gettering of metal impurities and annealing of intrinsic point defects. Time and injection dependent lifetime spectroscopy further reveals the mechanisms of lifetime improvement after Tabula Rasa treatment. Additionally, we report the efficacy of Tabula Rasa on n-type Cz-Si wafers and its dependence on oxygen concentration, correlated to position within the ingot.

14 SOLAR ENERGY↗
DOE OSTI2019

Theoretical Analysis of the EPR Spectrum Saturation Effect with Account for Spectral Diffusion in a System with a Gaussian Distribution of Spin Resonance Frequencies

We have analyzed the effect of saturation of the EPR spectrum for the model situation when a set of resonance frequencies of spin packets is described by a Gaussian distribution and spectral diffusion, which is a random process without correlation, exists. We have obtained the form of the spectrum for arbitrary values of magnetic–resonant parameters: relaxation times for the longitudinal and transverse components of the magnetization vector of the system, the spectral diffusion rate, the dispersion of the Gaussian distribution of resonance frequency, and the magnetic induction of the microwave field. We have derived an analytic expression for the saturation curve of the EPR spectrum. Detailed analysis of the dependence of the spectrum saturation curve on the spin system parameters is performed for the optimal case when the microwave field frequency coincides with the average frequency of the Gaussian distribution for resonance frequencies. The value of the microwave field induction at which the peak of the saturation curve is attained has been determined for arbitrary values of the paramagnetic relaxation time, spectral diffusion rate, and dispersion of the Gaussian spin frequency distribution. We have formulated the algorithm for determining the spin–lattice relaxation time using the experimentally determined microwave field induction for which the saturation curve passes through the peak. It should be noted that numerical calculations have been performed for parameters typical of EPR spectroscopy. The general analytic expression derived for the shape of the spectrum in the saturation conditions and qualitative conclusions concerning the effect of spectral diffusion on the saturation of stationary spectra obtained in this study can also be used in other fields of spectroscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Perona Malik anisotropic diffusion model using Peaceman Rachford scheme on digital radiographic image

In image processing, it is important to remove noise without affecting the image structure as well as preserving all the edges. Perona Malik Anisotropic Diffusion (PMAD) is a PDE-based model which is suitable for image denoising and edge detection problems. In this paper, the Peaceman Rachford scheme is applied on PMAD to remove unwanted noise as the scheme is efficient and unconditionally stable. The capability of the scheme to remove noise is evaluated on several digital radiography weld defect images computed using MATLAB R2009a. Experimental results obtained show that the Peaceman Rachford scheme improves the image quality substantially well based on the Peak Signal to Noise Ratio (PSNR). The Peaceman Rachford scheme used in solving the PMAD model successfully removes unwanted noise in digital radiographic image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Alpha

What is this parameter called alpha, and why is it so important? To answer these questions we consider the density of fission neutrons at any point $_r^→$ = (x, y, z) inside a volume of fissionable material. Ignoring neutron absorption, this time-varying density of neutrons is determined by two physical processes --- generation of more neutrons by fission chain reactions, and the loss of neutrons by diffusion to and through the surface. A simple, approximate mathematical description of these processes is the diffusion equation for the average neutron density n($_r^→$, t) including a fission source term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Interfacial solvation and slow transport of hydrated excess protons in non-ionic reverse micelles

This work employs molecular dynamics simulations to investigate the solvation and transport properties of hydrated excess protons (with a hydronium-like core structure) in non-ionic Igepal CO-520 reverse micelles of various sizes in a non-polar solvent. Multiscale Reactive Molecular Dynamics (MS-RMD) simulations were used to describe vehicular and hopping diffusion during the proton transport process. As detailed herein, an excess proton shows a marked tendency to localize in the interfacial region of micellar water pools. Slow proton transport was observed which becomes faster with increasing micellar size. Further analysis reveals that the slow diffusion of an excess proton is a combined result of slow water diffusion and the low proton hopping rate. This study also confirms that a low proton hopping rate in reverse micelles stems from the interfacial solvation of hydrated excess protons and the immobilization of interfacial water. The low water density in the interfacial region makes it difficult to form a complete hydrogen bond network near the hydrated excess proton, and therefore locks in the orientation of hydrated proton cations. The immobilization of the interfacial water also slows the relaxation of the overall hydrogen bond network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Facilitated transport membranes for H 2 purification from coal-derived syngas: A techno-economic analysis

A single-stage membrane process has been designed for using facilitated transport membranes (FTMs) to decarbonize the coal-derived syngas from an integrated gasification combined cycle (IGCC) power plant. The necessary process model and costing method have also been developed to assess the technical feasibility and process economics. In order to account for the carrier saturation phenomenon associated with FTMs, a homogeneous reactive diffusion model is integrated into the process model. The techno-economic study reveals that the mitigated carrier saturation upon bulk CO 2 removal can lead to appreciable increases in the CO 2 permeance and CO 2 /H 2 selectivity, which can be utilized to achieve 95% CO 2 purity and 95% H 2 recovery with a CO 2 /H 2 selectivity of 50 at the complete carrier saturation. FTMs with different facilitated transport characteristics can also be arranged in a hybrid membrane configuration to render a H 2 recovery of 99% and a cost of electricity of $\$118.5$ /MWh, which is 12.5% lower than that of the benchmark Selexol process.

42 ENGINEERING↗
DOE OSTI2019

First-principles study of hydrogen diffusion and self-clustering below tungsten surfaces

The diffusion and self-clustering nucleation behavior of hydrogen (H) without vacancies below tungsten (W) surfaces are important for understanding the retention of deuterium (D) in W crystals exposed to high-flux and low-energy D plasma. The H migration energy and binding energy of H to self-clusters near the W(100), W(110), and W(111) surfaces have been investigated by first-principles computer simulations using density functional theory. H diffusion from adsorption sites on the W(100), W(110), and W(111) surfaces into solute sites in the bulk requires energies of at least 1.21, 1.78, and 1.80 eV, respectively, while 0.27, 0.31, and 0.24 eV for the reverse process. In addition, the lateral diffusion of H between two subsurface layers below the W surfaces has been investigated. Two H atoms at a depth of 0.08 nm below the W(110) surface have the highest binding energy, followed by H atom pairs below the W(111) and W(110) surfaces. The nucleation and stability of H clusters depend on the surface orientation. A planar configuration between the first nearest neighbor {100} planes is energetically favorable for H self-clustering below the W(100), W(110), and W(111) surfaces. In conclusion, the thermal stability of a platelet containing 16 H atoms below the W surfaces at 300 and 600 K was also studied by ab initio molecular dynamics simulations, which indicate that the H platelet below a W(111) surface is more stable than that below either the W(100) or the W(110) surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

THE COSMIC BATTERY IN ASTROPHYSICAL ACCRETION DISKS

The aberrated radiation pressure at the inner edge of the accretion disk around an astrophysical black hole imparts a relative azimuthal velocity on the electrons with respect to the ions which gives rise to a ring electric current that generates large-scale poloidal magnetic field loops. This is the Cosmic Battery established by Contopoulos and Kazanas in 1998. In the present work we perform realistic numerical simulations of this important astrophysical mechanism in advection-dominated accretion flows, ADAFs. We confirm the original prediction that the inner parts of the loops are continuously advected toward the central black hole and contribute to the growth of the large-scale magnetic field, whereas the outer parts of the loops are continuously diffusing outward through the turbulent accretion flow. This process of inward advection of the axial field and outward diffusion of the return field proceeds all the way to equipartition, thus generating astrophysically significant magnetic fields on astrophysically relevant timescales. We confirm that there exists a critical value of the magnetic Prandtl number between unity and 10 in the outer disk above which the Cosmic Battery mechanism is suppressed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Vacuum ultra-violet damage and damage mitigation for plasma processing of highly porous organosilicate glass dielectrics

Porous organosilicate glass thin films, with k-value 2.0, were exposed to 147 nm vacuum ultra-violet (VUV) photons emitted in a Xenon capacitive coupled plasma discharge. Strong methyl bond depletion was observed, concomitant with a significant increase of the bulk dielectric constant. This indicates that, besides reactive radical diffusion, photons emitted during plasma processing do impede dielectric properties and therefore need to be tackled appropriately during patterning and integration. The detrimental effect of VUV irradiation can be partly suppressed by stuffing the low-k porous matrix with proper sacrificial polymers showing high VUV absorption together with good thermal and VUV stability. In addition, the choice of an appropriate hard-mask, showing high VUV absorption, can minimize VUV damage. Particular processing conditions allow to minimize the fluence of photons to the substrate and lead to negligible VUV damage. For patterned structures, in order to reduce VUV damage in the bulk and on feature sidewalls, the combination of both pore stuffing/material densification and absorbing hard-mask is recommended, and/or the use of low VUV-emitting plasma discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗