对称翼型近场尾流及后缘边界层流动特性的试验研究

本文简述了用热线风速仪测量“NACA 63-012”对称翼型近场尾流及后缘边界层流动特性,讨论了流场的时均特性及雷诺应力等湍流量的分布特性,提出了无因次速度分布的半经验公式、相关参数及相关曲线解析式,并得出了几点重要结论。     

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

1-3型压电纤维主动叠层板扭转特性的研究

对4类典型铺层的1-3型压电纤维复合材料主动板的变形模式进行了理论分析。在此基础上对具有扭转变形模式的主动叠层板进行了振动特性分析,得到了该类型主动叠层板在外界电场激励下的稳态解。最后对该类主动叠层板的扭转振动特性进行了深入的研究,包括作用电场强度的大小、频率对扭转特性的影响、压电纤维铺层角和压电纤维体积比对主动叠层板扭转特性的影响等。算例表明本文计算模型和有限元得到的结果符合良好,同时算例结果表明含1-3型压电纤维复合材料主动叠层板具有良好的主动扭转特性。     

三维舵面绕流的N-S方程计算研究

基于粘性非结构混合网格,采用格心格式的有限体积法,应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟。操纵面铰链力矩的计算结果与风洞试验结果的对比分析表明,计算方法对粘性流动具有良好的适应性。对舵面铰链力矩随迎角、马赫数、雷诺数变化的规律也进行了分析研究。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

不同当量比下氢气体积分数对甲烷-氢混合气预混火焰燃烧不稳定性的影响

在可调振幅的正弦波声场作用下,以甲烷-氢混合气预混钝体火焰为实验对象,通过火焰传递函数表征整个燃烧系统的燃烧不稳定性特征,借助CH基自发荧光图像描述火焰锋面运动及演化过程,研究了在声激励下不同当量比（0.8、1.0、1.2）下氢气体积分数的变化（0、10%、20%）对火焰燃烧不稳定的影响。结果表明: 对于当量比为0.8和1.0的预混火焰,氢气的加入使得火焰传递函数幅值增大,热释放波动变大,整体火焰的不稳定性增强;对于当量比为1.2的富燃预混火焰,随着氢气体积分数的增加,火焰传递函数幅值先减小后增加,火焰稳定性较强。      

半球谐振陀螺全角模式控制技术研究

介绍了半球谐振陀螺全角模式的工作机理及其技术特点。对全角模式下振型的控制进行了理论分析,采用乘法相干解调对陀螺振型参数进行提取,并采用参数激励方法对谐振子进行驱动和控制。设计了数字电路进行实验,结果表明该全角模式控制技术实现了半球谐振陀螺的角速率积分功能,能够0°～360°全角度敏感输出,测得的角速度测量值超过300(°)/s,角度测量线性度优于10~(-4),与传统的力平衡模式相比,半球谐振陀螺的动态特性得到大幅提升。     

高湍流度时全气膜涡轮叶片表面冷却和换热特性的实验研究

为获得高主流湍流度时全气膜涡轮叶片表面的冷却和换热特性,在跨声速风洞中实验研究了质量流量比(MFR)和主流雷诺数(Re)对叶片表面气膜冷却效率和换热系数比的影响。在叶片前缘布置了5排圆形孔,在吸力面和压力面分别布置了3排和6排圆形孔,实验结果由嵌入在叶片中截面的热电偶测得。实验中基于弦长的主流雷诺数的范围为3.0×105～9.0×105,叶栅出口马赫数Ma为0.8, MFR的范围是5.5%～12.5%,主流湍流度Tu为14.7%。实验结果表明:主流雷诺数升高显著增强了叶片表面的换热,使层流边界层到湍流边界层的转捩位置提前。对于吸力面S/C0.2的区域(S/C为当地弧长与弦长之比),气膜冷却效率受MFR影响明显,当MFR大于7.7%时提高MFR会导致气膜冷却效率降低;该区域的换热系数比在中低雷诺数时受MFR影响较小,在高雷诺数时随MFR升高而升高。压力面S/C-0.7区域的气膜冷却效率随MFR升高而升高,-0.7S/C-0.4区域的气膜冷却效率受MFR影响较小,对于整个压力面而言,MFR升高提高了叶片表面的换热系数。相对于叶片其它区域,压力面后半段区域和吸力面的气膜冷却效率受雷诺数影响较大。     

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。