QUANTITATIVE OH PLANAR LASER INDUCED FLUORESCENCE DIAGNOSTICS OF SYNGAS AND METHANE COMBUSTION IN A CAVITY COMBUSTOR

The current work reports quantitative OH species concentration in the cavity of a trapped vortex combustor (TVC) in the context of mixing and flame stabilization studies using both syngas and methane fuels. Planar laser induced fluorescence (PLIF) measurements of OH radical obtained using a Nd: YAG pumped dye laser are quantified using a flat flame McKenna burner. The momentum flux ratio (MFR), defined as the ratio of the cavity fuel jet momentum to that of the guide vane air stream, is observed to be a key governing parameter. At high MFRs similar to 4.5, the flame front is observed to form at the interface of the fuel jet and the air jet stream. This is substantiated by velocity vector field measurements. For syngas, as the MFR is lowered to similar to 0.3, the fuel-air mixing increases and a flame front is formed at the bottom and downstream edge of the cavity where a stratified charge is present. This trend is observed for different velocities at similar equivalence ratios. In case of methane combustion in the cavity, where the MFRs employed are extremely low at similar to 0.01, a different mechanism is observed. A fuel-rich mixture is now observed at the center of the cavity and this mixture undergoes combustion. On further increase of the cavity equivalence ratio, the rich mixture exceeds the flammability limit and forms a thin reaction zone at the interface with air stream. As a consequence, a shear layer flame at the top of the cavity interface with the mainstream is also observed. The equivalence ratio in the cavity also determines the combustion characteristics in the case of fuel-air mixtures that are formed as a result of the mixing. Overall, flame stabilization mechanisms have been proposed, which account for the wide range of MFRs and premixing in the mainstream as well.

Effect of Channel Confinement on Mixed Convective Flow Past an Equilateral Triangular Cylinder

The present work investigates the mixed convective flow and heat transfer characteristics past a triangular cylinder placed symmetrically in a vertical channel. At a representative Reynolds number, Re = 100, simulations are carried out for the blockage ratios beta = 1/3; 1/4; and 1/6. Effect of aiding and opposing buoyancy is brought about by varying the Richardson number in the range -1.0 <= Ri <= 1.0. At a blockage ratio of 1/3, suppression of vortex shedding is found at Ri = 1, whereas von Karman vortex street is seen both at beta = 1/4 and 1/6, respectively. This is the first time that such behavior of blockage ratio past a triangular cylinder in the present flow configuration is reported. Drag coefficient increases progressively with increasing Ri and a slightly higher value is noticed at beta = 1/3. For all b, heat transfer increases with increasing Ri. Flattening of Nu(avg)-Ri curve beyond Ri > 0: 75 is observed at beta = 1/3.

Intertwined vorticity and elastodynamics in flapping wing propulsion

We performed numerical experiments on a one-dimensional elastic solid oscillating in a two-dimensional viscous incompressible fluid with the intent of discerning the interplay of vorticity and elastodynamics in flapping wing propulsion. Perhaps for the first time, we have established the role of foil deflection topology and its influence on vorticity generation, through spatially and temporally evolving foil slope and curvature. Though the frequency of oscillation of the foil has a definite role, it is the phase relation between foil slope and pressure that determines thrust or drag. Similarly, the phase difference between flapping velocity, and pressure and inertial forces, determine the power input to the foil, and in turn drives propulsive efficiency. At low frequencies of oscillation, the sympathetic slope and curvature of deformation of the foil allow generation of leading-edge vortices that do not separate; they cause substantial rise in pressure between the leading edge and mid-chord. The circulatory component of pressure is determined primarily by the leading-edge vortex and therefore thrust too is predominantly circulatory in origin at low frequencies. In the intermediate and high-frequency range, thrust and drag on the foil spatially alternate and non-circulatory forces dominate over circulatory and viscous forces. For the mass ratios we simulated, thrust due to flapping varies quadratically as a function of Strouhal number or trailing-edge flapping velocity; further, the trailing edge flapping velocities peak at the same set of frequencies where the thrust is also a maximum. Propulsive efficiency, on the other hand, is roughly a mirror image of the thrust variation with respect to Strouhal number. Given that most instances of flapping propulsion in nature are primarily through distributed muscular actuation that enables precise control of deformation shape, leading to high thrust and efficiency, the results presented here are pointers towards understanding some of the mechanisms that drive thrust and propulsive efficiency.

Instabilities in unsteady boundary layers with reverse flow

Instabilities arising in unsteady boundary layers with reverse flow have been investigated experimentally. Experiments are conducted in a piston driven unsteady water tunnel with a shallow angle diffuser placed in the test section. The ratio of temporal (Pi(t)) to spatial (Pi(x)) component of the pressure gradient can be varied by a controlled motion of the piston. In all the experiments, the piston velocity variation with time is trapezoidal consisting of three phases: constant acceleration from rest, constant velocity and constant deceleration to rest. The adverse pressure gradient (and reverse flow) are due to a combination of spatial deceleration of the free stream in the diffuser and temporal deceleration of the free stream caused by the piston deceleration. The instability is usually initiated with the formation of one or more vortices. The onset of reverse flow in the boundary layer, location and time of formation of the first vortex and the subsequent flow evolution are studied for various values of the ratio Pi(x) (Pi(x) + Pi(t)) for the bottom and the top walls. Instability is due to the inflectional velocity profiles of the unsteady boundary layer. The instability is localized and spreads to the other regions at later times. At higher Reynolds numbers growth rate of instability is higher and localized transition to turbulence is observed. Scalings have been proposed for initial vortex formation time and wavelength of the instability vortices. Initial vortex formation time scales with convective time, delta/Delta U, where S is the boundary layer thickness and Delta U is the difference of maximum and minimum velocities in the boundary layer. Non-dimensional vortex formation time based on convective time scale for the bottom and the top walls are found to be 23 and 30 respectively. Wavelength of instability vortices scales with the time averaged boundary layer thickness. (C) 2015 Elsevier Masson SAS. All rights reserved.

Doping dependence of fluctuation diamagnetism in high T-c superconductors

Using a recently proposed Ginzburg-Landau-like lattice free energy functional due to Banerjee et al. (2011) we calculate the fluctuation diamagnetism of high -T-c superconductors as a function of doping, magnetic field and temperature. We analyse the pairing fluctuations above the superconducting transition temperature in the cuprates, ranging from the strong phase fluctuation dominated underdoped limit to the more conventional amplitude fluctuation dominated overdoped regime. We show that a model where the pairing scale increases and the superfluid density decreases with underdoping produces features of the observed magnetization in the pseudogap region, in good qualitative and reasonable quantitative agreement with the experimental data. In particular, we explicitly show that even when the pseudogap has a pairing origin the magnetization actually tracks the superconducting dome instead of the pseudogap temperature, as seen in experiment. We discuss the doping dependence of the `onset' temperature for fluctuation diamagnetism and comment on the role of vortex core -energy jn our model. (C) 2015 Elsevier Inc. All rights reserved.

Transitions and blowoff of unconfined non-premixed swirling flame

The present experimental work reports the first observations of primary and secondary transitions in the time-averaged flame topology in a non-premixed swirling flame as the geometric swirl number S-G (a non dimensional number used to quantify the intensity of imparted swirl) is varied from a magnitude of zero till flame blowout. First observations of two transition types viz. primary and secondary transitions are reported. The primary transition represents a transformation from yellow straight jet flame (at S-G = 0) to lifted flame with blue base and finally to swirling seated (burner attached) yellow flame. Time-averaged streamline plot obtained from 2D PIV in mid-longitudinal plane shows a recirculation zone (RZ) at the immediate vicinity of burner exit. The lifted flame is stabilized along the vortex core of this RZ. Further, when S-G similar to 1.4-3, the first occurrence of vortex breakdown (VB) induced internal recirculation zone (IRZ) is witnessed. The flame now stabilizes at the upstream stagnation point of the VB-IRZ, which is attached to the burner lip. The secondary transition represents a transformation from a swirling seated flame to swirling flame with a conical tailpiece and finally to a highly-swirled near blowout oxidizer-rich flame. This transition is understood to be the result of transition in vortex breakdown modes of the swirling flow field from dual-ring VB bubble to central toroidal recirculation zone (CTRZ). The physics of transition is described on the basis of modified Rossby number (Ro(m)). Finally, when the swirl intensity is very high i.e. SG similar to 10, the flame blows out due to excessive straining and due to entrainment of large amount of oxidizer due to partial premixing. The present investigation involving changes in flame topology is immensely important because any change in global flame structure causes oscillatory heat release that can couple with dynamic pressure and velocity fluctuations leading to unsteady combustion. In this light, understanding mechanisms of flame stabilization is essential to tackle the problem of thermo-acoustic instability. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.