超声速气流中液体横向射流雾化过程数值模拟

为了更加深入了解超燃冲压发动机燃烧室中的燃料雾化机理,对来流Mach数为1.94的超声速气流中液体横向射流的雾化过程进行了数值模拟研究.计算采用Euler-Lagrange方法,液滴二次破碎模型采用K-H/R-T模型.计算结果表明:考虑液滴二次破碎时,采用雾化锥模型获得的射流穿透深度以及液滴速度分布与实验结果符合得很好...     

超燃冲压发动机等截面隔离段中流动的LBM模拟

采用耦合的双分布函数格子Boltzmann方法模拟超燃冲压发动机中隔离段的流动.计算结果反映了隔离段内流场的基本结构,表明LBM在隔离段的模拟中是-种具有潜力的数值计算方法.     

进口条件对超燃冲压发动机氢气燃烧流场特性的影响分析

对不同进口条件下的超燃冲压发动机燃烧室内氢气喷流超声速燃烧流动特性进行了数值模拟与分析.宽范围超燃冲压发动机是吸气式高超声速飞行器推进系统设计中的热点问题之一,受实验设备硬件条件及实验技术限制,数值模拟技术仍然是超燃冲压发动机燃烧室内燃气燃烧特性及流场特性的主要研究手段.采用基于混合网格技术的多组元N-S方程有限体积方法求解器,在不同进口Mach数及压强条件下,对带楔板/凹腔结构的燃烧室模型氢气喷流燃烧流场进行了数值模拟,对比分析了氢气喷流穿透深度、喷口前后回流区结构、掺混效率及燃烧效率等流场结构与典型流场参数的变化特性及影响规律.研究成果可为宽范围超燃冲压发动机喷流燃烧流动特性分析提供参考.      

超燃冲压发动机先进性量化评估方法

超燃冲压发动机的评估评价对高超声速飞行技术的发展有重要影响。采用层次分析法并结合超燃冲压发动机的技术特点,利用各性能指标的无量纲化参数,研究建立了量化评估超燃冲压发动机先进性的方法与模型,并通过实际算例对其可行性和合理性进行了说明。该方法将复杂的超燃冲压发动机先进性量化评估问题简化为定量分析和数值模拟的过程,对超燃冲压发动机的评估方法和发展方向具有重要意义。     

JF12激波风洞高Mach数超燃冲压发动机实验研究

针对高Mach数（Ma ≥ 7）超燃冲压发动机高气动阻力下的燃烧组织问题,提出一种双突扩燃烧室结构方案.使用数值模拟方法考察了射流与双突扩燃烧室组合方式的混合燃烧特性.设计了双突扩超燃冲压发动机模型,在力学研究所JF12长试验时间激波风洞内,开展了Ma=7.0和Ma=9.5的氢燃料点火和燃烧试验对比.在风洞有效试验时间100 ms内,实现了Ma=7.0和Ma=9.5超燃冲压发动机的成功点火与稳定燃烧.在Ma=7.0情况下,进气道采用三维压缩,燃烧室入口设计Mach数Ma_c=2.5,壁面压力分布实验结果显示燃烧放热靠近燃烧室扩张段上游;在Ma=9.5情况下,进气道采用二维压缩,燃烧室入口设计Mach数Ma_c=3.5,由于燃烧室流动速度特别高,燃烧放热靠近燃烧室扩张段下游.      

一种新型的超燃冲压发动机闭式冷却循环

热防护技术是发展高超音速的关键技术之一.为了降低冷却用燃料需求,本文提出了基于闭式Brayton循环的超燃冲压发动机冷却循环,旨在提出一个缓解超燃冲压发动机热防护压力的新方法、新思路,同时可为高超音速飞行器提供足够的电力供应.文中介绍了冷却循环的组成和工作原理,定义了评价冷却循环性能的参数,在考虑系统内、外不可逆性的基础上,完成了系统的性能分析,分析表明该冷却循环可有效地节约冷却用燃料消耗,大幅提升超燃冲压发动机热防护系统性能.     

新建高焓激波风洞Ma=8飞行模拟条件的实现与超燃实验

针对高Mach数超燃冲压发动机实验能力空缺问题,基于航天十一院新建的FD-21高能脉冲风洞,进行了Ma=8超燃飞行条件的模拟能力设计与调试,获得了总焓2.9 MJ/kg、总压11.01 MPa实验条件,实现了Ma=8、高度31 km飞行条件的风洞模拟.在此基础上,研发了匹配的氢燃料供应及喷注时序控制系统,设计了超燃冲压发动机模型,开展了超燃冲压发动机模型自由射流应用性风洞实验,获得了氢气燃料与空气、氮气超声速气流耦合流动作用下的实验模型壁面压力数据.在当量比近似一致条件下,空气来流对应的燃烧室壁面压力明显高于氮气来流情况,表明氢气在1 ms有效实验时间内完成了与超声速空气来流的混合、点火与燃烧,获得燃烧释热特性,确认了在FD-21高能脉冲风洞开展高Mach数超燃实验是切实可行的,为后续研究奠定了良好的基础.      

斜爆轰发动机流动机理分析

为了研究高Mach数超燃冲压发动机和斜爆轰发动机的内流场燃烧流动机理，首先用CJ爆轰理论对超燃冲压发动机的内流场特性进行了理论分析，给出了燃烧室流场的气动规律，理论分析结果与现有实验结果吻合得非常好.其次，根据理论分析结果，提出了高Mach数超燃冲压发动机和斜爆轰发动机的气动设计原则.最后，根据提出的气动设计原则，设计了高Mach数斜爆轰发动机，飞行Mach数为9，对斜激波诱导燃烧机理开展了二维数值模拟研究.数值模拟结果表明，在高Mach数下，斜爆轰发动机燃烧室内可以得到稳定的燃烧流场.      

吸气式高超声速飞行器机体/推进一体化内外流非定常耦合方法

基于吸气式高超声速飞行器机体/推进一体化的气动布局设计方式,文章提出了一种内外流一体化流场的耦合求解方法,其中燃烧室内流场采用考虑有限速率化学反应动力学模型的一维非稳态方法求解,进气道和尾喷管外流场采用二维CFD软件计算,进气道与燃烧室在耦合界面处通过一维平均方法实现静温、静压和Mach数等参数传递.并分别以日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为内流场验证算例,以某典型高超声速飞行器一体化模型作为内外耦合流场验证算例.研究结果表明:有限速率化学反应准一维方法能较为准确地模拟燃烧室内燃烧流场,提出的内外流场耦合方法能够有效地计算出内外流耦合效应,计算后体压力分布与理论值较接近.该方法可为超燃冲压发动机的性能快速分析和吸气式高超声速飞行器机体/推进一体化的初步分析设计提供重要参考.      

换热器预冷的空气涡轮火箭性能分析研究

为了准确地分析以不同燃料为冷质的换热器预冷空气涡轮火箭发动机的性能,本文建立了可考虑工质组分化学平衡与工质热物性随温度和压力变化的热力循环模型。利用该模型分析了以LH2、LCH4和煤油为燃料的发动机循环性能,研究了预热温度、燃气发生器和主燃室出口温度对循环性能的影响。结果表明,换热器预冷能够降低压气机进口温度,减小相同增压比下的压缩功,提高燃料温度,减小推进剂流量,从而减小燃空比,增加比冲;燃料预热温度的提高可减小燃气发生器氧燃比,进而减小燃空比、增大比冲;较低的燃气发生器和主燃烧室出口温度可降低燃空比,提高比冲,但会造成单位推力的下降。同时,研究结果还表明,与LCH4和煤油相比,LH2的比热容与热值的比值最大,采用LH2预冷的发动机的飞行速度和比冲的增幅也最大,说明比热容与热值的比值越大,相应的预冷效果越好,即可依据燃料比热容与热值的比值定性衡量发动机的预冷效果。     

Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Dynamic behavior of diffusion flame interacting with a large-scale vortex by laser imaging techniques

Instantaneous and simultaneous measurements of two-dimensional temperature and OH-LIF profiles by combining Rayleigh scattering with laser induced fluorescence (LIF) were demonstrated in a nitrogen-diluted hydrogen (H₂ 30% + N₂ 70%) laminar normal diffusion flame interacting with a large-scale vortex by oscillating central fuel flow or in an inverse diffusion flame by oscillating central airflow. The dynamic behavior of the diffusion flame extinction and reignition during the flame–vortex interaction processes was investigated. The results obtained are described as follows. (1) The width of the reaction zone decreases remarkably, and a decrease in flame temperature and OH-LIF is seen with increasing central airflow in an inverse diffusion flame. OH-LIF increases, and temperature does not change with increasing central fuel flow in a normal diffusion flame. The computations predict the experimental results well, and it is revealed that flame temperature characteristics result from the preferential diffusion of heat and species, which induces excess enthalpy or on enthalpy deficit, and an increase or decrease in H₂ mole fraction in the flame. (2) When a large velocity fluctuation is given to the central flow, the temperature and the OH-LIF at the reaction zone become thin at the convex and circumferential part of the vortex where a high temperature layer exists, and the temperature at the reaction zone is lowered in the inverse flame and the normal flame. (3) The width and temperature of the reaction zone interacting with the vortex recover quickly to that of the laminar steady flame after the vortex passing in the normal flame, but the recovery to that of the steady flame after the vortex passing is delayed in the inverse flame. (4) When a remarkably large velocity fluctuation is given to the central airflow in the inverse flame, thinning of temperature and reaction zone starts at the convex and circumferential part of the vortex, resulting in a and flame extinction completely occurs at the tail part of the vortex and makes the pair of edge flames. The outside edge flame reignites and connects with the upstream reaction zone. The inside edge flame finally extinguishes as the supply of fuel is interrupted by the outside edge flame.     

双组分燃料的层流火焰速度模型

以往关于层流火焰速度的理论分析均只考虑单组分燃料,本文对双组分燃料的平面火焰进行了大活化能渐近理论分析。在理论分析中,将火焰结构分为预热区、化学反应区和平衡区,并在大活化能假设下对各个区域分别求解了关于温度与燃料质量分数的微分方程。根据每两个区域分界面上满足的结合条件,本文推导出了双组分燃料的层流火焰速度模型。该模型表明双组分燃料层流火焰速度的平方为各个单组分燃料层流火焰速度平方的加权平均。     

Behaviour of a confined fire located in an unventilated zone

The behaviour of a fire in an enclosure is studied for a configuration where the fuel source is located in the upper hot unventilated zone trapped by a soffit. The experimental study, undertaken in a laboratory-scale compartment with a fuel source above the level of a soffit, included the determination of the parameters (ventilation factor, rate of fuel supply) controlling the combustion or leading to extinction. Measurements (PIV, thermocouples, gas sampling and analysis) were performed to propose a hypothesis on the structure of the flame (flame stabilisation mechanisms, premixed or diffusion types). Video photography is used to determine the area covered by the flames. This information is used as a criterion to identify the combustion regimes. The results show that the gaseous fuel is diluted in the combustion products (CP) in the upper layer and that a recirculatory motion is formed, driven by buoyancy forces, which enhances the mixing of fuel and CP. These then travel horizontally towards the vent along the interface between the lower fresh air and upper zones, and are premixed with the convected air in the enclosure, before entering the reaction zone and being burnt. The flame stabilises at the interface between the upper hot and lower ventilated layers in the compartment. The observed “ghosting flame” is stabilised by a triple flame if the flame speed of the premixed flame is higher than the natural convection velocity induced in the compartment. The flame stability is quantified by a criterion based on the area of the horizontal flame. It has been observed that the combustion is controlled by the available mass fuel flux at the reaction zone if the ventilation is sufficient. This information is essential for the modelling of the phenomena involved in fires with such an underventilated fuel source.     

The extended IEM mixing model in the framework of the composition PDF approach: applications to diesel spray combustion

A steady, two-dimensional corner flame is established when fuel and oxidizer enter the reaction zone in mutually perpendicular directions. A model problem in which the velocity fields are linear functions of spatial position is utilized to study the resulting flame. The flame structure is comprised of a diffusion flame surrounded on either side by fuel-rich and fuel-lean partially premixed laminar flames, similar to, but distinct from, triple flames. Using suitable coordinate transformations and change of variables, the governing equations in the thermodiffusive approximation are recast into a form akin to classical triple flames, with the strain rate appearing as the eigenvalue. A new exact integral representation of the solution to the mixture fraction equation is then utilized and high activation energy asymptotics are applied to solve approximately for the resulting flame shape, the imposed strain rate and, most significantly, the position of flame stabilization. This theoretically predicted flame is computed numerically, and comparisons are made between theory and computation.     

Interaction of a premixed flame with a liquid fuel film on a wall

In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.     

Opposed flow oxy-flame synthesis of carbon and oxide nanostructures on molybdenum probes

The formation of carbon and metal-oxide nanostructures on molybdenum probes inserted in a counter-flow oxy-fuel flame is studied experimentally. Flame position and probe diameter were varied to achieve a controlled growth of carbon and metal-oxide nanostructures at fuel and oxygen-rich flame zones. Mo probes of 1-mm diameter were introduced in the flame at various heights, starting from the upper hydrocarbon-rich zone on the fuel side of the flame to the oxygen-rich zone on the oxidizer side. High density layers of carbon nanocoils (CNCs) and filamentous structures containing ribbon shapes and straight nanofibers were formed in the upper hydrocarbon-rich flame zone. The formation of carbon micro-trees was observed on the fuel side closer to the flame front. The structures formed in the oxidizer part of the flame were composed of molybdenum-oxides. MoO₂ micron-sized channel structures were formed on the oxidizer side in the vicinity of the flame front. The micro-channels had rectangular and square-framed shapes; they were completely hollow, closed, and semi-open with a small circular cavity at their tips. The application of probes with diameters of 0.75 and 0.25 mm resulted in the formation of spectacular 3-D structures with unique and distinct morphologies.     

Characteristics and structure of inverse flames of natural gas

Characteristics and structure of nominally non-premixed flames of natural gas are investigated using a burner that employs simultaneously two distinct features: fuel and oxidiser direct injection, and inverse fuel and oxidiser delivery. At low exit velocities, the result is an inverse diffusion flame that has been noted in the past for its low NO_x emissions, soot luminosity, and narrow stability limits. The present study aimed at extending the burner operating range, and it demonstrated that the inverse flame exhibits a varying degree of partial premixing dependent on the discharge nozzle conditions and the ratio of inner air jet and outer fuel jet velocities. These two variables affect the flame length, temperature distributions, and stability limits. Temperature measurements and Schlieren visualisation show areas of enhanced turbulent mixing in the shear region and the presence of a well-mixed reaction zone on the flame centreline. This reaction zone is enveloped by an outer diffusion flame, yielding a unique double-flame structure. As the fuel–air equivalence ratio is decreasing with an increase in the inner jet velocity, the well-mixed reaction zone extends considerably. These findings suggest a method for establishing a flame of uniform high temperature by optimising the coaxial nozzle geometry and flow conditions. The normalised flame length is decreasing exponentially with the air/fuel velocity ratio. Measurements demonstrate that the inverse flame stability limits change qualitatively with varying degree of partial premixing. At the low premixing level, the flame blow-out is a function of the inner and outer jet velocities and the nozzle conditions. The flame blow-out at high degree of partial premixing occurs abruptly at a single value of the inner air jet velocity, regardless of the fuel jet velocity and almost independent of the discharge nozzle conditions.     

Propagation and extinction of an unsteady spherical spray flame front

Experimental evidence seems to indicate that the life of a laminar spherical flame front propagating through a fresh mixture of air and liquid fuel droplets can be roughly split into three stages: (1) ignition, (2) radial propagation with a smooth flame front and (3) propagation with flame front cellularization and/or pulsation. In this work, the second stage is analysed using the slowly varying flame approach, for a fuel rich flame. The droplets are presumed to vaporize in a sharp front ahead of the reaction front. Evolution equations for the flame and evaporation fronts are derived. For the former the combined effect of heat loss due to droplet vaporization and radiation plays a dominant explicit role. In addition, the structure of the evaporation front is deduced using asymptotics based on a large parameter associated with spray vaporization. Numerical calculations based on the analysis point to the way in which the spray modifies conditions for flame front extinction. Within the framework of the present simplified model the main relevant parameters turn out to be the initial liquid fuel load in the fresh mixture and/or the latent heat of vaporization of the fuel.     

Laminar flame propagation on a horizontal fuel surface: Verification of classical Emmons solution

This work analyses the classical Emmons (1956) solution of flat plate laminar flame combustion on a film of liquid fuel. A two-dimensional (2D) numerical model developed for this purpose has been benchmarked with experimental results available in the literature for methanol. In the parametric study, numerical predictions have been compared with Emmons classical solution. The study shows that the Emmons solution is valid in a range of Reynolds numbers where flame anchors near the leading edge of the methanol pool and the combustion zone is confined around the hydrodynamic and thermal boundary layers. However, in cases of low free stream velocities the combustion zone is beyond the boundary layer zone and the Emmons solution deviates. In cases of very high free stream velocities, the flame moves away from the leading edge and anchors at a location downstream. The Emmons solution is not applicable in this case as well. For the fuel considered in this study (methanol), accounting for thermal radiation, employing an optically thin radiation model, allows better agreement between experimental and numerical temperature profiles but does not affect the mass burning rates.