Nonlinear analysis of hydrodynamic instability in laminar flames—I. Derivation of basic equations

An asymptotic nonlinear integrodifferential equation is derived for spontaneous instability of the plane front of a laminar flame.If the combustible mixture is deficient in the light component, spontaneous instability will lead to self-turbulization of the flame, and the flame front assumes a strongly nonstationary cellular structure.If there is an excess of the light component, spontaneous instability produces stationary, irregular wrinkles on the flame front, and the flame continues to propagate in a laminar regime.It is shown that in all cases spontaneous instability of the flame implies an increase in its propagation velocity.     

On a distorted flame front as a hydrodynamic discontinuity

The radius of curvature of a steady distorted flame often turns out to be considerably larger than the width of the thermal structure of the flame. Thus, even under finite deformations of the flame, its structure remains quasi-one-dimensional. This property enables the propagation velocity of a distorted flame front (relative to the gas) to be determined explicitly as a function of the hydrodynamic field and the physico-chemical parameters of the gaseous mixture. With the aid of this relationship one can try to determine the shape of the front in a model which is external with respect to the flame structure and treats the flame as a density jump in an ideal incompressible fluid. As applications, we consider (1) the structure of a Bunsen cone, (2) the extinction of a flame in a divergent flow and (3) the propagation velocity of a corrugated flame.     

An experimental study on the stability of jet diffusion flame

The stability behavior of jet diffusion flame developing in a coflowing high temperature air stream was studied experimentally, using city gas and hydrogen as fuel gases. The primary variables studied were temperature and velocity of the air stream. Two distinct types of stability limits were observed. The first is the blow-off of the rim-stabilized flame, which is familiar and has been studied by a number of investigators. The second is the break-off or extinction of the turbulent portion of the flame at the transition point from laminar to turbulent flow. The second limit is quite unfamiliar and is interesting in relation to the structure of turbulent diffusion flames. The important implications of the experimental findings concerning the stability and the structure of turbulent diffusion flames is discussed.     

同轴喷嘴自发激励高频燃烧不稳定性试验研究

为增进对液氧煤油火箭发动机同轴离心喷嘴燃烧不稳定性过程的理解,在大气环境下进行了同轴离心喷嘴的自发激励燃烧不稳定性试验。试验采用单喷嘴敞口模拟燃烧室,高温的氧气和空气混合物从同轴喷嘴的直流喷嘴喷注,高温的煤油蒸气从同轴喷嘴的离心喷嘴喷注。通过逐步改变氧化剂流量使模拟燃烧室内产生自发激励高频燃烧不稳定性,使用脉动压力传感器和黑白高速相机记录稳定和不稳定燃烧工况下的脉动压力和火焰。研究发现：气气同轴离心喷嘴的自发激励高频燃烧不稳定过程呈现“滞后”现象;不稳定工况下的火焰均为脱口火焰,火焰特征长度约等于喷嘴出口到脱口火焰团上沿的距离;气气同轴离心喷嘴燃烧不稳定性的发生原因可以被认为是因混合特征时间与声学特征时间相关。     

Combustion characteristics of the end burning hybrid rockets in laminar flow

In this study, we aim to clarify the blowoff mechanism for flame spreading in an opposed laminar flow in narrow solid fuel ducts. To clarify this mechanism we conducted two experiments. First, we observed the changes of the flame spread rate at various oxygen velocities, ambient pressures, and port diameters. For flame spreading in laminar flow, combustion modes could be classified into 3 distinct regimes based on the strength of the opposed flow, i.e., chemical regime, thermal regime, and stabilized regime. This result is consistent with the result in turbulent flow. In the stabilized regime, quenching distance is almost constant despite oxygen velocity. In order to investigate the effect of ambient pressure and port diameter of fuels on blowoff limit, transition oxygen velocity is observed. As a result, transition oxygen velocity is proportional to the logarithm of the ambient pressure and port diameter. This relation is applicable despite the flow condition. Furthermore, we calculated velocity gradient at the fuel surface to reveal the determining factor of the blowoff limit in laminar flow. Consequently, velocity gradient, which is considered to dominate flow separation in laminar flow, would not be constant. This is because the velocity gradient at the fuel surface could not be evaluated by only the assumption of Hagen–Poiseuille flow but other parameters, such as vaporized fuel gas and natural convection by buoyancy should be included.     

临近空间高超声速飞行器RCS特性研究

首次系统地对升力体外形临近空间飞行器雷达散射截面（RCS）特性开展了研究,主要研究内容包括本体及绕流RCS特性研究,亚密湍流尾迹RCS特性研究和层流尾迹的RCS特性研究等。结果显示：等离子体绕流将降低飞行器后向RCS;亚密湍流尾迹对低频段RCS影响明显,对高频段影响不明显。另外,双站雷达对临近空间高超声速飞行器也会有较好的探测效果。     

Prediction of turbulent diffusion flames with a four-equation turbulence model

For turbulent flame prediction a four-equation turbulence model consisting of the turbulent kinetic energy, the dissipation rate and the density-velocity correlations using unweighted statistics is suggested. The fluctuating concentrations are taken into account by using Spalding's (1971) equation for the variance of the mixture fraction f and by defining the probability density function of f as beta-function. The thermodynamic model for the combustion of H₂ with air assumes infinitely fast reaction in a global single step. Preliminary results are shown for the calculation of Kent and Bilger (1972) turbulent H₂-air diffusion flame. For this flame with fluctuating density the inclusion of the velocity-density correlations appears to be essential, because they have a strong effect especially on the development of the turbulent viscosity and the mean velocity and mass fraction profiles. Consequently the mean concentrations are in good agreement with the experimental results.     

场协同原理在高超声速化学非平衡流动中的推广

应用理论分析与数值模拟方法,将对流传热的场协同原理从不可压缩流动推广至高超声速化学非平衡流动中。结果表明,高超声速化学非平衡层流与湍流的热流密度取决于流动的当地单位体积的动量与单位质量总焓梯度的协同。用当地单位体积的动量与单位质量总焓梯度的协同研究高超声速化学非平衡流动的壁面传热问题,对层流流动下的对流传热,不但计及了高超声速化学非平衡流的密度变化对热流密度的影响,而且包括了静焓梯度、压力梯度、边界层内的分子黏性剪切效应对热流密度的作用;对湍流问题,除了上述层流流动各项对热流密度的影响外,还计及了雷诺剪切应力对热流密度的作用。考虑到高超声速化学非平衡流静焓的定义,高超声速化学非平衡层流及湍流的场协同同时计及所有组分的平动能、转动能、振动能及电子能等梯度的贡献。     

Stability of amino acids and their oligomerization under high-pressure conditions: implications for prebiotic chemistry

The polymerization of amino acids leading to the formation of peptides and proteins is a significant problem for the origin of life. This problem stems from the instability of amino acids and the difficulty of their oligomerization in aqueous environments, such as seafloor hydrothermal systems. We investigated the stability of amino acids and their oligomerization reactions under high-temperature (180-400°C) and high-pressure (1.0-5.5?GPa) conditions, based on the hypothesis that the polymerization of amino acids occurred in marine sediments during diagenesis and metamorphism, at convergent margins on early Earth. Our results show that the amino acids glycine and alanine are stabilized by high pressure. Oligomers up to pentamers were formed, which has never been reported for alanine in the absence of a catalyst. The yields of peptides at a given temperature and reaction time were higher under higher-pressure conditions. Elemental, infrared, and isotopic analyses of the reaction products indicated that deamination is a key degradation process for amino acids and peptides under high-pressure conditions. A possible NH(3)-rich environment in marine sediments on early Earth may have further stabilized amino acids and peptides by inhibiting their deamination.     

煤油/氧气同轴离心喷嘴燃烧数值模拟

为了更好地了解同轴离心喷嘴的工作特性,基于DDES模型研究了油气比分别为0.5、1、1.5下以煤油/氧气为推进剂的喷嘴的流体动力学特性与非预混燃烧特征.研究结果表明:由于旋流离心作用,在喷嘴出口轴心处和燃烧室顶部分别存在一个驻定涡和角涡,驻定涡径向分布在0.9 R～1.4 R,轴向尺寸在-1 R～14 R,随着燃料流量...     

An experimental study of flame turbulence

An experimental study of the turbulence generated by an enclosed, premixed, propane-air flame has been carried out in a combustion chamber of 25 × 20 cm cross section. Care was taken to reduce any effects of the axial pressure gradient. By suitable changes in the grid geometry, the turbulence intensity and scale of the approach flow were varied independently. The results of these experiments show that a strong link exists between the mechanisms of turbulent flame propagation and flame-generated turbulence. Thus three distinct regions may be identified, each having different characteristics in regard to the effects of turbulence scale on flame-generated turbulence. For each region, a physical mechanism for flame-generated turbulence is proposed. In particular, it is observed that over a wide range of intensity and scale of the approach turbulence, (a) the relative turbulence intensity in the flame zone varies in the range 1–2 times the relative turbulence intensity of the cold flow, (b) in the region of intermediate turbulence levels ( ) the flame-generated turbulence intensity reaches a minimum value which is equal to the approach stream turbulence intensity, (c) the flame-generated turbulence intensity reaches a maximum value when the rate of production of turbulent vorticity is equal to about half its rate of viscous dissipation.     

轴对称后向台阶不稳定流动及压强振荡数值研究

使用翼柱形装药的战术发动机,在翼槽烧完时会出现轴对称的后向台阶结构,易产生转角涡脱落引起的声涡耦合现象。根据战术发动机中后向台阶流动特点设计简化模型,从线性稳定性理论出发,应用声涡耦合模态预测方法进行分析。采用湍流流动的大涡模拟(LES)方法,进行了较宽来流Re数条件下后向台阶不稳定流动的数值模拟,获得了转角涡脱落产生的一般过程,分析了不同来流条件下的旋涡多尺度运动规律。通过分析流场中关键位置的压强振荡,研究了气动声学现象。结果表明,声涡耦合模态预测方法可对不同速度下可能诱发的声场固有频率进行合理预测。流场头部压强波动的振幅最大,频率与声场轴向声模态相耦合。再发展区出现低于轴向基频的宽峰振荡主要受大尺度旋涡运动的影响,该频率与基频接近时,会造成基频的多峰现象。当近壁面小尺寸旋涡形成频率接近声场某阶固有频率时,该频率及更低阶频率就有可能被激发。     

Alfven modes of solar coronal magnetic arches: Excitation of ballooning instability and modulation of flare emission

The effect of Alfven-type oscillations in a coronal magnetic arch on modulation of the gyrosynchrotron radiation and development of the ballooning instability in the arch is considered. On the basis of the energy method and the method of normal modes, the expressions are obtained for increments of ballooning instability at its swinging by natural oscillations of the arch. The conclusion is drawn that bending oscillations, which do not actually compress the plasma and, therefore, represent the Alfven-type modes, unlike the radial oscillations, are capable, under solar corona conditions, to effectively swing ballooning instability and, as a consequence, play a part of a trigger for solar flares. The ballooning instability of coronal arches is shown to be capable of causing formation of helmet-shaped structures in the lower solar corona. On the basis of calculations of the intensity modulation depth and the degree of circular polarization of non-thermal gyrosynchrotron radiation, under the assumption of excited Alfven oscillations of a coronal arch, the conclusion is drawn, that microwave observations at a frequency of > 10 GHz can be used for studying the conditions of excitation and propagation of Alfven modes in flare loops. The consequences of obtained results are discussed using the flare on April 15, 2002 as an example.     

RBCC燃烧室超声速反应混合层特性的大涡模拟

以飞行马赫数为4.5Ma的RBCC发动机典型工作状态为研究背景,采用大涡模拟研究了支板火箭射流和空气来流形成的超声速反应混合层的掺混燃烧过程,获得了燃烧室内详细的流场结构和流动特征,分析了强射流条件下超声速反应混合层的特性。结果表明由于速度梯度的存在,火箭射流进入燃烧室后与空气来流形成环形剪切层,剪切层内丰富的旋涡结构主导火箭射流和空气来流的掺混燃烧,随着湍流能量的串级输运,化学反应过程中释放的能量将被转化成细观尺度的湍流动能,大尺度旋涡将能量传递给小尺度旋涡并最终耗散,细小尺度的旋涡一方面能够促进燃烧反应物的掺混并强化燃烧过程,另一方面会给化学反应过程带来强烈的脉动,使得局部火焰淬灭,火焰结构表现出明显的非定常性。     

Dynamics of debye-scale nonstationary plasma structures in the region of auroral field-aligned currents

We consider the formation of small-scale nonstationary plasma structures in the region of relatively strong field-aligned electric currents. The formation of these structures has been shown to be associated with the density instability developed when the current velocity exceeds a critical value. The conditions for the development of this instability can be most favorable in the region of low-density plasma. Numerical calculations have been performed for the initial nonlinear stage of the structure development. The main parameters of the structure, i.e., the times of its formation and destruction, spatial scales, and electric field, have been estimated. The features of the structures are consistent with the existing data from space experiments in the region of auroral field-aligned currents of the Earth.     

多截面细径管结构参数对放气系统压降特性的影响

采用数值模拟方法研究了多截面细径管结构参数对放气系统压降特性的影响,得到了导管几何尺寸、导管内壁粗糙度和出口压力等参数对压降特性影响的变化规律。研究结果表明:缩短容器导管长度或扩大其内径对壅塞状态下的压降影响较小,但可以显著提高层流状态下气体的流动速度,从而减少放气时间;内壁粗糙度对壅塞状态下的压降特性影响较大,随着粗糙度的增大,压力下降速度有所降低,而对层流状态下的压降特性影响较小,在容器压力降至极低的情况下,粗糙度的影响可以忽略不计;出口压力对整个流动过程的放气速度影响较小,出口压力的进一步提高不会明显增加放气速度。该研究可为多截面细径管放气系统的优化设计提供依据。     

Hydrodynamic instabilities and transverse waves in propagation mechanism of gaseous detonations

The present study examines the role of transverse waves and hydrodynamic instabilities mainly, Richtmyer–Meshkov instability (RMI) and Kelvin–Helmholtz instability (KHI) in detonation structure using two-dimensional high-resolution numerical simulations of Euler equations. To compare the numerical results with those of experiments, Navier–Stokes simulations are also performed by utilizing the effect of diffusion in highly irregular detonations. Results for both moderate and low activation energy mixtures reveal that upon collision of two triple points a pair of forward and backward facing jets is formed. As the jets spread, they undergo Richtmyer–Meshkov instability. The drastic growth of the forward jet found to have profound role in re-acceleration of the detonation wave at the end of a detonation cell cycle. For irregular detonations, the transverse waves found to have substantial role in propagation mechanism of such detonations. In regular detonations, the lead shock ignites all the gases passing through it, hence, the transverse waves and hydrodynamic instabilities do not play crucial role in propagation mechanism of such regular detonations. In comparison with previous numerical simulations present simulation using single-step kinetics shows a distinct keystone-shaped region at the end of the detonation cell.     

A numerical study of the charge diameter effect on unsteady detonation

A numerical study for the unsteady detonation of an unconfined tetryl charge of small diameter, which is assumed to be homogeneous, was performed by using the two-dimensional Lagrangian hydrodynamic computer code, 2 DL, with the first order Arrhenius equation of reaction rate. Becker-Kistiakowsky-Wilson (BKW) and Kihara-Hikita (KH) equations of state have been applied to the detonation products.In the case of BKW, it is shown that the rarefaction waves propagating inward from the lateral surface make the reaction rate slow and give a curvature to the front. Then after an induction time, a strong initiation occurs in the reaction zone near the lateral surface and higher pressure zone moves to the axis. This higher pressure accelerates the detonation propagation near the lateral surface and the curvature of detonation front is reduced. Then, the reaction at the lateral surface again begins to decay by the rarefaction waves. Such a sequence of process is repeated periodically.The possibility of the occurrence of the strong initiation depends on the pressure and temperature in the shocked zone near the surface. In a small diameter charge, the delayed explosion becomes weaker near the surface, while its frequency increases. No shock interaction occurs because the direction of the particle flow is always divergent.In the case of KH equation of state, the temperature of detonation is higher than that obtained by BKW and the behaviour of instability becomes rather different from the previous result, i.e. in the axis the pressure oscillates repeating the overdriven and underdriven detonation similar with the case of BKW.     

高超声速条件下星光传输模拟方法与仿真

以高超声速飞行器为研究对象,建立了高速流场下的星光传输模拟仿真方法,定性、定量地研究平台气动光学效应对星光导航的影响。采用基于k-ω/SST两方程湍流模型的雷诺平均法计算星敏感器安装在飞行器不同位置处（头部、中部和尾部）的外流场密度分布,进而得到流场的空间折射率分布。在此基础上,利用几何光学法计算星光由高速层流流场引起的光程差、点扩散函数、像偏移等光学传输效应,得出经过层流流场的星光降晰图。同时运用统计光学理论计算高速湍流流场引起的星光出瞳面的波像差、密度均方差、相位均方差、Strehl比等光学传输效应,得出经过湍流流场的星光降晰图。仿真结果与理论分析吻合,从而为研究星光在高超声速流场中传输的机理提供了可靠的技术手段。     

Laminar dust flames in a reduced-gravity environment

The propagation of laminar dust flames in suspensions of iron in gaseous oxidizers was studied in a low-gravity environment onboard a parabolic flight aircraft. The reduction of buoyancy-induced convective flows and particle settling permitted the measurement of fundamental combustion parameters, such as the burning velocity and the flame quenching distance over a wide range of particle sizes and in different gaseous mixtures. Experimentally measured flame speeds and quenching distances were found in good agreement with theoretical predictions of a simplified analytical model that assumes particles burning in a diffusive mode. However, the comparison of flame speeds in oxygen–argon and oxygen–helium iron suspensions indicates the possibility that fine micron-sized particles burn in the kinetic mode. Furthermore, when the particle spacing is large compared to the scale of the reaction zone, a theoretical analysis suggests the existence of a new so-called discrete flame propagation regime. Discrete flames are strongly dependent on particle density fluctuations and demonstrate directed percolation behavior near flame propagation limits. The experimental observation of discrete flames in particle suspensions will require low levels of gravity over extended periods available only on orbital platforms.