The burning rate’s effect on the flame length of weak fire whirls

This paper discusses why the visibly-determined flame length of a weak fire whirl increases as compared with the corresponding pool fire without spin. Here, a fire whirl is called weak when the pure aerodynamic effect of flow circulation has a negligible influence on the flame length. Split cylinders were used to apply a flow circulation to a 3-cm-diameter methane burner flame and a 3-cm-diameter ethanol pool fire. After applying the flow circulation, the flame length of the ethanol pool fire increased about three times, while little change was observed in the flame length of the methane burner flame. The difference is explained by the fact that the burning rate of the methane burner flame was fixed constant, whereas that of the ethanol pool fire increased due to the increased heat input to the fuel surface caused by a change in flame shape pushed toward the fuel surface. The experimental observations thus demonstrate that the burning-rate effect can significantly increase the flame length even under a weak circulation condition. Computational fluid dynamics (CFD) simulations were conducted to understand the detailed flow structure of a fire whirl. An analytical model was then developed based on the experimental observations and CFD calculations; the predicted relationship between the flame height and the burning rate agreed with experimental data.     

Large eddy simulation of turbulent combustion systems

This paper reviews recent and ongoing work on numerical models for turbulent combustion systems based on a classical LES approach. The work is confined to single-phase reacting flows. First, important physico-chemical features of combustion-LES are discussed along with several aspects of overall LES models. Subsequently, some numerical issues, in particular questions associated with the reliability of LES results, are outlined. The details of chemistry, its reduction, and tabulation are not addressed here. Second, two illustrative applications dealing with non-premixed and premixed flame configurations are presented. The results show that combustion-LES is able to provide predictions very close to measured data for configurations where the flow is governed by large turbulent structures. To meet the future demands, new key challenges in specific modelling areas are suggested, and opportunities for advancements in combustion-LES techniques are highlighted. From a predictive point of view, the main target must be to provide a reliable method to aid combustion safety studies and the design of combustion systems of practical importance.     

Parametric studies on hydrocarbon fireball using large eddy simulations

Occurrences of fireball close to plant buildings due to the release of flammable hydrocarbon fuel caused by the formation of fuel vapour cloud poses severe safety concerns. On the availability of potential ignition source, the induced fireball would cause the damage to the structures of nuclear power plant by direct contact, radiation and/or convection of hot combustion products through the opening of air intakes and ducts. In the present paper, the accidental/ experimental observations and theoretical studies of fireball are summarised. Computational fluid dynamics (CFD) analyses have been carried out to study the behaviour of fireball using OpenFOAM CFD software. The parametric studies are conducted by varying the mass of fuel, inlet velocity and inlet diameter. The new correlations for fireball diameter and duration have been proposed based on the parametric studies using CFD simulations. The fireball with a larger amount of fuel releases the heat slower and for a longer duration. The high heat released rate (HRR) is observed in case of a larger inlet diameter used for the same mass. The incident radiation from the fireball is calculated at different locations to assess thermal hazard. Analysis performed show that various parameters like fireball diameter, duration and the radiative flux falling at different locations can be predicted well using CFD code.     

Large eddy simulation of swirling particle-laden flow in a model axisymmetric combustor

This paper focuses on the application of the large eddy simulation (LES) technique to a swirling particle-laden flow in a model combustion chamber. A series of calculations have been performed and compared directly with detailed experimental measurements. The computational domain identically matches the laboratory configuration, which effectively isolates effects related to dilute particle dispersion and momentum coupling. Results highlight the predictive capabilities of LES when implemented with the appropriate numerics, grid resolution (as dictated by the class of models employed) and well-defined boundary conditions. The case study provides a clearer understanding of the effectiveness and feasibility of current state-of-the-art models and a quantitative understanding of relevant modeling issues by analyzing the characteristic parameters and scales of importance. The novel feature of the results presented is that they establish a baseline level of confidence in our ability to simulate complex flows at conditions representative of those typically observed in gas-turbine (and similar) combustors.     

Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.     

Large eddy simulation of stably stratified turbulence

Stable stratification turbulence, as a common phenomenon in atmospheric and oceanic flows, is an important mechanism for numerical prediction of such flows. In this paper the large eddy simulation is utilized for investigating stable stratification turbulence numerically. The paper is expected to provide correct statistical results in agreement with those measured in the atmosphere or ocean. The fully developed turbulence is obtained in the stable stratification fluid by large eddy simulation with different...     

Large eddy simulation of city micro-atmospheric environment

Air quality is one of the important conditions for a better residence life in the populated urban area and it is closed related to the micro-atmospheric environment. Atmospheric environment is controlled by air motion with multi-scales in the city, while air flows in the residence area are of micro-scale atmospheric motion. This paper introduces a modern numerical simulation method, i.e. large eddy simulation (LES), for studying micro-atmospheric flows in the city residence area. For the complex flow features in the residence area, the proper application of LES is studied and various numerical methods are compared in order to investigate their effects on the prediction accuracy of micro-atmospheric flows, for instance, roughness elements and immersed boundary method for complex terrain, different subgrid models and so on. The wind field (including turbulence properties) and contaminant dispersion are computed by the proposed method for a model and a realistic residence area, and the numerical results are in good agreement with the experimental measurements. Supported by the National Natural Science Foundation of China (Grant No. 10572073) and Foundation for Development of Science and Technology in Macao (Grant No. 022/2006/A)     

Burst detection in turbulent channel flows based on large eddy simulation databases

As one of the important coherent structures in the near-wall region, turbulent burst is responsible for the production and transport of major turbulent kinetic energy and Rey- nolds stress[1]. Nearly half of turbulent kinetic energy or Reynolds stress is produced in the near-wall region, and 80% flows in outer region only contribute 20% of them. Both ejection and sweeping events contribute 60―70% of the turbulent shear stress respec- tively[2]. Recently, turbulent burst process has been foun…     

Evaluating the modulated gradient model in large eddy simulation of channel flow with OpenFOAM

Recently, a new family of subgrid-scale (SGS) models, termed as gradient-based models, has been introduced to calculate the SGS stresses in large eddy simulation (LES). In the present work, the modulated gradient model (MGM) was implemented in the OpenFOAM package, and the pimpleFoam solver was improved to be adopted with non-eddy viscosity models. The MGM is a new, nonlinear model that uses the local equilibrium hypothesis to assess the SGS kinetic energy and the velocity gradient tensor to calculate the relative weight of the different components of the SGS stress tensor. To evaluate the accuracy of the MGM along with the modified pimpleFoam solver, a turbulent channel flow was simulated at the three different frictional Reynolds numbers of 180, 395 and 590. Furthermore, the results were compared with direct numerical simulation data, as well as the numerical results obtained by the established SGS models such as the dynamic Smagorinsky model (DSM). A suitable accuracy for the first- and second-order turbulence parameters was reported. Moreover, it was demonstrated that MGM is computationally efficient compared to the DSM in treating channel flow.     

Large eddy simulation of turbulent horizontal buoyant jets

Large eddy simulations (LESs) of turbulent horizontal buoyant jets are carried out using a high-order numerical method and Sigma subgrid-scale (SGS) eddy-viscosity model, for a number of different Reynolds (Re) and Richardson (Ri) numbers. Simulations at previous experimental flow conditions (Re = 3200, 24, 000 and Ri = 0, 0.01) are carried out first, and the results are found to be qualitatively and quantitatively similar to the experimental results, thus validating the numerical methodology. The effect of varying Ri (values 2×10^?4, 0.001, 0.005, and 0.01) and Re (3200 and 24, 000) is studied next. The presence of stable stratification on one side and unstable stratification on the other side of the jet centreline leads to an asymmetric development of horizontal buoyant jets. It is found that this asymmetry, the total radial spread and the vertical deflection are significantly affected by Ri, while Re affects only the radial asymmetry. The need for developing improved integral models, accounting for this asymmetry, is pointed out. Turbulent production and dissipation rates are investigated, and are found to be symmetric in the horizontal plane, but asymmetric in the mid-vertical plane. A previously proposed model, for correlation between the vertical component of the fluctuating scalar flux vector and the vertical cross-correlation component of the Reynolds tensor, is modified based on the current LES results. Instantaneous scalar and velocity fields are analysed to reveal the structure of horizontal buoyant jets. Similar to the developed turbulent jet, the flow close to the nozzle too is found to be markedly different in the stable and unstable stratification regions. Persistent coherent vortex rings are found in the stable stratification region, while intermittent breakdown of vortex rings into small-scale structures is observed in the unstable stratification region. Similarities and differences between the flow structures in the horizontal buoyant jet configuration and those in the jet in crossflow configuration are discussed. Finally, a dynamic mode decomposition analysis is carried out, which indicates that the flow in the unstable stratification region is more energetic and prone to instabilities, as compared to the flow in the stable stratification region.     

A fractal flame-wrinkling large eddy simulation model for premixed turbulent combustion

Combustion plays an important role in a wide variety of industrial applications, such as gas-turbines, furnaces, spark-ignition engines, and various air-breathing engines. The ability to predict and understand the behavior of reacting flows in practical devices is fundamental to improved combustors with higher efficiency and reduced levels of emissions. At present, large eddy simulation is considered the most promising approach for premixed combustion modeling since the large-scale energy containing flow structures are resolved on the grid. However, the typically thin reaction zone cannot be resolved. To overcome this difficulty flamelet models, in which the reaction is assumed to take place in thin layers, wrinkled by the turbulence can sometimes be used. In these models, the turbulent flame speed can be represented as the product of the laminar flame speed, S_u, corrected for the effects of stretch (strain and curvature) and the flame-wrinkling, Ξ. In this study, we propose to model Ξ using fractal theory. This model requires sub-models for the fractal dimension, and the inner and outer cut-offs—the latter being set by the grid. A model is proposed for the inner cut-off, whereas an empirical parameterization is used to provide the fractal dimension. The proposed model is applied to flame kernel growth in homogeneous isotropic turbulence in a fan-stirred bomb and to a lean premixed flame in a plane symmetric dump combustor. Good qualitative and quantitative agreement with experimental data were obtained for the proposed model in both cases. Comparison with other well-known turbulent flame speed closure models shows that the proposed model behaves at least as good, or even better, than the reference models.     

Stochastic modeling of atomizing spray in a complex swirl injector using large eddy simulation

Large-eddy simulation of an atomizing spray issuing from a gas-turbine injector is performed. The filtered Navier–Stokes equations with dynamic subgrid scale model are solved on unstructured grids to compute the swirling turbulent flow through complex passages of the injector. The collocated grid, incompressible flow algorithm on arbitrary shaped unstructured grids developed by Mahesh et al. (J. Comp. Phys. 197 (2004) 215–240) is used in this work. A Lagrangian point-particle formulation with a stochastic model for droplet breakup is used for the liquid phase. Following Kolmogorov’s concept of viewing solid particle-breakup as a discrete random process, the droplet breakup is considered in the framework of uncorrelated breakup events, independent of the initial droplet size. The size and number density of the newly produced droplets is governed by the Fokker–Planck equation for the evolution of the pdf of droplet radii. The parameters of the model are obtained dynamically by relating them to the local Weber number and resolved scale turbulence properties. A hybrid particle-parcel is used to represent the large number of spray droplets. The predictive capability of the LES together with Lagrangian droplet dynamics models to capture the droplet dispersion characteristics, size distributions, and the spray evolution is examined in detail by comparing it with the spray patternation study for the gas-turbine injector. The present approach is computationally efficient and captures the global features of the fragmentary process of liquid atomization in complex configurations.     

Large eddy simulation of bluff-body stabilized swirling non-premixed flames

Large eddy simulations (LES) using a subgrid mixing and combustion model are carried out to study two bluff-body stabilized swirling non-premixed flames (SM1 and SMA2). The similarities and differences between the two flames are highlighted and discussed. Flow features, such as, the recirculation zone (RZ) size and the flame structure are captured accurately in both cases. The SM1 flame shows a toroidal RZ just behind the bluff body and a vortex breakdown bubble (VBB) downstream. In addition, a highly rotational non-recirculating region in-between the RZ and VBB is observed as well. On the other hand, the SMA2 shows a single elongated recirculation zone downstream the bluff body. Flame necking is observed downstream the bluff body for the SM1 flame but not for the SMA2 flame. The time-averaged velocity and temperature comparison also shows reasonable agreement. The study shows that the sensitivity of the flame structure to inflow conditions can be captured in the present LES without requiring any model changes.     

Large eddy simulation of the turbulent round jet dynamics

Large Eddy Simulation (LES) of turbulent free round jet submerged in water has been performed at Reynolds number Re = 25 000. LES results are compared against the PIV measurements data for the same flow configuration. The processes of the mixing layer formation in the jet initial region and its evolution down-stream have been studied with the analysts of turbulent power spectra. The effect of the subgrid Smagorinsky model coefficient C _s on the jet hydrodynamics has been studied. The work has been supported by the Russian Foundation for Basic Research (grants Nos. 03-02-16708-a and 04-02-16907-a).     

Large eddy simulations of a circular orifice jet with and without a cross-sectional exit plate

The effect of a cross-sectional exit plane on the downstream mixing characteristics of a circular turbulent jet is in- vestigated using large eddy simulation （LES）. The turbulent jet is issued from an orifice-type nozzle at an exit Reynolds number of 5 ×104. Both instantaneous and statistical velocity fields of the jet are provided. Results show that the rates of the mean velocity decay and jet spread are both higher in the case with the exit plate than without it. The existence of the plate is found to increase the downstream entrainment rate by about 10% on average over the axial range of 8-30de （exit diameter）. Also, the presence of the plate enables the formation of vortex rings to occur further downstream by 0.5-1 .Ode. A physical insight into the near-field jet is provided to explain the importance of the boundary conditions in the evolution of a turbulent jet. In addition, a method of using the decay of the centreline velocity and the half-width of the jet to calculate the entrainment rate is proposed.     

化学氧碘激光器内流场亚跨超音速混合的大涡模拟研究

 利用大涡模拟对化学氧碘激光器内的亚跨超音速混合过程进行模拟分析,其结果表明了大涡模拟对这种低压、低密度、亚跨超音速及夹杂多种介质的化学流场的可执行性。与传统的雷诺平均仿真结果相比较,大涡模拟能掌握更多的流场细节数据,能够对混合过程进行精准地判断和分析。在此基础上,提出了碘流反向45°入流的设计方案以增强混合程度,计算表明采用此种方案在相应出光面上平均小信号增益系数提高了5%。     

Experimental study and large eddy simulation of effect of terrain slope on marginal burning in shrub fuel beds

This paper presents a combined study of laboratory scale fire spread experiments and a three-dimensional large eddy simulation (LES) to analyze the effect of terrain slope on marginal burning behavior in live chaparral shrub fuel beds. Line fire was initiated in single species fuel beds of four common chaparral plants under various fuel bed configurations and ambient conditions. An LES approach was developed to model fire spreading through a fuel bed with a subgrid scale turbulent combustion model based on a flame surface density concept. By examining two fuel bed slope configurations, it was found that upslope fire spread depends not only on the increased radiant heat transfer but also on the aerodynamic effect created by the interaction of the flame with the inclined surface. Under certain conditions, the convective heat transfer induced by this interaction becomes the dominant mechanism in determining fire spread success. Seventy-three (or 42%) of 173 experimental fires successfully propagated for slopes ranging from −70% to 70%. It was found there exists a critical slope above which fire spread in these live fuel beds was successful, and below which fire spread was unsuccessful. This critical slope for marginal burning varied widely with fuel moisture content and fuel loading. A stepwise logistic regression model was developed from experimental data to predict the probability of successful fire spread. It is expected that this model may be helpful in providing guidelines for prescribed fire application.