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 613 records · Page 34

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↗

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↗

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↗

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.↗

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.↗

Mass-loading and parallel magnetized shocks

Recent observations at comets Giacobini-Zinner and Halley suggest that simple nonreacting gas dynamics or MHD is an inappropriate description for the bow shock. The thickness of the observed (sub)shock implies that mass-loading is an important dynamical process within the shock itself, thereby requiring that the Rankine-Hugoniot conditions possess source terms. This leads to shocks with properties similar to those of combustion shocks. The paper considers parallel magnetized shocks subjected to mass-loading, describes some properties which distinguish them from classical MHD parallel shocks, and establishes the existence of a new kind of MHD compound shock. These results will be of importance both to observations and numerical simulations of the comet-solar wind interaction.

Zank, G. P.↗

Shock capturing schemes for multidimensional flow

Progress made in the development of a genuinely multidimensional finite volume algorithm for problems in inviscid gas dynamics is presented. The approach entails: (1) the reconstruction of flowfield data using a planar wave pattern in which the strengths and orientations of the component waves are derived independently of the mesh geometry; and (2) the development of a flux formula which provides a numerical approximation to the flux at a finite volume cell face during the passage of waves which are in general oblique to the face. Several algorithms are outlined, and the most recent developments are included. The results of several numerical cases are also included.

Parpia, Ijaz H.↗

Composite load spectra for select space propulsion structural components

The objective of this program is to develop generic load models with multiple levels of progressive sophistication to simulate the composite load spectra that are induced in space propulsion system components, representative of Space Shuttle Main Engines (SSME), such as transfer ducts, turbine blades, and liquid oxygen (LOX) posts and system ducting. These models will be developed using two independent approaches. The first approach consists of using state-of-the-art probabilistic methods to describe the individual loading conditions and combinations of these loading conditions to synthesize the composite load spectra simulation. The methodology required to combine the various individual load simulation models (hot-gas dynamic, vibrations, instantaneous position, centrifugal field, etc.) into composite load spectra simulation models will be developed under this program. A computer code incorporating the various individual and composite load spectra models will be developed to construct the specific load model desired. The second approach, which is covered under the options portion of the contract, will consist of developing coupled models for composite load spectra simulation which combine the (deterministic) models for composite load dynamic, acoustic, high-pressure and high rotational speed, etc., load simulation using statistically varying coefficients. These coefficients will then be determined using advanced probabilistic simulation methods with and without strategically selected experimental data. This report covers the efforts of the third year of the contract. The overall program status is that the turbine blade loads have been completed and implemented. The transfer duct loads are defined and are being implemented. The thermal loads for all components are defined and coding is being developed. A dynamic pressure load model is under development. The parallel work on the probabilistic methodology is essentially completed. The overall effort is being integrated in an expert system code specifically developed for this project.

Newell, J. F.↗

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗

A general numerical model for wave rotor analysis

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Paxson, Daniel W.↗

Treating convection in sequential solvers

The treatment of the convection terms in the sequential solver, a standard procedure found in virtually all pressure based algorithms, to compute the flow problems with sharp gradients and source terms is investigated. Both scalar model problems and one-dimensional gas dynamics equations have been used to study the various issues involved. Different approaches including the use of nonlinear filtering techniques and adoption of TVD type schemes have been investigated. Special treatments of the source terms such as pressure gradients and heat release have also been devised, yielding insight and improved accuracy of the numerical procedure adopted.

Shyy, Wei↗

Swept shock/boundary layer interaction experiments in support of CFD code validation

Research on the topic of shock wave/turbulent boundary-layer interaction was carried out during the past three years at the Penn State Gas Dynamics Laboratory. This report describes the experimental research program which provides basic knowledge and establishes new data on heat transfer in swept shock wave/boundary-layer interactions. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angle ranging from 10 deg to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a k-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.

Settles, G. S.↗

The linear Boltzmann equation in slab geometry - Development and verification of a reliable and efficient solution

The linear Boltzmann equation can be cast in a form mathematically identical to the radiation-transport equation. A multigroup procedure is used to reduce the energy (or velocity) dependence of the transport equation to a series of one-speed problems. Each of these one-speed problems is equivalent to the monochromatic radiative-transfer problem, and existing software is used to solve this problem in slab geometry. The numerical code conserves particles in elastic collisions. Generic examples are provided to illustrate the applicability of this approach. Although this formalism can, in principle, be applied to a variety of test particle or linearized gas dynamics problems, it is particularly well-suited to study the thermalization of suprathermal particles interacting with a background medium when the thermal motion of the background cannot be ignored. Extensions of the formalism to include external forces and spherical geometry are also feasible.

Stamnes, K.↗

CW dye laser technique for simultaneous, spatially-resolved measurements of temperature, pressure, and velocity of NO in an underexpanded free jet

Gas dynamic quantities within an underexpanded free jet were measured nonintrusively using a rapid-tuning, CW ring dye laser. A nitrogen jet was seeded with 0.5 percent NO in N2, and the conditions were controlled such that a barrel shock formed. The frequency-doubled output of the dye laser was used to spectrally resolve rotational lines in the NO gamma band near 225 nm. With the rapid-tuning capability, these rotational spectra were acquired at a repetition rate of 4 kHz. Spatial resolution was afforded by monitoring the induced fluorescence via a lens and photomultiplier tube. Modeling the spectrally-resolved features with Voigt profiles permitted simultaneous measurements of NO velocity, rotational temperature, and pressure. Expansion of the jet was assumed to be isentropic, and agreement between measured and expected values was typically better than 10 percent over most of the Mach-number range encountered. At high Mach numbers, the measured rotational temperatures systematically departed from the isentropic temperature distribution. Such a measured departure could be ascribed to the onset of a non-Boltzmann distribution of NO rotational states.

Di Rosa, M. D.↗

Species and temperature measurement in H2/O2 rocket flow fields by means of Raman scattering diagnostics

Validation of CFD codes developed for prediction and evaluation of rocket performance is hampered by a lack of experimental data. Nonintrusive laser based diagnostics are needed to provide spatially and temporally resolved gas dynamic and fluid dynamic measurements. This paper reports the first nonintrusive temperature and species measurements in the plume of a 110 N gaseous hydrogen/oxygen thruster at and below ambient pressures, obtained with spontaneous Raman spectroscopy. Measurements at 10 mm downstream of the exit plane are compared with predictions from a numerical solution of the axisymmetric Navier-Stokes and species transport equations with chemical kinetics, which fully model the combustor-nozzle-plume flowfield. The experimentally determined oxygen number density at the centerline at 10 mm downstream of the exit plane is four times that predicted by the model. The experimental number density data fall between those numerically predicted for the exit and 10 mm downstream planes in both magnitude and radial gradient. The predicted temperature levels are within 10 to 15 percent of measured values.

De Groot, Wim A.↗

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 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 3D problems. The topical issues of reduced resolution procedures, stored lists, and variance reduction schemes are discussed in some detail.

Bird, G. A.↗

Thermomolecular effect on pressure measurements with orifices in transitional flow

Results are presented of a study carried out in order to verify a version of the method (presented at the 7th International Rarefied Gas Dynamics Symposium) for calculating the influence of thermomolecular flow on the pressures measured by means of orifices in walls exposed to rarefied gases. In these experiments, the 'orifice effect' is examined under conditions that are significantly different from those for which the semiempirical (SE) method of Kinslow and Potter (1971) was shown to be successful. Comparisons between the predicted results with flight data and the results obtained in the laboratory, combined with limited DSMC calculations, indicate that the SE method remains a useful tool which is relatively easy to apply and which gives good results when an appropriate energy accommodation coefficient is used.

Potter, J. L.↗

Calculating Shocks In Chemically Reacting Flows

Total-variation-diminishing dissipation operator results in crisp shocks. Report discusses method for numerical simulation of shocks in hypersonic, chemically reacting flows. Introduces mathematical model of gas of as many as three species for which species-conservation equations coupled with gas-dynamic equations.

Eberhardt, Scott↗