Shape Optimization for 3D Turbulent Flows Using Automatic Differentiation

Abstract

We present through this paper some new results with our approach for optimal shape design based on a CAD-free framework for shape and unstructured mesh deformations, automatic differentiation by program for the gradient computation and mesh adaption by metric control. Automatic differentiation allows for an easy, reliable and fast discrete adjoint computation. We managed to get the Jacobian of our Navier-Stokcs solver including a k-e turbulence model and wall-laws. The CAD-free framework is shown to be particularly convenient for optimization when the mesh connectivities and the shape discretization arc variable during optimization. Using these ingredients constrained optimization for turbulent transonic flows has been investigated in both 2 and 3D.     

DNS of Turbulent Flows in Ducts with Complex Shape

We carry out Direct Numerical Simulation (DNS) of flows in closed straight ducts with complex peripheral shape. To perform the simulations the Navier-Stokes equations in cylindrical coordinates are discretized by a second-order finite difference scheme, and the immersed-boundary technique is used to resolve the flow close to walls of complex shape. The basic geometry is a circular pipe of radius R, with imposed sinusoidal perturbations of the type \(\eta R \sin (N_{w}\theta )\). Simulations by varying N _w at fixed η were performed to investigate the effect of the perturbation wavenumber. Additional simulations by fixing N _w and varying η also allow to investigate the influence of the amplitude of the wall corrugations. The modifications of the near-wall structures due to change in the shape of the walls are well depicted through contour plots of the radial component of the vorticity. The presence of geometrical disturbances anchors the structures at the locations where curvature changes, and the shape of the structures is strongly linked to the amplitude of the wall corrugation. Our interest is also in understanding the influence of the shape of the surface on wall friction. We were expecting some changes in the profile of the total stress with respect to that of the circular pipe, which instead were not found. This is a first indication that changes in the near-wall region do not affect the outer region, and that Townsend’s similarity hypothesis holds.     

On the Shape and Dimensions of Three-Dimensional Cavities in Supercavitating Flows

The shape of long, trailing cavities behind three-dimensional headforms is discussed. The case of a flat elliptic wing is specifically treated. Three distinct shape regimes are found: quasi-planar, long-flat, spheroidal. These appear in successively higher speed ranges (lower cavitation numbers, ). It is argued that the cavities may be replaced by surrogates in the form of slender ellipsoids. The pressures on these are almost constant and correspond to a cavitation number equal to twice their longitudinal added mass coefficient, k₁. A heuristic theory based on kinetic energy fields is given, relating k₁ to the product of headform drag and cavity length. This theory correlates with an exact theory in the same form given by Garabedian for axi-symmetric cones and also with its extension to planar flows. Results are given here for the shape of the cavity behind an elliptic wing of any aspect ratio, given drag, and cavitation number. Specific formulae are given in the form, = f (C_D/AR), for the transition between the quasi-planar and long-flat regime, and the long-flat and spheroidal regime.     

On Some Background Flows for Tsunami Waves

With the aim to describe the state of the sea in a coastal region prior to the arrival of a tsunami, we show the existence of background flow fields with a flat free surface which model isolated regions of vorticity outside of which the water is at rest.     

Some Pitfalls in the Modelling of Convective Flows in Porous Media

Some modelling deficiencies in various recent papers, on convective flows in porous media, are pointed out and discussed. These deficiencies are related to the Forchheimer coefficient, mixed double-diffusive convection, magnetohydrodynamic mixed convection, convection in a rarefied gas, and geophysical phenomena.     

悬式绝缘子的形状优化

本文将带自适应运动极限的序列线性规划法应用于悬式绝缘子瓷头的形状优化问题中。建立了优化的数学模型;给出了相应的算法;分析了优化计算的结果;为悬式绝缘子的形状优化设计提供了可行的有效手段。     

Analytic Sensitivity Analysis for Shape Optimization

ＩｎｔｒｏｄｕｃｔｉｏｎＳｈａｐｅｏｐｔｉｍｉｚａｔｉｏｎｆｏｒｅｌａｓｔｏｍｅｒｉｓａｎｉｍｐｏｒｔａｎｔｍｅａｓｕｒｅｔｏｒｅｌｉｅｖｅｓｔｒｅｓｓｃｏｎｔｒｉｔｉｏｎ ,ｐｒｅｖｅｎｔｂｒｅａｋｉｎｇ ,ｉｍｐｒｏｖｅｌｏａｄ＿ｂｅａｒｉｎｇｃａｐａｃｉｔｙ .Ｓｔｒｕｃｔｕｒａｌｒｅｐｅｔｉｔｉｏｎａｎａｌｙｓｉｓ,ｓｅｎｓｉｔｉｖｉｔｙａｎａｌｙｓｉｓａｎｄｏｐｔｉｍａｌａｌｇｏｒｉｔｈｍａｒｅｔｈｒｅｅｉｍｐｏｒｔａｎｔｌｉｎｋｓｏｆｓｈａｐｅｏｐｔｉｍｉｚａｔｉｏｎ .Ｉｎｔｈｅｃｏｕｒｓ…     

弹性接触问题的形状优化设计方法

设计模型复杂和结构的高度非线性等问题使得目前结构优化技术中很少涉及接触问题。面向对象的有限元法利用了面向对象编程的优点，有效地提高了编程的效率，并便于程序的扩充。在Visual C 软件平台上，开发编制了面向对象的有限元法程序，并对含接触问题的结构进行形状优化设计，重点推导并给出了含接触问题有限元法的灵敏度全解析分析的表达式，包括位移和应力的灵敏度分析、接触条件的灵敏度分析及相应的边界支配方程等。最后通过算例，验证了所提出的优化方法的正确性。     

结构形状优化与智能模型化

阻碍形状优化软件广泛应用的障碍之一是依据自然设计度量描述和建立设计模型、分析模型、优化模型以及实现三个模型之间的转换.本文称这一困难为结构形状优优设计软件的适用性.本文提出了一种方法用来动态地确定平面连续体结构形状优化过程中的边界,应用基于人工智能方法的启发式规则与技术,自动生成由设计单元法表示的几何模型,也就是将一个结构自动剖分成若干个大的四边形映射单元.这些大单元对于进一步的网络生成是必要的,同时也是向全自动计算机辅助形状优化系统前进的重要一步.     

桁架结构形状与尺寸组合优化

提出了一种渐进优化方法 ,通过分离设计变量 ,分别采用渐进节点移动法和满应力法优化桁架结构的形状和尺寸。通过循环迭代 ,耦合形状和尺寸变量间的相互作用关系 ,以期得到结构的最小重量。结构受到应力、局部稳定和节点位移等多种约束。最优解由初始结构设计渐进得到。该方法的有效性通过二个算例得到证实     

Some Steady MHD Flows of the Second Order Fluid

The exact analytic solutions of two problems of a second order fluid in presence of a uniform transverse magnetic field are investigated. The governing equation is of fourth order ordinary differential equation and is solved using perturbation method. In the first problem we discuss the flow of a second order fluid due to non-coaxial rotations of a porous disk and a fluid at infinity. In second problem the flow of a second order conducting fluid between two infinite plates rotating about the same axis is investigated, with suction or blowing along the axial direction. For second order conducting fluid it is observed that asymptotic solution exists for the velocity both in the case of suction and blowing.     

Some Thermodynamical Aspects of the Solidification of Two-phase Flows

The solidification of a two-phase flow is analyzed. It is shown how fluid-to-fluid phase transformations and kinetic stresses due to relative diffusion between fluid components influence dissipation of energy, friction and heat flux at the interface solid–fluid.     

变截面悬臂梁的形状优化

讨论变截面悬臂梁的形状优化.围绕这一问题,对设计变量的选择做了有益的探讨,以避免无效的设计变量选取导致不真实或不可行的设计.     

Some Global Stability Results for Shear Flows of Viscoelastic Fluids

We consider plane shear flows of viscoelastic fluids. For a number of constitutive models, we prove stability of the rest state for perturbations of arbitrary size. We also consider stability of plane Poiseuille flow in a few special cases. This research was supported by the National Science Foundation under Grant DMS-0405810.     

多边形有限单元形函数的比较研究

多边形有限单元形函数有wachspress插值、Laplace插值和平均值插值三种类型.本文对三种多边形有限单元形函数的性质作了比较研究,给出了三种形函数各自的优点和局限性.Waclaspress和Laplace形函数是有理函数形式,而平均值形函数是无理函数形式.三种形函数均满足单位分解性、线性完备性,且在单元边界上呈线性.在三角形单元上,它们都等价于三角形面积坐标插值.在矩形单元上,Wachspress和Laplace形函数等价于双线性多项式插值形函数.Wachspress和平均值形函数适用于任意凸多边形单元,Laplace形函数更适用于圆内接多边形单元.Wachspress形函数不能推广到含有边节点的单元,平均值形函数可以直接推广到含有边节点的单元.数值试验,验证了本文理论分析的结论.     

A Discrete Adjoint Approach for the Optimization of Unsteady Turbulent Flows

In this paper we present a discrete adjoint approach for the optimization of unsteady, turbulent flows. While discrete adjoint methods usually rely on the use of the reverse mode of Automatic Differentiation (AD), which is difficult to apply to complex unsteady problems, our approach is based on the discrete adjoint equation directly and can be implemented efficiently with the use of a sparse forward mode of AD. We demonstrate the approach on the basis of a parallel, multigrid flow solver that incorporates various turbulence models. Due to grid deformation routines also shape optimization problems can be handled. We consider the relevant aspects, in particular the efficient generation of the discrete adjoint equation and the parallel implementation of a multigrid method for the adjoint, which is derived from the multigrid scheme of the flow solver. Numerical results show the efficiency of the approach for a shape optimization problem involving a three dimensional Large Eddy Simulation (LES).     

在静力与动力荷载作用下拱坝的体形优化

建立了拱坝体形优化的合理而实用的数学模型,提出了有效的求解方法,特别重视工程应用,目前已应用于近三十个实际工程,取得了巨大经济效益,并已应用于250米特高拱坝的设计。作者提出了一系列有效的求解方法,如内力展开法,设计变量指数变换法等等,并首次实现了地震荷载作用下拱坝的体形优化。     

预锻模具形状优化设计与有限元灵敏度分析

采用刚－粘塑性有限元灵敏度分析方法研究锻造过程的预锻模具形状优化设计问题。表示预锻模具形状的三次Ｂ样条曲线的控制系数用作设计变量。以实际终锻件与理想终锻件的形状差异为目标函数,给出了与设计变量有关的目标函数,节点坐标和节点速度等方面的灵敏度及其数学关系,对于典型的轴对称锻造过程,优化设计的预锻模具形状可获得理想的的终锻件形状,为实现净成形锻造提供了一种有效方法。     

结构形状优化设计数值方法的研究和应用

本文论述了连续体结构形状优化设计数值方法的研究和应用进展,讨论了结构模型化、灵敏度分析、优化方法改进、优化实用软件开发以及同CAD技术相结合等问题,介绍了在结构优化设计软件MCADS中采用的方法,并通过工程实例说明了结构形状优化设计的应用及其价值。     

高地应力下大型地下洞室拱形优化研究

高地应力区的大型洞室开挖会引起洞室拱周附近围岩应力集中和塑性区加大以及塑性应变的增加，但过高的应力集中会增加岩爆发生的几率，塑性区的加深会促使围岩失稳；洞室拱形形状会影响围岩应力集中和塑性区大小的分布，因此，以某高地应力下的大型发电洞岩体地质资料为依托，用弹塑性有限元分析高地应力地下大型发电洞室常采用的3种洞室拱形（即单心圆拱、三心圆拱和椭圆拱）对拱部围岩应力的影响，计算结果表明这种影响是非常明显的，有的应力集中系数达到3．33，无论采用哪种拱形，均避不开拱座附近应力集中现象。根据开挖后拱部围岩的塑性耗散能可以判断，无论初始地应力侧压系数是多少，椭圆拱是最优的，其次是三心圆拱。