Impact of Turbulent Flow and Mean Mixture Fraction Results on Mixing Model Behavior in Transported Scalar PDF Simulations of Turbulent Non-premixed Bluff Body Flames

Numerical simulation results are presented for three turbulent jet diffusion flames, stabilized behind a bluff body (Sydney Flames HM1-3). Interaction between turbulence and combustion is modeled with the transported joint-scalar PDF approach. The focus of the study is on the impact of the quality of simulation results in physical space on the behavior of two micro-mixing models in composition space: the Euclidean Minimum Spanning Tree (‘EMST’) model and the modified Curl coalescence dispersion (‘CD’) model. Profiles of conditional means and variances of thermo-chemical quantities, conditioned on the mixture fraction, are discussed in the recirculation region and in the neck zone behind. The impact of the flow and mixing fields in physical space on the mixing model behavior in composition space is strong for the CD model and increases as the turbulence – chemistry interaction becomes stronger. The EMST conditional profiles, on the contrary, are hardly affected.     

Conditional Source-Term Estimation as a Method for Chemical Closure in Premixed Turbulent Reacting Flow

A Conditional Source-term Estimation (CSE) model is used to close the mean reaction rates for a turbulent premixed flame. A product-based reaction progress variable is introduced as the conditioning variable for the CSE method. Different presumed probability density function (PDF) models are studied and a modified version of a laminar flame-based PDF model is proposed. Improved predictions of the variable distribution are obtained. The conditional means of reactive scalars are evaluated with CSE and compared to the direct numerical simulation (DNS). The mean reaction rates in a turbulent premixed flame are evaluated with the CSE model and the presumed PDFs. Comparison of the CSE closure method to DNS shows promising results. This paper was presented at the 2nd ECCOMAS Thematic Conference on Computational Combustion.     

Influence of Differential Diffusion on Super-Equilibrium Temperature in Turbulent Non-Premixed Hydrogen/Air Flames

In this paper we wish to investigate the occurrence of super-equilibrium temperature values, observed in many experimental configurations. We would like to understand the origin of this phenomenon. Previous authors have already shown that differential diffusion can lead to considerable changes in the temperature field and we would like to build on top of this observation. We investigate numerically super-equilibrium combustion by considering both laminar counter-flow and turbulent diluted hydrogen/air diffusion flames. These turbulent flames are computed using direct numerical simulations (DNS). A detailed reaction mechanism is employed and the transport properties are modeled using multicomponent diffusion velocities, including the Soret effect. Analyzing these results we introduce three complementary parameters (dilution-free mixture fraction, dilution excess and local enthalpy) to describe the local combustion conditions. Introducing a measure of dilution separately from the mixture fraction is necessary for a proper analysis. Using this set of parameters it becomes possible to explain super-equilibrium temperature levels as a consequence of differential diffusion.     

A Priori Evaluation of the Laminar Flamelet Decomposition Model for Turbulent Premixed Flames using DNS Data

Laminar flamelet decomposition (LFD) is a dynamic approach for modelling sub-filter scale turbulence-chemistry interactions in Large-Eddy Simulations using a stretched flamelet library. In this work, the performance of the LFD model – that was previously used only in non-premixed combustion—is investigated a priori for premixed combustion using positively-strained flamelets in the reactant-to-product configuration. For this purpose, a DNS database of methane-air premixed flames is utilized. The flames are propagating in a rectangular box under homogeneous isotropic turbulence conditions over a wide range of Karlovitz numbers. The results show that the LFD model can correctly account for the sub-filter scale turbulence-chemistry interactions to predict the filtered reaction rates and the filtered scalar field, provided that turbulent and laminar mixing are well predicted. The deviations from the DNS results are attributed to the shortcomings of the strained flamelet library and the non-flamelet effects. Finally, the LFD results are compared with a different sub-filter scale model using the same strained flamlelet library, and the relative performances of the two models are discussed.—

     

A Comparison of Strategies for Direct Numerical Simulation of Turbulence Chemistry Interaction in Generic Planar Turbulent Premixed Flames

Three different methods to introduce turbulence in the computational domain of Direct Numerical Simulations (DNS) of statistically planar turbulent premixed flame configurations have been reviewed and their advantages and disadvantages in terms of run time, natural flame development, control of turbulence parameters and convergence of statistics extracted from the simulations have been discussed in detail. It has been found that there is no method, which is clearly superior to the other two alternative methods. An analysis has been performed to explain why Lundgren’s physical space linear forcing results in an integral length scale which is, independent of the Reynolds number, a constant fraction of the domain size. Furthermore, an evolution equation for the integral length scale has been derived, and a scaling analysis of its terms has been performed to explain the evolution of the integral length scale in the context of Lundgren’s physical space linear forcing. Finally, a modification to Lundgren’s forcing approach has been suggested which ensures that the integral length scale settles to a predetermined value so that DNS of statistically planar turbulent premixed flames with physical space forcing can be conducted for prescribed values of Damköhler and Karlovitz numbers.     

A Study of Strain Rate Effects for Turbulent Premixed Flames with Application to LES of a Gas Turbine Combustor Model

Large-scale strain rate field, a resolved quantity which is easily computable in large-eddy simulations (LES), could have profound effects on the premixed flame properties by altering the turbulent flame speed and inducing local extinction. The role of the resolved strain rate has been investigated in a posterior LES study of GE lean premixed dry low-NOx emissions LM6000 gas turbine combustor model. A novel approach which is based on the coupling of the linear-eddy model with a one-dimensional counterflow solver has been applied to obtain the parameterizations of the resolved premixed flame properties in terms of the reactive progress variable, the local strain rate measure, and local Reynolds and Karlovitz numbers. The strain rate effects have been analyzed by comparing LES statistics for several models of the turbulent flame speed, i.e, with and without accounting for the local strain rate effects, with available experimental data. The sensitivity of the simulation results to the inflow velocity conditions as well as the grid resolution have been also studied. Overall, the results obtained demonstrate that the effects of the resolved strain rate are not dominant for the considered premixed flame configuration and the unstrained turbulent flame speed model is found to perform as well as the one that allows for the strain rate effects.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

基于能力谱法的土-桩-结构相互作用分析

本文基于能力谱方法对天然地基上带群桩基础的高层框架结构进行了平面静力弹塑性分析.文中首先应用子结构原理,将群桩基础用一组弹簧来模拟,并用Davies方法确定群桩刚度,建立了土-桩-高层框架结构相互作用体系Pushover分析模型;在此基础上,分析了均匀分布方式和依各层质量与第一弹性振型乘积比例分布等两种侧向加载模式下土-桩-结构相互作用模型和刚性地基模型的结构各性能点的加速度谱值Sa和位移谱值Sd、以及基底剪力和顶点位移等,探讨了不同地震烈度(7度、8度和9度)下两种不同分析模型间的结果差异,得到了一些有应用价值的结果.     

圆柱体形地基模型及与结构相互作用的样条函数解法

本文根据工程中大量遇到的圆形截面结构,建立了圆柱体形地基模型,并用B样条基函数来表征模型的位移。由于所设函数的合理性,使在一般荷载作用下圆心处的位移保持连续,应力应变不再奇异,且具有待定系数少、精度高等特点,使庞大的三维问题能在一般微机上求解。文中对地基-结构相互作用时地基人工边界的设置及范围进行了大量计算,得到了一些规律性的结果。     

管道流固耦合振动的行波方法研究

应用波动方法分析载流管道流固耦合振动问题。采用直梁模型，推导了管道流固耦合系统的轴向、横向波导方程及平面管系的坐标转换矩阵，通过节点位移连续和力平衡条件建立了耦合系统的散射模型。最终分析了两种耦合对载流管道振动的影响机理及流速变化对管道振动稳定性的影响。数值算例表明本文方法的正确性及有效性。     

A Numerical Study Promoting Algebraic Models for the Lewis Number Effect in Atmospheric Turbulent Premixed Bunsen Flames

This numerical investigation carried out on turbulent lean premixed flames accounts for two algebraic – the Lindstedt–Vaos (LV) and the classic Bray–Moss–Libby (BML) – reaction rate models. Computed data from these two models is compared with the experimental data of Kobayashi et al. on 40 different methane, ethylene and propane Bunsen flames at 1 bar, where the mean flame cone angle is used for comparison. Both models gave reasonable qualitative trend for the whole set of data, in overall. In order to characterize quantitatively, firstly, corrections are made by tuning the model parameters fitting to the experimental methane–air (of Le = 1.0) flame data. In case of the LV model, results obtained by adjusting the pre-constant, i.e., reaction rate parameter, C_R, from its original value 2.6 to 4.0, has proven to be in good agreement with the experiments. Similarly, for the BML model, with the tuning of the exponent n, in the wrinkling length scale, L_y = C_l⋅ l_x(s_L/u′)ⁿ from value unity to 1.2, the outcome is in accordance with the measured data. The deviation between the measured and calculated data sharply rises from methane to propane, i.e., with increasing Lewis number. It is deduced from the trends that the effect of Lewis number (for ethylene–air mixtures of Le = 1.2 and propane–air mixtures of Le = 1.62) is missing in both the models. The Lewis number of the fuel–air mixture is related to the laminar flame instabilities. Second, in order to quantify for its influence, the Lewis number effect is induced into both the models. It is found that by setting global reaction rate inversely proportional to the Lewis number in both the cases leads to a much better numerical prediction to this set of experimental flame data. Thus, by imparting an important phenomenon (the Lewis number effect) into the reaction rates, the generality of the two models is enhanced. However, functionality of the two models differs in predicting flame brush thickness, giving scope for further analysis.     

A Similarity Subgrid Model for Premixed Turbulent Combustion

The accuracy of large-eddy simulation (LES) of a turbulent premixed Bunsen flame is investigated in this paper. To distinguish between discretization and modeling errors, multiple LES, using different grid sizes h but the same filterwidth Δ, are compared with the direct numerical simulation (DNS). In addition, LES using various values of Δ but the same ratio Δ/h are compared. The chemistry in the LES and DNS is parametrized with the standard steady premixed flamelet for stochiometric methane-air combustion. The subgrid terms are closed with an eddy-viscosity or eddy-diffusivity approach, with an exception of the dominant subgrid term, which is the subgrid part of the chemical source term. The latter subgrid contribution is modeled by a similarity model based upon 2Δ, which is found to be superior to such a model based upon Δ. Using the 2Δ similarity model for the subgrid chemistry the LES produces good results, certainly in view of the fact that the LES is completely wrong if the subgrid chemistry model is omitted. The grid refinements of the LES show that the results for Δ = h do depend on the numerical scheme, much more than for h = Δ/2 and h = Δ/4. Nevertheless, modeling errors and discretization error may partially cancel each other; occasionally the Δ = h results were more accurate than the h ≤ Δ results. Finally, for this flame LES results obtained with the present similarity model are shown to be slightly better than those obtained with standard β-pdf closure for the subgrid chemistry.     

一种紊流润滑理论分析新方法—复合型紊流模式理论

在对主要用于高速轻载工况的常用紊理论进行简要分析和充分考察润滑流场的边界条件及内部结构的基础上，采用理论上比现有紊流润滑理论更为合理的复合型紊流模式理论，即在近壁区采用低紊流雷诺数的ｋ－ε模式，而在紊流核心区采用代数雷诺应力模式，对复杂流场的紊流润滑进行了分析，同时在计入惯性效应的情况下，推导出了一种适用于高压密封和高速重载轴承等计算，计算结果与实验数据十分吻合，验证了模型的有效性，可以应用于高压     

Flow Regimes of Interaction of a Turbulent Plane Jet into a Rectangular Cavity: Experimental Approach and Numerical Modelling

The present work is devoted to the experimental and numerical study ofthe interaction of a turbulent plane jet with a rectangular cavity.Several flow regimes have been found to occur: the non-oscillationregime, the stable oscillation regime and an unstable oscillationregime. The first two regimes have been particularly considered. Theexperimental study has been carried out using hot wire anemometry andsome visualisations. The numerical predictions based on statisticalmodelling have been made using on the one hand the standard k– model and on the other hand a two-scales split spectrum model. The structuralproperties of the flow have been described for the different situations.For the oscillatory regime, a parametrical study allowed to determinethe influence of the jet exit location and the Reynolds number on thefrequency of the jet flapping. The one point closures have been able topredict the oscillatory regime, and in particular the two-scales modelled to improved results because better account is taken of lag effectsin unsteady non-equilibrium situations.     

A Refined Asymptotic Model of Fluid-Structure Interaction in Scattering by Elastic Shells

A refined asymptotic model of fluid-structure interaction in scattering by elastic shells is proposed. The model takes into consideration transverse compression of a shell by a fluid and some other phenomena. As an illustration, scattering of a plane acoustic wave by a circular cylindrical shell is considered. Comparison of numerical data corresponding various approximate approaches is provided. This revised version was published online in July 2006 with corrections to the Cover Date.     

Algorithms for a Vehicle Dynamics Monitoring System,Based on Model Reference Structure

Vehicle control depends heavily on the knowledge of the vehicle operatingconditions. One of the most important parameters for its control is thetyre–road friction coefficient (µ). An appropriate way toestimate the vehicle operating conditions is the Model Reference Approach.This technique requires a model that provides estimated states which can becompared with the measured states, the difference is used to determine thereal operating conditions. This paper presents two different applications ofthe Model Reference Techniques to estimate tyre–road frictioncoefficient; these are based on the relation between tyre forces and slip,and on the vehicle lateral behaviour using an extended Kalman filter.Experimental data from the test vehicle confirms the good results obtainedin the friction estimation based on the tyre slip–force relation. Theestimation using an extended Kalman filter on lateral behaviour showsaccurate tracking. The next step to be taken is to integrate all thealgorithms in the test vehicle and to validate them for a wide range ofoperating conditions, in order to have reliable information for the activesystem control.     

Experimental Study on Stabilization of Lifted Swirl Flames in a Model GT Combustor

This study reports on experimental investigations on isothermal and reacting swirled non-premixed flows under varying pressure conditions. In this configuration, a central high speed fuel jet was surrounded by a heated swirling air flow. For the reacting case natural gas served as fuel whereas for isothermal conditions fuel was replaced by a mixture of helium and air to achieve Reynolds-similarity. The optically accessible combustor allowed for application of laser diagnostics. Here we report on Laser Doppler Anemometry and planar laser-induced fluorescence (PLIF) experiments used to characterize the flow field and visualize selected scalars, respectively. Acetone served as a fluorescence marker for mixture fraction investigations. The hydroxyl radical was used to provide general features of the reaction zone such as flame shape and mean stabilization. To expose the influence of pressure on the flame structure three different operating points were investigated varying the combustor pressure between 2 and 6 bar while the inflow bulk velocities remained the same. Striking features of the present configuration are a detached flame, multiple recirculation zones, and complex coherent flow structures.     

Vorticities in a LES Model for 3D Periodic Turbulent Flows

This paper is devoted to the study of a LES model to simulate turbulent 3D periodic flow. We focus our attention on the vorticity equation derived from this LES model for small values of the numerical grid size δ. We obtain entropy inequalities for the sequence of corresponding vorticities and corresponding pressures independent of δ, provided the initial velocity u₀ is in L_x² while the initial vorticity ω₀ = ∇ × u₀ is in L_x¹. When δ tends to zero, we show convergence, in a distributional sense, of the corresponding equations for the vorticities to the classical 3D equation for the vorticity.