Application of scalar dissipation rate modelling to industrial burners in partially premixed regimes

The objective of this paper is to test various available turbulent burning velocity models on an experimental version of Siemens small scale combustor using the commercial CFD code. Failure of burning velocity model with different expressions for turbulent burning velocity is observed with an unphysical flame flashback into the swirler. Eddy Dissipation Model/Finite Rate Chemistry is found to over-predict mean temperature and species concentrations. Solving for reaction progress equation with its variance using scalar dissipation rate modelling produced reasonably good agreement with the available experimental data. Two different turbulence models Shear Stress Transport (SST) and Scale Adaptive Simulation (SAS) SST are tested and results from transient SST simulations are observed to be predicting well. SAS-SST is found to under-predict with temperature and species distribution.

Application of scalar dissipation rate model to siemens DLE combustors

The standard design process for the Siemens Industrial Turbomachinery, Lincoln, Dry Low Emissions combustion systems has adopted the Eddy Dissipation Model with Finite Rate Chemistry for reacting computational fluid dynamics simulations. The major drawbacks of this model have been the over-prediction of temperature and lack of species data limiting the applicability of the model. A novel combustion model referred to as the Scalar Dissipation Rate Model has been developed recently based on a flamelet type assumption. Previous attempts to adopt the flamelet philosophy with alternative closure models have failed, with the prediction of unphysical phenomenon. The Scalar Dissipation Rate Model (SDRM) was developed from a physical understanding of scalar dissipation rate, signifying the rate of mixing of hot and cold fluids at scales relevant to sustain combustion, in flames and was validated using direct numerical simulations data and experimental measurements. This paper reports on the first industrial application of the SDRM to SITL DLE combustion system. Previous applications have considered ideally premixed laboratory scale flames. The industrial application differs significantly in the complexity of the geometry, unmixedness and operating pressures. The model was implemented into ANSYS-CFX using their inbuilt command language. Simulations were run transiently using Scale Adaptive Simulation turbulence model, which switches between Large Eddy Simulation and Unsteady Reynolds Averaged Navier Stokes using a blending function. The model was validated in a research SITL DLE combustion system prior to being applied to the actual industrial geometry at real operating conditions. This system consists of the SGT-100 burner with a glass square-sectioned combustor allowing for detailed diagnostics. This paper shows the successful validation of the SDRM against time averaged temperature and velocity within measurement errors. The successful validation allowed application of the SDRM to the SGT-100 twin shaft at the relevant full load conditions. Limited validation data was available due to the complexity of measurement in the real geometry. Comparison of surface temperatures and combustor exit temperature profiles showed an improvement compared to EDM/FRC model. Furthermore, no unphysical phenomena were predicted. This paper presents the successful application of the SDRM to the industrial combustion system. The model shows a marked improvement in the prediction of temperature over the EDM/FRC model previously used. This is of significant importance in the future applications of combustion CFD for understanding of hardware mechanical integrity, combustion emissions and dynamics of the flame. Copyright © 2012 by ASME.

Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

A novel theoretical approach to passive control of thermo-acoustic oscillations: Application to ducted heat sources

In this paper, we develop a linear technique that predicts how the stability of a thermo-acoustic system changes due to the action of a generic passive feedback device or a generic change in the base state. From this, one can calculate the passive device or base state change that most stabilizes the system. This theoretical framework, based on adjoint equations, is applied to two types of Rijke tube. The first contains an electrically-heated hot wire and the second contains a diffusion flame. Both heat sources are assumed to be compact so that the acoustic and heat release models can be decoupled. We find that the most effective passive control device is an adiabatic mesh placed at the downstream end of the Rijke tube. We also investigate the effects of a second hot wire and a local variation of the cross-sectional area but find that both affect the frequency more than the growth rate. This application of adjoint sensitivity analysis opens up new possibilities for the passive control of thermo-acoustic oscillations. For example, the influence of base state changes can be combined with other constraints, such as that the total heat release rate remains constant, in order to show how an unstable thermo-acoustic system should be changed in order to make it stable. Copyright © 2013 by ASME.

Application of fast oxygen sensors for investigations into air-path dynamics and EGR distribution in a diesel engine

The control of NOX emissions by exhaust gas recirculation (EGR) is of widespread application. However, despite dramatic improvements in all aspects of engine control, the subtle mixing processes that determine the cylinder-to-cylinder distribution of the recirculated gas often results in a mal-distribution that is still an issue for the engine designer and calibrator. In this paper we demonstrate the application of a relatively straightforward technique for the measurement of the absolute and relative dilution quantity in both steady state and transient operation. This was achieved by the use of oxygen sensors based on standard UEGO (universal exhaust gas oxygen) sensors but packaged so as to give good frequency response (∼ 10 ms time constant) and be completely insensitivity to the sample pressure and temperature. Measurements can be made at almost any location of interest, for example exhaust and inlet manifolds as well as EGR path(s), with virtually no flow disturbance. At the same time, the measurements yield insights into air-path dynamics. We argue that "dilution", as indicated by the deviation of the oxygen concentration from that of air, is a more appropriate parameter than EGR rate in the context of NOX control, especially for diesel engines. Experimental results are presented for the EGR distribution in a current production light duty 4-cylinder diesel engine in which significant differences were found in the proportion of the recirculated gas that reached each cylinder. Even the individual inlet runners of the cylinders exhibited very different dilution rates - differences of nearly 50% were observed at some conditions. An application of such data may be in the improvement of calibration and validation of CFD and other modelling techniques. Copyright © 2014 SAE International.