An evolutionary geometry parametrization is established to represent aerodynamic configurations. This geometry parametrization technique is constructed by integrating the classical B-spline formulation with the knot insertion algorithm. It is capable of inserting control points to a given parametrization without modifying its geometry. Taking advantage of this technique, a shape design problem can be solved as a sequence of optimizations from the basic parametrization to more refined parametrizations. Owing to the nature of the B-spline formulation, feasible parametrization refinements are not unique; guidelines based on sensitivity analysis and geometry constraints are developed to assist the automation of the proposed optimization sequence. Test cases involving airfoil optimization and induced drag minimization are solved adopting this method. Its effectiveness is demonstrated through comparisons with optimizations using uniform refined parametrizations.

An evolutionary geometry parametrization is established to represent aerodynamic configurations. This geometry parametrization technique is constructed by integrating the classical B-spline formulation with the knot insertion algorithm. It is capable of inserting control points to a given parametrization without modifying its geometry. Taking advantage of this technique, a shape design problem can be solved as a sequence of optimizations from the basic parametrization to more refined parametrizations. Owing to the nature of the B-spline formulation, feasible parametrization refinements are not unique; guidelines based on sensitivity analysis and geometry constraints are developed to assist the automation of the proposed optimization sequence. Test cases involving airfoil optimization and induced drag minimization are solved adopting this method. Its effectiveness is demonstrated through comparisons with optimizations using uniform refined parametrizations.

The first two steps in the development of an integrated multidisciplinary design optimization procedure capable of analyzing the nonlinear fluid flow about geometrically complex aeroelastic configurations have been accomplished in the present work. For the first step, a three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed. The advantage of unstructured grids, when compared with a structured-grid approach, is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited fro geometrically complex configurations of practical interest. In this work the time-dependent, nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional cases and a Gauss-Seidel algorithm for the three-dimensional; at steady-state, similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Various surface parameterization techniques have been employed in the current study to control the shape of the design surface. Once this surface has been deformed, the interior volume of the unstructured grid is adapted by considering the mesh as a system of interconnected tension springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR, an advanced automatic-differentiation software tool. To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for several two- and three-dimensional cases. In two-dimensions, an initially symmetric NACA-0012 airfoil and a high-lift multi-element airfoil were examined. For the three-dimensional configurations, an initially rectangular wing with uniform NACA-0012 cross-sections was optimized; in additions, a complete Boeing 747-200 aircraft was studied. Furthermore, the current study also examines the effect of inconsistency in the order of spatial accuracy between the nonlinear fluid and linear shape sensitivity equations. The second step was to develop a computationally efficient, high-fidelity, integrated static aeroelastic analysis procedure. To accomplish this, a structural analysis code was coupled with the aforementioned unstructured grid aerodynamic analysis solver. The use of an unstructured grid scheme for the aerodynamic analysis enhances the interactions compatibility with the wing structure. The structural analysis utilizes finite elements to model the wing so that accurate structural deflections may be obtained. In the current work, parameters have been introduced to control the interaction of the computational fluid dynamics and structural analyses; these control parameters permit extremely efficient static aeroelastic computations. To demonstrate and evaluate this procedure, static aeroelastic analysis results for a flexible wing in low subsonic, high subsonic (subcritical), transonic (supercritical), and supersonic flow conditions are presented. / Ph. D.

Increasing environmental concerns and fuel prices motivate the study of alternative, unconventional aircraft configurations. One such example is the blended-wing-body configuration, which has been shown to have several advantages over the conventional tube-and-wing aircraft configuration. In this thesis, a blended-wing-body aircraft is studied and optimized aerodynamically using a high-fidelity Euler-based flow solver, integrated geometry parameterization and mesh movement, adjoint-based gradient evaluation, and a sequential quadratic programming algorithm. Specifically, the aircraft is optimized at transonic conditions to minimize the sum of induced and wave drag. These optimizations are carried out with both fixed and varying airfoil sections. With varying airfoil sections and increased freedom, up to 52% drag reduction relative to the baseline geometry was achieved: at the target lift coefficient of 0.357, a drag coefficient of 0.01313 and an inviscid lift-to-drag ratio of 27.2 were obtained.

Increasing environmental concerns and fuel prices motivate the study of alternative, unconventional aircraft configurations. One such example is the blended-wing-body configuration, which has been shown to have several advantages over the conventional tube-and-wing aircraft configuration. In this thesis, a blended-wing-body aircraft is studied and optimized aerodynamically using a high-fidelity Euler-based flow solver, integrated geometry parameterization and mesh movement, adjoint-based gradient evaluation, and a sequential quadratic programming algorithm. Specifically, the aircraft is optimized at transonic conditions to minimize the sum of induced and wave drag. These optimizations are carried out with both fixed and varying airfoil sections. With varying airfoil sections and increased freedom, up to 52% drag reduction relative to the baseline geometry was achieved: at the target lift coefficient of 0.357, a drag coefficient of 0.01313 and an inviscid lift-to-drag ratio of 27.2 were obtained.

A set of techniques is developed that allows the incorporation of flight dynamics metrics as an additional discipline in a high-fidelity aerodynamic optimization. Specifically, techniques for including static stability constraints and handling qualities constraints in a high-fidelity aerodynamic optimization are demonstrated. These constraints are developed from stability derivative information calculated using high-fidelity computational fluid dynamics (CFD). Two techniques are explored for computing the stability derivatives from CFD. One technique uses an automatic differentiation adjoint technique (ADjoint) to efficiently and accurately compute a full set of static and dynamic stability derivatives from a single steady solution. The other technique uses a linear regression method to compute the stability derivatives from a quasi-unsteady time-spectral CFD solution, allowing for the computation of static, dynamic and transient stability derivatives. Based on the characteristics of the two methods, the time-spectral technique is selected for further development, incorporated into an optimization framework, and used to conduct stability-constrained aerodynamic optimization. This stability-constrained optimization framework is then used to conduct an optimization study of a flying wing configuration. This study shows that stability constraints have a significant impact on the optimal design of flying wings and that, while static stability constraints can often be satisfied by modifying the airfoil profiles of the wing, dynamic stability constraints can require a significant change in the planform of the aircraft in order for the constraints to be satisfied.

A set of techniques is developed that allows the incorporation of flight dynamics metrics as an additional discipline in a high-fidelity aerodynamic optimization. Specifically, techniques for including static stability constraints and handling qualities constraints in a high-fidelity aerodynamic optimization are demonstrated. These constraints are developed from stability derivative information calculated using high-fidelity computational fluid dynamics (CFD). Two techniques are explored for computing the stability derivatives from CFD. One technique uses an automatic differentiation adjoint technique (ADjoint) to efficiently and accurately compute a full set of static and dynamic stability derivatives from a single steady solution. The other technique uses a linear regression method to compute the stability derivatives from a quasi-unsteady time-spectral CFD solution, allowing for the computation of static, dynamic and transient stability derivatives. Based on the characteristics of the two methods, the time-spectral technique is selected for further development, incorporated into an optimization framework, and used to conduct stability-constrained aerodynamic optimization. This stability-constrained optimization framework is then used to conduct an optimization study of a flying wing configuration. This study shows that stability constraints have a significant impact on the optimal design of flying wings and that, while static stability constraints can often be satisfied by modifying the airfoil profiles of the wing, dynamic stability constraints can require a significant change in the planform of the aircraft in order for the constraints to be satisfied.

One key outcome of this research is the design for a 3-D tractor-trailer base-drag reduction device that predicts a 41% reduction in wind-averaged drag coefficient at 57 mph (92 km/h) and that is relatively insensitive to uncertain wind speed and direction and uncertain deflection angles due to mounting accuracy and static aeroelastic loading; the best commercial device of non-optimized design achieves a 12% reduction at 65 mph. Another important outcome is the process by which the optimized design is obtained. That process includes verification and validation of the flow solver, a less complex but much broader 2-D pathfinder study, and the culminating 3-D aerodynamic shape optimization under uncertainty (OUU) study. To gain confidence in the accuracy and precision of a computational fluid dynamics (CFD) flow solver and its Reynolds-averaged Navier-Stokes (RANS) turbulence models, it is necessary to conduct code verification, solution verification, and model validation. These activities are accomplished using two commercial CFD solvers, Cobalt and RavenCFD, with four turbulence models: Spalart-Allmaras (S-A), S-A with rotation and curvature, Menter shear-stress transport (SST), and Wilcox 1998 k-ω. Model performance is evaluated for three low subsonic 2-D applications: turbulent flat plate, planar jet, and NACA 0012 airfoil at α = 0°. The S-A turbulence model is selected for the 2-D OUU study. In the 2-D study, a tractor-trailer base flap model is developed that includes six design variables with generous constraints; 400 design candidates are evaluated. The design optimization loop includes the effect of uncertain wind speed and direction, and post processing addresses several other uncertain effects on drag prediction. The study compares the efficiency and accuracy of two optimization algorithms, evolutionary algorithm (EA) and dividing rectangles (DIRECT), twelve surrogate models, six sampling methods, and surrogate-based global optimization (SBGO) methods. The DAKOTA optimization and uncertainty quantification framework is used to interface the RANS flow solver, grid generator, and optimization algorithm. The EA is determined to be more efficient in obtaining a design with significantly reduced drag (as opposed to more efficient in finding the true drag minimum), and total predictive uncertainty is estimated as ±11%. While the SBGO methods are more efficient than a traditional optimization algorithm, they are computationally inefficient due to their serial nature, as implemented in DAKOTA. Because the S-A model does well in 2-D but not in 3-D under these conditions, the SST turbulence model is selected for the 3-D OUU study that includes five design variables and evaluates a total of 130 design candidates. Again using the EA, the study propagates aleatory (wind speed and direction) and epistemic (perturbations in flap deflection angle) uncertainty within the optimization loop and post processes several other uncertain effects. For the best 3-D design, total predictive uncertainty is +15/-42%, due largely to using a relatively coarse (six million cell) grid. That is, the best design drag coefficient estimate is within 15 and 42% of the true value; however, its improvement relative to the no-flaps baseline is accurate within 3-9% uncertainty. / Ph. D.

Strides in modern computational fluid dynamics and leaps in high-power computing have led to unprecedented capabilities for handling large aerodynamic problem. In particular, the emergence of adjoint design methods has been a break-through in the field of aerodynamic shape optimization. It enables expensive, high-dimensional optimization problems to be tackled efficiently using gradient-based methods in CFD; a task that was previously inconceivable. However, adjoint design methods are intended for gradient-based optimization; the curse of dimensionality is still very much alive when it comes to design space exploration, where gradient-free methods cannot be avoided. This research describes a novel approach for reducing dimensionality in large, computationally expensive design problems to a point where gradient-free methods become possible. This is done using an innovative application of Principal Component Analysis (PCA), where the latter is applied to the gradient distribution of the objective function; something that had not been done before. This yields a linear transformation that maps a high-dimensional problem onto an equivalent low-dimensional subspace. None of the original variables are discarded; they are simply linearly combined into a new set of variables that are fewer in number. The method is tested on a range of analytical functions, a two-dimensional staggered airfoil test problem and a three-dimensional Over-Wing Nacelle (OWN) integration problem. In all cases, the method performed as expected and was found to be cost effective, requiring only a relatively small number of samples to achieve large dimensionality reduction.

The development of a shape optimization solver using the existing Computational Fluid Dynamics (CFD) codes is taken up as topic of research in this thesis. A shape optimizer was initially developed based on Genetic Algorithm (GA) coupled with a CFD solver in an earlier work. The existing CFD solver is based on Kinetic Flux Vector Splitting and uses least squares discretization. This solver requires a cloud of points and their connectivity set, hence this CFD solver is a meshless solver. The advantage of a meshless solver is utilised in avoiding re-gridding (only connectivity regeneration is required) after each shape change by the shape optimizer. The CFD solver is within the optimization loop, hence evaluation of CFD solver after each shape change is mandatory. Although the earlier shape optimizer developed was found to be robust, but it was taking enoromous amount of time to converge to the optimum solution (details in Appendix). Hence a new evolving method, Ant Colony Optimization (ACO), is implemented to replace GA. A shape optimizer is developed coupling ACO and the meshless CFD solver. To the best of the knowledge of the present author, this is the first time when ACO is implemented for aerodynamic shape optimization problems. Hence, an exhaustive validation has become mandatory. Various test cases such as regeneration problems of (1) subsonic - supersonic nozzle with a shock in quasi - one dimensional flow (2) subsonic - supersonic nozzle in a 2-dimensional flow field (3) NACA 0012 airfoil in 2-dimensional flow and (4) NACA 4412 airfoil in 2-dimensional flow have been successfully demonstrated. A comparative study between GA and ACO at algorithm level is performed using the travelling salesman problem (TSP). A comparative study between the two shape optimizers developed, i.e., GA-CFD and ACO-CFD is carried out using regeneration test case of NACA 4412 airfoil in 2-dimensional flow. GA-CFD performs better in the initial phase of optimization and ACO-CFD performs better in the later stage. We have combined both the approaches to develop a hybrid GA-ACO-CFD solver such that the advantages of both GA-CFD and ACO-CFD are retained with the hybrid method. This hybrid approach has 2 stages, namely, (Stage 1) initial optimum search by GA-CFD (coarse search), the best members from the optimized solution from GA-CFD are segregated to form the input for the fine search by ACO-CFD and (Stage 2) final optimum search by ACO-CFD (fine search). It is observed that this hybrid method performs better than either GA-CFD or ACO- CFD, i.e., hybrid method attains better optimum in less number of CFD calls. This hybrid method is applied to the following test cases: (1) regeneration of subsonic-supersonic nozzle with shock in quasi 1-D flow and (2) regeneration of NACA 4412 airfoil in 2-dimensional flow. Two applications on shape optimization, namely, (1) shape optimization of a body in strongly rotating viscous flow and (2) shape optimization of a body in supersonic flow such that it enhances separation of binary species, have been successfully demonstrated using the hybrid GA-ACO-CFD method. A KFVS based binary diffusion solver was developed and validated for this purpose. This hybrid method is now in a state where industrial shape optimization applications can be handled confidently.