Gradient based aerodynamic shape optimization using the FIVER embedded boundary method

Aerodynamic shape optimisation technology is presented, comprising an efficient variable fidelity shape parameterisation method, an efficient andhigh quality mesh deformation scheme, and a parallel optimisation algorithm. The objective of the research presented here is the comparison of truly three-dimensional optimisations of aircraft wings in both aerodynamic and aeroelastic environments. The novel shape parameterisation technique allows various fidelities of design parameters, ranging from detailed surface changes to novel truly three-dimensional planform adjustments. An efficient interpolation scheme, using radial basis functions, transfers domain element movements into direct deformations of the design surface and corresponding CFD mesh, thus allowing total independence from the grid generation package and type (structured or unstructured). Optimisation is independent from the CFD flow solver by obtaining sensitivity information for an advanced parallel gradient-based optimiser by finite-differences. This ‘wrap-around’ optimisation technique is applied to a modern large transport aircraft wing in the cruise flight condition for minimum drag with stringent constraints in lift, volume, and two root moments. The objective of all optimisations is aerodynamic, however the static aeroelastic deflection provided by an aeroelastic solver will give that particular optimisation increased accuracy and real world relevance. The result of a constrained inviscid aerodynamic optimisation is presented and has a significant reduction in drag when compared to the initial wing with no violation of any constraints. The shape parameterisation method demonstrates that only a low number of design variables are necessary to achieve innovative planform and surface geometries with dramatically improved performance. Computational fluid dynamics (CFD) methods are now commonplace in aerospace industries, and at the forefront of analysis capabilities, providing a fast and effective method of predicting a design’s aerodynamic performance. However, with ever increasing complexity of designs, engineers can often struggle to interpret the intricacies of the CFD results sufficiently to be able to manually alter the geometry to improve performance. Hence, there has been an increase in demand for intelligent and automatic shape optimisation schemes. This requires combining geometry control methods with numerical optimisation algorithms, to provide a mechanism to mathematically seek improved and optimum designs, using CFD as the analysis tool. Optimisation requires consideration of three issues, each of which have numerous solutions: shape parameterisation including CFD surface and volume mesh deformation, computation of design variable derivatives, and effective use of these derivatives to improve design. Geometry parameterisation is critical for effective shape optimisation. This is the method of representing the design surface, and defines the degrees of freedom in which the geometry can be altered and, ideally, this should be linked with an effective method of deforming the CFD surface and volume mesh in a corresponding fashion. Parameterising complex shapes is a problem that remains a serious obstacle to both manual and automatic CFD-based optimisation. A wide variety of shape control and morphing methods have been developed, but many do not allow sufficiently

This paper presents an aerodynamic shape optimization method based on a genetic algorithm with the compressible Euler or Navier-Stokes equations. A new algorithm to represent aerodynamic configurations is proposed to find optimum aerodynamic configurations efficiently. We solve three aerodynamic shape optimization problems. First, an inverse problem is solved to understand a basic character and to examine the possibility of the genetic algorithm. Second, a drag minimization problem is solved to find a streamlined body. Finally, L/D(ratio of lift to drag) maximization problem in a transonic flow is solved and the optimized configuration is compared with a supercritical airfoil. The results show that genetic algorithm is powerful and robust in solving aerodynamic shape optimization problems.

An efficient approach to optimization of real-life 3D aerodynamic shapes for minimum drag is presented. The method allows for substitution of fictitious surfaces for those parts of aircraft which are not subject to modification in the process of design. The approach essentially reduces the geometrical complexity of computational models thus making the optimization of complete aerodynamic configurations practically feasible and, eventually, increasing the accuracy and reliability of optimization. The optimization framework is that of OPTIMAS, the in-house IAI optimization tool based on Genetic Algorithms and full Navier-Stokes computations. The method is illustrated by example of wing shape optimization for the full configuration of generic business jet with close coupled wing-body-nacelle, and vertical/horisontal tails. The results indicate the applicability of the method to practical aerodynamic design.KeywordsAerodynamic PerformanceWing ShapeAerodynamic DesignAerodynamic ShapeMinimum DragThese keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

Optimization problems in many industrial applications are very hard to solve. Many examples of them can be found in the design of aeronautical systems. In this field, the designer is frequently faced with the problem of considering not only a single design objective, but several of them, i.e., the designer needs to solve a Multi-Objective Optimization Problem (MOP). In aeronautical systems design, aerodynamics plays a key role in aircraft design, as well as in the design of propulsion system components, such as turbine engines. Thus, aerodynamic shape optimization is a crucial task, and has been extensively studied and developed. Multi-Objective Evolutionary Algorithms (MOEAs) have gained popularity in recent years as optimization methods in this area, mainly because of their simplicity, their ease of use and their suitability to be coupled to specialized numerical simulation tools. In this chapter, we will review some of the most relevant research on the use of MOEAs to solve multi-objective and/or multi-disciplinary aerodynamic shape optimization problems. In this review, we will highlight some of the benefits and drawbacks of the use of MOEAs, as compared to traditional design optimization methods. In the second part of the chapter, we will present a case study on the application of MOEAs for the solution of a multi-objective aerodynamic shape optimization problem.

Chapter 4 - Massive Parallelization as Prinicipal Technology for Constrained Optimization of Aerodynamic Shapes

Aerodynamic shape optimization of aircraft configurations often ignores stability considerations. To address this, a method for the computation of static, dynamic, and transient aircraft stability derivatives and their sensitivities for use in gradient-based optimization is introduced and evaluated. Computational fluid dynamics in the form of a three-dimensional structured-grid multiblock flow solver with both Euler and Reynolds-averaged Navier–Stokes equations is used. To compute the stability derivatives, a time-spectral formulation is used to compute an oscillating solution for the configuration of interest. From this oscillating solution, a series of linear regressions is performed to calculate the various stability derivatives. Because the solution is time dependent, it contains the information required to compute the transient, or “dot,” derivatives for the configuration. An adjoint method is used to compute the gradients of the stability derivatives of interest, enabling gradient-based stability-constrained aerodynamic shape optimization with respect to a large number of design variables. The computed stability derivatives are verified for an airfoil and validated for a generic unmanned combat aerial vehicle. The stability-constrained optimization of a wing demonstrates the viability and usefulness of the method for aircraft design optimization.

We present a new continuous adjoint method for aerodynamic shape optimization using the Euler equations, which reduces the computational cost of the gradients by reducing the volume integral part of the adjoint gradient formula to a surface integral. The savings are particularly significant for three-dimensional aerodynamic shape optimization problems on general unstructured and overset meshes. In order to validate the concept, the new gradient equations have been tested for various aerodynamic shape optimization problems, including an inverse problem for threedimensional wing configurations, and drag minimization problems of a single-element airfoil and a three-dimensional wing-fuselage configuration. In order to assess their accuracy, the results are compared with finite-difference gradients, complex-step gradients, and gradients calculated by the previous adjoint method which includes a volume integral. ∗Thomas V. Jones Professor of Engineering, Department of Aeronautics and Astronautics, Stanford University, AIAA Fellow †Postdoctoral Fellow, Department of Aeronautics and Astronautics, Stanford University, AIAA Member

ABSTRACTIn this article, a Virtual Stackelberg Game (VSG) is proposed for aerodynamic shape optimization, where the design variables are divided into two categories to optimize the same objective function, one acts as a leader, and the other ones as followers react independently and selfishly relative to the leader's strategy. During each Stackelberg strategy cycle, the Gradient-Based Method (GBM) with the adjoint method in Stanford University Unstructured (SU2) is applied in the optimization of each player. Firstly, parametric studies of VSG by two simple cases are conducted to assess the impact of critical parameters on aerodynamic shape optimization, including the design cycle, the split of design variables and role (leader and follower) assignment. Based on the criterion from parametric studies, two typical numerical cases under transonic flow are applied—the drag reduction design of a 2D airfoil and a 3D wing. It is found that, compared to the original GBM method, VSG can provide better optimization results with less computational cost.

Nowadays, one of the priorities of the European Commission is to reduce the environmental impact of aviation through the advanced design of novel aircraft configurations. This is of utmost importance in order to decrease the environmental footprint of aviation and to reduce fuel consumption and make airlines more profitable. This implies that new methods and tools for aerodynamic shape optimization will have to be developed, allowing aircraft configurations that cannot be obtained with traditional strategies. This paper focuses on the application of enhanced methods in aerodynamic shape design optimization to enable advanced aircraft configurations. In particular, this work aims to demonstrate the feasibility of the proposed strategy to reach optimal configurations that are far away from its baseline geometry. For this purpose, evolutionary algorithms are combined with support vector machines and applied to the optimization of a baseline geometry for different flow conditions. In particular, the selected application is based on the shape optimization problem of the landing gear master cylinder. Results pointed out the feasibility of the mentioned strategy to enable novel configurations within an aerodynamic shape optimization process.

A new method for aerodynamic shape optimization using a genetic algorithm with real number encoding is presented. The algorithm is used to optimize three different problems, a simple hill climbing problem, a quasi-one-dimensional nozzle problem using an Euler equation solver and a three-dimensional transonic wing problem using a nonlinear potential solver. Results indicate that the genetic algorithm is easy to implement and extremely reliable, being relatively insensitive to design space noise.

This paper presents an assessment of the potential fuel burn savings offered by the transonic strut-braced-wing configuration within the single-aisle class of aircraft relative to a modern conventional tube-and-wing aircraft through aerodynamic shape optimization based on the Reynolds-averaged Navier–Stokes equations. A representative strut-braced-wing aircraft is first developed through conceptual multidisciplinary design optimization based on the Airbus A320neo, with current technology levels assumed. A concept for the conventional tube-and-wing configuration is also developed to represent the Airbus A320neo as a performance baseline. Single-point aerodynamic shape optimization is then performed on wing–body–tail models of each aircraft to address aerodynamic design challenges and to provide more accurate performance estimates. Results indicate that shock formation can be mitigated from the wing–strut junction of the strut-braced wing at Mach 0.78 and a relatively high design lift coefficient of 0.750, providing an 8.2% reduction in block fuel over a 1000 n mile nominal mission when compared to the conventional tube-and-wing aircraft. Multipoint aerodynamic shape optimization is then performed to build toward a more credible estimate of fuel burn performance, with results showing a reduction in the fuel burn savings to 7.8% at the nominal design point relative to the conventional tube-and-wing aircraft to maintain a 7.6–8.0% improvement over the envelope of operating conditions, which includes design points at even higher Mach numbers and lift coefficients. These results demonstrate the viability of the transonic strut-braced-wing configuration for transport aircraft within the single-aisle class and its potential for reducing commercial fleet fuel burn.

In the aerodynamic shape design, the drag prediction has always been an extremely challenging mission for the exploration of a configuration. As for the more complex configurations, it is especially desired to the availability of a highly accurate and reliable aerodynamic numerical solution. For improving the drag prediction accuracy and promoting the aerodynamic shape designs, firstly, the characteristics of drag prediction based on far-field drag method and near-field drag method is analyzed and compared. Also, the merits and demerits of defining axial velocity defect with the current main far-field drag prediction approaches is summarized, which promotes the building of the improved method of axial velocity defect and the improved far-field drag prediction and decomposition approach. Moreover, during the establishment of the drag decomposition method, it is necessary to judge and decide on the selection of the drag region. Therefore, the discussions on the sensitivity of the relevant parameters are fulfilled. Furthermore, based on the far-field drag prediction and decomposition method constructed, the aerodynamic performance research of Common Research Model wing-body configuration is launched. The results show that it can effectively observe and analyze the changes in drag components, their impact on the total drag and the contribution percentage. Finally, combining the far-field drag prediction and decomposition method proposed in this paper with a gradient-based aerodynamic shape optimization design system, the aerodynamic shape optimization designs are studied with CRM wing-body configuration. The results can not only directly analyze the detailed change of the visualized drag region, but also can obtain the more accurate total drag and lift-to-drag ratio of the optimized configuration by removing the spurious drag.

The traditional multi-objective efficient global optimization (EGO) algorithms have been hybridized and adapted to solving the expensive aerodynamic shape optimization problems based on high-fidelity numerical simulations. Although the traditional EGO algorithms are highly efficient in solving some of the optimization problems with very complex landscape, it is not preferred to solve most of the aerodynamic shape optimization problems with relatively low-degree multi-modal design spaces. A new infill criterion encouraging more local exploitation has been proposed by hybridizing two traditional multi-objective expected improvements (EIs), namely, statistical multi-objective EI and expected hypervolume improvement, in order to improve their robustness and efficiency in aerodynamic shape optimization. Different analytical test problems and aerodynamic shape optimization problems have been investigated. In comparison with traditional multi-objective EI algorithms and a standard evolutionary multi-objective optimization algorithm, the proposed method is shown to be more robust and efficient in the tests due to its hybrid characteristics, easier handling of sub-optimization problems, and enhanced exploitation capability.

This paper presents a gradient-based aerodynamic shape optimization framework that utilizes a Graphics Processing Unit (GPU) for solving both flow and adjoint equations. It is built based on a GPU-accelerated flow solver that has been developed previously. Hence, the focus of this work is on how to solve the adjoint equations on the GPU and subsequently compute the gradients. The adjoint equations are right-preconditioned by a block Incomplete Lower Upper (ILU) preconditioner and solved by a restarted Generalized Minimum Residual (GMRES) method. The exact residual Jacobian matrix in the adjoint equations is computed using finite difference and a distance-2 graph coloring algorithm. With the adjoint-based gradients, the steepest descent method with momentum is employed for constrained aerodynamic shape optimization of a wing-body configuration at a transonic flow condition.

The present-day scenario of urbanization is mostly vertical densification that results in the building of super-tall structures in metropolitan areas with land shortage. In the process of going higher up, these buildings become part of the wind load rather than gravity acting on them. The combined action of turbulence wind flow and bluff bodies of high-rise buildings promotes the introduction of complex aerodynamic phenomena that include vortex shedding, buffeting, and galloping, all of which have a direct bearing on structural safety and occupant comfort. While traditional design techniques still focus on prismatic geometries, the recent architectural trends have moved to complex and irregular plan configurations, such as that of the Fish, Y, and Plus shapes, and also the globally twisted and tapered towers, which are mainly driven by aesthetic aspirations and functional optimization. Nevertheless, the current international design codes such as ASCE 7-16 and IS 875 (Part 3) do not provide many guidelines on these non-standard geometries and, consequently, performance-based wind engineering testing is commonly needed (ASCE, 2017; Bureau of Indian Standards, 2015). This is a review of the aerodynamic performance of complex plan-geometry high-rise buildings critically synthesized based on the peer-reviewed literature. It is conducted in a methodical assessment against the effectiveness of the Aerodynamic Shape Optimization (ASO) strategies by local geometric alterations including corner rounding, chamfering and recession against global morphological alterations including twisting and tapering. Additionally, the review examine the current analytical and numerical methods, find the drawbacks of Reynolds-Averaged Navier-Stokes (RANS) models regarding their inability to capture unsteady wake dynamics in comparison to Large Eddy Simulation (LES) and Delayed Detached Eddy Simulation (DDES) methods. The combined results illustrate the extent to which corner rounding and helical twisting lower the drag forces and the susceptibility to vibration caused by vortexes through disrupting the span-wise coherence of the vortices, while highly irregular configurations, like Fish-shaped plans, are very sensitive to the direction of the wind. Moreover, the review concludes with the identification of critical research gaps concerning the wind interference effects in dense urban clusters, and it positions the possible adoption of machine-learning-based surrogate models in the preliminary design phases as an attractive option to solve the problem of wind design codes currently having limitations.