Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete particle method”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A dynamic study of fragmentation and energy loss during high velocity impact

Research conducted under this contract can be divided into two main areas: hypervelocity (in the range up to 7 km/s) and high velocity (less than 1 km/s). Work in the former was performed at NASA-MSFC using the Light Gas Gun Facility. The lower velocity studies were conducted at Auburn University using the ballistic gun. The emphasis of the project was on the hypervelocity phenomenon especially in the characterization of the debris cloud formed by the primary impact events. Special devices were made to determine the angular distributions of momentum and energy of the debris cloud as a function of impact conditions. After several iteration processes, it was decided to concentrate on the momentum effort. Prototype devices were designed, fabricated, and tested. These devices were based on the conservation of momentum. Distributions of the debris cloud formed were measured by determining the amount of momentum transferred from the debris cloud to strategically placed pendulum measurement devices. The motion of the pendula was monitored using itegrated opto-interrupters. The distribution of momentum in the debris cloud was found to be a strong function of the impact condition. Small projectiles at high velocities were observed to produce finely dispersed debris whereas large projectiles generated discrete particles in the debris. Results also show that the momentum in the forward direction was enhanced due to the impact. This phenomenon of momentum multiplication was also observed in other studies and in computer simulations. It was initially planned to determine the energy distribution using deformation energy in a rod with strain gauges. Results from preliminary studies show that this technique is acceptable but too tedious. A new technique was explored based on measuring the heating effect of the debris cloud using an IR camera. The feasibility and sensitivity was established at Auburn University. This type of energy distribution measurement method can easily be adapted to the gas gun facility at MSFC. The objective of the lower velocity studies at Auburn was to simulate the damage produced in advanced materials by the lower energy debris cloud.

Zee, Ralph H.↗
NASA NTRS1995

Sensitivity analysis in multipole-accelerated panel methods for potential flow

In the design of an airframe, the effect of changing the geometry on resulting computations is necessary for design optimization. The geometry is defined in terms of a series of design variables, including design variables to define the wing planform, tail, canard, pylon, and nacelle. Design optimization in this research is based on how these design variable affect the potential flow. The potential flow is computed as a function of the geometry and location of a series of panels describing the airframe, which are in turn a function of the design variables. Multipole accelerated panel methods improve the computational complexity of the problem and thus are an attractive approach. To utilize the methods in design optimization, it was necessary to define the appropriate sensitivity derivatives. The overhead incurred from finding the sensitivity derivatives in conjunction with the original computation should be small. This research developed the background for multipole-accelerated panel methods and the framework for finding sensitivity derivatives in the methods. Potential flow panel codes are commonly used for powered-lift aerodynamic predictions for three dimensional geometries. Given an airframe which has been discretized into a series of panels to define the airframe geometry, potential is computed as a function of the influence of all panels on all other panels. This is a computationally intensive problem for which efficient solutions are desired to improve the computational time and to allow greater resolution by use of more panels. One such solution is the use of hierarchical multipole methods which entail approximations of the effects of far-field terms. Hierarchical multipole methods have become prevalent in molecular dynamics and gravitational physics, and have been introduced into the fields of capacitance calculations, computational fluid dynamics, and electromagnetics. The methods utilize multipole expansions to describe the effect of bodies (i.e. particles, astrophysical bodies, panels, etc.) within a sphere on points distant from the sphere, where the influence diminishes as a function of distance. The expansions are exact with infinite series, however, for practical computations, the series are truncated and accuracy is selected based on the number of terms retained in the expansions. A hierarchical tree structure groups bodies together based on proximity to allow definition of multipole expansions for each group. The multipole expansions are then used to compute the effect of the bodies in a group on distant bodies.

Leathrum, James F., Jr.↗
NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRS2017

Take-off Engine Particle Emission Indices for In-Service Aircraft at Los Angeles International Airport

We present ground‐based, advected aircraft engine emissions from flights taking off at Los Angeles International Airport. 275 discrete engine take‐off plumes were observed on 18 and 25 May 2014 at a distance of 400 m downwind of the runway. CO2 measurements are used to convert the aerosol data into plume‐average emissions indices that are suitable for modelling aircraft emissions. Total and non‐volatile particle number EIs are of order 1016‐1017 kg‐1 and 1014‐1016 kg‐1, respectively. Black‐carbon‐equivalent particle mass EIs vary between 175‐941 mg kg‐1 (except for the GE GEnx engines at 46 mg kg‐1). Aircraft tail numbers recorded for each take‐off event are used to incorporate aircraft‐ and engine‐specific parameters into the data set. Data acquisition and processing follow standard methods for quality assurance. A unique aspect of the data set is the mapping of aerosol concentration time series to integrated plume EIs, aircraft and engine specifications, and manufacturer‐reported engine emissions certifications. The integrated data enable future studies seeking to understand and model aircraft emissions and their impact on air quality.

Richard H Moore↗
NASA NTRS1994

The stability of perfect elliptic disks. 1: The maximum streaming case

Self-consistent distribution functions are constructed for two-dimensional perfect elliptic disks (for which the potential is exactly integrable) in the limit of maximum streaming; these are tested for stability by N-body integration. To obtain a discrete representation for each model, simulated annealing is used to choose a set of orbits which sample the distribution function and reproduce the required density profile while carrying the greatest possible amount of angular momentum. A quiet start technique is developed to place particles on these orbits uniformly in action-angle space, making the initial conditions as smooth as possible. The roundest models exhibit spiral instabilities similar to those of cold axisymmetric disks; the most elongated models show bending instabilities like those seen in prolate systems. Between these extremes, there is a range of axial ratios 0.25 approximately less than b/a approximately less than 0.6 within which these models appear to be stable. All the methods developed in this investigation can easily be extended to integrable potentials in three dimensions.

Levine, Stephen E.↗
NASA NTRS1993

Viscous analysis of three-dimensional rotor flows using a multigrid method

A three-dimensional code for rotating blade-row flow analysis was developed. The space discretization uses a cell-centered scheme with eigenvalues scaling for the artificial dissipation. The computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. An application is presented for the NASA rotor 67 transonic fan. Due to the blade stagger and twist, a zonal, non-periodic H-type grid is used to minimize the mesh skewness. The calculation is validated by comparing it with experiments in the range from the maximum flow rate to a near-stall condition. A detailed study of the flow structure near peak efficiency and near stall is presented by means of pressure distribution and particle traces inside boundary layers.

Arnone, A.↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRSDate not supplied

Using Polarized PCRTM-Solar Model to Validate EMIT Surface Reflectance L2A Product

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer scattering code, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate hyperspectral radiation fields from visible to shortwave infrared (VSWIR), i.e., 250 nm to 2500 nm. Since only limited accurate line-by-line calculation using VLIDORT is needed, this coupled PCRTM-VLIDORT approach is very efficient for simulating hyperspectral VSWIR spectra. These simulations incorporate accurately defined optical properties (gas absorptions, molecular and aerosol/cloud particle scattering) and detailed surface reflectance inputs. In this study we used the current L2A products of the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as the inputs to PCRTM-VLIDORT to simulate EMIT spectra. Comparing the simulated spectra with EMIT measurements provides an efficient method to validate the current EMIT product.

Xiaozhen Xiong↗
NASA NTRS2019

Retrieving Aerosol Characteristics From the PACE Mission, Part 1: Ocean Color Instrument

NASA’s Plankton, Aerosol, Clouds, ocean Ecosystem (PACE) satellite mission is scheduled to launch in 2022, with the Ocean Color Instrument (OCI) on board. For the first time reflected sunlight from the Earth across a broad spectrum from the ultraviolet (UV: 350 nm) to the short wave infrared (SWIR: 2260 nm) will be measured from a single instrument at 1 km spatial resolution. While seven discrete bands will represent the SWIR, the spectrum from 350 to 890 nm will be continuously covered with a spectral resolution of 5 nm. OCI will thus combine in a single instrument (and at an enhanced spatial resolution for the UV) the heritage capabilities of the Moderate resolution Imaging Spectroradiometer (MODIS) and the Ozone Monitoring Instrument (OMI), while covering the oxygen A-band (O2A). Designed for ocean color and ocean biology retrievals, OCI also enables continuation of heritage satellite aerosol products and the development of new aerosol characterization from space. In particular the combination of MODIS and OMI characteristics allows deriving aerosol height, absorption and optical depth along with a measure of particle size distribution. This is achieved by using the traditional MODIS visible-to-SWIR wavelengths to constrain spectral aerosol optical depth and particle size. Extrapolating this information to the UV channels allows retrieval of aerosol absorption and layer height. A more direct method to derive aerosol layer height makes use of O2A absorption methods, despite the relative coarseness of the nominal 5 nm spectral resolution of OCI. Altogether the PACE mission with OCI will be an unprecedented opportunity for aerosol characterization that will continue climate data records from the past decades and propel aerosol science forward toward new opportunities.

PACE↗
NASA NTRS2000

Sensitivity of Cirrus Bidirectional Reflectance at MODIS Bands to Vertical Inhomogeneity of Ice Crystal Habits and Size Distribution

A common assumption in satellite imager-based cirrus retrieval algorithms is that the radiative properties of a cirrus cloud may be represented by those associated with a specific ice crystal shape (or habit) and a single particle size distribution. However, observations of cirrus clouds have shown that the shapes and sizes of ice crystals may vary substantially with height within the clouds. In this study we investigate the sensitivity of the top-of-atmosphere bidirectional reflectances at two MODIS bands centered at 0.65 micron and 2.11 micron to the cirrus models assumed to be either a single homogeneous layer or three distinct but contiguous, layers. First, we define the single- and three-layer cirrus cloud models with respect to ice crystal habit and size distribution on the basis of in situ replicator data acquired during the First ISCCP Regional Experiment (FIRE-II), held in Kansas during the fall of 1991. Subsequently, fundamental light scattering and radiative transfer theory is employed to determine the single scattering and the bulk radiative properties of the cirrus cloud. Regarding the radiative transfer computations, we present a discrete form of the adding/doubling principle by introducing a direct transmission function, which is computationally straightforward and efficient an improvement over previous methods. For the 0.65 micron band, at which absorption by ice is negligible, there is little difference between the bidirectional reflectances calculated for the one- and three-layer cirrus models, suggesting that the vertical inhomogeneity effect is relatively unimportant. At the 2.11 micron band, the bidirectional reflectances computed for both optically thin (tau = 1) and thick (tau = 10) cirrus clouds show significant differences between the results for the one- and three-layer models. The reflectances computed for the three-layer cirrus model are substantially larger than those computed for the single-layer cirrus. Finally, we find that cloud reflectance is very sensitive to the optical properties of the small crystals that predominate in the top layer of the three-layer cirrus model. It is critical to define the most realistic geometric shape for the small "quasi-spherical" ice crystals in the top layer for obtaining reliable single-scattering parameters and bulk radiative properties of cirrus.

Yang, P.↗
NASA NTRS2005

Nonlinear Attitude Filtering Methods

This paper provides a survey of modern nonlinear filtering methods for attitude estimation. Early applications relied mostly on the extended Kalman filter for attitude estimation. Since these applications, several new approaches have been developed that have proven to be superior to the extended Kalman filter. Several of these approaches maintain the basic structure of the extended Kalman filter, but employ various modifications in order to provide better convergence or improve other performance characteristics. Examples of such approaches include: filter QUEST, extended QUEST, the super-iterated extended Kalman filter, the interlaced extended Kalman filter, and the second-order Kalman filter. Filters that propagate and update a discrete set of sigma points rather than using linearized equations for the mean and covariance are also reviewed. A two-step approach is discussed with a first-step state that linearizes the measurement model and an iterative second step to recover the desired attitude states. These approaches are all based on the Gaussian assumption that the probability density function is adequately specified by its mean and covariance. Other approaches that do not require this assumption are reviewed, including particle filters and a Bayesian filter based on a non-Gaussian, finite-parameter probability density function on SO(3). Finally, the predictive filter, nonlinear observers and adaptive approaches are shown. The strengths and weaknesses of the various approaches are discussed.

F Landis Markley↗
NASA NTRS2017

Kinetic Scale Structure of Low-Frequency Waves and Fluctuations

The dissipation of solar wind turbulence at kinetic scales is believed to be important for the heating of the corona and for accelerating the wind. The linear Vlasov kinetic theory is a useful tool for identifying various wave modes, including kinetic Alfven, fast magnetosonic/whistler, and ion-acoustic (or kinetic slow), and their possible roles in the dissipation. However, the kinetic mode structure in the vicinity of ion-cyclotron modes is not clearly understood. The present paper aims to further elucidate the structure of these low-frequency waves by introducing discrete particle effects through hybrid simulations and Klimontovich formalism of spontaneous emission theory. The theory and simulation of spontaneously emitted low-frequency fluctuations are employed to identify and distinguish the detailed mode structures associated with ion-Bernstein modes versus quasi-modes. The spontaneous emission theory and simulation also confirm the findings of the Vlasov theory in that the kinetic Alfven waves can be defined over a wide range of frequencies, including the proton cyclotron frequency and its harmonics, especially for high-beta plasmas. This implies that these low-frequency modes may play predominant roles even in the fully kinetic description of kinetic scale turbulence and dissipation despite the fact that cyclotron harmonic and Bernstein modes may also play important roles in wave-particle interactions.

waves↗
NASA NTRS2018

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗
NASA NTRS2015

Human Exploration of Phobos

This study developed, analyzed, and compared mission architectures for human exploration of Mars' Moons within the context of an Evolvable Mars Campaign. METHODS: All trades assumed conjunction class missions to Phobos (approximately 500 days in Mars system) as it was considered the driving case for the transportation architecture. All architectures assumed that the Mars Transit Habitat would remain in a High Mars Orbit with crewmembers transferring between HMO and Phobos in a small crew taxi vehicle. A reference science / exploration program was developed including performance of a standard set of tasks at 55 locations on the Phobos surface. Detailed EVA timelines were developed using realistic flight rules to accomplish the reference science tasks using exploration systems ranging from jetpacks to multi-person pressurized excursion vehicles combined with Phobos surface and orbital (L1, L4/L5, 20km Distant Retrograde Orbit) habitat options. Detailed models of propellant mass, crew time, science productivity, radiation exposure, systems and consumables masses, and other figures of merit were integrated to enable quantitative comparison of different architectural options. Options for pre-staging assets using solar electric propulsion (SEP) vs. delivering all systems with the crew were also evaluated. Seven discrete mission architectures were evaluated. RESULTS: The driving consideration for habitat location (Phobos surface vs. orbital) was radiation exposure, with an estimated reduction in cumulative mission radiation exposure of up to 34% (vs. Mars orbital mission) when the habitat is located on the Phobos surface, compared with only 3-6% reduction for a habitat in a 20km DRO. The exploration utility of lightweight unpressurized excursion vehicles was limited by the need to remain within 20 minutes of Solar Particle Event radiation protection combined with complex GN&C systems required by the non-intuitive and highly-variable gravitational environment. Two-person pressurized excursion vehicles as well as mobile surface habitats offer significant exploration capability and operational benefits compared with unpressurized EVA mobility systems at the cost of increased system and propellant mass. Mechanical surface translation modes (i.e. hopping) were modeled and offer potentially significant propellant savings and the possibility of extended exploration operations between crewed missions. Options for extending the utilization of the crew taxi vehicle were examined, including use as an exploration asset for Phobos surface exploration (when combined with an alternate mobility system) and as an EVA platform, both on Phobos and for contingency EVA on the Mars Transit Habitat. CONCLUSIONS: Human exploration of Phobos offers a scientifically meaningful first step towards human Mars surface missions that develops and validates transportation, habitation, and exploration systems and operations in advance of the Mars landing systems.

Abercromby, Andrew F. J.↗
NASA NTRS2004

Oh! I Slipped the Surly Bounds of Earth....and Ran into Space Weather!

Over the past decade the concept of space weather has been introduced and matured in both the scientific community and popular press. Likewise the concept of space climatology recently also is being advanced. Closely linked to these concepts are their impacts on ground- and space-based technological systems; one such system commonly mentioned is manned space flight exemplified by the Space Shuttle and International Space Station (ISS). From a manned space flight perspective, space weather and space climatology have significant effects on the amount of radiation exposure received by humans in space from the ambient high-energy charged particles present in interplanetary space and trapped in the geomagnetosphere. Whereas the impact of space weather for most technological systems is usually discrete and well correlated in time, the principle impact of space weather and space climatology is to increase the probability of latent cancer formation in thetraveler cohort. In this regard, while space weather may be the dominating factor for a given mission, over the life of a long-term program such as the Space Shuttle or ISS space climatology is the controlling factor of latent cancer risk. Human radiation exposure enhancements associated with space weather disturbances has been a concern among scientist and mission controllers since the inception of manned spaceflight nearly forty years ago. This led NASA to develop, in conjunction with the Environmental Science Services Administration s Space Disturbance Forecast Center and the USAF/AWS, the Solar Particle Alert Network (SPAN)-the foundation of an initial U.S. space weather monitoring and forecasting service. Since Apollo, routine space flight operations have evolved to include the use of space weather and climatology data provided through a world-wide network of operational space weather data services to predict and recommend actions to minimize astronaut radiation exposures. NASA Space Radiation Analysis Group (SRAG) flight controllers use real-time space weather data to detect and assess the impact of solar particle events, outer electron belt enhancements, the formation of pseudo-stable additional trapped radiation belts, and the solar cycle modulation of trapped radiation belts and galactic cosmic rays. Energetic particle data from GOES spacecraft are automatically ingested from NOAA Space Environment Center data servers and used to drive a model for the estimating the exposure to astronauts from solar particle events. While adequate for current manned space flight support, the existing operational space weather support system requires improvements to address the anticipated evolution in both the character of manned missions as well as space flight operations management. Necessary space weather data improvements include: reliably available (near) real-time space weather data on a fixed schedule via redundant access methods that support autonomous data acquisition by computer systems behind enterprise firewalls; and rapid transition of promising research sensors into operational systems.

Golightly, Michael J.↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
NASA NTRS2001

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗
NASA NTRS1973

On Lidar Sounding of the Atmosphere to Estimate Static and Dynamic Characteristics of Aerosol Inhomogeneities

1. A possible application of intensity fluctuations of a pulse light signal reflected by atmospheric aerosols is analyzed by the correlation method to evaluate static (medium sizes, shape) and dynamic (speed and direction of movement, lifetime) characteristics of aerosol inhomogeneities. The aerosol inhomogeneities are assumed to be expanded, pressed, disintegrated and originated constantly in accordance with random laws, the set of inhomogeneities as a whole traveling together with air masses and having predominant movement in wind direction. It is shown that the characteristics of aerosol inhomogeneities considered can be expressed by the coefficients of the correlation function expansion of the reflected signal fluctuation intensity in Tailor series. 2/ Correlation systems for evaluating static and dynamic characteristics of driving objects can be divided into two types according to the kind and quantity of used information: the systems with coordinates of the information removal "points" to be fixed in space, and the systems with a parallel simultaneous information removal at discrete moments of time. The systems for determination of wind direction considered in are the examples of the first type system. However, the operating information removal for two points is insufficient to estimate completely static and dynamic characteristics of inhomogenities, their quantity ought to be increased up to three of them for two-dimensional problem and up to four of them for three-dimensional problem as it is usually done in the ionospheric studies. The second type systems are used for the investigation of a medium shape and speed of the clouds according to photographs made from satellites. These systems are also used for solution of navigation problems. The use of optical quantum generators with a scanning beam is seen to increase greatly the working information removal in comparison with the first type systems. Nevertheless, scanning rate is not sufficient sometimes in order· to consider a general picture of aerosol inhomogeneities to be stationary. In this connection the use of the systems of second type treatment becomes a matter of essential difficulty. 3. Aerosol inhomogeneities simulation has been carried out on the basis of the digital computer experiments with the aim of estimating static and dynamic characteristics of inhomogeneities by an optical beam in the atmosphere at different scanning procedures. The dependence of determination accuracy of these characteristics on the type of chosen laws of aerosol particle distribution in the atmosphere, the parameters of inhomogeneities geometry, their speed and the law of scanning have been obtained.

Zuev, V. E.↗