基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

前缘吹气控制舵面流动分离

针对各种飞行器大舵偏下出现的流动分离问题,在北航D4风洞对旨在消除舵面流动分离的舵面前缘吹气技术进行了研究,为了降低控制分离所用的吹气量,吹气点设置在舵面前缘气流分离点处.应用粒子图像测速(PIV,Particle Image Velocimetry)技术,分析了舵面绕流在吹气量由小变大过程中所经历的3个不同演化阶段;由测压得到的舵面压力分布则显示,前缘吹气造成的引射作用使前缘吸力峰随吹气量增大而增大,这是前缘吹气能够使舵面升力增大的主要机理.实验结果还表明,前缘吹气可明显提高舵面升力,同时也可以显著降低舵面阻力.     

新型二维推力矢量喷管数值模拟

针对现有流体推力矢量控制方案的不足,提出利用喷流附壁效应的新型矢量喷管,借助于尾喷管射流对固壁延伸面的跟随作用控制尾喷流方向,实现推力转向.在此基础上采用限制流量的方法调节喷流的抽吸程度,产生不同的横向压力梯度,达到了矢量化控制推力转向的目的.运用这一概念设计了二维矢量喷管,用数值实验方式验证了喷管的推力转向效果,采用限制流量方法得到的最大矢量角度约13.3°,进一步结合射流控制可以使矢量偏角达到20°以上.通过对该喷管流场的数值计算研究,探讨了该矢量喷管内喷流转向形成的流动机理,从推力损失、转向效率上对喷管的性能特点进行了分析,为下一步开展实验研究奠定了基础.     

数字化组合测量辅助飞机装配质量检测技术

针对单独使用激光跟踪仪测量大型飞机壁板类组件装配质量时存在曲面数据细节信息不易采集、测量效率低等问题,提出基于激光跟踪仪和关节臂测量仪的大型飞机壁板类组件数字化复合测量方法.在保证测量精度的前提下,将关节臂测量仪(AACMM)融入到激光跟踪仪测量场中,使用两种设备进行集成检测,对大型飞机部件成形细节信息进行全面采集.同时根据飞机部件装配过程测量环节的测量特征,在Spatial Analyzer软件基础上开发了飞机壁板类组件数字化复合测量工具集,实现了对测量过程中所涉及到的设备、数据的集中管控.通过实例,验证了测量方法的有效性和效率.      

基于MapReduce的CME参数识别模型并行计算技术

日冕物质抛射（Coronal Mass Ejection,CME）参数识别模型是太阳风预报过程的重要组成部分.在空间环境预报业务中,为提高太阳风预报的准确率,需要提高CME参数识别的精度.模型以计算任务串行的方式运行,运算效率低导致模型运算时间长,不能满足这种需求.CME参数识别模型的物理运算过程相互不独立,其在单节点上的运行方式不能满足并行化要求.基于MapReduce的并行计算框架,改进了CME参数识别模型的计算流程,提出CDMR（CME detection under MapReduce）方法,实现了CME参数识别模型的并行计算,并对比分析CME参数识别模型在串行计算和MapReduce并行计算下的运行时间,提高了模型的识别精度和计算效率.      

平方阻尼在非线性缓冲系统中的作用

本文讨论平方阻尼(以参数a表征)在图(1)所示的缓冲系统中的作用,结果表明:平方阻尼对于降低相对位移峰值d_m总有良好的效果,但不一定降低加速度峰值a_m。在一定速度冲击V下,系统的平方阻尼有两种最佳值,使a_m最小的a_(am)和使缓冲效率γ最高的a_(cm),而平方阻尼最佳条件依赖于系统的缓冲弹簧的刚度变化类型: (1)线性弹簧情况:a_(am)=a_(cm),最佳条件为系统在首程中均匀减速; (2)渐硬弹簧情况:a_(am)=a_(cm),最佳条件为系统在首程中初始和终末加速度相等即; aV~2=ω_0~2q(dm) (3)渐软弹簧情况:a_(cm)>a_(am)。 本文的基础是所谓“首程显峰”,假定即系统在速度冲击V下的响应峰值d_m和a_m出现在第一行程中。当平方阻尼值不太小,缓冲弹簧刚度特性接近反对称时,该假定成立。     

一种基于承载效率准则的复合材料机翼结构布局优化方法(英文)

提出一种针对大型复合材料机翼结构的二级布局优化策略。系统级根据结构形变调整设计指标,子系统级优化结构布局以满足系统级约束。将机翼壁板中各种临界失效载荷的接近程度作为评判结构效率的标准,通过将结构效率作为优化目标,保证了所求解问题的连续性。加筋壁板通过等效的正交异性板模拟,并根据能量原理预测总体失稳载荷。讨论了加强筋支持弹性对于蒙皮抵抗局部失稳性能的影响,并通过神经网络代理模型对其进行逼近。采用解析模型对失效特性做出评估,减少了计算资源的占用。最后以一种大型机翼结构的综合优化作为算例,计算结果满足设计要求,验证了该方法的可行性。     

风力机组叶片的先进翼型族设计

大型风力机组叶片设计需要专用的翼型族,以满足风力机组运行的环境要求。高升阻比、低粗糙度影响是该翼型族的气动性能特点。为了发展我国具有自主知识产权的风力机专用翼型族,在科技部863项目的支持下,北京航空航天大学联合国内气动研究单位,通过理论分析、数值模拟和风洞试验等技术途径,开发了两个适用于大型风力机组叶片的先进翼型族。理论分析和风洞试验表明,设计结果基本达到了预先设定的气动性能指标。