变热线过热比可压缩流湍流度测量方法优化

开展了可压缩流中湍流度测量技术的优化研究，以满足对试验数据高精度评估的需求。在变热线过热比湍流度测量方法推导过程中，忽略了压力脉动项以简化湍流度求解过程。为更加准确评估高速风洞流场湍流度，引入了压力脉动项，以恒温热线风速仪响应关系式为基础，从理论上对可压缩流中湍流度的求解方法进行了优化。在马赫数0.3~0.7进行了湍流度测量试验，并分别利用优化前后的湍流度求解方法对试验数据进行了处理。结果表明两种求解方法所得的湍流度结果量值相近，但优化后的湍流度求解方法所得的湍流度结果随马赫数的变化趋势更加符合客观物理规律。利用蒙特卡洛模拟方法对湍流度的不确定度进行了求解，不确定度量值远小于湍流度量值，表明优化后的湍流度求解方法所得的湍流度结果基本能够代表真实值。试验结果证明了优化后湍流度测量方法的正确性及应用恒温热线风速仪对高速风洞流场湍流度进行测量的可行性。     

噪声对边界层转捩的影响

风洞流场的噪声和湍流度是影响层流边层转捩实验结果的两个重要因素。本文通过实验初步研究了噪声对转捩的影响，实验包括：输入声影响和风洞流场的噪声。结果表明，声频和声压在数值上都存在敏感区；风洞噪声是宽带噪声，以低频为主（至少低速风洞是如此），这对层流边界层转捩有重要作用；风洞噪声会随风速加大而增高，在稳定段用栅格改变湍流度时，实验流场的噪声会随湍流度升高而增加。     

风洞流场扰动和转捩区流场的相关性分析

本文通过实验研究了风洞实验区流场的扰动与边界层转捩区内的速度脉动及其强度、谱特性和两者之间的相关性,提出了实验区流场的扰动对转捩区流场起作用的条件和规律(即转捩对外界扰动响应的条件),这项工作是研究修正风洞背景流场噪声和湍流度对转捩雷诺数影响的基础。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

压缩性修正对γ-Re_θ转捩模型的影响研究

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

可压缩流体热线探针校准方法研究

在可压缩流体中利用热线技术进行湍流度测量时,其输出不仅与脉动速度有关,而且也与流体温度、密度紧密相关。因此,需要在与高速可压缩流体特征相似的校准装置中,在使用前对热线探针进行准确校准。本文提出了可压缩流体热线探针校准方法,建立了热线校准(工作)数学模型,并利用自建的热线校准风洞开展了热线校准实验和湍流度测量风洞试验验证,结果表明:热线校准方法可行,校准数学模型合理可靠,热线探针校准精度较高,湍流度测量结果可信,基本可满足高速可压缩流体湍流度测量的工程需要。     

可压缩流湍流度变热线过热比测量方法

开展了可压缩流中湍流度测量技术的研究,以满足高速风洞高精度试验能力的需求。以对流换热规律为基础,从理论上对可压缩流中热线金属丝热平衡关系式进行了推导,以此为基础,详细推导了恒温热线风速仪的响应关系式,得到了质量流量和总温灵敏度系数的显式表达式,建立了可压缩流中湍流度的求解方法。在马赫数为0.3~0.6范围内进行了湍流度测量试验,以响应关系式为数学模型,利用双曲线拟合方法对试验数据进行了拟合分析,求解得到了马赫数在0.3~0.6范围内流场湍流度约为0.3%~0.6%。对热线输出电压进行了频谱分析,根据频谱特性,利用低通滤波对频域信号进行了处理,有效降低了时域信号脉冲尖峰对湍流度求解的影响,滤波后求解得到马赫数在0.3~0.6范围内流场湍流度约为0.1%~0.3%,与前期测量结果相符。试验结果证明了所建立理论方法的正确性及利用恒温热线风速仪变过热比方法测量可压缩流湍流度的可行性。     

高对流Mach数三维混合层转捩特性分析及小激波结构模拟

采用空间大涡模拟方法对超/超混合、超/亚混合两类三维可压缩平面混合层转捩及其全湍流流场进行了研究,认为混合层动量厚度饱和点可作为流场转捩完成的标志.计算所得到的线性扰动波激励下流场转捩拟序结构与随机扰动下自然失稳结构以及文献结果进行了比对,其结果是一致的,表明了引入线性扰动激励来研究流场转捩结构是合理的.同时,本文还在较高对流Mach数流动下得到了三维流场动态小激波结构,其分布具有非对称特性,且形状与实验及直接数值模拟结果相似.不同条件混合层转捩计算表明:高对流Mach数下混合层转捩以Λ涡结构的形成和发展为主导机制,受扰动及对流条件的影响Λ涡结构不尽相同,某些情况下流场出现二维与三维涡结构共存现象.充分发展湍流区域,流场脉动速度分量量级相同,湍流压缩效应随着对流Mach数提高而明显增强.     

跨声速流场扰动模态与湍流度精细测量

开展了跨声速风洞可压流场扰动模态与湍流度精细测量研究，以满足高速飞行器的低湍流试验需求。基于变热线过热比测量方法，在理论上推导了跨声速流场扰动的一般模态和3种特殊模态特征方程（涡模态、熵模态和声模态），以及扰动模态对应的特征曲线，以此为基础，通过试验获得了流场一般模态、涡模态、声模态扰动图和较高精度的湍流度值，通过扰动图频域分析了马赫数1.50的流场声模态扰动机理，建立了一种跨声速可压流场扰动模态分析与低湍流度测量方法。在新型连续式跨声速风洞完成马赫数0.20～1.50流场测量试验，测得湍流度为0.037%～0.197%，数据非线性拟合优度为0.943～0.995，蒙特卡洛模拟不确定度为0.000 2%～0.004 1%。研究结果证明了所建立的可压流场扰动模态与低湍流度测量方法的有效性，对高性能跨声速风洞的流场评估、优化和飞行器低湍流试验有实际意义。     

Numerical study of compression corner flowfield using Gao-Yong turbulence model

A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM (advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.      

来流湍流度对超高负荷低压涡轮附面层流动控制手段的影响

实验研究了表面粗糙度耦合上游尾迹的流动控制技术,分析了来流湍流度（FSTI）在流动控制过程中对叶片吸力面附面层分离、转捩特性的影响.实验发现:在速度峰值点至分离点之间布置粗糙高度与弦长之比为1.05×10-4的粗糙条带可以在来流湍流度为0.4%与2.2%的低雷诺数范围内降低叶型损失.在雷诺数为85000的状态下,FSTI影响了尾迹通过区、尾迹诱导转捩区及自然转捩区的附面层动量厚度,造成了叶型损失的差异,但FSTI对抑制区的影响较小.      

高超声速钝头体边界层转捩试验

为研究高超声钝头体边界层转捩以及头部钝度对转捩的影响,在FD-14和FD-14A两座激波风洞中开展了热流、压力扰动和高速纹影显示等综合测量。试验结果表明,转捩雷诺数关于钝度雷诺数的变化显示出转捩反转的趋势。压力扰动的功率谱密度（PSD）分析结果以流向离散分布云图形式显示,边界层高速纹影图像显示了第二模态波的发展、湍流的生成和熵层对边界层结构的显著影响。大头部钝度带来的强熵梯度熵层流动对边界层压力扰动频谱特性和流动结构影响显著,在转捩反转机理中起到重要作用。此外,马赫数对转捩的影响不容忽视。     

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

湍流度对边界层转捩测量的影响

通过实验研究了气流湍流度对边界层转捩各种测量方法的影响,指出:随机信号均方根值法是检测边界层转捩较有效的方法。在转捩过程中随机信号均方根值变化较大且稳定。此外还研究了气流湍流度和不同频率的声波扰动对转捩区速度脉动均方根值分布的影响,指出:存在最大脉动均方根值,其位置在0.2δ高度附近,且不受湍流度和声波扰动影响。     

静叶时序对压气机叶片附面层流动影响的数值研究

采用数值方法对两级低速压气机中径处的非定常流场进行模拟,针对压气机第2排静叶两个典型周向位置对动静叶干扰下的叶片附面层流动进行研究.建立尾迹与附面层干扰分析模型,结合叶片壁面摩擦力和近壁面附面层湍动能,详细分析了尾迹和势流干扰下静叶时序改变对叶片附面层流动产生的影响.对第2排静叶附面层的研究结果表明:静叶时序改变了尾迹在其叶排中的输运特征,能够降低壁面摩擦力和近壁面湍动能及其非定常最大波动幅值,影响吸力面附面层内动叶尾迹后沉寂区的宽度.在非定常条件下,尾迹能够诱导静叶层流附面层在尾迹干扰的局部范围内转捩发展为湍流状态,同时高湍流度尾迹的干扰具有抑制逆压梯度下附面层分离的作用,并能够延长层流区的范围.     

扰动来流下声控机翼分离流动的研究

本文介绍在扰动来流下声控机翼流动的试验装置、低速大攻角时的声控效果以及附面层中流动参数的变化规律,并由此得出：随着来流湍流度的增加,因声激励而产生的附面层中大尺度有序结构涡（ＬＳＯＳＶ）的脉动速度幅值Ｖｌ、湍流度增量ｕ′ａ／Ｕ∞以及升力系数增益△ＣＬ均将迅速减小,而激励声压级有效阈值ＳＰＬｅｆｆ则随之逐渐增大。周期性脉动来流在沿机翼附面层流动中急剧衰减,并对机翼分离流动（包括有或无声控）基本上没有什么影响。     

伴流速度对平行喷口射流影响的数值研究

采用大涡模拟方法计算研究平行喷口出口马赫数为0.9,伴流速度比分别为0.1、0.3和0.5时的喷射流场特性。计算时采用高精度的数值模拟方法,并结合Smagorinsky亚网格尺度模型。考察了流场统计平均特性、脉动特性以及射流流场中涡结构的发展演变过程,结果表明:伴流速度的增大使得势流核长度变长,减缓了空间剪切层的发展,转捩延迟。喷射流场速度分布具有自相似性,而湍流强度的分布则不具有相似性。通过分析剪切层中轴向速度脉动、径向速度脉动、压力脉动在空间任意两点上的时空相关性,发现随着伴流速度的增大,脉动量在空间上的相关性减弱,而脉动量向下游的传递速率增加。该研究结果为进一步揭示伴流速度对喷流声场的影响提供基础。     

翼型尾缘噪声源空间分布与辐射特性关系

详细研究了翼型湍流边界层尾缘宽频噪声源空间分布与辐射特性的关系.采用基于雷诺平均流场的翼型尾缘宽频预测方法研究了NACA0012翼型湍流边界层尾缘宽频噪声在4种不同工况下的噪声源空间分布与辐射特性.首先计算了NACA0012翼型湍流边界层尾缘噪声源在不同频率下的空间分布.计算结果发现:边界层中湍流是翼型湍流边界层尾缘噪声声源.随着频率的增加,噪声源强度和噪声源空间尺寸都是先增加后减小,噪声源位置不断靠近翼型尾缘.同时也计算了边界层内不同位置处的噪声源对远场噪声的辐射特性,结果表明:边界层内层区域,其噪声频谱能量集中在高频;边界层外部区域,其噪声能量集中在中低频;攻角增大或者来流速度减小,噪声能量向低频转移.