Engineering topics
Publications and source records attributed to Yu, N. J..
The presentation will show some of our recent work in using the NASA-built Navier-Stokes solvers for various applications including airplane control surface effectiveness study, Reynolds number scaling, and high lift configuration analysis.
An iterative design method has been implemented into 2D and 3D Navier-Stokes codes for the design of airfoils or wings with given target pressure distributions. The method begins with the analysis of an initial geometry, and obtains the analysis pressure distributions of that geometry. The differences between analysis pressures and target pressures are used to drive geometry changes through the use of a streamline curvature method. This paper describes the procedure that makes the iterative design method work for Navier-Stokes codes. Examples of 2D airfoil design, and 3D wing design are included. It is demonstrated that the method is highly effective for airfoil or wing design at flow conditions where no substantial separation occurs. Problems encountered in the airfoil design with shock induced flow separations are discussed.
The feasibility of using boundary layer control (BLC) on a high speed civil transport (HSCT) high lift configuration for low speed performance improvement is studied. The possibility of using a part of the suction system previously designed for cruise laminar flow control (LFC) for low speed BLC mode was explored. A suction scheme for BLC was devised for the HSCT high lift configuration. The extent of the suction regions and levels were determined via panel method based inviscid flow analysis coupled with three-dimensional boundary layer analyses. The compatibility of the BLC suction requirements was assessed, and the modifications necessary to operate the system in both modes of operation, were determined. An assessment of the potential aerodynamic performance improvement resulting from an implementation of the BLC concept compared to an optimized simple flap configuration without BLC was made by an Euler code analysis of the simple flap configuration.
An Euler-based method for aerodynamic analysis of turboprop transport aircraft at transonic speeds has been developed. In this method, inviscid Euler equations are solved over surface-fitted grids constructed about aircraft configurations. Propeller effects are simulated by specifying sources of momentum and energy on an actuator disc located in place of the propeller. A stripwise boundary layer procedure is included to account for the viscous effects. A preliminary version of an approach to embed the exhaust plume within the global Euler solution has also been developed for more accurate treatment of the exhaust flow. The resulting system of programs is capable of handling wing-body-nacelle-propeller configurations. The propeller disks may be tractors or pushers and may represent single or counterrotation propellers. Results from analyses of three test cases of interest (a wing alone, a wing-body-nacelle model, and a wing-nacelle-endplate model) are presented. A user's manual for executing the system of computer programs with formats of various input files, sample job decks, and sample input files is provided in appendices.
The generalized minimal residual algorithm (GMRES) is a conjugate-gradient like method that applies directly to nonsymmetric linear systems of equations. In this paper, GMRES is modified to handle nonlinear equations characteristic of computational fluid dynamics. Attention is devoted to the concept of preconditioning and the role it plays in assuring rapid convergence. A formulation is developed that allows GMRES to be preconditioned by the solution procedures already built into existing computer codes. Examples are provided that demonstrate the ability of GMRES to greatly improve the robustness and rate of convergence of current state-of-the-art fluid dynamics codes. Theoretical aspects of GMRES are presented that explain why it works. Finally, the advantage GMRES enjoys over related methods such as conjugate gradients are discussed.
Euler codes for both axisymmetric and general three-dimensional nacelle-propeller flow analysis have been developed. Surface-fitted grids are generated either by a finite difference method or by an algebraic method. The propeller is represented by an actuator disk along a computational plane where proper boundary conditions are assigned to simulate the propeller power loading. Computed results for a NASA SR3 propeller as well as a NASA turboprop configuration are compared with test data. Good agreement has been achieved through the present simulation method.
An Euler code has been developed for the analysis of a wing-mounted propfan configuration. Surface-fitted grids are used to represent the wing, fuselage, and nacelle geometry. The propeller is simulated by an actuator disk along a computational plane. A grid embedding technique is employed to capture detailed flow field resolution in the vicinity of the engine exhaust plume. Results of a NASA turboprop configuration are compared with test data. Specific issues on grid embedding and methods of resolving them are discussed.
This paper presents the results of a study that used a three-dimensional transonic analysis method to analyze flow around a wing-mounted prop-fan configuration. The transonic analysis method employs a cell-oriented finite volume approach to solve the full potential equation in conservative form over a numerically generated surface-fitted grid. The propeller slipstream effects are simulated using a linearized transpiration boundary condition. The results of analyses of flow over clean wing/body, wing/body/nacelle, and wing/body nacelle with slipstream are compared to test data.
A suitable grid system for complex three dimensional configurations such as wing/body/nacelle shape for the solution of nonlinear transonic flow problems was constructed. Two approaches were explored based on Thompson's body fitted coordinate concept. The most general approach was to divide the computational domain into multiple rectangular blocks where the configuration itself was also represented by a set of blocks whose structure follows the natural lines of the configuration. The block structured grid system was adaptable to complex configurations and gives good grid quality near physical corners. However, it introduced algorithm issues for the flow solution concerning the treatment of nonanalytic grid block boundaries and nonstandard grid cells. These issues were explored in relation to the grid generation. A more limited approach treats a wing/body configuration with only a single rectangular block in computational space. The issues involving nonstandard cells were avoided, but other limitations on grid resolution appear. Both a linear and a nonlinear system of grid generation equations were developed including methods of grid control.
A new method for the design of shock-free supercritical airfoils, wings, and three-dimensional configurations is described. Results illustrating this procedure in two and three dimensions are given. They include modifications to part of the upper surface of an NACA 64A410 airfoil that will maintain shock-free flow over a range of Mach numbers for a fixed lift coefficient, and the modifications required on part of the upper surface of a swept wing with an NACA 64A410 root section to achieve shock-free flow. While the results are given for inviscid flow, the same procedures can be employed iteratively with a boundary layer calculation in order to achieve shock-free viscous designs. With a shock-free pressure field the boundary layer calculation will be reliable and not complicated by the difficulties of shock-wave boundary-layer interaction.
A method for the design of shock free supercritical airfoils, wings, and three dimensional configurations is described. Results illustrating the procedure in two and three dimensions are given. They include modifications to part of the upper surface of an NACA 64A410 airfoil that will maintain shock free flow over a range of Mach numbers for a fixed lift coefficient, and the modifications required on part of the upper surface of a swept wing with an NACA 64A410 root section to achieve shock free flow. While the results are given for inviscid flow, the same procedures can be employed iteratively with a boundary layer calculation in order to achieve shock free viscous designs. With a shock free pressure field the boundary layer calculation will be reliable and not complicated by the difficulties of shock wave boundary layer interaction.
The numerical procedures previously developed for computing nonlinear and time-linearized small-perturbation unsteady transonic flows are briefly reviewed, and the effects of unsteady modes of motion on two-dimensional transonic flows are evaluated. The numerical procedure used comprises an alternating-direction implicit scheme and treats shock waves as discontinuities in the flow. Comparison of the time-linearized results with fully nonlinear calculations delineates their range of applicability. The unsteady behavior due to harmonic pitching and flap oscillations of an NACA airfoil is also examined.
The alternating-direction implicit scheme developed by NASA Ames for unsteady transonic flows has been modified to include a shock-fitting algorithm as well as an analytically stretched coordinate system. The shock-fitting procedure treats shock waves as discontinuities normal to the free stream. Improvements in shock position and the unsteady pressure distributions are obtained by this modification. The various types of shock motion observed experimentally by Tijdeman are well simulated in calculations using the modified computational scheme. The method of detecting shock wave formation and the procedure for fitting a moving shock wave are illustrated. Results for a pulsating parabolic arc airfoil and for an NACA 64A006 airfoil with oscillating quarter-chord flap are presented and discussed.
The results of a systematic study of small perturbation transonic flows are presented. Both the flow over thin airfoils and the flow over wedges were investigated. Various numerical schemes were employed in the study. The prime goal of the research was to determine the efficiency of various numerical procedures by accurately evaluating the wave drag, both by computing the pressure integral around the body and by integrating the momentum loss across the shock. Numerical errors involved in the computations that affect the accuracy of drag evaluations were analyzed. The factors that effect numerical stability and the rate of convergence of the iterative schemes were also systematically studied.
First- and second-order numerical procedures for calculating two-dimensional transonic flows that treat shock waves as discontinuities are discussed. This short communication illustrates their application to a simple but non-trivial problem for which there are limited theoretical results.
First-and second-order numerical procedures are presented for calculating two-dimensional transonic flows that treat shock waves as discontinuities. Their application to a simple but nontrivial problem for which there are limited theoretical results is discussed.
This paper discusses the procedures we have developed to treat a canonical problem involving a mixed nonlinear equation with boundary data that imply a discontinuous solution. This equation arises in various physical contexts and is basic to the description of the nonlinear acoustic behavior of a shock wave near a caustic. The numerical scheme developed is of second order, treats discontinuities as such by applying the appropriate jump conditions across them, and eliminates the numerical dissipation and dispersion associated with large gradients. Our results are compared with the results of a first-order scheme and with those of a second-order scheme we have developed. The algorithm used here can easily be generalized to more complicated problems, including transonic flows with imbedded shocks.