可压缩气固混合层中离散相与连续相的相互作用研究

尽管已有许多文献采用数值模拟方法研究两相流问题,但主要是集中不可压流动方面.本文采用Eul-er-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流动进行了研究.气相流场采用非定常全Navier-Stokes方程描述,并应用具有空间三阶精度的WNND(Weighted Non-Oscillatory, Contai-ning No Free Parameters and Dissipative)格式进行数值高散.固相方程采用二阶单边三点差分离散.在考虑流场对固粒作用的同时,也计及颗粒对流场的反作用.主要研究混合层大尺度涡对颗粒扩散特性的影响及颗粒对流场结构的影响问题.在对流马赫数为0.5时,研究不同Stokes数颗粒在连续流场中的扩散特性,而在对流马赫数为0.8时研究了不同Stokes数颗粒对流场小激波结构的影响.     

气固两相混合层二维涡配对的数值研究

采用双向耦合模型中的速度耦合模型,数值模拟了气固两相混合层中涡的配对、合并过程,文中采用拟谱方法对流场进行了直接数值模拟,用Lagrange模型跟踪固粒,结果发现,在双向耦合过程中Stokes数仍然是重要的参数,但不是唯一影响流场的参数,流场的发展还与固粒的体积浓度、固粒的相对密度以及固粒大小等因素有关。固粒抑制流场的发展,阻碍涡的配对、合并,加快了涡量的扩散;小St数的固粒仍然跟随流体运动,大St数的固粒趋向于低涡量区的趋势减弱。     

可压缩混合层中细光束畸变效应的实验研究

利用双曝光CCD相机成像技术,对小直径光束穿越可压缩混合层流场存在的大尺度结构特性及其引起的气动光学效应进行了实验研究。实验研究结果表明:在不同对流马赫数(Mc=0.17,0.45)下,可压缩混合层中的大尺度结构形状不同,因而具有不同的折射率分布场,从而导致光束投影发生了不同程度的变形和偏移。利用光束质心的概念定量分析了不同对流马赫数下光束偏移和抖动的统计特性。统计结果表明在高对流马赫数下光束投影的偏移量减小,光束沿流场展向的平均偏移角度减小,但摆动幅度加大。     

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

悬浮固粒分区存在对射流稳定性的影响

本文利用气固两相耦合模型,理论推导出大雷诺数下悬浮固粒分区存在的射流稳定性方程,通过数值计算得到了固粒存在不同流域时流场扰动放大因子曲线和波速曲线。进而,在分析所得稳定性曲线的基础上,得到了关于固粒充满全流场、存在于射流区及存在于非射流区对流场中扰动增长及扰动传播影响的重要结论。这些结论对气固两相射流场的发展的认识和工程实际中实施对气固两相射流场的人工控制有指导意义。     

可压缩混合层流动近十年研究进展

从稳定性分析、实验研究和数值模拟3个方面回顾了过去近10年间关于可压缩混合层流动的研究.这些工作研究了可压缩性对混合层其它性质的影响,这些性质包括混合层增长率、不稳定模式、流场结构、平均速度和脉动速度、以及标志物的混合等.3种研究方法相互补充, 共同提供了认识可压缩混合层流动机制的大量信息.这些研究是湍流机理研究的重要组成部分,也是混合层流动在应用技术领域成功应用的前提.      

气固两相混合层固粒分布特性的分数维方法描述

本文对气固两相混合层中固粒的运动进行了数值模拟，得到了不同Ｓｔ数下，涡卷起和配对过程中固粒的分布特性，计算了不同时刻固粒分布的相关维数，说明分数维是一个可用来定量描述固粒分布均匀程度的重要参数，Ｓｔ数为１０的数量级时，对应最小的相关维数，固粒分布最不均匀，研究结果对控制固粒的均匀程度具有指导意义。     

一种减轻固粒对壁面冲蚀磨损的新方法

针对带固粒的近壁流场，提出了在壁面上开凿纵向浅沟槽以减轻固粒对壁面冲蚀磨损的方法．采用改进的κ-ε湍流模式和气-固两相双向耦合模型，计算了流场中颗粒的速度、轨迹以及固粒对壁面的冲蚀磨损，并进行了相应的试验对比．结果表明：在壁面上开凿纵向沟槽能减轻固粒对壁面的磨损；在一定的沟槽高度下，当沟槽宽度同沟槽问距离相等时，减轻磨损的效果最明显；在一定颗粒尺寸范围内，大颗粒导致磨损质量损失增大，当颗粒尺寸超过某一临界值后，颗粒尺寸对磨损的影响较小．进口处颗粒运动方向对壁面磨损具有直接影响．     

含悬浮固粒的旋转射流剪切层稳定性研究

本文在理想不可压旋转圆射流的运动方程中添加了固粒作用项，由此推得了时间增长率的表达式，进而得关于含悬浮固粒放置射流稳定性的修正瑞利稳定性准则，求出了不同固粒质量密度固－气脉动速度比值，固气脉动速度相位差及Ｓｔｏｋｅｓ数情况下旋转射流场的增长率与径向空间波数的关系曲线，在比较这些曲线的基础上，给出了关于固粒属性对旋转射流场稳定性影响的几个重要结论为控制旋转射流场和后续发展提供了依据。     

自由剪切湍流中颗粒-拟序结构相互作用研究进展

从实验和数值模拟两方面评述了颗粒-湍流拟序结构相互作用的近期研究进展.关于Stokes数不同对颗粒行为和拟序结构影响的试验研究,从单点激光多普勒测量到粒子图像全场测速,并与流场显示定性方法结合,揭示了不同Stokes数范围颗粒-拟序结构相互作用的规律.基于涡方法、直接数值模拟和大涡模拟等的模拟研究,进一步揭示颗粒-拟序结构的相互耦合作用和外界激励的调制作用等规律,同时推动了算法的发展.      

Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Research on the effect of particle of two-dimensional shear flow

ＩｎｔｒｏｄｕｃｔｉｏｎＴｕｒｂｕｌｅｎｔｆｌｏｗｌａｄｅｎｗｉｔｈｐａｒｔｉｃｌｅｓｉｓｃｏｍｍｏｎｉｎｂｏｔｈｎａｔｕｒｅａｎｄｉｎｄｕｓｔｒｙ .Ｔｈｅｍｉｘｉｎｇｄｅｇｒｅｅｏｆｐａｒｔｉｃｌｅｓｉｎｔｈｅｇａｓ＿ｓｏｌｉｄｆｌｏｗｉｓｓｉｇｎｉｆｉｃａｎｔｔｏｉｎｄｕｓｔｒｉａｌａｐｐｌｉｃａｔｉｏｎｓ.Ｔｈｅｒｅｃｅｎｔｒｅｓｅａｒｃｈｅｓｉｎｄｉｃａｔｅｔｈａｔｔｈｅｃｏｈｅｒｅｎｔｓｔｒｕｃｔｕｒｅｓｅｘｉｓｔｉｎｔｈｅｍｉｘｉｎｇｌａｙｅｒ.Ｃｒｏｗｅｅｔａｌ.[1]ｅｖａｌｕａｔｅｄ…     

双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况,本文以颗粒长径比为2,颗粒之间的初始距离ΔS_Py=4D为例,基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布,在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下,研究了颗粒的受力和运动的变化,探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明,双柱状颗粒在线性剪切流场中易向速度大的流体区域运动;前后两颗粒运动状态和轨迹不同,颗粒之间距离较近时,曳力会产生较大的波动;只有当颗粒在壁面附近时,滞后颗粒才能追上领先颗粒,两颗粒发生牵引、翻滚和分离过程。     

Effect of particle and fiber size on the morphology of deposits in fibrous filters

The capture of aerosol particles during their production and purification of air at workplace and atmospheric environment require an efficient separation method of particulate matter from the carrier gas. Many papers have been published recently on the deposition of particles on fibrous collectors, often using the classical continuum approach to describe the process. However, the approach is not convenient for studying the influence of particle deposition on filter performance (filtration efficiency and pressure drop) when one has to introduce nonsteady‐state boundary conditions. For the purpose of this work, the lattice Boltzmann model describes fluid dynamics, while the solid particle motion is modeled by the Brownian dynamics. The aim of this study is to model the influence of single particle loading (amount and morphology – e.g., porosity and fractal dimension) on its performance. Copyright © 2014 John Wiley & Sons, Ltd.     

垂直湍流液-固流中大颗粒的相对速度

通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式．用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数．实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。     

DRAG FORCE IN DENSE GAS—PARTICLE TWO—PHASE FLOW

Numerical simulations of flow over a stationary particle in a dense gas-particle two-phase flow have been carried out for small Reynolds numbers (less than 100). In order to study the influence of the particles interaction on the drag force, three particle arrangements have been tested: a single particle, two particles placed in the flow direction and many particles located regularly in the flow field. The Navier-Stokes equations are discretized in the three-dimensional space using finite volume method. For the first and second cases, the numerical results agree reasonably well with the data in literature. For the third case, i.e., the multiparticle case, the influence of the particle volume fraction and Reynolds numbers on the drag force has been investigated. The results show that the computational values of the drag ratio agree approximately with the published results at higher Reynolds numbers (from 34.2 to 68.4), but there is a large difference between them at small Reynolds numbers. The project supported by the Special Funds for Major State Basis Research Projects in China (G19990222).     

Importance of unsteady contributions to force and heating for particles in compressible flows. Part 2: Application to particle dispersal by blast waves

Particle dispersal by blast waves is an interesting phenomenon. A model problem, i.e., a sudden release of a compressed gas–particle mixture contained in a spherical container, is employed to investigate the fundamental physics of particle dispersal. The problem is simulated by the multiphase flow models proposed in Part 1 of this article that include unsteady contributions in momentum and energy exchange between gas and particles. At early times, when particles are accelerated in the expansion fan, unsteady force and heating contributions are much larger than the corresponding quasi-steady contributions. Consequently, neglecting unsteady contributions leads to significant errors in particle front location (the boundary of the particle cloud). The complex wave interactions in the flow have strong influence on the particle motion. As a result, the particle motion is a non-monotonic function of particle density or diameter and the evolution of particle concentration is non-uniform and unsteady.     

Discrete element method study of shear-driven granular segregation in a slowly rotating horizontal drum

Segregation and mixing of granular materials are complex processes and are not fully understood. Motivated by industrial need, we performed a simulation using the discrete element method to study size segregation of a binary mixture of granular particles in a horizontal rotating drum. Particles of two different sizes were poured into the drum until it was 50% full. Shear-driven segregation was induced by rotating the side-plates of the drum in the opposite direction to that of the cylindrical wall. We found that radial segregation diminished in these systems but did not completely vanish. In an ordinary rotating drum, a radial core of smaller particles is formed in the center of the drum, surrounded by larger revolving particles. In our system, however, the smaller particles were found to migrate toward the side-plates. The shear from anti-spinning side-plates reduces the voidage and increases the bulk density. As such, smaller particles in the mixer tend to move to denser regions. We varied the shear by changing the coefficient of friction on the side-plates to study the influence of shear rate on this migration. We also compared the extent of radial segregation with stationary side-plates and with side-plates moving in different angular directions.     

Three-dimensional simulation of a dust lifting process with varying parameters

This article presents modelling considerations and simulation results for a dust lifting process in a three-dimensional domain. The Eulerian–Lagrangian modelling technique is used. Multiple simulations with different values for the number of particles were performed. The results of the simulations are shown as snapshots of particle position at certain points in time after the passage of a shock wave. Statistical data for the particle positions and collisions are presented. These are: the average height of the particles, the mean square displacement of the particles and the cumulative number of recorded collisions plotted as functions of time. The particle averaged kinetic energy and the mechanical energy lost by particles during collisions are recorded as functions of time in order to study the motion of particles. The results show that simulations of an increasing number of particles render a less intense lifting effect and, more importantly, that the inter-particle and particle–wall collisions represent essential phenomena and need to be included in this type of model. Also, a comparison between two-dimensional and three-dimensional simulations was performed. It was found that, although 2D simulations are still useful, they overestimate the lifting process and therefore a 3D model is preferable. The influence of the magnitude of the restitution and friction coefficients on the process was also studied.