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.
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.
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.
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.
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.
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.
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.
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.
Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.
The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.
A three dimensional General Circulation Model (GCM) of Titan's atmosphere is described. Initial results obtained with an economical two dimensional (2D) axisymmetric version of the model presented a strong superrotation in the upper stratosphere. Because of this result, a more general numerical study of superrotation was started with a somewhat different version of the GCM. It appears that for a slowly rotating planet which strongly absorbs solar radiation, circulation is dominated by global equator to pole Hadley circulation and strong superrotation. The theoretical study of this superrotation is discussed. It is also shown that 2D simulations systemically lead to instabilities which make 2D models poorly adapted to numerical simulation of Titan's (or Venus) atmosphere.
The present capabilities offered by the combination of advanced special and general purpose numerical codes and supercomputers are discussed. Emphasis is placed on unique demands placed by the space environment on the computation of realistic space experiments and recommendations for improving the scientific return on the space investment. It is concluded that numerical simulations can provide timely and relevant answers which can guide and improve future space experiments.
The results of a numerical investigation to predict the flow in the wake regions of compressor cascades, and wakes and mixing in rotors are presented. Attention is given to the flow in compressor cascades including the effects of change in loading (incidence) and the inlet freestream turbulence intensity. The numerical analysis shows a slight increase in the total pressure loss coefficient through the cascade with increasing turbulence levels.
The results of a numerical investigation to predict the flow field including wakes and mixing in axial-flow compressor rotors are presented. The wake behavior in a moderately loaded compressor rotor is studied numerically using a 3D incompressible Navier-Stokes solver with a high Reynolds number form of a turbulence model. The equations are solved using a time dependent implicit technique. The agreement between the measured data and the predictions is good; including the blade boundary-layer profiles, wake mean-velocity profiles, and decay. The ability of the pseudocompressibility scheme to predict the entire flow field including the near and far wake profiles and its decay characteristics, effect of loading, and the viscous losses of a 3D rotor flow field are demonstrated. The mixing in the downstream regions away from the hub and annulus walls is dominated by wake diffusion. In regions away from the walls the radial mixing is predominantly caused by the transport of mass, momentum, and energy by the radial component of velocity in the wake.
The advancement of high-speed axial-flow multistage compressors is impeded by a lack of detailed flow-field information. Recent development in compressor flow modeling and numerical simulation have the potential to provide needed information in a timely manner. The development of a computer program is described to solve the viscous form of the average-passage equation system for multistage turbomachinery. Programming issues such as in-core versus out-of-core data storage and CPU utilization (parallelization, vectorization, and chaining) are addressed. Code performance is evaluated through the simulation of the first four stages of a five-stage, high-speed, axial-flow compressor. The second part addresses the flow physics which can be obtained from the numerical simulation. In particular, an examination of the endwall flow structure is made, and its impact on blockage distribution assessed.
Radiative MHD equations are used for an optically thin plasma to carry out a numerical experiment related to the formation of 'postflare' loops. The numerical experiment starts with a current sheet that is in mechanical and thermal equilibrium but is unstable to both tearing-mode and thermal-condensation instabilities. The current sheet is line-tied at one end to a photospheric-like boundary and evolves asymmetrically. The effects of thermal conduction, resistivity variation, and gravity are ignored. In general, reconnection in the nonlinear stage of the tearing-mode instability can strongly affect the onset of condensations unless the radiative-cooling time scale is much smaller than the tearing-mode time scale. When the ambient plasma is less than 0.2, the reconnection enters a regime where the outflow from the reconnection region is supermagnetosonic with respect to the fast-mode wave speed. In the supermagnetosonic regime the most rapidly condensing regions occur downstream of a fast-mode shock that forms where the outflow impinges on closed loops attached to the photospheric-like boundary. A similar shock-induced condensation might occur during the formation of 'postflare' loops.
The present numerical method for accurate and efficient solution of systems of linear equations proceeds by numerically developing a set of basis solutions characterized by slowly varying dependent variables. The solutions thus obtained are shown to have a computational overhead largely independent of the small size of the scale length which characterizes the solutions; in many cases, the technique obviates series solutions near singular points, and its known sources of error can be easily controlled without a substantial increase in computational time.
A numerical study is conducted of a reacting flowfield generated by a single hydrogen jet injected transversely behind a rearward facing step. A nonreacting mixing study was followed by a reacting simulation. The computation simulates a current experimental condition. The validation of the numerical solution was limited to a comparison of the wall-pressure distribution along the wall. The comparison indicates agreement between the experiment and the prediction with a maximum error of 7 percent. The addition of the reaction responsible for HO2 production (H + O2 + M = HO2 + M) to an initially simplified chemical mechanism reduced the exponential growth of the free radicals and reduced the amount of heat released by combustion; as a result, the flowfield was significantly changed. A reversed flow was obtained in limited regions of the flow.
Numerical simulations of an unstable supersonic slab-symmetric jet are described. The instabilities within the jet are characterized by growing internal body waves and their coupled surface waves that are also predicted in linear perturbation theory. The characteristic theory of fluid dynamics is used to help interpret the wave morphologies. It is demonstrated that these waves can be excited by imposing an arbitrary disturbance. From the numerical simulations, it is found that the sound waves propagating against the flow slow down as they propagate outward, and they grow in amplitude. These waves eventually disrupt the jet at a certain length. This disruption length is related to the jet Mach number and the perturbation intensity. Thus, the Mach number of a jet observed with a radio telescope can be estimated by measuring the disruption length and estimating the perturbation intensity. The jet Mach numbers in radio tailed sources determined in this way agree quite well with estimates from ram pressure bending arguments. The wiggles and flares observed in many extragalactic jets, especially in tailed radio sources, appear to be intimately related to instabilities and the jet disruption process.