CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

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.

Application of scalar dissipation rate modelling to industrial burners in partially premixed regimes

The objective of this paper is to test various available turbulent burning velocity models on an experimental version of Siemens small scale combustor using the commercial CFD code. Failure of burning velocity model with different expressions for turbulent burning velocity is observed with an unphysical flame flashback into the swirler. Eddy Dissipation Model/Finite Rate Chemistry is found to over-predict mean temperature and species concentrations. Solving for reaction progress equation with its variance using scalar dissipation rate modelling produced reasonably good agreement with the available experimental data. Two different turbulence models Shear Stress Transport (SST) and Scale Adaptive Simulation (SAS) SST are tested and results from transient SST simulations are observed to be predicting well. SAS-SST is found to under-predict with temperature and species distribution.

The propagation of acoustic waves in a slowly varying duct with radially sheared axial mean flow

We consider the propagation of acoustic waves along a cylindrical duct carrying radially sheared axial mean flow, in which the duct radius is allowed to vary slowly along the axis. In previous work [A.J. Cooper & N. Peake, Journal of Fluid Mechanics 445 (2001) 207-234.] radially sheared axial mean flow with nonzero swirl in a slowly varying duct was considered, but in this paper we set the swirl to zero, thereby allowing simplification of the calculations of both the mean and unsteady flows. In this approach the acoustic wavenumber and corresponding eigenfunction are determined locally, while the wave amplitude is found by solving an evolution equation along the duct. Sample results are presented, including a case in which, perhaps surprisingly, the number of cut-on modes increases as the duct radius decreases. © 2013 Elsevier Ltd. All rights reserved.

Laminar propagation of lean premixed flames ignited in stratified mixture

Combustion in stratified mixtures is envisaged in practical energy systems such as direct-injection spark-ignited (DISI) car engines, gas turbines, for reducing CO2 and pollutant emissions while protecting their efficiency. The mixture gradients change the fundamental properties of the flame, especially by a difference in temperature and composition between the burnt gases and those of a flame consuming a homogeneous mixture. This paper presents an investigation of the properties of the flame propagating in a lean homogeneous mixture after ignition in a richer mixture according to the magnitude of the stratification. Three magnitudes of stratification are investigated. The local flame burning velocity is determined by an original PIV algorithm developed previously. The local equivalence ratio in the fresh gases is measured from anisole PLIF. From the simultaneous PIV-PLIF measurements, the flame burning velocities conditioned on the local stretch rate and equivalence ratio in fresh gases are measured. The flame propagating through the homogeneous lean mixture has properties depending on the ignition conditions in the stratified layer. The flame propagating in the lean mixture is back-supported longer for ignition under the richer condition. The change of stretch sensitivity and burning velocity of the flame in the lean mixture is measured over time for the three magnitudes of mixture stratification investigated. The ignition in richer mixtures compensates for the nonequidiffusion effect of lean propane flame and sustains its robustness to stretch. The flame propagation in the lean homogeneous mixture is enhanced by ignition in a richer stratified layer, as much by their robustness to stretch as by an increase in the flame speed or the burning velocity. The decay time of this influence of the stratification, called memory effect, is determined. © 2013 The Combustion Institute.

Wave propagation in a slowly-varying duct with mean vortical swirling flow

The propagation of unsteady disturbances in a slowlyvarying cylindrical duct carrying mean swirling flow is investigated using a multiple-scales technique. This is applicable to turbomachinery flow behind a rotor stage when the swirl and axial velocities are of the same order. The presence of mean vorticity couples acoustic and vorticity equations which produces an eigenvalue problem that is not self-adjoint unlike that for irrotational mean flow. In order to determine the amplitude variation along the duct, an adjoint solution for the coupled system of equations is derived. The solution breaks down where a mode changes from cut on to cut off. In this region the amplitude is governed by a form of Airy's equation, and the effect of swirl is to introduce a small shift in the origin of the Airy function away from the turning-point location. The variation of axial wavenumber and amplitude along the duct is calculated. In hard-walled ducts mean swirl is shown to produce much larger amplitude variation along the duct compared with a nonswirling flow. Mean swirl also has a large effect in ducts with finite-impedance walls which differs depending on whether modes are co-rotating with the swirl or counter rotating. © 2001 by A.J. Cooper, Published by the American Institute of Aeronautics and Astronautics, Inc.

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.