Flame–wall interaction effects on the flame root stabilization mechanisms of a doubly-transcritical LO2/LCH4 cryogenic flame

High-fidelity numerical simulations are used to study flame root stabilization mechanisms of cryogenic flames, where both reactants (O₂ and CH₄) are injected in transcritical conditions in the geometry of the laboratory scale test rig Mascotte operated by ONERA (France). Simulations provide a detailed insight into flame root stabilization mechanisms for these diffusion flames: they show that the large wall heat losses at the lips of the coaxial injector are of primary importance, and require to solve for the fully coupled conjugate heat transfer problem. In order to account for flame–wall interaction (FWI) at the injector lip, detailed chemistry effects are also prevalent and a detailed kinetic mechanism for CH₄ oxycombustion at high pressure is derived and validated. This kinetic scheme is used in a real-gas fluid solver, coupled with a solid thermal solver in the splitter plate to calculate the unsteady temperature field in the lip. A simulation with adiabatic boundary conditions, an hypothesis that is often used in real-gas combustion, is also performed for comparison. It is found that adiabatic walls simulations lead to enhanced cryogenic reactants vaporization and mixing, and to a quasi-steady flame, which anchors within the oxidizer stream. On the other hand, FWI simulations produce self-sustained oscillations of both lip temperature and flame root location at similar frequencies: the flame root moves from the CH₄ to the O₂ streams at approximately 450 Hz, affecting the whole flame structure.     

A numerical simulation of transient ignition and ignition limit of a composite solid by a localised radiant source

An unsteady three-dimensional numerical model has been formulated, coded, and solved to study ignition and flame development over a composite solid fuel sample upon heating by a localised radiant beam in a buoyant atmosphere. The model consists of an unsteady gas phase and an unsteady solid phase. The gas phase formulation consists of full Navier-Stokes equations for the conservation of mass, momentum, energy, and species. A one-step, second-order overall Arrhenius reaction is adopted. Gas radiation is included by solving the radiation transfer equation. For the solid phase formulation, the energy (heat conduction) equation is employed to solve the transient solid temperature. A first-order in-depth solid pyrolysis relation between the solid fuel density and the local solid temperature is assumed. Numerical simulations provide time-and-space resolved details of the ignition transient and flame development and the existence of two types of ignition modes: one with reaction kernel initiated on the surface and the other with ignition kernel initiated in the gas phase. Other primary outputs of the computation are the minimum ignition energy (Joule) for the solid as a function of the external heating rate (Watt). Both the critical heat input for ignition and the optimal ignition energy are identified. Other parameters that were varied over the simulations include: sample thickness, ignition heat source spatial shape factor, and gravity level.     

振荡流共轭换热格子——Boltzmann模拟

振荡流共轭换热现象广泛存在于热声热机等工程应用中.基于双分布格子-Boltzmann模型,对平行平板间振荡流共轭换热进行了数值模拟.通过假定共轭界面处流体和固体的未知内能分布函数均为对应的平衡态滑移修正格式,提出了一种处理共轭换热边界的新方法.模拟结果表明,该方法可以保证共轭界面上温度连续和热流连续.分析了不同流体与固体导热系数比情况下振荡流共轭换热的速度场、温度场以及热流分布的特点.     

耦合传热和冷却流通道流动对气膜冷却的影响

本文通过数值模拟的方法同时考虑了耦合传热和冷却流通道流动对气膜冷却的影响.计算结果表明,在考虑耦合传热的情况下,冷却流通道流动的影响仍然存在,但随着壁面导热系数的增大,这种影响减弱;同时在考虑耦合传热的情况下,受保护壁面温度场分布更加均匀,冷却效果更好.计算结果还表明吹风比为0.5时的冷却效果优于吹风比为1.0的情况.     

A composite grid solver for conjugate heat transfer in fluid–structure systems

We describe a numerical method for modeling temperature-dependent fluid flow coupled to heat transfer in solids. This approach to conjugate heat transfer can be used to compute transient and steady state solutions to a wide range of fluid–solid systems in complex two- and three-dimensional geometry. Fluids are modeled with the temperature-dependent incompressible Navier–Stokes equations using the Boussinesq approximation. Solids with heat transfer are modeled with the heat equation. Appropriate interface equations are applied to couple the solutions across different domains. The computational region is divided into a number of sub-domains corresponding to fluid domains and solid domains. There may be multiple fluid domains and multiple solid domains. Each fluid or solid sub-domain is discretized with an overlapping grid. The entire region is associated with a composite grid which is the union of the overlapping grids for the sub-domains. A different physics solver (fluid solver or solid solver) is associated with each sub-domain. A higher-level multi-domain solver manages the entire solution process.     

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.     

A stable and high-order accurate conjugate heat transfer problem

This paper analyzes well-posedness and stability of a conjugate heat transfer problem in one space dimension. We study a model problem for heat transfer between a fluid and a solid. The energy method is used to derive boundary and interface conditions that make the continuous problem well-posed and the semi-discrete problem stable. The numerical scheme is implemented using 2nd-, 3rd- and 4th-order finite difference operators on Summation-By-Parts (SBP) form. The boundary and interface conditions are implemented weakly. We investigate the spectrum of the spatial discretization to determine which type of coupling that gives attractive convergence properties. The rate of convergence is verified using the method of manufactured solutions.     

Lattice Boltzmann simulation of viscous-fluid flow and conjugate heat transfer in a rectangular cavity with a heated moving wall

In the present work, conjugate heat transfer in a rectangular cavity with a heated moving lid is investigated using the lattice Boltzmann method (LBM). The simulations are performed for incompressible flow, with Reynolds numbers ranging from 100 to 500, thermal diffusivity ratios ranging from 1 to 100, and Prandtl numbers ranging from 0.7 to 7. A uniform heat flux through the top of the lid is assumed. Results show that LBM is suitable for the study of heat transfer in conjugate problems. Effects of the Reynolds number, the Prandtl number and the thermal diffusivity ratio on hydrodynamic and thermal characteristics are investigated and discussed. The streamlines and temperature distribution in flow field, dimensionless temperature and Nusselt number along the hot wall are illustrated. The results indicate that increase of thermal diffusivity yields the removal of a higher quantity of energy from lid and its temperature decreases when increasing the Reynolds and the Prandtl numbers.     

Stability of a spherical flame ball in a porous medium

Gaseous flame balls and their stability to symmetric disturbances are studied numerically and asymptotically, for large activation temperature, within a porous medium that serves only to exchange heat with the gas. Heat losses to a distant ambient environment, affecting only the gas, are taken to be radiative in nature and are represented using two alternative models. One of these treats the heat loss as being constant in the burnt gases and linearizes the radiative law in the unburnt gas (as has been studied elsewhere without the presence of a solid). The other does not distinguish between burnt and unburnt gas and is a continuous dimensionless form of Stefan's law, having a linear part that dominates close to ambient temperatures and a fourth power that dominates at higher temperatures.

Numerical results are found to require unusually large activation temperatures in order to approach the asymptotic results. The latter involve two branches of solution, a smaller and a larger flame ball, provided heat losses are not too high. The two radiative heat loss models give completely analogous steady asymptotic solutions, to leading order, that are also unaffected by the presence of the solid which therefore only influences their stability. For moderate values of the dimensionless heat-transfer time between the solid and gas all flame balls are unstable for Lewis numbers greater than unity. At Lewis numbers less than unity, part of the branch of larger flame balls becomes stable, solutions with the continuous radiative law being stable over a narrower range of parameters. In both cases, for moderate heat-transfer times, the stable region is increased by the heat capacity of the solid in a way that amounts, simply, to decreasing an effective Lewis number for determining stability, just as if the heat-transfer time was zero.     

Modelling nongrey gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modelling with respect to the treatment of nongrey radiation from both gas-phase species and soot particles, while radiation modelling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (M.F. Modest, 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured computational fluid dynamics (CFD) code for nongrey radiation modelling. The mixture full-spectrum k-distributions including nongrey absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongrey soot modelling is shown to be of greater importance than nongrey gas modelling in sooty flame simulations, with grey soot models producing large errors. The nongrey treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongrey treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongrey soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Modeling nongray gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modeling with respect to treatment of nongray radiation from both gas-phase species and soot particles, while radiation modeling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (Modest, M.F., 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured CFD code for nongray radiation modeling. The mixture full-spectrum k-distributions including nongray absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical-harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongray soot modeling is shown to be of greater importance than nongray gas modeling in sooty flame simulations, with gray soot models producing large errors. The nongray treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongray treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongray soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Comparisons of Heat Transfer Enhancement of an Internal Blade Tip with Metal or Insulating Pins

Cooling methods are needed for turbine blade tips to ensure a long durability and safe operation. A common way to cool a tip is to use serpentine passages with 180-deg turn under the blade tip-cap taking advantage of the three-dimensional turning effect and impingement like flow. Improved internal convective cooling is therefore required to increase the blade tip lifetime. In the present study, augmented heat transfer of an internal blade tip with pin-fin arrays has been investigated numerically using a conjugate heat transfer method. The computational domain includes the fluid region and the solid pins as well as the tip regions. Turbulent convective heat transfer between the fluid and pins, and heat conduction within pins and tip are simultaneously computed. The main objective of the present study is to observe the effect of the pin material on heat transfer enhancement of the pin-finned tips. It is found that due to the combination of turning, impingement and pin-fin crossflow, the heat transfer coefficient of a pin-finned tip is a factor of 2.9 higher than that of a smooth tip at the cost of an increased pressure drop by less than 10%. The usage of metal pins can reduce the tip temperature effectively and thereby remove the heat load from the tip. Also, it is found that the tip heat transfer is enhanced even by using insulating pins having low thermal conductivity at low Reynolds numbers. The comparisons of overall performances are also included.     

Simulation study of a hot metal cylinder cooling by gas-liquid flow

A mathematical model was developed for conjugate heat transfer in a heterogeneous system “solid body ? gas-liquid medium” with account for vapor generation at the surface of hot metal cylinder with cooling by a longitudinal water flow. Results are presented for numerical parametric calculations for influence of thermophysical and hydrodynamic characteristics on the pattern of vapor generation at the cooled cylinder surface.     

A study of the heat transfer characteristics of a laser rod in a concentric circular tube

An analytical solution was derived from a two-dimensional heat conduction model with non-uniform boundary conditions for a side-pumped, side-cooled cylindrical laser rod. The convective heat transfer coefficient and the spatially varying fluid temperature were determined from the theoretical solutions or experimental correlations of convective heat transfer in an annular passage with prescribed heat fluxes. The first term of the analytical solution coincides with the result of the one-dimensional model. The other terms indicate that the axial temperature-rise in a laser rod has positive correlations with the axial coolant temperature-rise, the radial Biot number and the length-to-radius ratio of the laser rod. Subsequently, a conjugate numerical simulation that couples up the fluid convection and the solid conduction was performed. Compared with the analytical solution, the conjugate numerical simulation better exhibits the entrance effects of flow and heat transfer; therefore, it may provide more accurate solution in specific cases. PACS 44.10.+i; 44.90.+c     

粗糙平板微通道流动和传热的数值模拟

以粗糙平行平板微通道为研究对象,以三角形锯齿状粗糙元模拟固体表面的粗糙度,采用CFD流体固体共轭传热技术数值研究了绝对粗糙度和相对粗糙度对平行平板微通道流动和传热特性的影响,着重分析了粗糙度和流体速度对水力入口段长度和热力入口段长度的影响规律,同时研究了相对粗糙度对微通道转捩雷诺数的影响,为进一步揭示微微通道的流动和传热机理提供了依据.     

Nonlinear equation for curved nonstationary flames and flame stability

A time-dependent nonlinear equation for a nonstationary curved flame front of an arbitrary expansion coefficient is derived under the assumptions of a small but finite flame thickness and weak nonlinearity. On the basis of the derived equation, stability of two-dimensional curved stationary flames propagating in tubes with ideally adiabatic and slip walls is studied. The stability analysis shows that curved stationary flames become unstable for sufficiently wide tubes. The obtained stability limits are in a good agreement with the results of numerical simulations of flame dynamics and with semiqualitative stability analysis of curved stationary flames. Possible outcomes of the obtained instability at the nonlinear stage are discussed. The instability may result in extra wrinkles at a flame front close to the stability limits and in self-turbulization of the flame far from the limits. The self-turbulization can also be interpreted as a fractal structure. The fractal dimension of a flame front and velocity of a self-turbulized flame are evaluated.     

On disintegration of near-limit cellular flames

A strongly non-linear geometrically-invariant model for the dynamics of near-limit cellular flame is proposed, where the flame evolution is governed by a system of equations for the flame interface and its temperature. The model generalizes its earlier weakly non-linear version pertinent to a mildly perturbed planar flame. Numerical simulations of the new model show that at sufficiently high levels of heat losses the cellular flame resulting from the diffusive instability exhibits a tendency toward self-fragmentation, quite in line with direct numerical simulations of the associated reaction-diffusion system.     

Combustion driven oscillations

Combustion driven oscillations can occur when a turbulent flame is enclosed in a tube or cavity. Interaction between heat fluctuations and the internal standing wave field at one of the natural frequencies of the air column produces strong organ pipe tones. The sound power emitted by this thermal-acoustic interaction depends on the impedance either side of the combustion zone and on a transfer function defining the response of the flame to sound wave disturbances. If this power exceeds the rate at which energy is dissipated at the cavity boundaries then there is a growth of the internal pressure field and an increase in the radiated sound. Plane wave theory is used to calculate the flame transfer function and adjacent impedances for a simple gas fired tube assembly. The predicted instability frequencies are then compared with experimental data. The results indicate that the flame transfer function plays a dominant role in determining the acoustic stability of the cavity and that insufficient data is available for accurately predicting unsteady flame front behaviour.     

Influences of heat flux on extinction characteristics of steady/unsteady premixed stagnation flames

This paper investigates the extinction characteristics of premixed stagnation flames (PSFs) with controlled heat losses and flow disturbances. The low-frequency air flow pulsations that imitate the operational transients in practical combustors were specially introduced. The tunable diode-laser absorption spectroscopy (TDLAS) measurement was applied to obtain the temperature profile and wall heat flux. It is found that, for steady flame with a fixed equivalence ratio, the extinction stretch rate dramatically increases as the wall heat flux decreases. The extinction criterion is summarized as a global Karlovitz number of 0.57 by establishing a relationship between the global and local stretch rates. Numerical simulations reveal that the local extinction Karlovitz number of steady PSFs is approximately 1.0 regardless of the conditions such as heat flux and equivalence ratio. Further experiments present that the air pulsations with a repetition of ～5 Hz significantly deteriorate the flame stability. Particularly, for unsteady perturbed flames, the extinction stretch rate exhibits a nonlinear trend, yielding two regimes with discrepant sensitivities to wall heat flux. The unsteady simulation then highlights a local stretch rate overshoot in the presence of pulsation. It is caused by the time delay between the inlet velocity and flame front movement that eventually leads to poor flame stability. Moreover, in the high heat-flux regime, a smaller local stretch rate overshoot results in the weak dependence of extinction limits on heat fluxes.     

Opposed flow flame spread over an array of thin solid fuel sheets in a microgravity environment

In this work a numerical study has been carried out to gain physical insight into the phenomena of opposed flow flame spread over an array of thin solid fuel sheets in a microgravity environment. The two-dimensional (2D) simulations show that the flame spread rates for the multiple-fuel configuration are higher than those for the flame spreading over a single fuel sheet. This is due to reduced radiation losses from the flame and increased heat feedback to the solid fuel. The flame spread rate exhibits a non-monotonic variation with decrease in the interspace distance between the fuel sheets. Higher radiation heat feedback primarily as gas/flame radiation was found to be responsible for the increase in the flame spread rate with the reduction of the interspace distance. It was noted that as the interspace distance between the fuel sheets was reduced below a certain value, no steady solution could be obtained. However, at very small interspace distances, steady state spread rates were obtained. Here, due to oxygen starvation the flame spread rate decreased and eventually at some interspace distance the flame extinguished. With fuel emittance (equal to absorptance) reduced to ‘0’ the flame spread rate was nearly independent of the interspace distance, except at very small distances where the flame spread rate dropped due to oxygen starvation. A flame extinction plot with the extinction oxygen level was constructed for the multiple-fuel configuration at various interspace distances. The default fuel with an emittance of 0.92 was found to be more flammable in the multiple-fuel configuration than in a single fuel sheet configuration. For a fuel emittance equal to zero, the extinction oxygen limit decreases for both the single and the multiple fuel sheet configurations. However, the two flammability curves cross over at a certain fuel separation distance. The multiple-fuel configurations become less flammable compared to the single fuel sheet configuration below a certain separation distance.