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

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

极高超声速稀薄气体原子辐射效应的p-DSMC方法

极高超声速流动激波层内的高温导致内能模态的激发并伴随热辐射发生, 过高的温度使得空气分子完全解离, 原子组分对辐射热的贡献将达到80%以上. 本文基于优化的原子辐射模型, 提出追踪光子?直接模拟蒙特卡罗(p-DSMC)方法, 研究了稀薄流区不同马赫数下的高超声速二维圆柱绕流的壁面辐射加热, 获得了有无激发辐射效应的壁面压力和热流以及沿驻点线变化的平动、振动和转动温度. 在不考虑激发辐射效应的情况下, 得到的壁面压力和热流与已有的模拟结果符合的非常好, 误差均在5%以内, 尤其是在驻点位置, 误差在1%以内; 获得的平动、振动以及转动温度均与文献结果符合的很好. 在相同的来流条件下, 考虑辐射效应后发现, 来流速度低于10 km/s时, 辐射加热不明显, 在驻点区域, 辐射加热占对流加热比重在7%左右; 来流速度大于10 km/s时, 在驻点区域, 辐射加热占对流加热比重将超过30%. 考虑辐射效应后, 对非平衡区的平动、转动和振动温度的最大值影响不大. 此外, 另一个重要结论是, 流场中原子的浓度是影响壁面辐射热流大小的一个重要因素.      

Measurements of non-equilibrium and equilibrium temperature behind a strong shock wave in simulated martian atmosphere

Non-equilibrium radiation measurements behind strong shock wave for simulated Martian atmosphere are presented in this paper. The shock wave is established in a hydrogen oxygen combustion driven shock tube. Timeresolved spectra of the Δv = 0 sequence of the B ²Σ⁺ → X ²Σ⁺ electronic transition of CN have been observed through optical emission spectroscopy (OES). A new method, which is based on fitting high resolution spectrum for rotational and vibrational temperatures measurement, is proposed to diagnose temperature distribution behind the shock wave. It is estimated that the current scheme has the maximum deviation less than 8% (1σ) for vibrational temperature measurement through detailed analysis of the influence of the uncertainties of spectroscopic constants and spectral resolution. Radiation structure of the shock layer, including induction, relaxation and equilibrium process, and corresponding rotational and vibrational temperatures are obtained through time gating OES diagnostics with sub-microsecond temporal resolution. The present extensive results will strongly benefit the reaction rate estimation and computational fluid dynamics (CFD) code validation in high enthalpy Mars reentry chemistry.     

Computation and experimental validation of N_{2}-CH_{4}-Ar mixture flows behind normal shock wave

Different vibration-dissociation-vibration coupling models have been used to compute the nonequilibrium N-CH-Ar mixture flow behind a normal shock wave. A three-temperature model was used and the diffuse nature of vibrational relaxation at high temperatures was accounted for. The numerical results obtained with the Treanor and Marrone model (preferential or non preferential) and the Park model of vibration-dissociation-vibration coupling are compared. These results show that the temperatures and the concentrations are mainly affected by the value of the characteristic temperature U in the preferential model of Marrone and Treanor. An assessment of the more realistic model was realized by comparing numerical results with shock tube experiments. The influence of argon addition on the nonequilibrium emission of CN behind the shock wave was also numerically studied and compared to experimental measurements. Received 1 September 1995 / Accepted 10 December 1996     

Experimental and numerical studies of rotational relaxation behind a strong shock wave in air

This paper describes the experimental and numerical investigations of unknown characteristics of the rotational nonequilibrium phenomena behind a strong shock wave in air. Experiments were carried out using a piston-driven shock tube with helium as driving gas and air as driven (test) gas, operated as a two-stage shock tube. In the experiments, emission spectra of NO were measured to evaluate the rotational temperature behind a strong shock wave. The numerical calculations use the computational code for the thermal and chemical nonequilibrium flow behind a strong shock wave developed by the present author's group, where 11 chemical species (N, O, NO, N, O, N, O, NO, N, O, e) and 48 chemical reactions of high-temperature air are considered. The thermal nonequilibrium is expressed by introducing an 8 temperature model composed of translational temperature, rotational and vibrational temperatures for N, O, NO, and electron temperature. The coupling of a rotation, vibration and dissociation (CRVD) model was incorporated to take sufficiently into account the rotational nonequilibrium. The calculations were conducted for the same conditions as the experimental ones. From the calculated flow properties, emission spectra were re-constructed using the code for computing spectra of high temperature air “SPRADIAN”. Furthermore, rotational and vibrational temperatures of NO (0,1) were determined from a curve fitting method and compared with the computed results. Received 12 September 2001 / Accepted 18 February 2002     

低稀释度条件下乙烯点火特性的激波管研究简

在激波管中利用反射激波后高温环境加热燃料,以燃料点火过程中氢氧自由基特征发射光谱强度的急剧变化作为点火发生的标志,在温度范围800 1 650K,压力0.2MPa, 0.7MPa, 1.2MPa,化学当量比为0.5, 1, 2,O₂ 浓度为空气含量20% 的条件下,进行了C₂H₄/O₂/Ar 混合气在低稀释度条件下点火特性的实验研究. 获得了乙烯点火延时随温度、压力、化学当量比、燃料以及氧化剂浓度等参数变化的拟合关系式. 对乙烯点火转爆轰现象进行了初步观察,考察了初始温度对乙烯点火特性以及点火转爆轰的影响.     

STUDIES ON THE AUTOIGNITION CHARACTERISTICS OF RP-3 AVIATION KEROSENE

在预加热到135℃的激波管反射激波后5区,以点火过程中OH自由基在306.5nm处特征发射光谱强度的急剧变化作为点火发生的标志,进行了RP-3航空煤油点火特性的实验研究.实验温度范围为800～1450K,当量比为0.5,1,1.5,压力为0.05,0.1,0.2MPa,O₂的摩尔浓度为空气含量20%.实验获得了低压条件下（0.05,0.1,0.2MPa）RP-3航空煤油点火延时与点火温度﹑压力﹑当量比以及煤油和氧气浓度的依赖关系.将低压实验结果与高压（0.55,1.1,2.2MPa）条件下煤油点火特性进行了对比.结果显示,当量比对煤油点火特性的影响存在一个临界温度.在临界温度以上的高温区,煤油点火延时随当量比增加而增长;在临界温度以下的低温区,煤油点火延时随当量比增加而缩短;这一临界温度随点火压力的降低而升高.采用3种煤油燃烧反应动力学机理对煤油点火过程进行了动力学数值模拟,并与实验结果进行了对比.结果显示,Honnet等提出的煤油反应机理在高压（2.2MPa）下与实验结果吻合得很好,而在低压下有一些差异.对不同压力条件下的点火过程进行敏感度分析表明,三体反应H+O₂+M=HO₂+M在高压时对煤油点火起轻微抑制作用,而在低压时对煤油点火起促进作用.     

Structure of the heating layer ahead of the front of a strong intensely emitting shock wave

The problem of the structure and brightness of strong shock waves arises in the investigation of such phenomena as the motion of large meteoroids in the atmosphere, optical and electrical discharges, the development of strong explosions, and other similar processes and in the creation of powerful radiation sources based on them. This problem also has a general physics interest. As the propagation velocity of a strong shock wave increases the gas temperature behind its front and the role of emission grow. Part of the radiation emitted by the gas heated and compressed in a shock wave is absorbed ahead of the front, forming the so-called heating layer. The quasisteady structure of a strong intensely emitting shock wave was studied in [1, 2]. In this case a diffusional approximation and the assumption of a gray gas were used to describe the radiation transfer. They introduced the concept of a wave of critical amplitude, when the maximum temperature T_- in the heating layer reaches the temperature T_a determined on the basis of the conservation laws, i.e., from the usual shock adiabat; it is shown that behind a compression shock moving through an already heated gas there is a temperature peak in which the maximum temperature T₊ exceeds T_a. The problem of the quasisteady structure of an emitting shock wave in air of normal density was solved numerically in [3]. The angular distribution of the radiation was approximately taken into account — it was assigned by a simple cosinusoidal law. The spectral effects were taken into account in a multigroup approximation. They introduced 38 spectral intervals, which is insufficient to describe a radiation spectrum with allowance for the numerous lines and absorption bands.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 86–92, September–October, 1978.     

STUDY ON THE RELATIONSHIP BETWEEN IGNITION DELAY TIME AND GASEOUS DETONATION CELL SIZE

讨论点火延迟时间和爆轰波胞格尺度的内在关系,将点火延迟时间作为特征参量来模拟胞格尺度. 分别对两个总包单步化学反应模型和一个基元反应模型的点火延迟时间进行了数值模拟研究. 对于满足当量比的氢气/空气混合气体,分析了不同初始压力下点火延迟时间随初始温度的变化关系. 研究表明:总包单步反应模型的点火延迟时间不随压力变化,且与初始温度呈线性关系. 基元反应模型的点火延迟时间随压力变化,而且存在理论上的S 型曲线,但是在拐点区域和低温区域与CHEMKIN 计算的结果相差1~3 个量级. 现有模型模拟的胞格尺度普遍偏小,其相应的点火延迟时间也偏小,胞格尺度与点火延迟时间具有正相关性. 入射激波后的气体的点火延迟时间与三波点的运动周期一致,是定量化模拟胞格的关键因素.     

Experimental investigation of CO vibrational deactivation in a supersonic cooling gas flow

Abstract. The vibrational relaxation time of CO molecules at collisions with H atoms was measured in shock tube experiments by means of the CARS-spectroscopy method. The measurements of the CO vibrational temperature at gas temperatures of 1800–3000 K were performed in a supersonic cooling gas flow. The gas was heated behind the incident and reflected shock wave and then flowed out of a wedge-shaped nozzle. H atoms were generated in the reflected shock wave because of dissociation of H and HO admixtures. The extremely high efficiency of H atoms in CO vibrational deactivation was confirmed. Received 1 February 2000 / Accepted 20 February 2000     

Modeling of vibration-dissociation oxygen kinetics at temperatures of 4,000-11,000 K

The time profiles of vibrational molecular oxygen temperature T _v measured earlier in experiments behind a strong shock wave were used for testing the theoretical and empirical models of thermal nonequilibrium dissociation of molecules. To do this, dissociating gas flows behind the strong shock wave front were calculated with account for these models. If the initial gas temperature behind the wave front T ₀ < 6.5×10³ K, the models well describe changing the temperature with time. However, for T ₀ > 7×10³ K neither of the models tested describes the measured temperature profiles satisfactorily. Using the empirical model proposed in the present study made it possible to satisfactorily describe the vibrational temperature evolution observed in experiments at temperatures up to 11×10³ K.     

ON SHOALING PROPERTY OF A SET OFFOURTH DISPERSIVE BOUSSINESQ MODEL

Boussinesq 类水波模型在港口、海岸以及海洋工程领域应用广泛,但以前对这类模型的变浅性能的研究不够充分. 针对Madsen 和Schäffer 提出的一组四阶Boussinesq 方程,从理论和数值两个方面对这一问题进行了探讨. 理论分析了其变浅性能,指出该文献中参数α₂ 和β₂ 的取值是不合理的,并重新确定其取值. 在交错网格下建立了基于混合4 阶Adams-Bashforth-Moulton 格式的预报-校正数值模型. 数值模拟了两个典型算例： 一是缓变平坡地形上波浪的传播变形,二是波浪在淹没梯形潜堤上的波浪演化过程. 计算结果分别与解析结果、物理模型实验结果进行了比较,发现变浅系数的取值对数值结果影响很大,新参数比原文参数模拟结果的吻合程度更高,这佐证了理论分析.