Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DISTRIBUTION FUNCTIONS”

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 NTRS1972

The solution of the relaxation problem for the Boltzmann equation by the integral iteration method

The Boltzmann equation is considered in terms of the problem of relaxation of some initial distribution function which depends only on velocities, to Maxwell's distribution function. The Boltzmann equation is given for the relaxation problem in which the distribution function f(t, u, v) is time dependent and is also dependent on two other variables u and v (the velocities of rigid spherical molecules). An iteration process is discussed in which the velocity space u, v is subdivided into squares, the distribution function in each square being approximated by the second-order surface from the values of the distribution function at nine points. The set of all of these points forms a network of u, v values at the nodes of which the distribution function can be found.

Limar, Y. F.↗
DOE OSTI2020

Renormalization and matching for the Collins-Soper kernel from lattice QCD

The Collins-Soper kernel, which governs the energy evolution of transverse- momentum dependent parton distribution functions (TMDPDFs), is required to accurately predict Drell-Yan like processes at small transverse momentum, and is a key ingredient for extracting TMDPDFs from experiment. Earlier we proposed a method to calculate this kernel from ratios of the so-called quasi-TMDPDFs determined with lattice QCD, which are defined as hadronic matrix elements of staple-shaped Euclidean Wilson line operators. Here we provide the one-loop renormalization of these operators in a regularization-independent momentum subtraction (RI'/MOM) scheme, as well as the conversion factor from the RI'/MOM-renormalized quasi-TMDPDF to the \( \overline{\mathrm{MS}} \) scheme. We also propose a procedure for calculating the Collins-Soper kernel directly from position space correlators, which simplifies the lattice determination.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Conditions for the validity of unmagnetized-plasma theory in describing weakly magnetized plasmas

The application of the unmagnetized-plasma theory to the derivation of the expression for the dielectric tensor for weakly magnetized plasmas is examined. The analysis of Baldwin et. al. (1969) is employed to determine the conditions under which the magnetic effects on the dielectric tensor of a plasma are small. The relation between wave and plasma parameters and the plasma distribution function is studied. Axisymmetric distribution functions in magnetized plasma are considered. The derived conditions are analyzed in terms of resonance ellipses of the cyclotron maser theory. Examples are presented in which the conditions are applied to a nonrelativistic Maxwellian distribution and the dispersion and instability of Langmuir and Bernstein waves in weakly magnetized plasma. The data reveal that the conditions under which the magnetic field can be completely neglected depend on the forms of the distribution function.

Robinson, P. A.↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗
NASA NTRS1966

X-ray analysis of fatigue damage in copper

A Warren-Averbach x-ray line profile analysis was applied to, broadened x-ray diffraction peaks from copper deformed in fatigue. The copper specimens were fatigued by four point bending at various strain amplitudes, and measurements were made at various fractions of the total fatigue life. The analysis results in an estimation of 1) an average coherently diffracting domain size in a direction normal to the diffracting planes; and 2) a strain distribution function. The strain distribution function is in terms of the root mean square strain, <ϵ^(2)> where ϵ is averaged over a given distance at all points in the diffracting crystals. This rms strain is expressed as a function of averaging distance for a direction normal to the diffracting planes. The rms strains were found to be a function of the fatigue strain amplitude level, increasing with amplitude and reaching a maximum limiting value during the first 2% of the fatigue life. The strains were normally larger in the 4lO0 direction compared to the <111> direction and the absolute values were quite small. These results may be explained by assuming the dislocations are present in the form of dipoles because the lattice strain around a dipole is much less than around single dislocations. Conclusions arrived at in this analysis are in agreement with current theories of fatigue damage in metals. The result that maximum hardening is achieved during the initial few percent of the fatigue life is confirmed by many other investigators using entirely different analytical techniques. It is felt that the structure of the material is changing during this initial hardening stage to an eventual steady state characteristic of the testing conditions. A unique substructure is developed which thereafter persists, and the remainder of the fatigue test is taken up by crack propagation through this substructure until ultimate failure occurs by fracture. The formation of this steady state structure is dependent on the maximum fatigue strain amplitude and on the frequency and temperature of cycling for a particular material.

COPPER↗
NASA NTRSDate not supplied

Assessing the Viability of Using GEOS-Forecast Product for Landslides Forecasting: A Step Toward Early Warning System

Landslides across the globe are mostly triggered by extreme rainfall events affecting infrastructure, transportation and livelihoods. The risks are rarely quantified due to lack of data, analytical skills and limited modeling techniques. Knowledge of local to global scale landslide risks provides communities and national agencies the ability to adapt disaster management practices to mitigate and recover from these hazards. In order to minimize the risks and improve characterization of community resilience to landslides, it is vital to have reliable information about the factors triggering landslides such as rainfall, well ahead in time. Forecasting potential landslide activity and impacts can be achieved through reliable precipitation forecast models. However, it is challenging because of the temporal and spatial variability of precipitation, an important factor in triggering landslides. Evaluation of the precipitation field, associated errors, and sampling uncertainties is integral for development of efficient and reliable landslide forecasting and early warning system. This study develops a methodology to assess the viability of using a precipitation field provided by a global model and its potential integration in the landslide forecasting system. The study focuses on the comparison between the IMERG (Integrated Multi-satellitE Retrievals for Global Precipitation Mission) and GEOS (NASA Goddard Earth Observing System)-Forecast product over contiguous United States (CONUS). GEOS model assimilates new observations every 6 hours, at 00, 06, 12, and 18 UTC. The framework is tested on the GEOS-Forecast Model initialized at 00 UTC using daily IMERG early product as reference using both categorical and continuous statistics. The categorical statistics includes the probability of detection (POD), success ratio (SR), critical success index (CSI), and the hit bias. Continuous statistics such as correlation, normalized standard deviation, and root-mean-square error are also evaluated. Overall, GEOS-Forecast precipitation field over the analysis period (~1 year) show underestimation with respect to IMERG early for the daily accumulated rainfall. However, the probability distribution function and cumulative distribution function of both show similar patterns. In terms of correlations, POD, SR, CSI, hit bias, the performance varies with respect to the rainfall threshold used.

Sana Khan↗
NASA NTRS2018

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

The triple Langmuir probe is an electrostatic probe in which three probe tips collect current when inserted into a plasma. The triple probe differs from a simple single Langmuir probe in the nature of the voltage applied to the probe tips. In the single probe, a swept voltage is applied to the probe tip to acquire a waveform showing the collected current as a function of applied voltage (I-V curve). In a triple probe three probe tips are electrically coupled to each other with constant voltages applied between each of the tips. The voltages are selected such that they would 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. Under the assumption of a Maxwellian plasma, one can determine the time-varying plasma temperature T(sub e)(t) and number density n(sub e)(t) from the applied voltage levels and the time-histories of the collected currents. In the present paper we examine the theory of triple probe operation, specifically focusing on the assumption of a Maxwellian plasma. Triple probe measurements have been widely employed for a number of pulsed and timevarying plasmas, including pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. While the equilibrium assumption may be justified for some applications, it is unlikely that it is fully justifiable for all pulsed and time-varying plasmas or for all times during the pulse of a plasma device. To examine a simple non-equilibrium plasma case, 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. We compare the distribution function that was assumed a priori with the distribution function one would infer when applying standard triple probe theory to analyze the collected currents. For the assumed class of non-Maxwellian distribution functions this serves to illustrate the effect a non-Maxwellian plasma would have on results interpreted using the equilibrium triple 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.↗
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.↗