反射激波作用下重气柱界面演化的PIV研究

在水平激波管中采用PIV方法研究了反射激波作用下SF6重气柱界面的发展演化。采用射流方法形成SF6无膜气柱界面,并以乙二醇作为示踪粒子。利用连续激光片光源结合高速摄影相机对流场进行显示,得到了反射激波作用下SF6气柱界面的发展过程。结果表明,入射激波的冲击会在界面上产生反向旋转的涡环结构,而反射激波的作用会在界面上产生与初始涡环旋转方向相反的次级涡环结构。此外,对反射激波作用后的流场图像进行PIV后处理,获得了流场连续的速度场和涡量场。获得的环量与已有的理论模型进行比较,取得了较好的一致性,也验证了本文实验方法的可行性。     

反射激波作用下两种重气柱界面不稳定性实验研究

在水平方形激波管中对两种无膜重气柱界面（分别是SF6和氩气）在反射激波作用下的不稳定性发展进行了实验研究。气柱界面采用射流技术形成,实验采用连续激光片光源照射流场,乙二醇作为示踪粒子,并用高速摄像机对流场进行拍摄,获得了入射激波以及反射激波共同作用下,两种不同气柱界面的演化过程。实验结果表明,两种气柱的Atwood数不同,界面演化速率不同,反射激波到达前后的界面形态不同。SF6气柱在入射激波作用下会产生两个比较明显的反向的涡环结构,而氩气柱界面上由于产生的涡量较少,涡环结构并不明显。在反射激波作用下,SF6气柱界面会出现明显的次级涡对,而且次级涡对的旋转方向与初始涡环结构的旋转方向相反。对于氩气柱而言,在反射激波作用下虽然也产生了与初始涡环方向相反的次级涡对,但次级涡对始终未充分发展。这是因为反射激波作用时氩气柱界面的Atwood数较小导致氩气柱界面上产生的反向涡量较少。实验结果充分表明了气体Atwood数对界面不稳定性的发展起到了较大的影响。     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

超声速横向射流三维流场结构特征

为了研究超声速燃烧室内燃料与空气快速掺混过程的流场特性，基于可压缩Navier-Stokes方程，采用大涡模拟方法和高精度WENO-TCD混合格式对来流马赫数为2.68，喷压比为36的超声速横向射流流场结构进行数值研究。数值结果清晰描述了超声速主流与横向射流相互作用过程的流场结构特征，得到了三维激波形态的演变规律以及它们在强化混合过程中的作用。另外，因桶形激波背风面低压区处的斜压效应，射流气体在桶形激波背风面形成一对螺旋向上的反向涡对，反向涡对的卷吸作用诱导进入壁面边界层的主流向上运动，形成冲击射流。冲击射流以v=557m/s的法向速度向上冲击桶形激波背风面，因而在桶形激波背风面留下类三角锥面凹痕。     

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

一种用于RM不稳定性研究的竖直环形激波管的设计与验证

设计并加工了一套竖直环形同轴无膜激波管,可用于环形汇聚激波诱导下的Richtmyer-Meshkov不稳定性实验研究。与前人工作相比,本文在流体界面的形成以及流场的观测方法上做了较大的改进。通过实验和数值方法,对该竖直激波管产生的环形柱状汇聚激波的参数进行测量和分析,验证了同轴激波管形成柱状汇聚激波方法的可行性和可靠性。在界面形成方面,采用细丝约束肥皂膜技术形成正八边形气体界面,并利用数值方法考察了细丝对界面发展的影响。结果表明在界面发展的前期,细丝的影响几乎可以忽略。利用连续激光片光结合高速摄影相机对流场进行观测,获得了正八边形air/SF6气体界面在环形汇聚激波及其反射激波冲击下的演化过程,并与数值结果进行了对比,获得了较好的一致性,进一步验证了汇聚激波的对称性以及细丝约束肥皂膜技术用于形成多边形气体界面的可靠性。     

冷气喷射法控制激波强度的数值研究

对冷气喷射时激波受到的影响进行了多方案数值研究.在高压级静叶吸力面反射激波生成点前后5个不同位置上采用相同总压及相同温度的冷气喷射,分析了不同位置的冷气喷射对激波强度和方向影响.结果表明:在吸力面激波折射点附近喷射冷气是一种行之有效的控制激波强度的方法.冷气喷射位置位于激波折射点附近时能够对激波的强度产生影响,在接近激波折射点前部位置注入冷气对流动有积极作用,减弱了激波.      

激波与氦气泡相互作用的实验与数值研究

提出了一种在激波管中实现运动激波与球形气一气界面相互作用的研究方案。以肥皂膜作为气体交界面,高压火花为光源,采用阴影法,在马赫数M=1．2的情况下,在矩形激波管中开展了一系列的实验研究,获得了激波作用下球形氦气泡界面发展变化的相关实验结果。还通过求解二维轴对称的Euler方程,结合网格自适应策略,对实验过程进行了数值模拟。通过与实验结果的对比分析,证实了涡量的产生和分布对于界面的变形和发展有着重要意义。     

射流涡发生器对激波边界层作用诱导的流体分离控制大涡模拟研究

为了研究射流涡发生器对激波边界层作用所诱导的流动分离控制机理及其流场特性,基于大涡模拟(Large eddy simulation)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=2.5情况下,平板上射流涡发生器对激波与边界层相互作用所诱导流场进行了数值模拟。结果表明,射流涡发生器对激波边界层的流体分离有一定的抑制作用,与无控制情况相比,射流作用下进出口总压恢复系数由85.9%提高到94.6%。射流尾涡主要集中于一环状区域内,在该区域内,入射激波与马蹄涡、桶形激波上方的涡管以及剪切涡相互作用,导致整体尾流被激波往下压缩。同时在激波的压缩下,各涡之间相互缠绕、挤压合并,形成多个流向小涡结构,将边界层外的高速流体卷入边界层内,从而增加边界层底层能量,达到抑制流动分离的目的。     

An integration method based on a novel combined flow for aerodynamic configuration of strutjet engine

In this paper a novel design method of aerodynamic configuration is proposed to integrate forebody, strut and inlet for strutjet engine, and a model at design point of Mach number 6 is generated to investigate the aerodynamic performance by both simulations and experiments. The basic flow field employed by proposed method is a combined flow named IBB, which is combined by Internal Conical Flow A (ICFA), truncated Busemann flow I (BI) for external section, and truncated Busemann flow II (BII) for internal section. The model configuration is generated by streamline tracing method from basic flow field, in which the forebody section is traced from ICFA and BI flows, and the inlet as well as strut section is traced from BII flow. The simulations in Mach number 4, 5, and 6 demonstrate uniform starting flow fields with relatively high total pressure recovery, which agree well with experiments in wind tunnel. Additionally, in low Mach number cases, this inlet could start at Mach number 3 while it is unstarted at Mach number 2.7; in high Mach number cases, a uniform flow could still exist in Mach number 6.5 while a relatively strong shock wave boundary layer interaction is found in cowl area of Mach number 7 case, indicating the inlet designed by proposed method works in a relatively wide Mach number range.     

跨音压气机叶栅的激波结构模型及损失

在分析讨论激波与附面层相互作用和栅后背压引起激波形状及强度变化的基础上, 给出了考虑激波/附面层相互作用及栅后背压的跨音叶栅激波结构的物理数学模型。应用本文所提出的模型分析了跨音叶栅的激波损失, 其结果和实验结果一致。激波损失的精确得出, 使得将激波与附面层相互作用所引起的流动分离损失从流动总损失中分离出来成为可能, 有助于了解激波与附面层相互作用引起流动分离的机理。      

流向涡与斜激波相互作用的试验研究和数值模拟

为研究流向涡与斜激波相互作用在超声速燃烧中的应用,进行了由翼产生的流向涡与楔块产生的二维斜激波相互作用的燃烧室冷流试验研究.在不同激波强度下,纹影仪捕捉到了强、中、弱不同的涡/波作用现象.仿真与试验结果符合得较好.试验结果表明:对马赫数为2.3流场,在翼攻角10°时,能产生强流向涡,此工况下,锲角越大,涡/波作用越强.仿真结果表明:马赫数对涡/波作用影响较大,总压影响不明显,总温可影响亚声速回流区的尺寸.      

平面激波绕矩形障碍物后的绕射与反射

本文通过激波管二维光测模型实验和采用矩形网格的FLIC方法进行数值计算两种手段针对平面激波绕矩形障碍物后的绕射与反射的问题进行了分析和研究。文中给出了绕射波的波形以及平面激波绕过矩形障碍物后与地面相互作用而产生的从规则反射到马赫反射的波系图案。本文考虑了入射激波马赫数M_i及不同矩形障碍物高度H的影响,并将二维光测模型实验同数值计算的结果进行了比较。     

Viscosity effects on weak irregular reflection of shock waves in steady flow

The viscosity effects on weak shock wave reflection are investigated with the Navier–Stokes and DSMC flow solvers. It is shown that the viscosity plays a crucial role in the vicinity of three-shock intersection. Instead of a singular triple point, in viscous flow there is a smooth three shock transition zone, where one-dimensional shock jump relations cannot be applied. At the flow parameters corresponding to the von Neumann reflection, when no inviscid three-shock solution exists, the computations predict an irregular shock-wave configuration similar to that observed previously in experiments. The existence of a viscous zone in the region of shock-wave interaction allows a continuous transition from the parameters behind the Mach stem to the parameters behind the reflected shock, which is impossible in the inviscid three-shock theory. Results of numerical simulations agree qualitatively with the conclusions of the approximate viscous theory, which describes the flow in the vicinity of the three-shock intersection.     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.     

爆轰驱动激波管缝合激波马赫数计算

采用数值求解的方法得出了氢氧爆轰驱动激波管的缝合状态参数.以空气为试验气体时,缝合激波马赫数随着H2摩尔浓度的增加而增加,H2摩尔浓度达到90%左右时达到最大.当缝合马赫数较高时,需要考虑高温真实气体效应的影响,缝合激波马赫数较理想气体的高.以氢空气混合气为试验气体时,缝合激波马赫数较以空气为试验气体的小.通过调整驱动气体与被驱动气体的初始参数,可以得到即能恰好消除Taylor波又能缝合的运行状态.     

压气机叶栅中的激波附面层相互作用

压气机叶栅中激波附面层相互作用是十分复杂的问题,由于相互作用引起分离是决定跨音速压气机性能的重要因素之一。然而迄今为止,尚未深入进行过压气机叶栅激波附面层相互作用的研究,发表的文献极少。作者经过计算和分析,说明压气机叶栅流中主要的相互作用形式是叶栅槽道中激波和湍流附面层的相互作用。作者通过分析指出,研究压气机叶栅激波附面层相互作用,不能直接应用Pearcey分离准则。作者并提出了适用于压气机叶栅的分离准则的函数关系为f(M_1,p_1/p_(L.E.),P_(r.E.)/p_2,Re_0)=O。     

AUSM+格式的改进

为了提高AUSM 格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM 格式的基础上经过数值研究分析,对AUSM 格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM 格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM 格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM 格式继承了原AUSM 格式优点,计算处理问题能力增强,可以有效地应用于亚声、跨声和高超声流动的工程实际问题中.     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.