Optimal shape design for the time‐dependent Navier–Stokes flow

This paper is concerned with the problem of shape optimization of two‐dimensional flows governed by the time‐dependent Navier–Stokes equations. We derive the structures of shape gradients for time‐dependent cost functionals by using the state derivative and its associated adjoint state. Finally, we apply a gradient‐type algorithm to our problem, and numerical examples show that our theory is useful for practical purposes and the proposed algorithm is feasible in low Reynolds number flows. Copyright © 2007 John Wiley & Sons, Ltd.     

Sensitivity of a general class of shape functionals to topological changes

The topological derivative represents the first term of the asymptotic expansion of a given shape functional with respect to the small parameter which measures the size of singular domain perturbations. The topological derivative has been successfully applied in the treatment of problems such as topology optimization, inverse analysis and image processing. In this paper, the calculation of the topological derivative for a general class of shape functionals is presented. In particular, we evaluate the topological derivative of a modified energy shape functional associated to the steady-state heat conduction problem, considering the nucleation of a small circular inclusion as the topological perturbation. Several methods were proposed to calculate the topological derivative. In this paper, the so-called topological-shape sensitivity method is extended to deal with a modified adjoint method, leading to an alternative approach to calculate the topological derivative based on shape sensitivity analysis together with a modified Lagrangian method. Since we are dealing with a general class of shape functionals, which are not necessarily associated to the energy, we will show that this new approach simplifies the most delicate step of the topological derivative calculation, namely, the asymptotic analysis of the adjoint state.     

Optimal design of paper machine headboxes

A shape optimization problem for incompressible flows within a stabilized finite element framework is studied. The goal is to develop and test numerical realizations of optimal shape design problems that could be applied to non‐trivial industrial problems. The resulting algorithm is applied to the optimization of the geometry of a tapered header in a paper machine headbox. Copyright © 2000 John Wiley & Sons, Ltd.     

Shape Design Sensitivity Analysis for Optimization Problems with Local and Global Functionals

ABSTRACT

In this paper, problems of sensitivity analysis and shape design for elastic solids are investigated. Optimization criteria are provided by integral and local functionals defined on the internal region and the boundary of the body. New sensitivity analysis relations are derived and implemented in successive optimization algorithms. Important aspects of effective shape optimization algorithms are discussed.     

Continuous and discrete adjoint approaches for aerodynamic shape optimization with low Mach number preconditioning

Discrete and continuous adjoint approaches for use in aerodynamic shape optimization problems at all flow speeds are developed and assessed. They are based on the Navier–Stokes equations with low Mach number preconditioning. By alleviating the large disparity between acoustic waves and fluid speeds, the preconditioned flow and adjoint equations are numerically solved with affordable CPU cost, even at the so‐called incompressible flow conditions. Either by employing the adjoint to the preconditioned flow equations or by preconditioning the adjoint to the ‘standard’ flow equations (under certain conditions the two formulations become equivalent, as proved in this paper), efficient optimization methods with reasonable cost per optimization cycle, even at very low Mach numbers, are derived. During the mathematical development, a couple of assumptions are made which are proved to be harmless to the accuracy in the computed gradients and the effectiveness of the optimization method. The proposed approaches are validated in inviscid and viscous flows in external aerodynamics and turbomachinery flows at various Mach numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of dual kriging to structural shape optimization based on the boundary contour method

Summary The paper presents an approach in which the coupling of dual kriging and the boundary contour method (BCM) is applied to structural shape optimization problems in mechanical engineering design. The problems consist of optimizing the shape of an elastic body, which requires minimizing an objective function subject to some given constraints, such as those of displacement, stress or manufacturing. The originality of the present work is involved with the use of two novel methods that are combined here to solve structural shape optimization problems. The first one, called dual kriging, is a general, versatile interpolation and geometric modeling tool. The second one is a new variant of the boundary element method (BEM), called the BCM, which achieves a further reduction in dimensionality of analysis problems. Based on the advantages of these two methods, the coupling approach presented here is expected to offer an effective as well as a straightforward manner for solving shape optimal design problems. Received 18 December 1997; accepted for publication 21 April 1998     

3-D Aerodynamic Shape Design Using Hessians Based on Incomplete Sensitivities

We investigate the practicability of an optimization algorithm based on the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method for 3D shape design problems, using approximate sensitivities of objective functions, where the contribution of the partial derivatives of the flow state with respect to the control variables is neglected. Therefore, it is worthwhile to investigate how an optimization method based on the Hessian behaves in this context. Indeed, the Hessian should be far from its real value if the gradient approximation is wrong. The optimization methodology is characterized by an unstructured CAD-free framework for shape and mesh deformations, an automatic differentiation of programs for the computation of the gradient of the cost function, and an unstructured flow solver. The redesign of transonic and supersonic wings has been considered and the performance of the BFGS method has been analyzed in comparison with a steepest descent method. Taking into account that a line search is too expensive to be carried out in such problems, a step size proportional to the gradient modulus has been employed for updating the control variables. Numerical results show that the BFGS method does not suffer from the approximation used in the evaluation of sensitivities, and leads to an effective improvement of the efficiency of the optimization methodology. These results can be then considered an a posteriori justification for incomplete sensitivities.     

Design Sensitivity Analysis in Structural Mechanics. III. Effects of Shape Variation

The dependence of static response and eigenvalues on the shape of plates and plane elastic solids is characterized. Shape of elastic bodies is taken as the design variable. The material derivative idea of continuum mechanics is used to obtain expressions for directional derivatives of displacement fields and eigenvalues with respect to a transformation function that defines a shape variation. The result is used to obtain explicit and computable expressions for variations of integral functionals that arise in structural optimization problems.     

One‐shot pseudo‐time method for aerodynamic shape optimization using the Navier–Stokes equations

This paper presents a numerical method for aerodynamic shape optimization problems in compressible viscous flow. It is based on simultaneous pseudo‐time stepping in which stationary states are obtained by solving the pseudo‐stationary system of equations representing the state, costate and design equations. The main advantages of this method are that it blends in nicely with previously existing pseudo‐time‐stepping methods for the state and the costate equations, that it requires no additional globalization in the design space, and that a preconditioner can be used for convergence acceleration which stems from the reduced SQP methods. For design examples of 2D problems, the overall cost of computation can be reduced to less than 2 times the forward simulation runs. Copyright © 2011 John Wiley & Sons, Ltd.     

TOPOLOGY OPTIMIZATION FOR GEOMETRIC STABILITY OF STRUCTURES WITH COMPENSATION DISPLACEMENTS

大量工程问题要求结构的局部区域在不同承载工况下保持位移响应的几何稳定性. 在结构的特定区域引入可以随承载工况调节的补偿位移是实现这一目标的有效手段. 在线弹性小变形范围内,通过最小的变形控制成本,研究了多工况下使结构特定位置的位移响应保持几何稳定性的拓扑优化问题. 设计目标为在保持结构位移响应几何稳定性的同时实现结构的最大刚度;优化模型包含两类设计变量:结构拓扑变量及补偿位移变量,两类变量采用分层寻优技术进行耦合. 采用伴随法分别推导了目标函数对两类设计变量的敏度求解格式. 结果表明,该优化模型能够在兼顾成本的同时较好地实现结构的变形控制目标.     

A biarc-based shape optimization approach to reduce stress concentration effects

In order to avoid stress concentration, the shape boundary must be properly designed via shape optimiza- tion. Traditional shape optimization approach eliminates the stress concentration effect by using free-form curve to present the design boundaries without taking the machin- ability into consideration. In most numerical control （NC） machines, linear as well as circular interpolations are used to generate the tool path. Non-circular curves, such as non- uniform rotational B-spline （NURBS）, need other more ad- vanced interpolation functions to formulate the tool path. Forming the circular tool path by approximating the opti- mal free curve boundary with arcs or biarcs is another op- tion. However, these two approaches are both at a cost of sharp expansion of program code and long machining time consequently. Motivated by the success of recent researches on biarcs, a reliable shape optimization approach is pro- posed in this work to directly optimize the shape boundaries with biarcs while the efficiency and precision of traditional method are preserved. Finally, the approach is validated by several illustrative examples.     

Optimization of the shapes of obstacles in jet-separation flow

The model of an ideal incompressible fluid is used to study the solvability of optimal control problems for the shape of a nozzle which discharges free-boundary fluid flow with and without accounting for gravity (internal aerodynamics) and shape optimization problems for an obstacle with jet separation (external aerodynamics). The qualitative properties of such flows are studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 30–39, May–June, 2007.     

Reducing memory requirements of unsteady adjoint by synergistically using check-pointing and compression

The unsteady adjoint method used in gradient-based optimization in 2D and, particularly, 3D industrial problems modeled by unsteady PDEs may have significant storage requirements and/or computational cost. The reason for this is that the backward in time integration of the adjoint equations requires the previously computed instantaneous flow fields to be available at each time-step. This article proposes remedies to this problem, by extending/upgrading relevant techniques proposed by the group of authors as well as other researchers. Their applicability is wide, even if these remedies are herein demonstrated in shape optimization problems in unsteady fluid mechanics. Check-pointing is in widespread use as it reduces the memory footprint and CPU cost of the optimization with a controllable computational overhead. Alternatively, flow field time-series can be stored in a lossless or lossly compressed form. The novelty of this article is the development of a Compressed Coarse-grained Check-Pointing strategy for second-order accurate schemes in time, by optimally combining check-pointing and lossy compression. The latter includes (a) the incremental Proper Generalized Decomposition (iPGD) algorithm and (b) a hybridization of the iPGD with the ZFP and Zlib algorithms. This is implemented within OpenFOAM, which is used to solve the flow and adjoint equations and conduct the optimization, and assessed in 2D/3D aerodynamic shape optimization problems on unstructured grids. Effectiveness in data reduction, computational cost, and reconstruction accuracy are compared, vis-à-vis also to the “standard” binomial check-pointing technique after adjusting it to second-order accurate schemes in time.     

On the incremental singular value decomposition method to support unsteady adjoint-based optimization

This paper proposes and evaluates an approximation model based on an incremental Singular Value Decomposition (iSVD) algorithm, for unsteady flow field reconstructions, needed for integrating the unsteady adjoint equations backward in time, within a gradient-based optimization loop. Due to the iSVD algorithm, the computational cost of solving the unsteady adjoint equations is reduced considerably, without practically affecting the accuracy of the computed gradient. Approximations to the unsteady flow fields are constructed while solving the time-varying flow equations (moving forward in time) and used to reconstruct these fields during the backward-in-time integration of the continuous adjoint equations. Optimization results obtained using the proposed method are compared to those computed using the binomial checkpointing technique, which acts as the reference method. Test cases for both flow control and shape optimization problems are presented.     

Insensitive Functionals,Inconsistent Gradients,Spurious Minima,and Regularized Functionals in Flow Optimization Problems

We use the simple context of Navier-Stokes flow in a channel with a bump to examine problems caused by the insensitivity of functionals with respect to design parameters, the inconsistency of functional gradient approximations, and the appearance of spurious minima in discretized functionals. We discuss how regularization can help overcome these problems. Along the way, we compare the discretize-then-differentiate and differentiate-then-discretize approaches to optimization, especially as they relate to the issue of inconsistent functional gradients. We close with a discussion of the implications that our observations have on more practical flow control and optimization problems.     

考虑动荷载响应的结构形状优化设计研究

以产生相同位移场为基础,将动荷载转化为一系列的等效静荷载,然后将这些等效静荷载作为多个载 荷工况进行有限元优化设计分析. 通过几个算例的验证,表明这种基于动荷载等效转换的算 法能够反映动荷载对结构的动力响应,对具有大规模自由度的结构形状优化问题是有效的.     

基于显式几何更新算法的地下管道自动形状优化

基于计算力学中的结构优化思想,应用一种新型的显式几何更新算法,自行编制C++程序,实现地下管道形状设计的自动优化。管道内的流体假设为牛顿不可压缩流,并考虑惯性项。优化区域主要为管道竖直方向和水平方向的过渡段。形状优化的设计变量是几何边界的有限元节点坐标,优化目标是实现流体黏性能耗散的最小化。优化过程基于形状梯度,即通过形状敏感度分析来求解目标函数相对于设计变量的偏导数。所使用的显式几何更新算法既可以通过网格清晰描述形状,也可以大范围地自动更新网格。详细介绍了地下管道自动形状优化过程的关键步骤。通过数值算例探讨了不同注入速度、密度和黏度对其最优形状的影响。     

基于虚荷载变量的形状优化和灵敏度分析

基于选择施加在结构“控制点”上的虚荷载作为优化设计变量,针对一种新的承受约束的形状优化数值方法进行了研究。借助于节点位移与虚荷载之间的线性关系,提出了一种新的计算灵敏度系数的解析方法。利用节点移动速度域概念构造了优化新形状产生的计算公式,以结构中节点的最大应力最小化作为优化目标,通过控制网格结点的最大位移量,较好地解决了单元网格在形状优化中的扭曲问题。对三个不同的实例成功地完成了形状优化。