Highly efficient parallel grand canonical simulations of interstitial-driven diffusion-deformation processes
Not provided.
Materials Science↗
SEARCH · Engineering Papers
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.
Entanglement and dynamics of diffusion-annihilation processes with Majorana defects
Not Available
Nahum, Adam↗
Processes of diffusion in a limited ensemble of charged particles in a static magnetic field
The dynamics of a limited ensemble of charged particles in a spatially uniform static magnetic field has been studied analytically and numerically. Calculations have been performed for particles with various masses and charges in a wide range of the parameters of the systems under study. It has been shown that the transverse diffusion coefficient in an ensemble of particles in the magnetic field for weakly nonideal systems is described by the Townsend formula. Analytical estimates have been proposed for the analysis of the effective diffusion radius of the ensemble of particles owing to their thermal motion. The proposed relations have been tested by the numerical simulation of the problem.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)
In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.
03 NATURAL GAS↗
Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics
This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Surface Self-Diffusion Induced Sintering of Nanoparticles
Despite the critical role of sintering phenomena in constraining the long-term durability of nanosized particles, a clear understanding of nanoparticle sintering has remained elusive due to the challenges in atomically tracking the neck initiation and discerning different mechanisms. Through the integration of in situ transmission electron microscopy and atomistic modeling, this study uncovers the atomic dynamics governing the neck initiation of Pt–Fe nanoparticles via a surface self-diffusion process, allowing for coalescence without significant particle movement. Real-time imaging reveals that thermally activated surface morphology changes in individual nanoparticles induce significant surface self-diffusion. Further, the kinetic entrapment of self-diffusing atoms in the gaps between closely spaced nanoparticles leads to the nucleation and growth of atomic layers for neck formation. This surface self-diffusion-driven sintering process is activated at a relatively lower temperature compared to the classic Ostwald ripening and particle migration and coalescence processes. The fundamental insights have practical implications for manipulating the morphology, size distribution, and stability of nanostructures by leveraging surface self-diffusion processes.
36 MATERIALS SCIENCE↗
Escape and finite-size scaling in diffusion-controlled annihilation
In this paper, we study diffusion-controlled single-species annihilation with a finite number of particles. In this reaction-diffusion process, each particle undergoes ordinary diffusion, and when two particles meet, they annihilate. We focus on spatial dimensions d>2 where a finite number of particles typically survive the annihilation process. Using scaling techniques we investigate the average number of surviving particles, M, as a function of the initial number of particles, N. In three dimensions, for instance, we find the scaling law M ~ N 1/3 in the asymptotic regime N»1. We show that two time scales govern the reaction kinetics: the diffusion time scale, T ~ N 2/3 , and the escape time scale, τ ~ N 4/3 . The vast majority of annihilation events occur on the diffusion time scale, while no annihilation events occur beyond the escape time scale.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights
The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Stochastic and fractal analysis of fracture trajectories
Analyses of fracture trajectories are used to investigate structures that fall between 'micro' and 'macro' scales. It was shown that fracture trajectories belong to the class of nonstationary processes. It was also found that correlation distance, which may be related to a characteristic size of a fracture process, increases with crack length. An assemblage of crack trajectory processes may be considered as a diffusive process. Chudnovsky (1981-1985) introduced a 'crack diffusion coefficient' d which reflects the ability of the material to deviate the crack trajectory from the most energetically efficient path and thus links the material toughness to its structure. For the set of fracture trajectories in AISI 304 steel, d was found to be equal to 1.04 microns. The fractal dimension D for the same set of trajectories was found to be 1.133.
Bessendorf, Michael H.↗
Fast Advective Water Flow through Nanochannels in Clay Interlayers: Implications for Moisture Transport in Soils and Unconventional Oil/Gas Production
Water flow in nanometer or sub-nanometer hydrophilic channels bears special importance in diverse fields of science and engineering. However, the nature of such water flow remains elusive. In this work, we report our molecular-modeling results on water flow in a sub-nanometer clay interlayer between two montmorillonite layers. We show that a fast advective flow can be induced by evaporation at one end of the interlayer channel, that is, a large suction pressure created by evaporation (~818 MPa) is able to drive the fast water flow through the channel (~0.88 m/s for a 46 Å-long channel). Scaled up for the pressure gradient to a 2 μm particle, the velocity of water is estimated to be about 95 μm/s, indicating that water can quickly flow through a μm-sized clay particle within seconds. The prediction seems to be confirmed by our thermogravimetric analysis of bentonite hydration and dehydration processes, which indicates that water transport at the early stage of the dehydration is a fast advective process, followed by a slow diffusion process. The possible occurrence of a fast advective water flow in clay interlayers prompts us to reassess water transport in a broad set of natural and engineered systems such as clay swelling/shrinking, moisture transport in soils, water uptake by plants, water imbibition/release in unconventional hydrocarbon reservoirs, and cap rock integrity of supercritical CO 2 storage.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
(Photo)Electrocatalytic CO 2 Reduction at the Defective Anatase TiO 2 (101) Surface
Excessive carbon dioxide (CO 2 ) emissions by combustion of fossil fuels are linked to global warming and rapid climate change. One promising route to lowering the concentration of CO 2 in the atmosphere is to reduce it to useful small molecules via photoelectrocatalytic hydrogenation, which would enable solar energy storage with a zero-carbon emission cycle and perform a more efficient separation of the photogenerated electron and hole pair than pure photocatalysis. Indeed, photoelectrocatalytic CO 2 reduction has been an intense focus of research. Using the density functional theory (DFT), we studied the CO 2 reduction reaction on the defective anatase TiO 2 (101) surface, at both the solvent/catalyst and the electrolyte/catalyst interfaces. The analysis of the electronic structure of the surface shows a contrast between the solvent/catalyst and the electrolyte/catalyst interfaces, which results in the two corresponding catalytic cycles being distinct. Furthermore, our study explains at the electronic and mechanistic levels why methanol is the main product in the presence of the electrolyte and why the overpotential is not only controlled by the reaction process but also by the diffusion process.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Galactic civilizations: Population dynamics and interstellar diffusion
The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found.
Newman, W. I.↗
Numerical experiments on the role of radiative processes in the development and maintenance of upper level clouds
A time-dependent, two-dimensional Eulerian model is presented whose purpose is to obtain more realistic parameterizations of extended high level cloudiness, and the results of a numerical experiment using the model are reported. The model is anelastic and the Bousinesque assumption is invoked. Unresolved subgrid scale processes are parameterized as eddy diffusion processes. Two phases of water are incorporated and equilibrium between them is assumed. The effects of infrared radiative processes are parametrically represented. Two simulations were conducted with identical initial conditions; in one of them, the radiation term was never turned on. The mean values of perturbation potential temperature at each level in the domain are plotted versus height after 15, 30, and 60 minutes of simulated time. The influence of the radiative term is seen to impose a cooling trend, leading to an increased generation of ice water and an increased generation of turbulent kinetic energy in the cloud layer.
Starr, D. O'C.↗
Diffusion mechanisms in chemical vapor-deposited iridium coated on chemical vapor-deposited rhenium
Radiation-cooled rocket thruster chambers have been developed which use CVD Re coated with CVD Ir on the interior surface that is exposed to hot combustion gases. The Ir serves as an oxidation barrier which protects the structural integrity-maintaining Re at elevated temperatures. The diffusion kinetics of CVD materials at elevated temperatures is presently studied with a view to the prediction and extension of these thrusters' performance limits. Line scans for Ir and Re were fit on the basis of a diffusion model, in order to extract relevant diffusion constants; the fastest diffusion process is grain-boundary diffusion, where Re diffuses down grain boundaries in the Ir overlayer.
Hamilton, J. C.↗
Low cost solar array project. Cell and module formation research area. Process research of non-CZ silicon material
Liquid diffusion masks and liquid dopants to replace the more expensive CVD SiO2 mask and gaseous diffusion processes were investigated. Silicon pellets were prepared in the silicon shot tower; and solar cells were fabricated using web grown where the pellets were used as a replenishment material. Verification runs were made using the boron dopant and liquid diffusion mask materials. The average of cells produced in these runs was 13%. The relationship of sheet resistivity, temperature, gas flows, and gas composition for the diffusion of the P-8 liquid phosphorus solution was investigated. Solar cells processed from web grown from Si shot material were evaluated, and results qualified the use of the material produced in the shot tower for web furnace feed stock.
Source record↗
Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas
Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.
Plasma kinetic theory↗
Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review
The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.
25 ENERGY STORAGE↗
Mass dependence of overshooting beneath convective envelopes
State of the art stellar evolution codes [Paxton et al., 2010, Demarque et al., 2004, Weiss and Schlattl, 2008, Siess et al., 2013, Christensen-Dalsgaard, 2008] evoke a diffusive process to model convective overshooting. The diffusion coefficient [Freytag et al., 1996, Pratt et al., 2017] used in the model can be set to change with the classification of the convective zone as non-burning, H-burning, He-burning, or metal-burning. It can also be set to change at defined evolutionary points, such as the bottom of the asymptotic giant branch, or during the third dredge up [Herwig, 2000, Lugaro et al., 2003]. Aside from these abrupt changes, the diffusion coefficient is typically locked to a percentage of the pressure scale height measured at the convective boundary. However, there is no theoretical reason for convective overshooting, or indeed other convective properties, to change in the same way that the pressure scale height changes as a star evolves.
79 ASTRONOMY AND ASTROPHYSICS↗