Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Diffusion of a multi-species component and its role in oxygen and water transport in silicates

The diffusion of a multispecies component is complicated by the different diffusion coefficient of each species and the interconversion reactions among the species. A diffusion equation is derived that incorporates the diffusive fluxes of all species contributing to the component's concentration. The effect of speciation on diffusion is investigated experimentally by measuring concentration profiles of all species developed during diffusion experiments. Data on water diffusion in rhyolitic glasses indicate that H2O molecules predominate over OH groups as the diffusing species at very low to high water concentrations. A simple theoretical relationship is drawn between the effective total oxygen diffusion coefficient and the total water concentration of silicates at low water content.

Zhang, Youxue↗
NASA NTRS1994

A unifying comparison of nearly scatter free transport models

Gombosi et al. (1993) recently derived a modified telegrapher's equation for charged particle transport under the influence of isotropic scattering. This equation obeys causality and disallows upstream diffusion for particles with random velocities smaller than the bulk flow velocity. The acausal diffusion equation was obtained to lowest order in the expansion of smallness prameters. The paper by Gombosi et al. (1993) prompted responses from Pauls et al. (1993) and Earl (1993). This paper is written to explain the differences between the methods, assumptions, and results of Gombosi et al. (1993), Pauls et al. (1993), and Earl (1993) and presents a new method of obtaining approximate solutions. It is shown that the assumptions used by Gombosi et al. (1993) and Pauls et al. (1993) are physically equivalent. In our solution method, the solution of the modified telegrapher's equation is obtained as the casual limit of solutions accurate to second order in the smallness parameter expansion. In order to investigate the coherent velocity, we have also developed `wavenumber eigenfunctions' which account for all the pitch angle dependence in our Boltzmann equation. Using truncation, Earl (1993) obtains approximations for the wavenumber dependence of the lowest two frequency modes, which correspond to two of the wavenumber eigenmodes. We find that a consequence of including only two wavenumber eigenmodes is that one obtains solutions which disobey causality at sufficiently short times. Furthermore, the coherent velocity of the two eigenmodes is strongly dependent on wavenumber and approaches the particle velocity in the limit of large wavenumber for both isotropic and anisotropic scattering processes. We conclude that Earl's (1993) solutions and solutions obtained using the new solution method implicitly assume weak acausality and reasonable behavior in the temporal regime, t less than 4 tau. The solutions are not strictly consistent with the behavior of the lowest two frequency modes but have similar behavior in the regime of low wavenumber.

Schwadron, N. A.↗
DOE OSTI2024

Theory of Beam Echoes

We develop the theory of beam echoes in circular accelerators under several different conditions. We derive detailed expressions for the echo amplitude and pulse width with nonlinear quadrupole and dipole kicks, first without and then with momentum spread. We use the theory with the linearized dipole and quadrupole kicks to solve the diffusion equation for different dependencies of the diffusion coefficient on the action. We then consider the use of multiple quadrupole kicks to increase the maximum echo amplitude. We have extended these calculations partially to the 2D case and we also have partial results for longitudinal echoes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Solvent and solute ingress into hydrogels resolved by a combination of imaging techniques

Using simultaneous neutron, fluorescence, and optical brightfield transmission imaging, the diffusion of solvent, fluorescent dyes, and macromolecules into a crosslinked polyacrylamide hydrogel was investigated. This novel combination of different imaging techniques enables us to distinguish the movements of the solvent and fluorescent molecules. Additionally, the swelling or deswelling of the hydrogels can be monitored. From the sequence of images, dye and solvent concentrations were extracted spatially and temporally resolved. Diffusion equations and different boundary conditions, represented by different models, were used to quantitatively analyze the temporal evolution of these concentration profiles and to determine the diffusion coefficients of solvent and solutes. Solute size and network properties were varied and their effect was investigated. Increasing the crosslinking ratio or partially drying the hydrogel was found to hinder solute diffusion due to the reduced pore size. By contrast, solvent diffusion seemed to be slightly faster if the hydrogel was only partially swollen and hence solvent uptake enhanced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Modelling of Trapped Radiation Near Jupiter

Energetic (62 to approx. 130 MeV) proton fluxes measured with the Galileo Probe inside Jupiter's main ring (radii 122,500 to 128,940 km) show a modest increase as Jupiter is approached. Solutions of a reduced (1-d) equatorial diffusion equation with-constant losses described by a lifetime tau match those data for tau greater than or equal to 10(exp 9) s for diffusion coefficient D(sub LL) = 10(exp -9) L(exp 4)/s, 10(exp -9) L(exp -3)/s, or 10(exp -10) L(exp 4)/s (this particle population may be undergoing nearly loss-free inward radial diffusion). Tau would increase as Jupiter is approached if its value were determined by pitch angle scattering due to EM wave-particle interactions. if the absolute amplitudes of the waves' magnetic fluctuations did not vary with radial distance. Exploration of numerical solutions of the diffusive transport equation for Jupiter's magnetospheric ions at die planetary ring and closer to Jupiter has been done. Explicit range-energy relationships for energy loss in ring matter (SiO2), modeled as a continuous disk, are incorporated. A spatial resolution of approx. 1000 km in the radial direction is used in the vicinity of the main ring; off-equatorial ion fluxes have been included (assuming the atmosphere is a perfect absorber). For protons the maximum effect of energy loss in the microscopic size range of ring particles is expected at energies approx. 0.1 Mev. Presumably there could be sufficiently large, unobserved bodies within the ring that would have to be modeled as discrete objects, as planetary satellites would be. Within limitations of the model and computational resources it's planned to characterize solutions for transport of magnetospheric ions past the planetary ring; the procedures that are used are fairly well known and may be applied generally to the simpler magnetospheric models.

Mihalov, John D.↗
DOE OSTI2018

Two-Level Least Squares Methods in Krylov Subspaces

Two-level least squares acceleration approaches are applied to the Chebyshev acceleration method and the restarted conjugate residual method in solving systems of linear algebraic equations with sparse unsymmetric coefficient matrices arising from finite volume or finite element approximations of boundary-value problems on irregular grids. Application of the proposed idea to other iterative restarted processes also is considered. The efficiency of the algorithms suggested is investigated numerically on a set of model Dirichlet problems for the convection-diffusion equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Self-similar and diffusive expansion of nonextensive plasmas

Exact analytical self-similar solution is presented for free collisionless expansion of a two-component plasma of inertial ions and nonextensive electrons into vacuum, using the generalized nonextensive velocity distribution for electrons. Furthermore, a hydrodynamic model of plasma expansion in the presence of the ambipolar diffusion caused by collisions among the plasma species, such as electrons and ions, is developed and a Fokker-Planck-like generalized diffusion equation for steady-state expansion of a nonextensive electron-ion plasma is derived. For the case of generalized statistics and in the absence of particle diffusion, the density, velocity, electric potential, and field of expansion profiles are exactly obtained and studied in terms of the self-similar parameter. It is found that superthermal electrons lead to an accelerated expansion of plasma compared to that of Maxwellian electrons. It is also revealed that the nonextensivity parameter plays a fundamental role on the density, velocity, electric potential, and field configuration of the expansion. Therefore, one is able to distinguish three different regimes q < 1, q = 1, and q > 1 for expansion corresponding to sub-nonextensive, extensive, and super-nonextensive statistical profiles for electrons, respectively. Current research can provide useful information and suggests techniques for investigation of the involved statistical mechanism on the role of the energetic electron fluid in the expansion of plasma in strong pulsed laser-matter interaction experiments. It is also shown that the particle diffusion expansion mechanism becomes more dominant for relatively large values of the nonextensivity parameter, q.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

SMALL-SCALE ANISOTROPIES OF COSMIC RAYS FROM RELATIVE DIFFUSION

The arrival directions of multi-TeV cosmic rays show significant anisotropies at small angular scales. It has been argued that this small-scale structure can naturally arise from cosmic ray scattering in local turbulent magnetic fields that distort a global dipole anisotropy set by diffusion. We study this effect in terms of the power spectrum of cosmic ray arrival directions and show that the strength of small-scale anisotropies is related to properties of relative diffusion. We provide a formalism for how these power spectra can be inferred from simulations and motivate a simple analytic extension of the ensemble-averaged diffusion equation that can account for the effect.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

rHARM: ACCRETION AND EJECTION IN RESISTIVE GR-MHD

Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity and respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests, and preliminary applications. The test simulations of the new code rHARM are compared to an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magnetorotational instability (MRI) in tori around black holes (BHs) for a range of magnetic diffusivities. We find an indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the BH is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild BHs that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Sagittarius A* as an origin of the Galactic PeV cosmic rays?

Supernova remnants (SNRs) have commonly been considered as a source of the observed PeV cosmic rays (CRs) or a Galactic PeV particle accelerator ('Pevatron'). In this work, we study Sagittarius A* (Sgr A*), which is the low-luminosity active galactic nucleus of the Milky Way Galaxy, as another possible canditate of the Pevatron, because it sometimes became very active in the past. We assume that a large number of PeV CRs were injected by Sgr A* at the outburst about 10{sup 7} yr ago when the Fermi bubbles were created. We constrain the diffusion coefficient for the CRs in the Galactic halo on the condition that the CRs have arrived on the Earth by now, while a fairly large fraction of them have escaped from the halo. Based on a diffusion-halo model, we solve a diffusion equation for the CRs and compare the results with the CR spectrum on the Earth. The observed small anisotropy of the arrival directions of CRs may be explained if the diffusion coefficient in the Galactic disk is smaller than that in the halo. Our model predicts that a boron-to-carbon ratio should be energy-independent around the knee, where the CRs from Sgr A* become dominant. It is unlikely that the spectrum of the CRs accelerated at the outburst is represented by a power-law similar to the one for those responsible for the gamma-ray emission from the central molecular zone (CMZ) around the Galactic center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Numerical simulations of chaotic and complex spatiotemporal patterns in fractional reaction–diffusion systems

The generalized fractional reaction–diffusion equations which exist in the form of noninteger order partial differential equations have now found wide application for illustrating important and useful physical phenomena, such as subdiffusive and superdiffusive scenarios. The space fractional derivatives are defined in the Riesz sense on the intervals 0 < α < 1 and 1 < α ≤ 2. We propose robust numerical techniques, such as a spectral representation of the fractional Laplacian operator in conjunction with the exponential time differencing method, and present the equivalent relationship between the Riesz fractional derivative and fractional Laplacian operator. We apply these techniques to numerically solve a range of chaotic processes, such as the Chua’s equations, Rössler system, Lorenz and Lorenz-type systems. Simulation results revealed various complex and spatiotemporal chaos, spiral chaos, intermittent chaos and spots patterns in two-dimensional space fractional reaction–diffusion problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Travelling fronts of the CO oxidation on Pd(111) with coverage-dependent diffusion

In this work, we study a surface reaction on Pd(111) crystals under ultra-high-vacuum conditions that can be modeled by two coupled reaction-diffusion equations. In the bistable regime, the reaction exhibits travelling fronts that can be observed experimentally using photo electron emission microscopy. The spatial profile of the fronts reveals a coverage-dependent diffusivity for one of the species. We propose a method to solve the nonlinear eigenvalue problem and compute the direction and the speed of the fronts based on a geometrical construction in phase-space. This method successfully captures the dependence of the speed on control parameters and diffusivities.

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

Entropically induced asymmetric passage times of charged tracers across corrugated channels

We analyze the diffusion of charged and neutral tracers suspended in an electrolyte embedded in a channel of varying cross section. Making use of systematic approximations, the diffusion equation governing the motion of tracers is mapped into an effective 1D equation describing the dynamics along the longitudinal axis of the channel where its varying-section is encoded as an effective entropic potential. This simplified approach allows us to characterize tracer diffusion under generic confinement by measuring their mean first passage time (MFPT). In particular, we show that the interplay between geometrical confinement and electrostatic interactions strongly affect the MFTP of tracers across corrugated channels hence leading to alternative means to control tracers translocation across charged pores. Finally, our results show that the MFPTs of a charged tracer in opposite directions along an asymmetric channel may differ We expect our results to be relevant for biological as well synthetic devices whose dynamics is controlled by the detection of diluted tracers.

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

Numerical model of the plasma formation at electron beam welding

The model of plasma formation in the keyhole in liquid metal as well as above the electron beam welding zone is described. The model is based on solution of two equations for the density of electrons and the mean electron energy. The mass transfer of heavy plasma particles (neutral atoms, excited atoms, and ions) is taken into account in the analysis by the diffusion equation for a multicomponent mixture. The electrostatic field is calculated using the Poisson equation. Thermionic electron emission is calculated for the keyhole wall. The ionization intensity of the vapors due to beam electrons and high-energy secondary and backscattered electrons is calibrated using the plasma parameters when there is no polarized collector electrode above the welding zone. The calculated data are in good agreement with experimental data. Results for the plasma parameters for excitation of a non-independent discharge are given. It is shown that there is a need to take into account the effect of a strong electric field near the keyhole walls on electron emission (the Schottky effect) in the calculation of the current for a non-independent discharge (hot cathode gas discharge). The calculated electron drift velocities are much bigger than the velocity at which current instabilities arise. This confirms the hypothesis for ion-acoustic instabilities, observed experimentally in previous research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A First-Passage Kinetic Monte Carlo method for reaction–drift–diffusion processes

Stochastic reaction–diffusion models are now a popular tool for studying physical systems in which both the explicit diffusion of molecules and noise in the chemical reaction process play important roles. The Smoluchowski diffusion-limited reaction model (SDLR) is one of several that have been used to study biological systems. Exact realizations of the underlying stochastic processes described by the SDLR model can be generated by the recently proposed First-Passage Kinetic Monte Carlo (FPKMC) method. This exactness relies on sampling analytical solutions to one and two-body diffusion equations in simplified protective domains. In this work we extend the FPKMC to allow for drift arising from fixed, background potentials. As the corresponding Fokker–Planck equations that describe the motion of each molecule can no longer be solved analytically, we develop a hybrid method that discretizes the protective domains. The discretization is chosen so that the drift–diffusion of each molecule within its protective domain is approximated by a continuous-time random walk on a lattice. New lattices are defined dynamically as the protective domains are updated, hence we will refer to our method as Dynamic Lattice FPKMC or DL-FPKMC. We focus primarily on the one-dimensional case in this manuscript, and demonstrate the numerical convergence and accuracy of our method in this case for both smooth and discontinuous potentials. We also present applications of our method, which illustrate the impact of drift on reaction kinetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Surface self-diffusion of silicon during high temperature annealing

The atomic-scale mechanisms driving thermally activated self-diffusion on silicon surfaces are investigated by atomic force microscopy. The evolution of surface topography is quantified over a large spatial bandwidth by means of the Power Spectral Density functions. We propose a parametric model, based on the Mullins-Herring (M-H) diffusion equation, to describe the evolution of the surface topography of silicon during thermal annealing. Usually, a stochastic term is introduced into the M-H model in order to describe intrinsic random fluctuations of the system. In this work, we add two stochastic terms describing the surface thermal fluctuations and the oxidation-evaporation phenomenon. Using this extended model, surface evolution during thermal annealing in reducing atmosphere can be predicted for temperatures above the roughening transition. A very good agreement between experimental and theoretical data describing roughness evolution and self-diffusion phenomenon is obtained. The physical origin and time-evolution of these stochastic terms are discussed. Finally, using this model, we explore the limitations of the smoothening of the silicon surfaces by rapid thermal annealing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

New discretization of Caputo–Fabrizio derivative

We derive a numerical approximation, namely L1–2 formula, to the Caputo–Fabrizio derivative by using a quadratic interpolation. Quadratic and cubic convergence rates are achieved for L1 and L1–2 formulas using Lagrange interpolation, respectively. We compute Caputo–Fabrizio derivatives of some known functions both theoretically and numerically. In addition, we solve non/linear sub-diffusion equations to test theoretical findings. Numerical results confirm the theoretically observed convergence rates.

97 MATHEMATICS AND COMPUTING↗