基元化学反应一维爆轰波的数值模拟

应用基元化学反应模型对一维氢氧混合气体爆轰波进行数值模拟。采用一端高温高压点火起爆的方式,模拟爆轰波的传播。对一维爆轰波的结构、传播特征以及各个化学组分的变化趋势进行了研究。通过将模拟的结果同国内外的研究结果进行对比,发现了一些规律,为多维爆轰波的数值模拟和脉冲爆轰发动机的研究奠定一定的基础。     

斜爆轰波系在受限空间内的演变及其临界条件的数值研究

作为一种面向未来高马赫数吸气推进技术的新概念发动机,斜爆轰发动机的运行依赖其稳定的宏观波系结构。本文针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析爆轰波系结构转变的临界条件。研究结果发现,随着楔面角度的增加,依次出现四种结构：激波诱导燃烧,斜爆轰波双规则反射,回流区马赫反射,楔面燃烧。对于稳定的波系结构,楔面压缩角同时存在下临界、亚临界以及超临界三种极限条件,在波系演变过程还伴随激波规则反射和马赫反射的转变。     

对数值模拟气相爆轰中二阶段化学反应模型再研究

采用改进的高精度时-空守恒元解元算法(the space-time conservation element and solution element method,CE/SE method)和考虑组分的二阶段化学反应模型(Sichel的二步模型)对气相爆轰问题的数值模拟进行了分析.分析发现采用Sichel的二步模型得到的数值结果虽然比早期二阶段化学反应模型(旧二步模型)更接近实验值,但是仍然不能得到爆轰过程准确气体动力学参数.为此通过修改组分的质量分数分布形式对Sichel的二步模型进行了改造,然后采用新的二步模型对平面爆轰波进行了数值模拟.数值结果表明采用新的二步模型计算得到气体动力学参数更接近于实验值和基元反应模型的计算值,在计算精度上有较大提高.     

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

旋转爆轰发动机工作过程的数值模拟

为了深入研究旋转爆轰发动机工作过程,采用高精度高分辨率迎风通量分裂格式(AUSMPW+)、氢氧7组分8反应化学动力学模型,求解三维带化学反应的Euler方程。对旋转爆轰发动机工作过程进行了数值模拟,分析了采用预爆轰管点火过程的流场特征以及整个发动机推力、推力偏心距、侧向力随时间的变化规律。计算结果表明:切向预爆轰管紧贴喷注入口时,点火过程仅产生一道单向爆轰波,能成功点燃旋转爆轰发动机;在本文给定的计算条件下,旋转爆轰发动机平均推力约180N,旋转爆轰波传播频率约为14285Hz;旋转发动机正常工作过程中,推力偏心矩、侧向力随时间周期性变化。     

准爆轰波传播机理的探索研究

本文分析了准爆轰波的传播过程与物理机制,提出了一个全新的准爆轰波物理模型,开展了理论分析与数值模拟研究.在该物理模型中考虑了超声速流动化学反应放热造成的动量损失,计算结果表明准爆轰波发展过程是热壅塞的.计算与实验结果的定性吻合表明了本文提出的理论模型是正确的,获得的计算结果是合理的,对于准爆轰波物理机理的讨论是有重要意义的.     

高密实颗粒药床燃烧转爆轰的数值模拟

在实验基础上建立起描述多孔固体推进剂中燃烧转爆轰过程的一维反应两相流模型。用Ｍａｃ Ｃｏｒｃｍａｋ有限差分格式和自适应网格技术数值求解了控制方程。数值结果表明，在燃烧转爆轰的各个阶段，颗粒床都存在不同程度的动态压缩，压缩波和燃烧波的相互作用是燃烧向爆轰转变的内在原因。数值预测与实验结果进行了比较，两者有相当的一致性。     

斜爆轰发动机燃烧机理试验研究

为了研究斜爆轰发动机的稳定燃烧机理,本文开展了飞行马赫数9的斜爆轰发动机的数值模拟研究和试验研究。首先,设计了全尺度斜爆轰发动机模型,发动机的总长度为2.8m。采用两级进气道压缩,每级压缩角度均为15°。利用三个小支板在进气道前缘主流核心区中进行氢气的喷射和混合。其次,对氢气混合过程和发动机燃烧过程进行了数值模拟研究。控制方程采用带化学反应的雷诺平均的N-S方程,湍流模型采用SST k-ω模型,化学反应模型采用9组份19步反应的基元反应模型。数值结果表明,氢气在进气道内混合得比较均匀,在燃烧室内获得了稳定的斜爆轰流场和正爆轰流场。接下来,在激波风洞中开展了马赫数9状态下的斜爆轰发动机稳定燃烧机理试验研究,在50 ms的风洞有效试验时间内获得了持续稳定的斜爆轰流场,试验结果与数值模拟结果吻合较好,表明在试验中形成了斜爆轰波。研究结果证明了斜爆轰发动机的技术可行性。     

气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.      

圆柱坐标系下连续旋转爆轰发动机的数值模拟

使用两步化学反应模型对连续旋转爆轰发动机(Continuous rotating detonation engine,CRDE)进行了二维数值模拟研究。数值计算获得了同轴圆管腔中间层曲面上连续爆轰的多个循环过程。分析了燃料入射、提前燃烧、爆轰波结构和波传播速度等几个关键问题。计算获得燃烧室内流场结构与之前实验研究结果定性符合。计算中燃料以最大1400m/s沿轴向入射,爆轰波可约以2400m/s沿周向连续循环传播。计算的燃烧室内爆轰波循环频率约为2300Hz。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

Research on Auto-detection for Remainder Particles of Aerospace Relay Based on Wavelet Analysis

Aerospace relay is one kind of electronic components which is used widely in national defense system and aerospace system. The existence of remainder particles induces the reliability declining, which has become a severe problem in the development of aerospace relay. Traditional particle impact noise detection (PIND) method for remainder detection is ineffective for small particles, due to its low precision and involvement of subjective factors. An auto-detection method for PIND output signals is proposed in this paper, which is based on direct wavelet de-noising (DWD), cross-correlation analysis (CCA) and homo-filtering (HF), the method enhances the affectiv-ity of PIND test about the small particles. In the end, some practical PIND output signals are analysed, and the validity of this new method is proved.     

Preparation of Highly Dispersed Antimony-doped Tin Oxide Nano-powder via Ion-exchange Hydrolysis of SnCl4 and SbCl3 and Azeotropic Drying

Antimony-doped tin hydroxide colloid precipitates have been synthesized by hydrolysis of SnCl4 and SbCl3 using: (1) an ion-exchange hydrolysis to remove chlorine ions, and (2) isoamyl acetate as an azeotropic solvent to obviate water. The obtained dried powder is of high dispersivity without any need for further grinding. The size and dispersivity of the final particles are investigated with the aid of TG-DTA, BET, XRD and TEM. After having calcined, the antimony-doped tin oxide nanopowder possesses a tetragonal rutile structure with high dispersivity, uniform particles and low hard agglomeration.     

Effects of Mo on the Microstructure and Hydrogen Sorption Properties of Ti-Mo Getters

The effects of Mo on the microstructure evolution, porosity and hydrogen sorption properties of Ti-Mo getters are investigated in this work. The results show that the addition of Mo prolongs the densification process of Ti-Mo getters and results in a significant amount of sintered pores. With the Mo content increasing, the porosity of getters firstly increases reaching the maximum value as it at- tains about 7.5wt.%, and then drops. At the room temperature, the hydrogen sorption property of getters increases progressively with the Mo content increasing, but the tendency is not very clear before its content lies below 2.5wt.%. When the Mo content achieves about 7.5wt.%, the hydrogen sorption property proves to be the best. The discussion is made about the above mentioned phenomena inclusive of hydrogen sorption properties of getters under different activation conditions (from 500-750 ℃).     

超声速燃烧数值模拟中的湍流与化学反应相互作用模型

高精度数值模拟有助于理解超声速湍流燃烧中湍流与化学反应的相互作用,可为发动机燃烧室等工程应用设计提供可靠的预测模型。除直接数值模拟外,目前在湍流燃烧应用中使用的大涡模拟和雷诺平均Navier-Stokes模拟均需要借助模型模化发生在湍流小尺度上的流动与化学反应过程对湍流大尺度运动的影响。现有的湍流与化学反应相互作用模型大致可分为:火焰面类模型和概率密度函数类模型,2类模型在不同的应用中各自具有优势和局限性。此外,现有模型大都基于低马赫数燃烧,而超声速燃烧中通常会伴随快速混合、局部熄火和再着火以及激波等复杂过程,这为发展其中的湍流与化学反应相互作用模型提出了更多的挑战。     

适于低轨卫星IP网络的单核共享树组播算法(英文)

为了解决低轨卫星IP网络中现有典型源组播算法的信道资源浪费问题,本文提出了一套单核共享树组播算法,即核心群合并共享树(CCST)和加权核心群合并共享树(w-CCST)算法。CCST 算法包括动态近似中心(DAC)选核方法和核心群合并组播路径构建方法。DAC方法专为周期、规律运动的低轨卫星网络提出,不需要复杂的星上计算。在核心群合并方法中,以核节点作为初始核心群,通过核心群和剩余组成员的最短路径方法逐步扩展直至整棵组播树构建完成,从而使得组播树的树代价最小,大大提高了网络的带宽利用率和组播传输效率。w-CCST 算法中所提出的加权因子可以调整树代价和端到端传播时延之间的折衷程度,因此,可以通过调整加权因子来适度增大树代价、降低端到端传播时延以支持某些端到端时延要求苛刻的实时组播业务。最后,与低轨卫星 IP 网络中典型算法进行了性能比较,仿真结果说明,CCST 算法的平均树代价比其它算法显著降低,w-CCST 算法的平均端到端传播时延小于 CCST 算法。     

Jet Vectoring Control Using a Novel Synthetic Jet Actuator

A primary air jet vectoring control system with a novel synthetic jet actuator (SJA) is presented and simulated numerically. The results show that, in comparison with an existing traditional synthetic jet actuator, which is able to perform the duty of either “push” or “pull”, one novel synthetic jet actuator can fulfill both “push” and “pull” functions to vector the primary jet by shifting a slide block inside it. Therefore, because the new actuator possesses greater efficiency, it has potentiality to replace the existing one in various appli- cations, such as thrust vectoring and the reduction of thermal signature. Moreover, as the novel actuator can fulfill those functions that the existing one can not, it may well be expected to popularize it into more flow control systems.     

Nonlinear Adaptive Slewing Motion Control of Spacecraft Truss Driven by Synchronous V-gimbaled CMG Precession

The slewing motion control of a truss arm driven by a V-gimbaled control-moment-gyro (CMG) is a nonlinear control problem. The V-gimbaled CMG consists of a pair of gyros that must precess synchronously. The moment of inertia of the system, the angular momentum of the gyros and the external disturbances are not exactly known. With the help of feedback linearization and recursive Lyapunov design method, an adaptive nonlinear controller is designed to deal with the unknown items. Performance of the proposed controller is verified by simulation.     

An SOR Implicit Time-accurate Scheme in Calculating Unsteady Flows

A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation (SOR) iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computa- tional tests are carried out by using the proposed SOR algorithm. Roe’s FDS scheme is used to discritize the inviscid flux terms. Un- steady computational results of SOR are compared with the experiment results and those of Gauss-Seidel. Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.