Triple-flame propagation in a parallel flow: an analytical study

We present an analytical study of triple-flame propagation in a two-dimensional mixing layer against a parallel flow. The problem is formulated within a constant density thermo-diffusive model, and solved analytically in the asymptotic limit of large activation energy of the chemical reaction for flames thin compared with their typical radius of curvature. Explicit expressions are obtained in this limit, describing the influence of the flow on the triple-flame. The results are expected to be applicable when the ratio between the flow-scale and the flame-front radius of curvature (which is mainly dictated by concentration gradients) is of order unity, or larger. When this ratio is large, as in the illustrative case of a Poiseuille flow in a porous channel considered here, the flow is found to negligibly affect the flame structure except for a change in its speed by an amount which depends on the stoichiometric conditions of the mixture. On the other hand, when this ratio is of order unity, the flow is able to significantly wrinkle the flame-front, modify its propagation speed, and shift its leading edge away from the stoichiometric line. The latter situation is investigated in the illustrative case of spatially harmonic flows. The results presented describe, in particular, how the leading-edge of the flame-front can be determined in terms of the flow amplitude A which is critical in determining the flame speed. The latter is found to depend linearly on A in the first approximation with a correction proportional to the flame thickness multiplied by , for |A| sufficiently large. The effect of varying the flow-scale on flame propagation in this context is also described, with explicit formulae provided, and interesting behaviours, such as non-monotonic dependence on the scale, identified.     

Hydrodynamic instability of premixed flame propagating in narrow planar channel in the presence of gas flow

The paper analyses the hydrodynamic instability of a flame propagating in the space between two parallel plates in the presence of gas flow. The linear analysis was performed in the framework of a two-dimensional model that describes the averaged gas flow in the space between the plates and the perturbations development of two-dimensional combustion wave. The model includes the parametric dependences of the flame front propagation velocity on its local curvature and on the combustible gas velocity averaged along the height of the channel. It is assumed that the viscous gas flow changes the surface area of the flame front and thereby affects the propagation velocity of the two-dimensional combustion wave. In the absence of the influence of the channel walls on the gas flow, the model transforms into the Darrieus–Landau model of flame hydrodynamic instability. The dependences of the instability growth rate on the wave vector of disturbances, the velocity of the unperturbed gas flow, the viscous friction coefficients and other parameters of the problem are obtained. It is shown that the viscous gas flow in the channel can lead, in some cases, to a significant increase in instability compared with a flame propagating in free space. In particular, the instability increment depends on the direction of the gas flow with respect direction of the flame propagation. In the case when the gas flow moves in the opposite direction to the direction of the flame propagation, the pulsating instability can appear.     

Flame balls dynamics in divergent channel

A three-dimensional reaction-diffusion model for lean low-Lewis-number premixed flames with radiative heat losses propagating in divergent channel is studied numerically. Effects of inlet gas velocity and heat-loss intensity on flame structure at low Lewis numbers are investigated. It is found that continuous flame front exists at small heat losses and the separate flame balls settled within restricted domain inside the divergent channel at large heat losses. It is shown that the time averaged flame balls coordinate may be considered as important characteristic analogous to coordinate of continuous flame stabilized in divergent channel.     

Molecular dynamic simulation of Copper and Platinum nanoparticles Poiseuille flow in a nanochannels

In this paper, simulation of Poiseuille flow within nanochannel containing Copper and Platinum particles has been performed using molecular dynamic (MD). In this simulation LAMMPS code is used to simulate three-dimensional Poiseuille flow. The atomic interaction is governed by the modified Lennard–Jones potential. To study the wall effects on the surface tension and density profile, we placed two solid walls, one at the bottom boundary and the other at the top boundary. For solid–liquid interactions, the modified Lennard–Jones potential function was used. Velocity profiles and distribution of temperature and density have been obtained, and agglutination of nanoparticles has been discussed. It has also shown that with more particles, less time is required for the particles to fuse or agglutinate. Also, we can conclude that the agglutination time in nanochannel with Copper particles is faster that in Platinum nanoparticles. Finally, it is demonstrated that using nanoparticles raises thermal conduction in the channel.     

Flame propagation in a small-scale parallel flow

We consider the propagation of laminar premixed flames in the presence of a parallel flow whose scale is smaller than the laminar flame thickness. The study addresses fundamental aspects with relevance to flame propagation in narrow channels, to the emerging micro-combustion technology, and to the understanding of the effect of small scales in a (turbulent) flow on the flame structure. In part, the study extends the results of a previous analytical study carried out in the thick flame asymptotic limit which has in particular addressed the validity of Damköhler's second hypothesis in the context of laminar steady parallel flows. Several new contributions are made here.

Analytical contributions include the derivation of an explicit formula for the effective speed of a premixed flame U _T in the presence of an oscillatory parallel flow whose scale ? (measured with the laminar flame thickness δ _L ) is small and amplitude A (measured with the laminar flame speed U _L ) is (1). The formula shows a quadratic dependence on both the amplitude and the scale of the flow. The validity of the formula is established analytically in two distinguished limits corresponding to (1) frequencies of oscillations (measured with the natural frequency of the flame U _L /δ _L ), and to higher frequencies of (A/?) (the natural frequency of the flow). The analytical study yields partial support of Damköhler's second hypothesis in that it shows that the flame behaves as a planar flame (to leading order) with an increased propagation speed which depends on both the scale and amplitude of the velocity fluctuation. However our formula for U _T contradicts the formula given by Damköhler in his original paper where U _T has a square root dependence on the scale and amplitude.

Numerical contributions include a significant set of two-dimensional calculations which determine the range of validity of the asymptotic findings. In particular, these account for volumetric heat loss and differential diffusion effects. Good agreement between the numerics and asymptotics is found in all cases, both for steady and oscillatory flows, at least in the expected range of validity of the asymptotics. The effect of the frequency of oscillation is also discussed. Additional related aspects such as the difference in the response of thin and thick flames to the combined effect of heat loss and fluid flow are also addressed. It is found for example that the sensitivity of thick flames to volumetric heat loss is negligibly affected by the parallel flow intensity, in marked contrast to the sensitivity of thin flames. Interestingly, and somewhat surprisingly, thin flames are found to be more resistant to heat loss when a flow is present, even for unit Lewis number; this ceases to be the case, however, when the Lewis number is large enough.     

Mathematical simulation of sound propagation in a flow channel with impedance walls

The paper considers the specifics of calculating tonal sound propagating in a flow channel with an installed sound-absorbing device. The calculation is performed on the basis of numerical integrating on linearized nonstationary Euler equations using a code developed by the authors based on the so-called discontinuous Galerkin method. Using the linear theory of small perturbations, the effect of the sound-absorbing lining of the channel walls is described with the modified value of acoustic impedance proposed by the authors, for which, under flow channel conditions, the traditional classification of the active and reactive types of lining in terms of the real and imaginary impedance values, respectively, remains valid. To stabilize the computation process, a generalized impedance boundary condition is proposed in which, in addition to the impedance value itself, some additional parameters are introduced characterizing certain fictitious properties of inertia and elasticity of the impedance surface.     

Slow hydromagnetic flow in a channel with porous walls in presence of a periodic magnetic field

An attempt has been made here to study the MHD effects on the slow motion of a viscous, incompressible and conducting fluid between two parallel porous infinite plane walls in presence of a transverse magnetic field varying periodically with time. The problem has been investigated firstly for the case of non-conducting walls and finally for the case of conducting walls.Notation biV velocity vector - U ₀ prescribed velocity - biH magnetic field vector - t time - density - coefficient of viscosity - coefficient of kinematic viscosity - conductivity of the medium - magnetic permeability - magnetic viscosity - _w conductivity of the wall - l thickness of the wall - L a constant characterising the thickness of the fluid slab - electric conductance ratio,= _w l/(L) - frequency of theH ₀ variations - P a constant, –1/p/x=P - (x, y, z) cartesion co-ordinates - Y dimensionless cartesian co-ordinate corresponding toy - M Hartmann number - R Reynolds number - R _m magnetic Reynolds number - R _c Reynolds number for cross-flow - S magnetic pressure number - P dimensionless pressure gradient - ¯t dimsnsionless time - Ps a parameter     

The flow of a micropolar fluid through a porous channel with expanding or contracting walls

The flow of a micropolar fluid in a porous channel with expanding or contracting walls is investigated. The governing equations are reduced to ordinary ones by using similar transformations. Homotopy analysis method (HAM) is employed to obtain the expressions for the velocity fields and microrotation fields. Graphs are sketched for the effects of some values of parameters, especially the expansion ratio, on the velocity and microrotation fields and associated dynamic characteristics are analyzed in detail.     

On combustion in a closed rectangular channel with initial vorticity

This article examines the detailed combustion process in a theoretical model with applicability to combustion in a wave rotor or wave disc engine. The model comprises a single channel into which an initial loading of methane and air is admitted and ignited after all inlet and exit ports have been closed. Combustion takes place at constant volume. However, the initial gaseous mixture in the channel is not at rest. The initial opening and closing of the ports generates significant vorticity which influences the evolution of the combustion process. Numerical evaluations are provided for the detailed flame shape for simplified (one-step) chemistry and a simulation using the detailed 235-step San Diego scheme. Quantities examined are the evolution of vorticity, pressure fluctuations, mass consumption rate, flame surface area and the influences on combustion of adiabatic and non-adiabatic channel walls. Combustion regimes are identified and compared with simpler model studies (no initial flow). Pointwise quantities are examined to describe the various stages of burning in the channel. The focus of the study is directed towards quantities that influence overall burning rate and completeness of combustion.     

Structural response of different Lewis number premixed flames interacting with a toroidal vortex

Simultaneous measurements of temperature, CH* and OH* chemiluminescent species are carried out to explore the impact of stretch rate and curvature on the structure of premixed flames. The configuration of an initially flat premixed flame interacting with a toroidal vortex is selected for the present study and reasons for this choice are discussed. Lewis number effects are assessed by comparing methane and propane flames. It is emphasized that the flame structure experiences very strong variations. In particular, the flame is shrunk (broadened) in the initial (final) period of the interaction with the vortex where strain rate (curvature) contribution of the stretch rate is predominant. By further analysing independently the thickness of the preheat and reaction zones, it is shown that for propane flames, not only the former but also the latter is significantly altered in zones where the flame curvature is negative. Changes in the reaction zone properties are further emphasized using CH* and OH* radicals. It is demonstrated that higher thermal diffusivity plays a significant role around curved regions, in which the enhanced diffusion of heat leads to a strong increase of CH* compared to OH* intensity. As an overall conclusion, this study suggests that it would be interesting to reassess the internal flame structure at lower and moderate Karlovitz numbers since changes might appear for a moderate vortex intensity with typical size much larger than the flame thickness.     

Effects of spanwise system rotation on turbulent stripes in a plane Poiseuille flow

In the transitional channel flow, the large-scale intermittent structure of localised turbulence, which is called the turbulent stripe pattern, can be found in the form of stripe arrangement. The structure of the turbulent stripe pattern is an oblique laminar–turbulent banded pattern and is inclined with respect to the streamwise direction. We performed direct numerical simulation at a transitional Reynolds number and very low-rotation numbers, and focused on the turbulent stripe pattern in the plane Poiseuille flow subjected to spanwise system rotation. We captured the turbulent stripe pattern in a rotating channel flow and found the augmentation and diminution of the turbulent stripe pattern were affected by the spanwise rotation. The contents of the discussion are the spatial size of the turbulent stripe pattern on the basis of the instantaneous flow fields, the energy spectra, and various statistics relating to the spanwise velocity component that characterise the turbulent stripe pattern. The turbulent stripe pattern was found to contain kinetic energy that was larger in very weakly rotating flows than in the static system. It was also found that the magnitude of the spanwise secondary flow increases, while the quasi-laminar region is wider at a very lowrotation number.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

Stationary combustion regimes and extinction limits of one-dimensional stretched premixed flames in a gap between two heat conducting plates

Stationary combustion regimes, their linear stability and extinction limits of stretched premixed flames in a narrow gap between two heat conducting plates are studied by means of numerical simulations in the framework of one-dimensional thermal-diffusion model with overall one-step reaction. Various stationary combustion modes including normal flame (NF), near-stagnation plane flame (NSF), weak flame (WF) and distant flame (DF) are detected and found to be analogous to the same-named regimes of conventional counterflow flames. For the flames stabilized in the vicinity of stagnation plane at moderate and large stretch rates (which are NF, NSF and WF) the effect of channel walls is basically reduced to additional heat loss. For distant flame characterized by large flame separation distance and small stretch rates intensive interphase heat transfer and heat recirculation are typical. It is shown that in mixture content / stretch rate plane the extinction limit curve has ε-shape, while for conventional counterflow flames it is known to be C-shaped. This result is quite in line with recent experimental findings and is explained by extension of extinction limits at small stretch rates at the expense of heat recirculation. Analysis of the numerical results makes possible to reveal prime mechanisms of flame quenching on different branches of ε-shaped extinction limit curve. Namely, two upper limits are caused by stretch and heat loss. These limits are direct analogs of the upper and lower limits on conventional C-shaped curve. Two other limits are related with weakening of heat recirculation and heat dissipation to the burner. Thus, the present study provides a satisfactory explanation for the recent experimental observations of stretched flames in narrow channel.     

Numerical simulation of heat transfer and fluid flow of Water-CuO Nanofluid in a sinusoidal channel with a porous medium

This study has compared the convection heat transfer of Water-based fluid flow with that of Water-Copper oxide (CuO) nanofluid in a sinusoidal channel with a porous medium. The heat flux in the lower and upper walls has been assumed constant, and the flow has been assumed to be two-dimensional, steady, laminar, and incompressible. The governing equations include equations of continuity, momentum, and energy. The assumption of thermal equilibrium has been considered between the porous medium and the fluid. The effects of the parameters, Reynolds number and Darcy number on the thermal performance of the channel, have been investigated. The results of this study show that the presence of a porous medium in a channel, as well as adding nanoparticles to the base fluid, increases the Nusselt number and the convection heat transfer coefficient. Also the results show that As the Reynolds number increases, the temperature gradient increases. In addition, changes in this parameter are greater in the throat of the flow than in convex regions due to changes in the channel geometry. In addition, porous regions reduce the temperature difference, which in turn increases the convective heat transfer coefficient.     

Strain rate and flame orientation statistics in the near-wall region for turbulent flame-wall interaction

Flame–wall interaction (FWI) in premixed turbulent combustion has been analysed based on a counter-flow like configuration at the statistically stationary state. For the present configuration, the two FWI sub-zones, i.e the influence zone and the quenching zone, can be identified from the DNS results. Detailed analysis of the important quantities, such as the flame temperature, flame–wall distance, wall heat flux, flame curvature and dilatation (including the flame normal and tangential strain rates), and some orientation relations between the flame normal and the principal strain rate directions, have been reported, together with the physical explanations. All these statistical results are determined by the relative strengths of the wall heat flux, thermal expansion and the flame–turbulence interaction.     

Allowable propagation of short pulse laser beam in a plasma channel and electromagnetic solitary waves

Nonparaxial and nonlinear propagation of a short intense laser beam in a parabolic plasma channel is analyzed by means of the variational method and nonlinear dynamics. The beam propagation properties are classified by five kinds of behaviors. In particularly, the electromagnetic solitary wave for finite pulse laser is found beside the other four propagation cases including beam periodically oscillating with defocussing and focusing amplitude, constant spot size, beam catastrophic focusing. It is also found that the laser pulse can be allowed to propagate in the plasma channel only when a certain relation for laser parameters and plasma channel parameters is satisfied. For the solitary wave, it may provide an effective way to obtain ultra-short laser pulse.