Modeling evaporation effects in conditional moment closure for spray autoignition

Simulations of an n-heptane spray autoigniting under conditions relevant to a diesel engine are performed using two-dimensional, first-order conditional moment closure (CMC) with full treatment of spray terms in the mixture fraction variance and CMC equations. The conditional evaporation term in the CMC equations is closed assuming interphase exchange to occur at the droplet saturation mixture fraction values only. Modeling of the unclosed terms in themixture fraction variance equation is done accordingly. Comparison with experimental data for a range of ambient oxygen concentrations shows that the ignition delay is overpredicted. The trend of increasing ignition delay with decreasing oxygen concentration, however, is correctly captured. Good agreement is found between the computed and measured flame lift-off height for all conditions investigated. Analysis of source terms in the CMC temperature equation reveals that a convective-reactive balance sets in at the flame base, with spatial diffusion terms being important, but not as important as in lifted jet flames in cold air. Inclusion of droplet terms in the governing equations is found to affect the mixture fraction variance field in the region where evaporation is the strongest, and to slightly increase the ignition delay time due to the cooling associated with the evaporation. Both flame propagation and stabilization mechanisms, however, remain unaffected. © 2011 Taylor & Francis.

Effects of turbulent reynolds number on the displacement speed statistics in the thin reaction zones regime of turbulent premixed combustion

The effects of turbulent Reynolds number on the statistical behaviour of the displacement speed have been studied using three-dimensional Direct Numerical Simulation of statistically planar turbulent premixed flames. The probability of finding negative values of the displacement speed is found to increase with increasing turbulent Reynolds number when the Damkhler number is held constant. It has been shown that the statistical behaviour of the Surface Density Function, and its strain rate and curvature dependence, plays a key role in determining the response of the different components of displacement speed. Increasing the turbulent Reynolds number is shown to reduce the strength of the correlations between tangential strain rate and dilatation rate with curvature, although the qualitative nature of the correlations remains unaffected. The dependence of displacement speed on strain rate and curvature is found to weaken with increasing turbulent Reynolds number when either Damkhler or Karlovitz number is held constant, but the qualitative nature of the correlation remains unaltered. The implications of turbulent Reynolds number effects in the context of Flame Surface Density (FSD) modelling have also been addressed, with emphasis on the influence of displacement speed on the curvature and propagation terms in the FSD balance equation. © 2011 Nilanjan Chakraborty et al.

Statistics and modelling of turbulent kinetic energy transport in different regimes of premixed combustion

The statistical behaviour of turbulent kinetic energy transport in turbulent premixed flames is analysed using data from three-dimensional Direct Numerical Simulation (DNS) of freely propagating turbulent premixed flames under decaying turbulence. For flames within the corrugated flamelets regime, it is observed that turbulent kinetic energy is generated within the flame brush. By contrast, for flames within the thin reaction zones regime it has been found that the turbulent kinetic energy decays monotonically through the flame brush. Similar trends are observed also for the dissipation rate of turbulent kinetic energy. Within the corrugated flamelets regime, it is demonstrated that the effects of the mean pressure gradient and pressure dilatation within the flame are sufficient to overcome the effects of viscous dissipation and are responsible for the observed augmentation of turbulent kinetic energy in the flame brush. In the thin reaction zones regime, the effects of the mean pressure gradient and pressure dilatation terms are relatively much weaker than those of viscous dissipation, resulting in a monotonic decay of turbulent kinetic energy across the flame brush. The modelling of the various unclosed terms of the turbulent kinetic energy transport equation has been analysed in detail. The predictions of existing models are compared with corresponding quantities extracted from DNS data. Based on this a-priori DNS assessment, either appropriate models are identified or new models are proposed where necessary. It is shown that the turbulent flux of turbulent kinetic energy exhibits counter-gradient (gradient) transport wherever the turbulent scalar flux is counter-gradient (gradient) in nature. A new model has been proposed for the turbulent flux of turbulent kinetic energy, and is found to capture the qualitative and quantitative behaviour obtained from DNS data for both the corrugated flamelets and thin reaction zones regimes without the need to adjust any of the model constants. © 2010 Springer Science+Business Media B.V.

Conditional Moment Closure for Turbulent Premixed Flames

The conditional moment closure (CMC) method has been successfully applied to various non-premixed combustion systems in the past, but its application to premixed flames is not fully tested and validated. The main difficulty is associated with the modeling of conditional scalar dissipation rate of the conditioning scalar, the progress variable. A simple algebraic model for the conditional dissipation rate is validated using DNS results of a V-flame. This model along with the standard k- turbulence modeling is used in computations of stoichiometric pilot stabilized Bunsen flames using the RANS-CMC method. A first-order closure is used for the conditional mean reaction rate. The computed non reacting and reacting scalars are in reasonable agreement with the experimental measurements and are consistent with earlier computations using flamelets and transported PDF methods. Sensitivity to chemical kinetic mechanism is also assessed. The results suggest that the CMC may be applied across the regimes of premixed combustion.

On the stabilisation of turbulent lifted flames

This study explores the stabilisation mechanisms of turbulent lifted flames by examining the scalar dissipation rate (SDR) of both passive and reactive scalars and their cross dissipation (CDR) in the stabilisation region. DNS results of a laboratory scale hydrogen turbulent lifted flame has been used for this analysis. Various definitions of the flame leading edge (FLE) has been compared and differences are illustrated. Time and spatial averaged statistic of SDR and CDR were examined. It was found that the averaged SDR for the mixture fraction at FLE was well below the reference quenching value for stoichiometric mixture. The averaged SDR for the progress variable is in the same order of the unstrained premixed laminar flame value. It was observed that the averaged CDR changed from negative to weakly positive at FLE. The change in sign was explained by a change in the relative alignment of the gradients of mixture fraction and progress variable. It was thus evident that the CDR was a good marker for stabilisation region and an important quantity in stabilisation mechanism.

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.

Simulation of spherical turbulent premixed flames

Outwardly-propagating spherical turbulent premixed flames are studied using Reynolds-Averaged Navier-Stokes methodology. The reaction rate closure is based on recently developed strained flamelets, where the flamelets are parametrised using the scalar dissipation rate. It is shown that the leading edge speed for the spherical flame is higher than for the planar case, for a given turbulence and thermo-chemical conditions. In addition, it is shown that including the mean curvature effects in the reaction rate closure do not influence the spherical flame speeds when compared with a model that excludes these effects.

Scalar dissipation rate transport in the corrugated flamelets and thin reaction zones regime for turbulent premixed flames

The characteristics of the scalar dissipation rate transport in the corrugated flamelets and the thin reaction zones regimes are studied based on two three-dimensional Direct Numerical Simulation (DNS) databases for freely propagating statistically planar turbulent premixed flames. The turbulent flame parameters are so chosen that the database which represents the corrugated flamelets regime has a global Damköhler number Da>1 whereas the database representing the thin reaction zones regime has Da <1. It is demonstrated that the terms originating from the correlation between fluctuating velocity and scalar gradient T1 shows strong Da dependence. The terms originating from dilatation T2, the scalar inner product of gradients of velocity and scalar fields T3 and the correlation between reaction rate and scalar gradients T4 and the dissipation term D2 remain important for both the flames. However, T3 dissipates scalar dissipation rate in the Da > 1 flame while it produces scalar dissipation rate in the Da < 1 flame. This difference is because of the change in the alignment between scalar and velocity gradients

The study of biological characteristics and heavy metal determination of Cd, Pb, Fe, V, Ni in the muscle, liver and gonad tissues of kingfish Scomberomorus commerson in the coastal waters of the Hormozagan Province (Persian Gulf)

Biological studies and heavy metal (Ni, v, Fe, Pb,Cd) determination in liver, gonad and muscle of Scomberomorus commerson were carried out from Oct 2006 to Sept 2007 in Hormozgan coastal waters. 599 Samples were gutted for reproduction and nutrition studies, fork length and weight were measured to nearest cm and g respectively. Meanwhile 40 samples were also investigated for heavy metal studies. All specimens were collected from two major landing sites (Bandar Lengeh & Bandar Abbas). Minimums & maximum fork length & weight were 29, 128 cm & 235 and 15350 g respectively.Isometric growth was shown according to our study and b was estimated 2.9 (overall), 2.91(male) & 2.89 (female). The average relative gut length was 0.52± 0.007 and it was determined that S. commerson is a carnivorous. More than 99 percent of gut content was different teleost fishes. Gastro somatic index had two peaks in Nov & Jan (before spawning) and with a decreased trend in July, the spawning period. Occurrence of empty stomach was estimated % 65.77. Maximum amount of condition factor was in Dec. Spawning season was started from June. The average of Absolute & relative fecundity (to weight unit) was 1217149±179315 and 178.2±15.58 respectively. Lm50% was estimated 75 cm for females. Sex ratio was 0.97: 1 (male: female). Chi- Square test showed no significant difference (p>0.05). Maximum amount of hepatosomatic index was estimated in March.Metal concentrations were determined using either Flame Atomic Absorption Spectroscopy (for Fe) or Graphite Furnace Atomic Absorption Spectroscopy (for Pb,Cd,Ni and V). The mean concentration (μg/g dry weight)of Pb,Cd,Ni,V and Fe in the liver were 0.0309, 0.0268, 0.0672, 0.0077, 2.5159 in the gonad 0.0440 ,0.0295, 0.1096, 0.0000, 1.4449 and in the muscle 0.0244, 0.0324, 0.0656, 0.0128, 1.6138 respectively. The maximum metal concentrations were below the maximum permissible limits for human consumption recommended by the USEPA, WHO and the UK. The results of Kendall's Tau-b correlation coefficient were as follows: The Liver tissue: There were significant positive linear relationships between accumulation of V, Fe, and Pb with Fork length, Pb and Fe with weight, GSI with Pb, Cd, V and 109 Fe, and a negative linear relationships between HSI with accumulation of V and Fe, Fork length, weight and GSI. The Gonad tissue: There were significant positive linear relationships between GSI with accumulation of Pb, Cd, Fe, Fork length and weight, a negative linear relationship between HSI with Fork length, weight and GSI. The Muscle tissue: There were significant positive linear relationships between accumulations of V, Fe with Fork length and weight factors and as well as GSI with Cd, V, Fe, Pb, Fork length and weight,a negative linear relationship between HSI with Fork length, weight, Cd, Fe and GSI. The results of Mann-Whitney U tests (P≤0.05) show that there were significant differences between summer and autumn from heavy metal contents in the studied tissues point of view. The only exceptions were for Ni in the liver, gonads and muscle and as well as there were significant differences between male and female from heavy metal contents in the studied tissues. The only exceptions were for Pb in muscle, Ni in liver, gonad and muscle, V in muscle, and Cd and Fe in gonads.

Investigation of self excited and acoustically forced combustion instabilities with laser diagnostics

Combustion oscillations in gas turbines can result in serious damage. One method used to predict such oscillations is to analyze the combustor acoustics using a simple linear model. Such a model requires a flame transfer function to describe the response of the heat release to flow perturbations inside the combustor. This paper reports on the application of Planar Laser Induced Fluorescence (PLIF) of OH radicals to analyze the response of a lean premixed flame to oncoming flow perturbations. Both self-excited oscillations and low amplitude forced oscillations at various frequencies are investigated in an atmospheric pressure model combustor rig. In order to visualize fluctuations of local fuel distribution, acetone-PLIF was also applied in non-reacting and acoustically forced flows at oscillation frequencies of 200 Hz and 510 Hz, respectively. OH-PLIF images were acquired over a range of operating parameters. The results presented in this paper originate from data sets acquired at fixed phase angles during the oscillation cycle. Comparative experiments in self excited and forced acoustic oscillations show that the flame and the combustion intensity develop similarly throughout the pressure cycle in both cases. Although the peak fluorescence intensities differ between self excited and the forced instabilities, there is a clear correspondence in the observed frequency and phase information from the two cases. This result encourages a comparison of the OH-PLIF and the acetone-PLIF results. Quantitative measurements of the equivalence ratio in specific areas of the measurement plane offer insight on the complex phenomena coupling acoustic perturbations, i.e. flow velocity fluctuations, to fluctuations in fuel distribution and combustion intensity, ultimately resulting in self excited combustion oscillations.

Study ectoparasite infestation of 10 common imported freshwater ornamental fish species

In the present study, ectoparasite infestation of common freshwater ornamental fish species imported into Iran including: Dwarf gourami (Colisa lalia), Angelfish (Pterophyllum scalare), Oscar (Astronatus ocellatus), Red eye Tetra ( Moenkhansia Sanctaefilomenae), Barb( Capoeta tetrazona), Arowana (Osteoglossum bicirrhosum), Goldfish(Carassius auratus), Red fin Shark (Lebeo elythrurus, Catfish (Hypostomus plecostomus) and Pangasius sutchi ( Pangasius hypophthalmus) from May until April was investigated.A total of 6 specimens sample of each fish species) randomly were obtained and taken alive to the lab. After observing gill and skin wet smear under the microscope , gill was dissected and examined carefully .Photo and films were taken from the isolated parasites and parasite identification was performed according to Yamaguti AP1), Bychowsky (■Y).From a total of examined fishes 1\ ► sample (YY.V.) were parasitized .Parasite groups that were observed consisting ciliated protozoan myxosporidian , monogenean Digenean metacercaria and crustacean .The only non-parasitized fish was Barb. Monogenea Trianchoratus sp. , Cleidodiscus sp. , Ancylodiscoides sp. , Thaparocleidus sp. , Centrocestus formosanus metacercaria and Myxobolus longisporus report for the first time in ornamental fish which imported from Southeast of Asia to Iran. According to the results, it seems that severe quarantine and sanitary rules are necessary.

Increasing the life of flare tips

Flare tips are essential for safety. Maintenance is difficult and costly. Flare tips are subjected to high combustion temperatures, thermal cycling, oxidation and marine corrosion. Following a number of flare tip failures an in depth study by Imperial College was carried out into the failure of a flare tip from a UK platform, looking for service life improvement. Materials selection and design solutions were considered. The study considered alternative materials and concluded that materials selection was the smaller part of the answer; design changes can double service life. This study used failure investigation, high temperature experimental and thermo-mechanical modelling analysis. The modelling process simulated two common flaring conditions and correctly predicted the observed failure of initiation and crack propagation from holes used to bolt on flame stabilizing plates to the top of the flare. The calculated thermal stress and strains enabled the low cycle fatigue life and minimum creep life to be predicted. It was concluded that service life could be improved by replacing Incoloy alloy 800HT (UNS N08800) with Inconel alloy 625 (UNS N06625), an alloy with attractive mechanical properties and improved high temperature corrosion resistance. Repositioning or eliminating bolt holes can double service life. Copyright 2008, Society of Petroleum Engineers.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.