密相液固两相湍流流动研究

在ＩＰＳＡＲ算法基础上推导了适用于考虑液度变化影响的密相液固两相湍流流动数值计算的ＤＩＰＳＡＲ算法，并采用低雷诺数模型，对竖直上升管中密相液固两相湍流分别采用ＩＰＳＡＲ算法和ＤＩＰＳＡＲ算法进行了数值计算，计算值与实测值较符合，计算结果的比较表明ＤＩＰＳＡＲ算法能更有效地预测密相液固两相湍流流动     

离心泵叶轮内密相液固两相湍流的数值模拟

利用作者建立的描述密相液固两相湍流的 KET模型和推导的基本控制方程组 ,在处理壁面边界条件时考虑了颗粒和叶片的相互碰撞作用 ,对离心泵叶轮内密相液固两相流动进行了数值模拟 ,得到了泵叶轮内两相流动的一些规律 ,为液固两相流泵的设计提供了一定的理论依据。     

高浓度固液两相流压力测量

本文介绍了确定水煤浆流动阻力的试验装置,应用膜片或传压器的传感器,可用来测定高浓度固液两相流压力。     

水平液─固流中颗粒抑制湍流的行为和条件

用激光多普勒分相测量方法，考察了水平含颗粒水流的主流区，亚毫米颗粒对水流湍流有抑制行为，湍流削弱程度与颗粒浓度、流动发展长度、水流雷诺数等因素有关．按照不出现尾涡脱落和至少能响应一部分大尺度涡的运动这两个基本条件，分析了湍流抑制与有关参量的依赖性．结果表明，流动尺度和颗粒尺寸之比，对于颗粒消耗而不增生湍流起着最主要的作用．在一定的水流雷诺数下，只要这一比值很大，亚毫米颗粒就会表现出对湍流的抑制作用     

气液两相瞬变流的流固耦合模型研究

传统的气液两相流瞬态分析和管道动力响应计算是分开的，存在一定的缺陷．针对石油工业中常见的多相混输问题，介绍了常见的气液两相流瞬态模型和流固之间存在的耦合机理，在不作薄壁管假设的前提下推导出了气液两相瞬变流的流固耦合模型．与现有相似模型的对比分析表明，这一模型比较全面地考虑了流体和管道的特性以及不同的耦合形式，可以适应实验和／或仿真研究的需要．     

固-液两相流黑水管道冲蚀磨损的数值模拟研究

煤气化黑水处理系统管道由于其流体介质高含固体颗粒和腐蚀性介质,且工作在高温、高压差环境中,极易受到冲蚀磨损和腐蚀的耦合作用而失效,影响其服役寿命.采用计算流体力学(CFD)方法数值模拟研究了煤气化黑水处理系统固-液两相流管道的冲蚀磨损行为和机理,以及流体介质速度和固体颗粒粒径对管道冲蚀磨损的影响规律,并分析了盲通管和涡室结构对弯管冲蚀磨损行为的优化改善效果.研究结果显示,煤气化黑水处理系统管线的冲蚀高危区主要分布在弯管外拱和变径管等结构突变区域;管道冲蚀磨损行为与其内部流体的运动和颗粒冲击特性有关;管道的冲蚀率均随着流体速度的增加而加剧,而粒径对弯管和变径管冲蚀率的影响并非单调关系,这与颗粒受力作用有关;弯管优化分析显示,涡室结构可以降低弯管的最大冲蚀率,减缓弯管的冲蚀磨损.     

瞬态通用固液两相流模型及其在动态井眼清洁模拟中的应用

为准确了解岩屑在井眼中的实时运移、积累过程及其变化规律,基于质量守恒和动量守恒方程,考虑岩屑运移流型和流型转换条件,建立了适应流型变化的瞬态通用固液两相流模型;并将其应用于一口窄密度窗口井的动态井眼清洁模拟过程中,实时模拟整个钻井过程中岩屑沿井筒分布情况.研究结果表明:岩屑对环空压力分布有重要影响;成功钻穿窄安全密度窗...     

两种颗粒湍流扩散模型数值模拟气液两相流泄漏扩散的比较

采用相间耦合的欧拉-拉格朗日方法,模拟了装有液化气（丙烷）的容器出现小孔或裂缝时,发生泄漏后的气液两相扩散过程。分别应用随机轨道模型及颗粒群模型来考察湍流对液滴扩散的影响,并与确定性轨道进行了比较。给出了数学物理模型,计算结果与实验数据进行了对比。结果表明：采用随机轨道模型能较好地描述液滴的湍流扩散,适用于有液相蒸发的两相流扩散问题。     

旋转直管内气固或液固两相流动浓度分布

本文针对旋转直管内气固或液固两相流动，建立了经过适当简经处理的颗粒无量纲运动方程，得到了稳定状态下直管内颗粒的浓度分布公式，从而为气固两相流风机或液固两相流泵的颗粒浓度分布研究打下了基础。     

沙粒粒径在水平水流中对湍流调制及减阻效应的影响作用

运用基于欧拉-欧拉方法的混合欧拉多相流模型结合雷诺应力湍流模型,对较大雷诺数下的水平固液两相湍流进行了数值计算,主要考察了以Kolmogorov尺度(约为2v/ur)为临界值的细小沙粒颗粒在5％的固相体积分数下对湍流场的调制及其减阻效应.研究发现,无量纲颗粒直径dp+≤2的小颗粒减小了固液两相湍流的雷诺应力,并且三个方向上的速度脉动也被不同程度地削弱;而dp+=4的大颗粒使缓冲层区域的雷诺应力稍增大,在增强法向速度脉动的同时对流向脉动有抑制作用,并且值得关注的是,较大颗粒的存在导致缓冲层中的部分区域出现了流变现象.在减阻方面,较小的颗粒(dp+≤2)有大致相同的减阻表现,而大颗粒(dp+=4)已经失去了减阻性能,总体上看,在所研究的情形下微颗粒的减阻性能随着其粒径的增大而降低.     

Particle/Turbulence Interactions, Mass Transfer and Gas/Solid Chemistry in a CFBC Riser

Gas/solid chemistry in the upper, dilute region of a circulating fluidised bed combustor (CFBC) riser is addressed. The limitations of turbulent mixing are illustrated by the example of the heterogeneous NO/CO/char reaction, relevant in CFB combustion of nitrogen-containing solid fuels. The mass transfer of the gaseous reactants to the char surface is determined, and how the conversion is influenced by the degree of mixing of the multiphase system by turbulent dispersion. Particle/turbulence interactions are taken into account by a (Lagrangian) frequency spectrum of velocity fluctuations, which determine the dispersion of the char particles described here with the Tchen–Hinze model. Chars from solid fuels characterised by fuel ratio (FR) ranging from 0.1 (wood) via 0.5 (peat) and 1 (coal A) to 2 (coal B) were considered. The effective rate of the NO/CO/char reaction is determined as a function of the size and type of the char particle, temperature, particle concentration, reactor dimensions and fluidization velocity, at atmospheric pressure. It was found that for this case the effective gas/solid conversion rate in the upper, dilute region of the riser is much lower than the gas/solid chemistry, mainly due to mass transfer limitations for char particles with sizes of typically 300 μm. The concentration of NO at the char particle surface is only a few % of that in the bulk gas phase. Strong influences were found for particle size and temperature, whilst the fluidisation velocity and the reactor size have only a small influence. It is concluded that for a typical CFB riser, for particles larger than approx. 20 μm, mass transfer has a stronger influence on the heterogeneous NO/CO/char reduction mechanism than the ``unmixedness' due to particle eddy dispersion limitations. It is recommended that this or a similar approach to turbulent dispersive mixing is implemented in CFD codes when these are used for boiler and furnace calculations. This revised version was published online in July 2006 with corrections to the Cover Date.     

一种紊流润滑理论分析新方法—复合型紊流模式理论

在对主要用于高速轻载工况的常用紊理论进行简要分析和充分考察润滑流场的边界条件及内部结构的基础上，采用理论上比现有紊流润滑理论更为合理的复合型紊流模式理论，即在近壁区采用低紊流雷诺数的ｋ－ε模式，而在紊流核心区采用代数雷诺应力模式，对复杂流场的紊流润滑进行了分析，同时在计入惯性效应的情况下，推导出了一种适用于高压密封和高速重载轴承等计算，计算结果与实验数据十分吻合，验证了模型的有效性，可以应用于高压     

Anisotropy Invariant Reynolds Stress Model of Turbulence (AIRSM) and its Application to Attached and Separated Wall-Bounded Flows

Numerical predictions with a differential Reynolds stress closure, which in its original formulation explicitly takes into account possible states of turbulence on the anisotropy-invariant map, are presented. Thus the influence of anisotropy of turbulence on the modeled terms in the governing equations for the Reynolds stresses is accounted for directly. The anisotropy invariant Reynolds stress model (AIRSM) is implemented and validated in different finite-volume codes. The standard wall-function approach is employed as initial step in order to predict simple and complex wall-bounded flows undergoing large separation. Despite the use of simple wall functions, the model performed satisfactory in predicting these flows. The predictions of the AIRSM were also compared with existing Reynolds stress models and it was found that the present model results in improved convergence compared with other models. Numerical issues involved in the implementation and application of the model are also addressed.     

AN EXTENDED k-ε MODEL FOR NUMERICAL SIMULATION OF WIND FLOW AROUND BUILDINGS

It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold’s normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.     

Computation of backward-facing step flows by a second-order Reynolds stress closure model

This paper scrutinizes the predictive ability of the differential stress equation model in complex shear flows. Two backward-facing step flows with different expansion ratios are solved by the LRR turbulence model with an anisotropic dissipation model and the near-wall regions of the separated side resolved by a near-wall model. The computer code developed for solving the transport equations is based on the finite-volume-finite-difference method. In the numerical solution of the time-averaged momenum equations the Reynolds stresses are treated partially as a diffusion term and partially as a source term to avoid numerical instability. Computational results are compared with experimental data. It is found that the near-wall region of the separated side resolved by the near-wall model, the LRR model with a simple modification of an anisotropic dissipation model can predict backward step flows well.     

固—液耦合Timoshenko管道的稳定性分析

根据Ｈａｍｉｌｔｏｎ原理的固－液耦合振同分方程用幂级数法计算了Ｔｉｍｏｓｈｅｎｋｏ管道的固有频率和临界流速。给出了管道前三阶固有频率－流速的关系曲线,分析了转动惯量对该输流管道的稳定特性的影响。计算结果表明,转动一对两端简支的固－液合Ｔｉｍｏｓｈｅｎｋｏ管道的静力失稳没有影响,但对其频率特性和动力失稳有影响。     

航空发动机轴承腔润滑的气液两相均匀流研究

基于轴承腔中润滑油的两相均匀流动模型,采用湍流模式和有限差分数值方法计算轴承腔内三维定常N-S方程,对腔内润滑油的气液两相均匀流动特性进行研究,以获得气液两相均匀流条件下润滑油流场、压力场和速度场在轴承腔内的分布情况,分析转子转速、含气率和润滑油进口速度对润滑油出口压力以及润滑油与壁面之间剪切力的影响,同时对单相流和两相均匀流润滑性能差异进行比较.结果表明,转子转速、含气率和润滑油进口速度对润滑油出口压力和腔内壁面与润滑油间的平均剪切力具有不同影响,而采用2种流动模型计算出的轴承腔润滑油出口压力的差异较大,同时支持了开展航空发动机轴承腔润滑两相流动分析的必要性.     

壁湍流猝发过程中速度分量的相位差对雷诺应力影响的实验研究

用IFA300恒温热线风速仪和×形二分量热线探针,以采样间隔小于最小湍流时间尺度的分辨率,精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、展向速度分量的时间序列信号.用子波分析辨识壁湍流相干结构猝发事件的能量最大准则,确定壁湍流相干结构猝发事件的时间尺度;用条件相位平均技术提取了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形,用互相关方法研究了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形的相位差关系及其对雷诺应力的影响,发现在缓冲层和对数律区,展向脉动速度与流向脉动速度的条件相位平均波形具有不同的相位;当两者相位基本一致时,雷诺应力达到正的最大值,此时湍流相干结构的产生非常活跃;当两者相位差分别集中在90°和270°附近时,雷诺应力的幅值减小并接近于零,此时湍流相干结构的产生和猝发都得到了抑制.     

Advanced Turbulence Modelling of Separated Flow in a Diffuser

The paper describes an investigation into the predictive performance of linear and non-linear eddy-viscosity models and differential stress-transport closures for separated flow in a nominally two-dimensional, asymmetric diffuser. The test case forms part of a broader collaborative exercise between academic and industrial partners. It is demonstrated that advanced turbulence models using strain-dependent coefficients and anisotropy-resolving closure offer tangible advantages in predictive capability, although the quality of their performance can vary significantly, depending on the details of closure approximations adopted. Certain features of the flow defy resolution by any of the closures investigated. In particular, no model resolves correctly the flow near the diffuser's inclined wall immediately downstream of the inlet corner, which may reflect the presence of a “flapping” motion associated with a highly-localised process of unsteady separation and reattachment. This revised version was published online in July 2006 with corrections to the Cover Date.     

两相流显微PIV/PTV系统的开发

开发了一个能同时测量两相流中两相速度和细颗粒尺寸分布的显微PIV/PTV系统,其硬件系统包括大功率连续激光器、显微镜、高速摄像机;软件系统由改进的球形颗粒图像识别算法、各种图像处理算法和各种先进的PIV/PTV算法组成。其中改进的圆弧识别算法能够进行更准确地进行曲线分割而能对充满噪音并相互重叠的颗粒图像给出较好的识别结果。应用该PIV系统,可以在微秒和微米数量级上捕获细颗粒/气泡图像,并能较准确地同时得到两相速度、颗粒尺寸和浓度分布。对焚香可吸入颗粒物进行了速度和尺寸的同时测量,得到了较满意的结果。