Experimental and Numerical Studies in a Compact Trapped Vortex Combustor: Stability Assessment and Augmentation

Fundamental studies on a compact trapped vortex combustor indicate that cavity injection strategies play a major role on flame stability. Detailed experiments indicate that blow-out occurs for a certain range of cavity air flow velocities. An unsteady RANS-based reacting flow simulation tool has been utilized to study the basic dynamics of cavity vortex for various flow conditions. The phenomenon of flame blow-out at certain intermediate cavity air velocities is explained on the basis of transition from a cavity-stabilized mode to an opposed flow stagnation mode. A novel strategy is proposed for achieving flame stability at all conditions. This involves using a flow guide vane in the path of the main flow to direct a portion of the main flow into the cavity. This seems to result in a desirable dual vortex structure, i.e., a small clockwise vortex behind the vane and large counterclockwise vortex in the cavity. Experimental results show stable flame at all flow conditions with the flow guide vane, and pressure drop is estimated to be within acceptable limits. Cold flow simulations show self-similar velocity profiles for a range of main inlet velocities, and high reverse velocity ratios (-0.3) are observed. Such a high-velocity ratio in the reverse flow shear layer profile leads to enhanced production of turbulence imperative to compact combustors. Reacting flow simulations show even higher reverse velocity ratios (above -0.7) due to flow acceleration. The flame is observed to be stable, even though minor shear layer oscillations are present in the form of vortex shedding. Self-similarity is also observed in reacting flow temperature profiles at combustor exit over the entire range of the mainstream velocity. This indicates that the present configuration holds a promise of delivering robust performance invariant of the flow operating conditions.

Analysis of Lined Ducts With Mean Flow, With Application to Dissipative Mufflers

Mufflers with at least one acoustically absorptive duct are generally called dissipative mufflers. Generally, for want of systems approach, these mufflers are characterized by transmission loss of the lined duct with overriding corrections for the terminations, mean flow, etc. In this article, it is proposed that dissipative duct should be integrated with other muffler elements, source impedance and radiation impedance, by means of transfer matrix approach. Towards this end, the transfer matrix for rectangular duct with mean flow has been derived here, for the least attenuated mode. Mean flow introduces a coupling between transverse wave numbers and axial wave number, the evaluation of which therefore calls for simultaneous solution of two or three transcendental equations. This is done by means of a Newton-Raphson iteration scheme, which is illustrated here for square ducts lined with porous ceramic tiles.

Transfer matrix modeling of hyperbolic and parabolic ducts with incompressible mean flow

A general differential equation for the propagation of sound in a variable area duct or nozzle carrying incompressible mean flow (of low Mach number) is derived and solved for hyperbolic and parabolic shapes. Expressions for the state variables of acoustic pressure and acoustic mass velocity of the shapes are derived. Self‐consistent expressions for the four‐pole parameters are developed. The conical, exponential, catenoidal, sine, and cosine ducts are shown to be special cases of hyperbolic ducts. Finally, it is shown that if the mean flow in computing the transmission loss of the mufflers involving hyperbolic and parabolic shapes was not neglected, little practical benefit would be derived.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 2. Comparison of fluctuating force simulations and experiments

The particle and fluid velocity fluctuations in a turbulent gas-particle suspension are studied experimentally using two-dimensional particle image velocimetry with the objective of comparing the experiments with the predictions of fluctuating force simulations. Since the fluctuating force simulations employ force distributions which do not incorporate the modification of fluid turbulence due to the particles, it is of importance to quantify the turbulence modification in the experiments. For experiments carried out at a low volume fraction of 9.15 x 10(-5) (mass loading is 0.19), where the viscous relaxation time is small compared with the time between collisions, it is found that the gas-phase turbulence is not significantly modified by the presence of particles. Owing to this, quantitative agreement is obtained between the results of experiments and fluctuating force simulations for the mean velocity and the root mean square of the fluctuating velocity, provided that the polydispersity in the particle size is incorporated in the simulations. This is because the polydispersity results in a variation in the terminal velocity of the particles which could induce collisions and generate fluctuations; this mechanism is absent if all of the particles are of equal size. It is found that there is some variation in the particle mean velocity very close to the wall depending on the wall-collision model used in the simulations, and agreement with experiments is obtained only when the tangential wall-particle coefficient of restitution is 0.7. The mean particle velocity is in quantitative agreement for locations more than 10 wall units from the wall of the channel. However, there are systematic differences between the simulations and theory for the particle concentrations, possibly due to inadequate control over the particle feeding at the entrance. The particle velocity distributions are compared both at the centre of the channel and near the wall, and the shape of the distribution function near the wall obtained in experiments is accurately predicted by the simulations. At the centre, there is some discrepancy between simulations and experiment for the distribution of the fluctuating velocity in the flow direction, where the simulations predict a bi-modal distribution whereas only a single maximum is observed in the experiments, although both distributions are skewed towards negative fluctuating velocities. At a much higher particle mass loading of 1.7, where the time between collisions is smaller than the viscous relaxation time, there is a significant increase in the turbulent velocity fluctuations by similar to 1-2 orders of magnitude. Therefore, it becomes necessary to incorporate the modified fluid-phase intensity in the fluctuating force simulation; with this modification, the mean and mean-square fluctuating velocities are within 20-30% of the experimental values.

Flow-acoustic Characterisation of a Cavity-based Combustor Configuration

This study concerns the flow-acoustic characterisation of a cavity-based combustor configuration. A well-validated numerical tool has been used to simulate the unsteady, two-dimensional reacting flow. Initially, a conventional flow over a cavity with dimensions and conditions corresponding to a compact cavity combustor was studied. Cavity mass injections in the form of fuel and air injections required for trapped vortex formation were then employed and the resonance features of this configuration were studied. The results indicate that the cavity depth mode resonance mechanism is dominant at the conditions studied in this work and that the oscillation frequencies do not change with cavity air injection. This observation is important since it implies that the only important variable which can alter resonant frequencies is the cavity depth. With combustion, the pressure oscillation amplitude was observed to increases significantly due to periodic entrainment of the cavity air jet and fluctuation of fuel-air mixture composition to produce highly fluctuating heat-release rates. The underlying mechanisms of the unsteady flow in the cavity combustor identified in this study indicate the strong dependence of the acoustics on the cavity injection strategies.

Indexing of decagonal quasicrystals I. The T8-phase

Indexing of a decagonal quasicrystal using the scheme utilizing five planar vectors and one perpendicular to them is examined in detail. A method for determining the indices of zone axes that a reciprocal vector would make in a decagonal phase of any periodicity has been proposed. By this method, the location of the zone axes made by any reciprocal vector can be predicted. The orthogonality condition has been simplified for the zone axes containing twofold vectors. The locations of zone axes have also been determined by an alternative method, utilizing spherical trigonometric calculations, which confirm the zone-axis locations given by the indices. The effect of one-dimensional periodicity on the indices and the accuracy of the zone-axis determination is discussed. Rules for the formation of zone axes between several reciprocal vectors and the prediction of all the reciprocal vectors in a zone are evolved.

3D Acoustic Analysis Of Spherical Chamber Having Single Inlet And Multiple Outlet: An Impedance Matrix Approach

The transmission loss (TL) performance of spherical chambers having single inlet and multiple outlet is obtained analytically through modal expansion of acoustic field inside the spherical cavity in terms of the spherical Bessel functions and Legendre polynomials. The uniform piston driven model based upon the impedance [Z] matrix is used to characterize the multi-port spherical chamber. It is shown analytically that the [Z] parameters are independent of the azimuthal angle (phi) due to the axisymmetric shape of the sphere; rather, they depend only upon the polar angle (theta) and radius of the chamber R(0). Thus, the effects of relative polar angular location of the ports and number of outlet ports are investigated. The analytical results are shown to be in good agreement with the 3D FEA results, thereby validating the procedure suggested in this work.

Effects of acoustic-streaming-induced flow in evaporating nanofluid droplets

We study the effect of acoustic streaming on nanoparticle motion and morphological evolution inside an acoustically levitated droplet using an analytical approach coupled with experiments. Nanoparticle migration due to internal recirculation forms a density stratification, the location of which depends on initial particle concentration. The time scale of density stratification is similar to that of perikinetic-driven agglomeration of particle flocculation. The density stratification ultimately leads to force imbalance leading to a unique bowl-shaped structure. Our analysis shows the mechanism of bowl formation and how it is affected by particle size, concentration, internal recirculation and fluid viscosity.

Surface pressure and heat transfer measurements in the unsteady separated hypersonic flow field over double cones

Reliable bench mark experimental database in the separated hypersonic flow regime is necessary to validate high resolution CFD codes. In this paper we report the surface pressure and heat transfer measurements carried out on double cones (first cone semi-apex angle = 15, 25 deg.; second cone semi-apex angle= 35, 68 deg.) at hypersonic speeds that will be useful for CFD code validation studies. The surface pressure measurements are carried out at nominal Mach number of 8.35 in the IISc hypersonic wind tunnel. On the other hand the surface heat transfer measurements are carried out at a nominal Mach number of 5.75 in the IISc hypersonic shock tunnel. The flow separation point on the first cone, flow reattachment on the second cone and the wild fluctuation of the transmitted shock on the second cone surface (25/68 deg. double cone) in the presence of severe adverse pressure gradient are some of the flow features captured in the measurements. The results from the CFD studies indicate good agreement with experiments in the attached flow regime while considerable differences are noticeable in the separated flow regime.

Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Effect of capillarity and heterogeneity in the numerical modelling of multiphase flow of fluids in unsaturated porous medium

The multiphase flow of fluids in the unsaturated porous medium is considered as a three phase flow of water, NAPL, and air simultaneously in the porous medium. The adaptive solution fully implicit modified sequential method is used for the numerical modelling. The effect of capillarity and heterogeneity effect at the interface between the media is studied and it is observed that the interface criteria has to be taken into account for the correct prediction of NAPL migration especially in heterogeneous media. The modified Newton Raphson method is used for the linearization and Hestines and Steifel Conjugate Gradient method is used as the solver.

Enhancing Speedup in Network Processing Applications by Exploiting Instruction Reuse with Flow Aggregation

Instruction reuse is a microarchitectural technique that improves the execution time of a program by removing redundant computations at run-time. Although this is the job of an optimizing compiler, they do not succeed many a time due to limited knowledge of run-time data. In this paper we examine instruction reuse of integer ALU and load instructions in network processing applications. Specifically, this paper attempts to answer the following questions: (1) How much of instruction reuse is inherent in network processing applications?, (2) Can reuse be improved by reducing interference in the reuse buffer?, (3) What characteristics of network applications can be exploited to improve reuse?, and (4) What is the effect of reuse on resource contention and memory accesses? We propose an aggregation scheme that combines the high-level concept of network traffic i.e. "flows" with a low level microarchitectural feature of programs i.e. repetition of instructions and data along with an architecture that exploits temporal locality in incoming packet data to improve reuse. We find that for the benchmarks considered, 1% to 50% of instructions are reused while the speedup achieved varies between 1% and 24%. As a side effect, instruction reuse reduces memory traffic and can therefore be considered as a scheme for low power.

An instantaneous control model for flow control in an ATM network

This paper is concerned with the optimal flow control of an ATM switching element in a broadband-integrated services digital network. We model the switching element as a stochastic fluid flow system with a finite buffer, a constant output rate server, and a Gaussian process to characterize the input, which is a heterogeneous set of traffic sources. The fluid level should be maintained between two levels namely b1 and b2 with b1

Flow from T-spout ladles

The ladle constitutes a crucial element in the pouring system for developing process consistency and quality in the manufacture of castings. Flow of molten metal from T-spout ladles, with the spout projecting at an angle from the shell of the ladle, was investigated by simulation, modeling and analysis. This was followed by experimental validation for water flow in translucent models, and verification in a cast iron foundry. Key parameters in the design of a ladle (for its geometry) have been identified as the ladle size, spout diameter and its angle. Velocity and flow rate from a T-spout ladle are governed by and increase with increase in angle of tilt apart from the aforementioned design parameters. Cross section and profile of the issuing jet displaying a stable twisting pattern, during its free-fall of model fluid water and molten metal iron, are notable aspects in the design of the ladle, particularly the spout exit.