Search NASA⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Smoothed particle hydrodynamics modeling and analysis of oxide reduction process for uranium oxides

A common kinetic feature for oxide reduction chemical/electrochemical processes is oxygen transport via a porous metallic layer, which has been considered as a rate-determining step for reducing uranium oxides to metallic uranium. Accounting this kinetic behavior must involve the resolution of the moving reactive interface between shrinking oxide and expanding metal phases. This study presents a numerical model using smoothed particles hydrodynamics (SPH) to effectively deal with the evolution of the shrinking core reaction interface and oxygen transport via mass transfer of lithium oxide (Li 2 O) species in multiple mass transfer domains. We successfully validated the proposed model against a theoretical derivation for the oxide reduction process handling a shrinking oxide core with molten salt and metal ash medium on a simple planar geometry. Armed with successful validation results, the model examined a realistic reactant geometry to extend the arguments beyond the one-dimensional analyses, allowing the proposed model to apply to general application scenarios with multi-dimensional geometry. Here, this study demonstrated that the proposed model could simulate and evaluate an arbitrary reaction basket design without iterative experimental trials, which is prohibitive for a scaled high-temperature molten salt study in an inert environment. The construct potentially provides not only deep insights on multiphysics behaviors governing the process dynamics but also a robust framework for evaluating and screening candidate basket designs in the most cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Study of the Effect of Doping with Iron on the Luminescence of Zinc-Selenide Single Crystals

The effect of doping with iron (thermal diffusion from a surface) on the luminescence of zinc-selenide single crystals in the wavelength range 0.44–0.72 μm and on the spatial distribution of luminescence centers are studied. By means of two-photon confocal microscopy, planar and volume maps of edge (exciton) and impurity–defect luminescence in the above-indicated spectral range are obtained for both doped and undoped crystals. It is shown that crystal regions containing a high iron concentration exhibit low-intensity luminescence in this range. It is found that, in the process of diffusion, several types of impurity–defect centers distributed in a complex way within the crystal bulk are formed. The nature of these centers is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Virtual Volumetric Additive Manufacturing (VirtualVAM)

Abstract Tomographic volumetric additive manufacturing (VAM) produces arbitrary 3D geometries by exposure of a rotating volume of photopolymer resin to tomographically‐patterned illumination. This enables high speed, layer‐less printing of parts from a wide range of photopolymers not amenable to layer‐by‐layer processes. Since the entire geometry is produced at once over the course of a few seconds to minutes, molecular diffusion length scales become significant to the printing process. Understanding these molecular reaction and diffusion processes is imperative for advancing VAM to a usable technology. These processes are experimentally very difficult to monitor and measure. Herein, VirtualVAM ‐ a simulation framework for modeling the tomographic VAM process, is developed and experimentally validated. VirtualVAM simulates reaction, diffusion, and heat generation processes over the course of a print with single‐voxel resolution. From a few experimentally‐determined input parameters and a set of images for projection, VirtualVAM is able to generate a large spatio‐temporal data set for any given tomographic VAM print. Using VirtualVAM, a number of experimentally‐unattainable aspects of the VAM process are investigated such as single‐voxel conversion profiles, effect of molecular oxygen, and stopping time determination. VirtualVAM also enables the optimization of exposure patterns to further improve contrast between in‐part and out‐of‐part delivered dose.

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

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗
DOE OSTI2020

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Structures and diffusion of two-dimensional dusty plasmas on one-dimensional periodic substrates

Using numerical simulations, we examine the structure and diffusion of a two-dimensional dusty plasma (2DDP) in the presence of a one-dimensional periodic substrate (1DPS) as a function of increasing substrate strength. Both the pair correlation function perpendicular to the substrate modulation and the mean-squared displacement (MSD) of dust particles are calculated. It is found that both the structure and dynamics of 2DDP exhibit strong anisotropic effects, due to the applied 1DPS. As the substrate strength increases from zero, the structure order of dusty plasma along each potential well of 1DPS increases first probably due to the competition between the interparticle interactions and the particle-substrate interactions, and then decreases gradually, which may be due to the reduced dimensionality and the enhanced fluctuations. The obtained MSD along potential wells of 1DPS clearly shows three processes of diffusion in our studied 2DDP. Finally, between the initial ballistic and finally diffusive motion, there is the intermediate subdiffusion discovered here, which may result from the substrate-induced distortion of the caging dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Percolation transition in carbon composite on the basis of fullerenes and exfoliated graphite

The electrical conductivity of a carbon composite on the basis of C{sub 60} fullerenes and exfoliated graphite is investigated in the range of relative contents of components from 0 to 100%. The samples are obtained by the thermal treatment of the initial dispersed mixtures in vacuum in the diffusion–adsorption process and their further cold pressing. The resistivity of the samples gradually increases with an increase in the fraction of fullerenes, and a sharp transition from the conductive state to the dielectric one is observed after achieving certain concentrations of C{sub 60}. The interpretation of the results within the percolation theory makes it possible to evaluate the percolation threshold (expressed as a relative content of graphite) as equal to 4.45 wt % and the critical conductivity index as equal to 1.85 (which is typical for three-dimensional two component disordered media including those having pores).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Fracture Toughness of Thick Boride Layers Estimated by the Cross-Sectioned Scratch Test

New results about the fracture toughness (K{sub c}) of thick boride layers estimated by the cross-sectioned scratch test are presented in this study. The FeB-Fe{sub 2}B layers developed at the surface of borided AISI 1018 and AISI 1045 steels and the Fe{sub 2}B layer formed on the borided AISI 1045 steel exposed to a diffusion annealing process (DAP) were used for this purpose. The cross-sectioned scratch tests were performed with a Vickers diamond stylus drawn across the thick boride layer under a constant load to produce a half-cone-shaped fracture near to the top surface of the borided steels. The height of the half-cone-shaped fracture as a function of the cross-sectioned scratch loads was used to determine the fracture toughness of the FeB and Fe{sub 2}B layers. The results showed a fracture resistance of ∼2.8 MPa√m for the FeB layer formed at the surface of borided AISI 1045 steel. Likewise, the effect of the DAP on the surface of the borided AISI 1045 steel promoted the formation of an exclusively Fe{sub 2}B layer, with an increase in the fracture toughness of the whole boride layer around 5 MPa√m. Finally, the principle of the technique can be used to minimize the influence of the anisotropic properties on the fracture toughness along the depth of boride layers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Structural evolution and magnetic hardness of (Sm,Zr)(Fe,Co,Ti) 12 alloy particles via reduction-diffusion

The quest for achieving high coercivity in Sm(Fe,Co,Ti) 12 alloys, despite their inherent strong magnetocrystalline anisotropy, has posed significant challenges. Recently, (Sm,Zr)(Fe,Co,Ti) 12 monocrystalline particles have exhibited coercivity μ 0 H c > 1.2 T, showcasing promising prospects and significant potential for both manufacturing and research endeavors. This study delves into the structural evolution of (Sm,Zr)(Fe,Co,Ti) 12 (1:12) alloy particles made via the calciothermic reduction-diffusion synthesis process as influenced by the molar ratios of Ca atoms to O 2- ions (Ca/O), annealing time and annealing temperature. Critical insight that informs conditions to optimize the magnetic response is gained via systematic experimentation and advanced electron microscopy. Complex structural features, including core-shell morphologies and intricate multiphase compositions within individual particles, are unveiled. An optimal Ca/O ratio of 1.30 produces particles with a coercivity up to μ 0 H c = 1.63 T, while higher Ca/O ratios induce the formation of a Sm-rich TbCu 7 -type (1: 7 ) phase, which only partially transforms into the desired 1:12 phase during annealing. Persistent remnants of the 1:7 phase locally impact atomic structure, particle morphology, and coercivity. Furthermore, these findings underscore the complex interplay between synthesis parameters, resulting structures, and magnetic properties, informing the design and optimization of high-performance permanent magnets comprised of the (1:12) compound.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiolytic degradation of 240 Plutonium and 242 Plutonium oxalates

Raman, FTIR, and diffuse reflectance spectroscopy were used to study the auto-radiolytic degradation of 240 Pu and 242 Pu oxalates. The significant differences in the lifetimes of 240 Pu and 242 Pu enabled the differentiation between environmental and radiolytic mechanisms. 240 Pu oxalates were observed to decompose to PuOCO 3 at intermediate times (~ 20 weeks) followed by partial conversion to PuO 2 at times greater than one year. Atmospheric oxidation was shown to be the primary decomposition mechanism for 242 Pu(IV) oxalate, and the alpha radiolysis of aquo and oxalate ligands serves as a secondary decomposition mechanism. In conclusion, this study offers a fresh perspective on radiolytic aging, which is crucial for long-term storage applications.

Analytical Techniques in Art Conservation↗
DOE OSTI2020

An Experimental Method to Measure Gaseous Diffusivity in Tight and Partially Saturated Porous Media via Continuously Monitored Mass Spectrometry

Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Finally, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Electrodeposition Process and Performance of CuIn(Se{sub x}S{sub 1−x}){sub 2} Film for Absorption Layer of Thin-Film Solar Cells

CuIn(Se{sub x}S{sub 1−x}){sub 2} thin film is prepared by the electrodeposition method for the absorption layer of the solar cell. The CuIn(Se{sub x}S{sub 1−x}){sub 2} films are characterized by cyclic voltammetry measurement for the reduction of copper, indium, selenium and sulfur in selenium and sulfur in aqueous solutions with sodium citrate and without sodium citrate. In the four cases, the defined reduction process for every single element is obtained and it is observed that sodium citrate changes the reduction potentials. A linear relationship between the current density of the reduction peak and (scan rate v){sup 1/2} for copper and indium is achieved, indicating that the process is diffusion controlled. The diffusion coefficients of copper and indium ions are calculated. The diffusional coefficient D value of copper is higher than that of indium, and this is the reason why the deposition rate of copper is higher. When four elements are co-deposited in the aqueous solution with sodium citrate, the quaternary compound of CuIn(Se{sub x}S{sub 1−x}){sub 2} is deposited together with Cu{sub 3}Se{sub 2} impure phases after annealing, as found by XRD spectra. Morphology is observed by SEM and AFM. The chemical state of the films components is analyzed by XPS. The UV-Visible spectrophotometer and electrochemistry workstation are employed to measure the photoelectric properties. The results show that the smooth, uniform and compact CuIn(Se{sub x}S{sub 1−x}){sub 2} film is a semiconductor with a band gap of 1.49 eV and a photovoltaic conversion efficiency of 0.45%.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Influence of Effective Physical Contact Area on Microstructure and Mechanical Properties of Diffusion-Bonded TC4/1060Al Joints

In the present study, dissimilar TC4 alloy and 1060Al alloy were successfully diffusion-bonded at a very low temperature about 410 °C by improving the effective physical contact area. A sound joint with a crack-free interface and a high shear strength (128 MPa) is obtained with a combination of TC4 surface roughness of 109.90 μm (S{sub a}) and 1060Al surface roughness of 101.43 μm (S{sub a}), and customized parallel surface scratches. Under this condition, the two parent metals share the largest contact area during the physical contact stage of diffusion bonding process, which is consistent with the theoretical calculation. The results show that no Ti-Al binary intermetallic compounds are observed at the interface of the joints. The maximum shear strength of the joint reaches 128 MPa with optimized bonding parameters (410 °C for 120 min with a pressure of 20 MPa), which is higher than the shear strength of 1060Al (80 MPa) base metal. The fractography analysis indicates the joints failed with a ductile fracture at the 1060Al side.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Pore-scale study of effects of macroscopic pores and their distributions on reactive transport in hierarchical porous media

For applications of reactive transport in porous media, optimal porous structures should possess both high surface area for reactive sites loading and low mass transport resistance. Hierarchical porous media with a combination of pores at different scales are designed for this purpose. In this paper, using the lattice Boltzmann method, pore-scale numerical studies are conducted to investigate diffusion-reaction processes in 2D hierarchical porous media generated by self-developed reconstruction scheme. Complex interactions between porous structures and reactive transport are revealed under different conditions. Simulation results show that adding macropores can greatly enhance the mass transport, but at the same time reduce the reactive surface, leading to complex change trend of the total reaction rate. Effects of gradient distribution of macropores within the porous medium are also investigated. It is found that a front-loose, back-tight (FLBT) hierarchical structure is desirable for enhancing mass transport, increasing total reaction rate, and improving catalyst utilization. Finally, on the whole, from the viewpoint of reducing cost and improving material performance, hierarchical porous structures, especially gradient structures with the size of macropores gradually decreasing along the transport direction, are desirable for catalyst application.

42 ENGINEERING↗
DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗