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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Scale Dependence of Multiplier Distributions for Particle Concentration, Enstrophy, and Dissipation in the Inertial Range of Homogeneous Turbulence

Turbulent flows preferentially concentrate inertial particles depending on their stopping time or Stokes number, which can lead to significant spatial variations in the particle concentration. Cascade models are one way to describe this process in statistical terms. Here, we use a direct numerical simulation (DNS) dataset of homogeneous, isotropic turbulence to determine probability distribution functions (PDFs) for cascade multipliers, which determine the ratio by which a property is partitioned into subvolumes as an eddy is envisioned to decay into smaller eddies. We present a technique for correcting effects of small particle numbers in the statistics. We determine multiplier PDFs for particle number, flow dissipation, and enstrophy, all of which are shown to be scale dependent. However, the particle multiplier PDFs collapse when scaled with an appropriately defined local Stokes number. As anticipated from earlier works, dissipation and enstrophy multiplier PDFs reach an asymptote for sufficiently small spatial scales. From the DNS measurements, we derive a cascade model that is used it to make predictions for the radial distribution function (RDF) for arbitrarily high Reynolds numbers, Re, finding good agreement with the asymptotic, infinite Re inertial range theory of Zaichik and Alipchenkov [New J. Phys. 11, 103018 (2009)]. We discuss implications of these results for the statistical modeling of the turbulent clustering process in the inertial range for high Reynolds numbers inaccessible to numerical simulations.

turbulent flows↗
NASA NTRS1995

The energy associated with MHD waves generation in the solar wind plasma

Gyrotropic symmetry is usually assumed in measurements of electron distribution functions in the heliosphere. This prevents the calculation of a net current perpendicular to the magnetic field lines. Previous theoretical results derived by one of the authors for a collisionless plasma with isotropic electrons in a strong magnetic field have shown that the excitation of MHD modes becomes possible when the external perpendicular current is non-zero. We consider then that any anisotropic electron population can be thought of as 'external', interacting with the remaining plasma through the self-consistent electromagnetic field. From this point of view any perpendicular current may be due to the anisotropic electrons, or to an external source like a stream, or to both. As perpendicular currents cannot be derived from the measured distribution functions, we resort to Ampere's law and experimental data of magnetic field fluctuations. The transfer of energy between MHD modes and external currents is then discussed.

delaTorre, A.↗
NASA NTRS2014

Wide-Field Hubble Space Telescope Observations of the Globular Cluster System in NGC 1399*

We present a comprehensive high spatial resolution imaging study of globular clusters (GCs) in NGC 1399, thecentral giant elliptical cD galaxy in the Fornax galaxy cluster, conducted with the Advanced Camera for Surveys(ACS) aboard theHubble Space Telescope(HST).Using a novel technique to construct drizzled point-spreadfunction libraries for HSTACS data, we accurately determine the fidelity of GC structural parameter measurementsfrom detailed artificial star cluster experiments and show the superior robustness of the GC half-light radius,rh,compared with other GC structural parameters, such as King core and tidal radius. The measurement ofrhfor themajor fraction of the NGC 1399 GC system reveals a trend of increasingrhversus galactocentric distance,Rgal,out to about 10 kpc and a flat relation beyond. This trend is very similar for blue and red GCs, which are found tohave a mean size ratio ofrh,redrh,blue0.820.11 at all galactocentric radii from the core regions of the galaxyout to40 kpc. This suggests that the size difference between blue and red GCs is due to internal mechanismsrelated to the evolution of their constituent stellar populations. Modeling the mass density profile of NGC 1399shows that additional external dynamical mechanisms are required to limit the GC size in the galaxy halo regionstorh2 pc. We suggest that this may be realized by an exotic GC orbit distribution function, an extended darkmatter halo, andor tidal stress induced by the increased stochasticity in the dwarf halo substructure at largergalactocentric distances. We compare our results with the GCrhdistribution functions in various galaxies and findthat the fraction of extended GCs withrh5 pc is systematically larger in late-type galaxies compared with GCsystems in early-type galaxies. This is likely due to the dynamically more violent evolution of early-type galaxies.We match our GCrhmeasurements with radial velocity data from the literature and split the resulting sample at themedianrhvalue into compact and extended GCs. We find that compact GCs show a significantly smaller line-of-sight velocity dispersion,cmp22525 km s1, than their extended counterparts,ext31721 km s1.Considering the weaker statistical correlation in the GCrhcolor and the GCrhRgalrelations, the more significantGC sizedynamics relation appears to be astrophysically more relevant and hints at the dominant influence of theGC orbit distribution function on the evolution of GC structural parameters.

Hubble↗
NASA NTRS1994

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗
NASA NTRS1991

A review of the kinetic detail required for accurate predictions of normal shock waves

Several aspects of the kinetic models used in the collision phase of Monte Carlo direct simulations have been studied. Accurate molecular velocity distribution function predictions require a significantly increased number of computational cells in one maximum slope shock thickness, compared to predictions of macroscopic properties. The shape of the highly repulsive portion of the interatomic potential for argon is not well modeled by conventional interatomic potentials; this portion of the potential controls high Mach number shock thickness predictions, indicating that the specification of the energetic repulsive portion of interatomic or intermolecular potentials must be chosen with care for correct modeling of nonequilibrium flows at high temperatures. It has been shown for inverse power potentials that the assumption of variable hard sphere scattering provides accurate predictions of the macroscopic properties in shock waves, by comparison with simulations in which differential scattering is employed in the collision phase. On the other hand, velocity distribution functions are not well predicted by the variable hard sphere scattering model for softer potentials at higher Mach numbers.

Muntz, E. P.↗
NASA NTRS1994

The effect of wave-particle interactions on the polar wind: Preliminary results

A Monte Carlo simulation was developed in order to study the effects of wave-particle interactions (WPI) on the plasma outflow in the polar wind. The simulation also considered the other mechanisms included in the `classical' polar wind studies such as gravity, the polarization electrostatic field, and the divergence of geomagnetic field lines. Although the plasma consisted of electrons, H(+) and O(+) ions, we emphasized the behavior of H(+) in this preliminary study. The ion distribution function, as well as the profiles of its moments (density, drift velocity, temperatures, etc.) were found for different levels of WPI, that is, for different values of normalized diffusion rates in the velocity space (normalized D(sub perpendicular) H(+)). We found that as the WPI strength increases: (1) the ion drift velocity increases and its density decreases; (2) the perpendicular temperature T(sub perpendicular) (H(+)) increases; (3) the parallel temperature T(sub parallel) (H(+)) first decreases and then increases due to the balance between the parallel adiabatic cooling and the transfer of the energy from the perpendicular to the parallel direction; and (4) the temperature anisotropy ((T(sub parallel) (H(+))/T(sub perpendicular) (H(+)) is reduced and even reversed in some cases. For strong WPI (normalized D(sub perpendicular) (H(+) much greater than 1), the ion distribution function shows weak conic features at high altitudes.

Barakat, A. R.↗
NASA NTRS1994

Electron collection by a highly positive satellite in the ionosphere: Test aprticle simulation

Collection of electrons by a satellite at a relatively high potential in low Earth orbit, including the effects of the satellite's orbital motion, remains a poorly understood phenomenon. Using a test particle simulation in which charged particle motion is tracked in prescribed electric fields, we calculate here the current collection and the nature of the energy distribution function of the electrons collected by the satellite, including the effects of the satellite orbital motion. Calculations of the collected current without the orbital motion show an excellent agreement with the current from the Parker-Murphy model, but with the orbital motion of the satellite the current collection is enhanced, the degree of enhancement depending on the size of the satellite sheath extending along the magnetic field line. In the latter case, the flow of electrons around the satellite shows some interesting behaviors including azimuthal E x B drift around the satellite, axial trapping along the magnetic field, and formation of field-aligned flow of electrons in the wake region. The total energy of the collected electrons is approximately equal e phi(sub 0), where phi(sub 0) is the satellite potential, but the partition of the energy into components parallel (W(sub parallel)) and perpendicular (W(sub perpendicular)) to the magnetic field shows interesting features. For the magnetic field along the polar axis Z, the energy distributions near the poles are perfectly field aligned, that is, the parallel energy (W(sub parallel) = e phi(sub 0). The perpendicular component (W(sub perpendicular)) progressively increases toward the 'equator' of the spherical satellite. For a sufficiently large bias potential phi(sub 0), the equatorial energy distributions f(W(sub parallel)) and f(W(sub perpendicular)) show fine structures having multiple peaks. These features of the distribution functions result from the chaotic motion of electrons and the associated coupling between the electron motions parallel and perpendicular to the magnetic field in strongly nonuniform electric fields near the satellite.

Singh, Nagendra↗
NASA NTRS2001

Langmuir Probe and Mass Spectroscopic Measurements in Inductively Coupled CF4 Plasmas

Abstract Electron and ion energy distribution functions and other plasma parameters such as plasma potential (V(sub p)) , electron temperature (T(sub e)), and electron and ion number densities (n (sub e) and n(sub i)) in low pressure CF4 plasmas have been measured. The experiments were conducted in a GEC cell using an inductively coupled plasma (ICP) device powered by a 13.56 MHz radio-frequency (rf) power source. The measurements were made at 300 W of input rf power at 10, 30 and 50 mTorr gas pressures. Langmuir probe measurements suggest that n(sub e), n(sub i) and V(sub p) remain constant over 60% of the central electrode area, beyond which they decrease. Within the limits of experimental error (+/- 0.25 eV), T(sub e) remains nearly constant over the electrode area. T(sub e) and V(sub p) increase with a decrease in pressure. n(sub e) and n(sub i) are not affected as significantly as T(sub e) or V(sub p) by variation in the gas pressure. The electron energy distribution function (EEDF) measurements indicate a highly non-Maxwellian plasma. CF3+ is the most dominant ion product of the plasma, followed by CF2+ and CF+. The concentrations of CF2+ and CF+ are much larger than that is possible from direct electron impact ionization of the parent gas. The cross-section data suggest that the direct electron impact ionization of fragment neutrals and negative ion production by electron attachment may be responsible for increase of the minor ions.

Rao, M. V. V. S.↗
NASA NTRS1994

Acceleration of interstellar pickup ions in the disturbed solar wind observed on Ulysses

Acceleration of interstellar pickup H(+) and He(+) as well as of solar wind protons and alpha particles has been observed on Ulysses during the passage of a corotating interaction region (CIR) at approximately 4.5 AU. Injection efficiencies for both the high thermal speed interstellar pickup ions (H(+) and He(+)) and the low thermal speed solar wind ions (H(+) and He(++) are derived using velocity distribution functions of protons, pickup He(+) and alpha particles from less than 1 to 60 keV/e and of ions (principally protons) above approximately 60 keV. The observed spatial variations of the few keV and the few hundred keV accelerated pickup protons across the forward shock of CIR indicate a two stage acceleration mechanism. Thermal ions are first accelerated to speeds of 3 to 4 times the solar wind speed inside the CIR, presumably by some statistical mechanism, before reaching higher energies by a shock acceleration process. Our results also indicate that (1) the injection efficiencies for pickup ions are almost 100 times higher than they are for solar wind ions, (2) pickup H(+) and He(+) are the two most abundant suprathermal ion species and they carry a large fraction of the particle thermal pressure, (3) the injection efficiency is highest for protons, lowest for He(+), and intermediate for alpha particles, (4) both H(+) and He(+) have identical spectral shapes above the cutoff speed for pickup ions, and (5) the solar wind frame velocity distribution function of protons has the form F(w) = F(sub o)w(sup -4) for 1 is less than w is less than approximately 5, where w is the ion speed divided by the solar wind speed. Above w approximately 5-10 the proton spectrum becomes steeper. These results have important implications concerning acceleration of ions by shocks and CIRs, acceleration of anomalous cosmic rays, and particle dynamics in the outer heliosphere.

Gloeckler, G.↗
NASA NTRS1978

Theory of flux anisotropies in a guiding center plasma

The one particle distribution function f on the scale of the bounce motion of particles in a magnetic field B is considered. The Vlasov equation is expanded through O(epsilon) in the adiabatic parameter which is the ratio of particle gyroradius to scale length of the magnetic field. Because f is directly proportional to particle flux differential in kinetic energy and solid angle, f is in principle measurable in space experiments, and the analysis is tailored to be explicitly applicable to space problems. To O(1), f is gyrotropic; its first velocity moment is (if non-vanishing) parallel to B, and hence macroscopic parallel flow is included in this term. The O(epsilon) contribution is non-gyrotropic and macroscopic flow parallel to B plus additional parallel flow results from these terms. The degree of non-gyrotropy and the amount of cross-field macroscopic flow depend on the perpendicular component of the electric field, on curvature and shear in the magnetic field, and on the spatial gradient, pitch angle derivative, and speed derivative of the lowest order distribution function.

Birmingham, T. J.↗
NASA NTRS1989

Theory of droplet. Part 1: Renormalized laws of droplet vaporization in non-dilute sprays

The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.

Chiu, H. H.↗
NASA NTRS1998

Profile Evolution Simulation in Etching Systems Using Level Set Methods

Semiconductor device profiles are determined by the characteristics of both etching and deposition processes. In particular, a highly anisotropic etch is required to achieve vertical sidewalls. However, etching is comprised of both anisotropic and isotropic components, due to ion and neutral fluxes, respectively. In Ar/Cl2 plasmas, for example, neutral chlorine reacts with the Si surfaces to form silicon chlorides. These compounds are then removed by the impinging ion fluxes. Hence the directionality of the ions (and thus the ion angular distribution function, or IAD), as well as the relative fluxes of neutrals and ions determines the amount of undercutting. One method of modeling device profile evolution is to simulate the moving solid-gas interface between the semiconductor and the plasma as a string of nodes. The velocity of each node is calculated and then the nodes are advanced accordingly. Although this technique appears to be relatively straightforward, extensive looping schemes are required at the profile corners. An alternate method is to use level set theory, which involves embedding the location of the interface in a field variable. The normal speed is calculated at each mesh point, and the field variable is updated. The profile comers are more accurately modeled as the need for looping algorithms is eliminated. The model we have developed is a 2-D Level Set Profile Evolution Simulation (LSPES). The LSPES calculates etch rates of a substrate in low pressure plasmas due to the incident ion and neutral fluxes. For a Si substrate in an Ar/C12 gas mixture, for example, the predictions of the LSPES are identical to those from a string evolution model for high neutral fluxes and two different ion angular distributions.(2) In the figure shown, the relative neutral to ion flux in the bulk plasma is 100 to 1. For a moderately isotropic ion angular distribution function as shown in the cases in the left hand column, both the LSPES (top row) and rude's string method (bottom row) predict tapered profiles. The LSPES uses an AD with a FWHM = 13.5 degrees, and rude's model uses a ratio of sheath voltage to ion temperature of 50. The more anisotropic IADs produce profiles with more vertical sidewalls and a wider bottom surface, as shown in the right hand column. Here, the LSPES has an AD with a FWHM = 2 degrees, and rude's model uses a ratio of 500 of sheath voltage to ion temperature. The agreement between the LSPES and rude's model is excellent in both cases. We will present etching profiles generated by the LSPES, including calculations of the re-emitted fluxes of both neutrals and ions off of the profile walls. In addition, we will show the effect of geometric structures (overhangs, etc.) on the etching profiles. Other physical aspects, such as surface diffusion, will also be included in the model.

Hwang, Helen H.↗
NASA NTRS1997

Nonguiding Center Motion and Substorm Effects in the Magnetotail

Thick and thin models of the middle magnetotail were developed using a consistent orbit tracing technique. It was found that currents carried near the equator by groups of ions with anisotropic distribution functions are not well approximated by the guiding center expressions. The guiding center equations fail primarily because the calculated pressure tensor is not magnetic field aligned. The pressure tensor becomes field aligned as one moves away from the equator, but here there is a small region in which the guiding center equations remain inadequate because the two perpendicular components of the pressure tensor are unequal. The significance of nonguiding center motion to substorm processes then was examined. One mechanism that may disrupt a thin cross-tail current sheet involves field changes that cause ions to begin following chaotic orbits. The lowest-altitude chaotic region, characterized by an adiabaticity parameter kappa approx. equal to 0.8, is especially important. The average cross-tail particle drift is slow, and we were unable to generate a thin current sheet using such ions. Therefore, any process that tends to create a thin current sheet in a region with kappa approaching 0.8 may cause the cross-tail current to get so low that it becomes insufficient to support the lobes. A different limit may be important in resonant orbit regions of a thin current sheet because particles reach a maximum cross-tail drift velocity. If the number of ions per unit length decreases as the tail is stretched, this part of the plasma sheet also may become unable to carry the cross-tail current needed to support the lobes. Thin sheets are needed for both resonant and chaotic orbit mechanisms because the distribution function must be highly structured. A description of current continuity is included to show how field aligned currents can evolve during the transition from a two-dimensional (2-D) to a 3-D configuration.

Kaufmann, Richard L.↗
NASA NTRS1973

A general relaxation theory of simple liquids

A relatively simple relaxation theory to account for the behavior of liquids under dynamic conditions was proposed. The general dynamical equations are similar in form to the phenomenological relaxation equations used in theories of viscoelasticity, however, they differ in that all the coefficients of the present equations are expressed in terms of thermodynamic and molecular quantities. The theory is based on the concept that flow in a liquid distorts both the radial and the velocity distribution functions, and that relaxation equations describing the return of these functions to their isotropic distributions, characterizing a stationary liquid, can be written. The theory was applied to the problems of steady and oscillatory shear flows and to the propagation of longitudinal waves. In all cases classical results are predicted for strain rates, and an expression for the viscosity of a liquid, simular to the Macedo-Litovitz equation, is obtained.

Merilo, M.↗
NASA NTRS1990

Cumulative Poisson Distribution Program

Overflow and underflow in sums prevented. Cumulative Poisson Distribution Program, CUMPOIS, one of two computer programs that make calculations involving cumulative Poisson distributions. Both programs, CUMPOIS (NPO-17714) and NEWTPOIS (NPO-17715), used independently of one another. CUMPOIS determines cumulative Poisson distribution, used to evaluate cumulative distribution function (cdf) for gamma distributions with integer shape parameters and cdf for X (sup2) distributions with even degrees of freedom. Used by statisticians and others concerned with probabilities of independent events occurring over specific units of time, area, or volume. Written in C.

Bowerman, Paul N.↗
NASA NTRS1989

Semiclassical approximations in the coherent-state representation

The semiclassical limit of the stationary Schroedinger equation in the coherent-state representation is analyzed simultaneously for the groups W1, SU(2), and SU(1,1). A simple expression for the first two orders for the wave function and the associated semiclassical quantization rule is obtained if a definite choice for the classical Hamiltonian and expansion parameter is made. The behavior of the modulus of the wave function, which is a distribution function in a curved phase space, is studied for the three groups. The results are applied to the quantum triaxial rotor.

Kurchan, J.↗
NASA NTRS2015

Wear-Out Sensitivity Analysis Project Abstract

During the course of the Summer 2015 internship session, I worked in the Reliability and Maintainability group of the ISS Safety and Mission Assurance department. My project was a statistical analysis of how sensitive ORU's (Orbital Replacement Units) are to a reliability parameter called the wear-out characteristic. The intended goal of this was to determine a worst case scenario of how many spares would be needed if multiple systems started exhibiting wear-out characteristics simultaneously. The goal was also to determine which parts would be most likely to do so. In order to do this, my duties were to take historical data of operational times and failure times of these ORU's and use them to build predictive models of failure using probability distribution functions, mainly the Weibull distribution. Then, I ran Monte Carlo Simulations to see how an entire population of these components would perform. From here, my final duty was to vary the wear-out characteristic from the intrinsic value, to extremely high wear-out values and determine how much the probability of sufficiency of the population would shift. This was done for around 30 different ORU populations on board the ISS.

Harris, Adam↗