基于控制理论的压气机叶型数值优化方法

将基于控制理论的气动优化方法应用于轴流压气机叶型设计.以Euler方程作为流动控制方程,具体推导得出了其相应的伴随方程,分析了边界条件,并给出求解方法.以给定压力分布作为目标函数,将参数化叶型作为设计变量,在求得目标函数对设计变量的梯度信息后,结合BFGS优化算法得到优化方向,更新设计变量完成叶型的优化设计.通过三个算例验证了该叶型优化设计方法的有效性.      

叶栅离散伴随气动优化设计

正叶栅气动优化设计是一个具有复杂约束环境的多变量、非线性、多目标优化问题。基于控制理论的气动优化设计方法,因其在求解目标函数对设计变量的梯度时引入了伴随系统,又被称为伴随方法。与传统的梯度求解方法相比,伴随方法的优势在于其梯度的计算量与设计变量数目无关,计算目标函数对所有设计变量的梯度只需计算一次流场和一次伴随场,且伴随方程是线性偏微分方程,远不及流动方程复杂。其中,离散伴随方法因其伴随方程及边界条件     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式，并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合，发展了一种跨声速机翼气动优化设计方法。为了提高计算效率，将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计，一般通过20~30次迭代就能得到满意的结果。     

基于控制理论的旋翼翼型优化设计

基于粘性Navier-Stokes方程和最优控制理论原理,研究了气动/几何约束条件下多设计变量的旋翼翼型气动优化设计问题。根据给定目标函数的表达形式,在计算坐标系下推导出了相应的伴随方程和边界条件,以及梯度方程的数学表达式,并对流动控制方程和伴随方程进行了有效数值求解。综合流动方程、伴随方程、目标函数敏感性导数和优化算法,发展了一种高效旋翼翼型气动优化方法。通过典型旋翼翼型的算例验证,表明方法的有效性。     

基于共轭方程法的跨音速机翼气动力优化设计

 设计状态的机翼气动力特性是设计人员最为关心的指标, 应用控制理论设计方法进行了有升力约束情形下跨音速机翼阻力优化设计研究, 根据给定的目标函数推导了相应的共轭方程和边界条件, 研究了共轭方程的数值求解方法, 以及计算目标函数对设计变量的敏感性导数时所涉及的度量矩阵变分求解问题, 研究了流场计算、共轭方程数值求解、敏感性导数求解和拟牛顿优化算法这几个主要方面的有效结合问题, 发展出了一种跨音速机翼气动力优化设计方法, 进行了跨音速机翼气动力优化设计研究验证, 优化后机翼气动力特性有一定程度的改善, 阻力系数能减少20%左右, 而升力系数有所增大, 说明所发展的设计方法是成功的, 该设计方法在跨音速及复杂外形气动设计方面比以往设计方法具有更好的适用性和优越性。     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

基于N-S方程的翼型气动优化设计

本文简要回顾了气动优化设计的最新发展，并采用连续伴随方法对粘性条件下的翼型气动外形进行了优化设计。设计过程中，采用B样条曲线的控制顶点作为设计变量，通过求解伴随方程得到目标函数关系设计变量的梯度，最后用最速下降法和BFGS变尺度法对二维翼型进行了外形优化，验证了本方法直接、快速的设计能力。     

考虑飞行器动力系统进排气效应的设计参数灵敏度分析研究

面向飞行器内外流一体化设计,基于自主研发的大规模并行化结构化网格RANS求解器以及离散伴随方程求解器,开展了考虑推进系统动力状态下进排气边界条件的变分研究,通过链式求导法则避免直接对守恒变量变分,进一步引入中间变量大幅度简化了进排气边界条件变分的难度,建立了考虑进排气效应的设计变量灵敏度高效分析方法,并通过TPS标准模型计算验证了进排气数值模拟精度,与有限差分对比验证了灵敏度计算精度,以翼上发动机气动布局进排气影响数值模拟为例,系统分析了低速、高速、定攻角、定升力状态,推进系统有无动力工况灵敏度的变化以及影响机理。     

基于伴随方程和自动微分的雷达散射截面表面灵敏度计算

反隐身技术的发展对军用飞行器的隐身性能提出更高要求，针对经典电磁场伴随方法无法获得表面灵敏度，难以为优化设计提供直观指导，且梯度求解时需反复填充阻抗矩阵，当设计变量、入射角度较多时梯度求解效率下降的问题，提出了基于伴随方程和自动微分的雷达散射截面表面灵敏度计算方法，并结合基函数特点提出了表面灵敏度的稀疏存储和稀疏矩阵-矢量相乘方法，避免了直接计算表面灵敏度时计算量、存储量无法承受的问题。提出的表面灵敏度计算方法可以通过一次阻抗矩阵偏导数求解得到所有网格点的表面灵敏度，避免反复填充阻抗矩阵的问题。当入射角度改变时，求解任意数量设计变量梯度的计算量约为8次矩阵-向量乘，显著提升基于伴随方程的梯度计算效率，为作战飞机气动隐身一体化优化设计提供有力技术支撑。     

基于控制理论和NS方程的气动设计方法研究

研究了基于控制理论和NS方程的气动设计方法,针对给定的目标函数表达形式,应用该设计理论在计算坐标下详细推导了相应的共轭方程及边界条件具体表达形式,以及梯度方程求解表达式,通过合理的数学变换,导出了共轭方程在笛卡尔坐标系下的直观表达形式,发展了有效的共轭方程数值求解方法,通过流动控制方程数值求解、共轭方程数值求解、目标函数对设计变量的梯度求解和优化算法等方面的有效结合,研究与发展了一种新的气动设计方法,以二维机翼气动设计为例,成功进行了亚、跨音速情形下的相关设计算例研究,研究结果表明应用控制理论和NS方程的气动设计方法在设计理论、适用性以及时间花费等方面都有着很好的特色和优点,且设计结果也更为可靠.     

一种基于共轭方程法求解黏性反问题的简化方法

 对于给定压力分布的黏性气动反问题,考虑到壁面微小扰动造成的压力变化主要由势流作用引起,因此可以简化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解,则其相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解,网格数减少,同时收敛速度加快,简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明,对于附着的边界层流动简化方法计算得到的敏感性导数具有较高的精度,能够有效完成反问题设计且减少总的计算耗时。     

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。     

一种三维弹性机翼气动优化新方法

基于控制理论方法和有限元分析软件,通过将网格生成、流场计算、有限元分析、粘性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,研究与发展了计及静气动弹性影响的三维机翼气动优化设计方法,其中,雷诺平均Navier-Stokes方程为主控方程,气动载荷和结构静弹性变形量由气动/结构方程的耦合迭代求解得到,目标函数梯度信息由共轭方程数值求解得到。典型大展弦比机翼气动减阻设计结果表明：研究及发展的计及静气动弹性影响的三维机翼气动优化设计方法是有效的,能够有效考虑静气动弹性的影响。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

Adjoint-based robust optimization design of laminar flow wing under flight condition uncertainties

It is an inherent uncertainty problem that the application of laminar flow technology to the wing of large passenger aircraft is affected by flight conditions. In order to seek a more robust natural laminar flow control effect, it is necessary to develop an effective optimization design method. Meanwhile, attention must be given to the impact of crossflow(CF) instability brought on by the sweep angle. This paper constructs a robust optimization design framework based on discrete adjoint methods ...     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。