Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion length”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Temperature Dependence of Volumetric and Dynamic Properties of Imidazolium-Based Ionic Liquids

Atomistically detailed molecular dynamics simulations were used to investigate the temperature dependence of the specific volume, dynamic properties, and viscosity of linear alkyl chain ([C n C 1 Im][NTf 2 ], n = 3–7) and branched alkyl chain ([(n – 2)mC n–1 C 1 Im][NTf 2 ]) ionic liquids (ILs). The trend of the glass transition temperature (T g ) values obtained in the simulations as a function of the alkyl chain length of cations was similar to the trend seen in experiments. In addition, the system relaxation behavior as determined from the temperature dependence of the diffusion coefficient, rotational relaxation time, and viscosity close to T g was observed to follow the Vogel–Fulcher–Tammann expression. Furthermore, the reciprocal of the diffusion coefficient of the anion and cation in both linear and branched IL systems showed a linear correlation with viscosity, thus confirming the validity of the Stokes–Einstein relationship for these systems. Similarly, the average rotational relaxation time of the ions was also found to correlate linearly with the viscosity of the ILs over a wide range of temperatures, thereby validating the Debye–Stokes–Einstein relationship for the ILs. These simulation findings suggest that the temperature dependence of the relaxation time of ILs is very similar to that of other glass-forming liquids.

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

Anderson Localization of Thermal Phonons Leads to a Thermal Conductivity Maximum

Our elastic model of ErAs disordered GaAs/AlAs superlattices exhibits a local thermal conductivity maximum as a function of length due to exponentially suppressed Anderson-localized phonons. By analyzing the sample-to-sample fluctuations in the dimensionless conductance, g, the transition from diffusive to localized transport is identified as the crossover from the multichannel to single-channel transport regime g ≈ 1. Single parameter scaling is shown to hold in this crossover regime through the universality of the probability distribution of g that is independent of system size and disorder strength.

42 ENGINEERING↗
DOE OSTI2023

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Communication: Polymer entanglement dynamics: Role of attractive interactions

The coupled dynamics of entangled polymers, which span broad time and length scales, govern their unique viscoelastic properties. To follow chain mobility by numerical simulations from the intermediate Rouse and reptation regimes to the late time diffusive regime, highly coarse grained models with purely repulsive interactions between monomers are widely used since they are computationally the most efficient. In this paper, using large scale molecular dynamics simulations, the effect of including the attractive interaction between monomers on the dynamics of entangled polymer melts is explored for the first time over a wide temperature range. Attractive interactions have little effect on the local packing for all temperatures T and on the chain mobility for T higher than about twice the glass transition T g . Finally, these results, across a broad range of molecular weight, show that to study the dynamics of entangled polymer melts, the interactions can be treated as pure repulsive, confirming a posteriori the validity of previous studies and opening the way to new large scale numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Effect of eddy diffusivity on wind-driven currents in a two-layer stratified lake

The steady state wind-driven circulation was numerically calculated in a rectangular stratified lake. The lake is composed of two layers having uniform but unequal densities and eddy diffusivities. The position in thermocline and the three-dimensional velocities in both layers calculated using shallow lake equations. The results show that, as the eddy diffusivity in the hypolimnion is increased, the thermocline tilt and hypolimnetic velocities increase. The effect of the other variables such as wind stress, density, basin length, and mean thermocline depth are also shown.

Gedney, R. T.↗
NASA NTRS1985

On the structure of the turbulent vortex

The trailing vortex generated by a lifting surface, the structure of its turbulent core and the influence of axial flow within the vortex on its initial persistence and on its subsequent decay are described. Similarity solutions of the turbulent diffusion equation are given in closed form and results are expressed in sufficiently simple terms that the influence of the lifting surface parameters on the length of persistence and the rate of decay of the vortex can be evaluated.

Roberts, L.↗
NASA NTRS1999

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗
NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗
NASA NTRSDate not supplied

Dependence of Lightning Cloud-Top Optical Emissions on Channel Characteristics Using Boltzmann Radiative Transfer Modeling

With an estimated global flash rate of 46 lightning per second, a significant amount of lightning nitrogen oxides (LNOx) are produced within the climate system. In turn, LNOx affect the concentrations of ozone (O 3 , an important greenhouse gas) as well as the very reactive hydroxyl radicals (OH, that affect all greenhouse gases). Thus, space-based estimation of LNOx is highly sought [Koshak 2017; Koshak et al. (2015)]. Total lightning channel length, L , is a proxy for the total flash energy (and hence LNOx production). Geostationary lightning mappers provide cloud-top 2D views of the diffuse (i.e., blurry) optical images of lightning, thereby making it difficult to determine L (and other channel features). Transient Boltzmann radiative transfer modeling is applied here to investigate the dependence of the spatio-temporal pattern of cloud-top lightning optical emissions on L and other channel characteristics (e.g., orientation, propagation direction, amplitude, location, ...).

William J Koshak↗
DOE OSTI2015

Correlated two-particle diffusion in dense colloidal suspensions at early times: Theory and comparison to experiment

The spatially resolved diffusive dynamic cross correlations of a pair of colloids in dense quasi-two-dimensional monolayers of identical particles are studied experimentally and theoretically at early times where motion is Fickian. In very dense systems where strong oscillatory equilibrium packing correlations are present, we find an exponential decay of the dynamic cross correlations on small and intermediate length scales. At large separations where structure becomes random, an apparent power law decay with an exponent of approximately –2.2 is observed. For a moderately dense suspension where local structural correlations are essentially absent, this same apparent power law decay is observed over all probed interparticle separations. A microscopic nonhydrodynamic theory is constructed for the dynamic cross correlations which is based on interparticle frictional effects and effective structural forces. Hydrodynamics enters only via setting the very short-time single-particle self-diffusion constant. No-adjustable-parameter quantitative predictions of the theory for the dynamic cross correlations are in good agreement with experiment over all length scales. The origin of the long-range apparent power law is the influence of the constraint of fixed interparticle separation on the amplitude of the mean square force exerted on the two tagged particles by the surrounding fluid. The theory is extended to study high-packing-fraction 3D hard sphere fluids. Here, the same pattern of an oscillatory exponential form of the dynamic cross correlation function is predicted in the structural regime, but the long-range tail decays faster than in monolayers with an exponent of –3.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1999

Noise and Bandwidth Measurements of Diffusion-Cooled Nb Hot-Electron Bolometer Mixers at Frequencies Above the Superconductive Energy Gap

Diffusion-cooled Nb hot-electron bolometer (HEB) mixers have the potential to simultaneously achieve high intermediate frequency (IF) bandwidths and low mixer noise temperatures for operation at THz frequencies (above the superconductive gap energy). We have measured the IF signal bandwidth at 630 GHz of Nb devices with lengths L = 0.3, 0.2, and 0.1 micrometer in a quasioptical mixer configuration employing twin-slot antennas. The 3-dB EF bandwidth increased from 1.2 GHz for the 0.3 gm long device to 9.2 GHz for the 0.1 gm long device. These results demonstrate the expected 1/L squared dependence of the IF bandwidth at submillimeter wave frequencies for the first time, as well as the largest EF bandwidth obtained to date. For the 0.1 gm device, which had the largest bandwidth, the double sideband (DSB) noise temperature of the receiver was 320-470 K at 630 GHz with an absorbed LO power of 35 nW, estimated using the isothermal method. A version of this mixer with the antenna length scaled for operation at 2.5 THz has also been tested. A DSB receiver noise temperature of 1800 plus or minus 100 K was achieved, which is about 1,000 K lower than our previously reported results. These results demonstrate that large EF bandwidth and low-noise operation of a diffusion-cooled HEB mixer is possible at THz frequencies with the same device geometry.

Wyss, R. A.↗
DOE OSTI2026

Controlling homogenization length scales and microstructure in additively manufactured Ti-Ta functionally graded materials

Materials with smooth compositional gradients or functionally grade materials (FGMs) produced via additive manufacturing (AM), enables joining dissimilar materials and optimizing multiple properties in advanced engineering applications. However, as-printed AM microstructures exhibit micro-segregation and solidification defects which, when combined with controlling macroscale gradient properties, complicates necessary post-processing. Here, we use CALPHAD-informed diffusion modelling to design post-processing heat treatments for lightweight to refractory FGMs. Ti-Ta (0 to 85 at. % Ta) FGMs were fabricated using laser-based directed energy deposition AM. Post-processing heat treatments at 1000° C and 1500° C were designed to promote homogenization across specific length scales and experimentally validated. Investigation of chemical segregation and microstructures demonstrated that the length scale of homogenization is controlled as a function of time, temperature, and local composition. Ta-rich regions exhibited incomplete homogenization compared to Ti-rich layers. Unmelted Ta particles were found to completely dissolve at 1500 °C. By controlling cooling rate (200 °C/min), martensitic structures were produced between 14–36 at. % Ta, consistent with martensite-start temperatures calculations, while furnace cooling (2 °C/min) produced α+β morphologies. This work establishes a validated predictive framework for designing post-processing to tailor microstructure and chemical architecture in AM FGMs, facilitating their deployment in demanding environments.

Materials science↗
DOE OSTI2020

Electrostatic Barriers to Nanoparticle Accessibility of a Porous Matrix

Translocation from one cavity to another through a narrow constriction (i.e., a “hole”) represents the fundamental elementary process underlying hindered mass transport of nanoparticles and macromolecules within many natural and synthetic porous materials, including intracellular environments. This process is complex and may be influenced by long-range (e.g., electrostatic) particle–wall interactions, transient adsorption/desorption, surface diffusion, and hydrodynamic effects. In this study, we used a three-dimensional (3D) tracking method to explicitly visualize the process of nanoparticle diffusion within periodic porous nanostructures, where electrostatic interactions were mediated via ionic strength. The effects of electrostatic interactions on nanoparticle transport were surprisingly large. For example, an increase in the Debye length of only a few nanometers (in a material with a hole diameter of ~100 nm) increased the mean cavity escape time 3-fold. A combination of computational and experimental analyses indicated that this hindered cavity escape was due to an electrostatic energy barrier in the region of the hole, which was quantitatively explained using DLVO theory. These findings explicitly demonstrate that the cavity escape process was barrier-limited and dominated by electrostatic effects.

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

Detailed Characterization of a Nanosecond-Lived Excited State: X-ray and Theoretical Investigation of the Quintet State in Photoexcited [Fe(terpy) 2 ] 2+

Theoretical predictions show that depending on the populations of the Fe 3dxy, 3dxz, and 3dyz orbitals two possible quintet states can exist for the high-spin state of the photoswitchable model system [Fe(terpy)2]2+. The differences in the structure and molecular properties of these 5B2 and 5E quintets are very small and pose a substantial challenge for experiments to resolve them. Yet for a better understanding of the physics of this system, which can lead to the design of novel molecules with enhanced photoswitching performance, it is vital to determine which high-spin state is reached in the transitions that follow the light excitation. The quintet state can be prepared with a short laser pulse and can be studied with cutting-edge time-resolved X-ray techniques. Here we report on the application of an extended set of X-ray spectroscopy and scattering techniques applied to investigate the quintet state of [Fe(terpy)2]2+ 80 ps after light excitation. High-quality X-ray absorption, nonresonant emission, and resonant emission spectra as well as X-ray diffuse scattering data clearly reflect the formation of the high-spin state of the [Fe(terpy)2]2+ molecule; moreover, extended X-ray absorption fine structure spectroscopy resolves the Fe-ligand bond-length variations with unprecedented bondlength accuracy in time-resolved experiments. With ab initio calculations we determine why, in contrast to most related systems, one configurational mode is insufficient for the description of the low-spin (LS)-high-spin (HS) transition. We identify the electronic structure origin of the differences between the two possible quintet modes, and finally, we unambiguously identify the formed quintet state as 5E, in agreement with our theoretical expectations.

Vankó, György↗
DOE OSTI2024

Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel

Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multicomponent mutual diffusion in the warm, dense matter regime

In this work, we present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗