内外函混合加力燃烧室扩压器内液雾湍流两相流场和浓度场数值模拟

 本文用κ-ε湍流模型和液雾轨道模型,对二维内外函混合加力燃烧室扩压器内有蒸发的液雾两相流动和混合进行了数值预报。其中对气相速度场和温度场的计算结果与实验符合令人满意,液雾两相流的计算结果在趋势上也是合理的,表明本文的模拟方法可用于加力燃烧室性能估算及优化设计。     

扫频式射流对涡轮叶栅间隙泄漏流动影响的数值研究

为控制涡轮叶栅中叶顶间隙泄漏流动和改善涡轮气动性能,将扫频式射流器(SJA)作为一种主动流动控制方法应用在涡轮叶栅的研究中。通过非定常数值计算,分析了SJA对涡轮叶栅叶顶间隙流动的作用过程以及作用机理,并且研究了不同工况下SJA对涡轮叶顶流场改善效果以及不同频率的SJA对叶顶流场的影响。结果表明:通过在涡轮叶栅上端壁增加单个SJA装置,可以有效地延迟上端壁的流动分离,其中最佳方案射流流量仅为进口总流量的0.35%,涡轮叶栅出口截面总压损失系数减少了11.48%。存在着最佳的频率284Hz,使SJA装置对流场的作用效果最佳,有效地改善了涡轮叶栅内的间隙流动。     

圆柱-平板角区湍流流动控制的风洞试验研究

角区流动中马蹄涡系的存在通常会造成不良影响。对圆柱-平板角区流动,在圆柱上游放置一倾斜的小圆棒能够改变角区流动结构。利用油流法和平板表面压力测量方法探讨了湍流流态下不同的小圆棒对平板表面的摩擦力线和压力分布的影响。油流实验揭示了倾斜棒能够改变角区的三维分离,新的分离线由倾斜棒和圆柱共同作用引起;倾斜棒对角区的作用可归类为两种不同的物理现象;倾斜棒能够引起圆柱侧面的分离线向下游发生极大的迁移,导致圆柱底部区域尾迹变窄。平板表面压力测量实验揭示附加的倾斜棒能够极大地改变压力分布情况,角区的逆压梯度相应减小;由此,逆压梯度引起的三维分离必然被削弱。     

湍流数据同化技术及应用

近年来数据同化（DA）被引入湍流动力学研究中,通过融合实验测量和数值计算,提高了实验测量的精度和广度,改善了数值模拟的预测性能。实验观测、预测模型和同化算法是数据同化的三要素,湍流研究中的实验观测包括热线风速仪、粒子图像测速法（PIV）、压力传感器等局部测量数据,预测模型主要指流动控制方程及湍流封闭方程,而同化算法包括贝叶斯推断、集合卡尔曼滤波（EnKF）、伴随等。稳态数据同化一般结合雷诺平均Navier-Stokes （RANS）模型方程,从重新标定模型常数、修正涡黏模型方程形式误差、修正雷诺应力项等方面着手。非稳态的数据同化一般包括四维变分（4DVar）等时间连续的数据同化方式以及顺序数据同化。4DVar通过时间正向和逆向积分迭代,存储量和计算量都非常大。顺序数据同化不需要时间逆向积分,可以在若干时刻对实验观测进行间断地植入,正向求解整个系统。另外,随着人工智能的飞速发展,湍流数据同化研究也向智能化迈进。对于纯数据驱动的湍流机器学习,其缺乏物理本质的约束,而基于物理信息的机器学习在物理本质上与数据同化是一致的。     

快响应PSP技术用于高超声速圆柱绕流的非定常压力测量

压敏涂料(PSP)技术是飞行器风洞实验大面积定量测压和流动显示的重要工具。为了发展非定常流动压力测量和脉冲风洞中全局压力测量能力,航天空气动力技术研究院开发了快速响应压敏涂料技术,与中科院化学所共同合作开发的压敏涂料采用PtTFPP作为发光基团,稳定性较强,持续光照下发光强度衰减为1.5%/h。采用自主研制的静态标定设备在标定腔内测量不同温度和压力条件下20mm×20mm涂料样片的表面发光强度来确定涂料的压力响应特性和Stern-Volmer公式,并设计制作了2套快响应时间动态标定设备,测得涂料典型响应时间在0.2ms。在中国航天空气动力技术研究院FD-03高超声速风洞中对平板圆柱装置进行了Ma=5的 PSP 试验,利用高速相机采集图像数据后,经过批量数据处理,采集频率为250Hz,得到了间隔时间4ms 的连续压力场数据。结合纹影图像和油流图对得到的PSP结果进行了分析,利用同时采集的测压孔数据对PSP 结果进行了比较。试验结果表明,快响应PSP技术可以更详细地显示流场结构,并能同时得到很好的定量压力数据。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

Friction moment analysis of space gyroscope bearing with ribbon cage under ultra-low oscillatory motion

This paper presents the model of calculating the total friction moment of space gyroscope ball bearings which usually work under ultra-low oscillatory motion and are very sensitive to the friction moment. The aim is to know the proportion of the friction moment caused by each frictional source in the bearing＇s total friction moment, which is helpful to optimize the bearing design to deduce the friction moment. In the model, the cage dynamic equations considering six degree-of-freedom and the balls dynamic equations considering two degree-of-freedom were solved.The good trends with different loads between the measured friction moments and computational results prove that the model under constant rate was validated. The computational results show that when the speed was set at 5 r/min, the bearing＇s maximum total friction moment when oscillation occurred was obviously larger than that occurred at a constant rate. At the onset of each oscillatory motion, the proportion of the friction moment caused by cage in the bearing＇s total friction moment was very high, and it increased with the increasing speed. The analyses of different cage thicknesses and different clearances between cage pocket and ball show that smaller thickness and clearance were preferred.     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

低雷诺数下50°后掠三角翼的旋涡流动

采用数值模拟和流动显示的方法研究了50°后掠角三角翼在低雷诺数下的旋涡流动,结果表明:低雷诺数下,非细长三角翼在5°攻角时就形成了稳定的前缘涡,较小攻角时前缘主涡就开始破裂,并观察到泡型和螺旋型两种旋涡破裂方式。另外,在一定的攻角范围内,前缘主涡的外侧又生成一对新的集中涡,构成双涡结构;随着攻角的增大,前缘涡涡核不断升高,主再附线向中心移动,二次分离区扩大。     

尖前缘翼干扰区的壁面压力和热流率分布

本文给出Ｍ＿∞＝７．８和６．７２，Ｒｅ＝３．５×１０￣７／ｍ和５．４×１０￣７／ｍ气流绕迎角为２０°、３０°和３５°尖前缘翼运动时，平板锥型干扰区的壁面压力和热流率分布。结果表明：（１）平板锥型干扰区的特征几何尺度与无粘激波角β＿０和翼迎角α相关，而壁面压力和热流率的峰值与法向马赫数Ｍ＿ｎ相关。（２）翼面压力和热流率分布由于受拐角涡影响，前者在翼根部呈波谷状，而后者呈波峰状，影响尺度与翼前缘处来流边界层厚度有关。