Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗
DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

In situ investigation of high-pressure hydrogen-induced swelling in elastomers and its correlation with material properties

The resistance of elastomeric materials to high-pressure hydrogen-induced damage is essential for ensuring the reliability of hydrogen infrastructure. Here, in this study, we systematically investigated the swelling behavior and hydrogen transport properties of four elastomer types – EPDM, NBR, FKM, and HNBR – using a custom in-situ view cell system capable of real-time monitoring during decompression from pressures up to 96.5 MPa. Each elastomer was formulated with and without fillers and plasticizers to assess the effects of formulation on swelling response. Thermal desorption analysis (TDA) was employed to determine equilibrium hydrogen content and diffusion coefficients, providing insight into gas uptake and mobility within each material. Correlation analyses using Pearson and Spearman coefficients revealed that the diffusion coefficient showed a stronger relationship with swelling behavior than hydrogen content, highlighting the dominant role of hydrogen mobility. Filled elastomers, particularly those with carbon black, consistently showed reduced swelling due to enhanced stiffness and reduced diffusivity. These results deepen our understanding of diffuso-mechanical interactions in elastomers and support the rational design of sealing materials for high-pressure hydrogen systems.

EPDM↗
DOE OSTI2022

Nanoscale Hybrid Electrolytes with Viscosity Controlled Using Ionic Stimulus for Electrochemical Energy Conversion and Storage

As renewable energy is rapidly integrated into the grid, the challenge has become storing intermittent renewable electricity. Technologies including flow batteries and CO 2 conversion to dense energy carriers are promising storage options for renewable electricity. To achieve this technological advancement, the development of next generation electrolyte materials that can increase the energy density of flow batteries and combine CO 2 capture and conversion is desired. Liquidlike nanoparticle organic hybrid materials (NOHMs) composed of an inorganic core with a tethered polymeric canopy (e.g., polyetheramine (HPE)) have a capability to bind chemical species of interest including CO 2 and redox-active species. In this study, the unique response of NOHM-I-HPE-based electrolytes to salt addition was investigated, including the effects on solution viscosity and structural configurations of the polymeric canopy, impacting transport behaviors. The addition of 0.1 M NaCl drastically lowered the viscosity of NOHM-based electrolytes by up to 90%, reduced the hydrodynamic diameter of NOHM-I-HPE, and increased its self-diffusion coefficient, while the ionic strength did not alter the behaviors of untethered HPE. This study is the first to fundamentally discern the changes in polymer configurations of NOHMs induced by salt addition and provides a comprehensive understanding of the effect of ionic stimulus on their bulk transport properties and local dynamics. These insights could be ultimately employed to tailor transport properties for a range of electrochemical applications.

25 ENERGY STORAGE↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Impact of Alloy-Disorder-Induced Localization on Hole Diffusion in Highly Excited c -Plane and m -Plane (In, Ga) N Quantum Wells

The diffusion coefficient of holes can provide knowledge about carrier localization in (In, Ga)N, where the carrier dynamics are altered by randomly fluctuating potential landscape. In group-III nitrides, the diffusivity of holes is difficult to measure by electrical methods but it can be studied using optical techniques. Here, we investigate the dependence of the hole diffusion coefficient on direction and carrier density in c-plane and m-plane (In, Ga)N structures by employing the light-induced transient-grating technique. We show that the hole diffusion coefficient is anisotropic in the m-plane structure, where it is several times larger along the a crystallographic direction than along the c direction. Such anisotropy is observed within the broad range of carrier densities from 10 18 to 10 20 cm –3 . The diffusivity changes nonmonotonously with increasing photoexcitation, this dependence being different in thick and thin layers. We argue that an unexpectedly high diffusion coefficient at low carrier densities in thick quantum wells can be a signature of efficient hole transport via percolative paths occurring due to compositional disorder. In turn, a decrease of diffusivity with the excitation can reflect the effect of Coulomb blockade of these paths. Lastly, we demonstrate that disorder impacts carrier diffusivity even at carrier densities above 10 19 cm –3 , where the overflow of localized states must be included to explain the observed increase of the diffusion coefficient with the carrier density.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock

Collisional plasma shocks generated from supersonic flows are an important feature in many astrophysical and laboratory high-energy-density plasmas. Compared to single-ion-species plasma shocks, plasma shock fronts with multiple ion species contain additional structure, including interspecies ion separation driven by gradients in species concentration, temperature, pressure, and electric potential. We present time-resolved density and temperature measurements of two ion species in collisional plasma shocks produced by head-on merging of supersonic plasma jets, allowing determination of the ion diffusion coefficients. Here our results provide the first experimental validation of the fundamental inter-ion-species transport theory. The temperature separation, a higher-order effect reported here, is valuable for advancements in modeling HED and ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects to be discovered by future direct imaging surveys. We present a low-resolution, R ∼ 300, M-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K {sub zz}) to estimate the strength of vertical mixing in each object. The K {sub zz} values of cooler gaseous objects are close to their theoretical maximum, and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K {sub zz} values imply that disequilibrium phosphine (PH{sub 3}) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH{sub 3} absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's M-band spectrum by including both CO and water clouds in the atmospheric model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Theoretical Analysis of Experimental Data of Sodium Diffusion in Oxidized Molybdenum Thin Films

In this work, the diffusion process of sodium (Na) in molybdenum (Mo) thin films while it was deposited on soda lime glass (SLG) was studied. A small amount of oxygen was present in the chamber while the direct-current (DC) magnetron sputtering was used for the deposition. The substrate temperatures were varied to observe its effect. Such molybdenum films, with or without oxidations, are often used in thin film solar cells, either as back contact or as hole transport layers. Secondary ion mass spectrometry (SIMS) was used to quantify the concentration of the species. A grain diffusion mechanistic model incorporating the effect of grain and grain boundary geometrical shape and size was developed. The model was used to provide an in-depth theoretical analysis of the sodium diffusion in molybdenum thin films that lead to the measured SIMS data. It was observed that not only diffusion coefficients should be considered when analyzing diffusion processes in thin films but also the ratio of grain boundary size to grain size. Both depend on substrate temperature and directly affect the amount of diffused species in the film. The data were analyzed under the light of the film growth speed versus diffusion front speed, the effect of oxygen content, and the effect of substrate temperature on the overall diffusion process. The temperature inversely affects the ratio of grain boundary size and grain size and directly affects the diffusion coefficient, which leads to a preferable temperature at which the highest amount of alkali can be found in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Halide Substitution Effects on Lithium-Ion Diffusion in Protonated Antiperovskites

Solid-state electrolytes (SSEs) for all-solid-state lithium-ion batteries are generating intense interest because these batteries can improve the safety and performance compared with devices fabricated with conventional, flammable liquid electrolytes. In these SSEs, it has been suggested that Li + ion diffusion in the grain boundaries is hindered and is a critical determinant of the overall ionic conductivity (σ). However, Li + ion diffusivities in the grain (D G ) and the grain boundary (D GB ) are difficult to determine experimentally, with few techniques capable of distinguishing the individual contributions. Here, we distinguished the D G and D GB for the protonated lithium antiperovskites (pLiAPs) SSEs: Li 2 OHCl, Li 2 OHBr, Li 2 OHF 0.1 Cl 0.9 , Li 2 OHF 0.1 Br 0.9 , and Li 2 OHCl 0.3 7Br 0.63 . The measurements were obtained directly from 7 Li pulsed-field gradient nuclear magnetic resonance (PFG-NMR) at 353 K. The 7 Li PFG-NMR echo profiles were composed of two primary components with additional secondary oscillatory components – the so-called NMR diffraction phenomenon. The length scale separating the two main components corresponds to a diffusion length of ~1.7 µm, which is thought to be the average grain size (by diameter). The short-range (≤ 1.7 µm) diffusion component associated with D G (≈10 -11 m 2 /s) varied minimally with halide substitution, while the long-range (≥ 1.7 µm) component D GB (≈10 -12 to 10 -15 m 2 /s) was highly sensitive to the substitution of halides and closely correlated with s. In addition, from the comparison of the ratio D GB /D G to D t (the Li + ion diffusion coefficient estimated from the rotational correlation time, t c ), it was determined that the contribution of D G to σ is negligible; 0.01 ~ 0.04 in the pLiAPs studied here. Finally, these insights provide fundamental understanding of the halide substitution effects on Li + ion grain versus grain boundary diffusion, and suggest that careful engineering of the grain boundaries at the microscopic level is necessary to achieve high-performance pLiAP SSEs.

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

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2024

Simulating Long–Term Dynamics of Radiation Belt Electrons Using DREAM3D Model

Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Extremely Slow Diffusion of Gold Nanoparticles under Confinement in Mesoporous Silica

Here, we report an X-ray photon correlation spectroscopy (XPCS) study of the mobility of colloidal gold nanoparticles with a diameter of approximately 4 nm in glycerol under confinement within the channels of SBA-15 mesoporous silica with pore diameters of approximately 6 and 12 nm. The XPCS correlation functions result from an effective heterodyne signal due to the mixing of the coherent scattering from the nanoparticles and the mesoporous silica. Over the range of wave vectors and hence length scales probed, the nanoparticle dynamics are well described by one-dimensional diffusion. The nanoparticle diffusivity varies with temperature in a manner expected based on the temperature-dependent viscosity of glycerol; however, the magnitudes of the diffusion coefficients are several orders of magnitude smaller than those of the nanoparticles in bulk glycerol. We consider mechanisms that might contribute to this reduction in diffusivity including enhancement in hydrodynamic drag under confinement, effects of nanoparticle adsorption to the pore walls, and slowing structural dynamics in glycerol due to the combination of the nanometer-scale confinement and the presence of the nanoparticles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Photocarrier dynamics at the interface between WS 2 and ozone-irradiated graphene

We investigate the photocarrier properties of heterostructures (HSs) formed by WS 2 and ozone-irradiated graphene (OI-graphene), with a focus on the impact of ozone irradiation. Photoluminescence and Raman spectroscopy reveal a pronounced doping effect induced by ozone irradiation on the graphene layer. Transient absorption measurements demonstrate a substantial suppression of charge transfer efficiency across the interface in the presence of ozone molecules. Spatially resolved pump-probe measurements further show a significant reduction in the diffusion coefficient of photocarriers in WS 2 on OI-graphene compared to that of individual WS 2 monolayers. The slower diffusion suggests enhanced spatial separation of electrons and holes in the WS 2 layer, induced by a built-in electric field at the interface. This effect arises from the combined influences of doping and dielectric property changes caused by ozone molecules. These findings highlight the critical influence of ozone irradiation on graphene and HS interfaces, offering valuable insights for the design of future graphene-based two-dimensional HSs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗