Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

The aeroacoustics of transonic helicopter blades

One of the most important issues facing the helicopter industry today is helicopter noise, in particular transonic rotor noise. It is the main factor limiting cruise speeds, and there is real demand for efficient and reliable prediction methods which can be used in the rotor design process. This paper considers the Ffowcs Williams-Hawkings equation applied to a permeable control surface. The surface is chosen to be as small as possible, while enclosing both the blade and any transonic flow regions. This allows the problematic quadrupole term to always be neglected, and requires only near field CFD input data. It is therefore less computationally intensive than existing prediction methods, and moreover retains the physical interpretation of the sources in terms of thickness, loading and shock-associated noise. A computer program has been developed which implements the permeable surface form of retarded time formulation. The program has been validated and subsequently used to validate an acoustic 2-D CFD code. It is fast and reliable for subsonic motion, but it is demonstrated that it cannot be used at high subsonic or supersonic speeds. A second computer program implementing a more general formulation has also been developed and is presently being validated. This general formulation can be applied at high subsonic and supersonic speeds, except under one specific condition. © 2002 by the author(s). Published by the American Institute of Aeronautics and Astronautics, Inc.

Discrete multicomponent model for biodiesel spray combustion simulation

An important first step in spray combustion simulation is an accurate determination of the fuel properties which affects the modelling of spray formation and reaction. In a practical combustion simulation, the implementation of a multicomponent model is important in capturing the relative volatility of different fuel components. A Discrete Multicomponent (DM) model is deemed to be an appropriate candidate to model a composite fuel like biodiesel which consists of four components of fatty acid methyl esters (FAME). In this paper, the DM model is compared with the traditional Continuous Thermodynamics (CTM) model for both diesel and biodiesel. The CTM model is formulated based on mixing rules that incorporate the physical and thermophysical properties of pure components into a single continuous surrogate for the composite fuel. The models are implemented within the open-source CFD code OpenFOAM, and a semi-quantitative comparison is made between the predicted spray-combustion characteristics and optical measurements of a swirl-stabilised flame of diesel and biodiesel. The DM model performs better than the CTM model in predicting a higher magnitude of heat release rate in the top flame brush region of the biodiesel flame compared to that of the diesel flame. Using both the DM and CTM models, the simulation successfully reproduces the droplet size, volume flux, and droplet density profiles of diesel and biodiesel. The DM model predicts a longer spray penetration length for biodiesel compared to that of diesel, as seen in the experimental data. Also, the DM model reproduces a segregated biodiesel fuel vapour field and spray in which the most abundant FAME component has the longest vapour penetration. In the biodiesel flame, the relative abundance of each fuel component is found to dominate over the relative volatility in terms of the vapour species distribution and vice versa in the liquid species distribution. © 2014 Elsevier Ltd. All rights reserved.

Application of fast oxygen sensors for investigations into air-path dynamics and EGR distribution in a diesel engine

The control of NOX emissions by exhaust gas recirculation (EGR) is of widespread application. However, despite dramatic improvements in all aspects of engine control, the subtle mixing processes that determine the cylinder-to-cylinder distribution of the recirculated gas often results in a mal-distribution that is still an issue for the engine designer and calibrator. In this paper we demonstrate the application of a relatively straightforward technique for the measurement of the absolute and relative dilution quantity in both steady state and transient operation. This was achieved by the use of oxygen sensors based on standard UEGO (universal exhaust gas oxygen) sensors but packaged so as to give good frequency response (∼ 10 ms time constant) and be completely insensitivity to the sample pressure and temperature. Measurements can be made at almost any location of interest, for example exhaust and inlet manifolds as well as EGR path(s), with virtually no flow disturbance. At the same time, the measurements yield insights into air-path dynamics. We argue that "dilution", as indicated by the deviation of the oxygen concentration from that of air, is a more appropriate parameter than EGR rate in the context of NOX control, especially for diesel engines. Experimental results are presented for the EGR distribution in a current production light duty 4-cylinder diesel engine in which significant differences were found in the proportion of the recirculated gas that reached each cylinder. Even the individual inlet runners of the cylinders exhibited very different dilution rates - differences of nearly 50% were observed at some conditions. An application of such data may be in the improvement of calibration and validation of CFD and other modelling techniques. Copyright © 2014 SAE International.

Influence of turbulence-chemistry interaction for n-heptane spray combustion under diesel engine conditions with emphasis on soot formation and oxidation

The influence of the turbulence-chemistry interaction (TCI) for n-heptane sprays under diesel engine conditions has been investigated by means of computational fluid dynamics (CFD) simulations. The conditional moment closure approach, which has been previously validated thoroughly for such flows, and the homogeneous reactor (i.e. no turbulent combustion model) approach have been compared, in view of the recent resurgence of the latter approaches for diesel engine CFD. Experimental data available from a constant-volume combustion chamber have been used for model validation purposes for a broad range of conditions including variations in ambient oxygen (8-21% by vol.), ambient temperature (900 and 1000 K) and ambient density (14.8 and 30 kg/m3). The results from both numerical approaches have been compared to the experimental values of ignition delay (ID), flame lift-off length (LOL), and soot volume fraction distributions. TCI was found to have a weak influence on ignition delay for the conditions simulated, attributed to the low values of the scalar dissipation relative to the critical value above which auto-ignition does not occur. In contrast, the flame LOL was considerably affected, in particular at low oxygen concentrations. Quasi-steady soot formation was similar; however, pronounced differences in soot oxidation behaviour are reported. The differences were further emphasised for a case with short injection duration: in such conditions, TCI was found to play a major role concerning the soot oxidation behaviour because of the importance of soot-oxidiser structure in mixture fraction space. Neglecting TCI leads to a strong over-estimation of soot oxidation after the end of injection. The results suggest that for some engines, and for some phenomena, the neglect of turbulent fluctuations may lead to predictions of acceptable engineering accuracy, but that a proper turbulent combustion model is needed for more reliable results. © 2014 Taylor & Francis.

集热栅对太阳能烟囱发电系统效率的影响

采用计算流体动力学（CFD）方法对太阳能烟囱发电装置进行数值模拟，得到装置内部的温度场、速度场、压力场等分布情况。对集热棚的各种几何和物理参数进行研究和分析。结果表明，集热棚直径、太阳辐照强度、覆盖材料的透明度等诸多参数对系统效率有直接而重要的影响。

迷宫流道转角对灌水器水力性能的影响

为研究齿形、梯形以及矩形流道转角变化对水力性能的影响,采用Fluent软件对不同形状下不同转角的流道进行了数值模拟。研究结果表明:当其他条件相同时,转角的变化与流量系数、流态指数呈负相关,其变化对梯形流道灌水器的流量系数影响最大,最多下降了19.03%,齿形流道次之,下降了10.14%,矩形流道是梯形流道转角角度增加的延伸,具有相同的水力性能变化规律;随着角度的增加,梯形流道总的局部水头损失系数最多增加了32.5%,而齿形流道总的局部水头损失系数最多增加了23.4%,变化都很明显;压力较高时,摩阻系数基本保持不变,流体为紊流状态。

基于混合发展的超声速段射流式COIL的数值优化

应用CFD技术,发展三维多组分化学反应流计算程序,对采用超声速段射流方式的氧碘化学激光进行数值模拟,考察分解率与增益等参考变量的详细三维分布.计算结果说明,超声速区域的高速流动以及混合效率降低使COIL无法在有限空间内完成整个运转流程;提高碘含量以加快反应速度的手段导致主流无法提供足够的载能介质,无法形成合理增益.在不改变喷管长度的前提下,提出主流无载气方式的探索,结果证明了超声速段射流方式采用无载气主流配置的可行性与优势,通过减小气流速度保证混合与化学反应在光腔上游完成,在合理的流量配比下光腔位置处可得到1.3%cm-1的增益峰值

气相规则胞格爆轰起爆与传播的统一框架理论

本文回顾了气相规则胞格爆轰起爆与传播研究的发展,根据中国科学院高温气体动力学重点实验室在爆轰现象方面的研究进展,提出了气相规则胞格爆轰波起爆和传播的统一框架理论。该框架理论由一个波传播/化学反应相互作用机制,两个基本物理过程(热点和化学反应带加速),三个临界状态(临界传播状态、临界起爆状态和临界胞格尺度)等六个关键要素构成。应用五个典型物理算例,界定了这些关键要素的物理机制、表现特性及其客观存在性。应用这个爆轰波统一框架理论,本文成功地解释了经典爆轰理论、应用CFD技术获得的多维爆轰波结果和实验观察到的胞格爆轰现象依据的关键要素。虽然爆轰波的三个临界状态需要进一步的定量化研究,但是应用框架理论解释已经获得的爆轰波研究结果表现的统一性及其对深入研究爆轰现象的预测性表明了该框架理论对于开展深入的爆轰物理具有重要意义。

污染空气对乙烯燃料超燃冲压发动机燃烧室流场影响的数值评估

由于纯净空气的超燃燃烧室试验条件比较缺乏，污染空气效应对试验结果的影响往往通过CFD计算来评估。本文通过数值模拟中国科学院力学研究所高温气动实验室的乙烯燃料超燃燃烧室试验，来数值评估来流空气污染效应对试验结果的影响。分别计算了燃烧室入口来流马赫数1.8和2.5的工况，计算结果表明：马赫1.8情况下，空气的污染成分有利于燃烧，纯净空气情况下由于凹腔内处于富油情况燃烧增压不够；马赫2.5情况下，空气的污染成分对燃烧影响不大。

超声速平板绕流的全流域阻力特性研究

超声速流动中摩擦阻力的准确预测一直是气动研究中的难点。本文以零攻角的平板绕流为对象，结合传统的CFD方法求解Navier-Stokes方程和直接模拟Monte-Carlo（DSMC）方法模拟粒子的微观运动，研究平板在超声速来流下的全流域阻力特性。研究发现，在超声速流中会产生激波-边界层相互干扰作用并随着来流马赫数的增加而加强，同时平板前缘存在非平衡效应显著的区域；努森数的增大使激波和边界层均增厚，粘性作用增强，非平衡区的范围增大，非平衡强度增强，而激波强度减弱；最后通过拟合连续流条件下的摩阻系数和修正自由分子流结果，结合桥函数，获得了全流域范围内均有效的平板阻力系数表达式。

交变流动中压力速度之间相位差对突变截面损失的影响分析

从弱非线性热声理论出发,给出交变流动中突变截面阻力系数的定义以及考察方法。通过PIV(粒子成像测速仪)测量与CFD计算结果的对比,验证了CFD计算结果的可靠性。进而利用CFD模拟考察了交变流动中压力与速度之间相位差对突变截面局部损失的影响,观察到产生这种影响的流场内部流动机制。分析表明,阻力特性确实对声场压力与速度之间的相位差存在依赖性,但这种依赖性会随着非线性的增强而减少。

Hydroelastic dynamic characteristics of a slender axis-symmetric body

The slender axis-symmetric submarine body moving in the vertical plane is the object of our investigation. A coupling model is developed where displacements of a solid body as a Euler beam (consisting of rigid motions and elastic deformations) and fluid pressures are employed as basic independent variables, including the interaction between hydrodynamic forces and structure dynamic forces. Firstly the hydrodynamic forces, depending on and conversely influencing body motions, are taken into account as the governing equations. The expressions of fluid pressure are derived based on the potential theory. The characteristics of fluid pressure, including its components, distribution and effect on structure dynamics, are analyzed. Then the coupling model is solved numerically by means of a finite element method (FEM). This avoids the complicacy, combining CFD (fluid) and FEM (structure), of direct numerical simulation, and allows the body with a non-strict ideal shape so as to be more suitable for practical engineering. An illustrative example is given in which the hydroelastic dynamic characteristics, natural frequencies and modes of a submarine body are analyzed and compared with experimental results. Satisfactory agreement is observed and the model presented in this paper is shown to be valid.

采用改进尾流振子模型的柔性海洋立管的涡激振动响应分析

建立了柔性杆件在非均匀流作用下的涡激振动响应预测模型, 考虑了涡激振动锁频阶段流体附加质量的变化, 以及振动响应和来流简缩速度的非线性关系. 该模型通过经验公式结合迭代求解的方式, 计算方便、速度快, 避免了数值计算(CFD)的繁琐, 较为适合于海洋工程实际应用. 与试验和数值结果的比较表明采用该文提出的计算模型, 可以更合理、准确地给出结构涡激振动响应. 最后, 结合实际平台参数, 进行了柔性立管在非均匀流场的作用下的涡激振动响应分析, 并研究了立管的预张力、流场分布等参数的影响. 分析结果表明: 随着立管张力和流场分布的改变, 各阶模态锁频区域发生了变化, 从而改变了结构的总体响应