Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical simulations of the nonlinear kink modes in linearly stable supersonic slip surfaces

High-resolution numerical simulations of supersonic slip surfaces were performed using the Cray 2 at the Minnesota Supercomputer Center to confirm earlier analytic nonlinear stability calculations of such structures. The present study presents a numerical solution for the problem solved analytically by Artola and Majda (1987) and follows growing modes well into the large-amplitude regime. The response of an equal-density Mach-4 slip surface to a variety of impinging nonlinear sound wave trains is computed using the piecewise-parabolic method (PPM). The nature of the nonlinear kink modes observed in the simulations was found to be similar to that discussed by Artola and Majda. Most of the speeds moved in either direction, with speeds near the predicted value. It is suggested that the stationary modes will eventually dominate the flow at much later times, and it is argued that the stationary modes are more disruptive than the propagating kink modes.

Pedelty, Jeffrey A.↗
NASA NTRS1990

Numerical simulation of actively controlled space structures

The numerical simulation of large, actively controlled space structures is a problem that has received considerable attention. A new computational method is presented in this paper, which simulates the vibration suppression of a large space structure by a linear compensator. The structure is assumed to be modeled by a large linear finite element model. Nested circularly linked lists and the trapezoidal rule are key aspects of this algorithm. There is an emphasis on generality, while establishing an equivalence between numerical stability and system stability.

Quan, Ralph↗
NASA NTRS1991

A detailed numerical investigation of Burke-Schumann gaseous and spray flames

The classical Burke-Schumann gaseous diffusion flame, and the effect of introducing fuel in the form of liquid sprays are investigated numerically. The time-dependent two-dimensional, axisymmetric conservation equations are solved for the gas phase. A Lagrangian approach is used for the dispersed phase. The chemical reactions are modeled through the one-step global reaction scheme. The numerical results show that the computed flame structure is significantly different from that given by the analytical solution. The computed flame is about 15 percent longer and 5 percent narrower than the classical Burke-Schumann flame. The predicted flame shape is, however, in better agreement with experimental observations than the Burke-Schumann flame. The effects of introducing liquid sprays is found to produce thinner and longer flames than its counterpart gaseous flame. A detailed comparison of the structures of gaseous and spray flames is made, and the differences are discussed.

Mawid, M. A.↗
NASA NTRS1991

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Cambier, Jean-Luc↗
NASA NTRS1991

Hypersonic waverider analysis - A comparison of numerical and experimental results

In this study two waverider configurations are analyzed numerically and experimentally. The geometries are designed using an inverse design and optimization code with the optimization procedure driven by maximized L/D and limited by geometric constraints. Euler and Navier-Stokes numerical simulations are performed with a partially flux-split finite-difference code developed at NASA Ames, F3D. Experimental data is obtained at NASA Langley in the Unitary Plan Wind Tunnel (UPWT) and the 20-Inch Mach 6 wind tunnel for flow visualization, force/moment, and surface pressure measurements at both on- and off-design conditions. Euler simulations are compared with analytic results obtained from conical flow theory, and Navier-Stokes simulations are compared with experimental results both qualitatively and quantitatively for on- and off-design conditions.

Jones, K. D.↗
NASA NTRS1991

A numerical study of the flow establishment time in hypersonic shock tunnels

This paper presents a numerical study of unsteady hypersonic viscous flows in 2D hypersonic shock tunnels. The numerical experiments on the combined nozzle/model flows, utilizing the configurations mounted in the test section of a 14.7-long shock tunnel, which has a design Mach number of Me = 6 and a stagnation temperature of T(0) = 3063 R, have been performed by solving the 2D Navier-Stokes equations with an ADI-type TVD scheme. The flow establishment times around aerodynamic models, such as a flat plate, a biconvex airfoil, thin and thick double wedges, and a circular cylinder, are calculated and compared. The analysis shows that, for the most severely separated case of a circular cylinder, the flow becomes stabilized within about 8 msec including the starting process, and its nondimensional flow establishment time is 46. The flow establishment time increases with increasing length of the separated region and decreases with increasing freestream velocity.

Lee, Jang-Yeon↗
NASA NTRS1991

Numerical analysis of rarefied slit flows. II - Navier-Stokes simulations

The model problem of pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is analyzed via full Navier-Stokes numerical simulation. Parametric solutions are generated for slit-height based Knudsen number ranging from continuum to transitional flow and for reservoir pressure ratios leading to subsonic and supersonic flow. The change in the structure of the flowfield near the slit as a function of pressure ratio and Knudsen number are quantified from a purely continuum standpoint. The choice of numerical domain size, boundary conditions and treatment of the slit are also discussed. As expected, comparison with a Direct Simulation Monte Carlo solution for a highly rarefied case shows large differences in the predicted mass flow. The cause of these differences can be quantified through detailed comparison of the local flowfield properties. For the larger pressure ratio cases qualitative trends with increasing rarefaction are discussed, including the change in the sonic line shape in the slit and the in total mass flow.

Wadsworth, D. C.↗
NASA NTRS1991

Numerical modeling of solidification and convection in a viscous pure binary eutectic system

The solidification and convection of the pure binary eutectic silicate system diopside-anorthite (Di-An) is numerically modeled. A mass-weighted enthalpy of fusion is used to account for the second solid phase (An) which crystallizes at the solidus temperature. Variable under-relaxation is used to aid convergence of the momentum equations in the implementation of the SIMPLER algorithm used to solve the two-dimensional continuum conservation equations. Numerical experiments of the solidification of Di80 melt show that a large temperature drop occurs across the solid and mush regions which decreases convective vigor in the liquid. Interesting compositional segregation patterns are produced during the solidification of Di80.

Oldenburg, Curtis M.↗
NASA NTRS1991

Numerical simulation of shrouded propellers

A numerical model was developed for the evaluation of the performance characteristics of a shrouded propeller. Using this model, a computational study was carried out to investigate the feasibility of improving the aerodynamic performance of a propeller by encasing it in a shroud. The propeller blade was modeled by a segmented bound vortex positioned along the span of the blade at its quarter-chord-line. The shroud was modeled by a number of discrete vortex rings. Due to the mutual dependence of shroud and propeller vortex strengths and the propeller vortex wake an iterative scheme was employed. Three shroud configurations were considered: a cylindrical and two conical shrouds. The computed performance of the shrouded propeller was compared with that of a free propeller of identical propeller geometry. The numerical results indicated that the cylindrical shroud outperformed the conical shroud configurations for the cases considered. Furthermore, when compared to the free propeller performance, the cylindrical shroud showed a considerable performance enhancement over the free propeller. However, the improvements were found to decrease with an increase in the advance ratio and to virtually diminish at advance ratios of about 2.5.

Afjeh, Abdollah A.↗
NASA NTRS1991

Numerical and literal aeroelastic-vehicle-model reduction for feedback control synthesis

The simplification of a high-order, literal model for large flexible aircraft is discussed. Areas of model fidelity that are critical if the model is to be used for control law synthesis are presented. Several simplification techniques, some new and some widely available, that can deliver the necessary model fidelity are presented and applied to a model from the literature. The techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory, is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle vibrational characteristics and rigid-body and aeroelastic stability derivatives, thus providing insight into the underlying causes for critical dynamic characteristics.

Newman, Brett↗
NASA NTRS1990

Numerical prediction of turbulent flame stability in premixed/prevaporized (HSCT) combustors

A numerical analysis of combustion instabilities that induce flashback in a lean, premixed, prevaporized dump combustor is performed. KIVA-II, a finite volume CFD code for the modeling of transient, multidimensional, chemically reactive flows, serves as the principal analytical tool. The experiment of Proctor and T'ien is used as a reference for developing the computational model. An experimentally derived combustion instability mechanism is presented on the basis of the observations of Proctor and T'ien and other investigators of instabilities in low speed (M less than 0.1) dump combustors. The analysis comprises two independent procedures that begin from a calculated stable flame: The first is a linear increase of the equivalence ratio and the second is the linear decrease of the inflow velocity. The objective is to observe changes in the aerothermochemical features of the flow field prior to flashback. It was found that only the linear increase of the equivalence ratio elicits a calculated flashback result. Though this result did not exhibit large scale coherent vortices in the turbulent shear layer coincident with a flame flickering mode as was observed experimentally, there were interesting acoustic effects which were resolved quite well in the calculation. A discussion of the k-e turbulence model used by KIVA-II is prompted by the absence of combustion instabilities in the model as the inflow velocity is linearly decreased. Finally, recommendations are made for further numerical analysis that may improve correlation with experimentally observed combustion instabilities.

Winowich, Nicholas S.↗
NASA NTRS1991

Practical design considerations and performance characteristics of high numerical aperture holographic lenses

The diffraction efficiency of interferometrically formed holographic lenses is influenced by the recording geometry and properties of the recording material. Variations in efficiency increase when attempting to make high numerical aperture elements. The factors which influence the diffraction efficiency of high numerical aperture holographic lenses are examined.

Kostuk, Raymond K.↗
NASA NTRS1991

A critical comparison of second order closures with direct numerical simulation of homogeneous turbulence

Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗
NASA NTRS1992

Gamma ray burst astrometry 2: Numerical tests

Since the announcement of the discovery of sources of gamma ray radiation in 1973, many more reports of such bursts have been published. Numerous artificial satellites have been equipped with gamma ray detectors including GRO. Unfortunately, almost no progress has been made in identifying the sources of this high energy radiation. Only one visible counterpart is known. It is suspected that this is a consequence of the methods currently used to define gamma ray burst source 'error boxes'. An alternative procedure was proposed in 1988 by Taff. Herein, Monte Carlo simulations are reported of the efficacy of this technique using realistic burst timing uncertainties and satellite location errors as well as a variety of satellite constellations. Since these are controlled numerical experiments, the dependence is examined of the statistics of the errors in the deduced burst wavefront normal as a function of the timing inconsistencies, detector location standard deviations, and especially the number and distribution of the detectors. The results clearly show that an arc minute prediction of a unique burst location is routinely obtainable once there are at least two interplanetary detectors.

Taff, L. G.↗
NASA NTRS1992

Numerical studies of the fluid and optical fields associated with complex cavity flows

Numerical solutions for the flowfield about several cavity configurations have been computed using the Reynolds averaged Navier-Stokes equations. Comparisons between numerical and experimental results are made in two dimensions for free shear layers and a rectangular cavity, and in three dimensions for the transonic aero-window problem of the Stratospheric Observatory for Infrared Astronomy (SOFIA). Results show that dominant acoustic frequencies and magnitudes of the self excited resonant cavity flows compare well with the experiment. In addition, solution sensitivity to artificial dissipation and grid resolution levels are determined. Optical path distortion due to the flow field is modelled geometrically and is found to match the experiment. The fluid field was computed using a diagonalized scheme within an overset mesh framework. An existing code, OVERFLOW, was utilized with the additions of characteristic boundary condition and output routines required for reduction of the unsteady data. The newly developed code is directly applicable to a generalized three dimensional structured grid zone. Details are provided in a paper included in Appendix A.

Atwood, Christopher A.↗
NASA NTRS1992

Study of turbulence structure through numerical simulations: The perspective of drag and reduction

The first part of this note concerns organized structures in the near-wall region of a turbulent boundary layer. A particular emphasis is given to those organized motions believed to be responsible for the turbulence production process in the wall region, and hence, indirectly responsible for the viscous drag. Examples are selected to illustrate how the analyses of numerically generated databases have contributed to improve our understanding of the organized structures of turbulent boundary layers. In the second part, results from an exploratory study based on the direct numerical simulation of the concepts for active control of turbulent boundary layers are presented. A significant drag reduction is achieved when the surface boundary condition is modified such that it could suppress the organized motion in the wall region. This drag reduction is accompanied by a significant reduction in the intensity of the organized structures and in the magnitude of Reynolds stresses throughout the flow.

Kim, J. J.↗
NASA NTRS1992

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗
NASA NTRS1990

A numerical method for prediction of compressible turbulent flows with closure models

A new computer code to solve the time averaged Navier-Stokes equations is developed. Many of the state-of-the-art numerical techniques and algorithms have been tested and implemented in the program in order to achieve a better numerical accuracy and code efficiency. Various turbulence models are tested for a wide range of flows. The initial focus has been on two-equation eddy-viscosity models, which are the most advanced available in current compressible flow codes. The long term goal will be to test Reynolds-Stress models and to explore their performance in the high Mach number range. Although testing and improvement of turbulence models for supersonic and hypersonic flows is the primary objective of this research, part of the effort has been devoted to analyzing the vortex breakdown phenomena using new computer programs. Some preliminary results on the breakdown of a vortex flow in a tube are reported. Present calculations are restricted to two dimensional flow geometry.

Huang, P. G.↗