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.
Finite-element algorithm devised to facilitate solution of two-dimensional Navier-Stokes equations governing kinematics and thermodynamics of variable-viscosity, compressible, multiple-species fluid. Algorithm has been implemented into existing computer program system.
A sequence of approximate solutions for the nonlinear, nonstationary Navier-Stokes equations for a two-dimensional domain, from which explicit error estimates and rates of convergence are obtained, is described. This sequence of approximate solutions is based primarily on the Newton-Kantorovich method.
This paper describes an efficient numerical method for solving the steady incompressible Navier-Stokes equations. The method is a fully implicit method based on the generalized Galerkin method, and the resulting system of equations is solved in a sweeping mode by iterative line relaxation. Results of the present method are compared with published results for separating and reattaching flows, and parametric studies showing the effects of step size, boundary locations, and Reynolds number are presented. The present method is substantially faster than previously published methods; typical run times range from 10 sec to 1 min of 7600 CPU time; and, based on results obtained to date, it is stable at any Reynolds number.
A time dependent, two dimensional Navier-Stokes code employing the method of body fitted coordinate technique was developed for supersonic flows past blunt bodies of arbitrary shapes. The bow shock ahead of the body is obtained as part of the solution, viz., by shock capturing. A first attempt at mesh refinement in the shock region was made by using the forcing function in the coordinate generating equations as a linear function of the density gradients. The technique displaces a few lines from the neighboring region into the shock region. Numerical calculations for Mach numbers 2 and 4.6 and Reynolds numbers from 320 to 10,000 were performed for a circular cylinder with and without a fairing. Results of Mach number 4.6 and Reynolds number 10,000 for an isothermal wall temperature of 556 K are presented in detail.
A time-dependent compressible turbulent Navier-Stokes analysis is applied to the oscillating airfoil flow field problem. The turbulence model is based upon a turbulence energy equation. Results obtained for an airfoil oscillating in pitch are compared to data.
Numerical solutions of the Navier-Stokes equations, with an algebraic turbulence model, for time-dependent two dimensional flow about multi-element airfoils were developed. Fundamental to these solutions was the use of numerically-generated boundary-conforming curvilinear coordinate systems to allow bodies of arbitrary shape to be treated. A general two dimensional grid generation code for multiple-body configuration was written as a part of this project and made available through the COSMIC code library.
A computer code for solving the Reynolds averaged full Navier-Stokes equations has been developed and applied using sheared H-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady state calculations. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications in less than 100 multigrid cycles on the finest mesh.
A vectorized 3D Navier-Stokes code has been implemented on CYBER 205 for solving the supersonic laminar flow over a swept fin/flat plate junction. The code extends MacCormack's predictor-corrector finite volume scheme to a generalized coordinate system in a locally one dimensional time split fashion. A systematic parametric study is conducted to examine the effect of fin sweep on the computed flow field. Calculated results for the pressure distribution on the flat plate and fin leading edge are compared with the experimental measurements of a right angle blunt fin/flat plate junction. The decrease in the extent of the separated flow region and peak pressure on the fin leading edge, and weakening of the two reversed supersonic zones with increase in fin sweep have been clearly observed in the numerical simulation.
Two efficient and robust finite-volume multigrid schemes for solving the Navier-Stokes equations are investigated. These schemes employ either a cell centered or a cell vertex discretization technique. An explicit Runge-Kutta algorithm is used to advance the solution in time. Acceleration techniques are applied to obtain faster steady-state convergence. Accuracy and convergence of the schemes are examined. Computational results for transonic airfoil flows are essentially the same, even for a coarse mesh. Both schemes exhibit good convergence rates for a broad range of artificial dissipation coefficients.
A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.
The fully three-dimensional, time-dependent compressible Navier-Stokes equations in cylindrical coordinates are presently used, in conjunction with the multistage Runge-Kutta numerical integration scheme for solution of the governing flow equations, to simulate complex flowfields within turbomechanical components whose pertinent effects encompass those of viscosity, compressibility, blade rotation, and tip clearance. Computed results are presented for selected cascades, emphasizing the code's capabilities in the accurate prediction of such features as airfoil loadings, exit flow angles, shocks, and secondary flows. Computations for several test cases have been performed on a Cray-YMP, using nearly 90,000 grid points.
A new technique has been developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure crrection method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. In two-dimensions a turbine flowfield, including heat transfer, has been computed with this method and the prediction of the cascade performance is presented. The extension of the pressure correction method for the solution of three-dimensional flows is also presented for laminar flow in an S-shaped duct and turbulent flow in the end-wall region of a turbine cascade.
The feasibility was determined of incorporating the Navier-Stokes computational code, CFL3D, into the supersonic wing design process. The approach taken is of two steps. The first step was to calibrate CFL3D against existing experimental data sets obtained on thin sharp edged delta wings. The experimental data identified six flow types which are dependent on the similarity parameters of Mach number and angle of attack normal to the leading edge. The calibration showed CFL3D capable of simulating these various separated and attached flow conditions. The second step was to use CFL3D to study the initial formation of leading edge separation over delta wings at supersonic speeds. This consisted of examining solutions obtained on a 65 deg delta wing at Mach number of 1.6 with varying cross sectional shapes. Reynolds number was held constant at 1000000 and the Baldwin-Lomax turbulence model was used. The study showed that through the use of leading edge radius and/or camber, the onset of leading edge separation can be delayed to a higher angle of attack than observed on a flat sharp edged wing. Based on the geometries studied, three wind tunnel models are being designed to verify these results.
A solution-adaptive grid procedure for the three-dimensional parabolized Navier-Stokes equations is developed on the basis of the two-dimensional line-by-line solution-adaptive technique of Harvey et al. (1990). It is shown that the new adaptive grid algorithm provides improved shock resolving characteristics over the conventional flow algorithm. The adaptation process has the ability of aligning grid lines with the existing flow-field structure, yielding increasing resolution of high gradient regions.
The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.
The three-dimensional Reynolds-averaged Navier-Stokes equations are numerically solved to obtain the pressure distribution and heat transfer rates on the endwalls and the blades of two linear turbine cascades. The TRAF3D code which has recently been developed in a joint project between researchers from the University of Florence and NASA Lewis Research Center is used. The effect of turbulence is taken into account by using the eddy viscosity hypothesis and the two-layer mixing length model of Baldwin and Lomax. Predictions of surface heat transfer are made for Langston's cascade and compared with the data obtained for that cascade by Graziani. The comparison was found to be favorable. The code is also applied to a linear transonic rotor cascade to predict the pressure distributions and heat transfer rates.
A fully elliptic, control volume solution of the two-dimensional incompressible Navier-Stokes equations for the prediction of cascade performance over a wide incidence range is presented. The numerical technique is based on a new pressure substitution method. A Poisson equation is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. The analysis of a double circular arc compressor cascade is presented, and the results are compared with the available experimental data at various incidence angles. Good agreement is obtained for the blade pressure distribution, boundary layer and wake profiles, skin friction coefficient, losses and outlet angles. Turbulence effects are simulated by the Low-Reynolds-Number version of the k-epsilon turbulence model.
A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.