MHD控制激波诱导边界层分离的机理

激波边界层相互作用是高超声速飞行器面临的重大问题,激波入射到平板引起的边界层分离是其中最具代表性的一种。用加权的3阶ENN格式计算了小磁雷诺数近似的MHD方程,研究了MHD控制层流边界层分离的机理。数值结果显示,通过局部电离空气并施加洛伦兹力,能使分离点向下游移动,分离区尺寸减小,从而抑制和缓解由于激波-边界层相互作用而引起的分离。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

洛仑兹力控制高超声速进气道边界层分离的数值模拟

采用空间HUE格式、时间LU-SGS推进、sst-kw湍流模型、多块结构网格程序,对磁流体动力学(Magnetohydrodynamic:MHD)控制高超声速二维进气道边界层分离进行了数值研究.研究发现,不施加控制时,数值模拟得到的壁面静压和实验结果符合良好,进气道喉道处分离区占据喉道高度的1/3左右.通过施加MHD控制,消除了进气道内部的边界层分离,总压恢复系数从0.502提高到0.56,喉道处流场畸变系数减小了18.6%.     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

基于DDES的跨声速导叶中激波与边界层干涉机理研究

随着透平负荷的逐渐提高,为了分析导叶中出现的激波与边界层干涉的流动机理,需要使用高精度、高分辨率的求解器和先进的湍流计算方法以获取精细的流场结构。采用自主开发的混合RANS/LES求解器,对高负荷透平环境中激波与边界层干涉的流动现象进行研究。结果表明:跨声速流场中同时存在强激波和弱激波,弱激波可在作用位置附近导致熵生成率增大;强激波会诱导边界层分离,而分离区的存在使激波出现马赫反射现象。进一步对湍动能和流动损失进行了机理研究,发现在分离泡前湍动能较大,且湍动能的生成和对流占主导机制;在强激波后由强激波引起的熵生成率较大,且强激波前后区域流动损失主导因素不同。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

确定低雷诺数翼型转捩分离泡位置的实验研究

在翼型模型的表面粘贴表面热膜,由其给出脉动电压的均方根值和波形图,可测出层流边界层分离点和湍流边界层再附点,转捩分离泡的位置也就确定了。     

NACA0012翼型低雷诺数绕流的实验研究

通过水槽氧气泡流动显示和PIV测速实验研究了NACA0012翼型在雷诺数为8200时的流动特性,重点炎注了翼型绕流结构随迎角的变化。研究发脱：分离点和分离翦切层形成旋涡的位置随迎角的增大而向上游移动,同时翼型上表面流动分离后形成的回流区尺寸随着翼利迎角的增加而增大。当流动再附于翼型上表面时,在再附点附近能够观测到展向涡的三维演化过程,并能观测到展向涡的局部配对现象。     

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.      

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

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

基于浸入式边界法的叶栅颤振数值模拟

基于浸入式边界法建立了求解振荡叶栅非定常流动的快速计算模型,并在模型中添加湍流模型,使之能够考虑更加接近实际流动的情况。具体的方法是求解目前为止最常用的雷诺时均Navier-Stokes方程并添加低雷诺数湍流模型(Lam-Bremhorst模型,k-ε模型的一种)进行计算,为了验证该方法的正确性,分别对层流边界层和湍流边界层进行了数值模拟,计算结果与布拉修斯解和壁面率吻合得很好,证明了湍流模型的可靠性。在此基础上,对高雷诺数条件下的振荡叶栅进行了数值模拟。结果表明,折合速度是影响叶栅振荡的重要因素,这与层流计算的结论类似。值得注意的是耦合过程没有生成贴体网格,减小了计算时间,可以准确快速地模拟真实的叶栅流动情况。     

Large-eddy simulation of shock-wave/turbulent boundary layer interaction with and without SparkJet control

The efficiency and mechanism of an active control device ‘‘Spark Jet" and its application in shock-induced separation control are studied using large-eddy simulation in this paper.The base flow is the interaction of an oblique shock-wave generated by 8° wedge and a spatially-developing Ma = 2.3 turbulent boundary layer.The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is20000.The detailed numerical approaches were presented.The inflow turbulence was generated using the digital filter method to avoid artificial temporal or streamwise periodicity.The numerical results including velocity profile,Reynolds stress profile,skin friction,and wall pressure were systematically validated against the available wind tunnel particle image velocimetry(PIV) measurements of the same flow condition.Further study on the control of flow separation due to the strong shock-viscous interaction using an active control actuator ‘‘Spark Jet" was conducted.The single-pulsed characteristic of the device was obtained and compared with the experiment.Both instantaneous and time-averaged flow fields have shown that the jet flow issuing from the actuator cavity enhances the flow mixing inside the boundary layer,making the boundary layer more resistant to flow separation.Skin friction coefficient distribution shows that the separation bubble length is reduced by about 35% with control exerted.     

Experimental investigation on aero-heating of rudder shaft within laminar/turbulent hypersonic boundary layers

The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m~(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test results is also validated by numerical simulations. Our analysis suggests that this phenomenon is due to the difference of boundary layer velocities caused by different thickness of boundary layer between laminar and turbulent flows, as well as the restricted flow within the rudder gap. When the turbulent boundary layer is more than three times thicker than that of the laminar boundary layer, the heat flux of the rudder shaft under the laminar flow condition is higher than that under the turbulent flow condition. Discovery of this phenomenon has great importance for guiding the design of the thermal protection system for the rudder shaft of hypersonic vehicles.