Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Friction of sliding surfaces carrying adsorbed lubricant layers

In concentrated contacts the behaviour of lubricants is much modified by the high local pressures: changes can arise both from molecular ordering within the very thin film lubricant layers present at the interface as well as from the deposition on the component surfaces of more solid-like polymeric boundary layers. These 'third bodies' separating the solid surfaces may have rheological or mechanical properties very different from those observed in the bulk. Classical elasto-hydrodynamic theory considers the entrapped lubricant to exhibit a piezo-viscous behaviour while the conventional picture of more solid boundary lubricant layers views their shear strength r as being linearly dependent on local pressure p, so that T = TO + ap where TO and a are constants. If TO is relatively small, then the coefficient of friction \i = T Ip ~ a and so Amonton's laws are recovered. However, the properties of adsorbed or deposited surface films, or indeed other third bodies such as debris layers, may be more complex than this. A preliminary study has looked quantitatively at the influence of the pressure dependence of the shear strength of any surface layer on the overall friction coefficient of a contact which is made up of an array of asperities whose height varies in a Gaussian manner. Individual contact points may be elastic or plastic. The analysis results in plots of coefficient of friction versus the service or load parameter PIH&NRa where P is the nominal pressure on the contact, HS the hardness of the deforming surface, N the asperity density, R the mean radius of curvature of the asperities, and a is the standard deviation of their height distribution. In principle, any variation oft withp can be incorporated into the model; however, in this initial study we have used data on colloidal suspensions from the group at the Ecole Centrale de Lyon as well as examining the effect of functional relationships of somewhat greater complexity than a simple linear form. Results of the analysis indicate that variations in fj. are possible as the load is varied which depend on the statistical spread of behaviour at individual asperity contacts. The value of this analysis is that it attempts to combine the behaviour of films on the molecular scale with the topography of real engineering surfaces and so give an indication of the effects at the full-size or macro-scale that can be achieved by chemical or molecular surface engineering.

Investigation of extinction in unsteady flames in turbulent combustion by 2D-LIF of OH radials and flamelet analysis

An experimental and theoretical investigation of premixed turbulent combustion in an engine simulator is presented. The distribution of hydroxyl radicals formed in the combustion of propane/air mixtures was visualized by 2D-LIF and used to monitor the progress of the combustion process. For stoichiometric mixtures, images showed a continuous wrinkled flame front, while in lean (λ=1.5) mixtures, local flame extinction was observed as discontinuities in the reaction zone. A bright active reaction zone was still observed in flame inlets and closed concave structures. The effects of self-absorption and of collisional quenching on the fluorescence signal are considered and appear to have only a minor net influence on the shape and width of the flame front. The images are evaluated and interpreted in terms of the Lewis number effect and the laminar flamelet model. Analysis was performed by determining the contour lines of the images (specifically, the ratios of average maximum to equilibrium OH concentration) and comparing with corresponding ratios from unstrained flame simulations. The results show that although the degree of turbulence is not high enough for straining effects to be important, flamelet curvature does play a significant role in the combustion of lean mixtures; this is manifested by a mean effective flame velocity that is less than the laminar burning velocity. © 1991 Combustion Institute.

Characterisation and simulation of asphalt road surfaces

Accurate simulation of rolling-tyre vibrations, and the associated noise, requires knowledge of road-surface topology. Full scans of the surface types in common use are, however, not widely available, and are likely to remain so. Ways of producing simulated surfaces from incomplete starting information are thus needed. In this paper, a simulation methodology based solely on line measurements is developed, and validated against a full two-dimensional height map of a real asphalt surface. First the tribological characteristics-asperity height, curvature and nearest-neighbour distributions-of the real surface are analysed. It is then shown that a standard simulation technique, which matches the (isotropic) spectrum and the probability distribution of the height measurements, is unable to reproduce these characteristics satisfactorily. A modification, whereby the inherent granularity of the surface is enforced at the initialisation stage, is introduced, and found to produce simulations whose tribological characteristics are in excellent agreement with the measurements. This method will thus make high-fidelity tyre-vibration calculations feasible for researchers with access to line-scan data only. In addition, the approach to surface tribological characterisation set out here provides a template for efficient cataloguing of road textures, as long as the resulting information can subsequently be used to produce sample realisations. A third simulation algorithm, which successfully addresses this requirement, is therefore also presented. © 2011 Elsevier B.V.

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.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Response of fishes to the underwater D.C. field

Response to external electric field (D. C.) of three different varieties of fish namely Puntius ticto, Heteropneustis fossilis and Tilapia mossambica having different anatomical and behavioural characteristics were studied. Clearly distinguished reactions occurred one after another m all the varieties of fish with the increase in field intensity with minor specific variations. These reactions can be broadly classified into initial start (first reaction), forced swimming (galvanotaxis), slackening of body muscle (galvanonarcosis) and state of muscular rigidity (tetanus). The orientation of the organism (projection of nervous element) to the surrounding field has been found to have important bearing on the behaviour reactions. Clearly differentiated anodic taxis and true narcosis set in when fish body axis was parallel to the lines of current conduction. But with greater angle between the body axis and the current lines, fish did not show well marked reactions. Fish body, when perpendicular to current lines responded for anodic curvature and off balance swimming followed by tetanus.

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.

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

COMPRESSOR BLADE

A blade (2) for a compressor comprising a pressure surface (6) and a suction surface (4), wherein the suction surface (4) comprises: a discontinuity in the chord-wise curvature of an intermediate portion of the suction surface between the blade leading edge and blade trailing edge; and a slot (20) arranged along at least a portion of the blade (2) in a substantially span-wise direction, the slot (2) being disposed at or near the said change in curvature, such that when in use lower momentum fluid near the suction surface (4) is aspirated into the slot (20).

Cinematographic OH-PLIF measurements of two interacting turbulent premixed flames with and without acoustic forcing

This paper describes an experimental investigation into the interactions that occur between two lean turbulent premixed flames stabilised on conical bluff-bodies when they are moved closer together. Cinematographic OH-PLIF measurements were acquired to investigate adjacent flame front interactions as a function of flame separation distance (S). Flame surface density (FSD) and curvature were determined to characterise the unforced flames. Acoustic forcing was then applied to explore the amplitude dependent thermo-acoustic response. Phase-averaged FSD and global heat release measurements in the form of OH * chemiluminescence were obtained for a range of forcing frequencies (f) and amplitudes (A) as a function of S. As the flames were brought closer together the adjacent annular jets were found to merge into a single jet structure. This caused adjacent flame fronts to merge above the wake region between the two flames at a location determined by the jet efflux (flame angle) and S. This region of flame-flame interaction we refer to as 'interacting region'. In the unforced flames, a trend of increasingly negative curvature for decreasing S produced a small net increase in flame surface area via cusp formation. When subjected to acoustic forcing, S-dependent regimes were found in the global heat release response as a function A. The overall trend showed that the occurrence of jet/flame merging reduces the value of A at which non-linear response occurs. In support of previous findings for flames stabilised along shear layers, the phase-averaged FSD showed that the flame dynamics that drive the thermo-acoustic response result from the roll-up of vortices which generate large-scale vortex-flame interactions. Compared with axisymmetric flames, the occurrence of jet merging alters the vortex-flame interactions resulting in an asymmetric contribution to the heat release between the wall and interacting regions. The majority of the heat release was found to occur in the interacting region through the rapid production and destruction of flame surface area. The occurrence of jet merging and large-scale interactions between adjacent flames result in different physical mechanisms that drive the thermo-acoustic response compared with single axisymmetric flames. © 2011.

The deposition of small particles from a turbulent air flow in a curved duct

The paper describes an experimental and theoretical study of the deposition of small spherical particles from a turbulent air flow in a curved duct. The objective was to investigate the interaction between the streamline curvature of the primary flow and the turbulent deposition mechanisms of diffusion and turbophoresis. The experiments were conducted with particles of uranine (used as a fluorescent tracer) produced by an aerosol generator. The particles were entrained in an air flow which passed vertically downwards through a long straight channel of rectangular cross-section leading to a 90° bend. The inside surfaces of the channel and bend were covered with tape to collect the deposited particles. Following a test run the tape was removed in sections, the uranine was dissolved in sodium hydroxide solution and the deposition rates established by measuring the uranine concentration with a luminescence spectrometer. The experimental results were compared with calculations of particle deposition in a curved duct using a computer program that solved the ensemble-averaged particle mass and momentum conservation equations. A particle density-weighted averaging procedure was used and the equations were expressed in terms of the particle convective, rather than total, velocity. This approach provided a simpler formulation of the particle turbulence correlations generated by the averaging process. The computer program was used to investigate the distance required to achieve a fully-developed particle flow in the straight entry channel as well as the variation of the deposition rate around the bend. The simulations showed good agreement with the experimental results. © 2012 Elsevier Ltd.

The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows

A series of flames in a turbulent methane/air stratified swirl burner is presented. The degree of stratification and swirl are systematically varied to generate a matrix of experimental conditions, allowing their separate and combined effects to be investigated. Non-swirling flows are considered in the present paper, and the effects of swirl are considered in a companion paper (Part II). A mean equivalence ratio of φ=0.75 is used, with φ for the highest level of stratification spanning 0.375-1.125. The burner features a central bluff-body to aid flame stabilization, and the influence of the induced recirculation zone is also considered. The current work focuses on non-swirling flows where two-component particle image velocimetry (PIV) measurements are sufficient to characterize the main features of the flow field. Scalar data obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution allow the behavior of key combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. Simultaneous cross-planar OH-PLIF is used to determine the orientation of the instantaneous flame normal in the scalar measurement window, allowing gradients in temperature and progress variable to be angle corrected to their three dimensional values. The relationship between curvature and flame thickness is investigated using the OH-PLIF images, as well as the effect of stratification on curvature.The main findings are that the behavior of the key combustion species in temperature space is well captured on the mean by laminar flame calculations regardless of the level of stratification. H 2 and CO are significant exceptions, both appearing at elevated levels in the stratified flames. Values for surface density function and by extension thermal scalar dissipation rate are found to be substantially lower than laminar values, as the thickening of the flame due to turbulence dominates the effect of increased strain. These findings hold for both premixed and stratified flames. The current series of flames is proposed as an interesting if challenging set of test cases for existing and emerging turbulent flame models, and data are available on request. © 2012 The Combustion Institute.

The structure of turbulent stratified and premixed methane/air flames II: Swirling flows

Experimental results are presented from a series of turbulent methane/air stratified flames stabilized on a swirl burner. Nine operating conditions are considered, systematically varying the level of stratification and swirl while maintaining a lean global mean equivalence ratio of φ̄=0.75. Scalar data are obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution, allowing the behavior of the major combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. The corresponding three-dimensional surface density function and thermal scalar dissipation rate are investigated, along with geometric characteristics of the flame such as curvature and flame thickness. Hydrogen and carbon monoxide levels within the flame brush are raised by stratification, indicating models with laminar premixed flame chemistry may not be suitable for stratified flames. However, flame surface density, scalar dissipation and curvature all appear insensitive to the degree of stratification in the flames surveyed. © 2012 The Combustion Institute.

Efficient truncation scheme for modeling deepwater mooring lines

This paper is aimed at enabling the confident use of existing model test facilities for ultra deepwater application without having to compromise on the widely accepted range of scales currently used by the floating production industry. Passive line truncation has traditionally been the preferred method of creating an equivalent numerical model at reduced depth; however, these techniques tend to suffer in capturing accurately line dynamic response and so reproducing peak tensions. In an attempt to improve credibility of model test data the proposed truncation procedure sets up the truncated model, based on line dynamic response rather than quasi-static system stiffness. The upper sections of each line are modeled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model that aims to simulate the remainder of the line. Stages 1 & 2 are used to derive a water depth truncation ratio. Here vibration decay of transverse elastic waves is assessed and it is found that below a certain length criterion, the transverse vibrational characteristics for each line are inertia driven, hence with respect to these motions the truncated model can assume a linear damper whose coefficient depends on the local line properties and vibration frequency. Stage 3 endeavors to match the individual line stiffness between the full depth and truncated models. In deepwater it is likely that taut polyester moorings will be used which are predominantly straight and have high axial stiffness that provides the principal restoring force to static and low frequency vessel motions. Consequently, it means that the natural frequencies of axial vibrations are above the typical wave frequency range allowing for a quasi-static solution. In cases of exceptionally large wave frequency vessel motions, localized curvature at the chain seabed segment and tangential skin drag on the polyester rope can increase dynamic peak tensions considerably. The focus of this paper is to develop an efficient scheme based on analytic formulation, for replicating these forces at the truncation. The paper will close with an example case study of a single mooring under extreme conditions that replicates exactly the static and dynamic characteristics of the full depth line. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).