CFD investigation of turbulent premixed flame response to transverse forcing

Screech is a high frequency oscillation that is usually characterized by instabilities caused by large-scale coherent flow structures in the wake of bluff-body flameholders and shear layers. Such oscillations can lead to changes in flame surface area which can cause the flame to burn unsteadily, but also couple with the acoustic modes and inherent fluid-mechanical instabilities that are present in the system. In this study, the flame response to hydrodynamic oscillations is analyzed in a controlled manner using high-fidelity Computational Fluid Dynamics (CFD) with an unsteady Reynolds-averaged Navier-Stokes approach. The response of a premixed flame with and without transverse velocity forcing is analyzed. When unforced, the flame is shown to exhibit a self-excitation that is attributed to the anti-symmetric shedding of vortices in the wake of the flameholder. The flame is also forced using two different kinds of low-amplitude out-of-phase inlet velocity forcing signals. The first forcing method is harmonic forcing with a single characteristic frequency, while the second forcing method involves a broadband forcing signal with frequencies in the range of 500 - 1000 Hz. For the harmonic forcing method, the flame is perturbed only lightly about its mean position and exhibits a limit cycle oscillation that is characteristic of the forcing frequency. For the broadband forcing method, larger changes in the flame surface area and detachment of the flame sheet can be seen. Transition to a complicated trajectory in the phase space is observed. When analyzed systematically with system identification methods, the CFD results, expressed in the form of the Flame Transfer Function (FTF) are capable of elucidating the flame response to the imposed perturbation. The FTF also serves to identify, both spatially and temporally, regions where the flame responds linearly and nonlinearly. Locking-in between the flame's natural self-excited frequency and the subharmonic frequencies of the broadband forcing signal is found to alter the dynamical behaviour of the flame. Copyright © 2013 by ASME.

On the application of LES to seal geometries

Five Large Eddy Simulation (LES) and hybrid RANS-NLES (Reynolds-Averaged Navier-Stokes-Numerical-LES) methods are used to simulate flow through a labyrinth seal geometry and are contrasted with RANS solutions. Results show that LES and RANS-NLES is capable of accurately predicting flow behaviour of two seal flows with a scatter of less than 5 %. RANS solutions show the potential to perform poorly for the turbulence models tested. LES and hybrid RANS-NLES are found to be consistent and in agreement with measurements, providing a flexible numerical platform for design investigations. It also allows greater flow physics insights. © Springer Science+Business Media Dordrecht 2013.

Prediction of combustion noise for an aeroengine combustor

Combustion noise may become an important noise source for lean-burn gas turbine engines, and this noise is usually associated with highly unsteady flames. This work aims to compute the broadband combustion noise spectrum for a realistic aeroengine combustor and to compare with available measured noise data on a demonstrator aeroengine. A low-order linear network model is applied to a demonstrator engine combustor to obtain the transfer function that relates to unsteadiness in the rate of heat release, acoustic, entropic, and vortical fluctuations. A spectral model is used for the heat release rate fluctuation, which is the source of the noise. The mean flow of the aeroengine combustor required as input data to this spectral model is obtained from Reynolds-averaged Navier-Stokes simulations. The computed acoustic field for a low-to-medium power setting indicates that the models used in this study capture the main characteristics of the broadband spectral shape of combustion noise. Reasonable agreement with the measured spectral level is achieved. © 2012 AIAA.

Global linear stability of the boundary-layer flow over a rotating sphere

We consider the linear global stability of the boundary-layer flow over a rotating sphere. Our results suggest that a self-excited linear global mode can exist when the sphere rotates sufficiently fast, with properties fixed by the flow at latitudes between approximately 55°-65° from the pole (depending on the rotation rate). A neutral curve for global linear instabilities is presented with critical Reynolds number consistent with existing experimentally measured values for the appearance of turbulence. The existence of an unstable linear global mode is in contrast to the literature on the rotating disk, where it is expected that nonlinearity is required to prompt the transition to turbulence. Despite both being susceptible to local absolute instabilities, we conclude that the transition mechanism for the rotating-sphere flow may be different to that for the rotating disk. © 2014 Elsevier Masson SAS. All rights reserved.

A flame surface density approach to large-eddy simulation of premixed turbulent combustion

The flame surface density approach to the modeling of premixed turbulent combustion is well established in the context of Reynolds-averaged simulations. For the future, it is necessary to consider large-eddy simulation (LES), which is likely to offer major advantages in terms of physical accuracy, particularly for unsteady combustion problems. LES relies on spatial filtering for the removal of unresolved phenomena whose characteristic length scales are smaller than the computational grid scale. Thus, there is a need for soundly based physical modeling at the subgrid scales. The aim of this paper is to explore the usefulness of the flame surface density concept as a basis for LES modeling of premixed turbulent combustion. A transport equation for the filtered flame surface density is presented, and models are proposed for unclosed terms. Comparison with Reynolds-averaged modeling is shown to reveal some interesting similarities and differences. These were exploited together with known physics and statistical results from experiment and from direct numerical stimulation in order to gain insight and refine the modeling. The model has been implemented in a combustion LES code together with standard models for scalar and momentum transport. Computational results were obtained for a simple three-dimensional flame propagation test problem, and the relative importance of contributing terms in the modeled equation for flame surface density was assessed. Straining and curvature are shown to have a major influence at both the resolved and subgrid levels.

Driving forces shaping phytoplankton assemblages in two subtropical plateau lakes with contrasting trophic status

1. We studied driving forces shaping phytoplankton assemblages in two subtropical plateau lakes with contrasting trophic status, the oligotrophic deep Lake Fuxian and the eutrophic shallow Lake Xingyun. 2. Phytoplankton samples were taken monthly for a year and phytoplankton species were sorted into the main taxonomic groups and associations proposed by Reynolds. A canonical correspondence analysis (CCA) was used to test the occurrence of these classification schemes and to determine their discriminatory power. 3. The results suggest that the major driving forces in Lake Fuxian were physical variables, and particularly the underwater light climate, whereas, nutrients were the important driving force in Lake Xingyun. 4. Top-down control through zooplankton grazing in Lake Fuxian was hardly ever a significant determinant itself, because of the scarcity of zooplankton and their low grazing efficiency of predation while a dominance of inedible cyanobacteria throughout the year rendered top-down controls ineffective failing in Lake Xingyun. Hence phytoplankton communities in both lakes appear to be regulated primarily by bottom-up controls.

Measurements of intake separation hysteresis in a model fan and nacelle RIG

A 1/20-scale, low speed model rig representing the fan and nacelle of a high bypass ratio jet engine has been tested under crosswind conditions. The flow conditions under which the intake flow separates and reattaches have been found to exhibit considerable hysteresis. This phenomenon has been examined by a careful test procedure in which the crosswind angle has been slowly increased and then decreased. Measurements of the hysteresis associated with separation and reattachment are presented for independent variations in stream-tube contraction ratio, ground clearance, fan operating point and Reynolds number. The results reveal that particular care must be taken to allow for any hysteresis when testing intakes under crosswind conditions. They also indicate that separation hysteresis is particularly sensitive to fan operating point and the position of the ground plane. These findings suggest that it is important for high Reynolds number intake tests and calculations to include a ground plane and a model of the downstream turbomachinery. © 2002 by the author(s).

Evidence of print gap airflow affecting web printing quality

The current study extends our earlier investigation on the real-time dynamics of print gap airflow around a single jetted drop over a moving substrate. In the present work, simulated web press printing was performed using a stationary grey-scale commercial inkjet print head to print full-width block of solid colour images onto a paper substrate with extended print gaps. The resultant printed images exhibit patterns or 'wood-graining' effects which become more prevalent as the relevant Reynolds number (Re) increases. The high-resolution scans of the printed images revealed that the patterns are created by oscillation and coalescence of neighboring printed tracks across the web. The phenomenon could be a result of drop stream perturbations caused by unsteady print gap airflow of the type similar to that observed in the previous study. ©2013; Society for Imaging Science and Technology.

On the laminar-turbulent transition of the rotating-disk flow: The role of absolute instability

© 2014 Cambridge University Press. This paper describes a detailed experimental study using hot-wire anemometry of the laminar-turbulent transition region of a rotating-disk boundary-layer flow without any imposed excitation of the boundary layer. The measured data are separated into stationary and unsteady disturbance fields in order to elaborate on the roles that the stationary and the travelling modes have in the transition process. We show the onset of nonlinearity consistently at Reynolds numbers, R, of ∼ 510, i.e. at the onset of Lingwood's (J. Fluid Mech., vol. 299, 1995, pp. 17-33) local absolute instability, and the growth of stationary vortices saturates at a Reynolds number of ∼ 550. The nonlinear saturation and subsequent turbulent breakdown of individual stationary vortices independently of their amplitudes, which vary azimuthally, seem to be determined by well-defined Reynolds numbers. We identify unstable travelling disturbances in our power spectra, which continue to grow, saturating at around R=585, whereupon turbulent breakdown of the boundary layer ensues. The nonlinear saturation amplitude of the total disturbance field is approximately constant for all considered cases, i.e. different rotation rates and edge Reynolds numbers. We also identify a travelling secondary instability. Our results suggest that it is the travelling disturbances that are fundamentally important to the transition to turbulence for a clean disk, rather than the stationary vortices. Here, the results appear to show a primary nonlinear steep-fronted (travelling) global mode at the boundary between the local convectively and absolutely unstable regions, which develops nonlinearly interacting with the stationary vortices and which saturates and is unstable to a secondary instability. This leads to a rapid transition to turbulence outward of the primary front from approximately R=565 to 590 and to a fully turbulent boundary layer above 650.

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.

Variational calculations of ionized-donor-bound excitons in GaAs-AlxGa1-xAs quantum wells

Using a two-parameter wave function, we calculate variationally the binding energy of an exciton bound to an ionized donor impurity (D+,X) in GaAs-AlxGa1-xAs quantum wells for the values of the well width from 10 to 300 Angstrom, when the dopant is located in the center of the well and at the edge of the well. The theoretical results confirm that the previous experimental speculation proposed by Reynolds tit al. [Phys. Rev. B 40, 6210 (1989)] is the binding energy of D+,X for the dopant at the edge of the well. in addition, we also calculate the center-of-mass wave function of the exciton and the average interparticle distances. The results are discussed in detail.

Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Study on methane hydration process in a semi-continuous stirred tank reactor

The methane hydration process is investigated in a semi-continuous stirred tank reactor. Liquid temperatures and reaction rates without stirrer are compared with those occurring with stirrer, while at the same time better stirring conditions of the methane hydration process are given by the experiments. Some basic data of fluid mechanics, for example, stirring Reynolds number, Froucle number and stirrer power, are calculated during the methane hydration process, which can be applied to evaluate stirrer capacity and provide some basic data for a scaled up reactor. Based on experiment and calculations in this work, some conclusions are drawn. First, the stirrer has great influence on the methane hydration process. Batch stirring is helpful to improve the mass transfer and heat transfer performances of the methane hydration process. Second, induction time can be shortened effectively by use of the stirrer. Third, in this paper, the appropriate stirring velocity and stirring time were 320 rpm and 30 min, respectively, at 5.0 MPa, for which the storage capacity and reaction time were 159.1 V/V and 370 min, respectively. Under the condition of the on-flow state, the initial stirring Reynolds number of the fluid and the stirring power were 12,150 and 0.54 W, respectively. Fourth, some suggestions, for example, the use of another type of stirrer or some baffles, are proposed to accelerate the methane hydration process. Comparing with literature data, higher storage capacity and hydration rate are achieved in this work. Moreover, some fluid mechanics parameters are calculated, which can provide some references to engineering application.

草本植物对土质路面径流水动力学特征及水沙过程的影响

通过室内人工模拟降雨,研究了草本植物不同覆盖度土质路面的径流水动力学特征以及路面的产流产沙过程。结果表明,路面径流的弗劳德数和雷诺数属层流的急流范畴,并且二者与径流流速及过水断面单位能量相同,均随草本植物覆盖度的增大而减小。而Darcy-weisbach阻力系数以及曼宁糙率系数则均随覆盖度的增大而增大。在不同覆盖度土质路面上,径流量随径流时间急剧增大后趋于稳定,且二者符合移轴双曲线关系,而输沙率则经历了先增大后减小的过程。随着植被覆盖度的增大,土质道路的水分入渗率逐渐增大,而径流量、输沙率及含沙率均逐渐减小。该研究表明土质路面种植的草本植物具有较好的水土保持功能,为阐明草本植物与土壤侵蚀的关系提供了理论依据。

薄层水流水力特性试验研究

对粘砂床面的薄层水流进行定床试验,研究了不同流量、坡度下薄层水流的水力学参数(流态、水深、流速及阻力系数)随坡度和流量的变化规律。结果表明:水流雷诺数Re随着坡度和流量的增大而增大,但Re主要是由流量决定的,坡度对Re的影响不大,水流流态基本上在过渡流区和紊流区;平均水深-单宽流量,流速-单宽流量,水流阻力系数-单宽流量均成很好的幂函数形式,平均水深和流速均随单宽流量的增大而增大,阻力系数随着单宽流量的增大而减小。