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.
The primary tasks performed were the continued development of inverse design procedures for the TAWFIVE code, the development of corresponding relofting and trailing edge closure procedures, and the testing of the methods for a variety of cases. The period from July 1, 1986 through December 31, 1986 is covered.
A fast, efficient and user friendly inverse design system for 3-D nacelles was developed. The system is a product of a 2-D inverse design method originally developed at NASA-Langley and the CFL3D analysis code which was also developed at NASA-Langley and modified for nacelle analysis. The design system uses a predictor/corrector design approach in which an analysis code is used to calculate the flow field for an initial geometry, the geometry is then modified based on the difference between the calculated and target pressures. A detailed discussion of the design method, the process of linking it to the modified CFL3D solver and its extension to 3-D is presented. This is followed by a number of examples of the use of the design system for the design of both axisymmetric and 3-D nacelles.
The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.
A model inverse design problem is used to investigate the effect of flow discontinuities on the optimization process. The optimization involves finding the cross-sectional area distribution of a duct that produces velocities that closely match a targeted velocity distribution. Quasi-one-dimensional flow theory is used, and the target is chosen to have a shock wave in its distribution. The objective function which quantifies the difference between the targeted and calculated velocity distributions may become non-smooth due to the interaction between the shock and the discretization of the flowfield. This paper offers two techniques to resolve the resulting problems for the optimization algorithms. The first, shock-fitting, involves careful integration of the objective function through the shock wave. The second, coordinate straining with shock penalty, uses a coordinate transformation to align the calculated shock with the target and then adds a penalty proportional to the square of the distance between the shocks. The techniques are tested using several popular sensitivity and optimization methods, including finite-differences, and direct and adjoint discrete sensitivity methods. Two optimization strategies, Gauss-Newton and sequential quadratic programming (SQP), are used to drive the objective function to a minimum.
Three optimization-based methods for solving aerodynamic design problems are compared. The Euler equations for one-dimensional duct flow was used as a model problem, and the three methods are compared for efficiency, robustness, and implementation difficulty. The smoothness of the design problem with respect to different shock-capturing finite difference schemes, and in the presence of grid refinement, is investigated.
An inverse wing design method was developed around an existing transonic wing analysis code. The original analysis code, TAWFIVE, has as its core the numerical potential flow solver, FLO30, developed by Jameson and Caughey. Features of the analysis code include a finite-volume formulation; wing and fuselage fitted, curvilinear grid mesh; and a viscous boundary layer correction that also accounts for viscous wake thickness and curvature. The development of the inverse methods as an extension of previous methods existing for design in Cartesian coordinates is presented. Results are shown for inviscid wing design cases in super-critical flow regimes. The test cases selected also demonstrate the versatility of the design method in designing an entire wing or discontinuous sections of a wing.
A novel method for implementing aerodynamic data dispersion analysis is herein introduced. A general mathematical approach combined with physical modeling tailored to the aerodynamic quantity of interest enables the generation of more realistically relevant dispersed data and, in turn, more reasonable flight simulation results. The method simultaneously allows for the aerodynamic quantities and their derivatives to be dispersed given a set of non-arbitrary constraints, which stresses the controls model in more ways than with the traditional bias up or down of the nominal data within the uncertainty bounds. The adoption and implementation of this new method within the NASA Ares I Crew Launch Vehicle Project has resulted in significant increases in predicted roll control authority, and lowered the induced risks for flight test operations. One direct impact on launch vehicles is a reduced size for auxiliary control systems, and the possibility of an increased payload. This technique has the potential of being applied to problems in multiple areas where nominal data together with uncertainties are used to produce simulations using Monte Carlo type random sampling methods. It is recommended that a tailored physics-based dispersion model be delivered with any aerodynamic product that includes nominal data and uncertainties, in order to make flight simulations more realistic and allow for leaner spacecraft designs.
The work presented in this paper describes the application of a multiblock gridding strategy to the solution of aerodynamic design optimization problems involving complex configurations. The design process is parallelized using the MPI (Message Passing Interface) Standard such that it can be efficiently run on a variety of distributed memory systems ranging from traditional parallel computers to networks of workstations. Substantial improvements to the parallel performance of the baseline method are presented, with particular attention to their impact on the scalability of the program as a function of the mesh size. Drag minimization calculations at a fixed coefficient of lift are presented for a business jet configuration that includes the wing, body, pylon, aft-mounted nacelle, and vertical and horizontal tails. An aerodynamic design optimization is performed with both the Euler and Reynolds Averaged Navier-Stokes (RANS) equations governing the flow solution and the results are compared. These sample calculations establish the feasibility of efficient aerodynamic optimization of complete aircraft configurations using the RANS equations as the flow model. There still exists, however, the need for detailed studies of the importance of a true viscous adjoint method which holds the promise of tackling the minimization of not only the wave and induced components of drag, but also the viscous drag.
A new inverse method for aerodynamic design of airfols is presented for subcritical flows. The pressure distribution in this method can be prescribed as a function of the arc length of the as-yet unknown body. This inverse problem is shown to be mathematically equivalent to solving only one nonlinear boundary value problem subject to known Dirichlet data on the boundary. The solution to this problem determines the airfoil, the freestream Mach number, and the upstream flow direction. The existence of a solution to a given pressure distribution is discussed. The method is easy to implement and extremely efficient. A series of results for which comparisons are made with the known airfoils is presented.
This paper deals with selected linearized-aerodynamic and vortex-flow methods as applied to aircraft design problems at high subsonic speeds. In particular, the NASA Vortex Lattice and Modified Multhopp methods are the linearized techniques employed, and the suction analogy is used to provide estimates associated with vortex-flow aerodynamics. Many examples are given as to how researchers at Langley have used these methods to design the high subsonic, wing-mean-camber shapes for various configurations such as a supersonic transport, high-aspect-ratio transport, trapezoidal fighter wing, strake wing, tandem wing, joined wing, delta wing, and slender cranked wing. Many of these have been built, tested, and have had their data compared with theory. In addition, a technique for defining efficiently performing strake planforms for use in strake-wing combinations is discussed, and further improvements in wing design are outlined. The latter may be obtained by using higher-ordered linear panel methods as well as nonlinear-transonic methods.
This paper describes an efficient numerical approach for the design and optimization of aerodynamic bodies. As in classical optimal control methods, the present approach introduces a cost function and a costate variable (Lagrange multiplier) in order to achieve a minimum. High efficiency is achieved by using a multigrid technique to solve for all the unknowns simultaneously, but restricting work on a design variable only to grids on which their changes produce nonsmooth perturbations. Thus, the effort required to evaluate design variables that have nonlocal effects on the solution is confined to the coarse grids. However, if a variable has a nonsmooth local effect on the solution in some neighborhood, it is relaxed in that neighborhood on finer grids. The cost of solving the optimal control problem is shown to be approximately two to three times the cost of the equivalent analysis problem. Examples are presented to illustrate the application of the method to aerodynamic design and constraint optimization.
Progress in the direct-inverse wing design method in curvilinear coordinates has been made. This includes the remedying of a spanwise oscillation problem and the assessment of grid skewness, viscous interaction, and the initial airfoil section on the final design. It was found that, in response to the spanwise oscillation problem that designing at every other spanwise station produced the best results for the cases presented, a smoothly varying grid is especially needed for the accurate design at the wing tip, the boundary layer displacement thicknesses must be included in a successful wing design, the design of high and medium aspect ratio wings is possible with this code, and the final airfoil section designed is fairly independent of the initial section.
This paper describes a formal optimization procedure for helicopter rotor blade designs which minimizes hover horsepower while assuring satisfactory forward flight performance. The approach is to couple hover and forward flight analysis programs with a general purpose optimization procedure. The resulting optimization system provides a systematic evaluation of the rotor blade design variables and their interaction, thus reducing the time and cost of designing advanced rotor blades. The paper discusses the basis for and details of the overall procedure, describes the generation of advanced blade designs for representative Army helicopters, and compares designs and design effort with those from the conventional approach which is based on parametric studies and extensive cross-plots.
Described is a formal optimization procedure for helicopter rotor blade design which minimizes hover horsepower while assuring satisfactory forward flight performance. The approach is to couple hover and forward flight analysis programs with a general-purpose optimization procedure. The resulting optimization system provides a systematic evaluation of the rotor blade design variables and their interaction, thus reducing the time and cost of designing advanced rotor blades. The paper discusses the basis for and details of the overall procedure, describes the generation of advanced blade designs for representative Army helicopters, and compares design and design effort with those from the conventional approach which is based on parametric studies and extensive cross-plots.
A new inverse method for aerodynamic design of subcritical airfoils is presented. The pressure distribution in this method can be prescribed in a natural way, i.e. as a function of arclength of the as yet unknown body. This inverse problem is shown to be mathematically equivalent to solving a single nonlinear boundary value problem subject to known Dirichlet data on the boundary. The solution to this problem determines the airfoil, the free stream Mach number M(sub x) and the upstream flow direction theta(sub x). The existence of a solution for any given pressure distribution is discussed. The method is easy to implement and extremely efficient. We present a series of results for which comparisons are made with the known airfoils.
A numerical method based on the conservation form of the full potential equation has been applied to the problem of three-dimensional supersonic flows with embedded subsonic regions. The governing equation is cast in a nonorthogonal coordinate system, and the theory of characteristics is used to accurately monitor the type-dependent flow field. A conservative switching scheme is employed to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. The newly developed computer program can handle arbitrary geometries with fuselage, canard, wing, flow through nacelle, vertical tail and wake components at combined angles of attack and sideslip. Results are obtained for a variety of configurations that include a Langley advanced fighter concept with fuselage centerline nacelle, Rockwell's Advanced Tactical Fighter (ATF) with wing mounted nacelles, and the Shuttle Orbiter configuration. Comparisons with available experiments were good.
Progress in the development of computational methods for steady and unsteady aerodynamics has perennially paced advancements in aeroelastic analysis and design capabilities. Since these capabilities are of growing importance in the analysis and design of high-performance aircraft, considerable effort has been directed toward the development of appropriate aerodynamic methodology. The contributions to those efforts from the integral-equations research program at the NASA Langley Research Center is reviewed. Specifically, the current scope, progress, and plans for research and development for inviscid and viscous flows are discussed, and example applications are shown in order to highlight the generality, versatility, and attractive features of this methodology.
Since the project is rapidly nearing conclusion, the status of the tasks outlined in the original proposal are briefly outlined. These tasks include: viscous interation and wake curvature effects; code optimization and design methodology studies; methods for the design of isolated regions; program improvement efforts; and validation, testing, and documentation.