Search results

The purpose of this paper is to propose the use of the continuous adjoint method as a tool to identify the appropriate location and “type” (suction or blowing) of steady jets used in active flow control systems.

The method is based on continuous adjoint and covers both internal and external aerodynamics. The adjoint equations, including the adjoint to the SpalartAllmaras turbulence model and their boundary conditions are formulated. At the cost of solving the flow and adjoint equations just once, the sensitivity derivatives of the objective function with respect to hypothetical (normal) jet velocities at all wall nodes are computed. Comparisons of the computed sensitivities with finite differences and parametric studies to assess the present method are included.

Though the sensitivities are computed for zero jet velocities, they adequately support decision making on: the recommended location of jet(s), at boundary nodes with high absolute valued sensitivities; and the selection between suction or blowing jets, based on the sign of the computed sensitivities. Regarding adjoint methods, two important findings of this work are: the role of the adjoint pressure which proves to be an excellent sensor in flow control problems; and the prediction accuracy of the proposed adjoint method compared to the commonly made assumption of “frozen turbulence”.

First use of the continuous adjoint method using full differentiation of the turbulence model, in flow control optimization. A low‐cost design tool for recommending some of the most important jet characteristics.

This paper aims to conduct the optimization of the multi-stage gas turbine with the effect of the cooling air injection based on the adjoint method.

Continuous adjoint method is combined with the S2 surface code.

The optimization of the stagger angles, stacking lines and the passage can improve the attack angles and restrain the development of the boundary, reducing the secondary flow loss caused by the cooling air injection.

The aerodynamic performance of the gas turbine can be improved via the optimization of blade and passage based on the adjoint method.

The results of the first study on the adjoint method applied to the S2 surface through flow calculation including the cooling air effect are presented.

The goal of this paper is to give a comprehensive and short review on how to compute the first- and second-order topological derivatives and potentially higher-order topological derivatives for partial differential equation (PDE) constrained shape functionals.

The authors employ the adjoint and averaged adjoint variable within the Lagrangian framework and compare three different adjoint-based methods to compute higher-order topological derivatives. To illustrate the methodology proposed in this paper, the authors then apply the methods to a linear elasticity model.

The authors compute the first- and second-order topological derivatives of the linear elasticity model for various shape functionals in dimension two and three using Amstutz' method, the averaged adjoint method and Delfour's method.

In contrast to other contributions regarding this subject, the authors not only compute the first- and second-order topological derivatives, but additionally give some insight on various methods and compare their applicability and efficiency with respect to the underlying problem formulation.

In magnetostatics, topology optimization (TO) addresses the problem of finding the distributions of both current densities and ferromagnetic materials to comply with fixed magnetic specifications. The purpose of this paper is to develop TO in order to design Hall-effect Thrusters (HETs).

In fact, TO problems are known to be large-scale optimization problems. The authors therefore adopt the adjoint method to reduce the computation time required to obtain the gradient information. In this paper, they illustrate the continuous variant of the adjoint method in the context of magnetostatics TO. Herein, the authors propose an implementation of the adjoint method then use it within a gradient-based optimization solver fmincon-MATLAB to solve a HET TO design problem.

By comparison with finite difference method, the authors validate the accuracy of the suggested implementation of the adjoint method. Then, they solve a large-scale HET TO design problem. The resultant design of TO is distinctly original and not intuitive.

In this paper, the authors introduce TO as a tool that has allowed them to explore new and innovative design of a HET. However, although the design presented is original, its manufacture is not feasible. Thus, a discussion section has been included at the end of paper to suggest a possible way to concretize topological solutions.

TO helps to explore more original design possibilities. In this paper, the authors present an implementation of the adjoint method that makes it possible to solve efficiently and in less central processing unit time large-scale TO design problem.

An easy implementation of the adjoint method is presented in magnetostatics TO. This implementation was first validated by comparison with the finite difference method and then used to solve a large-scale design problem. The result of the TO design problem is distinctly original and non-intuitive.

An aerodynamic shape optimization algorithm is presented, which includes all aspects of the design process: parameterization, flow computation and optimization. The purpose of this paper is to show that the Class‐Shape‐Refinement‐Transformation method in combination with an Euler/adjoint solver provides an efficient and intuitive way of optimizing aircraft shapes.

The Class‐Shape‐Transformation method was used to parameterize the aircraft shape and the flow was computed using an in‐house Euler code. An adjoint solver implemented into the Euler code was used to compute the required gradients and a trust‐region reflective algorithm was employed to perform the actual optimization.

The results of two aerodynamic shape optimization test cases are presented. Both cases used a blended‐wing‐body reference geometry as their initial input. It was shown that using a two‐step approach, a considerable improvement of the lift‐to‐drag ratio in the order of 20‐30 per cent could be achieved. The work presented in this paper proves that the CSRT method is a very intuitive and effective way of parameterizating aircraft shapes. It was also shown that using an adjoint algorithm provides the computational efficiency necessary to perform true three‐dimensional shape optimization.

The novelty of the algorithm lies in the use of the Class‐Shape‐Refinement‐Transformation method for parameterization and its coupling to the Euler and adjoint codes.

In the shadow of X‐ray tomography and nuclear magnetic resonance tomography a new tomographic technique based on low‐frequency electric currents has been successfully developed in the past decade. Impedance computed tomography (ICT), although it gives poor spatial resolution images, is unmatched in certain cases. Two different approaches to the sensitivity analysis are presented in this paper. The direct differentiation method has been applied, but the adjoint variable method has not been used in ICT reconstruction algorithms until now. Some problems associated with the adjoint variable method of sensitivity analysis are discussed. The new algorithm is compared with the direct differentiation approach.

Topology optimization is a method used for developing optimized geometric designs by distributing material pixels in a given design space that maximizes a chosen quantity of interest (QoI) subject to constraints. The purpose of this study is to develop a problem-agnostic automatic differentiation (AD) framework to compute sensitivities of the QoI required for density distribution-based topology optimization in an unstructured co-located cell-centered finite volume framework. Using this AD framework, the authors develop and demonstrate the topology optimization procedure for multi-dimensional steady-state heat conduction problems.

Topology optimization is performed using the well-established solid isotropic material with penalization approach. The method of moving asymptotes, a gradient-based optimization algorithm, is used to perform the optimization. The sensitivities of the QoI with respect to design variables, required for optimization algorithm, are computed using a discrete adjoint method with a novel AD library named residual automatic partial differentiator (Rapid).

Topologies that maximize or minimize relevant quantities of interest in heat conduction applications are presented. The efficacy of the technique is demonstrated using a variety of realistic heat transfer applications in both two and three dimensions, in conjugate heat transfer problems with finite conductivity ratios and in non-rectangular/non-cuboidal domains.

In contrast to most published work which has either used finite element methods or Cartesian finite volume methods for transport applications, the topology optimization procedure is developed in a general unstructured finite volume framework. This permits topology optimization for flow and heat transfer applications in complex design domains such as those encountered in industry. In addition, the Rapid library is designed to provide a problem-agnostic pathway to automatically compute all required derivatives to machine accuracy. This obviates the necessity to write new code for finding sensitivities when new physics are added or new cost functions are considered and permits general-purpose implementations of topology optimization for complex industrial applications.

The purpose of this study is to develop and compare two methods of determining the total field, including phase information, when only field amplitudes have been measured on a set of planes in the near field of a complex electromagnetic source.

The first method is a gradient‐based optimization algorithm, based on the adjoint fields. The second method employs an optimization algorithm based on the phase angle gradients of a functional.

The first method, the adjoint field method, is functioning well for a 2D test case. The second method, the phase angle gradient method, gives very good results for 3D test cases.

The next step is to test the methods with results from real measurement data.

The developed methods are intended for use in dosimetry studies and other applications, where the field distribution from electromagnetic sources are needed.

The methods extend previously made constant phase approximations. The present methods are useful in situations where the electromagnetic source is hard to model.

The purpose of this paper is to present a method to calculate the derivative of a functional that depends on the shape of an object. This functional depends on the solution of a linear acoustic problem posed in an unbounded domain. We rewrite this problem in terms of another one posed in a bounded domain using the Dirichlet‐to‐Neumann (DtN) map or the modified DtN map. Using a classical method in shape sensitivity analysis, called the adjoint method, we are able to calculate the derivative of the functional using the solution of an auxiliary problem. This method is particularly efficient because the cost of calculating the derivatives is independent of the number of parameters used to approximate the shape of the domain. The resulting variational problems are discretized using the finite‐element method and solved using an efficient Krylov‐subspace iterative scheme. Numerical examples that illustrate the efficacy of our approach are presented.

This paper aims to present a dynamically adjusted deflated restarting procedure for the generalized conjugate residual method with an inner orthogonalization (GCRO) method.

The proposed method uses a GCR solver for the outer iteration and the generalized minimal residual (GMRES) with deflated restarting in the inner iteration. Approximate eigenpairs are evaluated at the end of each inner GMRES restart cycle. The approach determines the number of vectors to be deflated from the spectrum based on the number of negative Ritz values, k∗.

The authors show that the approach restores convergence to cases where GMRES with restart failed and compare the approach against standard GMRES with restarts and deflated restarting. Efficiency is demonstrated for a 2D NACA 0012 airfoil and a 3D common research model wing. In addition, numerical experiments confirm the scalability of the solver.

This paper proposes an extension of dynamic deflated restarting into the traditional GCRO method to improve convergence performance with a significant reduction in the memory usage. The novel deflation strategy involves selecting the number of deflated vectors per restart cycle based on the number of negative harmonic Ritz eigenpairs and defaulting to standard restarted GMRES within the inner loop if none, and restricts the deflated vectors to the smallest eigenvalues present in the modified Hessenberg matrix.