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

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

钝体边界层失稳特性分析

推导了适合可压缩二维/轴对称边界层的线性抛物化稳定性方程。在法向采用4阶中心格式,流向采用1阶向后差分,对线性抛物化稳定性方程进行离散。通过来流马赫数为4.5的平板与曲面边界层稳定性问题验证了程序的正确性。讨论了来流马赫数6的球钝锥边界层内,扰动波沿端头向下游演化的失稳特性。发现了一种新的失稳机制。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

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

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

一种新PSO混合算法在直升机配平中的应用

直升机配平计算是动力学分析的基础,其实质是求解高维复杂的非线性方程组。针对经典算法与智能算法的特点与不足,提出了一种求解非线性方程组的新粒子群方法。在粒子群(Particle swarm optimization, PSO)算法的基础上,根据模拟退火(Simulated annealing, SA)思想,引入了嵌入式LM (Levenberg-marquardt)优化 算子。该方法充分发挥了3种算法的优势,克服了LM算法初值敏感性,PSO算法易陷入局部极值等问题。通过UH-60A直升机实例配平计算,验证了本文算例模型的准确性。在此基础 上,针对某一前飞状态下的配平算例,在收敛可靠性和计算效率上通过与其他算法进行对比,表明该算法具有可靠的收敛性和较高的计算效率,进一步验证了该算法在配平问题上的可信度与实用性,为直升机飞行动力学问题的处理提供了一种新的有效方法。     

风力机翼型阵风响应计算与分析

基于网格速度方法,对阵风条件下风力机翼型运动的气动特性进行数值计算。中心格式有限体积法、双时间推进法求解预处理后的非定常N-S方程。计算结果表明:在弦向阵风作用下,翼型的运动能有效降低响应波动峰值,即弦向阵风对运动影响较小;而在法向阵风作用下,气动响应波动很大。计算翼型在不同参数阵风作用下的响应过程,阵风的频率变化对响应影响不大。弦向阵风的幅值增大,叶片运动对阵风响应峰值减缓效果增加;法向阵风的幅值增大1倍,响应峰值也增加1倍。     

自适应翼型的计算和分析

计算并讨论了不同马赫数（Ma）和攻角(α）下简化的自适应翼型的舵面偏转角的规律。计算中采用代数方法生成计算网络，用有限全积法离散二维可压缩流动的Fuler方程，采用了矢通量分裂Van-leer格式离散无粘通量，用隐式的时间积分和多重网格加速收敛。用数值优化方法取得自适应翼型最优的舵面偏角，并计算其气动特性。为验证自适应翼型增益的效果，计算了双目标优化翼型的气动特性，本文计算了对称圆弧翼型和NACA65A006翼型。计算结果表明，自适应翼型比双目标优化的翼型有更好的气动特性。     

二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

数值模拟二维跨音速叶栅可压粘性流动

用有限体积法数值模拟了任意回转流面跨音速压气机叶栅中的粘性流动，四步龙格－库塔法用于N－S方程的时间推进。给出了一种抑制叶栅计算中“数值失速”的方法，并采用隐式残差光顺技术加速收敛。用所发展的方法计算了三种叶栅在宽广来流Ma数条件下的流动，计算结果与实验结果吻合良好，流场激波具有较高的分辨率。