圆转矩形喷管射流掺混特性试验研究

本文采用热线风速仪对圆转矩形收敛喷管的射流流场与掺混特性进行了试验研究,测量了喷口下游不同截面上,射流宽、窄对称面上射流的速度、雷诺剪应力、湍流强度,得到其沿径向和轴向的分布规律。发现在径向方向上,雷诺剪应力先增大后减小,最大值出现在射流与外流的交界面附近,湍流强度则逐渐减小;在轴向上,雷诺剪应力在近喷口区存在小幅波动,而后逐渐减小,湍流度则先略有增大而后变化逐渐减慢。射流特性在宽对称面与窄对称面上的分布规律相同。     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

减阻表面换热特性的实验

采用实验方法对带"V"型沟槽表面通道在不同雷诺数下的流动与换热特性进行了研究.研究结果表明:宽与深均为0.5 mm的沟槽板通道在实验范围内具有明显的减阻效果, 减阻效果随着雷诺数的增大而减小;相比平板通道流动, 低雷诺数范围下(10000雷诺数范围下(30000雷诺数的增加而增加.      

基于大涡模拟的圆转矩喷管尾喷流强化掺混机理研究

为揭示2元圆转矩喷管尾喷流强化掺混的内在机制,应用大涡模拟(LES)方程对2种相同进、出口直径的喷管模型(轴对称、2元圆转矩)在Ma=0.8、高雷诺数(2×10~5)条件下进行了数值模拟计算。结果表明:与轴对称喷管相比,圆转矩喷管射流掺混效应增强,速度衰减快,核心区长度和高温区域面积减小。同时尾喷流拟序结构变化说明:2种喷管主要拟序结构均包含涡环、涡辫、发卡涡、螺旋涡等相似结构;但圆转矩喷管在射流近场诱导出的涡旋更丰富,边角剪切涡发展更快,形成明显的CVP结构,导致其射流柱失稳时刻提前、距离缩短;同时,喷管形式的改变使得射流剪切层内雷诺应力增大,速度脉动增强。拟序结构发展及雷诺剪切应力变化说明在射流流场中涡旋发展耗散速度增大、速度边界层脉动增强、射流柱易失稳是导致射流掺混增强的本质因素,为异形喷管的强化掺混机理提供了依据。     

圆转矩形喷管出口宽高比对射流雷诺剪应力的影响

采用热线风速仪,利用单斜丝,对宽高比W/H分别为1,4,8,12,16的5个圆转矩形收敛喷管和一个轴对称喷管的射流对称面上雷诺剪应力分布特性进行了实验研究.研究发现:在喷口下游不同截面上,射流宽、窄对称面上的雷诺剪应力沿径向均先缓慢增大,到达射流边界后迅速减小,射流边界逐渐沿径向外移.矩形喷管射流相比轴对称射流具有较强的旋流,雷诺剪应力较大,且随着宽高比增大,旋流强度增大,剪应力也逐渐提高,导致了射流与外流掺混增强.宽高比大于8以后,增大幅度逐渐减小.射流宽、窄对称面上的分布规律相同.     

硅微陀螺平板结构的压膜阻尼研究

压膜阻尼对硅微陀螺的精度具有重要影响。平板结构是硅微陀螺中一种典型结构,本文通过对雷诺方程的推导,得出平板在正弦规律小振幅刚性运动时的阻尼力解析式,并给出低频时阻尼系数的近似算法。利用Matlab和CoventorWare软件对运动平板的阻尼进行仿真计算,得出运动平板的阻尼系数和刚度系数,并分析平板参数对阻尼系数及刚度系数的影响。建立不同类型的带孔平板模型,通过比较可以看出,当阻尼孔按一定方式排列时,可以降低结构的阻尼。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

轴流风机中雷诺应力与应变滞后涡黏模型 的PIV实验

采用粒子成像测速(PIV)技术,精细测量了低速轴流风机内部叶片60%跨度处压力面与吸力面的切向和轴向速度,以及风机出口径向平面的切向和径向速度,基于叶片的周缘速度和叶片弦长的雷诺数为370000.应用周期相位平均技术分析雷诺剪切应力与平均速度剪切应变的时空演化过程,并应用互相关技术研究了两者之间的相位关系.结果表明雷诺剪切应力与平均速度剪切应变的极值发生在不同时空相位,存在相位滞后.这一现象广泛存在于三维非平衡复杂湍流中,显然,工业领域广泛应用的线性k-ε模型具有一定的缺陷,应考虑雷诺应力张量与平均速度应变张量之间的时空滞后现象,从而得到符合物理规律的计算结果.因此,雷诺应力张量与平均速度应变张量滞后的涡黏模型将成为一个很有发展前景的封闭模型,来更加精准地预测工业领域中广泛存在的非平衡复杂湍流.      

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...