SEARCH · Engineering Papers
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.
Solution of Boltzmann and rate equations for electron distribution function and state populations in nonequilibrium MHD plasmas
Comparison of monte carlo, percus-yevick integral equations and convolution-hypernetted-chain methods for computing radial distribution functions of gaussian model
Estimating accurate Time-of-Failure (ToF) of a system is key in making the decisions that impact operational safety and optimize cost. In this context, it is interesting to note that different approaches have been explored to tackle the problem of estimating ToF. The difference is in part characterized by different definitions of the hazard zones. The conventional definition for the cumulative distribution function (CDF) calculation is assumed to have well-defined hazard zones, that is, hazard zones defined as a function of the system state trajectory. An alternate method suggests the use of hazard zones defined as a function of the system state at time , instead of hazard zones defined as a function of system state up to and including time k (Acuna and Orchard 2018, 2017). This paper explores these differences and their impact on ToF estimation. Results for the conventional CDF definition indicated that, (i) the cumulative distribution function is always an increasing function of time, even when realizations of the degradation process are not monotonic, (ii) the sum of all probabilities is always 1 and does not need to be normalized, and (iii) all probabilities are positive and less than or equal to 1. Similar results are not observed for CDF calculation with hazard zones defined as a function only of the system state at time . Results for ToF estimation using Acuna's definition differ, suggesting that there is an underlying assumption of independence in the hazard zone definition. Therefore, we present an alternate definition of hazard zone which guarantees the properties of a well-defined CDF with a more straightforward ToF definition.
Effect of ionization and recombination on the electron distribution function of a slightly ionized gas in the presence of an external magnetic and electric field
Our present understanding of the distribution functions of the solar wind ion and electron thermal/suprathermal populations will be reviewed covering selected theoretical and observational topics. Roughly half the review will be devoted to recent discoveries (since Solar Wind 7). Among recent results are those of the Ulysses mission. The SWICS experiment observed pick-up protons and alpha particles, and acceleration of these particles at interplanetary shocks. Positive ion phenomena observed by SWOOPS include enhanced proton-alpha particle differential streaming both: (1) at high latitudes; and (2) in the ecliptic beyond 1 AU in the vicinity of shocks. The SWOOPS positive ion observations also demonstrate relative constancy of the entropy per proton at high latitudes beyond about 2.5 AU. Double beaming in both protons and alphas is observed by SWOOPS on either side of the current sheet, but not within; a possible explanation is reconnection at the edge of the current sheet. SWOOPS observed bi-directional streaming of electrons beyond 2 AU not only within Coronal Mass Ejections, but also upstream of corotating interplanetary shocks. Latitudinal and radial gradients of Ulysses electron core and halo temperatures will be discussed, as will work in progress on the three dimensional (in velocity space) properties of ion distributions.
Electron distribution function for steady homogeneous plasma of arbitrary ionization degree under strong electric fields, discussing kinetic equations for elastic and nonelastic encounters
Results of an extensive investigation of probability distribution functions (pdfs) for Rayleigh-Benard convection, in hard turbulence regime, are presented. It is shown that the pdfs exhibit a high degree of internal universality. In certain cases this universality is established within two Kolmogorov scales of a boundary. A discussion of the factors leading to the universality is presented.
Results of an extensive investigation of probability distribution functions (pdf's) for Rayleigh-Benard convection, in the hard turbulence regime, are presented. It is seen that the pdf's exhibit a high degree of internal universality. In certain cases this universality is established within two Kolmogorov scales of a boundary. A discussion of the factors leading to universality is presented.
Our present understanding of the distribution functions of the solar wind ion and electron thermal/suprathermal populations are reviewed, covering selected theoretical and observational topics (half of material new since Solar Wind 7).
We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.
We study electron distribution functions in a diffusion region of antiparallel asymmetric reconnection by means of particle-in-cell simulations and analytical theory. At the electron stagnation point, the electron distribution comprises a crescent-shaped population and a core component. The crescent-shaped distribution is due to electrons coming from the magnetosheath toward the stagnation point and accelerated mainly by electric field normal to the current sheet. Only a part of magnetosheath electrons can reach the stagnation point and form the crescent-shaped distribution that has a boundary of a parabolic curve. The penetration length of magnetosheath electrons into the magnetosphere is derived. We expect that satellite observations can detect crescent-shaped electron distributions during magnetopause reconnection.
Ytterbium-doping in all-inorganic lead-halide perovskites (CsPb(Cl1-xBrx)3) generates novel properties including quantum cutting and narrow line emission, making these materials attractive spectral down-converters for solar photovoltaics. The relationship between this optical efficiency and the defect structure(s) associated with Yb3+ dopants within perovskites is not well understood. Various charge-neutral doping motifs have previously been proposed and studied computationally, including clusters involving two substitutional Yb3+ ions charge-compensated by a single local Pb2+ vacancy. Near-band-edge defect states associated with such motifs are believed to play an important mechanistic role in quantum cutting itself. Here, we report the results of X-ray absorption and X-ray total-scattering measurements on ytterbium-doped CsPbCl3. XANES shows that the dopant oxidation state is exclusively Yb3+, and a combination of Yb L3 and Pb L3 EXAFS shows that this Yb3+ substitutes exclusively at Pb2+ sites, where it adopts a pseudo-octahedral [YbCl6]3- coordination environment. Shell-by-shell fits to the data show a short Yb-Cl bond distance of 2.58 Å compared to the Pb-Cl bond distance of 2.83 Å. We confirm this finding by X-ray pair distribution function analysis, which also shows evidence of additional Pb2+ vacancy formation induced by Yb3+ doping. We evaluate whether this is the primary mechanism of charge compensation using simulated EXAFS and pair distribution function data for several computed defect structures. Together, these results resolve the local dopant structures and charge-compensation mechanisms in lanthanide-doped all-inorganic lead-halide perovskites, and, as such, significantly advance the understanding of structure-function relationships in this important class of materials.
An instrument which can rapidly measure plasma particle distribution functions has been developed based upon recent innovations in electrostatic analyzer design and position sensitive particle detection. The new analyzer uses a quadrispherical geometry, but has a completely uniform 360 deg fan-shaped field of view. The polar angular distribution of entering particles is spatially imaged onto a position sensitive detector at the annular exit aperture after a deflection through 90 deg. Several methods of position sensitive detection have been successfully used in conjunction with this analyzer. The simplest is individual channel multipliers spaced around the annular exit. Microchannel plate electron multipliers permit greater position resolution to be obtained, and a detector using microchannel plates followed by a resistive anode image converter obtains angular resolution of about one degree, i.e., 360 individual angle pixels. Instruments of this type were flown on a sounding rocket in early 1982 and will be included on the Giotto comet mission and the AMPTE ion release module (IRM).
Electron energy distribution functions were calculated in a U235 plasma at 1 atmosphere for various plasma temperatures and neutron fluxes. The distributions are assumed to be a summation of a high energy tail and a Maxwellian distribution. The sources of energetic electrons considered are the fission-fragment induced ionization of uranium and the electron induced ionization of uranium. The calculation of the high energy tail is reduced to an electron slowing down calculation, from the most energetic source to the energy where the electron is assumed to be incorporated into the Maxwellian distribution. The pertinent collisional processes are electron-electron scattering and electron induced ionization and excitation of uranium. Two distinct methods were employed in the calculation of the distributions. One method is based upon the assumption of continuous slowing and yields a distribution inversely proportional to the stopping power. An iteration scheme is utilized to include the secondary electron avalanche. In the other method, a governing equation is derived without assuming continuous electron slowing. This equation is solved by a Monte Carlo technique.
Here, we measure the Voronoi density probability distribution function (PDF) for both dark matter and halos in N-body simulations. For the dark matter, Voronoi densities represent the matter density field smoothed on a uniform mass scale, which approximates the Lagrangian density field. For halos, the Voronoi densities contain information about the local environment of each halo. We measure the halo virial masses, the total amount of dark matter within each halo Voronoi cell, and the halo Voronoi cell volumes, and we show how halo abundances depend on these three quantities. We then study the position-dependent Voronoi density PDF, measured within finite subregions of the Universe, using separate universe simulations. We demonstrate that the spatial variation of the position-dependent PDF is due to large-scale density fluctuations, indicating that the position-dependent PDF is a biased tracer of large-scale structure. We measure this bias for the dark matter, and interpret it as the bias of regions of the Lagrangian density field that are selected based on density. For the halos, this bias can be interpreted as a form of assembly bias. We present the mapping from late-time to early-time Voronoi density for each simulation dark matter particle, which is highly stochastic. We compare the median of this stochastic map with spherical collapse calculations and discuss challenges involved in modeling the evolution of the density field on these scales.
Transport properties and distribution function dynamics in one dimensional plasma models
Non-Maxwellian ion velocity distribution functions have been theoretically predicted and confirmed by observations, to occur at high latitudes. These distributions deviate from Maxwellian due to the combined effect of the E x B drift and ion-neutral collisions. At high altitude and/or for solar maximum conditions, the ion-to-neutral density ratio increases and, hence, the role of ion self-collisions becomes appreciable. A Monte Carlo simulation was used to investigate the behavior of O(+) ions that are E x B-drifting through a background of neutral O, with the effect of O(+) (Coulomb) self-collisions included. Wide ranges of the ion-to-neutral density ratio n(sub i)/n(sub n) and the electrostatic field E were considered in order to investigate the change of ion behavior with solar cycle and with altitude. For low altitudes and/or solar minimum (n(sub i)/n(sub n) less than or equal to 10(exp -5)), the effect of self-collisions is negligible. For higher values of n(sub i)/n(sub n), the effect of self-collisions becomes significant and, hence, the non-Maxwellian features of the O(+) distribution are reduced. The Monte Carlo results were compared to those that used simplified collision models in order to assess their validity. In general, the simple collision models tend to be more accurate for low E and for high n(sub i)/n(sub n).
In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.