Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

MODEL OF THE FIELD LINE RANDOM WALK EVOLUTION AND APPROACH TO ASYMPTOTIC DIFFUSION IN MAGNETIC TURBULENCE

The turbulent random walk of magnetic field lines plays an important role in the transport of plasmas and energetic particles in a wide variety of astrophysical situations, but most theoretical work has concentrated on determination of the asymptotic field line diffusion coefficient. Here we consider the evolution with distance of the field line random walk using a general ordinary differential equation (ODE), which for most cases of interest in astrophysics describes a transition from free streaming to asymptotic diffusion. By challenging theories of asymptotic diffusion to also describe the evolution, one gains insight on how accurately they describe the random walk process. Previous theoretical work has effectively involved closure of the ODE, often by assuming Corrsin's hypothesis and a Gaussian displacement distribution. Approaches that use quasilinear theory and prescribe the mean squared displacement ({Delta}x {sup 2}) according to free streaming (random ballistic decorrelation, RBD) or asymptotic diffusion (diffusive decorrelation, DD) can match computer simulation results, but only over specific parameter ranges, with no obvious 'marker' of the range of validity. Here we make use of a unified description in which the ODE determines ({Delta}x {sup 2}) self-consistently, providing a natural transition between the assumptions of RBD and DD. We find that the minimum kurtosis of the displacement distribution provides a good indicator of whether the self-consistent ODE is applicable, i.e., inaccuracy of the self-consistent ODE is associated with non-Gaussian displacement distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SIMPLE ANALYTICAL FORMS OF THE PERPENDICULAR DIFFUSION COEFFICIENT FOR TWO-COMPONENT TURBULENCE. I. MAGNETOSTATIC TURBULENCE

We explore perpendicular diffusion based on the unified nonlinear transport theory. We derive simple analytical forms for the perpendicular mean free path and investigate the influence of different model spectra. We show that for cases where the field line random walk is normal diffusive, the perpendicular diffusion coefficient consists of only two transport regimes. Details of the spectral shape are less important, especially those of the inertial range. Only the macroscopic properties of the turbulence spectrum control the perpendicular diffusion coefficient. Simple formulae for the perpendicular diffusion coefficient are derived which can easily be implemented in solar modulation or shock acceleration codes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE COLLAPSE OF TURBULENT CORES AND RECONNECTION DIFFUSION

In order for a molecular cloud clump to form stars, some transport of magnetic flux is required from the denser internal regions to the outer regions; otherwise, this can prevent the gravitational collapse. Fast magnetic reconnection, which takes place in the presence of turbulence, can induce a process of reconnection diffusion that has been elaborated on in earlier theoretical works. We have named this process turbulent reconnection diffusion, or simply RD. This paper continues our numerical study of this process and its implications. In particular, we extend our studies of RD in cylindrical clouds and consider more realistic clouds with spherical gravitational potentials (from embedded stars); we also account for the effects of the gas self-gravity. We demonstrate that, within our setup reconnection, diffusion is efficient. We have also identified the conditions under which RD becomes strong enough to make an initially subcritical cloud clump supercritical and induce its collapse. Our results indicate that the formation of a supercritical core is regulated by a complex interplay between gravity, self-gravity, the magnetic field strength, and nearly transonic and trans-Alfvénic turbulence; therefore, it is very sensitive to the initial conditions of the system. In particular, self-gravity helps RD and, as a result, the magnetic field decoupling from the collapsing gas becomes more efficient compared with the case of an external gravitational field. Our simulations confirm that RD can transport magnetic flux from the core of collapsing clumps to the envelope, but only a few of them become nearly critical or supercritical sub-Alfvénic cores, which is consistent with the observations. Furthermore, we have found that the supercritical cores built up in our simulations develop a predominantly helical magnetic field geometry that is also consistent with recent observations. Finally, we have also evaluated the effective values of the turbulent RD coefficient in our simulations and found that they are much larger than the numerical diffusion coefficient, especially for initially trans-Alfvénic clouds, thus ensuring that the detected magnetic flux removal is due to the action of turbulent RD rather than numerical diffusivity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The diffusion approximation versus the telegraph equation for modeling solar energetic particle transport with adiabatic focusing. I. Isotropic pitch-angle scattering

The diffusion approximation to the Fokker-Planck equation is commonly used to model the transport of solar energetic particles in interplanetary space. In this study, we present exact analytical predictions of a higher order telegraph approximation for particle transport and compare them with the corresponding predictions of the diffusion approximation and numerical solutions of the full Fokker-Planck equation. We specifically investigate the role of the adiabatic focusing effect of a spatially varying magnetic field on an evolving particle distribution. Comparison of the analytical and numerical results shows that the telegraph approximation reproduces the particle intensity profiles much more accurately than does the diffusion approximation, especially when the focusing is strong. However, the telegraph approximation appears to offer no significant advantage over the diffusion approximation for calculating the particle anisotropy. The telegraph approximation can be a useful tool for describing both diffusive and wave-like aspects of the cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Interstellar simulations using a unified microscopic-macroscopic Monte Carlo model with a full gas-grain network including bulk diffusion in ice mantles

We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach that includes photon penetration and bulk diffusion. To determine the significance of these two processes, we simulate the chemistry with three different models. In Model 1, we use an exponential treatment to follow how photons penetrate and photodissociate ice species throughout the grain mantle. Moreover, the products of photodissociation are allowed to diffuse via bulk diffusion and react within the ice mantle. Model 2 is similar to Model 1 but with a slower bulk diffusion rate. A reference Model 0, which only allows photodissociation reactions to occur on the top two layers, is also simulated. Photodesorption is assumed to occur from the top two layers in all three models. We found that the abundances of major stable species in grain mantles do not differ much among these three models, and the results of our simulation for the abundances of these species agree well with observations. Likewise, the abundances of gas-phase species in the three models do not vary. However, the abundances of radicals in grain mantles can differ by up to two orders of magnitude depending upon the degree of photon penetration and the bulk diffusion of photodissociation products. We also found that complex molecules can be formed at temperatures as low as 10 K in all three models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ANALYSIS OF MAGNETOROTATIONAL INSTABILITY WITH THE EFFECT OF COSMIC-RAY DIFFUSION

We present the results obtained from the linear stability analysis and 2.5 dimensional magnetohydrodynamic (MHD) simulations of magnetorotational instability (MRI), including the effects of cosmic rays (CRs). We took into account the CR diffusion along the magnetic field but neglected the cross-field-line diffusion. Two models are considered in this paper: the shearing box model and differentially rotating cylinder model. We studied how MRI is affected by the initial CR pressure (i.e., energy) distribution. In the shearing box model, the initial state is uniform distribution. Linear analysis shows that the growth rate of MRI does not depend on the value of the CR diffusion coefficient. In the differentially rotating cylinder model, the initial state is a constant angular momentum polytropic disk threaded by a weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both results are confirmed by MHD simulations. The MHD simulation results show that the outward movement of matter by the growth of MRI is not impeded by the CR pressure gradient, and the centrifugal force that acts on the concentrated matter becomes larger. Consequently, the growth rate of MRI is increased. On the other hand, if the initial CR pressure is uniform, then the growth rate of the MRI barely depends on the value of the CR diffusion coefficient.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Dust transport in MRI turbulent disks: Ideal and non-ideal MHD with ambipolar diffusion

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified MHD simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers (Sc) for radial and vertical turbulent diffusion are Sc{sub r}∼1 and Sc{sub z}≳3, but in the AD dominated regime both Sc{sub r} and Sc{sub z} are ≲1. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Generalized theory of diffusion based on kinetic theory

Here, we propose to use spin hydrodynamics, a two-fluid model of spin propagation, as a generalization of the diffusion equation. We show that in the dense limit spin hydrodynamics reduces to Fick's law and the diffusion equation. In the opposite limit spin hydrodynamics is equivalent to a collisionless Boltzmann treatment of spin propagation. Spin hydrodynamics avoids unphysical effects that arise when the diffusion equation is used to describe to a strongly interacting gas with a dilute corona. We apply spin hydrodynamics to the problem of spin diffusion in a trapped atomic gas. We find that the observed spin relaxation rate in the high-temperature limit [Sommer et al., Nature (London) 472, 201 (2011)] is consistent with the diffusion constant predicted by kinetic theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Odd Diffusivity of Chiral Random Motion

Diffusive transport is characterized by a diffusivity tensor which may, in general, contain both a symmetric and an antisymmetric component. Although the latter is often neglected, we derive Green-Kubo relations showing it to be a general characteristic of random motion breaking time-reversal and parity symmetries, as encountered in chiral active matter. In analogy with the odd viscosity appearing in chiral active fluids, we term this component the odd diffusivity. Furthermore, we show how odd diffusivity emerges in a chiral random walk model, and demonstrate the applicability of the Green-Kubo relations through molecular dynamics simulations of a passive tracer particle diffusing in a chiral active bath.

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

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Phosphorus Doping of Polycrystalline CdTe by Diffusion

Phosphorus diffusion in single crystal and polycrystalline CdTe material is explored using various methods. Dynamic secondary ion mass spectroscopy (SIMS) is used to determine 1D P diffusion profiles. A 2D diffusion model is used to determine the expected cross-sectional distribution of P in CdTe after diffusion anneals. Time of flight SIMS and cross-sectional cathodoluminescence corroborates expected P distributions. Devices fabricated with diffused P exhibit hole concentrations up to low 1015 cm-3, however a subsequent activation anneal enabled hole concentrations greater than 1016 cm-3. CdCl2 treatments and Cu based contacts were also explored in conjunction with the P doping process.

14 SOLAR ENERGY↗
DOE OSTI2016

Hindered cytoplasmic diffusion of inositol trisphosphate restricts its cellular range of action

The range of action of intracellular messengers is determined by their rates of diffusion and degradation. Previous measurements in oocyte cytoplasmic extracts indicated that the Ca 2+ -liberating second messenger inositol trisphosphate (IP 3 ) diffuses with a coefficient (~280 μm 2 s -1 ) similar to that in water, corresponding to a range of action of ~25 μm. Consequently, IP 3 is generally considered a “global” cellular messenger. We also reexamined this issue by measuring local IP 3 -evoked Ca 2+ puffs to monitor IP 3 diffusing from spot photorelease in neuroblastoma cells. Fitting these data by numerical simulations yielded a diffusion coefficient (≤10 μm 2 s -1 ) about 30-fold slower than that previously reported. Here, we propose that diffusion of IP 3 in mammalian cells is hindered by binding to immobile, functionally inactive receptors that were diluted in oocyte extracts. The predicted range of action of IP 3 (<5 μm) is thus smaller than the size of typical mammalian cells, indicating that IP 3 should better be considered as a local rather than a global cellular messenger.

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

Effect of Ion Streaming on Diffusion of Dust Grains in Dissipative System

The presence of strong electric fields in the sheath region of laboratory complex plasma induces an ion drift and perturbs the field around dust grains. The downstream focusing of ions leads to the formation of oscillatory kind of attractive wake potential which superimpose with the normal Debye-Hückel (DH) potential. The structural properties of complex plasma and diffusion coefficient of dust grains in the presence of such a wake potential have been investigated using Langevin dynamics simulation in the subsonic regime of ion flow. The study reveals that the diffusion of dust grains is strongly affected by the ion flow, so that the diffusion changes its character in the wake potential to the DH potential dominant regimes. The dependence of the diffusion coefficient on the parameters, such as the neutral pressure, dust grain size, ion flow velocity, and Coulomb coupling parameter, have been calculated for the subsonic regime by using the Green-Kubo expression, which is based on the integrated velocity autocorrelation function. It is found that the diffusion and the structural property of the system is intimately connected with the interaction potential and significantly get affected in the presence of ion flow in the subsonic regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Numerical Solution of Moving Phase Boundary and Diffusion-Induced Stress of Sn Anode in the Lithium-Ion Battery

Here, we have previously observed a large transient stress in Sn film anodes at the beginning of the Sn-Li 2 Sn 5 phase transformation. To understand this behavior, we use numerical modeling to simulate the kinetics of the 1-D moving boundary and Li diffusion in the Sn anodes. A mixture of diffusion-controlled and interface-controlled kinetics is found. The Li concentration in the Li 2 Sn 5 phase remains near a steady-state profile as the phase boundary propagates, whereas the Li diffusion in Sn is more complicated. Li continuously diffuses into the Sn layer and produces a supersaturation; the Li can then diffuse toward the Sn/Li 2 Sn 5 interface and contribute to further phase transformation. Finally, the evolution of Li concentration in the Sn induces strain which involves rate-dependent plasticity and elastic unloading, resulting in the complex stress evolution that is observed. In the long term, the measured stress is dominated by the stress in the growing Li 2 Sn 5 phase.

25 ENERGY STORAGE↗
DOE OSTI2015

Determination of hydrogen diffusion coefficients in F82H by hydrogen depth profiling with a tritium imaging plate technique

Hydrogen diffusion coefficients in a reduced activation ferritic/martensitic steel (F82H) and an oxide dispersion strengthened F82H (ODS-F82H) have been determined from depth profiles of plasma-loaded hydrogen with a tritium imaging plate technique (TIPT) in the temperature range from 298 K to 523 K. Data on hydrogen diffusion coefficients, D, in F82H, are summarized as D [m{sup 2}*s{sup -1}] =1.1*10{sup -7}exp(-16[kJ mol{sup -1}]/RT). The present data indicate almost no trapping effect on hydrogen diffusion due to an excess entry of energetic hydrogen by the plasma loading, which results in saturation of the trapping sites at the surface and even in the bulk. In the case of ODS-F82H, data of hydrogen diffusion coefficients are summarized as D [m{sup 2}*s{sup -1}] =2.2*10{sup -7}exp(-30[kJ mol{sup -1}]/RT) indicating a remarkable trapping effect on hydrogen diffusion caused by tiny oxide particles (Y{sub 2}O{sub 3}) in the bulk of F82H. Such oxide particles introduced in the bulk may play an effective role not only on enhancement of mechanical strength but also on suppression of hydrogen penetration by plasma loading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI

The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λ L ), transversal diffusivity (λ T ), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.

60 APPLIED LIFE SCIENCES↗