Active Subspaces for Shape Optimization

Design and optimization benefit from understanding the dependence of a quantity of interest (e.g., a design objective or constraint function) on the design variables. A low-dimensional active subspace, when present, identifies important directions in the space of design variables; perturbing a design along the active subspace associated with a particular quantity of interest changes that quantity more, on average, than perturbing the design orthogonally to the active subspace. This low-dimensional structure provides insights that characterize the dependence of quantities of interest on design variables. Airfoil design in a transonic flowfield with a parameterized geometry is a popular test problem for design methodologies. Two particular airfoil shape parameterizations are examined (parametric sections and class-shape transformation), and the active subspaces present in two common design quantities of interest (transonic lift and wave-drag coefficients) are studied under each shape parameterization. The two parameterizations are mathematically related by a simple change of variable and subsequent odd-power polynomial series. The active subspaces enable low-dimensional approximations of lift and drag that relate to physical perturbations of the airfoil. In particular, a two-dimensional approximation of both transonic lift and drag is obtained and interpreted, where 1) each dimension is a particular linear combination of the design variables resulting in intuitive airfoil perturbations, and 2) the linear combination weights are related between lift and drag. It is shown how these approximations along with similarities between the lift and drag active subspaces inform a multi-objective design problem.

Shared low-dimensional subspaces for propagating kinetic uncertainty to multiple outputs

Global sensitivity analysis and uncertainty quantification of soot formation in an n-dodecane spray flame

From regional sensitivity to intra-sensitivity for parametric analysis of free-form shapes: Application to ship design

Active subspaces for sensitivity analysis and dimension reduction of an integrated hydrologic model

Gradient free active subspace construction using Morris screening elementary effects

Groundwater models always involve high‐dimensional parameters, which makes computationally tractable data assimilation using surrogate models very challenging. To address this issue, one common practice is to employ dimension reduction (DR) techniques. Nevertheless, traditional DR methods are usually implemented based on prior parameter statistics, that is, without considering the inherent system dynamics. Here, we show that when significant difference in parameter sensitivity exists, further efficiency can be achieved by adopting a supervised DR method, that is, the active subspace (AS) method. To avoid non‐trivial efforts in calculating the gradient information needed in the standard AS method, a cluster‐based gradient‐free AS (GFAS) method is developed in this study. By combining GFAS with Gaussian process regression, a surrogate model for the CPU‐demanding groundwater model can be adaptively constructed to accelerate data assimilation. Furthermore, a compensation scheme is proposed to cope with uncertainty underestimation caused by DR. The developed approach is tested with numerical experiments and field cases, which illustrated that the new approach is more efficient than the previously developed unsupervised ones by incorporating sensitivity information. Although an iterative ensemble smoother is employed in this study, the proposed method can also be used in other data assimilation approaches, such as Markov chain Monte Carlo and ensemble Kalman filter.

We present a distributed active subspace method for training surrogate models of complex physical processes with high-dimensional inputs and function valued outputs. Specifically, we represent the model output with a truncated Karhunen–Loeve (KL) expansion, screen the structure of the input space with respect to each KL mode via the active subspace method, and finally form an overall surrogate model of the output by combining surrogates of individual output KL modes. To ensure scalable computation of the gradients of the output KL modes, needed in active subspace discovery, we rely on adjoint-based gradient computation. The proposed method combines benefits of active subspace methods for input dimension reduction and KL expansions used for spectral representation of the output field. We provide a mathematical framework for the proposed method and conduct an error analysis of the mixed KL active subspace approach. Specifically, we provide an error estimate that quantifies errors due to active subspace projection and truncated KL expansion of the output. We demonstrate the numerical performance of the surrogate modeling approach with an application example from biotransport.

We present a distributed active subspace method for training surrogate models of complex physical processes with high-dimensional inputs and function valued outputs. Specifically, we represent the model output with a truncated Karhunen-Lo\`eve (KL) expansion, screen the structure of the input space with respect to each KL mode via the active subspace method, and finally form an overall surrogate model of the output by combining surrogates of individual output KL modes. To ensure scalable computation of the gradients of the output KL modes, needed in active subspace discovery, we rely on adjoint-based gradient computation. The proposed method combines benefits of active subspace methods for input dimension reduction and KL expansions used for spectral representation of the output field. We provide a mathematical framework for the proposed method and conduct an error analysis of the mixed KL active subspace approach. Specifically, we provide an error estimate that quantifies errors due to active subspace projection and truncated KL expansion of the output. We demonstrate the numerical performance of the surrogate modeling approach with an application example from biotransport.

Fuel octane number prediction based on topological indices and active subspace method

Reliability analysis and design of a single diode solar cell model using polynomial chaos and active subspace

A physics-guided machine learning for multifunctional wave control in active metabeams

PurposeAerodynamic shape optimisation is a complex problem usually governed by transonic non-linear aerodynamics, a high dimensional design space and high computational cost. Consequently, the use of a numerical simulation approach can become prohibitive for some applications. This paper aims to propose a computationally efficient multi-fidelity method for the optimisation of two-dimensional axisymmetric aero-engine nacelles.Design/methodology/approachThe nacelle optimisation approach combines a gradient-free algorithm with a multi-fidelity surrogate model. Machine learning based on artificial neural networks (ANN) is used as the modelling technique because of its ability to handle non-linear behaviour. The multi-fidelity method combines Reynolds-averaged Navier Stokes and Euler CFD calculations as high- and low-fidelity, respectively.FindingsRatios of low- and high-fidelity training samples to degrees of freedom of nLF/nDOFs = 50 and nHF/nDOFs = 12.5 provided a surrogate model with a root mean squared error less than 5% and a similar convergence to the optimal design space when compared with the equivalent CFD-in-the-loop optimisation. Similar nacelle geometries and aerodynamic flow topologies were obtained for down-selected designs with a reduction of 92% in the computational cost. This highlights the potential benefits of this multi-fidelity approach for aerodynamic optimisation within a preliminary design stage.Originality/valueThe application of a multi-fidelity technique based on ANN to the aerodynamic shape optimisation problem of isolated nacelles is the key novelty of this work. The multi-fidelity aspect of the method advances current practices based on single-fidelity surrogate models and offers further reductions in computational cost to meet industrial design timescales. Additionally, guidelines in terms of low- and high-fidelity sample sizes relative to the number of design variables have been established.

Design space exploration can reveal the underlying structure of design problems of interest. In a set-based approach, for example, exploration can identify sets of designs or regions of the design space that meet specific performance requirements. For some problems, promising designs may cluster in multiple regions of the design space, and the boundaries of those clusters may be irregularly shaped and difficult to predict. Visualizing the promising regions can clarify the design space structure, but design spaces are typically high-dimensional, making it difficult to visualize the space in three dimensions. Techniques have been introduced to map high-dimensional design spaces to low-dimensional, visualizable spaces. Before the promising regions can be visualized, however, the first task is to identify how many clusters of promising designs exist in the high-dimensional design space. Unsupervised machine learning methods, such as spectral clustering, have been utilized for this task. Spectral clustering is generally accurate but becomes computationally intractable with large sets of candidate designs. Therefore, in this paper a technique for accurately identifying clusters of promising designs is introduced that remains viable with large sets of designs. The technique is based on spectral clustering but reduces its computational impact by leveraging the Nyström Method in the formulation of self-tuning spectral clustering. After validating the method on a simplified example, it is applied to identify clusters of high performance designs for a high-dimensional negative stiffness metamaterials design problem.

In this paper, we present an optimal shape design method for a dielectric microlens which is used to focus an incoming infrared plane wave in wideband, by exploiting the finite difference time domain (FDTD) technique and the topology optimization technique. Topology optimization is a scheme to search an optimal shape by adjusting the material properties, which are design variables, within the design space. And by introducing the adjoint variable method, we can effectively calculate a derivative of the objective function with respect to the design variable. To verify the proposed method, a shape design problem of a dielectric microlens is tested when illuminated by a transverse electric (TE)-polarized infrared plane wave. In this problem, the design variable is the dielectric constant within the design space of a dielectric microlens. The design objective is to maximally focus the incoming magnetic field at a specific point in wideband.