SEARCH · Engineering Papers
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.
Presented here is a comprehensive review of the following aerodynamics elements: computational methods and applications, computational fluid dynamics (CFD) validation, transition and turbulence physics, numerical aerodynamic simulation, drag reduction, test techniques and instrumentation, configuration aerodynamics, aeroacoustics, aerothermodynamics, hypersonics, subsonic transport/commuter aviation, fighter/attack aircraft and rotorcraft.
A description is given of the hardware and software tools and techniques in use at NASA's Numerical Aerodynamic Simulation Facility for visualization of computational fluid dynamics. The hardware consists of high-performance graphics workstations connected to the supercomputer with high-bandwidth lines, a frame buffer connected to the supercomputer with UltraNet, a digital video recording system, and film recorders. The software permits the scientist to view the three-dimensional scenes dynamically, to zoom into a region of interest, and to rotate his viewing position to study any region of interest in more detail. The software also provides automated animation and video recording of the scenes. The digital effects unit on the video system facilitates comparison of computer simulations with flight or wind tunnel experiments.
An unsymmetric block Lanczos algorithm has been employed for the dynamic analysis of a large system which has arbitrary damping and/or repeated (or closely spaced) eigenvalues. In the algorithm development, the right and left Lanczos vectors are all theoretically biorthogonal to each other. However, these vectors may lose the biorthogonality owing to cancellation and roundoff errors. For the unsymmetric case there can be a breakdown, even without numerical errors. This paper describes computational techniques which have led to a robust unsymmetric block Lanczos algorithm.
A differential second-moment closure is applied to calculate nonswirling turbulent particle-laden jet flow. The closure is effected by extending constant-density single-phase turbulence models to the continuous (carrier) phase of two-phase flow. The dispersed (discrete) phase is treated by stochastic-Lagrangian technique. To reduce the effects of numerical (false) diffusion on the predicted results, a higher order differencing scheme, namely the flux-spline, is employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The experiment consists of a confined coaxial jet flow of air with glass beads. The beads' diameter is 105 microns and is injected from the inner jet with a mass loading ratio of 0.9. The computed results are in good agreement with experimental data, especially for the carrier phase. However, the predicted turbulence fluctuations associated with the discrete phase are lower than the data at the downstream region of the flow field.
Two recent theories of nonlinear mechanical response are quantitatively compared and related to experimental data. Computer techniques are formulated to handle the numerical integration and iterative procedures needed to solve the associated sets of coupled nonlinear differential equations. Problems encountered during these formulations are discussed and some open questions described. Bearing in mind these cautions, the consequences of changing parameters that appear in the formulations on the resulting engineering properties are discussed. Hence, engineering approaches to the analysis of thermoplastic matrix material can be suggested.
Phase space structure and chaos near a neutral line are studied with numerical surface-of-section (SOS) techniques and analytic methods. Results are presented for a linear neutral line model with zero crosstail electric field. It was found that particle motion can be divided into three regimes dependening on the value of the conserved canonical momentum, Py, and the conserved Hamiltonian, h. The phase space structure, using Poincare SOS plots, is highly sensitive to bn = Bn/B0 variations, but not to h variations. It is verified that the slow motion preserves the action, Jz, as evaluated by Sonnerup (1971), when the period of the fast motion is smaller than the time scale of the slow motion. Results show that the phase space structure and particle chaos depend sensitively upon Py and bn, but are independent of h.
Measurements were taken of the path of a simulated typical tissue segment or 'particle' within a rotating fluid as a function of gravitational strength, fluid rotation rate, particle sedimentation rate, and particle initial position. Parameters were examined within the useful range for tissue culture in the NASA rotating wall culture vessels. The particle moves along a nearly circular path through the fluid (as observed from the rotating reference frame of the fluid) at the same speed as its linear terminal sedimentation speed for the external gravitational field. This gravitationally induced motion causes an increasing deviation of the particle from its original position within the fluid for a decreased rotational rate, for a more rapidly sedimenting particle, and for an increased gravitational strength. Under low gravity conditions (less than 0.1 G), the particle's motion through the fluid and its deviation from its original position become negligible. Under unit gravity conditions, large distortions (greater than 0.25 inch) occur even for particles of slow sedimentation rate (less than 1.0 cm/sec). The particle's motion is nearly independent of the particle's initial position. Comparison with mathematically predicted particle paths show that a significant error in the mathematically predicted path occurs for large particle deviations. This results from a geometric approximation and numerically accumulating error in the mathematical technique.
The application of a two-step hybrid technique to the geometrically nonlinear analysis of curved beams is used to demonstrate the potential of hybrid analytical techniques in nonlinear structural mechanics. The hybrid technique is based on successive use of the perturbation method and a classical direct variational procedure. The functions associated with the various-order terms in the perturbation expansion of the fundamental unknowns, and their sensitivity derivatives with respect to material and geometric parameters of the beam, are first obtained by using the perturbation method. These functions are selected as coordinate functions (or modes) and the classical direct variational technique is then used to compute their amplitudes. The potential of the proposed hybrid technique for nonlinear analysis of structures is discussed. The effectiveness of the hybrid technique is demonstrated by means of numerical examples. The symbolic computation system Mathematica is used in the present study. The tasks performed on Mathematica include: (1) generation of algebraic expressions for the perturbation functions of the different response quantities and their sensitivity derivatives: and (2) determination of the radius of convergence of the perturbation series.
We examine the characteristics of the classical radiation emitted by a relativistic electron beam that propagates perpendicularly through a large amplitude relativistic plasma wave. Such a study is useful for evaluating the feasibility of using relativistic plasma waves as extremely short wavelength undulators for generating short wavelength radiation. The electron trajectories in a plasma wave undulator are obtained using perturbation techniques and are then compared to numerical simulation results. The frequency spectrum and angular distribution of the spontaneous radiation emitted by a single electron and the stimulated radiation gain are obtained analytically, and are then compared to 3-D numerical simulations. The characteristics of the plasma wave undulator are compared to the ac FEL undulator and the conventional FEL.
The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.
In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.
The effectiveness of a four-dimensional variational data assimilation (4-DVAR) technique for creating "bogus" vortices in numerical simulations of hurricanes is evaluated in this study. A series of numerical experiments is conducted to generate initial vortices for Hurricane Georges and Bonnie (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and wind information into a mesoscale numerical model (MM5). Several different strategies are tested for investigating the sensitivity of the initial vortex representation to the type of bogus information.
We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.
The main objective of this study is to validate the jet noise reduction potential of a concept associated with distributed exhaust nozzles. Under this concept the propulsive thrust is generated by a larger number of discrete plumes issuing from an array of small or mini-nozzles. The potential of noise reduction of this concept stems from the fact that a large number of small jets will produce very high frequency noise and also, if spaced suitably, they will coalesce at a smaller velocity to produce low amplitude, low frequency noise. This is accomplished through detailed acoustic and fluid measurements along with a Computational Fluidic Dynamic (CFD) solution of the mean (DE) Distributed Exhaust nozzle flowfield performed by Northrop-Grumman. The acoustic performance is quantified in an anechoic chamber. Farfield acoustic data is acquired for a DE nozzle as well as a round nozzle of the same area. Both these types of nozzles are assessed numerically using Computational Fluid Dynamic (CFD) techniques. The CFD analysis ensures that both nozzles issued the same amount of airflow for a given nozzle pressure ratio. Data at a variety of nozzle pressure ratios are acquired at a range of polar and azimuthal angles. Flow visualization of the DE nozzle is used to assess the fluid dynamics of the small jet interactions. Results show that at high subsonic jet velocities, the DE nozzle shifts its frequency of peak amplitude to a higher frequency relative to a round nozzle of equivalent area (from a S(sub tD) = 0.24 to 1. 3). Furthermore, the DE nozzle shows reduced sound pressure levels (as much as 4 - 8 dB) in the low frequency part of the spectrum (less than S(sub tD) = 0.24 ) compared to the round nozzle. At supersonic jet velocities, the DE nozzle does not exhibit the jet screech and the shock-associated broadband noise is reduced by as much as 12 dB.
The advent of inexpensive, high performance computer and new efficient algorithms have made possible the automatic recognition of numerically computed constants. In other words, techniques now exist for determining, within certain limits, whether a computed real or complex number can be written as a simple expression involving the classical constants of mathematics. In this presentation, some of the recently discovered techniques for constant recognition, notably integer relation detection algorithms, will be presented. As an application of these methods, the author's recent work in recognizing "Euler sums" will be described in some detail.
The Stardust Interstellar Dust Collector (SIDC) was intended to capture and return contemporary interstellar dust. The approx.0.1 sq m collector was composed of aerogel tiles (85% of the collecting area) and aluminum foils and was exposed to the interstellar dust stream for a total exposure factor of 20 sq m day. The Stardust Interstellar Preliminary Examination (ISPE) is a consortium-based project to characterize the collection using nondestructive techniques. Sandford et al. recently assessed numerous potential sources of organic contaminants in the Stardust cometary collector. These contaminants could greatly complicate the analysis and interperetation of any organics associated with interstellar dust, particularly because signals from these particles are expected to be exceedingly small. Here, we present a summary of FTIR analyses of over 20 aerogel keystones, many of which contained candidates for interstellar dust.
Clouds are highly effective in obscuring optical images of the Space Shuttle taken during its ascent by ground-based and airborne tracking cameras. Because the imagery is used for quick-look and post-flight engineering analysis, the Columbia Accident Investigation Board (CAIB) recommended the return-to-flight effort include an upgrade of the imaging system to enable it to obtain at least three useful views of the Shuttle from lift-off to at least solid rocket booster (SRB) separation (NASA 2003). The lifetimes of individual cloud elements capable of obscuring optical views of the Shuttle are typically 20 minutes or less. Therefore, accurately observing and forecasting cloud obscuration over an extended network of cameras poses an unprecedented challenge for the current state of observational and modeling techniques. In addition, even the best numerical simulations based on real observations will never reach "truth." In order to quantify the risk that clouds would obscure optical imagery of the Shuttle, a 3D model to calculate probabilistic risk was developed. The model was used to estimate the ability of a network of optical imaging cameras to obtain at least N simultaneous views of the Shuttle from lift-off to SRB separation in the presence of an idealized, randomized cloud field.
A simulation framework based on the Memory-Mapped-Files technique was created to operate multiple numerical processes in locked time-steps and send I/O data synchronously across to one-another to simulate system-dynamics. This simulation scheme is currently used to study the complex interactions between inlet flow-dynamics, variable-geometry actuation mechanisms, and flow-controls in the transition from the supersonic to hypersonic conditions and vice-versa. A study of Mode-Transition Control for a high-speed inlet wind-tunnel model with this MMF-based framework is presented to illustrate this scheme and demonstrate its usefulness in simulating supersonic and hypersonic inlet dynamics and controls or other types of complex systems.