A new grid deformation technology with high quality and robustness based on quaternion

Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the orthogonality of deformed grid, the displacement of grid points is divided into rotational and translational parts in this paper, and inverse distance weighted interpolation is used to transfer the changing location from boundary grid to the spatial grid. Moreover, the deformation of rotational part is implemented in combination with the exponential space mapping that improves the certainty and stability of quaternion interpolation. Furthermore, the new grid deformation technique named ‘‘layering blend deformation' is built based on the basic quaternion technique, which combines the layering arithmetic with transfinite interpolation(TFI) technique. Then the proposed technique is applied in the movement of airfoil, parametric modeling, and the deformation of complex configuration, in which the robustness of grid quality is tested. The results show that the new method has the capacity to deal with the problems with large deformation, and the ‘‘layering blend deformation' improves the efficiency and quality of the basic quaternion deformation method significantly.     

超临界机翼民机全模腹撑干扰数值模拟研究

本文研究了超临界机翼民机全模型的腹撑支臂干扰，通过数值计算模拟了腹撑支臂位于民机模型腹部的两种支撑位置以及后位支撑时腹撑支臂不同支撑角度的支臂干扰，计算结果可为民机模型试验腹支撑方式位置选择及支臂角度选择提供参考依据。     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

扩散作用对y-Re优模型转捩预测的影响

Gamma．theta转捩模型是Menter及其合作者提出的一个基于经验关系的转捩模型,其经验关系式和模型参数是通过来流米雷诺数均在150万以下的仍系列平板实验在CFX程序中标定的。首先通过33A平板修正了转捩经验关系式．然后对S＆K平板实验进行了数值模拟,研究了较高雷诺数下模型控制方程扩散项系数以及与数值扩散相关的流向网格尺度对gamma．theta模型转捩预测结果的影响。模拟的结果表明,这些因素对转捩位置都有明显影响,并且gamma方程的扩散项系数还影响转捩完成后全湍区的摩阻值。三种因素对转捩预测影响的趋势分别为：gamma方程的扩散项系数越大,theta方程的扩散项系数越小,流向网格尺度越大,则转捩位置越靠前。     

累积复合轧制TZM合金板的组织与性能

采用粉末冶金法制备出TZM合金烧结坯,并将其轧制成2mm厚板材,在此基础上采用累积复合轧制技术(ARB),完成TZM板四次轧制复合,并对不同道次TZM复合板进行性能测试和组织分析。结果表明:随累积变形量增加,复合板晶粒显著细化、组织分布更加均匀,经过三次累积复合后,材料的抗拉强度和界面结合状态最佳。试样晶粒在轧制过程中被拉长、展宽,晶粒断面直径200～500nm,抗拉强度提高50%,2mm TZM复合板最高抗拉强度可达968MPa;材料断裂特征仍为脆性断裂,伸长率变化不明显,轧制复合板伸长率仅为2.0%。     

FL-23风洞级间分离与网格测力试验系统

近年来,随着飞行器研制不断高速化发展,一些型号要求在超声速条件下实现级间分离与网格测力试验。为了实现超声速飞行器级间分离与网格测力风洞试验,利用FL-23跨超声速风洞独有的投放机构,通过对上、下支撑及其控制系统进行改造升级,实现了X、Y两个方向的复合运动,建立了马赫数0.3～4.0的级间分离系统。经过多期型号试验验证,该系统对模型定位控制精度达到要求,满足飞行器高马赫数下开展级间分离与网格测力风洞试验需求。     

构建非结构网格高精度有限体积方法的新途径

综述了笔者所在研究团队在发展非结构网格紧致模板高精度有限体积方法方面的研究进展。非结构网格二阶精度有限体积方法在各类商用和自研计算流体力学（CFD）软件中得到了广泛应用。当进一步提高精度时,遇到的主要困难是高阶有限体积方法重构模板过大的问题。这已成为发展非结构网格高精度有限体积方法的主要技术瓶颈之一。近年来,为解决此问题开展了系统研究。基于首先提出的操作紧致性概念,先后提出了3种紧致模板高精度重构方法,包括紧致最小二乘重构、变分重构和多步重构。这些重构方法的共同特点是可在只包含面相邻单元的紧致模板上实施,并达到任意高阶精度。综述了这3种方法,对这些方法的构造思路、实施策略和进一步发展做了概要的阐述。其中变分重构方法将作为非结构网格高精度有限体积方法的方案之一,在国家数值风洞（NNW）工程中得到发展及应用。     

带冲击-气膜双层壁结构的涡轮导叶冷却特性源项法数值研究

为了降低由于大量冷却孔导致的计算资源消耗,采用源项法对带冲击-气膜双层壁冷却结构的Mark II涡轮导叶传热进行了简化数值模拟,并与不采用源项法的结果进行了对比。分析了冷却孔网格密度以及流量系数经验关系式对源项法模拟精度的影响。结果表明,采用源项法能够基本准确地模拟叶片表面温度分布并且缩短75%的计算时间。当冷却孔进出口面网格最小尺度小于孔径的1/10时,数值模拟结果将不随面网格密度改变。影响模拟精度的主要因素是经验关系式能否准确模拟总冷气流量以及横流对流量系数的影响。     

三维椭圆型方程网格生成中的新源项修正法

构造一种新的源项修正方法用于三维椭圆型方程网格生成，可以实现对网格与边界正交及网格与边界间距的双重直接控制．     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算