Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARTICLE TRAJECTORY”

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 505 records · Page 28

Observations of Subvisible Cirrus Clouds and Gravity Waves at the Tropical Tropopause

Thin, subvisible cirrus (SVC) clouds at the tropical tropopause have been observed by a number of methods in a variety of observational programs, including in situ sampling and aircraft and space-based lidar. Modeling studies suggest that these clouds play an important role in dehydrating tropospheric air as it enters the stratosphere. This is because particles large enough to have significant fall speeds can form under the conditions of slow cooling that are implied by the large horizontal extent of the SVC sheets. The IR radiation that these clouds absorb, and the upward vertical motion this implies, also makes them candidates for a tropical troposphere-to-stratosphere mass transfer mechanism. They may also play a role in the earth's radiation budget. These sheets were observed on five flights during the Tropical Ozone Transport Experiment (TOTE) by the NASA Langley DIAL lidar aboard NASA's DC-8 research aircraft operating during December 1995 and February 1996 south of Hawaii. The relationship of the SVC's observed during TOTE to convection was not a simple one. One class of SVC's are within 1000 km of the persistent strong convection near 15S (the SPCZ). Trajectory analyses indicated that the SVC air masses have in fact passed through the SPCZ within a few days of observation. These clouds are very close to the tropopause, with maximum potential temperatures not much higher than 370K, consistent with in situ water and total water measurements near the tropopause made during the Stratosphere Troposphere Exchange Project in January 1987 at Darwin, Australia. A second class of SVC's are not immediately downstream of convection. These clouds tend to be higher, reaching potential temperatures of 390K or more. Trajectory analyses indicate that the air in these SVC's originates either in the equatorial western Pacific or along the subtropical jet. In any case, the warm temperatures the SVC air masses encounter just prior to the observation time along the back trajectory imply that the clouds cannot be residual particles from cirrus blowoff, but must form locally as the air move upward and equatorward south of Hawaii. Since all the parcels have encountered colder temperatures than those at the time of observation early in their history, subsynoptic scale temperatures colder than the analysis temperatures appear to be required to explain the formation of ice particles. In fact, the sloping shapes of the SVC's do suggest that they are gravity or inertia-gravity waves. In situ meteorological measurements made by the ER-2 within a day of the DC-8 remote lidar observations show a gravity wave structure near the equator with an estimated period of about 30 hours. This is sufficiently long to allow large particles to form and fall out (thus allowing dehydration). Other ER-2 flights south of Hawaii at other times of year show gravity and inertia-gravity waves with a poleward wavenumber component and significant (5 degrees peak to peak) temperature perturbation.

Pfister, Leonhard↗

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

Ion kinematics in a plasma focus.

The results of numerical integrations of three-dimensional equations of motion of ions subject to given electric and magnetic fields are presented. The fields represent those which may exist in the pinch phase of the plasma focus, although here they depend only on the radial coordinate. The ions initially have Maxwellian velocity distributions, and their trajectories are interpreted in terms of single-particle constants of the motion. Two models of the axial electric field Ez are considered. For strong Ez away from the axis, there is a cyclotron acceleration which leads to ion heating. For positive Ez on the axis, ions within a Larmor radius of the axis undergo very efficient acceleration; the results for this new model are in general agreement with experimental results.

Gary, S. P.↗

Current status of the impact theory for the terminal Cretaceous extinction

An iridium abundance anomaly at the Cretaceous/Tertiary (C/T) boundary has been disclosed at 36 sites throughout the globe. The impact of a large extraterrestrial object is now widely accepted as the best explanation of this anomaly. Au/Ir and Pt/Ir ratios from two C/T boundary clays indicate a type I carbonaceous chondrite composition for the impacting object, for which a sea floor impact with subsequent subduction has been calculated to have a 20 percent probability. The extinction's 'killing mechanism' has not yet been established, but both temperature changes and darkness due to atmospheric dust are probable contributors. In view of the rapid dispersal of dust in ballistic trajectories, the rapid settling of heavy, coagulated dust particles, the effects of darkness on phytoplankton, and plant records, it is suggested that darkness would have lasted a few months, rather than the originally hypothesized few years.

Alvarez, W.↗

International cometary explorer solar array performance during the encounter with comet Giacobini-Zinner

The International Cometary Explorer (ICE) passed through the tail of comet Giacobini-Zinner (G-Z) on September 11, 1985. During the comet encounter, the ICE solar array was exposed to an environment consisting of hypervelocity cometary dust particles. The pre-encounter calculations, which provided verification that the ICE solar array would survive the comet encounter, are presented. The potential for damage to the ICE solar array by impacts from hypervelocity particles was examined for prospective closest-approach trajectories in the range of 1000 km to 10,000 km. The actual closest-approach trajectory was 8000 km from the nucleus of the comet. The ICE solar array experienced negligible degradation in performance during the encounter with comet G-Z.

Day, John H.↗

Effects of low frequency waves and spiky electric fields in the magnetotail

Intense low frequency waves and large amplitude spiky electric fields are commonly observed in the active magnetotail. These fields are of dynamic significance for substorms because the amplitude of the waves is large enough to provide the dissipation necessary for reconnection. The waves modify the trajectories of ions in the magnetotail, resulting in enhanced energization and pitch angle scattering in comparison with the trajectories obtained without waves. The spikes may represent the nonlinear evolution of the waves. Observations are presented of low frequency waves and spiky electric fields measured in the plasma sheet. The relationship between the waves and the spikes is compared to wave evolution theories. The Lundquist number is evaluated for the magnetotail based on the wave observations and previously published studies of the current sheet scale sizes. The effects of the waves on the ion trajectories are discussed using the results from a particle tracing code.

Cattell, C.↗

Strategies for Radiation Hardness Testing of Power Semiconductor Devices

Plans on the drawing board for future space missions call for much larger power systems than have been flown in the past. These systems would employ much higher voltages and currents to enable more powerful electric propulsion engines and other improvements on what will also be much larger spacecraft. Long term human outposts on the moon and planets would also require high voltage, high current and long life power sources. Only hundreds of watts are produced and controlled on a typical robotic exploration spacecraft today. Megawatt systems are required for tomorrow. Semiconductor devices used to control and convert electrical energy in large space power systems will be exposed to electromagnetic and particle radiation of many types, depending on the trajectory and duration of the mission and on the power source. It is necessary to understand the often very different effects of the radiations on the control and conversion systems. Power semiconductor test strategies that we have developed and employed will be presented, along with selected results. The early results that we have obtained in testing large power semiconductor devices give a good indication of the degradation in electrical performance that can be expected in response to a given dose. We are also able to highlight differences in radiation hardness that may be device or material specific.

Soltis, James V.↗

Wall effects observed in tissue-equivalent proportional counters from 1.05 GeV/nucleon iron-56 particles

Tissue-equivalent proportional counters (TEPCs) have been used to measure energy deposition in simulated volumes of tissue ranging in diameter from 0.1 to 10 microm. There has been some concern that the wall used to define the volume of interest could influence energy deposition within the sensitive volume because it has a density significantly greater than that of the cavity gas. These effects become important for high-velocity heavy ions. Measurements of energy deposition were made for 1 GeV/nucleon iron particles in a TEPC simulating a 1-microm-diameter sphere of tissue. The TEPC was nested within a particle spectrometer that provided identification and flight path of individual particles. Energy deposition was studied as a function of pathlength through the TEPC. Approximately 30% of the energy transfer along trajectories through the center of the detector escapes the sensitive volume. The response of the TEPC, for trajectories through the detector, is always larger than calculations for energy loss in a homogeneous medium. This enhancement is greatest for trajectories near the cavity/wall interface. An integration of the response indicates that charged-particle equilibrium is essentially achieved for a wall thickness of 2.54 mm. However, estimates of the linear energy transfer for the incident particles are influenced by these wall effects.

Non-NASA Center↗

Stardust Interstellar Preliminary Examination X: Impact Speeds and Directions of Interstellar Grains on the Stardust Dust Collector

On the basis of an interstellar dust model compatible with Ulysses and Galileo observations, we calculate and predict the trajectories of interstellar dust (ISD) in the solar system and the distribution of the impact speeds, directions, and flux of ISD particles on the Stardust Interstellar Dust Collector during the two collection periods of the mission. We find that the expected impact velocities are generally low (less than 10 km per second) for particles with the ratio of the solar radiation pressure force to the solar gravitational force beta greater than 1, and that some of the particles will impact on the cometary side of the collector. If we assume astronomical silicates for particle material and a density of 2 grams per cubic centimeter, and use the Ulysses measurements and the ISD trajectory simulations, we conclude that the total number of (detectable) captured ISD particles may be on the order of 50. In companion papers in this volume, we report the discovery of three interstellar dust candidates in the Stardust aerogel tiles. The impact directions and speeds of these candidates are consistent with those calculated from our ISD propagation model, within the uncertainties of the model and of the observations.

interstellar↗

NAIRAS Ionizing Radiation Model: Extension from Atmosphere to Space

The Nowcast of Aerospace Ionizing RAdiation System (NAIRAS) model is a real-time, global, physics-based model originally developed to predict exposure from cosmic radiation to air travelers from both galactic and solar sources. A prototype operational NAIRAS model has provided tabular and graphical data products via its public web site for about ten years. A new version of the NAIRAS model has been developed that incorporates an extension of the model domain from the atmospheric ionizing radiation environment to the space radiation environment, with the addition of the trapped inner belt proton source and altitude-dependent and rigidity-dependent geomagnetic shielding of the galactic cosmic rays (GCR) and solar energetic particle (SEP) protons. New output products of differential and integral particle flux have been developed for the characterization of single-event effects (SEE), expanding the application of NAIRAS from human radiation exposure assessment to allowing end-users to quantify radiation environment risks to aviation and spacecraft microelectronic systems. The NAIRAS model has transitioned to prototype operations at the Community Coordinated Modeling Center (CCMC) where the model now operates in two modes: (1) real-time global predictions of the atmospheric radiation environment and (2) a run-on-request (RoR) service allowing the user to select a specific time period for the global dosimetric calculations, or to upload an aircraft, balloon, or spaceflight trajectory file to provide predictions of the dosimetric and particle flux quantities along the flight path. The new features of NAIRAS version 3.0 are described in this paper and example results of the new output products for low-Earth orbit (LEO), medium-Earth orbit (MEO), and free-space radiation environments are presented

Christopher J. Mertens↗

Motion about the stable libration points in the linearized, restricted three-body problem

The motion of a point particle in the neighborhood of a triangular libration point (L sub 4 or L sub 5) in the linearized, restricted problem of three bodies in the plane is described. The derivation of the equations of motion is standard. From these equations, three invariants of the motion are obtained; the Jacobi integral is expressed linearly in terms of two of these. The trajectories for varied initial conditions are drawn, and a complete geometric description of the particle motion is given in elementary terms. Each trajectory has an exterior boundary curve; its equation is found. An approximation to this boundary curve was known; the two curves are compared graphically. For certain initial conditions, there is an interior region from which the trajectory is excluded; the equation of the boundary of this region is given.

Mittleman, D.↗

Putting quasilinear theory to the test: A comparison of analytical theory with results from a trajectory code

Parallel and perpendicular diffusion coefficients were computed numerically by following particle orbits in a simulated magnetic field. The simulated field was chosen to have delta B/B(sub o) small, so as to provide a test of quasilinear theory in a regime where the theory should be most accurate. The simulation space is large enough to contain many magnetic field correlation lengths, so that effects of field line random walk can be studied. After presenting results for parallel diffusion, we will focus on two controversial issues relating to perpendicular diffusion: (1) Do quasilinear descriptions of perpendicular diffusion retain any validity for particles whose Larmor radius is smaller than a correlation length? (2) Does field line random walk lead to particle diffusion in the usual sense, or does it produce 'compound' diffusion for which particles spread out proportionally to t(exp 1/4) instead of t(exp 1/2)?

Bieber, J. W.↗

Impingement of Droplets in 60 Deg Elbows with Potential Flow

Trajectories were determined for water droplets or other aerosol particles in air flowing through 600 elbows especially designed for two-dimensional potential motion. The elbows were established by selecting as walls of each elbow two streamlines of a flow field produced by a complex potential function that establishes a two-dimensional flow around. a 600 bend. An unlimited number of elbows with slightly different shapes can be established by selecting different pairs of streamlines as walls. Some of these have a pocket on the outside wall. The elbows produced by the complex potential function are suitable for use in aircraft air-inlet ducts and have the following characteristics: (1) The resultant velocity at any point inside the elbow is always greater than zero but never exceeds the velocity at the entrance. (2) The air flow field at the entrance and exit is almost uniform and rectilinear. (3) The elbows are symmetrical with respect to the bisector of the angle of bend. These elbows should have lower pressure losses than bends of constant cross-sectional area. The droplet impingement data derived from the trajectories are presented along with equations so that collection efficiency, area, rate, and distribution of droplet impingement can be determined for any elbow defined by any pair of streamlines within a portion of the flow field established by the complex potential function. Coordinates for some typical streamlines of the flow field and velocity components for several points along these streamlines are presented in tabular form. A comparison of the 600 elbow with previous calculations for a comparable 90 elbow indicated that the impingement characteristics of the two elbows were very similar.

Hacker, Paul T.↗

Sources, Transport and Visibility Impact of Ambient Submicrometer Particle Size Distributions in an Urban Area of Central Taiwan

This study applied positive matrix factorization (PMF) to identify the sources of size-resolved submicrometer (10–1000 nm) particles and quantify their contributions to impaired visibility based on the particle number size distributions (PNSDs), aerosol light extinction (b p ), air pollutants (PM 10 , PM 2.5 , SO 2 , O 3 , and NO), and meteorological parameters (temperature, relative humidity, wind speed, wind direction, and ultraviolet index) measured hourly over an urban basin in central Taiwan between 2017 and 2021. The transport of source-specific PNSDs was evaluated with wind and back trajectory analyses. The PMF revealed six sources to the total particle number (TPN), surface (TPS), volume (TPV), and b p . Factor 1 (F1), the key contributor to TPN (35.0 %), represented nucleation (<25 nm) particles associated with fresh traffic emission and secondary new particle formation, which were transported from the west-southwest by stronger winds (>2.2 m s -1 ). F2 represented the large Aitken (50–100 nm) particles transported regionally via northerly winds, whereas F3 represented large accumulation (300–1000 nm) particles, which showed elevated concentrations under stagnant conditions (<1.1 m s −1 ). F4 represented small Aitken (25–50 nm) particles arising from the growth and transport of the nucleation particles (F1) via west-southwesterly winds. F5 represented large Aitken particles originating from combustion-related SO 2 sources and carried by west-northwesterly winds. F6 represented small accumulation (100–300 nm) particles emitted both by local sources and by the remote SO 2 sources found for F5. Overall, large accumulation particles (F3) played the greatest role in determining the TPV (66.4 %) and TPS (34.8 %), and their contribution to bp increased markedly from 17.3 % to 40.7 % as visibility decreased, indicating that TPV and TPS are better metrics than TPN for estimating b p . Furthermore, slow-moving air masses—and therefore stagnant conditions—facilitate the build-up of accumulation mode particles (F3 + F6), resulting in the poorest visibility.

source apportionment↗

Modeling the Superstorm in November 2003

The superstorm on 20.21 November 2003 was the largest geomagnetic storm in solar cycle 23 as measured by Dst, which attained a minimum value of .422 nT. We have simulated this storm to understand how particles originating from the solar wind and ionosphere get access to the magnetosphere and how the subsequent transport and energization processes contribute to the buildup of the ring current. The global electromagnetic configuration and the solar wind H+ distribution are specified by the Lyon-Fedder-Mobarry (LFM) magnetohydrodynamics model. The outflow of H+ and O+ ions from the ionosphere are also considered. Their trajectories in the magnetosphere are followed by a test-particle code. The particle distributions at the inner plasma sheet established by the LFM model and test-particle calculations are then used as boundary conditions for a ring current model. Our simulations reproduce the rapid decrease of Dst during the storm main phase and the fast initial phase of recovery. Shielding in the inner magnetosphere is established at early main phase. This shielding field lasts several hours and then breaks down at late main phase. At the peak of the storm, strong penetration of ions earthward to L shell of 1.5 is revealed in the simulation. It is surprising that O+ is significant but not the dominant species in the ring current in our calculation for this major storm. It is very likely that substorm effects are not well represented in the models and O+ energization is underestimated. Ring current simulation with O+ energy density at the boundary set comparable to Geotail observations produces excellent agreement with the observed symH. As expected in superstorms, ring current O+ is the dominant species over H+ during the main to mid-recovery phase of the storm.

Fok, Mei-Ching↗

Model exospheres of the ringed planets

The theory of rotating ion exospheres in dipolar magnetic field geometries is used to predict certain effects of the rings on the inner plasmaspheres of the ringed planets Jupiter, Saturn, and Uranus. The analysis here assumes a purely ionospheric source. It is found that if rings affect only the trapped plasma, the cold plasma density at Saturn will be either reduced or increased by a factor of approximately 2 at L = 1.65, depending on whether the ring populates or empties the trapped trajectories. It is found that if all cold plasma particles on trapped orbits are absorbed and the rings are not a source of plasma, this simple exospheric theory can explain the ionospheric density profile inferred from the Pioneer 11 Saturn radio occultation experiment, in which a localized peak in the profile was observed on field lines threading the Guerin division between the C and D rings.

Luhmann, J. G.↗

Cloud Condensation Nuclei in Fire-3

The centerpiece of this research was the cloud condensation nuclei (CCN) measurements of the Desert Research Institute (DRI) CCN spectrometers on board the NCAR C-130 aircraft during the Arctic Cloud Experiment (ACE) in May, 1998. These instruments operated successfully throughout all eight 10-hour research flights based in Fairbanks and the two ferry flights between Colorado and Fairbanks. Within a few months of completion of ACE the CCN data was edited and put into the archives. A paper was completed and published on the CCN climatology during the previous two FIRE field projects-FIRE 1 based in San Diego in June and July, 1987 and ASTEX based in the Azores Islands in June, 1992. This showed distinct contrasts in concentrations and spectra between continental and maritime CCN concentrations, which depended on air mass trajectories. Pollution episodes from Europe had distinct influences on particle concentrations at low altitudes especially within the boundary layer. At higher altitudes concentrations were similar in the two air mass regimes. Cloudier atmospheres showed lower concentrations especially below the clouds, which were a result mostly of coalescence scavenging.

Source record↗