Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Analysis of hydrogen H-alpha observations of the coma of Comet P/Halley

The Monte Carlo Particle Trajectory Model of Combi and Smyth (1988) is used here to analyze observations of the H-alpha coma of Comet Halley. The solar excitation mechanism for the H-alpha emissions line is described. The H2O production rates derived for the H-alpha brightness measurements are shown to be very consistent with the H2O production rates determined from other Comet Halley observations of the H, O, and OH comae. Revised H2O production rates determined from 6300 A brightness measurements are presented.

Smyth, William H.↗
NASA NTRS1995

Observations and analysis of O(1D) and NH2 line profiles for the coma of comet P/Halley

A set of high-resolution Fabry-Perot measurements of the coma of comet P/Halley was acquired in the (O I) 6300 A and NH2 6298.62 A emission lines. These high-resolution measurements provide the first optical observations capable of studying directly the photochemical kinetics and dynamic outflow of the coma. The observations were analyzed by a Monte Carlo Particle Trajectory Model. The agreement of the model and observed line profiles was excellent and verified the underlying dynamics, exothermic photodissociative chemistry, and collisional thermalization in the coma. The somewhat wider intrinsic line profile width for the O(1D) emission in 1986 January compared to 1986 May, is, for example, produced by the larger outflow speeds and gas temperatures nearer perihelion in January. The January O(1D) profile, which is wider than the January NH2 profile, is indicative of the photochemical kinetics in the dissociation of the parent molecules H2O and OH in the coma. The absolute calibration of the observations in 1986 January allowed the production rates for H2O and the NH2-parent molecules to be determined. The average daily water production rates derived from the O(1D) emission data for January 16 and 17 are presented. These very large water production rates are consistent with the extrapolated (and 7.6 day time variable) water production rates determined from the analysis of lower spectral resolution observations for O(1D) and H-alpha emissions that covered the time period up to January 13. The large production rates on January 16 and 17 establish that the maximum water production rate for comet Halley accurred pre-perihelion in January. Implications drawn from comparison with 18 cm radio emission data in January suggest that the peak water production rate was even larger. The average production rate for NH3 determined from the NH2 emission data for January 17 was (1.48 +/- 0.10) x 10(exp 28) molecules/s, yielding an NH3/H2O production rate ratio of 0.55%.

Smyth, William H.↗
NASA NTRS1995

Analysis of hydrogen Lyman-alpha observations of the coma of Comet P/Halley near the perihelion

The pioneer Venus Orbiter Ultraviolet Spectrometer measurements of the Lyman-alpha intensity of atomic hydrogen excited by solar resonance scattering in the coma of Comet P/Halley acquired from December 28, 1985, to January 6, 1986, and from January 31, 1986, to March 6, 1986, are simulated with the Monte Carlo Particle Trajectory Model corrected for optical depth effects. Spatially detailed comparisons between data and model show excellent agreement and are used to infer that the highest cometary activity may not be at perihelion, but about 2 1/2 weeks before. An improved set of H2O production rates is presented for the period of time that the spectrometer was observing and is found to be consistent with the rates from other types of measurements. The apparent discrepancy between Stewart (1987) in early March and International Ultraviolet Explorer OH derived rates is resolved. The problem with the conversion of 18-cm OH radio brightness to H2O production rates is also discussed.

Smyth, William H.↗
NASA NTRS1995

Geomagnetic transmission of solar energetic protons during the geomagnetic disturbances of October 1989

Orbit-averaged geomagnetic transmission measurements during the large solar energetic particle events of October 1989 are presented using proton data from the NOAA-10 and GOES-7 satellies. The measurements are compared to geomagnetic transmission calculations determined by tracing particle trajectories through the combination of the International Geomagnetic Reference Field (IGRF) model and the 1989 Tsyganenko magnetospheric magnetic field model. The effective 'ring current' parameter in the 1989 Tsyganenko model based on the Dst data. Results are compared to calculations employing only the IGRF and to a parameterization of geomagnetically quiet-time cutoff rigidities derived from Cosmos/intercosmos observations. The 3-hour orbit-averaged results have approximately 15% accuracy during the October 1989 events.

Boberg, P. R.↗
NASA NTRS1995

Interactive data exploration and particle tracking for general circulation models

The SCENE environment for interactive visualization of complex data sets is discussed. This environment is used to create tools for graphical exploration of atmospheric flow models. These tools may be extended by the user in a seamless manner, so that no programming is required. A module for accurately tracing field lines and particle trajectories in SCENE is presented. This is used to examine the flowfield qualitatively with streamlines and pathlines and to identify critical points in the velocity field. The paper also describes a visualization tool for general circulation models on which the primary features of the environment are demonstrated.

Rosenbaum, R. I.↗
NASA NTRS1998

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗
NASA NTRS1998

Analysis of IUE Observations of Hydrogen in Comets

The 15-years worth of hydrogen Lyman-alpha observations of cometary comae obtained with the International Ultraviolet Explorer (IUE) satellite had gone generally unanalyzed because of two main modeling complications. First, the inner comae of many bright (gas productive) comets are often optically thick to solar Lyman-alpha radiation. Second, even in the case of a small comet (low gas production) the large IUE aperture is quite small as compared with the immense size of the hydrogen coma, so an accurate model which properly accounts for the spatial distribution of the coma is required to invert the infrared brightnesses to column densities and finally to H atom production rates. Our Monte Carlo particle trajectory model (MCPTM), which for the first time provides the realistic full phase space distribution of H atoms throughout the coma has been used as the basis for the analysis of IUE observations of the inner coma. The MCPTM includes the effects of the vectorial ejection of the H atoms upon dissociation of their parent species (H2O and OH) and of their partial collisional thermalization. Both of these effects are crucial to characterize the velocity distribution of the H atoms. This combination of the MCPTM and spherical radiative transfer code had already been shown to be successful in understanding the moderately optically thick coma of comet P/Giacobini-Zinner and the coma of comet Halley that varied from being slightly to very optically thick. Both of these comets were observed during solar minimum conditions. Solar activity affects both the photochemistry of water and the solar Lyman-alpha radiation flux. The overall plan of this program here was to concentrate on comets observed by IUE at other time during the solar cycle, most importantly during the two solar maxima of 1980 and 1990. Described herein are the work performed and the results obtained.

Combi, Michael R.↗
NASA NTRS1999

Three-Dimensional Simulations of Electron Beams Focused by Periodic Permanent Magnets

A fully three-dimensional (3D) model of an electron beam focused by a periodic permanent magnet (PPM) stack has been developed. First, the simulation code MAFIA was used to model a PPM stack using the magnetostatic solver. The exact geometry of the magnetic focusing structure was modeled; thus, no approximations were made regarding the off-axis fields. The fields from the static solver were loaded into the 3D particle-in-cell (PIC) solver of MAFIA where fully 3D behavior of the beam was simulated in the magnetic focusing field. The PIC solver computes the time-integration of electromagnetic fields simultaneously with the time integration of the equations of motion of charged particles that move under the influence of those fields. Fields caused by those moving charges are also taken into account; thus, effects like space charge and magnetic forces between particles are fully simulated. The electron beam is simulated by a number of macro-particles. These macro-particles represent a given charge Q amounting to that of several million electrons in order to conserve computational time and memory. Particle motion is unrestricted, so particle trajectories can cross paths and move in three dimensions under the influence of 3D electric and magnetic fields. Correspondingly, there is no limit on the initial current density distribution of the electron beam, nor its density distribution at any time during the simulation. Simulation results including beam current density, percent ripple and percent transmission will be presented, and the effects current, magnetic focusing strength and thermal velocities have on beam behavior will be demonstrated using 3D movies showing the evolution of beam characteristics in time and space. Unlike typical beam optics models, this 3D model allows simulation of asymmetric designs such as non- circularly symmetric electrostatic or magnetic focusing as well as the inclusion of input/output couplers.

Kory, Carol L.↗
NASA NTRS1994

Computational Investigation of an F-18 Aircraft in the High-Alpha Regime

One of the goals of NASA's High Alpha Technology Program is to provide flight-validated design methods for the high-angle-of-attack regime. This is an integrated effort utilizing computational simulations, wind tunnel experiments, and flight tests using the F-18 High Alpha Research Vehicle (HARV). The dominant physics of the aircraft flows in the high alpha regime changes as the angle of attack is increased. At moderate angle of attack the flow is characterized by boundary layer separation and the formation of tight vortices. As the angle of attack is increased, these vortices break down producing unsteady wakes. With further increase in angle of attack, the, vortex breakdown moves progressively upstream until the entire flowfield becomes dominated by the unsteady wake. Previous computational work has demonstrated the ability to simulate flows about the F-18 HARV in the medium-to-high angle of attack range, where the flowfield is characterized by the vortex formation and subsequent breakdown. This paper extends the previous computations to include conditions of 45 degree angle of attack where the flowfield becomes dominated by the unsteady wake shed from the Leading Edge Extension (LEX), and regions of laminar and transitional flow appear on the fuselage forebody. A more complete surface geometry is utilized, which includes the features of the engine nacelle, inlet diffuser, and the boundary layer diverter duct. A volume grid sensitivity study was also performed to extend the accuracy of the results, most notably in the prediction of the LEX vortex breakdown position. This paper includes comparisons of computational results with both in-flight surface pressure measurements, and flow visualizations of the surface and off-surface particle trajectories.

Murman, Scott M.↗
NASA NTRS2001

Investigation of Nitride Morphology After Self-Aligned Contact Etch

Self-Aligned Contact (SAC) etch has emerged as a key enabling technology for the fabrication of very large-scale memory devices. However, this is also a very challenging technology to implement from an etch viewpoint. The issues that arise range from poor oxide etch selectivity to nitride to problems with post etch nitride surface morphology. Unfortunately, the mechanisms that drive nitride loss and surface behavior remain poorly understood. Using a simple langmuir site balance model, SAC nitride etch simulations have been performed and compared to actual etched results. This approach permits the study of various etch mechanisms that may play a role in determining nitride loss and surface morphology. Particle trajectories and fluxes are computed using Monte-Carlo techniques and initial data obtained from double Langmuir probe measurements. Etched surface advancement is implemented using a shock tracking algorithm. Sticking coefficients and etch yields are adjusted to obtain the best agreement between actual etched results and simulated profiles.

Hwang, Helen H.↗
NASA NTRS1995

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗
NASA NTRS2003

Preliminary Results of a Microgravity Investigation to Measure Net Charge on Granular Materials

Accurate characterization of the electrostatic charge on granular materials has typically been limited to materials with diameters on the order of 10 microns and below due to high settling velocities of larger particles. High settling velocities limit both the time and the acceptable uncertainty with which a measurement can be made. A prototype device has been developed at NASA Glenn Research Center (GRC) to measure coulombic charge on individual particles of granular materials that are 50 to 500 microns in diameter. This device, a novel extension of Millikan's classic oil drop experiment, utilizes the NASA GRC 2.2 second drop tower to extend the range of electrostatic charge measurements to accommodate moderate size granular materials. A dielectric material with a nominal grain diameter between 1.06 and 250 microns was tribocharged using a dry gas jet, suspended in a 5x10x10 cm enclosure during a 2.2 second period of microgravity and exposed to a known electric field. The response was recorded on video and post processed to allow tracking of individual particles. By determining the particle trajectory and velocity, estimates of the coulombic charge were made. Over 30 drops were performed using this technique and the analysis showed that first order approximations of coulombic charge could successfully be obtained, with the mean charge of 3.4E-14 coulombs measured for F-75 Ottawa quartz sand. Additionally, the measured charge showed a near-Gaussian distribution, with a standard deviation of 2.14E -14 coulombs.

Green, Robert D.↗
NASA NTRS2004

NASCAP-2K: An Overview

Nascap-2k is the modern replacement for the older 3-D charging codes NASCAP/GEO, NASCAP/LEO, POLAR, and DynaPAC. Built on the DynaPAC kernel and incorporating surface charging, environment and space potential models from the older codes, Nascap-2k performs charging calculations for a wide variety of space environments under control of a unified graphical interface. In this paper we illustrate the use of Nascap-2k for spacecraft charging calculations. We touch on some of the unique physical and mathematical models on which the code is based. Examples/demos include the use of Object Toolkit, charging calculations in geosynchronous substorm, solar wind, low earth orbit, and auroral environments, and display and analysis of surface potentials, space potentials and particle trajectories.

Mandell, M. J.↗
NASA NTRS2004

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗
NASA NTRS2004

Accuracy of the Cloud Integrating Nephelometer

Potential error sources for measurements with the Cloud Integrating Nephelometer (CIN) are discussed and analyzed, including systematic errors of the measurement approach, flow and particle-trajectory deviations at flight velocity, ice-crystal breakup on probe surfaces, and errors in calibration and developing scaling constants. It is concluded that errors are minimal, and that the accuracy of the CIN should be close to the systematic behavior of the CIN derived in Gerber et al (2000). Absolute calibration of the CIN with a transmissometer operating co-located in a mountain-top cloud shows that the earlier scaling constant for the optical extinction coefficient obtained by other means is within 5% of the absolute calibration value, and that the CIN measurements on the Citation aircraft flights during the CRYSTAL-FACE study are accurate.

Gerber, Hermann E.↗
NASA NTRS2005

Flow Visualization Proposed for Vacuum Cleaner Nozzle Designs

In 1995, the NASA Lewis Research Center and the Kirby Company (a major vacuum cleaner company) began negotiations for a Space Act Agreement to conduct research, technology development, and testing involving the flow behavior of airborne particulate flow behavior. Through these research efforts, we hope to identify ways to improve suction, flow rate, and surface agitation characteristics of nozzles used in vacuum cleaner nozzles. We plan to apply an advanced visualization technology, known as Stereoscopic Imaging Velocimetry (SIV), to a Kirby G-4 vacuum cleaner. Resultant data will be analyzed with a high-speed digital motion analysis system. We also plan to evaluate alternative vacuum cleaner nozzle designs. The overall goal of this project is to quantify both velocity fields and particle trajectories throughout the vacuum cleaner nozzle to optimize its "cleanability"--its ability to disturb and remove embedded dirt and other particulates from carpeting or hard surfaces. Reference

Source record↗
NASA NTRS2005

Plasma Heating Simulation in the VASIMR System

The paper describes the recent development in the simulation of the ion-cyclotron acceleration of the plasma in the VASIMR experiment. The modeling is done using an improved EMIR code for RF field calculation together with particle trajectory code for plasma transport calculat ion. The simulation results correlate with experimental data on the p lasma loading and predict higher ICRH performance for a higher density plasma target. These simulations assist in optimizing the ICRF anten na so as to achieve higher VASIMR efficiency.

Ilin, Andrew V.↗