Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Significant initial results from the environmental measurements experiment on ATS-6

The Applications Technology Satellite (ATS-6), launched into synchronous orbit on 30 May 1974, carried a set of six particle detectors and a triaxial fluxgate magnetometer. The particle detectors were able to determine the ion and electron distribution functions from 1 to greater than 10 to the 8th power eV. It was found that the magnetic field is weaker and more tilted than predicted by models which neglect internal plasma and that there is a seasonal dependence to the magnitude and tilt. ATS-6 magnetic field measurements showed the effects of field-aligned currents associated with substorms, and large fluxes of field-aligned particles were observed with the particle detectors. Encounters with the plasmasphere revealed the existence of warm plasma with temperatures up to 30 eV. A variety of correlated waves in both the particles and fields were observed: pulsation continuous oscillations, seen predominantly in the plasmasphere bulge; ultralow frequency (ULF) standing waves; ring current proton ULF waves; and low frequency waves that modulate the energetic electrons. In additon, large scale waves on the energetic-ion-trapping boundary were observed, and the intensity of energetic electrons was modulated in association with the passage of sector boundaries of the interplanetary magnetic field.

Fritz, T. A.↗
NASA NTRS1979

Accuracy of RGD approximation for computing light scattering properties of diffusing and motile bacteria

The study tests the accuracy of the Rayleigh-Gans-Debye (RGD) approximation against a rigorous scattering theory calculation for a simplified model of E. coli (about 1 micron in size) - a solid spheroid. A general procedure is formulated whereby the scattered field amplitude correlation function, for both polarized and depolarized contributions, can be computed for a collection of particles. An explicit formula is presented for the scattered intensity, both polarized and depolarized, for a collection of randomly diffusing or moving particles. Two specific cases for the intermediate scattering functions are considered: diffusing particles and freely moving particles with a Maxwellian speed distribution. The formalism is applied to microorganisms suspended in a liquid medium. Sensitivity studies revealed that for values of the relative index of refraction greater than 1.03, RGD could be in serious error in computing the intensity as well as correlation functions.

Kottarchyk, M.↗
DOE OSTI2021

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. In this work, we present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, M ν ≲ 0.6eV, and applicable to redshifts z ≲ 2. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in ~ 25 steps, with a starting redshift of z=99. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between ~ 5-20% for runs with 1024 3 CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction.

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 NTRS2001

Orbits of Saturn's Inner Moons and Other Observations Connected within the 1995-1996 Saturnian Ring Plane Crossing

We reduced and analyzed radio images of Saturn taken in 1990-1995 at 0.35, 2.0, 3.6 and 6.1 cm wavelength. The latitudinal brightness distribution varies substantially over time. Detailed modeling of these data shows that NH4SH condenses at the thermochemical equilibrium temperature of 235.5 K only at Northern mid-latitudes. Over most of the planet, condensation does not occur until the gas cools to 190 +/- 5 K. Supersaturation may also cause the dark equatorial region seen in 1995 at 6.1 cm. Observations of the rings show that the West (dusk) ansa is brighter than the East (dawn) ansa by factors up to 2. The magnitude of the asymmetry increases with increasing wavelength and with decreasing distance to the planet. This East-West asymmetry may be due to multiple scattering in gravitational wakes. We improved the Monte Carlo scattering code used to model the microwave emission from Saturn's rings at multiple wavelengths. The basic code outlined in Dunn, Molnar, and Fix describes the rings as a uniform, infinite layer of ice spheres with a finite vertical extent. The particle sizes follow a power law distribution and scatter light according to a linear combination of Mie and isotropic phase functions. These parameters may be adjusted to reflect the physical properties in the scattering properties of the ring particles. This is used to construct a model which can be directly compared to Very Large Array data. The code has been extended by including anisotropic particle distributions, primarily in the form of the wake structures that are likely to exist in the enhancements. The vertical extent of the wakes equals that of the ring layer. Between the wakes are regions of diminished density so that the optical depth integrated over many wakes matches that of the Voyager occultation experiments. This allows us to model the brightness temperature as a function of azimuth and look for asymmetric scattering behavior in the rings. By changing the width and relative density of the wakes, we can alter the nature of the gross scattering properties of the rings to match what is seen in our data.

Lissauer, Jack J.↗
NASA NTRS1976

Particles from comet Kohoutek detected by the micrometeoroid experiment on HEOS 2

HEOS B measurements on particles ejected from comet Kohoutek reflect average particle rate as a function of particle speed and mass in relation to random distribution with known speed from the interplanetary region. The micrometeoroid experiment detector onboard the satellite passed through the orbital plane of the comet and encountered ejected particles for approximately two months.

Hoffmann, H. J.↗
NASA NTRS1980

A simulation study of cold electron heating by loss cone instabilities

Results are presented of a computer simulation study of electron cyclotron harmonic waves; an electrostatic finite size particle simulation code is used. The initial electron velocity distribution function is modeled as a hot ring perpendicular to the ambient magnetic field and a cold Maxwellian component. Nonlinear cyclotron resonance is found to be the saturation mechanism of the unstable waves and the heating mechanism of the cold electrons. The results can be related to certain ionospheric phenomena.

Ashour-Abdalla, M.↗
NASA NTRS1989

Particle distributions in a two-dimensional reconnection field geometry

In order to understand space phase density variations near the magnetopause and in the geomagnetic tail, particle trajectories in a two-dimensional steady state reconnection field geometry are studied by tracing particles from a spacecraft position back in time. The distribution function at the spacecraft is calculated using Liouville's theorem and known distribution functions in the source regions. It is found, in particular, that multiple ion beams in the plasma sheet boundary layer are a natural consequence of reconnection. In the center of the central plasma sheet, the omnidirectional flux of keV ions should be at a minimum. The depth of this minimum is very large near the neutral point but decreases as the distance from the neutral point increases.

Gurran, D. B.↗
NASA NTRS2011

Uncertainty Analysis in Space Radiation Protection

Space radiation is comprised of high energy and charge (HZE) nuclei, protons, and secondary radiation including neutrons. The uncertainties in estimating the health risks from galactic cosmic rays (GCR) are a major limitation to the length of space missions, the evaluation of potential risk mitigation approaches, and application of the As Low As Reasonably Achievable (ALARA) principle. For long duration space missio ns, risks may approach radiation exposure limits, therefore the uncertainties in risk projections become a major safety concern and methodologies used for ground-based works are not deemed to be sufficient. NASA limits astronaut exposures to a 3% risk of exposure induced death (REID) and protects against uncertainties in risks projections using an assessment of 95% confidence intervals in the projection model. We discuss NASA s approach to space radiation uncertainty assessments and applications for the International Space Station (ISS) program and design studies of future missions to Mars and other destinations. Several features of NASA s approach will be discussed. Radiation quality descriptions are based on the properties of radiation tracks rather than LET with probability distribution functions (PDF) for uncertainties derived from radiobiology experiments at particle accelerators. The application of age and gender specific models for individual astronauts is described. Because more than 90% of astronauts are never-smokers, an alternative risk calculation for never-smokers is used and will be compared to estimates for an average U.S. population. Because of the high energies of the GCR limits the benefits of shielding and the limited role expected for pharmaceutical countermeasures, uncertainty reduction continues to be the optimal approach to improve radiation safety for space missions.

Cucinotta, Francis A.↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS2018

Electron Bulk Acceleration and Thermalization at Earth's Quasiperpendicular Bow Shock

Electron heating at Earth's quasiperpendicular bow shock has been surmised to be due to the combined effects of a quasistatic electric potential and scattering through wave-particle interaction. Here we report the observation of electron distribution functions indicating a new electron heating process occurring at the leading edge of the shock front. Incident solar wind electrons are accelerated parallel to the magnetic field toward downstream, reaching an electron-ion relative drift speed exceeding the electron thermal speed. The bulk acceleration is associated with an electric field pulse embedded in a whistler-mode wave. The high electron-ion relative drift is relaxed primarily through a nonlinear current-driven instability. The relaxed distributions contain a beam traveling toward the shock as a remnant of the accelerated electrons. Similar distribution functions prevail throughout the shock transition layer, suggesting that the observedacceleration and thermalization is essential to the cross-shock electron heating.

Chen, L.-J.↗
NASA NTRS2022

Assessing the Sampleability of Bennu’s Surface for the OSIRIS-REx Asteroid Sample Return Mission

NASA’s first asteroid sample return mission, OSIRIS-REx, collected a sample from the surface of near-Earth asteroid Bennu in October 2020 and will deliver it to Earth in September 2023. Selecting a sample collection site on Bennu’s surface was challenging due to the surprising lack of large ponded deposits of regolith particles exclusively fine enough (≤ 2 cm diameter) to be ingested by the spacecraft’s Touch-and-Go Sample Acquisition Mechanism (TAGSAM). Here we describe the Sampleability Map of Bennu, which was constructed to aid in the selection of candidate sampling sites and to estimate the probability of collecting sufficient sample. “Sampleability” is a numeric score that expresses the compatibility of a given area’s surface properties with the sampling mechanism. The algorithm that determines sampleability is a best fit functional form to an extensive suite of laboratory testing outcomes tracking the TAGSAM performance as a function of four observable properties of the target asteroid. The algorithm and testing were designed to measure and subsequently predict TAGSAM collection amounts as a function of the minimum particle size, maximum particle size, particle size frequency distribution, and the tilt of the TAGSAM head off the surface. The sampleability algorithm operated at two general scales, consistent with the resolution and coverage of data collected during the mission. The first scale was global and evaluated nearly the full surface. Due to Bennu’s unexpected boulder coverage and lack of ponded regolith deposits, the global sampleability efforts relied heavily on additional strategies to find and characterize regions of interest based on quantifying and avoiding areas heavily covered by material too large to be collected. The second scale was site-specific and used higher-resolution data to predict collected mass at a given contact location. The rigorous sampleability assessments gave the mission confidence to select the best possible sample collection site and directly enabled successful collection of hundreds of grams of material.

Asteroid exploration↗
NASA NTRS1988

A predicted energetic ion signature of a neutral line in the geomagnetic tail

Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.

Martin, R. F., Jr.↗
NASA NTRS1981

Analytical description of charged particles transport along arbitrary guiding-field configurations

A new description of focused transport was developed which is valid for arbitrary spatial dependences of both the scattering mean free path and the focusing length. In particular, it describes the supercoherent transition from the coherent modes that occurs when focusing is strong compared to scattering to the diffusive mode that occurs when focusing is weak. Because it specifies the pitch angle distribution in terms of relatively simple functions, the new formulation facilitates not only the detailed interpretation of solar particle events, but also the qualitative understanding of focused transport. In the steady state, it predicts that the angular distribution consists of a highly collimated 'strahl' which is superimposed on an anisotropic 'halo'. This picture provides new insight into observed configurations of steady-state particle fluxes.

Earl, J. A.↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

A Radiation Chemistry Code Based on the Greens Functions of the Diffusion Equation

Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.

Plante, Ianik↗