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 937 records · Page 52

A horseshoe-shaped ring of diffuse emission detected at 1.4 GHz

We identify a horseshoe-shaped ring (HSR) of diffuse emission in J1407+0453 from the Faint Images of Radio Sky at Twenty-cm (FIRST) survey using the Very Large Array telescope at 1.4 GHz. An optical galaxy, SDSS J140709.01+045302.1, is present near the limb of the HSR of J1407+0453, with a spectroscopic redshift ofz = 0.13360. The total extent of the source, including the diffuse emission of J1407+0453, is 65 arcsec (with a physical extent of 160 kpc), whereas the diameter of the HSR is approximately 10 arcsec (25 kpc). The flux density of the HSR is ∼47 mJy at 1400 MHz, whereas the flux densities of the whole diffuse emission of J1407+0453 at 1400 MHz and 150 MHz are 172 mJy and 763 mJy, respectively. We measure the radio luminosity of HSR J1407+0453 as 1.94 × 10 24 W Hz −1 , with a spectral index,α 150 1400 = −0.67. The black hole mass of J1407+0453 is 5.8 × 10 8 M ⊙ . We compare the HSR of diffuse emission of J1407+0453 with other discovered diffuse circular sources. The possible formation scenarios for J1407+0453 are discussed, so as to understand the nature of the source. We present a spectral index map of source J147+0453 to study the spectral properties of the source.

Astronomy & Astrophysics↗

Analysis and experimental study on formation conditions of large-scale barrier-free diffuse atmospheric pressure air plasmas in repetitive pulse mode

Atmospheric air diffuse plasmas have enormous application potential in various fields of science and technology. Without dielectric barrier, generating large-scale air diffuse plasmas is always a challenging issue. This paper discusses and analyses the formation mechanism of cold homogenous plasma. It is proposed that generating stable diffuse atmospheric plasmas in open air should meet the three conditions: high transient power with low average power, excitation in low average E-field with locally high E-field region, and multiple overlapping electron avalanches. Accordingly, an experimental configuration of generating large-scale barrier-free diffuse air plasmas is designed. Based on runaway electron theory, a low duty-ratio, high voltage repetitive nanosecond pulse generator is chosen as a discharge excitation source. Using the wire-electrodes with small curvature radius, the gaps with highly non-uniform E-field are structured. Experimental results show that the volume-scaleable, barrier-free, homogeneous air non-thermal plasmas have been obtained between the gap spacing with the copper-wire electrodes. The area of air cold plasmas has been up to hundreds of square centimeters. The proposed formation conditions of large-scale barrier-free diffuse air plasmas are proved to be reasonable and feasible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diffusivity anomaly in modified Stillinger-Weber liquids

By modifying the tetrahedrality (the strength of the three body interactions) in the well-known Stillinger-Weber model for silicon, we study the diffusivity of a series of model liquids as a function of tetrahedrality and temperature at fixed pressure. Previous work has shown that at constant temperature, the diffusivity exhibits a maximum as a function of tetrahedrality, which we refer to as the diffusivity anomaly, in analogy with the well-known anomaly in water upon variation of pressure at constant temperature. We explore to what extent the structural and thermodynamic changes accompanying changes in the interaction potential can help rationalize the diffusivity anomaly, by employing the Rosenfeld relation between diffusivity and the excess entropy (over the ideal gas reference value), and the pair correlation entropy, which provides an approximation to the excess entropy in terms of the pair correlation function. We find that in the modified Stillinger-Weber liquids, the Rosenfeld relation works well above the melting temperatures but exhibits deviations below, with the deviations becoming smaller for smaller tetrahedrality. Further we find that both the excess entropy and the pair correlation entropy at constant temperature go through maxima as a function of the tetrahedrality, thus demonstrating the close relationship between structural, thermodynamic, and dynamical anomalies in the modified Stillinger-Weber liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnonic band gaps in two-dimension magnonic crystals with diffuse interfaces

In this paper, the plane wave method is extended to include the diffuse interface in the calculation of the dispersion of spin waves in two-dimension magnonic crystals. The diffuse interfaces with linear and sinusoidal profiles of variation in the spontaneous magnetization and exchange constant are considered and the effects of the thicknesses and profiles of diffuse interfaces on the magnonic band gaps are investigated. The results show that the thicknesses and profiles of diffuse interfaces are clearly seen to play a significant role in determining the size and position of the magnonic band gaps in the both square and triangular lattices in the exchange interaction regime. The smooth (linear or sinusoidal) interface does not lead to disappearance of the band gaps, instead it may lead to larger band gaps than those in the model with sharp (infinitely thin) diffuse interface under certain conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Generation of large-scale, barrier-free diffuse plasmas in air at atmospheric pressure using array wire electrodes and nanosecond high-voltage pulses

This paper introduces a method to generate large-scale diffuse plasmas by using a repetition nanosecond pulse generator and a parallel array wire-electrode configuration. We investigated barrier-free diffuse plasmas produced in the open air in parallel and cross-parallel array line-line electrode configurations. We found that, when the distance between the wire-electrode pair is small, the discharges were almost extinguished. Also, glow-like diffuse plasmas with little discharge weakening were obtained in an appropriate range of line-line distances and with a cathode-grounding cross-electrode configuration. As an example, we produced a large-scale, stable diffuse plasma with volumes as large as 18 × 15 × 15 cm{sup 3}, and this discharge region can be further expanded. Additionally, using optical and electrical measurements, we showed that the electron temperature was higher than the gas temperature, which was almost the same as room temperature. Also, an array of electrode configuration with more wire electrodes had helped to prevent the transition from diffuse discharge to arc discharge. Comparing the current waveforms of configurations with 1 cell and 9 cells, we found that adding cells significantly increased the conduction current and the electrical energy delivered in the electrode gaps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-dimensional numerical simulation of boron diffusion for pyramidally textured silicon

Multidimensional numerical simulation of boron diffusion is of great relevance for the improvement of industrial n-type crystalline silicon wafer solar cells. However, surface passivation of boron diffused area is typically studied in one dimension on planar lifetime samples. This approach neglects the effects of the solar cell pyramidal texture on the boron doping process and resulting doping profile. In this work, we present a theoretical study using a two-dimensional surface morphology for pyramidally textured samples. The boron diffusivity and segregation coefficient between oxide and silicon in simulation are determined by reproducing measured one-dimensional boron depth profiles prepared using different boron diffusion recipes on planar samples. The established parameters are subsequently used to simulate the boron diffusion process on textured samples. The simulated junction depth is found to agree quantitatively well with electron beam induced current measurements. Finally, chemical passivation on planar and textured samples is compared in device simulation. Particularly, a two-dimensional approach is adopted for textured samples to evaluate chemical passivation. The intrinsic emitter saturation current density, which is only related to Auger and radiative recombination, is also simulated for both planar and textured samples. The differences between planar and textured samples are discussed.

36 MATERIALS SCIENCE↗

Diffusion anisotropy of poor metal solute atoms in hcp-Ti

Atom migration mechanisms influence a wide range of phenomena: solidification kinetics, phase equilibria, oxidation kinetics, precipitation of phases, and high-temperature deformation. In particular, solute diffusion mechanisms in α-Ti alloys can help explain their excellent high-temperature behaviour. The purpose of this work is to study self- and solute diffusion in hexagonal close-packed (hcp)-Ti, and its anisotropy, from first-principles using the 8-frequency model. The calculated diffusion coefficients show that diffusion energy barriers depend more on bonding characteristics of the solute rather than the size misfit with the host, while the extreme diffusion anisotropy of some solute elements in hcp-Ti is a result of the bond angle distortion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Indium and gallium diffusion through zirconia in the TiN/ZrO{sub 2}/InGaAs stack

Angle-resolved X-ray Photoelectron Spectroscopy (ARXPS) was applied to the TiN/ZrO{sub 2}/InGaAs stack to assess its thermal stability. Through a robust ARXPS analysis, it was possible to observe subtle effects such as the thermally induced diffusion of substrate atomic species (In and Ga) through the dielectric layer. The detailed characterization of the film structure allowed for assessing the depth profiles of the diffused atomic species by means of the scenarios-method. Since the quantification for the amount of diffused material was done at different temperatures, it was possible to obtain an approximate value of the activation energy for the diffusion of indium through zirconia. The result is very similar to the previously reported values for indium diffusion through alumina and through hafnia.

36 MATERIALS SCIENCE↗

Reaction rates for reaction-diffusion kinetics on unstructured meshes

The reaction-diffusion master equation is a stochastic model often utilized in the study of biochemical reaction networks in living cells. It is applied when the spatial distribution of molecules is important to the dynamics of the system. A viable approach to resolve the complex geometry of cells accurately is to discretize space with an unstructured mesh. Diffusion is modeled as discrete jumps between nodes on the mesh, and the diffusion jump rates can be obtained through a discretization of the diffusion equation on the mesh. Reactions can occur when molecules occupy the same voxel. Here in this paper, we develop a method for computing accurate reaction rates between molecules occupying the same voxel in an unstructured mesh. For large voxels, these rates are known to be well approximated by the reaction rates derived by Collins and Kimball, but as the mesh is refined, no analytical expression for the rates exists. We reduce the problem of computing accurate reaction rates to a pure preprocessing step, depending only on the mesh and not on the model parameters, and we devise an efficient numerical scheme to estimate them to high accuracy. We show in several numerical examples that as we refine the mesh, the results obtained with the reaction-diffusion master equation approach those of a more fine-grained Smoluchowski particle-tracking model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasilinear diffusion coefficients in a finite Larmor radius expansion for ion cyclotron heated plasmas

In this study, a reduced model of quasilinear velocity diffusion by a small Larmor radius approximation is derived to couple the Maxwell’s equations and the Fokker Planck equation self-consistently for the ion cyclotron range of frequency waves in a tokamak. The reduced model ensures the important properties of the full model by Kennel-Engelmann diffusion, such as diffusion directions, wave polarizations, and H-theorem. The kinetic energy change (Wdot ) is used to derive the reduced model diffusion coefficients for the fundamental damping (n = 1) and the second harmonic damping (n = 2) to the lowest order of the finite Larmor radius expansion. The quasilinear diffusion coefficients are implemented in a coupled code (TORIC-CQL3D) with the equivalent reduced model of the dielectric tensor. We also present the simulations of the ITER minority heating scenario, in which the reduced model is verified within the allowable errors from the full model results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗

DIFFUSE HARD X-RAY EMISSION IN STARBURST GALAXIES AS SYNCHROTRON FROM VERY HIGH ENERGY ELECTRONS

The origin of the diffuse hard X-ray (2-10 keV) emission from starburst galaxies is a long-standing problem. We suggest that synchrotron emission of 10-100 TeV electrons and positrons (e {sup {+-}}) can contribute to this emission, because starbursts have strong magnetic fields. We consider three sources of e {sup {+-}} at these energies: (1) primary electrons directly accelerated by supernova remnants, (2) pionic secondary e {sup {+-}} created by inelastic collisions between cosmic ray (CR) protons and gas nuclei in the dense interstellar medium of starbursts, and (3) pair e {sup {+-}} produced between the interactions between 10 and 100 TeV {gamma}-rays and the intense far-infrared (FIR) radiation fields of starbursts. We create one-zone steady-state models of the CR population in the Galactic center (R {<=} 112 pc), NGC 253, M82, and Arp 220's nuclei, assuming a power-law injection spectrum for electrons and protons. We consider different injection spectral slopes, magnetic field strengths, CR acceleration efficiencies, and diffusive escape times, and include advective escape, radiative cooling processes, and secondary and pair e {sup {+-}}. We compare these models to extant radio and GeV and TeV {gamma}-ray data for these starbursts, and calculate the diffuse synchrotron X-ray and inverse Compton (IC) luminosities of these starbursts in the models which satisfy multiwavelength constraints. If the primary electron spectrum extends to {approx}PeV energies and has a proton/electron injection ratio similar to the Galactic value, we find that synchrotron emission contributes 2%-20% of their unresolved, diffuse hard X-ray emission. However, there is great uncertainty in this conclusion because of the limited information on the CR electron spectrum at these high energies. IC emission is likewise a minority of the unresolved X-ray emission in these starbursts, from 0.1% in the Galactic center to 10% in Arp 220's nuclei, with the main uncertainty being the starbursts' magnetic field. We also model generic starbursts, including submillimeter galaxies, in the context of the FIR-X-ray relation, finding that anywhere between 0% and 16% of the total hard X-ray emission is synchrotron for different parameters, and up to 2% in the densest starbursts assuming an E {sup -2.2} injection spectrum and a diffusive escape time of 10 Myr (E/3 GeV){sup -1/2} (h/100 pc). Neutrino observations by IceCube and TeV {gamma}-ray data from HESS, VERITAS, and CTA can further constrain the synchrotron X-ray emission of starbursts. Our models do not constrain the possibility of hard, second components of primary e {sup {+-}} from sources like pulsars in starbursts, which could enhance the synchrotron X-ray emission further.

79 ASTRONOMY AND ASTROPHYSICS↗

The diffuse gamma-ray flux associated with sub-PEV/PEV neutrinos from starburst galaxies

One attractive scenario for the excess of sub-PeV/PeV neutrinos recently reported by IceCube is that they are produced by cosmic rays in starburst galaxies colliding with the dense interstellar medium. These proton-proton (pp) collisions also produce high-energy gamma rays, which finally contribute to the diffuse high-energy gamma-ray background. We calculate the diffuse gamma-ray flux with a semi-analytic approach and consider that the very high energy gamma rays will be absorbed in the galaxies and converted into electron-positron pairs, which then lose almost all of their energy through synchrotron radiation in the strong magnetic fields in the starburst region. Since the synchrotron emission goes into energies below GeV, this synchrotron loss reduces the diffuse high-energy gamma-ray flux by a factor of about two, thus leaving more room for other sources to contribute to the gamma-ray background. For an E{sub ν}{sup −2} neutrino spectrum, we find that the diffuse gamma-ray flux contributes about 20% of the observed diffuse gamma-ray background in the 100 GeV range. However, for a steeper neutrino spectrum, this synchrotron loss effect is less important, since the energy fraction in absorbed gamma rays becomes lower.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffuse γ-ray emission from galactic pulsars

Millisecond pulsars (MSPs) are old fast-spinning neutron stars that represent the second most abundant source population discovered by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi). As guaranteed γ-ray emitters, they might contribute non-negligibly to the diffuse emission measured at high latitudes by Fermi-LAT (i.e., the Isotropic Diffuse γ-Ray Background (IDGRB)), which is believed to arise from the superposition of several components of galactic and extragalactic origin. Additionally, γ-ray sources also contribute to the anisotropy of the IDGRB measured on small scales by Fermi-LAT. In this manuscript we aim to assess the contribution of the unresolved counterpart of the detected MSPs population to the IDGRB and the maximal fraction of the measured anisotropy produced by this source class. To this end, we model the MSPs' spatial distribution in the Galaxy and the γ-ray emission parameters by considering observational constraints coming from the Australia Telescope National Facility pulsar catalog and the Second Fermi-LAT Catalog of γ-ray pulsars. By simulating a large number of MSP populations through a Monte Carlo simulation, we compute the average diffuse emission and the anisotropy 1σ upper limit. We find that the emission from unresolved MSPs at 2 GeV, where the peak of the spectrum is located, is at most 0.9% of the measured IDGRB above 10° in latitude. The 1σ upper limit on the angular power for unresolved MSP sources turns out to be about a factor of 60 smaller than Fermi-LAT measurements above 30°. Our results indicate that this galactic source class represents a negligible contributor to the high-latitude γ-ray sky and confirm that most of the intensity and geometrical properties of the measured diffuse emission are imputable to other extragalactic source classes (e.g., blazars, misaligned active galactic nuclei, or star-forming galaxies). Nevertheless, because MSPs are more concentrated toward the galactic center, we expect them to contribute significantly to the γ-ray diffuse emission at low latitudes. Because, along the galactic disk, the population of young pulsars overcomes in number that of MSPs, we compute the γ-ray emission from the whole population of unresolved pulsars, both young and millisecond, in two low-latitude regions: the inner Galaxy and the galactic center.

79 ASTRONOMY AND ASTROPHYSICS↗

A Herschel/HIFI legacy survey of HF and H{sub 2}O IN the galaxy: Probing diffuse molecular cloud chemistry

We combine Herschel observations for a total of 12 sources to construct the most uniform survey of HF and H{sub 2}O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-infrared (HIFI) allows us to compare the HF and H{sub 2}O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H{sub 2}O velocity distributions follow each other almost perfectly and establish that HF and H{sub 2}O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H{sub 2} in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H{sub 2} in our sample, we find that the N(H{sub 2}O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H{sub 2}O might be used as a tracer of H{sub 2}—within a factor of 2.5—in the diffuse interstellar medium (ISM). We show that the measured factor of ∼2.5 variation around the median is driven by true local variations in the H{sub 2}O abundance relative to H{sub 2} throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H{sub 2}O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse ISM (n{sub H} ⩽ 1000 cm{sup −3}).

79 ASTRONOMY AND ASTROPHYSICS↗