混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

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

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

基于光流信息的圆弧偏置比例导引规律

提出了一种改进的用于被动寻的导弹的基于光流带落角约束的导引规律,其具有更强的适应性和更小的落角误差,而且该导引律不要求弹目距离信息,克服了被动寻的导弹不能测距的约束.基于采用光学传感器和光流算法的测量模型,受昆虫导航的启发,利用光流信息进行被动寻的导弹的导引和控制.同时借鉴偏置比例导引律的结构,重新设计偏置项,实现导弹以期望的落角命中目标.仿真表明,拦截地面运动目标时,该导引律在保证小脱靶量的同时,基本达到期望的落角,对于测量噪声也具有较好的鲁棒性.      

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

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

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.      

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

湍流数据同化技术及应用

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

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

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

Topological Studies of Three-dimensional Flows in a High Pressure Compressor Stator Blade Row without and with Boundary Layer Aspiration

This paper presents a numerical study of the flow topologies of three-dimensional (3D) flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blade is representative of the first stage operating under transonic inlet conditions and the blade design encourages development of highly complex 3D flows. The blade has a small tip clearance. The computational fluid dynamics (CFD) studies show progressive increase of hub corner stall with the increase in incidence. Aspiration is implemented on the hub wall via a slot in the corner between the hub wall and the suction surface. The CFD studies show aspiration to be sensitive to the suction flow rate; lower rate leads to very complex flow structures and increased level of losses whereas higher rate renders aspiration effective for control of hub corner separation. The flow topologies are studied by trace of skin friction lines on the walls. The nature of flow can be explained by the topological rules of closed separation. Furthermore, a deeper analysis is done for a particular case with advanced criterion to test the non-degeneracy of critical points in the flow field.     

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA（W）-1多段翼型的大迎角分离流动及其气动特性。