振动激发对高超声速气动力/热影响

随着飞行马赫数的不断提高,空气的高温气体效应越来越明显,对高超声速飞行器的气动力/热特性产生重要影响.高温气体效应对气动力/热的影响机理复杂,影响参数众多,迄今为止国内外尚未完全研究清楚.发生高温气体效应时,多个非线性物理过程耦合在一起,地面试验和数值模拟无法将这些过程解耦,无法给出关键物理机理.为了解决这一问题,文章提出一种理论分析与数值模拟相结合的两步渐进新方法:先通过牛顿迭代法得到发生振动激发过程的斜激波无黏解;再将该无黏解的结果作为边界条件,求解边界层的黏性解.利用该方法研究了振动激发过程对二维斜劈的气动力/热特性的影响规律.研究结果表明,振动激发过程对斜激波后的温度、密度、马赫数、雷诺数和斜激波角影响较大,而对压力和速度影响较小.斜激波波后的无黏流动与边界层流动是耦合在一起的.发生振动激发后,斜激波波后雷诺数的增大会导致边界层厚度减小,结合多个物理量的变化,如速度增大和温度减小,共同对边界层内的摩擦阻力和气动热产生影响.对比完全气体的结果发现,振动激发使壁面摩阻升高,而使壁面热流降低.分别通过影响激波层和边界层,振动激发对摩阻的影响是弱耦合的,而对热流的影响则是强耦合的.     

六相交流放电再入飞行器热环境地面模拟研究

再入飞行器高速飞行过程中,其表面受到强烈的气动加热作用,所产生的复杂高温气体环境会破坏飞行器材料,影响飞行器结构的可靠性.因此,基于地面装置实现高速飞行器再入过程中表面热环境的模拟,对于再入飞行器的热防护测试具有十分重要的意义.文章基于数值模拟,分析了工作气压的变化对等离子体中非平衡能量输运过程以及等离子体气体温度等参数的影响规律,提出了通过改变工作气压来调节等离子体冲击壁面的热流密度的方法.基于此,首先以表面热流密度和加热时间与真实飞行条件下一致为原则,基于六相交流电弧放电等离子体实验平台,产生了大体积、高气体温度,且壁面热流密度可调的等离子体电弧射流.然后,对采用酚醛浸渍基碳热防护材料的烧蚀体进行了地面烧蚀实验,在壁面热流密度为1.07～3.95 MW/m2范围内获得了与文献报道吻合较好的实验结果,初步验证了该方法的可行性.对高速再入飞行器典型部件进行了烧蚀实验,在壁面最高热流密度为5 MW/m2的实验条件下,获得了与空间飞行实验吻合良好的地面模拟实验结果.这表明在不采用高成本风洞的前提下,本论文所提出的地面模拟实验方法可在一定程度上模拟飞行器再入过程中的表面热环境.     

等离子体激励气动力学探索与展望

等离子体激励气动力学是研究等离子体激励与流动相互作用下, 绕流物体受力和流动特性以及管道内部流动规律的科学, 属于空气动力学、气体动力学与等离子体动力学交叉前沿领域. 等离子体激励是等离子体在电磁场力作用下运动或气体放电产生的压力、温度、物性变化, 对气流施加的一种可控扰动. 局域、非定常等离子体激励作用下, 气流运动状态会发生显著变化, 进而实现气动性能的提升. 国际上对介质阻挡放电等离子体激励、等离子体合成射流激励及其调控附面层、分离流动、含激波流动等开展了大量研究. 等离子体激励调控气流呈现显著的频率耦合效应, 等离子体冲击流动控制是提升调控效果的重要途径. 发展高效能等离子体激励方法, 通过等离子体激励与气流耦合, 激发和利用气流不稳定性, 揭示耦合机理、提升调控效果, 是等离子体激励气动力学未来的发展方向.      

带有中间裂纹载流薄板放电瞬间耦合场的数值模拟

应用耦合场理论采用数值分析的方法计算了具有中间裂纹的导电薄板在放电瞬间裂纹端区域附近的温度场、温度梯度场的分布状态。计算结果表明：由于电流产生的焦耳热源的作用，裂纹尖端处温度瞬时急剧升高，能够在很小的范围内熔化形成微小的口，裂纹前缘的曲率半径增加了几个数量级，显著地减少了应力集中，阻止了干线裂纹源的形成，有效地遏制了裂纹的扩展，数值模拟采用了热－电耦合的分析方法，考虑了材料电阻随温度变化和电流流过生成的内热源间的相互作用，同时考虑了导热系数和比热随温度变化所产生的影响，对于具有中间裂纹的载流反，根据结构、材质、边界条件及通电电流的对称性，计算得到了两个裂纹尖端完全对称的结果。     

磁流体动力学斜激波控制数值模拟分析

高超声速飞行器MHD(磁流体动力学)斜激波控制应用的关键在于理解等离子体斜激波流场与磁场的相互作用规律,这里发展了全MHD数值模拟方法对其进行研究,数值方法基于八波方程附加源项形式,进行有限体积离散,采用了Roe求解器、OC-TVD空间格式和LU-SGS方法,且采用投影方法降低磁场伪散度误差.考察外加均匀磁场的马赫10无粘导电拐角流动,压缩角为10°.结果中散度误差较低,并且通过激波参数验证了结果的准确性.流场显示,磁场使得激波角增大,部分情况下出现了快、慢激波结构,其中快激波变化更明显;壁面压强根据磁场的不同出现了不同程度的降低.最后采用群速度图方法进行了快慢激波形式分析,解释了磁场影响下流场形式变化机理.     

JP-10裂解的激波管实验与动力学模拟

利用单脉冲激波管对碳氢燃料JP-10在1150~1300K条件下的高温热裂解特性进行了实验研究,采用气相色谱法分析热裂解产物并获得了热裂解速率系数.主要裂解产物有乙烯、乙炔、丙烯、丁烯、1,3-丁二烯、环戊二烯、环戊烯、苯、甲苯,以及少量的甲烷、乙烷、二甲苯和甲基环戊烯.将每次激波管实验后所有产物浓度累加, JP-10裂解速率系数由实验测定.为了消除激波运行中非理想性和边界层效应导致反应温度确定的误差,采用对比速率法确定裂解温度,即在反应物中加入少量热解速率已知的内标物,根据内标物在相同的激波管实验条件下的裂解程度确定反应温度.根据内标物裂解量确定的激波管裂解反应温度通常小于采用传统测量激波速度由激波关系计算的反射激波后5区温度.在1200~1300K之间两种方法得到的温度吻合得较好,差异在20K以内,随着温度升高,两者差异增大.在实验研究的基础上,依据San Diego Mechanism对JP-10高温裂解过程进行了动力学模拟.结果显示:主要裂解产物中乙烯、乙炔和1,3-丁二烯产量随温度变化的实验值与San Diego Mechanism的模拟结果有很好的一致性,但环戊烯产量的实验值比模拟值高很多,预示JP-10裂解中完全开环和部分开环反应都是重要的裂解通道.      

JP-10裂解的激波管实验与动力学模拟

利用单脉冲激波管对碳氢燃料JP-10在1150～1300 K条件下的高温热裂解特性进行了实验研究,采用气相色谱法分析热裂解产物并获得了热裂解速率系数.主要裂解产物有乙烯、乙炔、丙烯、丁烯、1,3-丁二烯、环戊二烯、环戊烯、苯、甲苯,以及少量的甲烷、乙烷、二甲苯和甲基环戊烯.将每次激波管实验后所有产物浓度累加, JP-10裂解速率系数由实验测定.为了消除激波运行中非理想性和边界层效应导致反应温度确定的误差,采用对比速率法确定裂解温度,即在反应物中加入少量热解速率已知的内标物,根据内标物在相同的激波管实验条件下的裂解程度确定反应温度.根据内标物裂解量确定的激波管裂解反应温度通常小于采用传统测量激波速度由激波关系计算的反射激波后5区温度.在1200～1300 K之间两种方法得到的温度吻合得较好,差异在20K以内,随着温度升高,两者差异增大.在实验研究的基础上,依据San Diego Mechanism对JP-10高温裂解过程进行了动力学模拟.结果显示:主要裂解产物中乙烯、乙炔和1,3-丁二烯产量随温度变化的实验值与San Diego Mechanism的模拟结果有很好的一致性,但环戊烯产量的实验值比模拟值高很多,预示JP-10裂解中完全开环和部分开环反应都是重要的裂解通道.     

非线性热环境下高温合金蜂窝板隔热性能研究

金属蜂窝板结构在高温热环境下的隔热特性是高速飞行器热防护设计的重要参数. 使用自行研制的高速飞行器瞬态气动热试验模拟系统, 对高温合金蜂窝平板结构在高达800℃的非线性热环境下的隔热性能进行实验研究, 获得了蜂窝板结构的瞬态和稳态传热特性以及在多种不同温度下金蜂窝平板结构隔热效果的实验数据. 在考虑结构内部蜂窝芯壁面间辐射、金属结构的传热以及蜂窝腔内空气传热的多重热交换条件下, 采用三维有限元计算方法对蜂窝板的隔热特性进行了数值模拟, 计算结果和试验结果的吻合性良好, 验证了数值模拟方法的可信性和有效性, 并为数值模拟方法能够在一定程度上较好地替代价格昂贵的气动热模拟试验打下了基础. 讨论了在复杂非线性高温环境下金属蜂窝板隔热效率的变化, 加热面温度的升降速度与隔热效率的关联性以及金属蜂窝板表面发射率的选取等问题, 对高速飞行器金属蜂窝结构的热防护研究具有重要的参考价值.      

可压缩流向涡与激波轴对称干扰的数值模拟1)

用NS方程数值模拟了可压缩流向涡和激波轴对称相互作用现象.数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较.结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征.提出了判断流向涡与运动激波相互作用中旋涡破碎的准则.     

固定网格下的特征线法求解溃坝问题

溃坝问题是典型的非线性双曲方程的Riemann问题,其数值求解的难点在于对间断面的捕捉以及避免间断面处在数值计算过程中产生数值色散,因而为求解此问题所产生的各种数值计算方法的优劣也体现在这两个方面。本文针对溃坝问题提出一种新的计算方法。该方法基于对偶变量推导的浅水波方程,根据方程的特点,从方程的特征值和黎曼不变量出发,采用高精度的激波捕捉方法计算黎曼不变量的位置随时间的变化,然后映射至不随时间变化的固定网格。根据黎曼不变量的位置,采用保形分段三次Hermite插值将物理量映射至网格节点。计算结果显示,该方法不仅操作简单,计算量小,而且结果准确。     

Quantification of the effect of surface heating on shock wave modification by a plasma actuator in a low-density supersonic flow over a flat plate

This paper describes experimental and numerical investigations focused on the shock wave modification induced by a dc glow discharge. The model is a flat plate in a Mach 2 air flow, equipped with a plasma actuator composed of two electrodes. A weakly ionized plasma was created above the plate by generating a glow discharge with a negative dc potential applied to the upstream electrode. The natural flow exhibited a shock wave with a hyperbolic shape. Pitot measurements and ICCD images of the modified flow revealed that when the discharge was ignited, the shock wave angle increased with the discharge current. The spatial distribution of the surface temperature was measured with an IR camera. The surface temperature increased with the current and decreased along the model. The temperature distribution was reproduced experimentally by placing a heating element instead of the active electrode, and numerically by modifying the boundary condition at the model surface. For the same surface temperature, experimental investigations showed that the shock wave angle was lower with the heating element than for the case with the discharge switched on. The results show that surface heating is responsible for roughly 50 % of the shock wave angle increase, meaning that purely plasma effects must also be considered to fully explain the flow modifications observed.     

Plasma mitigation of shock wave: experiments and theory

Two types of plasma spikes, generated by on-board 60 Hz periodic and pulsed dc electric discharges in front of two slightly different wind tunnel models, were used to demonstrate the non-thermal plasma techniques for shock wave mitigation. The experiments were conducted in a Mach 2.5 wind tunnel. (1) In the periodic discharge case, the results show a transformation of the shock from a well-defined attached shock into a highly curved shock structure, which has increased shock angle and also appears in diffused form. As shown in a sequence with increasing discharge intensity, the shock in front of the model moves upstream to become detached with increasing standoff distance from the model and is eliminated near the peak of the discharge. The power measurements exclude the heating effect as a possible cause of the observed shock wave modification. A theory using a cone model as the shock wave generator is presented to explain the observed plasma effect on shock wave. The analysis shows that the plasma generated in front of the model can effectively deflect the incoming flow; such a flow deflection modifies the structure of the shock wave generated by the cone model, as shown by the numerical results, from a conic shape to a curved one. The shock front moves upstream with a larger shock angle, matching well with that observed in the experiment. (2) In the pulsed dc discharge case, hollow cone-shaped plasma that envelops the physical spike of a truncated cone model is produced in the discharge; consequently, the original bow shock is modified to a conical shock, equivalent to reinstating the model into a perfect cone and to increase the body aspect ratio by 70%. A significant wave drag reduction in each discharge is inferred from the pressure measurements; at the discharge maximum, the pressure on the frontal surface of the body decreases by more than 30%, the pressure on the cone surface increases by about 5%, whereas the pressure on the cylinder surface remains unchanged. The energy saving from drag reduction is estimated to make up two-thirds of the energy consumed in the electric discharge for the plasma generation. The measurements also show that the plasma effect on the shock structure lasts much longer than the discharge period.

---

     

The formation of a secondary shock wave behind a shock wave diffracting at a convex corner

This paper deals with the formation of a secondary shock wave behind the shock wave diffracting at a two-dimensional convex corner for incident shock Mach numbers ranging from 1.03 to 1.74 in air. Experiments were carried out using a 60 mm 150 mm shock tube equipped with holographic interferometry. The threshold incident shock wave Mach number () at which a secondary shock wave appeared was found to be = 1.32 at an 81° corner and = 1.33 at a 120° corner. These secondary shock waves are formed due to the existence of a locally supersonic flow behind the diffracting shock wave. Behind the diffracting shock wave, the subsonic flow is accelerated and eventually becomes locally supersonic. A simple unsteady flow analysis revealed that for gases with specific heats ratio the threshold shock wave Mach number was = 1.346. When the value of is less than this, the vortex is formed at the corner without any discontinuous waves accompanying above the slip line. The viscosity was found to be less effective on the threshold of the secondary shock wave, although it attenuated the pressure jump at the secondary shock wave. This is well understood by the consideration of the effect of the wall friction in one-dimensional duct flows. In order to interpret the experimental results a numerical simulation using a shock adaptive unstructured grid Eulerian solver was also carried out. Received 1 May 1996 / Accepted 12 September 1996     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

飞行Ma12条件超燃发动机流场及燃烧特征分析

为提升针对高马赫数发动机的模拟能力,对计算方法进行了可压缩性修正,并针对飞行Ma12条件下超燃冲压发动机进行了多状态三维数值模拟,分析了发动机内波系、参数以及燃烧性能特征.研究结果表明:(1)修正后的方法计算所得激波位置及强度与试验值吻合,在激波串模拟、高马赫数发动机模拟上均展现了更优的能力.(2)发动机内形成激波与反...     

高压捕获翼构型亚跨超流动特性数值研究

为研究高压捕获翼布局在亚跨超条件下的流动特性, 选取圆锥?圆台机体组合捕获翼概念构型, 在马赫数0.3 ~ 3速域范围内, 选取典型状态点, 采用数值模拟在 0°攻角条件下进行了计算和分析. 结果表明, 在整个速域范围内, 由于机体与捕获翼在对称面附近的垂向距离最小, 因此二者之间的气动干扰最为明显, 且沿展向逐渐减弱. 同时, 随马赫数增大, 机体与捕获翼间的流场结构明显不同, 具体表现为: 当Ma<0.5时, 未出现流动分离现象, 当Ma>0.5时, 机体后段开始出现明显的流动分离, 由于捕获翼与机体形成先收缩后扩张的等效通道, 捕获翼下表面和机体上表面的压力均先减小后增大; 进入跨声速速域后, 在捕获翼的影响下, 流动分离更加明显, 机体与捕获翼之间开始出现激波, 并且与分离区相互作用, 同时出现激波串, 捕获翼下表面产生明显的压力波动现象, Ma=1.5时, 通道内激波位置基本到达机体尾部, 分离区基本消失; 当Ma>2以后, 整个流场呈现以激波为主导的结构形式, 捕获翼下表面和机体上表面的压力分布逐渐趋于平缓.      

Numerical investigation of axisymmetric shock wave focusing over paraboloidal reflectors

The problem of a plane shock wave incident to a paraboloidal reflector is numerically investigated. The numerical solver used is developed by an improved, implicit, upwind total variation diminishing scheme in a finite-volume approach. The real-gas effect is taken into account if high temperature occurs. The solver is validated on four test problems. The complicated flow fields of axisymmetric shock wave focusing for different-depth reflectors at various incident shock Mach numbers are studied. An interesting result of a maximum pressure happening at the reflector center is found. This is due to the occurrence of an implosion phenomenon. A maximum temperature might occur at the reflector center or at other locations, depending on the incident shock Mach number and the reflector depth. Moreover, vortical flows induced by shock wave focusing and their formation mechanism are explored. It was found that the vortices near the reflector are caused by a ring-shaped shock/slipline interaction. Owing to the slipline on the symmetry axis, a jet flow is induced, resulting in the formation of vortices near the symmetry axis. Received 13 January 1998 / Accepted 10 November 1998     

Shock wave interaction with area changes in ducts

Whitham's approximation for handling shock wave propagation in area changes (reductions) in a duct was checked in comparison with a numerical solution. Also the Whitham approximation for shock wave propagation from a constant cross-sectional duct to a duct of a smaller cross-sectional area was studied and compared with a numerical solution. It was found that for modest incident shock Mach numbers and modest area reductions the Whitham approximation provided a fair solution for the shock Mach number and for the post-shock pressure. For higher shock Mach numbers and/or area reductions, large discrepancies exit between the approximate and exact solutions. A wider range of applicability of the Whitham approximation is found for the monotonical area reduction case; it is quite narrow for the passage of a shock wave from a wider to a narrower duct case. In addition, the effect of the extent of the area change region on the time required for reaching a quasi-steady flow was studied. It was shown that the longer the area change segment is, the longer it takes to reach a quasi-steady flow.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gasdynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel containing a gas in an initial state of homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface and also the hub-to-tip radius ratio. The nature of the rarefaction and shock wave propagation is one-dimensional, although the acoustic waves are diffracted due to the radially varying propagation speed. Under conditions of initially homentropic solid-body rotation, a degree of similarity exists between rotating and stationary shock-tube flows. This similarity is exploited to arrive at an approximate analytical solution to the Riemann problem in a rotating shock-tube.     

Experimental comparison between the effect of electrical discharge and the effect of surface heating on Mach 2 rarefied airflow over a flat plate

This study describes the experimental investigation of the effect of a negative DC glow discharge on a Mach 2 rarefied airflow around a flat plate. More precisely, we will show a comparison between two experiments. In the first one we will observe the effect of discharge by Pitot probe measurement. This discharge is created by applying negative DC potential difference between two electrodes flush mounted on the surface of a quartz flat plate placed in Mach 2 rarefied airflow. The electrodes are arranged in the spanwise direction. In the second experiment, electrodes are removed and replaced with a surface heater. The pressure profiles obtained by a glass Pitot tube are presented, and a comparison between the plasma effect and the surface heater effect is made, for the same surface temperature and in thermal equilibrium, with the aim of identifying the origin of the observed effect. For both experiments, surface heating causes a decrease in the boundary layer stagnation pressure, while increasing the boundary layer thickness, with the effects becoming larger for higher mean surface temperature. The effects due to the plasma actuator seem to be larger over the active electrode.