强旋湍流气－固两相流动的颗粒随机轨道法模拟

应用颗粒随机轨道模型，并与一种新的代数Ｒｅｙｎｏｌｄｓ应力模型相结合，对新型煤粉涡旋燃烧炉内强旋湍流气。固两相流动进行了数值模拟。得到了与实验相符合的颗粒相密度分布和质量流分布。计算结果表明，在涡旋燃烧炉内的强旋湍流流场中，外壁附近颗粒浓度最高，颗粒停留时间加长，气－固两相间滑移速度增大。     

强旋湍流气－固两相流动的颗粒随机轨道法模拟

应用颗粒随机轨道模型，并与一种新的代数Ｒｅｙｎｏｌｄｓ应力模型相结合，对新型煤粉涡旋燃烧炉内强旋湍流气。固两相流动进行了数值模拟。得到了与实验相符合的颗粒相密度分布和质量流分布。计算结果表明，在涡旋燃烧炉内的强旋湍流流场中，外壁附近颗粒浓度最高，颗粒停留时间加长，气－固两相间滑移速度增大。     

稠密气固两相湍流流动的实验和数值模拟

基于气固两相流动模型计算循环流化床内稠密气固两相流湍流动,颗粒动理学方法模拟颗粒相湍动能,ＳＧＳ模型模拟气相湍流,采用γ－射线密度计和非等速取样管测量局部颗粒浓度和流率,利用ＦＦＴ方法计算颗粒浓度功率谱密度。模拟计算得到上升管内气相和固相速度和浓度分布等。同时数值模拟与Ｔｓｕｊｉ等和Ｋｎｏｗｌｔｏｎ等试验结果进行了比较,结果表明数值模拟计算与实验结果相吻合。     

气体—颗粒激波松弛流动特性研究

提出了一种解气体－颗粒粘性流动方程组的方法。在此方法中，计算颗粒阻力时采用Ｈｅｎｄｅｒｓｏｎ阻力系数公式，得到两相激波松弛流动数值结果，并与实验结果作了比较，详细分析和讨论了计算结果，揭示了两相激波松弛流动特性随颗粒容积比分的变化规律     

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

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

透平叶栅跨音速流动计算中的新型有限元法

将Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法和多级有限元的思想结合起来，构成了在收敛速度和稳定性两方面均较好的新型有限元算法：多级广义有限元。利用这一方法，分别基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和Ｅｕｌｅｒ方程，研究了透平跨音速叶栅无粘流动和粘性流动，并将计算结果与实验结果作了比较。计算结果表明，本方法是透平机械内部跨音速流动计算的强有力的手段。     

管内湍流旋流的数值计算

本文用有限差分法对直管内的湍流旋流进行了数值模拟。计算中采用Ｂｏｕｓｓｉｎｅｓｑ湍流涡粘性假设的基本思想和Ｋ－ε双方程模型来求解雷诺应力各分量。为了反映旋流中湍流转输的非均匀性和各向异性特征，对雷诺应力各分量及与之相主尖的各湍流粘性系数分别进行计算。计算结果表明该模型能较好地反映直管内湍流旋流的流动结构。     

激波点燃粉尘的数值模拟研究

建立了用于模拟入射激波后可燃粉尘颗粒点火的一维非定常两相化学反应流模型，该模型考虑了气固两相间的相互作用、粉尘颗粒的加速、加热和化学反应。粉尘颗粒着火前的化学反应用发生在颗粒外表面和内孔表面的非均相反应描述，颗粒内部的温度变化用一含有化学反应源项的非稳态热传导方程来描述，以颗粒外表面温度的突跃上升作为可燃粉尘颗粒点燃的着火条件。我们用该模型和ＰＳＩＣ方法，对由中等强度激波从纯气相传入煤粉－氧气混合物而引起的非定常两相流动现象，包括气固两相间的相互作用、粉尘颗粒的加速、加热以及点火过程进行了数值研究，计算了对应于不同载荷比、马赫数为４～５的入射激波后煤尘颗粒的点火延迟时间，分析了由于可燃粉尘颗粒的存在，入射激波及波后气固两相流动参数的变化规律。数值计算结果与实验数据符合较好。文中建立的模型和所用的基于ＰＳＩＣ算法的数值方法，用最自然的方式描述气固两相流动，即用连续流模型（欧拉方程）描述输运相（气相）的流动，用轨道颗粒模型（拉格朗日方程）描述分散相（颗粒相）的运动。用这种方法模拟含尘介质中激波后颗粒的点火是很有效的，它可以清楚地确定哪一个颗粒群最先着火，它的初始位置以及在整个点火延迟时间内     

统一二阶矩模型用于模拟旋流湍流两相流动

用统一二阶矩模型（ＵＳＭ）模拟了旋流数为０４７和１５的气粒两相流动，并和实验结果以及ｋ ε ｋｐ模型的模拟结果进行了对比．研究结果表明，提高旋流数减小了轴向速度反流区，增大了切向速度似固核区．ＵＳＭ和ｋ ε ｋｐ模型预报旋流数为０４７时的两相速度场差别不大，并都和实验结果接近，但前者预报的旋流数为１５的两相速度场比后者有改进，在两种情况下，前者都能揭示出后者无法预报的两相湍流各向异性规律．     

平面突扩流动非稳定性的大涡模拟

用ＴＧ有限元法求解二维大涡模拟的非定常Ｎ－Ｓ方程，数值模拟了平面突扩流动，计算结果展示了回流区内流动的多涡结构及其不稳定的周期振荡过程，而该振荡过程的周期时间历程平均，恰好与诸多湍流模型的计算结果及实验结果相吻合，证明了数值模拟的正确性。     

Numerical simulation of incompressible two‐phase flows with a Boussinesq–Scriven interface stress tensor

We consider the numerical simulation of a three‐dimensional two‐phase incompressible flow with a viscous interface. The simulation is based on a sharp interface Navier–Stokes model and the Boussinesq–Scriven constitutive law for the interface viscous stress tensor. In the recent paper [Soft Matter 7, 7797–7804, 2011], a model problem with a spherical droplet in a Stokes Poiseuille flow with a Boussinesq–Scriven law for the surface viscosity has been analyzed. In that paper, relations for the droplet migration velocity are derived. We relate the results obtained with our numerical solver for the two‐phase Navier–Stokes model to these theoretical relations. Copyright © 2013 John Wiley & Sons, Ltd.     

明渠气二相流的双流体模型

从解决平均二相流基本方程的封闭问题出发,分析了水利工程中气水二相流的特点,介绍了据此提出的可用于水利工程稀疏气泡流计算的双流体模型,并对明渠气水二相流进行了算例验证。     

水平管内油——水两相流动压降规律的实验研究

设计和建造了内径为26．1mm,长30m的水平不锈钢多相流实验环道,利用白油与水进行了油－水两相流流型和压降实验。本文针对各种流型,分析了油－水两相流动的压降规律和油－水混合液有效粘度,指出有效粘度法只适用于油－水分散流型的压降预测,对于分层流型式其它混合流型使用合适该种流型的压降预测模型来计算压降。研究结论对油田现场的油水混输管路的经济运行具有较大的指导意义。     

单裂缝中携砂液流动规律研究

裂缝中携砂液流动是一种固液两相流,携砂液的运移与支撑剂的铺置是水力压裂裂缝保持导流能力的关键.本文基于FLUENT流体计算软件,采用双流体模型,将颗粒看作拟流体,携砂液按照牛顿流体处理,分析了支撑剂体积分数α_s、阿基米德数Ar、颗粒雷诺数Re以及裂缝入口边界对流动规律的影响.研究结果表明:携砂液在裂缝中的流动过程中,发展成为支撑剂体积分数不同的四个区域,包括砂堤区、颗粒悬浮区、颗粒滚流区和无砂区;支撑剂的沉降程度随着支撑剂体积分数和阿基米德数的增加而增加,而随着雷诺数增加而降低;入口为网眼型时,进入裂缝后过流面积的增加导致流速突降,使得支撑剂更容易在入口处产生堆积,在同一入口流速下,较均匀入口的工况铺砂高度大.     

Gas-particle two-phase turbulent flow in a vertical duct

Two-phase gas-phase turbulent flows at various loadings between the two vertical parallel plates are analyzed. A thermodynamically consistent turbulent two-phase flow model that accounts for the phase fluctuation energy transport and interaction is used. The governing equation of the gas-phase is upgraded to a two-equation low Reynolds number turbulence closure model that can be integrated directly to the wall. A no-slip boundary condition for the gas-phase and slip-boundary condition for the particulate phase are used. The computational model is first applied to dilute gas-particle turbulent flow between two parallel vertical walls. The predicted mean velocity and turbulence intensity profiles are compared with the experimental data of Tsuji et al. (1984) for vertical pipe flows, and good agreement is observed. Examples of additional flow properties such as the phasic fluctuation energy, phasic fluctuation energy production and dissipation, as well as interaction momentum and energy supply terms are also presented and discussed.

Applications to the relatively dense gas-particle turbulent flows in a vertical channel are also studied. The model predictions are compared with the experimental data of Miller & Gidaspow and reasonable agreement is observed. It is shown that flow behavior is strongly affected by the phasic fluctuation energy, and the momentum and energy transfer between the particulate and the fluid constituents.     

The flow structure of dilute gas-particle suspensions behind a shock wave moving along a flat surface

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented.The Saffman lift force acting on a particle in a shear flow istaken into account.It is shown that particle migration across the boundary layer leads tointersections of particle trajectories.The corresponding modification of dusty gas model isproposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at aconstant speed are obtained by using the method of matched asymptotic expansions.Themethod of the calculation of particle phase parameters in Lagrangian coordinates isdescribed in detail.Some numerical results for the case of small particle concentration aregiven.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd.     

A THREE-FLUID MODEL OF THE SAND-DRIVEN FLOW

ＡＴＨＲＥＥ－ＦＬＵＩＤＭＯＤＥＬＯＦＴＨＥＳＡＮＤ－ＤＲＩＶＥＮＦＬＯＷ￥(刘大有，董飞)ＬｉｕＤａｙｏｕ；ＤｏｎｇＦｅｉ（ＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＡｃａｄｅｍｉａＳｉｎｉｃａ，Ｂｅｉｊｉｎｇ１０００８０，Ｐ．Ｒ．Ｃｈｉｎａ）Ａｂｓ...     

基于结构两相阻力分析的泥石流堆积形态研究

针对泥石流的结构两相阻力特征,采用结构两相阻力对泥石流的堆积形态进行了静力学分析,推导出在一般条件下泥石流堆积形态的公式,并讨论了在不同条件下泥石流堆积形态的特点。分析表明采用结构两相阻力对泥石流堆积形态分析可以较合理的模化阻力项,能够表现泥石流体复杂的成分结构组成,较真实的反映实际泥石流的堆积形态特征。在特定材料条件下,结构两相阻力项转化为单一材料阻力项,所得到的相应结论与已有对单一材料的分析结论一致。     

Modeling of unidirectional blood flow in microvessels with effects of shear-induced dispersion and particle migration

A cell-free layer, adjacent to microvessel walls, is present in the blood flow in the microcirculation regime. This layer is of vital importance for the transport of oxygen-saturated red cells to unsaturated tissues. In this work, we first discuss the physics of formation of this cell-free layer in terms of a balance between the shear-induced dispersion and particle migration. To this end, we use high-viscosity drops as prototypes for cells, and discuss our results in terms of physical parameters such as the viscosity ratio and the capillary number. We also provide a short-time analysis of the transient drift-dispersion equation, which helps us better explain the formation process of the cell-free layer. Moreover, we present models for investigating the blood flow in two different scales of microcirculation. For investigating the blood flow in venules and arterioles, we consider a continuous core-flow model, where the core-flow solution is considered to be a Casson fluid, surrounded by a small annular gap of Newtonian plasma, corresponding to the cell-free layer. We also propose a simple model for smaller vessels, such as capillaries, whose diameters are of a few micrometers. In this lower-bound limit, we consider a periodic configuration of aligned, rigid, and axi-symmetric cells, moving in a Newtonian fluid. In this regime, we approximate the fluid flow using the lubrication theory. The intrinsic viscosity of the blood is theoretically predicted, for both the lower and upper-bound regimes, as a function of the non-dimensional vessel diameter, in good agreement with the previous experimental works. We compare our theoretical predictions with the experimental data, and obtain qualitatively good agreement with the well-known Fåhræus-Lindqvist effect. A possible application of this work could be in illness diagnosis by evaluating changes in the intrinsic viscosity due to blood abnormalities.