非等温液滴蒸发燃烧的传热传质特性

对无重力、相对静止热环境中碳氢燃料单液滴的蒸发燃烧过程,建立了非等温液滴蒸发-火焰面燃烧耦合模型。考虑了液滴内导热、液滴表面蒸发、火焰面化学反应以及气相区对流传热传质的瞬态耦合作用。将整个耦合系统分解为液滴、火焰面内气相区、火焰面外气相区三个相对独立的计算域,采用固定计算域法与随动坐标系,结合耦合界面条件,进行数值求解。分析了液滴蒸发燃烧过程中的传热传质特性以及环境压力的影响。结果表明,液滴蒸发与燃烧强烈耦合,过程后期的蒸发和火焰面动态特性与初期明显不同。     

高速数字全息测量乙醇喷雾火焰中液滴三维位置、粒径及蒸发率

喷雾蒸发燃烧的研究对指导发动机燃烧系统设计具有重要意义。本文搭建了高速数字全息系统,在线测量乙醇喷雾火焰中液滴的粒径、三维位置、速度及蒸发率。对喷雾火焰中的液滴进行了统计分析,得到液滴粒径及三维空间分布。燃烧喷雾场液滴的平均粒径为68μm;非燃烧火焰测试区液滴数量多且较密集,燃烧火焰测试区液滴数量少且稀疏.追踪单液滴并处理得到湍流火焰中液滴的运动轨迹及速度。通过研究粒径的平方D2随停留时间ts的变化,测得液滴平均蒸发率为-3.343×10-7 m2/s.     

液体火箭有机凝胶喷雾液滴蒸发模型及仿真研究

凝胶推进剂虽然兼具有液体推进剂流量可控和固体推进剂长期可储存等优点, 但凝胶喷雾液滴蒸发燃烧问题却一直困扰着凝胶推进剂研制及燃烧室设计工作, 阻碍了凝胶推进剂实际工程应用.设计实现了凝胶单液滴蒸发燃烧实验系统, 通过某型有机凝胶偏二甲肼(UDMH)单液滴在四氧化二氮蒸气中的蒸发燃烧实验现象, 进一步深入分析了凝胶液滴蒸发燃烧机理.根据实验中凝胶单液滴在不同阶段的蒸发特性, 建立了有机凝胶喷雾液滴在胶凝剂膜形成、膨胀、破裂三个不同蒸发阶段的多组分蒸发模型, 采用初步选定的模型参数及物性参数对凝胶单液滴在高温气体环境中的蒸发全过程进行了仿真计算, 并与常规液体液滴的仿真结果进行了对比分析.结果表明,凝胶喷雾液滴表面胶凝剂含量在蒸发初期增加比较缓慢, 但在某临界时刻后的极短时间内迅速升高至形成胶凝剂膜的质量分数95%, 导致表面质量流率迅速下降至0,表面温度则快速上升至UDMH推进剂沸点.胶凝剂膜形成后, 液滴半径及表面UDMH蒸气质量分数出现了实验现象中凝胶液滴反复膨胀-破裂的震荡现象, 液滴表面温度维持在略高于沸点的某温度范围内,凝胶液滴内部的沸腾蒸发明显强于液体液滴表面稳态蒸发流率, 使得凝胶喷雾液滴生存时间小于常规液体液滴.     

加热基底上附壁液滴蒸发动力学实验研究

为了了解液滴表面热流型演变规律及其对蒸发速率的影响,采用稳态蒸发实验研究了加热基底上乙醇和水滴的蒸发过程。液滴半径固定在R=2.5 mm,高度变化范围为0.4至1.4 mm。结果表明,随着液滴高度降低,乙醇液滴表面出现了中心不均匀的温度波动、热流体波及Bénard-Marangoni流胞,而水滴表面则没有明显温度波动。随着基底温度升高,液滴内部热对流增强;随着液滴高度降低,乙醇液滴蒸发速率先减小、后增大,而水滴蒸发速率则无明显变化。实验结果证实热对流对液滴蒸发速率具有明显促进作用。     

旋流喷雾燃烧的直接数值模拟

对三维旋流喷雾燃烧进行了初步的直接数值模拟,其中液滴蒸发采用无限热传导蒸发模型描述,气相燃烧采用自适应单步反应机理,液滴的跟踪在拉格朗日框架中进行。模拟结果表明,喷雾燃烧的火焰结构十分复杂,仅采用传统的非预混燃烧模型是不够的。在蒸发和燃烧共存并存在强烈相互作用的区域,组分与混合物份额之间存在着复杂的关联,这给发展精确的湍流喷雾燃烧模型带来了挑战。     

气液两相旋流燃烧的时均统计特性

通过直接数值模拟的方法研究了旋流燃烧器中三维正庚烷喷雾燃烧,为实际喷雾蒸发和燃烧问题提供参考。气相燃烧模型采用自适应单步反应机理,液相采用拉格朗日方法跟踪,液滴蒸发采用无限传导蒸发模型。本文研究了气相和液相时均特性.结果发现燃烧和中心回流区(CRZ)之间有相互促进作用;同向旋流导致更强的中心回流区(CRZ),但是会有更小的外部回流区(ORZ);富预混燃烧会有更高的化学反应速率,并且蒸发冷却的影响更为明显;湍流入流会导致更高的液滴散布和更小的液滴直径。     

高焓横向射流尾迹区掺混燃烧的数值特性

尾迹区作为横向射流流场的重要结构受到广泛关注,其掺混和燃烧特性对近壁面区域的流场特性有重要影响.文章在对仿真充分验证的基础上,采用Reynolds平均模拟方法对Ma=8飞行条件下高焓横向射流尾迹区中的掺混和燃烧特性进行了数值研究.探究了冷热流场尾迹区中的氢气掺混特性,冷态流场尾迹区中的激波结构对氢气分布产生一定影响,热态流场尾迹区中存在多种氢气掺混路径.V形回流区中的高浓度氢气对燃烧产生了一定的阻碍作用.定量测量了尾迹区中的火焰结构,尾迹区火焰的顶点位置随高度增加向下游线性移动,受射流主流影响,尾迹区火焰的展向宽度在距离壁面一定高度后开始增大.对冷热流场中的主要参数进行了对比,燃烧消耗了氢气使温度升高,但是尾迹区中的流动速度没有明显增加,燃烧放出的热量没有完全转化为流体的动能.      

不同工况下旋流液雾燃烧的直接数值模拟

采用有限差分方法对不同工况下三维旋流液雾燃烧进行了直接数值模拟,其中液滴的跟踪在拉格朗日框架中进行,液滴的蒸发相变采用无限热传导蒸发模型描述,气相燃烧采用自适应单步反应机理,模拟中采用的模型燃烧器尽可能逼近真实的燃气轮机旋流燃烧器.结果发现,旋流液雾燃烧流动和火焰结构受到旋流方式和当量比的影响,流场中出现了反平行排列的...     

竖直表面液滴运动的数值模拟

为研究具有不同润湿特性的表面上液滴运动的规律,本文依据液滴接触角与基底方位角的关系以及运动液滴的动态接触角模型,运用Fluent软件,实现了对竖直表面上静止和运动液滴的数值模拟,得到了液滴的形貌以及内部的速度分布。模拟所得竖直表面上静止液滴和运动液滴的接触角和形貌均与实验结果吻合较好。液滴内部速度分布表明,裸铝表面(亲水)上液滴的运动方式为滑动;114°疏水表面上的液滴内部出现不完整的环流,但仍然以滑动运动为主;而145°疏水表面上液滴以滚动方式向下运动,液滴运动较快。     

钝体驻定湍流扩散火焰的数值研究——燃烧模型比较

分别采用标量联合的概率密度函数方法、稳态火焰面模型、Euler非稳态火焰面模型和基于有限体积/Monte Carlo混合算法的完备PDF模型对钝体驻定的Sydney湍流扩散火焰HM1进行数值模拟,以比较不同燃烧模型的性能,并比较标量联合的概率密度函数方法和Euler非稳态火焰面模型对氮氧化物排放预测的差异.计算结果和实验数据的比较表明,采用概率密度函数方法计算化学反应可以得到更好的结果但计算量较大,而用火焰面模型求解计算量较小,在接近完全燃烧的情形下,其计算结果比较合理.     

Transient convective burning of interactive fuel droplets in single-layer arrays

The transient convective burning of n-octane droplets interacting within single-layer arrays in a hot gas flow perpendicular to the layer is studied numerically, with considerations of droplet surface regression, deceleration due to the drag of the droplets, internal liquid motion, variable properties, non-uniform liquid temperature and surface tension. Infinite periodic arrays, semi-infinite periodic arrays with one row of droplets (linear array) or two rows of droplets, and finite arrays with nine droplets with centers in a plane are investigated. All arrays are aligned orthogonal to the free stream direction. This paper compares the behavior of semi-infinite periodic arrays and finite arrays with the behavior of previously studied infinite periodic arrays. Furthermore, it identifies the critical values of the initial Damköhler number for bifurcations in flame behavior at various initial droplet spacing for all these arrays. The initial flame shape is either an envelope flame or a wake flame as determined by the initial Damköhler number, the array configuration and the initial droplet spacing. The critical initial Damköhler number separating initial wake flames from initial envelope flames decreases with increasing interaction amongst droplets at intermediate droplet spacing (when the number of rows in the array increases or the initial droplet spacing decreases for a specific number of rows in the array). In the transient process, an initial wake flame has a tendency to develop from a wake flame to an envelope flame, with the moment of wake-to-envelope transition advanced for the increasing interaction amongst droplets at intermediate droplet spacing. For the array with nine droplets with centers in a plane, the droplets at different types of positions have different critical initial Damköhler number and different wake-to-envelope transition time for initial wake flame.     

Transient convective burning of a periodic fuel-droplet array

The transient convective burning of fuel-droplets interacting within 3-D infinite periodic arrays in a hot gas stream is numerically studied for the first time, with considerations of droplet regression, deceleration due to the drag of the droplets, internal liquid motion, variable properties, non-uniform liquid temperature, surface tension, and n-octane one-step oxidation kinetics. Depending upon the initial conditions and other constraints, a flame is established early as either a wake flame or an envelope flame. An initial envelope flame remains an envelope flame, and an initial wake flame has a tendency to develop from a wake flame to an envelope flame. The flame shows no strong tendency to modify significantly the standoff distance during the lifetime of the droplet. For an initial wake flame, the moment of wake-to-envelope transition is advanced as the initial droplet spacing (intermediate) is decreased, but tends to be postponed as the initial droplet spacing is further reduced. The burning rate at smaller initial droplet spacing or smaller initial Reynolds number might be greater for some period during the lifetime because of an earlier wake-to-envelope transition which elevates the average surface temperature. Lower ambient temperature yields a later wake-to-envelope transition time and smaller mass burning rate. At the lower ambient pressure with the same initial relative stream velocity, the average surface temperature is reduced, the wake-to-envelope transition is later, and the mass burning rate is smaller. Validation of our analysis is made by comparing with the results of an isolated droplet Wu and Sirignano [11].     

Large Eddy Simulation of a swirl stabilized spray flame

Large Eddy Simulations (LES) of kerosene spray combustion in an axial-swirl combustor have been carried out focusing on the effect of the evaporating droplets on the flame temperature and species concentrations. The LES-PDF methodology is used for both dispersed (liquid) and gas phases. The liquid phase is described using a Lagrangian formulation whilst an Eulerian approach is employed for the gas phase. The predictive capability of LES with sub-grid scale models for spray dispersion and evaporation is assessed placing emphasis on the effect of the unresolved velocity and temperature fields on the droplet evaporation rate. The results of the fully coupled LES formulation exhibit good agreement between the measured and simulated mean velocity fields. The global behaviour of the spray combustion, such as droplet dispersion and evaporation, are captured reasonably well in the simulations. It was found that the large velocity fluctuations observed in the shear layer strongly affect the evaporation rate and thus the temperature distributions. The present work also demonstrated the feasibility of LES to study complex flow features which are typical of gas-turbine combustion chambers.     

Influence of water mist on propagation and suppression of laminar premixed flame

The combustion of premixed gas mixtures containing micro droplets of water was studied using one-dimensional approximation. The dependencies of the burning velocity and flammability limits on the initial conditions and on the properties of liquid droplets were analyzed. Effects of droplet size and concentration of added liquid were studied. It was demonstrated that the droplets with smaller diameters are more effective in reducing the flame velocity. For droplets vaporizing in the reaction zone, the burning velocity is independent of droplet size, and it depends only on the concentration of added liquid. With further increase of the droplet diameter the droplets are passing through the reaction zone with completion of vaporization in the combustion products. It was demonstrated that for droplets above a certain size there are two stable stationary modes of flame propagation with transition of hysteresis type. The critical conditions of the transition are due to the appearance of the temperature maximum at the flame front and the temperature gradient with heat losses from the reaction zone to the products, as a result of droplet vaporization passing through the reaction zone. The critical conditions are similar to the critical conditions of the classical flammability limits of flame with the thermal mechanism of flame propagation. The maximum decrease in the burning velocity and decrease in the combustion temperature at the critical turning point corresponds to predictions of the classical theories of flammability limits of Zel'dovich and Spalding. The stability analysis of stationary modes of flame propagation in the presence of water mist showed the lack of oscillatory processes in the frames of the assumed model.     

Determination of temperature and concentration of a vapor–gas mixture in a wake of water droplets moving through combustion products

Characteristic temperatures and concentrations of a vapor–gas mixture in a wake of water droplets moving through combustion products (initial temperature 1170 K) were determined using the Ansys Fluent mathematical modeling package. We investigated two variants of motion: motion of two droplets (with sizes from 1 mm to 3 mm), consecutive and parallel, and motion of five staggered droplets. The influence of the relative position of droplets and also of distances between them (varied from 0.01 mm to 5 mm) on temperatures and concentrations of water vapor was established. The distances determine the relation between the evaporation areas and the total volume occupied by a droplet aggregate in the gas medium. The results of modeling for conditions that take into account vaporization on the droplet surface at average constant values of evaporation rate and also with consideration of the change in the latter, depending on the droplet temperature field, are compared. We determined conditions under which the modeling results are comparable for the assumption of a constant vaporization rate and with regard to the dependence of the latter on temperature. The earlier hypothesis on formation of a buffer vapor layer (“thermal protection”) around a droplet, which decreases the thermal flow from the external gas medium, was validated.     

Investigation of droplet combustion in strained counterflow diffusion flames using planar laser-induced fluorescence

Experimental investigation of an isolated droplet burning in a convective flow is reported. Acetone droplets were injected in a steady laminar diffusion counterflow flame operating with methane. Planar laser-induced fluorescence measurements applied to OH radical and acetone was used to measure the spatial distribution of fuel vapour and the structure of the flame front around the droplet. High-magnification optics was used in order to image flow areas with a ratio of 1:1.2. The different combustion regimes of an isolated droplet could be observed from the configuration of the envelope flame to that of the boundary-layer flame, and occurrence of these regimes was found to depend on the droplet Reynolds number. Experimental results were compared with 1D numerical simulations using detailed chemistry for the configuration of the envelope flame. Good agreement was obtained for the radial profile of both OH radical and fuel vapour. Influence of droplet dynamics on the counterflow flame front was also investigated. Results show that the flame front could be strongly distorted by the droplet crossing. In particular, droplets with high velocity led to local extinction of the flame front whereas droplets with low velocity could ignite within the flame front and burn on the oxidiser side. PACS 33.50.-j; 42.62.-b; 47.55.D-; 47.70.Pq; 47.80.Jk     

Vaporization and combustion in three-dimensional droplet arrays

The quasi-steady vaporization and combustion of multiple-droplet arrays is studied numerically. Utilizing the Shvab–Zeldovich formulation, a transformation of the governing equations to a three-dimensional Laplace’s equation is performed, and the solution to Laplace’s equation is obtained numerically to find the effects of droplet interactions in symmetric, multiple-droplet arrays. Vaporization rates, flame surface shapes, and flame locations are found for different droplet array configurations and fuels. The number of droplets, the droplet arrangement within the arrays, and the droplet spacing within the arrays are varied to determine the effects of these parameters. Computations are performed for uniformly spaced three-dimensional arrays of up to 216 droplets, with center-to-center spacing ranging from 3 to 25 droplet radii. As a result of the droplet interactions, the number of droplets and relative droplet spacing significantly affect the vaporization rate of individual droplets within the array, and consequently the flame shape and location. For small droplet spacing, the individual droplet vaporization rate decreases below that obtained for an isolated droplet by several orders of magnitude. A similarity parameter which correlates vaporization rates with array size and spacing is identified. Individual droplet flames, internal group combustion, and external group combustion can be observed depending on the droplet geometry and boundary conditions.     

Conditional statistics of nonreacting and reacting sprays in turbulent flows by direct numerical simulation

Conditional statistics concerning evaporation and combustion of a spray are investigated in homogeneous, isotropic, and decaying two-dimensional (2D) turbulence. Randomly distributed, polydisperse droplets of n-heptane go through single-step combustion chemistry. Attention is focused on parametric effects of initial Sauter mean radius (SMR), turbulence level and droplet velocity in both reacting and nonreacting cases. A simple linear model for the conditional evaporation rate is proposed and validated against DNS data. A conventional β-probability density function (pdf) is shown to be valid with no peak occurring on the fuel side. The amplitude mapping closure (AMC) model works well for the conditional scalar dissipation rate with evaporating and reacting sprays. Parametric study shows that initial SMR and droplet velocity are major factors affecting conditional flame structures, whereas the effect of reaction is not significant except during autoignition.     

On the behaviour of organic gel multi-size spray diffusion flames

The recently reported, experimentally observed, unusual behaviour of organic gellant-based fuel droplets which, under appropriate ambient thermal conditions, evaporate and burn in an oscillatory fashion is incorporated in a phenomenological manner in a model of a two-dimensional arbitrary multi-size spray diffusion flame. Non-unity Lewis numbers are permitted for the fuel vapour and oxidant. A combined analytical/numerical solution of the governing equations is presented and used to investigate how a spray's initial polydispersity and the frequency of oscillatory evaporation influence the combustion field. It is demonstrated that the initial droplet size distribution and the frequency of evaporation of the burning gel droplets can have an acute impact both on the homogeneous diffusion flame shape, height and width and on the thermal field downstream of the flame front. Hot spots of individual (or clusters of) burning droplets can be created and under certain operating conditions can lead to hotter temperatures than experienced in the main homogeneous flame. The intensity of these hotspots, their number and location are sensitive to spray related parameters. In realistic combustion chambers there is a danger inherent in the existence of hotspots in undesirable regions as they can damage the structural integrity. Other computed results demonstrate that, in relation to the spray diffusion flames obtained using an equivalent purely liquid fuel spray, the use of a gel fuel spray can lead, under certain operating conditions, to a reduction in flame height and temperature. The latter effect is critical when considering flame extinction.