Comment and reply on the 'assessment of combustion and submodels for turbulent nonpremixed hydrocarbon flames'

In deriving the flamelet model for nonpremixed combustion certain terms, but not the unsteady term, are assumed to be negligible. This results in a relation between all reacting scalars and the mixture fraction as independent variable. An ideal test of the flamelet assumption can be based on direct numerical simulation (DNS) data, if all reacting scalars are conditioned on mixture fraction and conditional moments are evaluated. The fundamental assumption of the flamelet model are unwillingly justified. The unsteady and steady formulations of the same equations are compared and found that unsteadiness is important in an unsteady simulation.

Analytic model for turn off in the silicon-on-insulator LIGBT

Lateral insulated gate bipolar transistors (LIGBTs) in silicon-on-insulator (SOI) show a unique turn off characteristic when compared to junction-isolated RESURF LIGBTs or vertical IGBTs. The turn off characteristic shows an extended `terrace' where, after the initial fast transient characteristic of IGBTs due to the loss of the electron current, the current stays almost at the same value for an extended period of time, before suddenly dropping to zero. In this paper, we show that this terrace arises because there is a value of LIGBT current during switch off where the rate of expansion of the depletion region with respect to the anode current is infinite. Once this level of anode current is approached, the depletion region starts to expand very rapidly, and is only stopped when it reaches the n-type buffer layer surrounding the anode. Once this happens, the current rapidly drops to zero. A quasi-static analytic model is derived to explain this behaviour. The analytically modelled turn off characteristic agrees well with that found by numerical simulation.

Vortex Transport and Blade Interactions in High Pressure Turbines

This paper presents a study of the three-dimensional flow field within the blade rows of a high-pressure axial flow steam turbine stage. Half-delta wings were fixed to a rotating hub to simulate an upstream rotor passage vortex. The flow field is investigated in a Low-Speed Research Turbine using pneumatic and hot-wire probes downstream of the blade row. The paper examines the impact of the delta wing vortex transport on the performance of the downstream blade row. Steady and unsteady numerical simulations were performed using structured 3D Navier-Stokes solver to further understand the flow field. The loss measurements at the exit of the stator blade showed an increase in stagnation pressure loss due to the delta wing vortex transport. The increase in loss was 21% of the datum stator loss, demonstrating the importance of this vortex interaction. The transport of the stator viscous flow through the rotor blade row is also described. The rotor exit flow was affected by the interaction between the enhanced stator passage vortex and the rotor blade row. Flow underturning near the hub and overturning towards the mid-span was observed, contrary to the classical model of overturning near the hub and underturning towards the mid-span. The unsteady numerical simulation results were further analysed to identify the entropy producing regions in the unsteady flow field.

The time-resolved flow within an unsteady ejector

The flow field within an unsteady ejector has been investigated using experimental and computational techniques. The experimental results show a peak thrust augmentation of 1.4; numerical simulation gives a value of 1.37. It is shown that the vortex ring dominates the flow field. At optimal thrust augmentation the vortex ring acts like a fluid piston accelerating the fluid inside the ejector. A model is proposed for the operation of unsteady ejectors, based on the vortex ring acting like a fluid piston. Control volume analysis is presented showing that mass entrainment is responsible for thrust augmentation. It is proposed that the spacing of successive vortex rings determines the mass entrainment and therefore thrust augmentation. The efficiency of unsteady ejectors was found to vary between 28% and 32% depending on the L/D ratio of the unsteady jet source. Copyright © 2008 by J H Heffer.

On computing the physical sources of jet noise

We derive a closed system of equations that relates the acoustically radiating flow variables to the sources of sound for homentropic flows. We use radiating density, momentum density and modified pressure as the dependent variables which leads to simple source terms for the momentum equations. The source terms involve the non-radiating parts of the density and momentum density fields. These non-radiating components are obtained by removing the radiating wavenumbers in the Fourier domain. We demonstrate the usefulness of this new technique on an axi-symmetric jet solution of the Navier-Stokes equations, obtained by direct numerical simulation (DNS). The dominant source term is proportional to the square of the non-radiating part of the axial momentum density. We compare the sound sources to that obtained by an acoustic analogy and find that they have more realistic physical properties. Their frequency content and amplitudes are consistent with. We validate the sources by computing the radiating sound field and comparing it to the DNS solution. © 2010 by S. Sinayoko, A. Agarwal.

RANS-URANS in axial compressor, a design methodology

Turbomachinery flows are inherently unsteady. Until now during the design process, unsteadiness has been neglected, with resort merely to steady numerical simulations. Despite the assumption involved, the results obtained with steady simulations have been used with success. One of the questions arising in recent years is can unsteady simulations be used to improve the design of turbomachines? In this work the numerical simulation of a multi-stage axial compressor is carried out. Comparison of Reynolds averaged Navier-Stokes (RANS) and unsteady Reynolds averaged Navier-Stokes (URANS) calculation shows that the unsteadiness affects pressure losses and the prediction of stall limit. The unsteady inflow due to the wake passing mainly modifies the losses and whirl angle near the endwalls. The computational cost of the fully unsteady compared with a steady simulation is about four times in terms of mesh dimension and two orders of magnitude as number of iterations. A mixed RANS-URANS solution has been proposed to give the designer the possibility to simulate an unsteady stage embedded in a steady-state simulation. This method has been applied to the simulation of a four-stage axial compressor rig. The mixed RANS-URANS approach has been developed using sliding and mixing planes as interface conditions. The rotor-stator interaction has been captured physically while reducing the computational time and mesh size.

Effects of soil fabric on behaviors of granular soils: Microscopic modeling

The effects of initial soil fabric on behaviors of granular soils are investigated by using Distinct Element Method (DEM) numerical simulation. Soil specimens are represented by an assembly of non-uniform sized spheres with different initial contact normal distributions. Isotropically consolidated triaxial compression loading and extension unloading in both undrained and drained conditions are simulated for vertically- and horizontally-sheared specimens. The numerical simulation results are compared qualitatively with the published experimental data and the effects of initial soil fabric on resulting soil behaviors are discussed, including the effects of specimen reconstitution methods, effects of large preshearing, and anisotropic characteristics in undrained and drained conditions. The effects of initial soil fabric and mode of shearing on the quasi-steady state line are also investigated. The numerical simulation results can systematically explain that the observed experimental behaviors of granular soils are due principally to their conditions of the initial soil fabric. This outcome provides insights into the observed phenomena in microscopic view. © 2011 Elsevier Ltd.

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.

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.

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.

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

Design and optimization of a 250nm SOI LDMOSFET

This work is aimed at optimising the static performance of a high voltage SOI LDMOSFET. Starting with a conventional LDMOSFET, 2D and 3D numerical simulation models, able to accurately match datasheet values, have been developed. Moving from the original device, several design techniques have been investigated with the target of improving the breakdown voltage and the ON-state resistance. The considered design techniques are based on the modification of the doping profile of the drift region and the Superjunction design technique. The paper shows that a single step doping within the drift region is the best design choice for the considered device and is found to give a 24% improvement in the breakdown voltage and a 17% reduction of the ON-state resistance. © 2011 IEEE.

Effect of flow-geometry on turbulence-scalar interaction in premixed flames

Turbulent combustion of stoichiometric hydrogen-air mixture is simulated using direct numerical simulation methodology, employing complex chemical kinetics. Two flame configurations, freely propagating and V-flames stabilized behind a hot rod, are simulated. The results are analyzed to study the influence of flame configuration on the turbulence-scalar interaction, which is critical for the scalar gradient generation processes. The result suggests that this interaction process is not influenced by the flame configuration and the flame normal is found to align with the most extensive strain in the region of intense heat release. The combustion in the rod stabilized flame is found to be flamelet like in an average sense and the growth of flame-brush thickness with the downstream distance is represented well by Taylor theory of turbulent diffusion, when the flame-brushes are non-interacting. The thickness is observed to saturate when the flame-brushes interact, which is found to occur in the simulated rod stabilized flame with Taylor micro-scale Reynolds number of 97. © 2011 American Institute of Physics.

Assessment of computational fluid dynamics and experimental data for shock boundary-layer interactions

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop, numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric, and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and, in general, it was difficult to discern clear trends in the data. For the Reynolds-averaged Navier-Stokes (RANS) methods, the choice of turbulence model appeared to be the largest factor in solution accuracy. Scale-resolving methods, such as large-eddy simulation (LES), hybrid RANS/LES, and direct numerical simulation, produced error levels similar to RANS methods but provided superior predictions of normal stresses. Copyright © 2012 by Daniella E. Raveh and Michael Iovnovich.

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.