Geometric perspectives on MDO and MDO architectures

Two main perspectives have been developed within the Multidisciplinary Design Optimization (MDO) literature for classifying and comparing MDO architectures: a numerical point of view and a formulation/data flow point of view. Although significant work has been done here, these perspectives have not provided much in the way of a priori information or predictive power about architecture performance. In this report, we outline a new perspective, called the geometric perspective, which we believe will be able to provide such predictive power. Using tools from differential geometry, we take several prominent architectures and describe mathematically how each constructs the space through which it moves. We then consider how the architecture moves through the space which it has constructed. Taken together, these investigations show how each architecture relates to the original feasible design manifold, how the architectures relate to each other, and how each architecture deals with the design coupling inherent to the original system. This in turn lays the groundwork for further theoretical comparisons between and analyses of MDO architectures and their behaviour using tools and techniques derived from differential geometry. © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

On the application of differential geometry to MDO

Multidisciplinary Design Optimization (MDO) is a methodology for optimizing large coupled systems. Over the years, a number of different MDO decomposition strategies, known as architectures, have been developed, and various pieces of analytical work have been done on MDO and its architectures. However, MDO lacks an overarching paradigm which would unify the field and promote cumulative research. In this paper, we propose a differential geometry framework as such a paradigm: Differential geometry comes with its own set of analysis tools and a long history of use in theoretical physics. We begin by outlining some of the mathematics behind differential geometry and then translate MDO into that framework. This initial work gives new tools and techniques for studying MDO and its architectures while producing a naturally arising measure of design coupling. The framework also suggests several new areas for exploration into and analysis of MDO systems. At this point, analogies with particle dynamics and systems of differential equations look particularly promising for both the wealth of extant background theory that they have and the potential predictive and evaluative power that they hold. © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Acoustic analogy in swirling mean flow applied to predict rotor trailing-edge noise

The present study aims at accounting for swirling mean flow effects on rotor trailing-edge noise. Indeed, the mean flow in between the rotor and the stator of the fan or of a compressor stage is highly swirling. The extension of Ffowcs-Williams & Hawkings' acoustic analogy in a medium at rest with moving surfaces and of Goldstein's acoustic analogy in a circular duct with uniform mean flow to a swirling mean flow in an annular duct is introduced. It is first applied to tonal noise. In most cases, the swirl modifies the pressure distribution downstream of the fan. In several configurations, when the swirl is rather close to a solid body swirl, it is often sufficient to apply a simple Doppler effect correction when predicting the duct modes in uniform mean flow in order to predict accurately the noise radiated with swirl. However, in other realistic configurations, the swirling mean-flow effect cannot be addressed using this simple Doppler effect correction. Second, a rotor trailing-edge noise model accounting for both the effects of the annular duct and the swirling mean flow is developed and applied to a realistic fan rotor with different swirling and sheared mean flows (and as a result different associated blade stagger angles). The benchmark cases are built from the Boeing 18-inch Fan Rig Broadband Noise Test. In all cases the swirling mean flow has an effect. In some cases the a simple Doppler effect may address it, but, in other realistic configurations our acoustic analogy with swirl is needed. © 2012 by the authors. Published by the American Institute of Aeronautics and Astronautics, Inc.

A canonical normal SBLI flow relevant to external compression inlets

The normal shock wave / boundary layer interaction (normal SBLI) is important to the operation and performance of a supersonic inlet, and the normal SBLI is particularly prominent in external compression inlets. To improve our understanding of such interactions, it is helpful to make use of fundamental flows which capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental fiow-fleld configurations have been considered as possible test cases to represent the normal SBLI aspects found in typical external compression inlets, and it was found that the spillage-diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flow-fleld allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to be all held approximately constant mid independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage-diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, the width around twice or three times the height, and with the area expansion just downstream of the shock on the conservative side of the stall limit for incompressible diffusers. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Enhancing and developing the practical optimisation capabilities and intelligence of automatic design software

In the modern engineering design cycle the use of computational tools becomes a neces- sity. The complexity of the engineering systems under consideration for design increases dramatically as the demands for advanced and innovative design concepts and engineering products is expanding. At the same time the advancements in the available technology in terms of computational resources and power, as well as the intelligence of the design software, accommodate these demands and make them a viable approach towards the chal- lenge of real-world engineering problems. This class of design optimisation problems is by nature multi-disciplinary. In the present work we establish enhanced optimisation capabil- ities within the Nimrod/O tool for massively distributed execution of computational tasks through cluster and computational grid resources, and develop the potential to combine and benefit from all the possible available technological advancements, both software and hardware. We develop the interface between a Free Form Deformation geometry manage- ment in-house code with the 2D airfoil aerodynamic efficiency evaluation tool XFoil, and the well established multi-objective heuristic optimisation algorithm NSGA-II. A simple airfoil design problem has been defined to demonstrate the functionality of the design sys- tem, but also to accommodate a framework for future developments and testing with other state-of-the-art optimisation algorithms such as the Multi-Objective Genetic Algorithm (MOGA) and the Multi-Objective Tabu Search (MOTS) techniques. Ultimately, heav- ily computationally expensive industrial design cases can be realised within the presented framework that could not be investigated before. © 2012 by the authors. Published by the American Institute of Aeronautics and Astronautics, Inc.

A computational fluid-structure interaction study of inflation cycles of a tension cone decelerator

Inflatable aerodynamic decelerators present potential advantages for planetary entry in missions of robotic and human exploration. The design of these structures face many engineering challenges, including complex deformable geometries, anisotropic material response, and coupled shockturbulence interactions. In this paper, we describe a comprehensive computational fluid-structure interaction study of an inflation cycle of a tension cone decelerator in supersonic flow and compare the simulations with earlier published experimental results. The aeroshell design and flow conditions closely match recent experiments conducted at Mach 2.5. The structural model is a 16-sided polygonal tension cone with seams between each segment. The computational model utilizes adaptive mesh refinement, large-eddy simulation, and shell mechanics with self-contact modeling to represent the flow and structure interaction. This study focuses on the dynamics of the structure as the inflation pressure varies gradually, and the behavior of forces experienced by the flexible and rigid (the payload capsule) structures. © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Manufacturing dish shaped rings on radial-axial ring rolling mills

Ring rolling is an established method to produce seamless rings of different cross-sectional geometries. For dish shaped rings, there are applications in different areas such as offshore, aeronautics or the energy sector. At the moment, dish shaped rings are produced by machining of rings with rectangular shaped cross section, by (open die) hollow forging on a conical mandrel or by using shaped ring rolling tools. These ways of manufacturing have the disadvantage of high material waste, additional costs for special tools, long process time and limited or inflexible geometries. Therefore, the manufacturing of dish shaped rings on conventional radial-axial ring rolling mills would expand the range of products for ring producers. The aim of this study is to investigate the feasibility of an alternative to the current manufacturing processes, without requiring additional tooling and material costs. Therefore, the intended formation of dish shaped rings-previously regarded as a form error-is investigated. Based on an analysis of geometrical requirements and metal flow mechanisms, a rolling strategy is presented, causing dishing and ring climbing by a large height reduction of the ring. Using this rolling strategy dish shaped rings with dishing angles up to 18° were achieved. In addition to the experiments finite element method (FEM)-simulations of the process have been successfully conducted, in order to analyze the local strain evolution. However, when the contact between ring and main roll is lost in the process the ring starts to oscillate around the mandrel and neither dishing nor ring climbing is observed. © 2013 German Academic Society for Production Engineering (WGP).

Multispacecraft refueling optimization considering the J2 perturbation and window constraints

The optimization of a near-circular low-Earth-orbit multispacecraft refueling problem is studied. The refueling sequence, service time, and orbital transfer time are used as design variables, whereas the mean mission completion time and mean propellant consumed by orbital maneuvers are used as design objectives. The J2 term of the Earth's nonspherical gravity perturbation and the constraints of rendezvous time windows are taken into account. A hybridencoding genetic algorithm, which uses normal fitness assignment to find the minimum mean propellant-cost solution and fitness assignment based on the concept of Pareto-optimality to find multi-objective optimal solutions, is presented. The proposed approach is demonstrated for a typical multispacecraft refueling problem. The results show that the proposed approach is effective, and that the J2 perturbation and the time-window constraints have considerable influences on the optimization results. For the problems in which the J2 perturbation is not accounted for, the optimal refueling order can be simply determined as a sequential order or as the order only based on orbitalplane differences. In contrast, for the problems that do consider the J2 perturbation, the optimal solutions obtained have a variety of refueling orders and use the drift of nodes effectively to reduce the propellant cost for eliminating orbital-plane differences. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Joint experimental and numerical approach to three-dimensional shock control bump research

Previous studies of transonic shock control bumps have often been either numerical or experimental. Comparisons between the two have been hampered by the limitations of either approach. The present work aims to bridge the gap between computational fluid dynamics and experiment by planning a joint approach from the outset. This enables high-quality validation data to be produced and ensures that the conclusions of either aspect of the study are directly relevant to the application. Experiments conducted with bumps mounted on the floor of a blowdown tunnel were modified to include an additional postshock adverse pressure gradient through the use of a diffuser as well as introducing boundary-layer suction ahead of the test section to enable the in-flow boundary layer to be manipulated. This has the advantage of being an inexpensive and highly repeatable method. Computations were performed on a standard airfoil model, with the flight conditions as free parameters. The experimental and computational setups were then tuned to produce baseline conditions that agree well, enabling confidence that the experimental conclusions are relevant. The methods are then applied to two different shock control bumps: a smoothly contoured bump, representative of previous studies, and a novel extended geometry featuring a continuously widening tail, which spans the wind-tunnel width at the rear of the bump. Comparison between the computational and experimental results for the contour bump showed good agreement both with respect to the flow structures and quantitative analysis of the boundary-layer parameters. It was seen that combining the experimental and numerical data could provide valuable insight into the flow physics, which would not generally be possible for a one-sided approach. The experiments and computational fluid dynamics were also seen to agree well for the extended bump geometry, providing evidence that, even though thebumpinteracts directly with the wind-tunnel walls, it was still possible to observe the key flow physics. The joint approach is thus suitable even for wider bump geometries. Copyright © 2013 by S. P. Colliss, H. Babinsky, K. Nubler, and T. Lutz. Published by the American Institute of Aeronautics and Astronautics, Inc.

Adaptive combustion instability control with saturation

The control of a class of combustion systems, suceptible to damage from self-excited combustion oscillations, is considered. An adaptive stable controller, called Self-Tuning Regulator (STR), has recently been developed, which meets the apparently contradictory challenge of relying as little as possible on a particular combustion model while providing some guarantee that the controller will cause no harm. The controller injects some fuel unsteadily into the burning region, thereby altering the heat release, in response to an input signal detecting the oscillation. This paper focuses on an extension of the STR design, when, due to stringent emission requirements and to the danger of flame extension, the amount of fuel used for control is limited in amplitude. A Lyapunov stability analysis is used to prove the stability of the modified STR when the saturation constraint is imposed. The practical implementation of the modified STR remains straightforward, and simulation results, based on the nonlinear premixed flame model developed by Dowling, show that in the presence of a saturation constraint, the self-excited oscillations are damped more rapidly with the modified STR than with the original STR. © 2001 by S. Evesque. Published by the American Institute of Aeronautics and Astronautics, Inc.

Wave propagation in a slowly-varying duct with mean vortical swirling flow

The propagation of unsteady disturbances in a slowlyvarying cylindrical duct carrying mean swirling flow is investigated using a multiple-scales technique. This is applicable to turbomachinery flow behind a rotor stage when the swirl and axial velocities are of the same order. The presence of mean vorticity couples acoustic and vorticity equations which produces an eigenvalue problem that is not self-adjoint unlike that for irrotational mean flow. In order to determine the amplitude variation along the duct, an adjoint solution for the coupled system of equations is derived. The solution breaks down where a mode changes from cut on to cut off. In this region the amplitude is governed by a form of Airy's equation, and the effect of swirl is to introduce a small shift in the origin of the Airy function away from the turning-point location. The variation of axial wavenumber and amplitude along the duct is calculated. In hard-walled ducts mean swirl is shown to produce much larger amplitude variation along the duct compared with a nonswirling flow. Mean swirl also has a large effect in ducts with finite-impedance walls which differs depending on whether modes are co-rotating with the swirl or counter rotating. © 2001 by A.J. Cooper, Published by the American Institute of Aeronautics and Astronautics, Inc.

The absorption of sound by helmholtz resonators with and without flow

Helmholtz resonators are commonly used as absorbers of incident acoustic power. Theoretical and experimental investigations have been performed in the four cases of no mean flow, grazing mean flow, bias mean flow and a combination of grazing and bias mean flows. In the absence of a mean flow, the absorption coefficient (deflned as the proportion of incident energy absorbed) is a non-linear function of the acoustic pressure and high incident acoustic pressures are required before the absorption becomes signiflcant. In contrast, when there is a mean flow present, either grazing or bias, the absorption is linear and thus absorption coefficient is independent of the magnitude of the acoustic pressure, and absorption is obtained over a wider range of frequencies. Non-linear effects are only discernible very close to resonance and at very-high amplitude. With grazing mean flow, there is the undesirable effect that sound can be generated over a range of frequencies due to the interaction between the unsteadily shed vorticity waves and the downstream edge of the aperture. This production is not observed when there is a bias flow because here the vorticity is shed all around the rim of the aperture and swept away by the mean flow. When there is both a grazing mean flow and a mean bias flow, we flnd that only a small amount of bias mean flow, compared with grazing mean flow, is required to destroy the production of acoustic energy. © 2002 by the author(s). Published by the American Institute of Aeronautics and Astronautics, Inc.

The aeroacoustics of transonic helicopter blades

One of the most important issues facing the helicopter industry today is helicopter noise, in particular transonic rotor noise. It is the main factor limiting cruise speeds, and there is real demand for efficient and reliable prediction methods which can be used in the rotor design process. This paper considers the Ffowcs Williams-Hawkings equation applied to a permeable control surface. The surface is chosen to be as small as possible, while enclosing both the blade and any transonic flow regions. This allows the problematic quadrupole term to always be neglected, and requires only near field CFD input data. It is therefore less computationally intensive than existing prediction methods, and moreover retains the physical interpretation of the sources in terms of thickness, loading and shock-associated noise. A computer program has been developed which implements the permeable surface form of retarded time formulation. The program has been validated and subsequently used to validate an acoustic 2-D CFD code. It is fast and reliable for subsonic motion, but it is demonstrated that it cannot be used at high subsonic or supersonic speeds. A second computer program implementing a more general formulation has also been developed and is presently being validated. This general formulation can be applied at high subsonic and supersonic speeds, except under one specific condition. © 2002 by the author(s). Published by the American Institute of Aeronautics and Astronautics, Inc.

A systematic approach to materials information management for aerospace design

The materials information requirements of the aerospace sector are considered, specifically 'consolidation' (management of raw test data), 'analysis' (investigation of material trade-offs) and 'dissemination (secure distribution of data throughout an organization). An information architecture that satisfies the complex requirements of the aerospace materials industry is discussed and a case-study is presented. © 2003 by Granta Design Limited. Published by the American Institute of Aeronautics and Astronautics, Inc.