Search NASA⌕ Search

SEARCH · Search NASA

Results for “gas dynamics”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A model of the energetic ion environment of Mars

Because Mars has a weak intrinsic magnetic field and a substantial atmosphere, instruments on orbiting spacecraft should detect a population of energetic heavy planetary ions which result from comet-like ion pickup in the solar wind and magnetosheath convection electric fields, in addition to those that might result from processes internal to a Martian 'magnetosphere.' Although this ion exosphere has been previously discussed in the literature, detailed predictions that might be directly applied to the interpretation of data are not available. Here a test particle model is used to construct a global picture of Martian pickup ions in the Mars environment. The model makes use of the recent Nagy and Cravens (1988) model of the Martian exosphere and Spreiter and Stahara's (1980) gas dynamic model of the magnetosheath. The pickup of ions originating at Phobos is also considered. Notable properties of the resulting ion distributions include their near-monoenergetic spectra, pancake pitch angle distributions, and large gyroradii compared to the planetary scale.

Luhmann, J. G.↗
NASA NTRS1990

A collision-selection rule for a particle simulation method suited to vector computers

A theory is developed for a selection rule governing collisions in a particle simulation of rarefied gas-dynamic flows. The selection rule leads to an algorithmic form highly compatible with fine grain parallel decomposition, allowing for efficient utilization of supercomputers having vector or massively parallel single instruction multiple data architectures. A comparison of shock-wave profiles obtained using both the selection rule and Bird's direct simulation Monte Carlo (DSMC) method show excellent agreement. The equation on which the selection rule is based is shown to be directly related to the time-counter procedure in the DSMC method. The results of several example simulations of representative rarefied flows are presented, for which the number of particles used ranged from 10 to the 6th to 10 to the 7th demonstrating the greatly improved computational efficiency of the method.

Baganoff, D.↗
NASA NTRS1990

Control-oriented modeling of combustion and flow processes in liquid propellant rocket engines

This paper presents a control-oriented model of the flow, reaction, and transport processes in liquid propellant rocket combustion chambers, based on the multicomponent conservation laws of gas dynamics. This model provides a framework for the inclusion of detailed chemical kinetic relations, viscous and other dissipative effects, a variety of actuators and sensors, as well as process and measurement disturbances. In addition to its potential usefulness to the designer in understanding the dynamical complexity of the system and the sources of model uncertainty, the model provides a rigorous basis for control system design. An appraisal of current and feasible actuators and sensors, and their mathematical representation are included.

Bentsman, Joseph↗
NASA NTRS1990

Assessment of shock capturing schemes for resonant flows in nonlinear instability analysis

The paper presents computational assessment of advanced numerical schemes for nonlinear acoustic problems related to combustion instabilities in liquid rocket engines. Several time-accurate, shock capturing schemes have been evaluated on a benchmark, closed-end resonant pipe flow problem. It involves the numerical solution of inviscid, compressible gas dynamics equations to predict acoustic wave propagation, wave steepening, formation of shocks, acoustic energy dissipation and wave-wall reflection for several hundred wave cycles. It was demonstrated that high accuracy TVD type schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time. The selected scheme was then applied to analyze the acoustic responses of resonant pipe-resonator, radial acoustic modes and hub-baffle configurations. Interesting observations of wave shape and damping characteristics have been drawn from presented computational studies.

Przekwas, A. J.↗
NASA NTRS1990

Efficiency and discrepancy in the direct simulation methods

The paper reviews the factors that affect the versatility, computational efficiency, and accuracy of the various direct simulation methods that are being applied to problems in rarefied gas dynamics. The quality of the simulation can be affected by the underlying motivation and objectives. The influence of computer hardware developments on the cost-effectiveness is discussed, as are the software factors that influence the speed with which simulation programs can be applied to new problems. Particular attention is given to a novel method that improves the efficiency of the simulation for complex three-dimensional problems. The topical issues of reduced resolution procedures, stored lists, and variance reduction schemes are discussed in some detail.

Bird, G. A.↗
NASA NTRS1990

Least-squares finite element methods for compressible Euler equations

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L-sq-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L-sq method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H1-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.

Jiang, Bo-Nan↗
NASA NTRS1991

High-Resolution Numerical Simulation Of Shock Waves

NASA technical memorandum compares results of upwind and symmetric shock-capturing methods in numerical simulation of gas-dynamic flows. Both methods find shocks as sharp variations in fluid properties over few grid points. Methods differ in type of artificial viscosity introduced to stabilize computations. Symmetric method shown as accurate as upwind method, but with fewer and simpler time steps in transient case or fewer iteration steps in steady-state case.

Yee, H. C.↗
NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗
NASA NTRS1990

Gravitational star formation thresholds and gas density in three galaxies

It has long been held that the star formation rate (SFR) may be described as a power law of the gas density, p(exp n), as given by Schmidt (1959). However, this relation has as yet remained poorly defined and is likewise poorly understood. In particular, most studies have been investigations of global gas and star formation properties of galaxies, due to lack of adequate high-resolution data for detailed studies of individual galaxies. The three spiral galaxies in this study have published maps of both H2 (as traced by CO), and HI, thereby enabling the authors to investigate the relationship between total gas surface density and SFR. The purpose of the present investigation is the comparison of spatially-resolved total surface gas density in three galaxies (NGC 6946, M51, and M83) to sigma sub c as given by the above model. CO, HI and H alpha data for NGC 6946 were taken from Tacconi-Garman (1988), and for M51 and M83 from Lord (1987). The authors used a CO-H2 conversion of N(H2)/I sub CO(exp cos i = 2.8 x 10(exp 20) atoms cm(-2)/(K kms(-1), and summed the H2 and HI data for each galaxy to obtain the total hydrogen gas density. This total was then multiplied by a factor of 1.36 to include the contribution of helium to the total surface gas density. The authors assumed distances to NGC 6946, M51, and M83 to be 6.0, 9.6, and 8.9 Mpc respectively, with inclination angles of 30, 20, and 26 degrees. H alpha flux was used as the measure of SFR for NGC 6946, and SFR for the remaining two galaxies was taken directly from Lord as computed from H alpha measurements. The results of these full-disk studies thus show a remarkable correlation between the total gas density and the threshold densities given by the gravitational stability criterion. In particular, the threshold density appears to mark a lower boundary to the range of gas densities in these galaxies, which may have consequence in determining appropriate models for star formation and gas dynamics. More evidence is required to verify this result, and the authors are currently undertaking a high-resolution study of the nearby spiral M33 and other galaxies to further investigate this problem.

Oey, M. S.↗
NASA NTRS1990

An analysis of CCD images of the coma of Comet Halley

The analysis of selected CCD images of the coma of comet P/Halley is presented. The images were taken using specially designed filters that isolate regions of a comet's spectrum such that only sunlight which has been scattered by the dust in the coma is recorded. The modeling analysis objective is to make use of the skills developed in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae and to use those models as a basis for a new dust coma model. This model will include a self-consistant picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. The model is intended to be used as a tool to analyze selected images from the two sets of CCD images with the hope that it will help the understanding of the effects of a number of important processes on the spatial morphology of the observed dust coma. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observe scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; and (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gasdynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology on the outer dust coma.

Combi, Michael↗
NASA NTRS1991

Distributed Minimal Residual (DMR) method for acceleration of iterative algorithms

A new method for enhancing the convergence rate of iterative algorithms for the numerical integration of systems of partial differential equations was developed. It is termed the Distributed Minimal Residual (DMR) method and it is based on general Krylov subspace methods. The DMR method differs from the Krylov subspace methods by the fact that the iterative acceleration factors are different from equation to equation in the system. At the same time, the DMR method can be viewed as an incomplete Newton iteration method. The DMR method was applied to Euler equations of gas dynamics and incompressible Navier-Stokes equations. All numerical test cases were obtained using either explicit four stage Runge-Kutta or Euler implicit time integration. The formulation for the DMR method is general in nature and can be applied to explicit and implicit iterative algorithms for arbitrary systems of partial differential equations.

Lee, Seungsoo↗
NASA NTRS1991

Uniform high order spectral methods for one and two dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗
NASA NTRS1991

Detonation duct gas generator demonstration program

The feasibility of the generation of detonation waves moving periodically across high speed channel flow is experimentally demonstrated. Such waves are essential to the concept of compressing requirements and increasing the engine pressure compressor with the objective of reducing conventional compressor requirements and increasing the engine thermodynamic efficiency through isochoric energy addition. By generating transient transverse waves, rather than standing waves, shock wave losses are reduced by an order of magnitude. The ultimate objective is to use such detonation ducts downstream of a low pressure gas turbine compressor to produce a high overall pressure ratio thermodynamic cycle. A 4 foot long, 1 inch x 12 inch cross-section, detonation duct was operated in a blow-down mode using compressed air reservoirs. Liquid or vapor propane was injected through injectors or solenoid valves located in the plenum or the duct itself. Detonation waves were generated when the mixture was ignited by a row of spark plugs in the duct wall. Problems with fuel injection and mixing limited the air speeds to about Mach 0.5, frequencies to below 10 Hz, and measured pressure ratios of about 5 to 6. The feasibility of the gas dynamic compression was demonstrated and the critical problem areas were identified.

Wortman, Andrew↗
NASA NTRS1991

Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗
NASA NTRS1991

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially nonoscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1991

Theoretical description of spectral line profiles of parent molecules in cometary comae

The present overview of cometary spectral-line profiles obtainable through advancements in high-resolution spectroscopic studies, which allow the retrieval of coma kinematic properties from velocity-resolved spectral-line profiles, incorporates the most important gas dynamic processes into an outflow model which is tailored to the interpretation of spectroscopic observations of parent molecules. The model is then used to study the influence on parent-molecule spectral line profile formation of the field-of-view, the expansion velocity, the kinetic temperature, and the anisotropic outflow distributions.

Hu, Hong-Yao↗
NASA NTRS1991

A fully coupled flow simulation around spacecraft in low earth orbit

The primary objective of this investigation is to provide a full flow simulation of a spacecraft in low earth orbit (LEO). Due to the nature of the environment, the simulation includes the highly coupled effects of neutral particle flow, free stream plasma flow, nonequilibrium gas dynamics effects, spacecraft charging and electromagnetic field effects. Emphasis is placed on the near wake phenomenon and will be verified in space by the Wake Shield Facility (WSF) and developed for application to Space Station conditions as well as for other spacecraft. The WSF is a metallic disk-type structure that will provide a controlled space platform for highly accurate measurements. Preliminary results are presented for a full flow around a metallic disk.

Justiz, C. R.↗
NASA NTRS1991

The solar wind interaction with Mars - Mariner 4, Mars 2, Mars 3, Mars 5, and Phobos 2 observations of bow shock position and shape

An aggregate Mars bow shock data set using Mariner 4, Mars 2, Mars 3, Mars 5, and Phobos 2 observations has been analyzed. The results support the earlier conclusion that the mean distance to the subsolar shock at Mars is nearly 1.5 planetary radii, from which gas dynamic models predict an obstacle altitude of 500 km. The Martian bow shock does not appear to vary significantly in shape or altitude with the phase of the solar cycle. The unusually distant dayside bow shock crossings reported by Mars 2 and 3 also appear in the Phobos 3 observations, suggesting that the dayside obstacle can on rare occasions reach altitudes over 1000 km. The Martian bow shock differs from that of Venus in that its mean altitude is greater, it lacks a strong solar cycle variation, and its location is far more variable, including the occurrence of strong bow shocks over the dayside hemisphere at distances at least as great as the orbit of Phobos 2, i.e., 2.8 Mars radii.

Slavin, J. A.↗