Search NASA⌕ Search

SEARCH · Search NASA

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Update on the Fire (solar probe) mission study

Since mid-1994 the U.S. and Russia have been studying the technical feasibility of a joint solar probe mission as part of the 'Fire and Ice' concept to explore close to the Sun, and Pluto, together. In the current concept of the 'Fire' mission, separate spacecraft built by each country would be launched together, fly by Jupiter to shed orbital angular momentum and achieve a solar polar orbit, and arrive 3.6 years later at 4 and 10 R(sub s). The Fire mission would measure basic parameters of the modes of energy and momentum flow and transfer to the coronal plasma that are not observable remotely. Specifically, measurement of magnetic fields, waves, suprathermal particles, and critical features of the plasma particle composition and distribution function would be made from 4 to 30 R(sub s) where the solar wind is known to be accelerated. In addition, the Fire spacecraft should image coronal structures unambiguously and relate the underlying and flown-through structures to plasma characteristics measured in situ. Each country is developing a backup plan to pursue the solar probe objectives alone if the other side is unable to carry out its mission.

Jones, W. Veron↗
NASA NTRS2005

Microwave Signatures of Melting/Refreezing Snow: Observations and Modeling Using Dense Medium Radiative Transfer Theory

Microwave brightness temperatures of snow covered terrains can be modeled by means of the Dense Radiative Transfer Medium Theory (DMRT). In a dense medium, such as snow, the assumption of independent scattering is no longer valid and the scattering of correlated scatterers must be considered. In the DMRT, this is done considering a pair distribution function of the particles position. In the electromagnetic model, the snowpack is simulated as a homogeneous layer having effective permittivity and albedo calculated through the DMRT. In order to account for clustering of snow crystals, a model of cohesive particles can be applied, where the cohesion between the particles is described by means of a dimensionless parameters called stickiness (z), representing a measure of the inversion of the attraction of the particles. The lower the z the higher the stickiness. In this study, microwave signatures of melting and refreezing cycles of seasonal snowpacks at high altitudes are studied by means of both experimental and modeling tools. Radiometric data were collected 24 hours per day by the University of Michigan Tower Mounted Radiometer System (TMRS). The brightness temperatures collected by means of the TMRS are simulated by means of a multi-layer electromagnetic model based on the dense medium theory with the inputs to the model derived from the data collected at the snow pits and from the meteorological station. The paper is structured as follows: in the first Section the temperature profiles recorded by the meteorological station and the snow pit data are presented and analyzed; in the second Section, the characteristics of the radiometric system used to collect the brightness temperatures are reported together with the temporal behavior of the recorded brightness temperatures; in the successive Section the multi-layer DMRT-based electromagnetic model is described; in the fourth Section the comparison between modeled and measured brightness temperatures is discussed. We dedicate the last Section to the conclusions and future works.

Tedesco, Marco↗
NASA NTRS2007

Space Radiation Risk Assessment for Future Lunar Missions

For lunar exploration mission design, radiation risk assessments require the understanding of future space radiation environments in support of resource management decisions, operational planning, and a go/no-go decision. The future GCR flux was estimated as a function of interplanetary deceleration potential, which was coupled with the estimated neutron monitor rate from the Climax monitor using a statistical model. A probability distribution function for solar particle event (SPE) occurrence was formed from proton fluence measurements of SPEs occurred during the past 5 solar cycles (19-23). Large proton SPEs identified from impulsive nitrate enhancements in polar ice for which the fluences are greater than 2 10(exp 9) protons/sq cm for energies greater than 30 MeV, were also combined to extend the probability calculation for high level of proton fluences. The probability with which any given proton fluence level of a SPE will be exceeded during a space mission of defined duration was then calculated. Analytic energy spectra of SPEs at different ranks of the integral fluences were constructed over broad energy ranges extending out to GeV, and representative exposure levels were analyzed at those fluences. For the development of an integrated strategy for radiation protection on lunar exploration missions, effective doses at various points inside a spacecraft were calculated with detailed geometry models representing proposed transfer vehicle and habitat concepts. Preliminary radiation risk assessments from SPE and GCR were compared for various configuration concepts of radiation shelter in exploratory-class spacecrafts.

Kim, Myung-Hee Y.↗
NASA NTRS2008

Sensor Technologies for Particulate Detection and Characterization

Planned Lunar missions have resulted in renewed attention to problems attributable to fine particulates. While the difficulties experienced during the sequence of Apollo missions did not prove critical in all cases, the comparatively long duration of impending missions may present a different situation. This situation creates the need for a spectrum of particulate sensing technologies. From a fundamental perspective, an improved understanding of the properties of the dust fraction is required. Described here is laboratory-based reference instrumentation for the measurement of fundamental particle size distribution (PSD) functions from 2.5 nanometers to 20 micrometers. Concomitant efforts for separating samples into fractional size bins are also presented. A requirement also exists for developing mission compatible sensors. Examples include provisions for air quality monitoring in spacecraft and remote habitation modules. Required sensor attributes such as low mass, volume, and power consumption, autonomy of operation, and extended reliability cannot be accommodated by existing technologies.

Greenberg, Paul S.↗
NASA NTRS2021

Polarimetric sensitivity of light-absorbing carbonaceous aerosols over ocean: A theoretical assessment

Visible-light-absorbing carbonaceous aerosols within the boundary layer affect the radiance and polarization states of the radiation at the top of the atmosphere. Remote sensing from suborbital and satellite-based platforms utilizes these radiance and polarization signals to retrieve the key properties of these aerosols. Recent retrieval algorithms have shown a progressive trend toward including multi-angular and multi-spectral polarimetric measurements to produce better retrieval accuracy in comparison to those using measurements based on a single viewing angle. Here, we perform a theoretical investigation of the top of atmosphere (TOA) radiance-related reflectance factor (bidirectional reflectance factor (BRF)) and the two types of polarimetry-related factors (polarized bidirectional reflectance factor (pBRF) and the degree of linear polarization (DoLP)) for different types of atmospheric light-absorbing carbonaceous aerosols as a function of particle size distribution. We selected three polarimetric bands corresponding to those utilized by NASA's Airborne Multiangle SpectroPolarimetric Imager (AirMSPI)—near-UV (470 nm), visible (660 nm), and near-infrared (865 nm)—for our simulations which were performed over ocean surface using the successive order of scattering (SOS) algorithm coupled to a Lorenz-Mie aerosol optics model. The analysis of particle phase matrix elements indicates a close relationship between the angular dependencies of DoLP and associated phase matrix components at the shortest polarimetric band (470 nm). Using Jacobian analysis, we find that the radiance- and polarimetry-related reflectance factors of weakly light-absorbing aerosols, such as brown carbon, are more sensitive to changes in particle size and imaginary refractive index in comparison with those of black carbon, which is strongly light-absorbing. Our results suggest that the DoLP data could be used by future retrieval algorithms for reliably estimating microphysical properties of absorbing carbonaceous aerosols with imaginary refractive index less than 0.4.

refractive index↗
NASA NTRS2026

Recommended Practices in Laser-Induced Fluorescence (LIF) Diagnostics for Electric Propulsion

Laser-induced fluorescence (LIF) spectroscopy uses Doppler-shifted laser photons to excite a bound electron transition in an atom or ion, with subsequent fluorescence emission detected in order to measure the local velocity distribution function (VDF) and/or particle density. Due to its non-invasive nature, high spatial resolution, and acceptable difficulty of implementation, LIF has become a favored diagnostic for electric propulsion (EP) plasmas, particularly in the study of Hall thrusters. This paper presents recommended best practices for LIF measurements in electric thrusters, drawn from over three decades of implementation heritage in the EP community. Focusing on both single-photon LIF (typically used to measure VDFs) and two-photon LIF (TALIF, typically used to determine densities), the paper covers selection of atomic transitions to target, experimental setups, and interpretation of data. A number of special topics and advanced applications are discussed, including uncertainty analysis, techniques for high-speed measurements, and applications to a variety of propellants.

plasma plume↗
NASA NTRS1981

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2022

A Generalized Dual-Frequency Ratio (DFR) Approach for Rain Retrievals

The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.

Robert Meneghini↗
DOE OSTI2023

Electrostatic upper-hybrid mode instability driven by a ring electron distribution

Quasi electrostatic fluctuations in the upper-hybrid frequency range are commonly detected in the planetary magnetospheric environment. The origin of such phenomena may relate to the instability driven by a loss-cone feature associated with the electrons populating the dipole-like magnetic field. The present paper carries out a one-dimensional electrostatic particle-in-cell simulation accompanied by a reduced quasilinear kinetic theoretical analysis to investigate the dynamics of the upper-hybrid mode instability driven by an initial ring electron distribution function, which is a form of loss-cone distribution. A favorable comparison is found between the two approaches, which shows that the reduced quasilinear theory, which is grounded in the concept of a model of the particle distribution function that is assumed to maintain a fixed mathematical form except that the macroscopic parameters that define the distribution are allowed to evolve in time, can be an effective tool in the study of plasma instabilities, especially if it is guided by and validated against the more rigorous simulation result.

Physics↗
NASA NTRS1981

Observations of the ionization states of energetic particles accelerated in solar flares

The ionization states and spectra of 0.3 to 2.4 MeV/nuc He, C, O and Fe are measured in a survey of ten solar flare particle events. He(plus) is found to be present in all events, which indicates the common presence of energetic He(plus) in the source material from which solar particles are accelerated. The distribution functions of He(plus), He(plus plus) and heavier elements are represented by simple exponentials of the particle speed times its rigidity to a power n, where n is between 0 and 1, and equal e-folding values. Results are consistent with a model where ions are accelerated in the corona by multi-dimensional shocks out of a population taken from both hot and cold coronal regions.

Gloeckler, G.↗
NASA NTRS2019

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

A model of current collection in Langmuir probes is used to investigate the effect representative non- equilibrium plasmas under various conditions have on the electron temperature and number density that would be calculated through analysis of the probe collection characteristics. The model uses the distribution function to calculate the charged particle flux to a probe and then, for fixed applied voltages, the current values that satisfy continuity in the probes are determined. The triple probe is not scanned in voltage, so there is no practical way using experimental triple probe data to determine if the plasma is in equilibrium. As a consequence, a triple probe analysis typically relies on the assumption that the plasma is in equilibrium. Proceeding from this point, the numerically-generated non-equilibrium triple probe data are analyzed assuming that the plasma is in equilibrium, with the data compared to the initial distribution function inputs of plasma temperature and number density to determine the effect the non-equilibrium distribution has on plasma measurements. The temperature and number density are both significantly affected when a fraction of the particles in the distribution are shifted from the equilibrium configuration into the non-equilibrium part of the distribution function. For all instances studied, the computed electron temperature and number density are extremely sensitive to small deviations from equilibrium ( ≤ 5% of the plasma shifted into the non-equilibrium function). Shifting more of the plasma into the non-equilibrium distribution beyond this initial level does not produce a significant additional shift in the computed plasma properties.

Polzin, Kurt A.↗
NASA NTRS2001

The Space Debris Environment for the ISS Orbit

With thirty-five planned missions over the next five years, the International Space Station (ISS) will be the focus for manned space activity. At least 6 different vehicles will transport crew and supplies to and from the nominally 400 km, 51.6 degree orbit. When completed, the ISS will be the largest space structure ever assembled and hence the largest target for space debris. Recent work at the Johnson Space Center has focused on updating the existing space debris models. The Orbital Debris Engineering Model, has been restructured to take advantage of state of the art desktop computing capability and revised with recent measurements from Haystack and Goldstone radars, additional analysis of LDEF and STS impacts, and the most recent SSN catalog. The new model also contains the capability to extrapolate the current environment in time to the year 2030. A revised meteoroid model based on the work of Divine has also been developed, and is called the JSC Meteoroid Model. The new model defines flux on the target per unit angle per unit speed, and for Earth orbit, includes the meteor showers. This paper quantifies the space debris environment for the ISS orbit from natural and anthropogenic sources. Particle flux and velocity distributions as functions of size and angle are be given for particles 10 microns and larger for altitudes from 350 to 450 km. The environment is projected forward in time until 2030.

Theall, Jeff↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗