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.
Papers and discussions on surface representation and its integration with aerodynamics, computers, graphics, wind tunnel model fabrication, and flow field grid generation are presented. Surface definition is emphasized.
Digitally acquired and processed results from an experimental investigation of grid generated turbulence of various scales through and downstream of nine matched cubic contour contractions ranging in area ratio from 2 to 36, and in length to inlet diameter ratio from 0.25 to 1.50 are reported. An additional contraction with a fifth order contour was also utilized for studying the shape effect. Thirteen homogeneous and nearly isotropic test flow conditions with a range of turbulence intensities, length scales and Reynolds numbers were generated and used to examine the sensitivity of the contractions to upstream turbulence. The extent to which the turbulence is altered by the contraction depends on the incoming turbulence scales, the total strain experienced by the fluid, as well as the contraction ratio and the strain rate. Varying the turbulence integral scale influences the transverse turbulence components more than the streamwise component. In general, the larger the turbulence scale, the lesser the reduction in the turbulence intensity of the transverse components. Best agreement with rapid distortion theory was obtained for large scale turbulence, where viscous decay over the contraction length was negligible, or when a first order correction for viscous decay was applied to the results.
A method for constructing a direct functional relation between the computational and physical domains such that boundaries in the computational domain map into boundaries in the physical domain is presented. This method which is called the 'two-boundary technique' is based on generating functions that represent two disconnected boundaries in the physical domain and an independent function which connects one boundary with another to generate grids. The suggested connecting functions are simple linear and cubic parametric polynomials with parameters based on position and derivatives at the disconnected boundaries. The boundaries can be defined by analytical, approximate, and tension spline functions; the grid density in specified regions is controlled by embedding contraction functions in either the boundary representation or the connecting function.
An implicit finite-difference code is developed to solve either inviscid or viscous flow about two-dimensional cascade blade elements. General coordinate transformations are used so that boundaries can coincide with coordinate lines, and an automatic grid generation routine based on elliptic partial differential equations is employed to mesh arbitrary cascade elements. Characteristic combinations of the differential equations are used at inflow and outflow boundaries. Computed results for both inviscid and viscous flow are compared with other existing cascade solutions and experimental data.
Results of hot-wire measurements in an incompressible partially confined jet issuing from an array of rectangular nozzles, equally spaced with their small dimensions aligned are presented. The quantities measured include mean velocity and the Reynolds stress in the two central planes of the jet at stations covering up to 115 widths (small dimension of a nozzle) downstream of the nozzle exit. For downstream distances greater than 60 widths, the flowfield is observed to be nearly homogenous and the turbulence appears to be quite similar to that of a grid generated turbulence.
An evaluation of the transonic-wing-analysis computer code TWING is given. TWING utilizes a fully implicit approximate factorization iteration scheme to solve the full potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations were analyzed, and the limits of applicability of this code was evaluated. Comparisons of computed results were made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative full potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.
A fast, fully implicit, approximate factorization algorithm designed to solve the conservative full-potential equation is used to compute lifting-wing flow fields with embedded transonic flow. The computer program (TWING - transonic wing analysis) uses an elliptic-solver numerical grid-generation routine and is capable of analyzing arbitrary wing geometries with sweep, twist, taper, and section variation. Computed flowfield results for several wing geometries are presented; they demonstrate substantial improvement in convergence speed relative to the FL028 computer code.
A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.
The Lagrangian dispersion theory of Durbin (1980) is used to analyze experiments by Warhaft and Lumley (1978) and by Sreenivasan et al. (1980) on temperature fluctuations in grid-generated turbulence. Both theory and experiment show that the decay exponent m depends on the ratio of the initial length scales of velocity and temperature, although when this ratio is greater than 2.5 such dependence is negligible. The theory shows that m is not truly constant, but within the range covered by the experiments it is nearly so. The agreement between theory and experiment lends credence to the idea that the decay of fluctuations is controlled largely by turbulent relative dispersion.
Recent developments with finite-difference techniques are emphasized. The quotation marks reflect the fact that any finite discretization procedure can be included in this category. Many so-called finite element collocation and galerkin methods can be reproduced by appropriate forms of the differential equations and discretization formulas. Many of the difficulties encountered in early Navier-Stokes calculations were inherent not only in the choice of the different equations (accuracy), but also in the method of solution or choice of algorithm (convergence and stability, in the manner in which the dependent variables or discretized equations are related (coupling), in the manner that boundary conditions are applied, in the manner that the coordinate mesh is specified (grid generation), and finally, in recognizing that for many high Reynolds number flows not all contributions to the Navier-Stokes equations are necessarily of equal importance (parabolization, preferred direction, pressure interaction, asymptotic and mathematical character). It is these elements that are reviewed. Several Navier-Stokes and parabolized Navier-Stokes formulations are also presented.
The present effectiveness and future promise of computational fluid dynamics (CFD) are discussed. The CFD ideal of computer-designed flight vehicles and the progress that has been made toward its achievement, expressed in terms of primary and secondary pacing items, are also addressed. The primary pacing items listed by Chapman - grid generation, turbulence modeling, solution methodology development, and mainframe computer methods and architecture - are modified and reviewed. Important secondary pacing items in CFD development are also examined: algorithm development, complex geometry definition, and predata and postata processing. In addition, user requirements placed on CFD computer codes are defined and some problems associated with the dissemination of CFD technology are identified.
A viscous-inviscid interactive calculation procedure is developed for application to flow in cascades of two-dimensional airfoils. This procedure has essentially three components. First, a numerical solution of the Euler equations which can accommodate an arbitrarily specified cascade geometry is carried out on a nonorthogonal curvileanear grid mesh that is fitted to the geometry of the cascade. A method of grid generation has been used which relies in part on a succession of conformal mappings. Second, a viscous solution for use in boundary layer and wake regions has been programmed. Finally, an interactive scheme which takes the form of a source-sink distribution along the blade surface and wake centerline is employed. Results have been obtained with this procedure for several cascade flow situations, and some comparisons with experiment are presented.
Recent progress in computational methods for time dependent fluid dynamics is presented. The emphasis is on advances applicable to large scale systems with the connection between the numerics and the physics of the code stressed. All aspects of a working code are discussed including such topics as initialization, boundary conditions, grid generation in addition to algorithmic advances. One sometimes uses a time dependent method as an iteration procedure to reach a steady state solution. Work in accelerating the convergence to the steady state is also surveyed.
A viscous-inviscid interactive calculation procedure is developed for application to flow in cascades of two-dimensional airfoils. This procedure has essentially three components. First, a numerical solution of the Eulers equations which can accommodate an arbitrarily specified cascade geometry of the cascade. A method of grid generation has been used which relics in part on a succession of conformal mappings. Second, a viscous solution for use in boundary layers and wake regions was programmed. Finally, an interactive scheme which takes the form of a source-sink distribution along the blade surface and wake centerline is employed. Results were obtained with this procedure for several cascade flow situations, and some comparisons with experiment are presented.
A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.
Basic concepts associated with the numerical solution of elliptic partial differential equations are introduced, and procedures used to solve the full potential equation for transonic flow fields are discussed. Governing equations, classical relaxation schemes and concepts regarding transonic, full potential equation algorithms are covered. The equation transformation and grid generation procedures; full potential spatial differencing schemes; full potential iteration schemes, emphasizing convergence acceleration; and three dimensional applications are presented.
Seven low cost multi-100 kW planar solar array modules were fabricated and tested. Two different designs were used, demonstrating advanced solar array construction practices. Both module types utilized second generation gridded back cells featuring high efficiency and IR transparency. A silicon dioxide AR coating optimized for transmission at gamma = 1.7 microns was applied to the back surface. Two interconnect types, a single sheet printed circuit and a roll type, with alternate approaches to increasing transparency and reducing cost were designed and fabricated. Hinge stress and electrical power optimization were also examined. Two point designs were studied. The first design used a coilage longeron mast and is autonomously deployable. The second design used a Stac Beam for high natural frequency response and required astronaut assistance and assembly on orbit. It was conclusively demonstrated that planar arrays are the most cost effective design for use on the space station or other high power applications.
The CMC fluid mechanics program system was developed to transmit the theoretical solution of finite element numerical solution methodology, applied to nonlinear field problems into a versatile computer code for comprehensive flow field analysis. Data procedures for the CMC 3 dimensional Parabolic Navier-Stokes (PNS) algorithm are presented. General data procedures a juncture corner flow standard test case data deck is described. A listing of the data deck and an explanation of grid generation methodology are presented. Tabulations of all commands and variables available to the user are described. These are in alphabetical order with cross reference numbers which refer to storage addresses.