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.

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

在高超声速飞行条件下，流入冲压发动机燃烧室并降至低速的空气的温度随飞行马赫数增加而愈来愈高。燃料与高温空气混合燃烧释放的化学能中的一部分将转化为解离能。这些解离能在长度受限的尾喷管中难以充分复合形成推力，使冲压发动机推力在高超声速范围内随飞行马赫数增大而下降，难以满足高超声速飞行器的推进要求。 与亚燃冲压发动机相比，流入超燃冲压发动机燃烧室的空气的温度在同样飞行马赫数条件下将明显降低，上述困难可大大缓解。然而目前超燃冲压发动机还存在关键性难点有待克服。若保持现有亚燃冲压发动机的吸气与燃烧方式，通过催化促进燃气解离组分在尾喷管膨胀过程中复合，可以增大冲压发动机的推力，满足高超声速飞行器的推进要求，为高超声速飞行器推进提供新的选择。 本论文主要研究内容如下： （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)在边界层外层,越靠近壁面,作用参数越小;而在边界层近壁区黏性底层,虽然惯性力减小,但黏性力却迅速增加,因此越靠近壁面,作用参数反而越大,加速低速流的代价增加.

激励射流在超声速流场中的混合特性研究

将Hartmann-Sprenger(H-S)管集成到传统射流装置中形成激励射流,研究其在超声速流场中的混合特性.设计了3种频率的C02激励射流,将其横向注入Ma=2.5的均匀超声速来流当中.采用平面激光散射技术对其进行瞬态可视化成像,利用组分采样进行时均分析,结果表明: H-S激励射流可以有效提高射流穿透深度,形成较多大尺度结构,从而获得较好的混合效果

Study on the Mechanism of Detonation to Quasi-detonation Transition

In this paper, the mechanism of detonation to quasi-detonation transition was discussed, a new physical model to simulate quasi-detonation was proposed, and one-dimensional theoretical and numerical simulation was conducted. This study firstly demonstrates that the quasi-detonation is of thermal choking. If the conditions of thermal choking are created by some disturbances, the supersonic flow is then unable to accept additional thermal energy, and the CJ detonation becomes the unstable quasi-detonation precipitately. The kinetic energy loss caused by this transition process is firstly considered in this new physical model. The numerical results are in good agreement with previous experimental observations qualitatively, which demonstrates that the quasi-detonation model is physically correct and the study are fundamentally important for detonation and supersonic combustion research.

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.

Modeling study to compare the flow and heat transfer characteristics of low-power hydrogen, nitrogen and argon arc-heated thrusters

A modelling study is performed to compare the plasma °ow and heat transfer char- acteristics of low-power arc-heated thrusters (arcjets) for three di®erent propellants: hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equa- tions, which take into account the e®ects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, veloc- ity and Mach number distributions calculated within the thruster nozzle obtained with di®erent propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the °ow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appear- ing near the cathode tip; the °ow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant °ows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, speci¯c enthalpies and thermal conductivities, are di®erent, there are appreciable di®er- ences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest speci¯c impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.