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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

A strictly Markovian expansion for plasma turbulence theory

The collision operator that appears in the equation of motion for a particle distribution function that has been averaged over an ensemble of random Hamiltonians is non-Markovian. It is non-Markovian in that it involves a propagated integral over the past history of the ensemble averaged distribution function. All formal expansions of this nonlinear collision operator to date preserve this non-Markovian character term by term yielding an integro-differential equation that must be converted to a diffusion equation by an additional approximation. In this note we derive an expansion of the collision operator that is strictly Markovian to any finite order and yields a diffusion equation as the lowest non-trivial order. The validity of this expansion is seen to be the same as that of the standard quasi-linear expansion.

Jones, F. C.↗
NASA NTRS1988

Molecular tori in Seyfert galaxies - Feeding the monster and hiding it

The principal properties of the tori of gas which surround Seyfert nuclei are discussed. The internal state of the clouds and their size distribution function are examined, and it is shown that the Jeans mass scale results in clouds which are individually sufficiently opaque to block out the nucleus, and that the balance of processes which controls their size distribution function also forces the covering factor to be of the order of or greater than unity. Where the gravitational potential is dominated by stars, cloud-cloud collisions keep the molecular clouds close to the equatorial plane. Stirring by stellar processes is never strong enough to compete with collisional losses. The position of the inner edge of the torus is determined by a balance between the inward flow of clouds and the rate at which the nuclear continuum can evaporate them.

Krolik, Julian H.↗
NASA NTRSDate not supplied

Inter-Calibrating the CERES SW Instrument Fluxes Using Deep Convective Clouds(DCC) Within the Same Sun-Synchronous Orbit

The NASA CERES EBAF product provides the scientific community observed TOA fluxes to monitor the Earth’s energy imbalance and to validate climate models. To provide a seamless EBAF 24-year record, the CERES instrument calibration must be stable and consistent across satellite records. Once the Terra and Aqua spacecraft are decommissioned at the end of 2025, the CERES project will rely on CERES instrument observations on the SNPP and NOAA-20 1:30 PM sun-synchronous orbits. The SNPP and NOAA-20 satellite orbit placement will prevent any time matched observations for inter-calibration efforts. The future Libera instrument, which will continue the CERES record and scheduled for launch in 2028 onboard the NOAA-22 satellite, will also need to be inter-calibrated without the aid of time-matched observations. Deep convective clouds (DCC) are the most Lambertian, brightest, tropical Earth invariant targets located at the tropopause making them ideal to radiometrically scale the CERES SW observed radiance to a common reference. The SNPP, NOAA-20 and future NOAA-22 satellites are in the same 16-day repeatable orbits, allowing the DCC targets to be observed with the same angular configuration. The empirically derived SNPP CERES SW channel Bidirectional Reflectance Distribution Function (BRDF) can be easily be applied to the NOAA-20 CERES SW channel radiances since they observe the nearly the same DCC systems. The DCC BRDF corrected radiances are analyzed collectively into radiance probability distribution functions (PDF). By comparing the PDF statistics, the SNPP and NOAA-20 CERES SW channel can be inter-calibrated. The inter-calibration coefficients are validated with the CERES instrument calibration team coefficients to optimize the methodology.

David Doelling↗
NASA NTRS1977

Nonlinear Alfven waves in high-speed solar wind streams

A nonlinear proton distribution function that is an exact stationary solution of the nonlinear Vlasov equation and Maxwell's equations and which supports a single nonlinear transverse Alfven (ion cyclotron) wave that is circularly polarized and nondispersive is proposed for most of the observations during high-speed solar wind streams. This nonlinear distribution removes the strong Alfven wave instability, inconsistent with the persistence of the observed proton distribution functions in high-speed streams, found by the linear stability analysis. Model temperature anisotropies and drift velocities of the two spatially inhomogeneous bi-Maxwellian components are consistent with typical proton velocity distributions measured in high-speed streams at 1 AU. Two derived relations for each of the wave number and the phase velocity of the wave are obeyed within experimental uncertainties by two typical proton measurements. Our model also predicts that the alpha particle bulk flow velocity exceeds the proton particle bulk flow velocity, as is observed.

Abraham-Shrauner, B.↗
NASA NTRS1991

Variations in plasma wave intensity with distance along the electron foreshock boundary at Venus

Plasma waves are observed in the solar wind upstream of the Venus bow shock by the Pioneer Venus Orbiter. These wave signatures occur during periods when the interplanetary magnetic field through the spacecraft position intersects the bow shock, thereby placing the spacecraft in the foreshock region. Wave intensity is analyzed as a function of distance along the electron foreshock boundary. It is found that the peak wave intensity may increase along the foreshock boundary from the tangent point to a maximum value at several Venus radii, then decrease in intensity with subsequent increase in distance. These observations could be associated with the instability process: the instability of the distribution function increasing with distance from the tangent point to saturation at the peak. Thermalization of the beam for distances beyond this point could reduce the distribution function instability resulting in weaker wave signatures.

Crawford, G. K.↗
NASA NTRS1990

Neuromorphic learning of continuous-valued mappings from noise-corrupted data. Application to real-time adaptive control

The ability of feed-forward neural network architectures to learn continuous valued mappings in the presence of noise was demonstrated in relation to parameter identification and real-time adaptive control applications. An error function was introduced to help optimize parameter values such as number of training iterations, observation time, sampling rate, and scaling of the control signal. The learning performance depended essentially on the degree of embodiment of the control law in the training data set and on the degree of uniformity of the probability distribution function of the data that are presented to the net during sequence. When a control law was corrupted by noise, the fluctuations of the training data biased the probability distribution function of the training data sequence. Only if the noise contamination is minimized and the degree of embodiment of the control law is maximized, can a neural net develop a good representation of the mapping and be used as a neurocontroller. A multilayer net was trained with back-error-propagation to control a cart-pole system for linear and nonlinear control laws in the presence of data processing noise and measurement noise. The neurocontroller exhibited noise-filtering properties and was found to operate more smoothly than the teacher in the presence of measurement noise.

Troudet, Terry↗
NASA NTRS1972

Resonant diffusion in strongly turbulent plasmas.

The effect of turbulent fluctuations on plasma particles is considered, and equations are derived which describe the evolution of macroscopic properties such as temperature and flow speed of the turbulent plasma. Initially, a diffusion equation for a single-particle distribution function averaged over an ensemble of plasmas is derived for an unmagnetized plasma. For the resonant diffusion in strongly turbulent plasmas, an ensemble of three-dimensional plasmas is considered with an approximately homogeneous and stationary distribution of random electromagnetic fluctuations. For each realization, the single-particle distribution function satisfies the Vlasov equation.

Birmingham, T. J.↗
NASA NTRS1977

The signature of parallel electric fields in a collisionless plasma

The presence of a parallel electric field in the magnetosphere should be apparent through its effect on particles. An analysis is made of particle trajectories in phase space for a steady state collisionless plasma. The analysis utilizes the concepts of distribution function, total energy, and magnetic moment because all three quantities are conserved in a steady state situation in the absence of collisions. Two methods are presented for using particle data to gain information on the potential distribution. The first method involves the use of measurements of both ion and electron distribution functions and invokes the principle of quasi-neutrality. The second method involves the observation of discontinuities in velocity space which can be interpreted in terms of the effect of the electric field on the apparent loss or source cone.

Whipple, E. C., Jr.↗
NASA NTRS1981

Characteristics of reflected and diffuse ions upstream from the earth's bow shock

The distinction between two types of upstream ion populations is made on the basis of pronounced differences in their distribution functions. The reflected ions represent a fast beam with temperatures typically one-million to five-million K and speeds up to five times the solar wind speed. An important feature of the reflected ion distributions is their strong temperature anisotropy, with perpendicular temperature exceeding parallel temperature by a factor of two or three. In contrast, the diffuse ions occupy a much larger region of phase space, both in energy and angle; their distribution function generally has the form of a circular ridge in two dimensions and a spherical shell in three dimensions. Accordingly, their temperature is much larger (not less than about 10-million K) and their bulk speed is typically smaller than the solar wind speed.

Paschmann, G.↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
NASA NTRS1981

Magnetic pumping of particles in the outer Jovian magnetosphere

The mechanism of magnetic pumping consists of two processes, the adiabatic motion of charged particles in a time-varying magnetic field and their pitch angle diffusion. The result is a systematic increase in the energy of charged particles trapped in mirror (and particularly, magnetospheric) magnetic fields. A numerical model of the mechanism is constructed, compared with analytic theory where possible, and, is used to predict the consequences of the process for cases that are not tractable by analytical means. The model is applied to the outer Jovian magnetosphere for two purposes; to find magnetospheric regions in which the mechanism may energize trapped particles, and to generate distribution functions involving pitch angle diffusion caused by wave-particle interactions. Beyond 20 Jupiter radii in the outer magnetosphere particles may be magnetically pumped to energies of the order of 1-2 MeV and two-temperature distribution functions with 'break points' at 1-4 keV for electrons and 8-35 keV for ions are predicted.

Borovsky, J. E.↗
NASA NTRS1986

A second-order accurate kinetic-theory-based method for inviscid compressible flows

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Deshpande, Suresh M.↗
NASA NTRS1996

Dust and Sputtered Particle Streams in Cometary Atmosphere

Dynamics of dust particles in a comet atmosphere was investigated using Monte-Carlo simulation method taking into account radiation and gravitational forces applied to the particles. The model distribution function found out in VEGA experiments was chosen for calculations. Dust particles with masses more then 10(exp -10) g were considered. The delay effect for the big mass particles is shown. In this paper the boundary distribution function for sputtered particles drifting off the surface was determined using Nochilla's model. General dependences for the normal component of the relative mass velocity of sputtered particles are given as the function of particle energy normalized by Tomas-Fermi energy.

Svirschevsky, Stanislav B.↗
DOE OSTI2025

Cosmic Reionization on Computers: Statistical Properties of the Distributions of Mean Opacities

Quasar absorption lines provide a unique window to the relationship between galaxies and the intergalactic medium during the Epoch of Reionization. In particular, high redshift quasars enable measurements of the neutral hydrogen content of the universe. However, the limited sample size of observed quasar spectra, particularly at the highest redshifts, hampers our ability to fully characterize the intergalactic medium during this epoch from observations alone. In this work, we characterize the distributions of mean opacities of the intergalactic medium in simulations from the Cosmic Reionization on Computers (CROC) project. We find that the distribution of mean opacities along sightlines follows a non-trivial distribution that cannot be easily approximated by a known distribution. When comparing the cumulative distribution function of mean opacities measurements in subsamples of sample sizes similar to observational measurements from the literature, we find consistency between CROC and observations at redshifts $z\lesssim 5.7$. However, at higher redshifts ($z\gtrsim5.7$), the cumulative distribution function of mean opacities from CROC is notably narrower than those from observed quasar sightlines implying that observations probe a systematically more opaque intergalactic medium at higher redshifts than the intergalactic medium in CROC boxes at these same redshifts. This is consistent with previous analyses that indicate that the universe is reionized too early in CROC simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Acceleration of low-energy protons and alpha particles at interplanetary shock waves

The low-energy protons and alpha particles in the energy range 30 keV/charge to 150 keV/charge associated with three different interplanetary shock waves in the immediate preshock and postshock region are studied using data obtained by the ISEE 3. The spatial distributions in the preshock and postshock medium are presented, and the dependence of the phase space density at different energies on the distance from the shock and on the form of the distribution function of both species immediately at the shock is examined. It is found that in the preshock region the particles are flowing in the solar wind frame of reference away from the shock and in the postshock medium the distribution is more or less isotropic in this frame of reference. The distribution function in the postshock region can be represented by a power law in energy which has the same spectral exponent for both protons and alpha particles. It is concluded that the first-order Fermi acceleration process can consistently explain the data, although the spectra of diffuse bow shock associated particles are different from the spectra of the interplanetary shock-associated particles in the immediate vicinity of the shock. In addition, the mean free path of the low energy ions in the preshock medium is found to be considerably smaller than the mean free path determined by the turbulence of the background interplanetary medium.

Scholer, M.↗
NASA NTRS2019

The Characterization of Secondary Interstellar Neutral Oxygen Beyond the Heliopause: A Detailed Analysis of the IBEX-Lo Oxygen Observations

In this study, we analyze the directional distribution of the secondary interstellar neutral (ISN) O population observed by the IBEX-Lo neutral atom camera on the Interstellar Boundary EXplorer (IBEX) via the comparison with simulated ISN O intensity maps produced by an analytical model. In the analytical model, we assume that there are primary and secondary ISN populations at the heliopause. We further assume that each population is represented by a Maxwellian velocity distribution function with its own flow parameters. For the viewing directions of IBEX-Lo, we compute the incoming atom speeds at the heliopause with a Keplerian equation of motion in the solar gravity field. Then, we calculate analytically the distribution function to obtain the ISN intensities at Earth’s orbit. We compare the simulated O intensity maps with the IBEX-Lo O sky map to determine the most likely flow parameters of the secondary ISN O population. Using this method, we find the most likely flow parameters of the secondary ISN O population: V(sub SecISNO) = 11 ± 2.2 km s(exp -1), λ(sub SecISNO) = 67° ± 1°.5, β(sub SecISNO) = -12° ± 1°.6, and T(sub SecISNO) = 10,000 ± 1500 K. The results indicate that the secondary ISN O flow direction is deflected toward lower ecliptic longitude and higher negative ecliptic latitude from the ISN gas flow direction at the heliopause. The secondary ISN O flow direction is more deflected from the ISN gas flow direction than the secondary ISN He flow direction.

Interstellar matter clouds↗