Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Some axisymmetric self-similar galaxy models

Distribution functions are presented for some axisymmetric self-similar galaxy models with flat rotation curves. In particular, we record velocity distribution functions for the logarithmic potentials studied by Richstone. These distribution functions depend on two integrals only. Other distribution functions (possibly involving more or other integrals) associated with the same mass model differ only by an additive function that integrates over velocities to yield zero density everywhere. However, a given distribution function can be associated with different mass models.

Miller, R. H.↗
NASA NTRS1992

Accelerated cometary ions observed downstream of the Comet Halley bow shock

Results are presented of fluxes of energetic ions with energies exceeding 100 keV that were detected upstream of the bow shock of Comet Halley by the Tuende instrument on board the VEGA 1 spacecraft. Downstream of the shock, ion fluxes in the energy range 100 to 180 keV were observed. The measured ion fluxes were transformed into distribution functions in the solar wind frame using a variety of assumptions concerning the energy dependence of the distribution function and the identity of the ion species. The derived distribution function upstream of the shock falls off steeply with energy between 100 and 150 keV, with an effective temperature of about 7 keV or spectral index of about -15. The distribution function increases with decreasing cometocentric distance, on average, reaching a maximum at the bow shock. The measured distribution functions are compared with those obtained by similar instruments on Giotto and ICE as well as with the predictions of several theoretical models that employ different acceleration mechanisms.

Kecskemety, K.↗
DOE OSTI2022

Mapping variations of redshift distributions with probability integral transforms

We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Representation of the Geosynchronous Plasma Environment in Spacecraft Charging Calculations

Historically, our ability to predict and postdict spacecraft surface charging has been limited by the characterization of the plasma environment. One difficulty lies in the common practice of fitting the plasma data to a Maxwellian or Double Maxwellian distribution function, which may not represent the data well for charging purposes. We use electron and ion flux spectra measured by the Los Alamos National Laboratory (LANL) Magnetospheric Plasma Analyzer (MPA) to examine how the use of different spectral representations of the charged particle environment in computations of spacecraft potentials during magnetospheric substorms affects the accuracy of the results. We calculate the spacecraft potential using both the measured fluxes and several different fits to these fluxes. These measured fluxes have been corrected for the difference between the measured and calculated potential. The potential computed using the measured fluxes and the best available material properties of graphite carbon, with a secondary electron escape fraction of 81%, is within a factor of three of the measured potential for 87% of the data. Potentials calculated using a Kappa function fit to the incident electron flux distribution function and a Maxwellian function fit to the incident ion flux distribution function agree with measured potentials nearly as well as do potentials calculated using the measured fluxes. Alternative spectral representations gave less accurate estimates of potential. The use of all the components of the net flux, along with spacecraft specific average material properties, gives a better estimate of the spacecraft potential than the high energy flux alone.

Davis, V. A.↗
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Probabilistic Life and Reliability Analysis of Model Gas Turbine Disk

In 1939, W. Weibull developed what is now commonly known as the "Weibull Distribution Function" primarily to determine the cumulative strength distribution of small sample sizes of elemental fracture specimens. In 1947, G. Lundberg and A. Palmgren, using the Weibull Distribution Function developed a probabilistic lifing protocol for ball and roller bearings. In 1987, E. V. Zaretsky using the Weibull Distribution Function modified the Lundberg and Palmgren approach to life prediction. His method incorporates the results of coupon fatigue testing to compute the life of elemental stress volumes of a complex machine element to predict system life and reliability. This paper examines the Zaretsky method to determine the probabilistic life and reliability of a model gas turbine disk using experimental data from coupon specimens. The predicted results are compared to experimental disk endurance data.

Holland, Frederic A.↗
NASA NTRS1972

A method for direct numerical integration of the Boltzmann equation

The principal difficulties in numerical solution of the Boltzmann equation are considered. The study is aimed at formulating a numerical solution in such a manner that it contains a minimum amount of excess information at the distribution function level. It is pointed out that the accurate calculation of the distribution function at each point in phase space requires a tremendous number of operations, due to the necessity of solving five-fold quadratures in the collision integral. This results in the operational memory of the digital computer being insufficient to store all the data on the distribution functions at the necessary points in phase space. An algorithm is constructed involving successive iterations of the Boltzmann equation which does not require storage of each step of the new distribution function.

Cheremisin, F. G.↗
NASA NTRS1984

Correlation of auroral hiss and upward electron beams near the polar cusp

Data were obtained from the DE-1 high-altitude plasma instrument (HAPI) and plasma wave instrument (PWI) during outbound passes through the polar cusp near local noon. The observed distribution functions of electron beams are fitted by drifting Maxwellian functions and the observed distribution functions of hot background electrons by isotropic Maxwellian functions. In addition, the cold plasma density is inferred from knowledge of the electron plasma frequency and the measured density of the warm plasma, including the electron beam distribution. The empirically fitted plasma parameters, including density, temperature and drifting energy, are used to solve the linear dispersion equation for the resulting whistler mode emissions. Because the whistler mode becomes quasi-electrostatic for wave-normal angles near the resonance cone, the electrostatic approximation is used for the whistler mode dispersion relation. The results of wave instability analyses are then compared with the wave observations. A ray tracing of cusp hiss emission is conducted to locate the wave source region (at about one earth-radius).

Lin, C. S.↗
NASA NTRS1993

The Generation of Oblique Magnetosonic Waves

One of the outstanding issues regarding the excitation of magnetosonic waves has been the observational evidence that obliquely propagating waves are dominant at comets and planetary foreshocks despite the predictions of linear theory that maximum growth occurs at parallel propagation. To address this issue, we have conducted a detailed linear theory using a beam-ring distribution function. The results have shown that such distribution functions are associated with four separate instabilities. Two of these instabilities are similar to the right hand resonant ion/ion and the non-resonant instabilities which are also present in the case of a field aligned beam. The other two instabilities are associated with the presence of ring part of the distribution function. One of these instabilities excites magnetosonic waves with maximum growth in the oblique directions. The other excites Alfven waves with maximum growth in the oblique directions. The importance of the former instability is that it may explain the oblique nature of the magnetosonic waves observed at planetary foreshocks and comets. In order to understand the nonlinear properties of these various instabilities, we have also conducted 2-D hybrid (particle ions, fluid electrons) simulations. We have found that when the beam density is sufficiently large so that the non-resonant instability has the largest growth rate, these waves dominate the initial wave power in the system. At later times, however, the obliquely propagating magnetosonic waves become dominant. Another important finding was that when the two instabilities which excite magnetosonic waves have the largest growth rates, the system is dominated by the obliquely propagating magnetosonic waves. This is despite the fact that the largest growth rate for one of the instabilities occurs in the parallel direction. The exact cause of this is not currently understood and is under investigation.

Source record↗
NASA NTRS1976

An evaluation of the method for determining the Whitham F-function using distributions of downwash and sidewash angles

The method of computing the Whitham F function using distributions of downwash and sidewash angles was evaluated with two different models. F functions which were calculated for a half angle cone cylinder at M infinites = 2.01, using theoretically and experimentally derived flow angles, show that the method is sensitive to small inaccuracies in the measured flow angles. An oblique wing transport model was tested at 0 deg angle of attack at M infinitely = 2.01. In this test, two different probes were used at two different distances from the model. The pressure signature derived from the F function was extrapolated and compared to the pressure signature measured at the distance of 0.87 body lengths with the static pressure probe. The agreement between the two pressure signatures was poor due to the many inaccuracies involved in using a probe designed to measure flow angularity.

Mendoza, J. P.↗
NASA NTRS2018

Behavior of Langmuir Probes in Non-Equilibrium Plasmas

Langmuir probes are diagnostic tools used to determine electron temperature, number density, and plasma potential. Irving Langmuir first used an electrostatic probe in the 1920s to find these characteristics in ionized gases. Single, double, and triple Langmuir probes are commonly used in plasma diagnostics because of their relative simplicity. In the single probe, a swept voltage is applied between the probe tip and circuit common to acquire a waveform showing the collected current as a function of applied voltage. A double Langmuir probe consists of two tips, both inserted into the plasma, with a voltage applied between them. As this voltage is swept, a current-voltage characteristic is measured. In a triple probe three probe tips are electrically coupled to each other with constant non-swept voltages applied between each of the tips. The voltages are selected to represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Triple Langmuir probe measurements have been widely employed for various types of plasmas, including pulsed and time-varying plasmas such as those seen in pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. The typical Langmuir probe analysis for determining electron temperature and number density of the plasma (for a single, double, or triple Langmuir probe) includes an assumption that the plasma is in thermal equilibrium. While the this assumption may be justified for some applications, it is unlikely that it is fully justifiable for pulsed and time-varying plasmas or for the entire time a plasma device is in use. In the present work, we model the responses of Langmuir probes as they are inserted into a range of simple equilibrium and non-equilibrium plasmas. We return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. Comparing the distribution function that was assumed a priori with the plasma density and temperature one would infer when applying standard probe theory to analyze the collected currents serves to illustrate the effect a non- Maxwellian plasma would have on results interpreted using the equilibrium probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗
DOE OSTI2020

ReaxFF molecular dynamics simulations of electrolyte–water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal–chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

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

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Propagation Pattern for Moment Representation of the Lattice Boltzmann Method

A propagation pattern for the moment representation of the regularized lattice Boltzmann method (LBM) in three dimensions is presented. Using effectively lossless compression, the simulation state is stored as a set of moments of the lattice Boltzmann distribution function, instead of the distribution function itself. An efficient cache-aware propagation pattern for this moment representation has the effect of substantially reducing both the storage and memory bandwidth required for LBM simulations. This article extends recent work with the moment representation by expanding the performance analysis on central processing unit (CPU) architectures, considering how boundary conditions are implemented, and demonstrating the effectiveness of the moment representation on a graphics processing unit (GPU) architecture.

42 ENGINEERING↗
NASA NTRS1973

TLC determination of functionality in prepolymers

Application of thin-layer chromatographics provides rapid qualitative determination of functional distribution in experimental prepolymer. Functionality distribution is of fundamental importance for it determines; (1) manner in which given carboxyl-terminated prepolymer will cure and (2) physical properties of resulting product.

Potts, J. E., Jr.↗
NASA NTRS1991

Particle motion in the tail current sheet

Theory of particle motion in current sheets is reviewed. For small, approximately constant normal magnetic field, Bz, particles oscillate about the current sheet and 'live' within the sheet for one-half gyroperiod based on Bz. This lifetime replaces the mean collision time in the Lorentzian conductivity and thus gives rise to the concept of an inertial (or gyro-) conductivity. A substorm model by Coroniti utilizes this conductivity to allow reconnection to proceed without anomalous processes, due to wave-particle interactions. Chaotic particle orbits may at times be important to the dynamics, depending on parameters such as particle energy, current sheet thickness, and field line curvature. A current sheet model with neutral line predicts a ridge structure and asymmetries in the distribution function. Ion distributions near the plasma sheet boundary layer, during the CDAW 6 interval, are consistent with the model predictions. In recent studies by Mitchell et al. and Williams et al., the major current carriers during the growth phase of a substorm were found to be adiabatic electrons not more than 1 keV, but just before a current disruption event, the tail current was mainly carried by energetic ions undergoing current sheet oscillation.

Speiser, T. W.↗
NASA NTRS1991

The modified plasma dispersion function

The modified plasma dispersion function (MPDF), based on the generalized Lorentzian (kappa) particle distribution function, is introduced, and a comprehensive set of graphs of the real and imaginary parts of the MPDF is presented. For any positive integral value of kappa, MPDF is calculated in closed form as a finite series. It is demonstrated how the MPDF approaches the plasma dispersion function in the limit as kappa yields infinity, a result to be expected since the kappa distribution function formally approaches the Maxwellian as kappa yields infinity. It is concluded that the MPDF can provide a tool in studying microinstabilities in plasmas when the particle distribution function is not only the standard generalized Lorentzian, but also of the Lorentzian type, including the loss-cone, bi-Lorentzian, and product bi-Lorentzian distributions.

Summers, Danny↗