SEARCH · Search NASA
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.
An explicit-implicit and an implicit two-dimensional Navier-Stokes code along with various grid generation capabilities were developed. A series of classical benckmark cases were simulated using these codes.
A thin-layer Navier-Stokes code capable of predicting steady-state viscous flows is applied to the transonic flow over a Space Shuttle configuration. The code is written in the generalized coordinate system, and the grid-generation code of Fujii (1983) is used for the discretization of the flow field. The flow-field computation is done using the CRAY 1S computer at NASA Ames. The computed result is physically reasonable, even though no experimental data is available for the comparison purpose.
The numerical simulation of flow about ice accretion shapes has been accomplished by solving the Navier-Stokes equations using MacCormack's explicit finite difference scheme. The computations were performed on a CRAY-XMP computer. The influence of turbulence is taken into account by means of an algebraic eddy-viscosity model. In order to optimize the grid spacing and to achieve near orthogonality at the surface of the complex ice shapes, a hyperbolic grid generation scheme is utilized. Particular attention is given to the heat transfer process for which good agreement between the numerical and experimental results is achieved. In addition, liquid water droplet trajectories are coupled within the flowfield along with the resulting collection efficiencies using a parabolized Navier-Stokes formulation.
Analysis of a NACA0012 airfoil with leading edge ice has been performed using a Navier-Stokes code coupled with a grid generation code. The computed results were compared to experimental information obtained for an airfoil with a well defined artificial ice shape. The computations were performed at angles of attack ranging from zero to ten degrees. This range is sufficient to show the development of the separation bubble aft of the ice shape on both the upper and lower surfaces. Velocity profile plots in the separation bubble are examined in order to determine if recirculation patterns are predicted properly and if separation and reattachment points are found within the resolution of the experimental information. Also, the massive separation near the point of stall is examined in order to more accurately evaluate the lift coefficient curve in that region. Lift, drag, and moment coefficients are computed and compared to experiment.
Two optimization methods for solving a variational problem in grid generation are described and evaluated. The smoothness, cell volumes, and orthogonality of the variational integrals are examined. The Jacobi-Newton iterative method is compared to the Fletcher-Reeves conjugate gradient method. It is observed that a combination of the Jacobi-Newton iteration and the direct solution of the variational problem produces an algorithm which is easy to program and requires less storage and computer time/iteration than the conjugate gradient method.
An approach to analyzing CFD knowledge-based systems is proposed which is based, in part, on the concept of knowledge-level analysis. Consideration is given to the expert cooling fan design system, the PAN AIR knowledge system, grid adaptation, and expert zonal grid generation. These AI/CFD systems demonstrate that current AI technology can be successfully applied to well-formulated problems that are solved by means of classification or selection of preenumerated solutions.
Report describes calculations of transonic airflows about wing in wind tunnel. Basic equations of flow used in study are Reynolds-averaged Navier-Stokes equations in strong conservation-law form. Equations of flow incorporated into finite-difference computer code called TNS (Transonic Navier-Stokes). Computational grid generated by solution of partial differential equations yielding smooth meshes conforming to surfaces of wing and wind tunnel.
A turbulence generating grid was designed and installed in the turbine inlet which produced the target nominal value of 10 percent free stream turbulence. Aerodynamic documentation of the rotor and stator midspan surface pressure distributions were obtained. Midspan heat transfer data were obtained on the rotor and stator for variations in inlet turbulence, rotor-stator axial spacing, and rotor incidence.
An existing Navier-Stokes code (PARC2D) was used to compute the nozzle flow field. Grids were generated by the interactive grid generator codes TBGG and GENIE. All computations were made on the NASA/MSFC CRAY X-MP computer. Comparisons were made between the computations and MSFC in-house wall pressure measurements for CO2 flow through a conical nozzle having an area ratio of 40. Satisfactory agreements exist between the computations and measurements for different stagnation pressures of 29.4, 14.7, and 7.4 psia, at stagnation temperature of 1060 R. However, agreements did not match precisely near the nozzle exit. Several reasons for the lack of agreement are possible. The computational code assumes a constant gas gamma, whereas the gamma i.e. the specific heat ratio for CO2 varied from 1.22 in the plenum chamber to 1.38 at the nozzle exit. The computations also assumes adiabatic and no-slip walls. Both assumptions may not be correct. Finally, it is possible that condensation occurs during the nozzle expansion at the low stagnation pressure. The next phase of the work will incorporate variable gamma and slip wall boundary conditions in the computational code and develop a more accurate computer code.
An adaptive method to solve partial differential equations in fluid mechanics is presented. The approach requires internal boundary conditions that must be conservative, data structures for keeping track of several layers of fine grid patches, error estimation, and heuristics for automatic grid generation. In practical calculations gains in computer efficiency up to 10 over nonadaptive methods are observed. The whole procedure takes 3000 lines of FORTRAN code.
A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.
A computer code called NCOREL (for Nonconical Relaxation) has been developed to solve for supersonic full potential flows over complex geometries. The method first solves for the conical at the apex and then marches downstream in a spherical coordinate system. Implicit relaxation techniques are used to numerically solve the full potential equation at each subsequent crossflow plane. Many improvements have been made to the original code including more reliable numerics for computing wing-body flows with multiple embedded shocks, inlet flow through simulation, wake model and entropy corrections. Line relaxation or approximate factorization schemes are optionally available. Improved internal grid generation using analytic conformal mappings, supported by a simple geometric Harris wave drag input that was originally developed for panel methods and internal geometry package are some of the new features.
Automation flow field zoning in two dimensions is an important step towards easing the three-dimensional grid generation bottleneck in computational fluid dynamics. A knowledge based approach works well, but certain aspects of flow field zoning make the use of such an approach challenging. A knowledge based flow field zoner, called EZGrid, was implemented and tested on representative two-dimensional aerodynamic configurations. Results are shown which illustrate the way in which EZGrid incorporates the effects of physics, shape description, position, and user bias in a flow field zoning.
A transonic equivalent strip (TES) method was further developed for unsteady flow computations of arbitrary wing planforms. The TES method consists of two consecutive correction steps to a given nonlinear code such as LTRAN2; namely, the chordwise mean flow correction and the spanwise phase correction. The computation procedure requires direct pressure input from other computed or measured data. Otherwise, it does not require airfoil shape or grid generation for given planforms. To validate the computed results, four swept wings of various aspect ratios, including those with control surfaces, are selected as computational examples. Overall trends in unsteady pressures are established with those obtained by XTRAN3S codes, Isogai's full potential code and measured data by NLR and RAE. In comparison with these methods, the TES has achieved considerable saving in computer time and reasonable accuracy which suggests immediate industrial applications.
Topics addressed include: wind tunnel and flight experiments; computational fluid dynamics (CFD) applications, industry overviews; and inviscid methods and grid generations.
In order to assess the state of the art in transonic flow disciplines and to glimpse at future directions, NASA-Langley held a Transonic Symposium. Emphasis was placed on steady, three dimensional external, transonic flow and its simulation, both numerically and experimentally. The symposium included technical sessions on wind tunnel and flight experiments; computational fluid dynamic applications; inviscid methods and grid generation; viscous methods and boundary layer stability; and wind tunnel techniques and wall interference. This, being volume 1, is unclassified.
The results of an investigation into the feasibility of using the MPP for direct and large eddy simulations of the Navier-Stokes equations is presented. A major part of this study was devoted to the implementation of two of the standard numerical algorithms for CFD. These implementations were not run on the Massively Parallel Processor (MPP) since the machine delivered to NASA Goddard does not have sufficient capacity. Instead, a detailed implementation plan was designed and from these were derived estimates of the time and space requirements of the algorithms on a suitably configured MPP. In addition, other issues related to the practical implementation of these algorithms on an MPP-like architecture were considered; namely, adaptive grid generation, zonal boundary conditions, the table lookup problem, and the software interface. Performance estimates show that the architectural components of the MPP, the Staging Memory and the Array Unit, appear to be well suited to the numerical algorithms of CFD. This combined with the prospect of building a faster and larger MMP-like machine holds the promise of achieving sustained gigaflop rates that are required for the numerical simulations in CFD.
Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.