考虑传质传热效应的翼型积冰计算

本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

翼型表面明冰的数值模拟

提出一种翼型表面明冰的数值模拟方法.基于欧拉两相流理论,建立空气、水滴两相流控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击;采用考虑了粗糙度影响的边界层积分法求解结冰表面的对流换热特性,基于经典Messinger模型提出明冰数值模拟的热力学模型;利用冰层时间推进法模拟积冰过程,采用积冰法向生长假设生成积冰外形.对NACA 0012翼型在不同环境温度下的明冰冰形进行了预测,并与实验结果进行了对比.      

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程；四阶龙格一库塔法求解水滴运动轨迹方程；逐次二分法寻找水滴撞击极限，然后计算局部水滴收集率；分析翼型表面的传质、传热项，建立相应质量、能量方程并迭代求解，得到翼型表面结冰质量；根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好，结冰过程中的特征参数与冰型及相应冰形特征吻合。     

多段翼构型的积冰数值模拟

基于非结构网格提出一种针对多段翼构型的积冰数值模拟方法.通过求解气体欧拉方程得到多段翼的绕流流场,基于空气绕流流场的气流速度分布,采用迭代法对水滴控制方程进行数值求解,得到多段翼各翼面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对某典型三段翼霜冰和明冰的结冰情况进行了预测研究,根据多段翼的流场特点提出一种明冰冻结模型的求解思路,并分析了积冰对多段翼性能的影响.      

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

单边非对称螺杆齿形曲线及铣刀计算公式的统一形式

将文献［１］，［２］中推导出的各段公式归纳为统一的形式。介绍一种求解接触方程的方法，并附有一个计算实例图。     

飞行器再入机动制导的一种新方法

基于飞行力学的原理，建立了飞行器攻击地面目标的运动方程，进而运用计算结构力学与最优控制的相似性原理，将求解Ｒｉｃｃａｔｉ方程的一种新方法引入飞行器再入机动制导问题中。由于采用的是结构静力学中的子结构消元法，从而避免了本征值问题的求解。该方法大幅度地提高了计算速度，同时具有很好的数值稳定性，计算结果表明该方法的有效性。     

连续广义系统的镇定性分析

本文主要研究奇异系统的稳定性和镇定性问题，通过Lyapunov议程2和Ricaati方程，给出了奇异系统的渐近稳定性和镇定性的判据。     

基于数值模拟的GH536高温合金铣削力模型

为了研究GH536高温合金的铣削力来指导生产,利用DEFORM-3D对GH536的铣削过程进行了建模和仿真;设计铣削正交仿真方案,对仿真结果进行极差分析。结果表明:铣削深度、铣削宽度和每齿进给量对F_x、F_y和F_z均成正相关、而铣削速度对F_z正相关,其余负相关,其中铣削深度对铣削力影响程度最大,而铣削速度影响程度最小;利用SPSS对正交仿真数据进行处理,建立GH536铣削力预测模型,并对模型进行回归分析;设计铣削测力试验,将其结果与铣削力预测模型进行对比,并用铣削力理论模型求解结果作为参照对比,结果表明铣削力预测模型能够比理论模型更准确的预测GH536切削过程中的铣削力。     

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

无波动、无自由参数的耗散差分格式

本文通过对“模型问题”的研究发现,在激波附近差分解的虚假波动与差分方程的修正方程式的三阶色散项相关。如果适当地调节色散项的系数使之过激波时该系数按要求变号,则可完全抑制激波附近虚假的波动。 在以上发现的基础上,提出了一个有效的无自由参数、无波动的差分格式。文中证明,这个格式就是“TVD格式”,亦可视为格式的推广。 文中对许多二维、三维流动作了计算,结果表明,本文的格式应用方便且可得到满意的结果,对激波有较高的分辨能力。     

基于影响矩阵的故障诊断方法及其在航空航天领域的应用

提出了系统影响矩阵的概念,在此基础上提出了故障诊断的一种新方法--基于影响矩阵的诊断方法。该方法的核心是系统影响方程和系统测量方程。系统测量方程描述系统的输出（测量）特性;影响方程描述系统的各部件对各输出（测量）的影响特性。该方法能将测量信号的分析处理和故障的定位很好地结合起来。     

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

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

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