Inviscid,laminar and turbulent opposed flows

This paper attempts to reproduce numerically previous experimental findings with opposed flows and extends their range to quantify the effects of upstream pipes and nozzles with inviscid, laminar and turbulent flows. The choice of conservation equations, boundary conditions, algorithms for their solution, the degree of grid dependence, numerical diffusion and the validity of numerical approximations are justified with supporting calculations where necessary. The results of all calculations on the stagnation plane show maximum strain rates close to the annular exit from the nozzles and pipes for lower separations and it can be expected that corresponding reacting flows will tend to extinguish in this region with the extinction moving towards the axis. With laminar flows, the maximum strain rate increased with Reynolds number and the maximum values were generally greater than with inviscid flows and smaller than with turbulent flows. With large separations, the strain rates varied less and this explains some results with reacting flows where the extinction appeared to begin on the axis. The turbulent‐flow calculations allowed comparison of three common variants of a two‐equation first‐moment closure. They provided reasonable and useful indications of strain rates but none correctly represented the rms of velocity fluctuations on the axis and close to the stagnation plane. As expected, those designed to deal with this problem produced results in better agreement with experiment but were still imperfect. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of local variations of the laminar flame speed on the global finger-flame acceleration scenario

Numerous formulations describing the dynamics and morphology of corrugated flames, including the scenarios of flame acceleration, are based on a “geometrical consideration”, where the wrinkled-to-planar flame velocities ratio, S_w /S_L , is evaluated as the scaled flame surface area, while the entire combustion chemistry is immersed into the planar flame speed S_L , which is assumed to be constant. However, S_L may experience noticeable spatial/temporal variations in practice, in particular, due to pressure/temperature variations as well as non-uniform distribution of the equivalence ratio and/or that of combustible or inert dust impurities. The present work initiates the systematic study of the impact of the local S_L -variations on the global flame evolution scenario. The variations are assumed to be imposed externally, in a manner being a free functional of the formulation. Specifically, the linear, parabolic and hyperbolic spatial S_L -distributions are incorporated into the formulations of finger flame acceleration in pipes, and they are compared to the case of constant S_L . Both two-dimensional channels and cylindrical tubes are considered. The conditions promoting or moderating flame acceleration are identified, and the revisited equations for the flame shape, velocity and acceleration rate are obtained for various S_L -distributions. The theoretical findings are validated by the computational simulations of the reacting flow equations, with agreement between the theory and modelling demonstrated.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

Construction of a Chemical Kinetic Model of Five-Component Gasoline Surrogates under Lean Conditions

The requirements for improving the efficiency of internal combustion engines and reducing emissions have promoted the development of new combustion technologies under extreme operating conditions (e.g., lean combustion), and the ignition and combustion characteristics of fuels are increasingly becoming important. A chemical kinetic reduced mechanism consisting of 115 species and 414 elementary reactions is developed for the prediction of ignition and combustion behaviors of gasoline surrogate fuels composed of five components, namely, isooctane, n-heptane, toluene, diisobutylene, and cyclohexane (CHX). The CHX sub-mechanism is obtained by simplifying the JetSurF2.0 mechanism using direct relationship graph error propagating, rate of production analysis, and temperature sensitivity analysis and CHX is mainly consumed through ring-opening reactions, continuous dehydrogenation, and oxygenation reactions. In addition, kinetic parameter corrections were made for key reactions R14 and R391 based on the accuracy of the ignition delay time and laminar flame velocity predictions. Under a wide range of conditions, the mechanism’s ignition delay time, laminar flame speed, and the experimental and calculated results of multi-component gasoline surrogate fuel and real gasoline are compared. The proposed mechanism can accurately reproduce the combustion and oxidation of each component of the gasoline-surrogate fuel mixture and real gasoline.     

Quenching of Flame Propagation with Heat Loss

The steady propagation of a planar laminar premixed flame, with a one-step exothermic reaction and linear heat loss, is studied. The corresponding travelling wave equations are solved numerically. The dependence of the flame velocity on the heat loss parameter is determined and compared with known results obtained by asymptotic expansion and other approximations. Due to the introduction of an ignition temperature the problem can be reduced to a bounded interval (of length L) and the graph of flame speed versus heat loss parameter can be parametrised by L. The numerical method is tested in the case of a step function nonlinearity when the exact solution of the differential equations can also be calculated.     

Printed Circuit Board-Based Electrolyte Regulator for Laminar Microfluidics

The trench on a printed circuit board was reconstructed to fabricate a microfluidic framework that allows low-cost production for small quantities and integration with multifunctional elements. An on-chip electrolyte regulator was thus proposed on this platform to analyze diffusion properties in laminar flow. A numerical model was developed, highlighting the interplay between the electrolyte migration and hydrodynamic properties. Solutions with dissolved sodium chloride were simulated and experimentally tested for the regulation of electrical conductivity under the guidance of the normalized Nernst-Planck equation. The diffusion mechanism and the resulting concentration field were demonstrated in detail. This approach provides a satisfactory manufacturing method and a useful tool for integrated microfluidic systems.     

“蜡烛火焰划分”实验的局限成因分析及改进研究

人教版初中化学首个探究性实验"蜡烛及其燃烧"明确指出火焰的划分方法，但在实际教学中肉眼观察法和火柴梗燃烧法对蜡烛火焰的划分存在一定的局限性。采用文献研究法和实验验证法分析其局限性的成因——燃料的析炭能力强弱、材料结构的不同等，并提出相应的改进建议。     

Laminar Flow Emulsification Process to Control the Viscosity Reduction of Heavy Crude Oils

The formation of heavy crude oil in water (O/W) emulsion by a low energy laminar controlled flow has been investigated. The emulsion was prepared in an eccentric cylinder mixer. Its geometry allows the existence of chaotic flows that are able to mix well highly viscous fluids. This new mixer design is used to produce high internal phase ratio emulsions for three oils: castor oil and two heavy crude oils of different initial viscosity (Zuata and Athabasca crude oils). The influence of the stirring conditions, geometrical parameters, and water volume fraction on the rheological properties of the resulting O/W emulsion is studied.     

精确计算Bingham流体圆管层流压降及速度分布规律的新方法

Bingham(宾汉)模型情况下,多采用通用公式进行圆管层流压降的解析计算,即将Bingham模型本构方程代入粘性流体圆管层流流动通用公式进行计算,仅能得到压降的解析解.新方法结合Bingham流体本构方程与运动方程,建立有关力学平衡方程,并运用代数方程的根式解理论对圆管层流流动时的非线性方程进行求解,可直接求得Bingham流体圆管层流压降及速度流核区半径的解析解,进一步可求得圆管层流速度解析解;Bingham流体圆管层流速度的直接影响因素为流量、塑性粘度和屈服值,研究发现速度流核宽度与屈服值成正比,与流量及塑性粘度成反比,且流核的宽度越大,流核区的速度越小.