带主动冷却的超声速燃烧室传热分析

介绍了流体、固体传热耦合的一维分析方法,对带主动冷却系统的马赫数2.5超声速燃烧室进行了传热分析.该分析以实验测量的燃烧室壁面静压以及超临界煤油换热特性数据为基础,考虑了燃气的高温离解效应,燃烧特性以及碳氢燃料的高温热物理特性,对不同燃烧状态、冷却条件下的主动冷却过程进行了分析.结果表明有燃烧时壁面热流可高达1 MW/m~2以上,是无燃烧时的2～3倍.当煤油流量较小时(当量比为0.45),冷却后的壁面温度仍偏高,而且冷却壁内温度分布不均匀.随着煤油流量的增加,冷却效果明显提高,冷却壁内温度分布趋于均匀;并且煤油的出口温度也显著减小.

一种含运动固壁超声速流动的Descartes网格算法

提出一种Descartes网格算法,用于数值求解含任意复杂及运动固壁的超声速流动问题.采用位标集函数确定和跟踪流-固界面.引入虚网格技术处理流-固边界条件,并沿法向和切向分别进行计算.该算法简单、稳健,可与高阶有限差分格式并用.选取一组一维/二维静止或运动物体绕流算例,验证其有效性.

Catalytic Cracking and Heat Sink Capacity of Aviation Kerosene Under Supercritical Conditions

Catalytic cracking of China no. 3 aviation kerosene using a zeolite catalyst was investigated under supercritical conditions. A three-stage heating/cracking system was specially designed to be capable of heating 0.8 kg kerosene to a temperature of 1050 K and pressure of 7.0 MPa with maximum mass flow rate of 80 g/s. Sonic nozzles of different diameters were used to calibrate and monitor the mass flow rate of the cracked fuel mixture. With proper experiment arrangements, the mass flow rate per unit throat area of the cracked fuel mixture was found to well correlate with the extent of fuel conversion. The gaseous products obtained from fuel cracking under different conditions were also analyzed using gas chromatography. Composition analysis showed that the average molecular weight of the resulting gaseous products and the fuel mass conversion percentage were a strong function of the fuel temperature and were only slightly affected by the fuel pressure. The fuel conversion was also shown to depend on the fuel residence time in the reactor, as expected. Furthermore, the heat sink levels due to sensible heating and endothermic cracking were determined and compared at varying test conditions. It was found that at a fuel temperature of similar to 1050 K, the total heat sink reached similar to 3.4 MJ/kg, in which chemical heat sink accounted for similar to 1.5 MJ/kg.

Aerodynamic characteristics research on wide-speed range waverider configuration

Waverider generated from a given flow field has a high lift-to-drag ratio because of attached bow shock on leading edge. However, leading edge blunt and off-design condition can make bow shock off leading edge and have unfavorable influence on aerodynamic characteristics. So these two problems have always been concerned as important engineering science issues by aeronautical engineering scientists. In this paper, through respectively using low speed and high speed waverider design principles, a wide-speed rang vehicle is designed, which can level takeoff and accelerate to hypersonic speed for cruise. In addition, sharp leading edge is blunted to alleviated aeroheating. Theoretical study and wind tunnel test show that this vehicle has good aerodynamic performance in wide-speed range of subsonic, transonic, supersonic and hypersonic speeds.

Configuration optimization of hypersonic vehicles under transitional flow conditions

The rarefied gas effects on several configurations are investigated under hypersonic flow conditions using the direct simulation Mont Carlo method. It is found that the Knudsen number, the Mach number, and the angle of attack all play a mixed role in the aerodynamics of a flat plate. The ratio of lift to drag decreases as the Knudsen number increases. Studies on 3D delta wings show that the ratio of lift to drag could be increased by decreasing the wing thickness and/or by increasing the wing span. It is also found that the waveriders could produce larger ratio of lift to drag as compared with the delta wing having the same length, wing span, and cross section area.

冲压发动机喷气催化复合效应研究

在高超声速飞行条件下，流入冲压发动机燃烧室并降至低速的空气的温度随飞行马赫数增加而愈来愈高。燃料与高温空气混合燃烧释放的化学能中的一部分将转化为解离能。这些解离能在长度受限的尾喷管中难以充分复合形成推力，使冲压发动机推力在高超声速范围内随飞行马赫数增大而下降，难以满足高超声速飞行器的推进要求。 与亚燃冲压发动机相比，流入超燃冲压发动机燃烧室的空气的温度在同样飞行马赫数条件下将明显降低，上述困难可大大缓解。然而目前超燃冲压发动机还存在关键性难点有待克服。若保持现有亚燃冲压发动机的吸气与燃烧方式，通过催化促进燃气解离组分在尾喷管膨胀过程中复合，可以增大冲压发动机的推力，满足高超声速飞行器的推进要求，为高超声速飞行器推进提供新的选择。 本论文主要研究内容如下： （1） 研究了亚燃冲压发动机燃烧室内燃气解离能与飞行马赫数的关系。通过对冻结流、平衡流和有限化学反应速率的流动的数值计算，确定了回收解离能增大推力的潜力。 （2） 以双爆轰技术为基础，建立起一套地面燃气产生装置。所产生的燃气的组分、温度和压力均与冲压发动机在高空飞行时燃气完全相同。调试出总温3200K、总压20Bar（对应来流马赫数6）和试验时间17.5ms以及总温4000K、总压5Bar（对应来流马赫数8）和试验时间12.5ms两种状态参数的试验用燃气。 （3） 建立了基于动量守恒原理的通过皮托管测压力换算推力的测量方法。对催化复合增大推力的实验而言，一般要进行特定流动条件下喷水与未喷水两种情况下推力大小的比较，其精度可以达到2%甚至更高。 （4） 完成了尾喷管喉道下游管壁喷水试验，成功释放出高温燃气中的解离能，有效增大了推力，证实了催化增推的想法是可行的。在来流马赫数6的条件下获得了11.0%的推力增量；在来流马赫数8.0的条件下也获得了11.7%的推力增加。

Studies of ion energy and yield from argon clusters heated by intense femtosecond laser pulses

Using time-of-flight spectrometry, the interaction of intense femtosecond laser pulses with argon clusters has been studied by measuring the energy and yield of emitted ions. With two different supersonic nozzles, the dependence of average ion energy (E) over bar on cluster size (n) over bar in a large range of (n) over bar approximate to 3 x 10(3) similar to 3 x 10(6) has been measured. The experimental results indicate that when the cluster size (n) over bar <= 3 x 10(5), the average ion energy (E) over bar proportional to (n) over bar (0.5), Coulomb explosion is the dominant expansion mechanism. Beyond this size, the average ion energy gets saturated gradually, the clusters exhibit a mixed Coulomb-hydrodynamic expansion behavior. We also find that with the increasing gas backing pressure, there is a maximum ion yield, the ion yield decreases as the gas backing pressure is further increased.

飞秒强激光场中大尺寸氩团簇爆炸机理的实验研究

By using time-of-flight spectroscopy, the ionization and explosion of large argon clusters ( (n) over bar = 3 x 10(3) - 3 x 10(6)) in the intense femtosecond pulsed laser field (60 fs,2 x 10(16) W/cm(2)) has been studied, and the dependence of average energy of ions emitted from argon clusters on the gas backing pressure has been measured. By comparing the average ion energies obtained with two different supersonic conical nozzles and considering the Hagena's scaling law of clusters, we have found that the average ion energy is determined by the cluster size when the laser parameters are kept unchanged. The experimental results indicated that when the cluster size is less than 3 x 10(5) atoms per cluster, the Coulomb repulsion force is the dominating factor in the expansion mechanism. Beyond this size, for 3 x 10(5) < (n) over barn < 3 x 10(6), the expansion is the result of the combined effect of both the Coulomb repulsion force and the hydrodynamic force, and the latter will play the dominating role for increasing cluster size.

Nuclear fusion from Coulomb explosions of deuterated methane clusters subjected to ultraintense femtosecond laser pulses

This paper reports that Coulomb explosions taken place in the experiment of heteronuclear deuterated methane clusters ((CD4)(n)) in a gas jet subjected to intense femtosecond laser pulses (170 mJ, 70 fs) have led to table-top laser driven DD nuclear fusion. The clusters produced in supersonic expansion had an average energies of deuterons produced in the laser-cluster interaction were 60 and 1.5 KeV, respectively. From DD collisons of energetic deuterons, a yield of 2.5(+/-0.4)x10(4) fusion neutrons of 2.45 MeV per shot was realized, giving rise to a neutron production efficiency of about 1.5 x 10(5) per joule of incident laser pulse energy. Theoretical calculations were performed and a fairly good agreement of the calculated neutron yield with that obtained from the present experiment was found.

Au@Si-n: Growth behavior, stability and electronic structure

The configurations, stability, and electronic structure of AuSin (n = 1-16) clusters have been investigated within the framework of the density functional theory at the B3PW91/LanL2DZ and PW91/DNP levels. The results show that the Au atom begins to occupy the interior site for cages as small as Si-11 and for Si-12 the Au atom completely falls into the interior site forming Au@Si-12 cage. A relatively large embedding energy and small HOMO-LUMO gap are also found for this Au@Si-12 structure indicating enhanced chemical activity and good electronic transfer properties. All these make Au@Si-12 attractive for cluster-assembled materials.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

Nonlinear elastic waves in a monatomic chain with nonlinear interaction

Nonlinear wave equation for a one-dimensional anharmonic crystal lattice in terms of its microscopic parameters is obtained by means of a continuum approximation. Using a small time scale transformation, the nonlinear wave equation is reduced to a combined KdV equation and its single soliton solution yields the supersonic kink form of nonlinear elastic waves for the system.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

基于混合发展的超声速段射流式COIL的数值优化

应用CFD技术,发展三维多组分化学反应流计算程序,对采用超声速段射流方式的氧碘化学激光进行数值模拟,考察分解率与增益等参考变量的详细三维分布.计算结果说明,超声速区域的高速流动以及混合效率降低使COIL无法在有限空间内完成整个运转流程;提高碘含量以加快反应速度的手段导致主流无法提供足够的载能介质,无法形成合理增益.在不改变喷管长度的前提下,提出主流无载气方式的探索,结果证明了超声速段射流方式采用无载气主流配置的可行性与优势,通过减小气流速度保证混合与化学反应在光腔上游完成,在合理的流量配比下光腔位置处可得到1.3%cm-1的增益峰值

MHD控制超声速边界层的理论研究和数值分析

对MHD(mechanisms of magnetohy drodynamics)控制超声速平板湍流边界层的机理进行了理论研究和数值模拟.理论上,采用等离子体低频近似碰撞频率模型,建立等离子体中电子和离子的力平衡方程,得到等离子体速度、极化电场以及边界层速度.数值上,通过空间HLLE格式、LU-SGS时间推进求解时均磁流体动力学湍流方程,其中湍流模型采用sst-kω双方程模型.研究结果表明:(1)边界层速度的理论结果和数值结果误差在7%范围内;(2)只有磁场而电场为零时,洛仑兹力起到减小摩阻的作用.施加电场后,洛仑兹力能够加速边界层低速区流体;(3)在边界层外层,越靠近壁面,作用参数越小;而在边界层近壁区黏性底层,虽然惯性力减小,但黏性力却迅速增加,因此越靠近壁面,作用参数反而越大,加速低速流的代价增加.