三元二次样条函数及其计算（Ⅰ）——高维数据拟合

半节点二次样条是Ｃ１类连续的，其二阶导函数是阶梯函数，在半节点处产生跳跃。鉴于此，本文利用最小二乘法，获得了一种半节点二次插值样条边界条件确定方法，该方法可以保证二次插值样条在半节点处二阶导数的变化最小，这相当于保证了曲率的变化最小。为了适应四维数据插值的需要，给出了三元二次插值样条的定义及基样条表示，提出了一种边界条件，证明了其存在唯一性。最后把二次插值样条边界条件确定方法推广到三元二次插值样条上去，获得了三元二次插值样条边界条件的确定方法。该方法的意义是，可以直接从插值条件获得边界条件，从而克服了多元样条边界条件难以确定的困难     

孔式抽吸对带间隙高负荷压气机性能的影响

实验研究了端壁孔式附面层抽吸对带间隙的高负荷直列叶栅流动特性和气动性能的影响,通过三种附面层抽吸方案与原型方案的对比分析,探讨了附面层抽吸抑制间隙流动、减小损失的机理。对叶片表面和下端壁进行了墨迹流动显示,并利用五孔探针对叶栅出口气动参数进行了测量,对比分析了叶栅流道内的流场特征和损失分布。结果表明,在间隙内布置抽吸孔能有效降低间隙流动动能,削弱间隙流与主流之间的相互掺混作用,减小因间隙流动引起的端壁区域流动分离,从而达到对三维流动分离的抑制,有效降低损失,且最大降幅达16.7%;间隙流动引起的流动分离和损失在流道中占主导地位,尽管在端壁沿周向合理布置抽吸孔能在一定程度上抑制近端壁的附面层或二次流发展,但这种局部效应仍不能显著改善叶栅的整体性能;端壁上不合理的抽吸控制策略反而可能影响主流的正常流动,导致损失增加,其中方案4损失增加了约3%。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

大展弦比飞机模型边界层转捩模拟技术

针对大展弦比飞机模型高速风洞试验需求,进行了边界层转捩模拟技术研究.研究表明:与传统的金刚砂粗糙带相比,柱状转捩带具有粘贴牢靠、可操作性强等优势,已经推广应用于大展弦比飞机模型高速风洞试验;边界层转捩对大展弦比飞机模型气动特性有重要影响,不仅影响阻力,而且对升力与俯仰力矩均有明显影响;在目前国内低雷诺数风洞条件下建议在边界层固定转捩条件下进行风洞试验,以提高风洞试验精度.     

Scaling of interaction lengths for hypersonic shock wave/turbulent boundary layer interactions

The interaction length induced by Shock Wave/Turbulent Boundary-Layer Interactions (SWTBLIs) in the hypersonic flow was investigated using a scaling analysis, in which the interaction length normalized by the displacement thickness of boundary layer was correlated with a corrected non-dimensional separation criterion across the interaction after accounting for the wall temperature effects. A large number of hypersonic SWTBLIs were compiled to examine the scaling analysis over a wide range of Mach numbers, Reynolds numbers, and wall temperatures. The results indicate that the hypersonic SWTBLIs with low Reynolds numbers collapse on the supersonic SWTBLIs, while the hypersonic cases with high Reynolds numbers show a more rapid growth of the interaction length than that with low Reynolds numbers. Thus, two scaling relationships are identified according to different Reynolds numbers for the hypersonic SWTBLIs. The scaling analysis provides valuable guidelines for engineering prediction of the interaction length, and thus, enriches the knowledge of hypersonic SWTBLIs.     

基于时域有限差分法平面光波的一维吸收边界条件

吸收边界条件的选取在使用时域有限差分法分析设计衍射光学元件的过程中有着举足轻重的作用。文章从麦克斯韦方程组出发,推导出了一维吸收边界条件的解析解及该吸收边界对TEM波的反射系数,用数值模拟验证了在激励源分别处于讨论区域右边、中间,且在TEM波传播的同时一直存在的情况下,一维Mur一阶吸收边界条件和波延迟法的优质吸收效果。对两种边界条件的计算误差进行了对比和分析,结果表明Neum ann边界条件波形稳定早,数值精度高,Mur一阶边界条件波形稳定稍迟,有一定误差,但吸收效果很好。     

斜出口合成射流激励器横流输运特性与边界层控制

研制了一种斜出口合成射流(BSJ)激励器,应用粒子图像测速(PIV)技术获得了合成射流的瞬态和时均流动结构,分析比较了斜出口和常规平直出口激励器流场的异同。该斜出口激励器能将周围气体进行有方向的能量和质量输送控制,时均流场呈现沿壁面的横向流动输运特性。解释了横向流动输运的形成机理,并应用该激励器进行了边界层控制实验。研究结果表明:通过改变激励器的工作电压和频率,可方便有效地实现对平板边界层速度型的调控,使受控边界层速度型更"饱满",这对控制推迟边界层分离非常有利。     

基于原子内聚力与表面能等效的内聚裂纹模型

 从非局部连续介质理论出发,采用一种新的理性力学方法对裂纹前缘内聚应力分布规律进行研究。首先,在内聚裂纹表面引入非局部应力边界条件,从而将内聚区内表面诱发张力(非局部表面残余)与内聚应力等价联系起来;然后,利用能量平衡关系,得到仅与表面能密度相关的I型裂纹内聚力新的本构方程。最后,在推导结果的基础上,计算一个具体的脆性断裂算例研究内聚区内表面能与内聚应力随裂纹张开位移(COD)变化的分布规律。由计算结果发现,裂纹尖端应力奇异性消除,且应力最大值不一定出现在裂纹尖端,而是发生在裂纹尖端周围的内聚区内。     

介质阻挡放电等离子体气动激励的动量特性

从诱导气流速度和体积力两个方面开展了介质阻挡放电等离子体气动激励的动量特性研究,对不同激励参数下的诱导气流速度进行了测量与结果分析;通过对体积力的实验测量与理论计算,对比研究了体积力的变化规律与主要影响因素.结果表明:等离子体气动激励可以增大气流速度,但随着气流速度的增大,等离子体气动激励的加速效果减弱;激励电压或激励频率增大,体积力均表现为线性增加,但激励电压增大时可以更好的增大体积力;与增大电极内间距相比,增大下层电极宽度可以产生更大的体积力.      

一种基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

本文发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型。该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在二者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解。这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率。该模型先得到了Levine and Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响。