Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle Drag”

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 271 records · Page 15

The role of particle collisions in pneumatic transport

A model of dilute gas-solid flow in vertical risers is developed in which the particle phase is treated as a granular material, the balance equations for rapid granular flow are modified to incorporate the drag force from the gas, and boundary conditions, based on collisional exchanges of momentum and energy at the wall, are employed. In this model, it is assumed that the particle fluctuations are determined by inter-particle collisions only and that the turbulence of the gas is unaffected by the presence of the particles. The model is developed in the context of, but not limited to, steady, fully developed flow. A numerical solution of the resulting governing equations provides concentration profiles generally observed in dilute pneumatic flow, velocity profiles in good agreement with the measurements of Tsuji, et al. (1984), and an explanation for the enhancement of turbulence that they observed.

Mastorakos, E.↗

Radial pressure in the solar nebula as affecting the motions of planetesimals

Growing planetesimals and a range of drag laws depending on the Reynolds number and on the ratio of particle size to mean free path are considered. Particles spiral in the direction of positive gradient, thus being concentrated toward toroidal concentrations of gas. The effect increases with decreasing rates of particle growth, i.e., with increasing time scales of planet formation by accretion. In the outer regions, where evidence suggests that comets were formed and Uranus and Neptune were so accumulated, the effect of the pressure gradient is to clear the forming comets from those regions. The large mass of Neptune may have developed because of this effect, perhaps Neptune's solar distance was reduced from Bode's law, and perhaps no comet belt exists beyond Neptune. In the asteroid belt, on a slow time scale, the effect may have spiraled planetesimals toward Mars and Jupiter, thus contributing to the lack of planet formation in this region.

Whipple, F. L.↗

Planetesimal Formation in the Outer Solar Nebula

A numerical investigation of the orbital trajectories of individual particles in the turbulent outer solar nebula has been performed. The (spherical) particle consists of an unchanging mm-sized 'dust' core surrounded by an H2O ice mantle; the density of both core and mantle is 0.5 g/cm(exp 3). The simulations include the effects of H2O condensation from the gas phase, H2O sublimation from the particle surface, and collisional growth via particle collisions with a background distribution of small H2O grains. The model nebula is an azimuthally symmetric minimum-mass nebula of solar composition with a vertical (and radial) temperature gradient. Particle evolution follows a pattern. A particle starting out in a cool region grows via condensation and collisional accretion until it is large enough (decimeter- to meter-sized) to decouple somewhat from the turbulence. (This growth occurs on a timescale of several thousand years at 10 AU; at 30 AU, the timescale is approx. 104 years.) The particle then moves rapidly inward toward the sun due to secular gas drag forces, sublimates much of its icy mantle, and slows its inward migration as it gets caught up in the turbulence again (due to its now-smaller size) at the 'sublimation boundary,' where the ambient gas temperature is approx. 150 K. Such a process could, on a short timescale (i.e., a timescale much shorter than the nebular gas lifetime of approx. 106 yr), generate a population of decimeter- to meter-sized bodies which would then collisionally accrete to form planetesimals.

Supulver, K. D.↗

Lorentz forces on the dust in Jupiter's ring

The paths of dust particles in the Jovian ring are investigated using a numerical integration program, including the acceleration due to gravity and the Lorentz and drag accelerations arising from the motions of the charged dust through the Jovian plasma. It is determined that the orbit of a 2.5 micron radius spherical dust particle with a density of 2 g/cu cm -10V will become significantly perturbed. The ring will tend to warp northwards near 130-160 deg longitude, with the maximum excursion of the Jupiter ring grains equalling about 0.1 deg (consistent with a distance of 220 km above the equatorial plane). It is found that either the particles are larger or the voltages on them less than what has been determined by previous investigators, while the plasma near the ring may be considerably cooler than was estimated. Calculations show that particles of 0.3 micron with -10 V potentials are spread from 1.68-1.98 of the radius of Jupiter and inclined up to 7 deg out of the equatorial plane. The paths of these particles do not follow Keplerian orbits, and the particle positions are not symmetric about the equatorial plane. Particles of 0.4 micron radius have less asymmetric orbits than 0.3 micron particles, while particles less than 0.2 micron are perturbed into Jupiter cloudtops within a few tens of hours.

Consolmagno, G. J.↗

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Goulds Belt, Interstellar Clouds, and the Eocene Oligocene Helium-3 Enhancement

Drag from hydrogen in the interstellar cloud which formed Gould's Belt may have sent interplanetary dust particle (IDPs) and small meteoroids with embedded helium to the Earth, perhaps explaining part the helium-3 flux increase seen in the sedimentary record near the Eocene-Oligocene transition. Assuming the Solar System passed through part of the cloud, IDPs in the inner Solar System may have been dragged to Earth, while dust and small meteoroids in the asteroid belt up to centimeter size may have been dragged to the resonances, where their orbital eccentricities were pumped up into Earth-crossing orbits; however, this hypotheses does not explain the Popigai and Chesapeake Bay impacts.

Rubincam, David Parry↗

The rings of Uranus - Nature and origin

Effects of a resonance between a set of ring particles and individual satellites are considered. It is suggested that the resonances are 1:1 for the rings of Uranus and that each ring must contain an as-yet unseen satellite. A model is proposed in which the rings of Uranus were caused by small satellites that entered the Roche zone due to tidal drag. According to this model, the satellites gradually lost material from their surfaces in this zone, the particles left a satellite with relatively low speeds and must have entered orbits similar to that of the satellite, and the gravitational force of the satellite compelled these particles to remain in a narrow sharply defined ring.

Dermott, S. F.↗

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

Optimization of a Remote Sensing Energy Balance Method over Different Canopy Applied at Global Scale

Parameterization methods which calculate turbulent heat and water fluxes with thermal remote sensing data were evaluated in the revised remote sensing surface energy balance system (SEBS) model (Chen et al., 2013). The model calculates sensible heat (H) based on the Monin-Obukhov similarity theory (MOST) and determines latent heat (LE) as the residual of energy balance. We examined the uncertainties of H and LE in the SEBS model due to five key parameters at the local station point scale. Observations at 27 flux towers located in seven land cover types (needle-leaf forest, broad leaf forest, shrub, savanna, grassland, cropland, and sparsely vegetated land) and an artificial intelligence particle swarm optimization (PSO) algorithm was combined to calibrate the five parameters (leaf drag coefficient, leaf heat transfer coefficients, roughness length for soil, and two parameters for ground heat calculation) in the SEBS model. The root-mean-square error at the site scale was reduced by 9 W/sq.m for H, and 92 W/sq.m for LE, and their correlation coefficients were increased by 0.07 (H) and 0.11 (LE) after using the calibrated parameters. The updated model validation was further conducted globally for the remotely sensed evapotranspiration (ET) calculations. Overestimation of SEBS global ET was significantly improved by using the optimized values of the parameters. The results suggested PSO was able to consistently locate the global optimum of the SEBS model, and appears to be capable of solving the ET model optimization problem.

Chen, Xuelong↗

Creeping flow of a conducting fluid past axisymmetric bodies in the presence of an aligned magnetic field

The use of strong magnetic fields for the control of particle settling in metallic systems is investigated by altering the fluid mechanics in the melt. The fluid mechanism of particle settling is analyzed for the motion around single, axisymmetric particles in the limit of creeping flow for a fluid with a large electrical conductivity. The drag is found to increase proportionately to the intensity of the magnetic field or the Hartmann number Ha. The flowfield forms boundary layers, which thin out with increasing Ha, along the surfaces parallel to the flow. For axisymmetric bodies, the boundary layer separates as the poles of the surface are approached and encloses regions of almost stagnant fluid. These regions spread upstream and downstream along the body with increasing Ha, thereby trapping the particle in a column of stagnant fluid.

Kyrlidis, A.↗

Erosion of circumstellar particle disks by interstellar dust

Circumstellar particle disks appear to be a common phenomenon; however, their properties vary greatly. Models of the evolution of such systems focus on internal mechanisms such as interparticle collisions and Poynting-Robertson drag. Herein it is shown that 'sandblasting' by interstellar dust can be an important and even dominant contributor to the evolution of circumstellar particle disks. Stars spend up to about 3 percent of their main-sequence lifetimes within atomic clouds. Among an IRAS sample of 21 nearby main-sequence A stars, beta Pictoris has the brightest disk; it also possesses the smallest random velocity and therefore the slowest predicted erosion rate.

Lissauer, Jack J.↗

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Particle trajectory computer program for icing analysis of axisymmetric bodies - A progress report

Aircraft exposed to an atmospheric icing environment can accumulate ice, resulting in a sharp increase in drag, a reduction in lift, control surface fouling, and engine damage all of which result in a hazardous flight situation. NASA Lewis Research Center (LeRC) has conducted a program to examine, with the aid of high-speed computer codes, how the trajectories of particles contribute to the ice accumulation on airfoils and engine inlets. For this effort, a computer code was developed to calculate icing particle trajectories and impingement limits for axisymmetric inlets. The original research-oriented NASA code was upgraded and modified to meet the requirements of the design engineer. The improved code is capable of performing trajectory calculations for any atmospheric conditions and droplet sizes. It can handle single droplets or a distribution of various droplet sizes. The four programs that comprise the code are described and the results of a test case using flight conditions for a Fokker F100 icing tunnel test are presented.

Maltezos, Dimitrios G.↗

Stratospheric Collection of Dust from Comet 73P/Schwassmann-Wachmann 3

Interplanetary dust particles (IDPs) collected in the stratosphere are unique materials that are compositionally distinct from meteorites. Astronomical observations and dynamical models indicate that both asteroids and short-period comets are significant sources of IDPs. IDPs having fragile, porous structures, unequilibrated, anhydrous mineralogy, and high atmospheric entry velocities are thought to derive from comets, whereas asteroidal IDPs are identified by their compact structure, hydrated mineralogy and low atmospheric entry velocities. Uncertainty remains in the classification of asteroidal and cometary IDPs owing to our limited sampling of comets and the asteroid belt and the complex dynamical histories of most IDPs in space. Most IDPs spend thousands of years in space prior to being accreted by the Earth. During this time, dust particles undergo orbital evolution, including gradual reduction in their perihelion and eccentricity as a result of Poynting-Robertson drag. Planetary encounters may also significantly change their orbital parameters. Consequently, it is generally not possible to identify the specific parent body of a given IDP. However, it has been proposed that it is possible to identify dust from comets that have formed Earth-crossing dust trails. In this case, the dust particles have been in space for such a short period of time (a few decades or less) that their orbits have not significantly changed. Furthermore, these fresh IDPs could be identified in the laboratory from their short space-exposure histories (low solar noble gas abundance and lack of solar flare tracks). NASA flew several dedicated IDP collection missions attempting to collect dust from comet 26P/Grigg-Skjellerup, the best candidate identified. Remarkably, many particles from those collectors exhibit unusual properties, including low abundances of solar noble gases and high abundances of presolar grains. These observations are consistent with the dust particles originating from comet Grigg-Skjellerup (hereafter G-S). This study considers the prospects for collection of dust from comet 73P/Schwassmann-Wachmann 3 (hereafter SW3). SW3 is a small (2 km diameter) Jupiter family comet whose perihelion is close to and just inside the Earth's orbit. The orbit of SW3 is suitable for producing a low-velocity Earth-crossing dust stream and is the likely parent of the Tau Herculid meteor stream. This study complements a previously published model of the SW3 meteor stream that predicted a very low level of activity for grains 100 micron -- 100 mm in size.

Messenger, Scott R.↗

Regimes of Pulsar Pair Formation and Particle Energetics

We investigate the conditions required for the production of electron-positron pairs above a pulsar polar cap (PC) and the influence of pair production on the energetics of the primary particle acceleration. Assuming space-charge limited flow acceleration including the inertial frame-dragging effect, we allow both one-photon and two-photon pair production by either curvature radiation (CR) photons or photons resulting from inverse-Compton scattering of thermal photons from the PC by primary electrons. We find that,, while only the younger pulsars can produce pairs through CR, nearly all known radio pulsars are capable of producing pairs through non-resonant inverse-Compton scatterings. The effect of the neutron star equations of state on the pair death lines is explored. We show that pair production is facilitated in more compact stars and more a massive stars. Therefore accretion of mass by pulsars in binary systems may allow pair production in most of the millisecond purser population. We also find that two-photon pair production may be important in millisecond pursers if their surface temperatures are above approx. or equal to three million degrees K. Pursers that produce pairs through CRT wilt have their primary acceleration limited by the effect of screening of the electric field. In this regime, the high-energy luminosity should follow a L(sub HE) proportional to dot-E(sup 1/2, sub rot) dependence. The acceleration voltage drop in pursers that produce pairs only through inverse-Compton emission will not be limited by electric field screening. In this regime, the high-energy luminosity should follow a L(sub HE) proportional to dot-E(sub rot) dependence. Thus, older pursers will have significantly lower gamma-ray luminosity.

Harding, Alice K.↗