Lagrangian conditional moment closure model with flame group interaction for lifted turbulent spray jet flames

A new Lagrangian conditional moment closure (CMC) model is developed for multiple Lagrangian groups of sequentially evaporating fuel in turbulent spray combustion. Flame group interaction is taken into account as premixed combustion by the eddy breakup (EBU) model in terms of the probability of finding flame groups in the burned and the unburned state. Evaporation source terms are included in the two phase conditional transport equations, although they turn out to have negligible influence on the mean temperature field during combustion. The Lagrangian CMC model is implemented in OpenFOAM [1 H.G. Weller, G. Tabor, H. Jasak, and C. Fureby, A tensorial approach to computational continuum mechanics using object-oriented techniques, Comput. Phys. 12 (1998), pp. 620–631.[Crossref] [Google Scholar]] and validated for test cases in the Engine Combustion Network (ECN) [2 Engine Combustion Network. (2011). Available at http://www.sandia.gov/ecn. [Google Scholar],3 L.M. Pickett, C.L. Genzale, G. Bruneaux, L. Malbec, L. Hermant, C. Christian, and J. Schramm, Comparison of diesel spray combustion in different high-temperature, high-pressure facilities, SAE technical paper 2010-01-2106 (2010). [Google Scholar]]. Similar ignition delays and lift-off lengths are predicted by the incompletely stirred reactor (ISR) and the Eulerian CMC models due to relatively uniform conditional flame structure in the domain. The improved Lagrangian CMC model shows no abrupt reaction or oscillatory behaviour with an appropriate model constant K and gives results in better agreement with measurements lying between the predictions by ISR and Lagrangian CMC without flame group interaction.     

Prediction of NO in turbulent diffusion flames using Eulerian particle flamelet model

The combustion characteristics for the turbulent diffusion flames using the unsteady flamelet concept have been numerically investigated. The Favre-averaged Navier–Stokes equations are solved by a finite volume method of SIMPLE type that incorporates the laminar flamelet concept with a modified k ? ε turbulence model. The NO formation is estimated by solving the Eulerian particle transport equations in a postprocessing mode. Two test problems are considered: CH₄/H₂/N₂ jet flame and CH₄/H₂ stabilised bluff body flame. The temperature and species profiles are well captured by the flamelet model. Two different chemical mechanisms (GRI 2.11 and 3.0) give nearly identical results for temperature and species except NO. The GRI 3.0 gives significantly higher NO levels compared to the GRI 2.11. This is mainly attributed to the difference in NO formation by the prompt mechanism. The NO formation is sensitive to the number of flamelet particles. The NO levels for two test flames do not change when the flamelet particle number exceeds six.     

Comparison of RANS/CMC modeling of Flame D with conventional and with Presumed Mapping Function statistics

Numerical simulations of Sandia Flame D are presented using Reynolds-averaged formulations plus a two equation turbulence model for the flow and mixing fields and a first order Conditional Moment Closure for the flame model. The distributions of probability and Conditional Scalar Dissipation (CSD) in mixture fraction space are modeled first using a Beta PDF (Probability Distribution Function) plus a theoretical model for CSD and then using the consistent, Presumed Mapping Function-based formulation of Mortensen for both PDF and CSD, in both two and three stream mixing modes. It is shown that there is an improvement in predictions compared with experiment when using the consistent models.     

Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k–? model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC.     

Simulation of spray development and turbulent combustion processes in low and high speed diesel engines by the CMC-ISR model

Simulation is performed to analyse the characteristics of turbulent spray combustion in conventional low and high speed diesel engine conditions. Turbulence–chemistry interaction is resolved by the Conditional Moment Closure (CMC) model in the spatially integrated form of an Incompletely Stirred Reactor (ISR). After validation against measured pressure traces, characteristic length and time scales and dimensionless numbers are estimated at the locations of sequentially injected fuel groups. Conditional flame structures are calculated for sequentially evaporated fuel groups to consider different available periods for ignition chemistry. Injection overlaps the combustion period in the high rpm engine, while most combustion occurs after injection and evaporation are complete in the low rpm engine. Ignition occurs in rich premixture with the initial peak temperature at the equivalence ratio around 2–4 as observed in Dec [2]. It corresponds to the most reactive mixture fraction of the minimum ignition delay for the given mixture states. Combustion proceeds to lean and rich sides in the mixture fraction space as a diffusion process by turbulence. The mean scalar dissipation rates (SDRs) are lower than the extinction limit to show stability of diffusion flames throughout the combustion period.     

Effect of fuel mixture fraction and velocity perturbations on the flame transfer function of swirl stabilized flames

A novel methodology is developed to decompose the classic Flame Transfer Function (FTF) used in the thermo-acoustic stability analysis of lean premix combustors into contributions of different types. The approach is applied, in the context of Large Eddy Simulation (LES), to partially-premixed and fully-premixed flames, which are stabilized via a central recirculation zone as a result of the vortex breakdown phenomenon. The first type of decomposition is into contributions driven by fuel mixture fraction and dynamic velocity fluctuations. Each of these two contributions is further split into the components of turbulent flame speed and flame surface area. The flame surface area component, driven by the pure dynamic velocity fluctuation, which is shown to be a dominant contribution to the overall FTF, is also additionally decomposed over the coherent flow structures using proper orthogonal decomposition. Using a simplified model for the dynamic response of premixed flames, it is shown that the distribution of the FTF, as obtained from LES, is closely related to the characteristics of the velocity field frequency response to the inlet perturbation. Initially, the proposed method is tested and validated with a well characterized laboratory burner geometry. Subsequently, the method is applied to an industrial gas turbine burner.     

Pulverized coal swirl flames in oxy-fuel conditions: Effects of oxidizer O2 concentration on flow field and local gas composition

In an experimental study the effects of varied oxygen concentrations in the oxidizer gas on resulting flow fields, combustion products and general behavior of pulverized coal swirl flames under oxy-fuel conditions have been investigated. Experiments were carried out in a small scale down-fired cylindrical combustion chamber equipped with an annular swirl burner. Studied flames had a constant power output of 40 kW_th and O₂/CO₂ oxidizer gas mixtures with O₂ concentrations ranging from 23 to 33 vol%. Detailed two-dimensional flow field measurements are obtained from laser Doppler anemometry (LDA). Velocity profiles (Mean and RMS) have been obtained for all conditions investigated and serve as basis for identification of flow field characteristics. Velocity RMS values are provided as supplementary material. To complement flow field measurements, in-flame gas composition measurements were also conducted using a sampling probe combined with infrared gas absorption analysis via Fourier-transform infrared (FTIR) spectrometry. The results obtained show increased velocities, particularly along the main vortex for flames with increased oxygen contents, while lower velocities are found to occur inside the recirculation regions. The opposite occurs with lower O₂ concentrations, showing significantly reduced velocities in the main vortex, but stronger recirculation than the high oxygen counterparts. This effect is attributed to a modification of the swirl level introduced by the expansion of product gases. Measured NO and CO in-flame concentrations showed significant variations under different O₂ concentrations in the oxidizer.     

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.     

Modeling and measurements of size distributions in premixed ethylene and benzene flames

This study demonstrates the major differences in the evolution of the particle size distributions (PSDs), both measured and modeled, of soot in premixed benzene and ethylene flat flames. In the experiments, soot concentration and PSDs were measured by using a scanning mobility particle sizer (SMPS, over the size range of 3-80 nm). The model employed calculations of gas phase species coupled with a discrete sectional approach for the gas-to-particle conversion. The model includes reaction pathways leading to the formation of nano-sized particles and their coagulation to larger soot particles. The particle size distribution, both experimental and modeled, evolved from a single particle mode (the nucleation mode) to a bimodal size distribution. An important distinction between the results for the ethylene and benzene flames is the behavior of the nucleation mode which persists at all heights above the burner (HAB) for ethylene whereas it was greatly suppressed at greater HAB for the benzene flames. The explanation for the decreased nucleation mode at higher elevations in the benzene flame is that the aromatics are consumed in the oxidation zone of the flame. Fair predictions of particle-phase concentrations and particle sizes in the two flames were obtained with no adjustments to the kinetic scheme. In agreement with experimental data, the model predicts a higher formation of particulate in the benzene flame as compared with the ethylene flame.     

Modelling of the correlation between velocity and reactive scalar gradients in turbulent premixed flames based on DNS data

The interaction between turbulence and reactive scalar fields is discussed for the wrinkled flamelets regime of turbulent premixed combustion. Emphasis is placed on the effects associated with the turbulent straining term. In the regime of turbulent combustion under consideration, which corresponds to Karlovitz and Damköhler numbers such that Ka < 1 and Da > 1, a clear and simple formulation is proposed to explain and to model the influence of the correlation between velocity and reactive scalar gradients. This formulation is based on the conservative variables budget across one-dimensional premixed laminar flamelets. The analysis firmly confirms the dependence on both the Damköhler number and the expansion factor, a feature already foreseen in recent studies. Nevertheless, in contrast with previous work, (i) the scaling arguments used in the present contribution are different from those used in other recent proposals, and (ii) the proposed closures are not only deduced from dimensional arguments but also from the consideration of conservative variable budgets across laminar flamelets. The resulting functional dependence on the expansion factor is found to be influenced by the underlying one-dimensional flamelet representation and two possible closures are put forward to take this dependence into account. (iii) The two closures do not exhibit a proportionality to the mean scalar dissipation rate as suggested in previous studies but to the square of this quantity. This results in the presence of a second contribution proportional to in the modelled transport equation for the mean scalar dissipation rate, in addition to the modelled molecular dissipation term. (iv) Since previous Direct Numerical Simulation (DNS) studies have been essentially devoted to the influence of the Damköhler number, the present DNS validation step is focused on the effects of the expansion rate. To this purpose, the proposed models are validated against three available DNS databases obtained for turbulent premixed flames with different values of the density ratio between unburned and fully burned gases.     

Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine

In the framework of Reynolds-averaged Navier–Stokes simulation, supersonic turbulent combustion flows at the German Aerospace Centre (DLR) combustor and Japan Aerospace Exploration Agency (JAXA) integrated scramjet engine are numerically simulated using the flamelet model. Based on the DLR combustor case, theoretical analysis and numerical experiments conclude that: the finite rate model only implicitly considers the large-scale turbulent effect and, due to the lack of the small-scale non-equilibrium effect, it would overshoot the peak temperature compared to the flamelet model in general. Furthermore, high-Mach-number compressibility affects the flamelet model mainly through two ways: the spatial pressure variation and the static enthalpy variation due to the kinetic energy. In the flamelet library, the mass fractions of the intermediate species, e.g. OH, are more sensible to the above two effects than the main species such as H₂O. Additionally, in the combustion flowfield where the pressure is larger than the value adopted in the generation of the flamelet library or the conversion from the static enthalpy to the kinetic energy occurs, the temperature obtained by the flamelet model without taking compressibility effects into account would be undershot, and vice versa. The static enthalpy variation effect has only little influence on the temperature simulation of the flamelet model, while the effect of the spatial pressure variation may cause relatively large errors. From the JAXA case, it is found that the flamelet model cannot in general be used for an integrated scramjet engine. The existence of the inlet together with the transverse injection scheme could cause large spatial variations of pressure, so the pressure value adopted for the generation of a flamelet library should be fine-tuned according to a pre-simulation of pure mixing.     

时域有限差分法在建筑声学中的应用及前景

介绍了声波方程的基本差分格式及稳定条件、数值色散、吸收边界条件等数值计算理论，例举了前人用时域有限差分法对噪声传播过程的模拟和室内声学中座椅吸声低谷效应模拟的模拟结果。本文指出，由于时域有限差分法的特点使其具有在模拟脉冲响应方而的特别优势，因而，应用这一技术研究厅堂的声学特性，尤其是低频特性，将有重要意义。     

细菌觅食优化算法在光散射法颗粒粒度反演中的应用

针对函数约束算法中传统的智能算法反演时存在鲁棒性差和易陷入局部最优的缺点,提出了将正则化理论与细菌觅食优化算法相结合应用在颗粒粒度的测量中。引入Tikhonov平滑泛函来构建算法的目标函数,采用L曲线法确定正则化参数;再利用细菌觅食优化算法通过趋向、聚群、复制和迁徙等四种智能行为,迭代计算来搜寻函数的最优解。实验仿真结果表明:利用细菌觅食优化算法实现了在不同程度的随机噪声下的服从J-SB分布的单峰分布的均匀球形颗粒粒度分布反演,其反演结果更稳定,反演精度高,对于实现稳定、快速、准确的颗粒粒度在线测量具有重要的意义。     

Analysis on the efficiency of parallel FDTD method and its application in two-dimensional photonic crystal

A new kind of parallel finite-difference time-domain (FDTD) method is used to calculate and analyze the transmission properties of photonic crystal. The effectiveness and the accuracy of parallel FDTD method were proved by the typical two-dimensional sine wave. The result shows that parallel FDTD method, which can save time effectively in electromagnetic problems with large size and long time, performs as accurate as the series one. This method, whose results shows that the photonic band gap normalized frequency moves to the low-frequency direction with the enlargement of the section area of the square cylinder and that the band gap gets thinner, is applied to simulate and analyze the two-dimensional photonic crystal and gives calculating method and analyzing accordance to calculate and analyze large-sized photonic crystal structure.     

表面增强拉曼光谱对免疫球蛋白IgG分子与银基底作用的研究

采用高灵敏度的表面增强拉曼光谱(SERS)技术,以具有强SERS信号的金纳米粒子标记抗体,以此SERS标记免疫金溶胶为探针,结合扫描电镜技术,研究免疫球蛋白羊抗小鼠IgG分子与银基底的相互作用。我们发现,羊抗小鼠IgG分子可直接与银基底通过疏水作用或形成Ag-S键而牢固结合。为消除这种非特性吸附,本文以小牛血清白蛋白(BSA)封闭银基底,取得了较好的效果。     

小波变换和连续投影算法在火龙果总酸无损检测中的应用

应用可见/近红外光谱技术、小波变换(WT)和连续投影算法(SPA),对火龙果总酸含量(TA)进行精确、快速的无损检测,为火龙果内部品质无损检测提供科学依据。利用Maya2000光纤光谱仪采集380～1 099 nm范围的火龙果漫反射光谱数据,通过WT消噪、SPA优选波长和偏最小二乘回归(PLSR)分析方法,建立了火龙果总酸的定量预测模型。试验结果表明：经过WT消噪联合SPA优选波长压缩光谱变量后建立的WT-SPA-PLSR模型,预测精度都高于全谱PLSR模型。由全部样本的原始光谱变量作为输入变量建立PLSR模型的预测相关系数(R_p)为0.851 394, 预测均方根误差(RMSEP)为0.086 848;全部样本的原始光谱数据使用dbN(N=2,3,…,10)小波进行分解消噪,其中消噪效果最优的是db4小波2层分解(db4-2),WT-PLSR模型的R_p为0.915 635,RMSEP为0.066 752,小波变换消噪后的光谱预测模型精度明显提高;原始光谱经过db10-3小波消噪联合SPA算法,从570个光谱变量中优选出530,545,604,626,648,676,685,695,730,897,972,1 016 nm共12个变量作为输入变量,建立WT-SPA-PLSR预测模型,模型的R_P为0.882 83, RMSEP为0.077 39。SPA算法适合火龙果TA模型的光谱变量选择,能够有效提取与总酸相关度高的波长变量,增加了预测模型的精度和稳定性。研究结果表明小波变换技术联合连续投影算法的漫反射近红外光谱无损检测火龙果总酸含量具有可行性。     

非分离状态下FI-KR-FAAS测定铁的价态

建立了流动注射(FI)编结反应器(KR)在线预富集与火焰原子吸收光谱(FAAS)联用非分离状态下测定铁的价态的新方法。进样流速6.0 mL.min-1;进样时间60s,测定0.04 mg.L-1的Fe(Ⅲ)和Fe(Ⅱ),Fe(Ⅲ)和Fe(Ⅱ)的浓集系数(EF)分别为41和9;检出限分别为2.5和14.3μg.L-1;相对标准偏差(RSD,n=11)分别为2.3%和3.1%。以0.1‰φ的三乙醇胺为掩蔽剂,Fe(Ⅲ)和Fe(Ⅱ)在水样中的回收率分别为97%～101%和96%～100%。     

光谱分析在碲化镉晶体制备和晶体评价中的应用

籍助于原子吸收法及X射线光电子能谱(XPS)的分析,比较了CdTe晶体制备的不同工艺条件。根据红外光谱和光致发光光谱的测量结果,对不同方法制备的CdTe晶体质量作了评价。     

时域有限差分法在光波导分析中的应用及改进

对时域有限差分算法在圆柱坐标下弱导光纤中应用时吸收边界条件的设置及其对计算精度的影响进行了讨论，提出了改善的方法。并且将这种方法应用于一般光纤和光电子晶体光纤传输特性的计算，得到了有用结果。     

近红外光谱快速分析技术及其在动物饲料和产品品质检测中的应用

近红外光谱是近年来发展最快、最引人注目的光谱学技术。主要介绍了近红外光谱技术的基本原理和发展, 包括近红外光谱预处理技术如微分处理、信号平滑等技术的发展和近红外光谱分析模型如多元线性回归、主成分分析、偏最小二乘法和人工神经网络等的发展。综述了国内外近几年来此技术在动物饲料和产品品质检测中的应用。文献调查显示, 近红外光谱分析技术以其快速、无损、不污染环境等诸多优点在国内外饲料和动物产品检测方面得到广泛应用。在饲料分析方面,近红外不仅能用于其常量成分干物质、粗蛋白、粗纤维、粗脂肪等的测定,而且能用于微量成分、有毒有害成分的测定。在动物产品分析方面,该技术已用于禽蛋、牛肉、羊肉、猪肉等的各种物理和化学指标的测定。文中详细给出了已经报道的利用近红外光谱技术测定饲料和动物产品测定指标和光谱处理以及模型建立的情况, 并讨论了近红外光谱快速检测技术的在饲料分析和动物产品分析领域的应用新趋势和局限性。