Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle distributions”

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

An analysis of whistler waves at interplanetary shocks

We present an analysis of whistler wave magnetic and electric field amplitude ratios from which we compute wave propagation angles and energies of electrons in resonance with the waves. To do this analysis, we compute the theoretical dependence of ratios of wave components on the whistler wave propagation angle Theta for various combinations of orthogonal wave components. Ratios of wave components that would be observed by a spinning spacecraft are determined, and the effects of arbitrary inclinations of the spacecraft to the ambient magnetic field and to the whistler wave vector are studied. This analysis clearly demonstrates that B/E, the ratio of magnetic to electric field amplitudes, cannot be assumed to be the wave index of refraction, contrary to assumptions of some earlier studies. Therefore previous interpretations of whistler wave observations based on this assumption must be reinvestigated. B/E ratios derived using three orthogonal wave components can be used to unambiguously determine Theta. Using spin plane observations alone, a significant uncertainty occurs in the determination of Theta. Nevertheless, for whistler waves observed downstream of several interplanetary shocks by the Ulysses plasma wave experiment we find that Theta is highly oblique. We suggest that the analysis of wave amplitude ratios used in conjunction with traditional stability analyses provide a promising tool for determining which particle distributions and resonances are likely to be dominant contributors to wave growth.

Lengyel-Frey, D.↗

Ulysses field and plasma observations of magnetic holes in the solar wind and their relation to mirror-mode structures

The term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

Winterhalter, Daniel↗

VLF waves in the foreshock

Plasma waves observed in the VLF range upstream of planetary bow shocks not only modify the particle distributions, but also provide important information about the acceleration processes that occur at the bow shock. Electron plasma oscillations observed near the tangent field line in the electron foreshock are generated by electrons reflected at the bow shock through a process that has been referred to as Fast Fermi acceleration. Fast Fermi acceleration is the same as shock-drift acceleration, which is one of the mechanisms by which ions are energized at the shock. We have generated maps of the VLF emissions upstream of the Venus bow shock, using these maps to infer properties of the shock energization processes. We find that the plasma oscillations extend along the field line up to a distance that appears to be controlled by the shock scale size, implying that shock curvature restricsts the flux and energy of reflected electrons. We also find that the ion acoustic waves are observed in the ion foreshock, but at Venus these emissions are not detected near the ULF forshock boundary. Through analogy with terrestrial ion observations, this implies that the ion acoustic waves are not generated by ion beams, but are instead generated by diffuse ion distributions found deep within the ion foreshock. However, since the shock is much smaller at Venus, and there is no magnetosphere, we might expect ion distributions within the ion foreshock to be different than at the Earth. Mapping studies of the terrestrial foreshock similar to those carried out at Venus appear to be necessary to determine if the inferences drawn from Venus data are applicable to other foreshocks.

Strangeway, R. J.↗

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗

Thermal emission from particulate surfaces: A comparison of scattering models with measured spectra

Emissivity spectra of particulate mineral samples are highly dependent on particle size when that size is comparable to the wavelength of light emitted (5-50 micrometers for the midinfrared). Proper geologic interpretation of data from planetary infrared spectrometers will require that these particle size effects be well understood. To address this issue, samples of quartz powders were produced with narrow, well-characterized particle size distributions. Mean particle diameters in these samples ranged from 15 to 227 micrometers. Emission spectra of these powders allow the first detailed comparison of the complex spectral variations with particle size observed in laboratory data with the predictions of radiative transfer models. Four such models are considered here. Hapke's relectance theory (converted to emissivity via Kirchoff's law) is the first model tested. Hapke's more recently published emission theory is also employed. The third model, the 'Mie/Conel' model, uses Mie single scattering with a two-stream approximation for multiple scattering. This model, like the first, is a converted reflec- tance model. Mie scattering assumes particles are both spherical and well separated, which is not true for the quartz powders, but includes diffraction effects. The fourth model uses the Mie solution for single scattering by spheres and inputs those results into the multiple scattering formalism of Hapke's emission theory. The results of the four models are considered in relation to the values of the optical constants n and k. We have grouped these as class 1 (k large), class 2 (k moderate, n is approximately 2), class 3 (k small, n is approximately 2), and class 4 (k small, n is approximately 1). In general, the Mie/Hapke hybrid model does best at predicting variations with grain size. In particular, it predicts changes of the correct pattern, although incorrect magnitude, for class 1 bands, where large increases in emissivity with decreasing grain size are observed. This model also does an excellent job on moderate (class 2) and very weak and intraband (class 3) regions, and correctly predicts the emission maximum and its invariance with grain size near the Christiansen frequency (class 4). The Mie/Hapke hybrid model also has the fewest free parameters of the four models examined, while maintaining the most physical treatment of the radiative transfer. The Mie/Conel model performs as well as the Mie/Hapke hybrid model in strong bands (class 1) but does not accurately model the behavior of moderate (class 2) and very weak (class 3) bands.

Moersch, J. E.↗

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗

Heat flux and viscosity of ions in the collisionless solar wind

Between 1 and 2 solar radii, the Coulomb-collision mean free path for thermal ions exceeds the scale height of the solar atmosphere. The expanding solar plasma becomes collisionless and the kinetics of the solar wind are no longer dominated by thermalizing collisions. The usual Braginskii-type expressions for solar wind ion heat flux and viscosity are no longer valid. However, another microscale still exists in the solar wind, dictated by the gyro-radius of ions in the turbulent embedded solar wind magnetic field. Wave-particle interactions will act to isotropize (but not thermalize) particle distributions, and the relevant microscale for this process is the ion gyro-radius. The ion distribution can be modelled as undergoing isotropizing 'collisions,' with the relevant mean free path scaling with gyro-radius. Here, the author presents the heat flux and viscosity expected for solar wind protons which are relaxing to isotropy on a microscale that scales with gyro-radius. The collisionless viscosity and heat flux have a functional dependence different than their collisional analogs. The collisional expressions for ion viscosity and heat flux drastically overestimate the efficiency of diffusive energy and momentum transport actually operative in the solar wind.

Williams, L. L.↗

Studies of Interstellar Pickup Ions in the Solar Wind

The work under this grant involves studies of the interaction of interstellar pickup ions with the solar wind, with the goal of a comprehensive model of the particle distributions and wave intensities to be expected throughout the heliosphere, as well as the interactions of those distributions with the solar wind termination shock. In the past year, we have completed a number of projects, including observations and modeling of the effects of a large scattering mean free path on the pickup He(+) seen at AMPTE, an analytical model of anisotropic pickup tons in a steady radial magnetic field, and a derivation of a reduced solar wind Mach number due to increased estimates on the inflowing hydrogen density allowing for a weak termination shock. In the next year, we plan to investigate in more detail the correspondence between our models of anisotropic pickup ions and the data on spectra, variations, and proton-He(+) correlation provided by AMPTE, Ulysses, and our instrument on SOHO. We will model the time-dependent pickup ion density resulting from finite periods of radial magnetic field. We will also incorporate the effects of a large mean free path into our analysis of the He(+) focusing cone, leading to more accurate parameter values for the interstellar helium gas. This progress report also includes a discussion of our Space Physics Educational Outreach activities in the past year and plans for the next year.

Isenberg, Philip A.↗

A Gas-Kinetic Scheme for Reactive Flows

In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

Lian,Youg-Sheng↗

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗

The Development of an Edible Peanut Protein Film

The peanut is one of the crops chosen for use in NASA's Advanced Life Support Program (ALS). The peanut is a source of both oil and protein. After oil is extracted from the peanut, a protein rich flour remains. An edible peanut protein film is one use for this flour. Two types of film are developed for this study, one set of film contains 10% fat while the other set contains no fat. For film without fat the defatting of the peanut by the Soxhlet method is the first step in the manufacturing process of the film. Secondly, the protein is precipitated at its isoelectric point (pH 4.5) and centrifuged to separate the protein from the non-protein. After freeze-drying the protein it is milled in a ceramic ball mill to decrease particle size and sifted through a series of sieves to determine particle size distribution. Those particles retained on the 100 mesh sieves are utilized for film formation. Larger particles are re-ground and sifted. Five grams of protein is mixed with 50 mL of distilled water, 70 mL of 80% ethanol, 15 mL of 6N ammonium hydroxide and a plasticizer. This mixture is heated for 30 minutes until the temperature reaches 70 C. The mixture is then poured onto a level Teflon coated glass surface. After allowing the film to form overnight under a ventilation hood, it is manually removed from the plate. The processes and methods adopted have created flexible films of uniform thickness that are free of air bubbles. Thickness of films made from defatted peanut protein and partially defatted peanut protein were 0.10 Lm and 0.13 Lm respectively. Films with natural peanut fat are approximately three times as flexible and almost four times as strong as the films made without fat. Further research will be performed to evaluate its mechanical properties. This paper will greatly contribute to food preservation and waste management. Potential applications of this film are edible/biodegradable containers, wrapping for food preservation (against water, oxygen and oil), storage packets for seasonings or other ingredients and encapsulation for pharmaceutical use.

Patrick, N.↗

Laser Transmission Measurements of Soot Extinction Coefficients in the Exhaust Plume of the X-34 60K-lb Thrust Fastrac Rocket Engine

A four-channel laser transmissometer has been used to probe the soot content of the exhaust plume of the X-34 60k-lb thrust Fastrac rocket engine at NASA's Marshall Space Flight Center. The transmission measurements were made at an axial location approximately equal 1.65 nozzle diameters from the exit plane and are interpreted in terms of homogeneous radial zones to yield extinction coefficients from 0.5-8.4 per meter. The corresponding soot mass density, spatially averaged over the plume cross section, is, for Rayleigh particles, approximately equal 0.7 microgram/cc, and alternative particle distributions are briefly considered. Absolute plume radiance at the laser wavelength (515 nm) is estimated from the data at approximately equal 2,200 K equivalent blackbody temperature, and temporal correlations in emission from several spatial locations are noted.

Dobson, C. C.↗

Laser Transmission Measurements of Soot Extinction Coefficients in the Exhaust Plume of the X-34 60k-lb Thrust Fastrac Rocket Engine

A four-channel laser transmissometer has been used to probe the soot content of the exhaust plume of the X-34 60k-lb thrust Fastrac rocket engine at NASA's Marshall Space Flight Center. The transmission measurements were made at an axial location about equal 1.65 nozzle diameters from the exit plane and are interpreted in terms of homogeneous radial zones to yield extinction coefficients from 0.5-8.4 per meter. The corresponding soot mass density, spatially averaged over the plume cross section, is, for Rayleigh particles, approximately equal to 0.7 micrograms/cubic cm and alternative particle distributions are briefly considered. Absolute plume radiance at the laser wavelength (515 nm) is estimated from the data at approximately equal to 2.200 K equivalent blackbody temperature, and temporal correlations in emission from several spatial locations are noted.

Dobson, C. C.↗

Photochemical Phenomenology Model for the New Millennium

The "Photochemical Phenomenology Model for the New Millennium" project tackles the issue of reengineering and extension of validated physics-based modeling capabilities ("legacy" computer codes) to application-oriented software for use in science and science-support activities. While the design and architecture layouts are in terms of general particle distributions involved in scattering, impact, and reactive interactions, initial Photochemical Phenomenology Modeling Tool (PPMT) implementations are aimed at construction and evaluation of photochemical transport models with rapid execution for use in remote sensing data analysis activities in distributed systems. Current focus is on the Composite Infrared Spectrometer (CIRS) data acquired during the CASSINI flyby of Jupiter. Overall, the project has stayed on the development track outlined in the Year 1 annual report and most Year 2 goals have been met. The issues that have required the most attention are: implementation of the core photochemistry algorithms; implementation of a functional Java Graphical User Interface; completion of a functional CORBA Component Model framework; and assessment of performance issues. Specific accomplishments and the difficulties encountered are summarized in this report. Work to be carried out in the next year center on: completion of testing of the initial operational implementation; its application to analysis of the CASSINI/CIRS Jovian flyby data; extension of the PPMT to incorporate additional phenomenology algorithms; and delivery of a mature operational implementation.

Bishop, James↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

Force Evaluation in the Lattice Boltzmann Method Involving Curved Geometry

The present work investigates two approaches for force evaluation in the lattice Boltzmann equation: the momentum- exchange method and the stress-integration method on the surface of a body. The boundary condition for the particle distribution functions on curved geometries is handled with second order accuracy based on our recent works. The stress-integration method is computationally laborious for two-dimensional flows and in general difficult to implement for three-dimensional flows, while the momentum-exchange method is reliable, accurate, and easy to implement for both two-dimensional and three-dimensional flows. Several test cases are selected to evaluate the present methods, including: (i) two-dimensional pressure-driven channel flow; (ii) two-dimensional uniform flow past a column of cylinders; (iii) two-dimensional flow past a cylinder asymmetrically placed in a channel (with vortex shedding); (iv) three-dimensional pressure-driven flow in a circular pipe; and (v) three-dimensional flow past a sphere. The drag evaluated by using the momentum-exchange method agrees well with the exact or other published results.

Mei, Renwei↗

Characterization of Magnetospheric Spacecraft Charging Environments Using the LANL Magnetospheric Plasma Analyzer Data Set

An improved specification of the plasma environment has been developed for use in modeling spacecraft charging. It was developed by statistically analyzing a large part of the LANL Magnetospheric Plasma Analyzer (MPA) data set for ion and electron spectral signature correlation with spacecraft charging, including anisotropies. The objective is to identify a relatively simple characterization of the full particle distributions that yield an accurate predication of the observed charging under a wide variety of conditions.

Hardage, Donna↗

Quenching Combustible Dust Mixtures Using Electric Particulate Suspensions (EPS): A New Testing Method For Microgravity

The Electric Particulate Suspension (EPS) is a combustion ignition system being developed at Iowa State University for evaluating quenching effects of powders in microgravity (quenching distance, ignition energy, flammability limits). Because of the high cloud uniformity possible and its simplicity, the EPS method has potential for "benchmark" design of quenching flames that would provide NASA and the scientific community with a new fire standard. Microgravity is expected to increase suspension uniformity even further and extend combustion testing to higher concentrations (rich fuel limit) than is possible at normal gravity. Two new combustion parameters are being investigated with this new method: (1) the particle velocity distribution and (2) particle-oxidant slip velocity. Both walls and (inert) particles can be tested as quenching media. The EPS method supports combustion modeling by providing accurate measurement of flame-quenching distance as a parameter in laminar flame theory as it closely relates to characteristic flame thickness and flame structure. Because of its design simplicity, EPS is suitable for testing on the International Space Station (ISS). Laser scans showing stratification effects at 1-g have been studied for different materials, aluminum, glass, and copper. PTV/PIV and a leak hole sampling rig give particle velocity distribution with particle slip velocity evaluated using LDA. Sample quenching and ignition energy curves are given for aluminum powder. Testing is planned for the KC-135 and NASA s two second drop tower. Only 1-g ground-based data have been reported to date.

Colver, Gerald M.↗