来流条件对叶片表面附面层转捩影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明:来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。      

来流条件对叶片表面附面层转捩 影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。     

前体涡诱导机翼摇滚的人工转捩技术研究

针对翼身组合体由前体绕流主控的机翼摇滚运动现象,在常规低速风洞中为了模拟其前体流动的过临界Re数绕流流态,研究并提出了一种分布式多转捩丝人工转捩技术。该技术是在模型前体的两侧分别以模型周向角±56.25°为中心周向对称分布多对转捩丝,其间隔为22.5°。转捩丝的轴向位置是从x/D=0.306开始延伸到机翼前缘翼根所在截面。在这样的布局下,当模型摇滚时,前体两侧绕流的转捩有效区内就能始终保持至少各有一根转捩丝。这样,在常规低速风洞中进行模型滚转振荡实验时,利用上述人工转捩技术,只要摇滚的减缩频率足够低,在任一滚转角下就能使得前体的边界层流动转捩到过临界流态。最后,在北航D4常规低速风洞中通过该人工转捩技术,得到了在前体过临界Re数绕流流态下,翼身组合体的机翼摇滚运动形态。     

压气机叶栅叶片表面附面层流态变化影响因素探讨

以平面叶栅中的二元叶栅模型为试验对象,测量了在不同来流条件下叶片表面流场分布情况及栅后气流参数,分析了不同来流条件下叶片表面附面层流动状态的变化。并借助数值模拟手段重点研究了在不同来流马赫数和冲角下,叶片表面压力梯度对层流附面层向紊流附面层转捩过程的影响,通过利用实验数据分析研究来流条件对转捩过程的影响,为从机理上更深刻地认识叶片表面粘性附面层转捩机制提供了科学参考依据。     

超高负荷低压涡轮非定常特性实验

针对某后加载超高负荷低压涡轮叶型的定常与非定常气动特性进行了实验研究.结果表明:由于上游尾迹的作用,在低来流雷诺数与来流湍流度下,尾迹诱导转捩叠加自然转捩可以抑制分离、降低叶型损失.在一个尾迹通过周期,由于尾迹诱导转捩的作用,附面层分离周期性地放大、缩小.同时在吸力面速度峰值点下游观测到了尾迹被割裂为两部分的现象,形成了一个主尾迹与一个副尾迹.副尾迹同样可以诱导转捩,但强度较低,对分离抑制效果有限.      

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

CKDO一方程模型计算高超声速横流转捩

高超声速横流转捩集湍流-激波作用、可压缩、横流、转捩等难题于一身，基于雷诺平均N-S方程（RANS）的转捩建模是一项挑战。本文采用CKDO-tran计算高超声速横流转捩的经典标模—HIFiRE-5，分可压缩性效应、雷诺数效应和迎角效应进行评估，发现：不考虑可压缩效应的KDO-tran模型无法给出准确的预测，而不引入特定转捩机理，CKDO-tran仍然可以预测一系列雷诺数工况的高超声速横流转捩。迎角0°时，CKDO-tran可以较好地预测出双肺叶转捩图像，迎角为4°时，CKDO-tran产生的转捩线提前，但给出了与实验相似的转捩图像，具有进一步开发的潜力。对4°工况进行来流湍流度敏感性的研究发现，随着来流湍流度减小，CKDO-tran产生的转捩线逐渐后移，并且转捩图像与实验的相似性逐渐减弱直至消失。通过对迎角效应进行的分析，推测出来流-激波相互作用的建模是消除迎角效应的关键。     

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前。因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ωSST两方程湍流模式。采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力。     

受周向扰动的涡环的直接数值模拟

采用作者自己发展的柱坐标系中的谱方法对受周向扰动后的涡环的非线性演变过程以及向湍流转捩的初期行为进行了直接数值模拟。结果表明，在自诱导与整体应变场的作用下涡线逐渐前倾，与z轴成45°角，并逐渐拉伸，扰动增长的结果是小尺度运动以能量级串方式逐渐生长，最终向湍流转捩；与稳定性预测一致，扰动波不沿周向传播。     

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

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

Accurately predicting hypersonic transitional flow on cone via a symmetry approach

A new algebraic transition model is proposed based on a Structural Ensemble Dynamics (SED) theory of wall turbulence, for accurately predicting the hypersonic flow heat transfer on cone. The model defines the eddy viscosity in terms of a two-dimensional multi-regime distribution of a Stress Length (SL) function, and hence is named as SED-SL. This paper presents clear evidence of precise predictions of transition onset location and peak heat flux of a wide range of hypersonic Transitional Boundary Layers (TrBL) around straight cone at zero incidence, to an unprecedented accuracy as validated by over 70 measurements for varying five crucial influential factors (Mach number, temperature ratio, cone half angle, nose Reynolds number and noise level). The results demonstrate the universality of the postulated multi-regime similarity structure, in characterizing not only the spatial non-uniform distribution of the eddy viscosity in hypersonic TrBL on cone, but also the dependence of the transition onset location on the five influential factors. The latter yields a novel correlation formula for transition center Reynolds number which takes similar functional form as the SL function within the symmetry approach. It is concluded that the SED-SL model simulates TrBL around cone with uniformly high accuracy, and then points out to an optimistic alternative way to construct hypersonic transition model.     

电弧加热流场湍流度对尖锥边界层转捩影响的研究

电弧加热流场的热环境特性直接影响热防护系统的地面试验数据,由于电弧加热器高温气流和参数波动的原因,直接测量湍流度非常困难。为研究电弧加热流场湍流度对于边界层转捩的影响,采用红外热图热像仪,在电弧加热流场中进行了5°尖锥模型边界层转捩研究。结合数值计算,将试验结果与常规风洞的尖锥边界层转捩结果进行了比较。结果表明：马赫数影响的雷诺数转捩判别准则可以用于计算电弧加热流场的转捩雷诺数;电弧加热流场的尖锥边界层转捩雷诺数显著小于常规风洞的转捩雷诺数,表明在该试验条件下,电弧加热流场的湍流度显著大于常规风洞。     

圆锥高超声速湍流边界层转捩试验研究

主要介绍了一种用于高超声速湍流边界层表面参数测量的试验方法,利用边界层的表面热流值,作为判断边界层转捩与否的依据,开展边界层攻角效应转捩特性试验研究.试验在我院FD-20炮风洞中进行,马赫数为6,试验模型是半锥度5°的圆锥模型.开展了对于圆锥模型,头部半径对转捩位置的影响,小攻角情况下对圆锥模型转捩位置的影响,并且获得了在2°和3°攻角情况下圆锥模型的攻角效应转捩全局曲线.试验结果证实了,在小攻角下,随着攻角的增加,边界层迎风面转捩位置向后移动,背风面转捩位置向前移动.     

高超声速边界层转捩对旋转钝锥自由飞运动的影响

通过在钝锥模型表面上布置人工绊线促使边界层强迫转捩,采用运动自由度不受约束的风洞模型自由飞试验技术研究边界层转捩对高超声速旋转钝锥自由飞行运动特性和气动特性的影响规律,并与自然转捩的旋转钝锥风洞模型自由飞试验结果作对比分析,试验马赫数为5.0,以模型长为特征尺寸的自由流雷诺数为1.68×106。研究结果表明:有人工绊线的旋转钝锥在自由飞行过程中有"激励稳定"的绕流流场,产生动态稳定的自由飞运动(动稳定导数系数小于0),而无转捩绊线的旋转钝锥在自由飞行中则有"激励不稳定"的绕流流场,产生动态不稳定的自由飞运动(动稳定导数系数大于0)。     

激波风洞边界层转捩测量技术及应用

高超声速边界层转捩对摩阻、传热等有重要影响。在高超声速飞行器研制中，迫切希望能精确预测和控制边界层转捩。激波风洞作为高超声速气动热环境试验的主要地面模拟设备，是研究高超声速边界层转捩的重要设备。但激波风洞原有测量技术适用于工程型号试验，需要依据高超声速边界层转捩特点进行适应性改造和升级。依据高超声速边界层转捩过程中的热流、压力、密度等物理参数变化，发展了薄膜热流传感器测热技术、温敏热图测量技术、高频脉动压力测量技术、高清晰度纹影显示技术等适用于激波风洞的边界层转捩测量技术。并针对头部钝度0.05 mm的半锥角7°尖锥模型，在中国空气动力研究与发展中心Ø2 m激波风洞（FD-14A）马赫数10、单位雷诺数1.2×10⁷/m的流场条件下开展了边界层转捩试验。采用多种转捩测量技术同时进行测量，获得尖锥模型表面边界层转捩情况、边界层脉动压力频谱特征、边界层内清晰的第2模态波和湍流斑纹影图像，不同测量技术获取的试验结果可相互印证，线性稳定性理论分析结果与试验结果相吻合。     

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.     

激波绕射触发爆震波的试验研究

在单爆震试验平台上对激波绕射触发爆震波的特性进行了试验研究.针对不同的激波反射器结构尺寸,测量了爆震管内气体压力变化历程及火焰传播速度,获得了不同工况下激波绕射触发爆震波的特性.结果表明:激波反射器对爆震波的形成有促进作用,并取决于反射器中的反射锥及反射孔板的结构与尺寸;在反射锥堵塞比一定的情况下,存在最佳反射锥锥角使得爆震波的形成距离最小;而反射锥与反射孔板的间距增大,爆震波形成距离增大,一旦形成稳定的爆震波后,反射器对爆震波压力的影响不大.      

三维高超声速边界层中eN方法的应用

针对e N方法中的鞍点法(SPM),结合内波理论,提出了一种简便算法.并用该算法计算了一个由平板、圆锥、圆柱3部分组成的模型以马赫数为10的速度在高度位于30~45km之间、攻角为0°和10°的工况下N值分布情况.结果表明:高度的增加会使转捩位置向下游移动,转捩面积会变小甚至消失;与攻角为0°相比,攻角为10°会使迎风的平板部分转捩发生的可能性增大,转捩面积也会变大;对于背风处的圆锥部分,攻角为10°却会减弱转捩发生可能性,使得转捩面积变小,甚至消失,而且转捩位置也会变化.     

基于小波变换的高超声速边界层动态图像分析

利用小波变换技术,研究高超声速下动态尖锥周期俯仰摆动的边界层转捩图像。通过基于二维小波变换的图像消噪和图像增强算法,优化动态实验图像,根据小波变换图像的奇异性提取沿时间历程变化的尖锥边界线和边界层图像,得到尖锥的俯仰摆动周期和迎角变化范围。结果显示边界层厚度的变化周期与尖锥的俯仰摆动周期存在相位差,随着迎角增大,背风面转捩点向尖锥前缘移动。实验结果表明,基于小波变换的图像处理技术不仅能够作为高超声速尖锥实验流动显示进行定性观察,还可以实现对高超声速边界层转捩实验的精确观测和定量分析,成为一种重要的量化和自动化的实验技术和分析手段。     

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

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