The effect of clearance on shrouded and unshrouded turbines at two levels of reaction

In this paper, the effect of seal clearance on the efficiency of a turbine with a shrouded rotor is compared with the effect of the tip clearance when the same turbine has an unshrouded rotor. The shrouded versus unshrouded comparison was undertaken for two turbine stage designs one having 50% reaction the other having 24% reaction. Measurements for a range of clearances, including very small clearances, showed three important phenomena. Firstly, as the clearance is reduced, there is a "break-even clearance" at which both the shrouded turbine and the unshrouded turbine have the same efficiency. If the clearance is reduced further, the unshrouded turbine performs better than the shrouded turbine, with the difference at zero clearance termed the "offset loss". This is contrary to the traditional assumption that both shrouded and unshrouded turbines have the same efficiency at zero clearance. The physics of the break-even clearance and the offset loss are discussed. Secondly, the use of a lower reaction had the effect of reducing the tip leakage efficiency penalty for both the shrouded and the unshrouded turbines. In order to understand the effect of reaction on the tip leakage, an analytical model was used and it was found that the tip leakage efficiency penalty should be understood as the dissipated kinetic energy rather than either the tip leakage mass flow rate or the tip leakage loss coefficient. Thirdly, it was also observed that, at a fixed flow coefficient, the fractional change in the output power with clearance was approximately twice the fractional change in efficiency with clearance. This was explained by using an analytical model. © 2010 by ASME.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.

The effect of clearance on shrouded and unshrouded turbines at two levels of reaction

In this paper, the effect of seal clearance on the efficiency of a turbine with a shrouded rotor is compared with the effect of the tip clearance when the same turbine has an unshrouded rotor. The shrouded versus unshrouded comparison was undertaken for two turbine stage designs one having 50% reaction, the other having 24% reaction. Measurements for a range of clearances, including very small clearances, showed three important phenomena. Firstly, as the clearance is reduced, there is a "break-even clearance" at which both the shrouded turbine and the unshrouded turbine have the same efficiency. If the clearance is reduced further, the unshrouded turbine performs better than the shrouded turbine, with the difference at zero clearance termed the "offset loss." This is contrary to the traditional assumption that both shrouded and unshrouded turbines have the same efficiency at zero clearance. The physics of the break-even clearance and the offset loss are discussed. Secondly, the use of a lower reaction had the effect of reducing the tip leakage efficiency penalty for both the shrouded and the unshrouded turbines. In order to understand the effect of reaction on the tip leakage, an analytical model was used and it was found that the tip leakage efficiency penalty should be understood as the dissipated kinetic energy rather than either the tip leakage mass flow rate or the tip leakage loss coefficient. Thirdly, it was also observed that, at a fixed flow coefficient, the fractional change in the output power with clearance was approximately twice the fractional change in efficiency with clearance. This was explained by using an analytical model. © 2014 by ASME.

Spatial correlation of heat release rate and sound emission from turbulent premixed flames

The two-point spatial correlation of the rate of change of fluctuating heat release rate is central to the sound emission from open turbulent flames, and a few attempts have been made to address this correlation in recent studies. In this paper, the two-point correlation and its role in combustion noise are studied by analysing direct numerical simulation (DNS) data of statistically multi-dimensional turbulent premixed flames. The results suggest that this correlation function depends on the separation distance and direction but, not on the positions inside the flame brush. This correlation can be modelled using a combination of Hermite-Gaussian functions of zero and second order, i.e. functions of the form (1-Ax2)e-Bx2 for constants A and B, to include its possible negative values. The integral correlation volume obtained using this model is about 0.2δL3 with the length scale obtained from its cube root being about 0.6δ L, where δ L is the laminar flame thermal thickness. Both of the values are slightly larger than the values reported in an earlier study because of the anisotropy observed for the correlation. This model together with the turbulence-dependent parameter K, the ratio of the root-mean-square (RMS) value of the rate of change of reaction rate to the mean reaction rate, derived from the DNS data is applied to predict the far-field sound emitted from open flames. The calculated noise levels agree well with recently reported measurements and show a sensitivity to K values. © 2012 The Combustion Institute.

Effects of Lewis number on scalar variance transport in premixed flames

The influences of differential diffusion rates of heat and mass on the transport of the variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been studied using three-dimensional simplified chemistry based Direct Numerical Simulation (DNS) data of statistically planar turbulent premixed flames with global Lewis number ranging from Le = 0.34 to 1.2. The Lewis number effects on the statistical behaviours of the various terms of the transport equations of variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been analysed in the context of Reynolds Averaged Navier Stokes (RANS) simulations. It has been found that the turbulent fluxes of the progress variable and temperature variances exhibit counter-gradient transport for the flames with Lewis number significantly smaller than unity whereas the extent of this counter-gradient transport is found to decrease with increasing Lewis number. The Lewis number is also shown to have significant influences on the magnitudes of the chemical reaction and scalar dissipation rate contributions to the scalar variance transport. The modelling of the unclosed terms in the scalar variance equations for the non-unity Lewis number flames have been discussed in detail. The performances of the existing models for the unclosed terms are assessed based on a-priori analysis of DNS data. Based on the present analysis, new models for the unclosed terms of the active scalar variance transport equations are proposed, whenever necessary, which are shown to satisfactorily capture the behaviours of unclosed terms for all the flames considered in this study. © 2010 Springer Science+Business Media B.V.

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding.

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.

Kinetics of the chemical looping oxidation of H2 by a co-precipitated mixture of CuO and Al2O3

Chemical-looping combustion (CLC) has the inherent property of separating CO2 from flue gases. Instead of air, it uses an oxygen-carrier, usually in the form of a metal oxide, to provide oxygen for combustion. When used for the combustion of gaseous fuels, such as natural gas, or synthesis gas from the gasification of coal, the technique gives a stream of CO2 which, on an industrial scale, would be sufficiently pure for geological sequestration. An important issue is the form of the metal oxide, since it must retain its reactivity through many cycles of complete reduction and oxidation. Here, we report on the rates of oxidation of one constituent of synthesis gas, H2, by co-precipitated mixtures of CuO+Al2O3 using a laboratory-scale fluidised bed. To minimise the influence of external mass transfer, and also of errors in the measurement of [H2], particles sized to 355-500μm were used at low [H2], with the temperature ranging from 450 to 900°C. Under such conditions, the reaction was slow enough for meaningful measurements of the intrinsic kinetics to be made. The reaction was found to be first order with respect to H2. Above ∼800°C, the reaction of CuO was fast and conformed to the shrinking core mechanism, proceeding via the intermediate, Cu2O, in: 2CuO+H2→Cu2O+H2O, ΔH1073 K0=- 116.8 kJ/mol; Cu2O+H2→2Cu+H2O, ΔH1073 K0-80.9 kJ/mol. After oxidation of the products Cu and Cu2O back to CuO, the kinetics in subsequent cycles of chemical looping oxidation of H2 could be approximated by those in the first. Interestingly, the carrier was found to react at temperatures as low as 300°C. The influence of the number of cycles of reduction and oxidation is explored. Comparisons are drawn with previous work using reduction by CO. Finally, these results indicate that the kinetics of reaction of the oxygen carrier with gasifier synthesis gases is very much faster than rates of gasification of the original fuel. © 2010 The Institution of Chemical Engineers.

An investigation of the kinetics of CO2 uptake by a synthetic calcium based sorbent

This study examines the kinetics of carbonation by CO2 at temperatures of ca. 750 °C of a synthetic sorbent composed of 15 wt% mayenite (Ca12Al14O33) and CaO, designated HA-85-850, and draws comparisons with the carbonation of a calcined limestone. In-situ XRD has verified the inertness of mayenite, which neither interacts with the active CaO nor does it significantly alter the CaO carbonation–calcination equilibrium. An overlapping grain model was developed to predict the rate and extent of carbonation of HA-85-850 and limestone. In the model, the initial microstructure of the sorbent was defined by a discretised grain size distribution, assuming spherical grains. The initial input to the model – the size distribution of grains – was a fitted parameter, which was in good agreement with measurements made with mercury porosimetry and by the analysis of SEM images of sectioned particles. It was found that the randomly overlapping spherical grain assumption offered great simplicity to the model, despite its approximation to the actual porous structure within a particle. The model was able to predict the performance of the materials well and, particularly, was able to account for changes in rate and extent of reaction as the structure evolved after various numbers of cycles of calcination and carbonation.

An investigation of the kinetics of CO 2 uptake by a synthetic calcium based sorbent

This study examines the kinetics of carbonation by CO 2 at temperatures of ca. 750°C of a synthetic sorbent composed of 15wt% mayenite (Ca 12Al 14O 33) and CaO, designated HA-85-850, and draws comparisons with the carbonation of a calcined limestone. In-situ XRD has verified the inertness of mayenite, which neither interacts with the active CaO nor does it significantly alter the CaO carbonation-calcination equilibrium. An overlapping grain model was developed to predict the rate and extent of carbonation of HA-85-850 and limestone. In the model, the initial microstructure of the sorbent was defined by a discretised grain size distribution, assuming spherical grains. The initial input to the model - the size distribution of grains - was a fitted parameter, which was in good agreement with measurements made with mercury porosimetry and by the analysis of SEM images of sectioned particles. It was found that the randomly overlapping spherical grain assumption offered great simplicity to the model, despite its approximation to the actual porous structure within a particle. The model was able to predict the performance of the materials well and, particularly, was able to account for changes in rate and extent of reaction as the structure evolved after various numbers of cycles of calcination and carbonation. © 2011 Elsevier Ltd.

Increased gas production from hydrates by combining depressurization with heating of the wellbore

To extract gas from hydrate reservoirs, it has to be dissociated in situ. This endothermic dissociation process absorbs heat energy from the formation and pore fluid. The heat transfer governs the dissociation rate, which is proportional to the difference between the actual temperature and the equilibrium temperature. This study compares three potential gas production schemes from hydrate-bearing soil, where the radial heat transfer is governing. Cylindrical samples with 40% pore-filling hydrate saturation were tested. The production tests were carried out over 90 min by dissociating the hydrate from a centered miniature wellbore, by either lowering the pressure to 6, 4, or 6 MPa with simultaneous heating of the wellbore to 288 K. All tests were replicated by a numerical simulation. With additional heating at the same wellbore pressure, the gas production from hydrates could, on average, be increased by 1.8 and 3.6 times in the simulation and experiments, respectively. If the heat influx from the outer boundary is limited, a simulation showed that the specific heat of the formation is rapidly used up when the wellbore is only depressurized and not heated. © 2012 American Chemical Society.

The effect of addition of ZrO2 to Fe2O3 for hydrogen production by chemical looping

In this paper, a synthetic mixture of ZrO2 and Fe 2O3 was prepared by coprecipitation for use in chemical looping and hydrogen production. Cycling experiments in a fluidized bed showed that a material composed of 30 mol % ZrO2 and 70 mol % Fe 2O3 was capable of producing hydrogen with a consistent yield of 90 mol % of the stoichiometric amount over 20 cycles of reduction and oxidation at 1123 K. Here, the iron oxide was subjected to cycles consisting of nearly 100% reduction to Fe followed by reoxidation (with steam or CO 2 and then air) to Fe2O3. There was no contamination by CO of the hydrogen produced, at a lower detection limit of 500 ppm, when the conversion of Fe3O4 to Fe was kept below 90 mol %. A preliminary investigation of the reaction kinetics confirmed that the ZrO2 support does not inhibit rates of reaction compared with those observed with iron oxide alone. © 2012 American Chemical Society.

Application of axial discontinuity factors for high conversion water reactors

The paper shows that generating cross sections using three-dimensional geometry and application of axial discontinuity factors are essential requirements for obtaining accurate prediction of criticality and zone average reaction rates in highly heterogeneous RBWR-type systems using computer codes based on diffusion theory approximation. The same methodology as presented here will be used to generate discontinuity factors for each axial interface between fuel assembly zones to ensure preservation of reaction rates in each zone and global multiplication factor. The use of discontinuity factors and three-dimensional cross sections may allow for a coarser energy group structure which is desirable to simplify and speed up transient calculations.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.

Reynolds number effects on scalar dissipation rate transport and its modeling in turbulent premixed combustion

The effects of turbulent Reynolds number, Ret, on the transport of scalar dissipation rate of reaction progress variable in the context of Reynolds averaged Navier-Stokes simulations have been analyzed using three-dimensional simplified chemistry-based direct numerical simulation (DNS) data of freely propagating turbulent premixed flames with different values of Ret. Scaling arguments have been used to explain the effects of Ret on the turbulent transport, scalar-turbulence interaction, and the combined reaction and molecular dissipation terms. Suitable modifications to the models for these terms have been proposed to account for Ret effects, and the model parameters include explicit Ret dependence. These expressions approach expected asymptotic limits for large values of Ret. However, turbulent Reynolds number Ret does not seem to have any major effects on the modeling of the term arising from density variation. Copyright © Taylor and Francis Group, LLC.

Scalar gradient behaviour in MILD combustion

The results of three-dimensional Direct Numerical Simulation (DNS) of Moderate, Intense Low-oxygen Dilution (MILD) and conventional premixed turbulent combustion conducted using a skeletal mechanism including the effects of non-unity Lewis numbers and temperature dependent transport properties are analysed to investigate combustion characteristics using scalar gradient information. The DNS data is also used to synthesise laser induced fluorescence (LIF) signals of OH, CH2O, and CHO. These signals are analysed to verify if they can be used to study turbulent MILD combustion and it has been observed that at least two (OH and CH2O) LIF signals are required since the OH increase across the reaction zone is smaller in MILD combustion compared to premixed combustion. The scalar gradient PDFs conditioned on the reaction rate obtained from the DNS data and synthesised LIF signals suggests a strong gradient in the direction normal to the MILD reaction zone with moderate reaction rate implying flamelet combustion. However, the PDF of the normal gradient is as broad as for the tangential gradient when the reaction rate is high. This suggests a non-flamelet behaviour, which is due to interaction of reaction zones. The analysis of the conditional PDFs for the premixed case confirms the expected behaviour of scalar gradient in flamelet combustion. It has been shown that the LIF signals synthesised using 2D slices of DNS data also provide very similar insights. These results demonstrate that the so-called flameless combustion is not an idealised homogeneous reactive mixture but has common features of conventional combustion while containing distinctive characteristics. © 2013 The Combustion Institute.