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 objective is to develop efficient numerical techniques for the study of aeroelastic response of a propfan in an unsteady transonic flow. A three dimensional unsteady Euler solver is being modified to address this problem.
The fundamentals and the state-of-the-art technology for unstructured grids and solvers are highlighted. Algorithms and techniques pertinent to mesh generation are discussed. It is shown that grid generation and grid manipulation schemes rely on fast multidimensional searching. Flow solution techniques for the Euler equations, which can be derived from the integral form of the equations are discussed. Sample calculations are also provided.
Euler solutions for steady transonic flow (Free-stream Mach 0.63-0.8, alpha = 0-2 deg) over NACA 0012 and supercritical airfoils with solid as well as porous surfaces suggest porosity as a means to realize multipoint design for transonic airfoils. The porous surfaces extend over at least 90 percent of the chord. The porosity distribution is described by a modified sine wave with several amplitudes. Either connected or separated cavities are assumed to lie underneath the upper and lower surfaces. Applied to an NACA 0012 airfoil, porosity generally increases lift, in some instances by up to 65 percent. Porous NACA 0012 airfoils in supercritical flow yield reductions of an order of magnitude in wave drag at constant lift, compared to their solid counterpart. Making the surface of a supercritical airfoil permeable also leads to sizeable reductions in wave drag at constant lift for overspeed conditions. The discussion of the computed results addresses issues such as grid sensitivity and checks for systematic errors.
The capability of generating 2-D unstructured triangular meshes about arbitrary geometries is demonstrated. This work uses a distribution of boundary points and triangulates the computational domain using a Delaunay triangulation algorithm. Typically, initial cells are added based on cell aspect ratios or cell areas. A resulting mesh can then be used along with the connectivity of the cells to solve either a Euler or Navier-Stokes flow problem.
The research described in this paper is a numerical investigation of the effects of unsteady flow on gas turbine heat transfer, particularly on a rotor blade surface. The unsteady flow in a rotor blade passage and the unsteady heat transfer on the blade surface as a result of wake/blade interaction are modeled by the inviscid flow/boundary layer approach. The Euler equations which govern the inviscid flow are solved using a time-accurate marching scheme. The unsteady flow in the blade passage is induced by periodically moving a wake model across the passage inlet. Unsteady flow solutions in the passage provide pressure gradients and boundary conditions for the boundary-layer equations which govern the viscous flow adjacent to the blade surface. Numerical solutions of the unsteady turbulent boundary layer yield surface heat flow values which can then be compared to experimental data. Comparisons with experimental data show that unsteady heat flux on the blade suction surface is well predicted, but the predictions of unsteady heat flux on the blade pressure surface do not agree.
A multigrid diagonal implicit algorithm has been developed to solve the three-dimensional Euler equations of inviscid compressible flow on block-structured grids. An improved method of advancing the multigrid cycle has been examined with respect to convergence rates, accuracy, and efficiency. In this method, the multigrid cycle is advanced independently in each of the blocks, and the information exchange between the blocks is done using buffer arrays, allowing for the asynchronous updating of interface boundary conditions. This updating scheme is used to eliminate the convergence problems found in a previous implementation of the algorithm while retaining its potential for efficient parallel execution. Results are computed for transonic flows past wings and include pressure distributions to verify the accuracy of the scheme and convergence histories to demonstrate the efficiency of the method. Efficiencies that were obtained using a modest number of processors in parallel are also presented and discussed.
The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.
A symmetric Total-Variation-Diminishing (TVD) formulation of the numerical dissipation terms has been incorporated into a diagonalized alternating direction implicit multigrid algorithm to solve the Euler equations of inviscid compressible flow. The new treatment of the dissipation makes is possible to capture both very strong and very weak shocks, virtually without oscillation for the steady flows of interest here. In addition, the TVD constraint fixes one of the two previously arbitrary constants in the formulation of the dissipation, and results in both converged solutions and convergence rates which are relatively insensitive to the choice of the remaining dissipation parameter.
A conference was held on parallel computational fluid dynamics and produced related papers. Topics discussed in these papers include: parallel implicit and explicit solvers for compressible flow, parallel computational techniques for Euler and Navier-Stokes equations, grid generation techniques for parallel computers, and aerodynamic simulation om massively parallel systems.
Observations regarding the use of advanced computational fluid dynamics (CFD) analysis, sensitivity analysis (SA), and design codes in gradient-based multidisciplinary design optimization (MDO) reflect our perception of the interactions required of CFD and our experience in recent aerodynamic design optimization studies using CFD. Sample results from these latter studies are summarized for conventional optimization (analysis - SA codes) and simultaneous analysis and design optimization (design code) using both Euler and Navier-Stokes flow approximations. The amount of computational resources required for aerodynamic design using CFD via analysis - SA codes is greater than that required for design codes. Thus, an MDO formulation that utilizes the more efficient design codes where possible is desired. However, in the aerovehicle MDO problem, the various disciplines that are involved have different design points in the flight envelope; therefore, CFD analysis - SA codes are required at the aerodynamic 'off design' points. The suggested MDO formulation is a hybrid multilevel optimization procedure that consists of both multipoint CFD analysis - SA codes and multipoint CFD design codes that perform suboptimizations.
We solve the problem of determining airfoils that approximate, in a least square sense, given surface pressure distributions in transonic flight regimes. The flow is modeled by means of the Euler equations and the solution procedure is an adjoint- based minimization algorithm that makes use of the inverse Theodorsen transform in order to parameterize the airfoil. Fast convergence to the optimal solution is obtained by means of the pseudo-time method. Results are obtained using three different pressure distributions for several free stream conditions. The airfoils obtained have given a trailing edge angle.
Several examples from the past decade of success stories involving the design and flight test of three true X-planes will be described: in particular, X-plane design techniques that relied heavily upon computational fluid dynamics (CFD). Three specific examples chosen from the author s personal experience are presented: the X-36 Tailless Fighter Agility Research Aircraft, the X-45A Unmanned Combat Air Vehicle, and, most recently, the X-48B Blended Wing Body Demonstrator Aircraft. An overview will be presented of the uses of CFD analysis, comparisons and contrasts with wind tunnel testing, and information derived from the CFD analysis that directly related to successful flight test. Some lessons learned on the proper application, and misapplication, of CFD are illustrated. Finally, some highlights of the flight-test results of the three example X-planes will be presented. This overview paper will discuss some of the author s experience with taking an aircraft shape from early concept and three-dimensional modeling through CFD analysis, wind tunnel testing, further refined CFD analysis, and, finally, flight. An overview of the key roles in which CFD plays well during this process, and some other roles in which it does not, are discussed. How wind tunnel testing complements, calibrates, and verifies CFD analysis is also covered. Lessons learned on where CFD results can be misleading are also given. Strengths and weaknesses of the various types of flow solvers, including panel methods, Euler, and Navier-Stokes techniques, are discussed. The paper concludes with the three specific examples, including some flight test video footage of the X-36, the X-45A, and the X-48B.
Examples of the design and flight test of three true X-planes are described, particularly X-plane design techniques that relied heavily on computational fluid dynamics(CFD) analysis. Three examples are presented: the X-36 Tailless Fighter Agility Research Aircraft, the X-45A Unmanned Combat Air Vehicle, and the X-48B Blended Wing Body Demonstrator Aircraft. An overview is presented of the uses of CFD analysis, comparison and contrast with wind tunnel testing, and information derived from CFD analysis that directly related to successful flight test. Lessons learned on the proper and improper application of CFD analysis are presented. Highlights of the flight-test results of the three example X-planes are presented. This report discusses developing an aircraft shape from early concept and three-dimensional modeling through CFD analysis, wind tunnel testing, further refined CFD analysis, and, finally, flight. An overview of the areas in which CFD analysis does and does not perform well during this process is presented. How wind tunnel testing complements, calibrates, and verifies CFD analysis is discussed. Lessons learned revealing circumstances under which CFD analysis results can be misleading are given. Strengths and weaknesses of the various flow solvers, including panel methods, Euler, and Navier-Stokes techniques, are discussed.
A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.
The present consideration of procedures for the definition of boundary conditions for the Navier-Stokes equations emphasizes the derivation of boundary conditions that are compatible with nondissipative algorithms applicable to direct simulations of turbulent flows. A novel formulation for the Euler equations is derived on the basis of characteristic wave relations through boundaries; this formulation is generalized to the Navier-Stokes equations. The method, which applies to both sub- and supersonic flows, is used in reflecting and nonreflecting boundary-condition treatments. Attention is given to practical implementations involving inlet and outlet boundaries and slip and nonslip walls, as well as the test cases of a ducted shear layer, vortices propagating through boundaries, and Poiseuille flow.
Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.
A comparison of flow separation in transonic flows is made using various computational schemes which solve the Euler and the Navier-Stokes equations of fluid mechanics. The flows examined are computed using several simple two-dimensional configurations including a backward facing step and a bump in a channel. Comparison of the results obtained using shock fitting and flux vector splitting methods are presented and the results obtained using the Euler codes are compared to results on the same configurations using a code which solves the Navier-Stokes equations.
The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.