不同雷诺数下二维椭圆柱绕流的数值模拟研究

利用计算流体力学软件Fluent对不同雷诺数(Re=100,3900,3.5×10~6)下二维椭圆柱绕流进行了数值模拟研究,分析了不同轴长比(2b/2a=cosθ,θ=0°,15°,30°,45°,60°)下椭圆柱绕流的特性。通过对比尾部涡流情况、升力系数C_L、阻力系数C_D以及斯特劳哈尔数St初步发现:由于椭圆形截面偏流线型的特点,在三种雷诺数下随着θ的增大椭圆柱绕流尾涡强度减小,流场的变化使圆柱表面的压力系数减小,最终导致圆柱的升力系数幅值与阻力系数均值减小。而斯特劳哈尔数St在三种雷诺数下的变化不同,随着θ的增大,层流雷诺数(Re=100)下St值减小;亚临界雷诺数(Re=3900)下St值在45°处轻微上扬,在60°处明显减小;超临界雷诺数(Re=3.5×10~6)下St值增大。     

稳态条件下气溶胶粒子在口腔模型中的沉淀率

以实际人体解剖学数据为基础,利用流体力学软件构建了口腔模型。由于口腔形状不规则,流场中会有湍流产生,所以本文选用的是k-ω方程。在此基础上,假设吸入空气流量为一个常量,改变吸入气体、气溶胶粒子的物理参数,以调整St数的变化,模拟不同情况下气溶胶粒子运动情况。由模拟结果可知:气溶胶粒子在口腔中的沉淀率受粒子半径、密度、及流速的影响,并随这些量的增加而增加;沉淀率可以视作St数的函数,并随St数的增加而增加;气溶胶粒子在口腔中的沉淀分布受模型几何形状和湍流的影响,在喉部附近沉淀较为明显;粒子在口腔中将沿怎样的轨迹运动则取决于粒子在入口处的初始位置。     

壁面约束对柱状粒子在牛顿流体中沉降影响的研究

格子Boltzmann方法(LBM)直接数值模拟了在有壁面约束的流场中柱状粒子的沉降。结果说明粒子在壁面的影响下，会离开管壁向管中心移动，移动速度先是增大，然后随着与管壁距离的增加，逐渐减小，粒子中心最后并不一定停留在管中心。随着粒子与壁面距离的增加，粒子的沉降速度将增大，可见壁面对粒子沉降有阻碍作用，平行于粒子轴线的壁面对粒子沉降的影响比垂直于粒子轴线壁面的影响大。     

涡旋波流场的涡量测试与计算及特征参数的影响

从实验和模拟计算两个方面,对涡旋波流场的速度矢量图和涡量进行了定量测试和综合分析,并对涡旋波流动与其特征影响参数之间的关系进行了研究。结果表明:涡旋波流动只在二维层流状态下产生,涡旋波形成后,Re数对旋涡尺寸的影响较小,涡量随Re数的增大而增大,但涡量不随时间的增加而单调增加。随着St数的减小,旋涡尺寸明显增大,涡量却随之减小。本文也对不同槽道内的涡旋波流场进行了数值模拟计算,进而确定了涡旋波流场的形成条件。     

扩散沉积导致的管内气溶胶穿透率简明公式

针对气溶胶粒子在长直扁平方管(通道)及圆管内受布朗扩散导致的沉积损失问题,从对流扩散方程出发,采用分离变量法,通过严格的理论求解得到了这两类问题中粒子穿透率的无穷级数解。利用Python语言中高精度计算功能确定了相应本征问题的本征值、本征函数及无穷级数展开系数。根据所得的相对严格意义上的穿透率数据,回归拟合得到了扩散沉积条件下粒子通过方管及圆管时的穿透率简明公式,并与已有研究进行了详细的比较。结果表明:本文得到的两个简明公式更加简单和准确,对方管与圆管内粒子穿透率估计的误差最大值分别仅为3.7×10~(-4)和8.3×10~(-4),可以快速有效地估计层流条件下和受扩散沉积作用时,气溶胶粒子通过两种管(通)道的穿透率。     

纳米AlOOH及纳米Fe_3O_4粒子在液体石蜡中的摩擦学性能研究

采用四球摩擦磨损试验机对表面修饰的纳米AlOOH粒子及纳米Fe3O4粒子在液体石蜡中的摩擦学性能进行了对比研究.结果表明:这2种纳米粒子均能提高液体石蜡的减摩耐磨性能和PB值,纳米AlOOH粒子因具有层状结构,表现出更好的减摩耐磨性能,随着纳米粒子粒径的增大,其减摩耐磨的最佳浓度出现升高的趋势.对磨斑表面的SEM、AFM和XPS分析结果表明,纳米粒子能沉积在摩擦副表面,减少摩擦副表面微凸体的直接接触,从而减少微凸体之间的犁削和黏着.     

大气边界层内羽流扩散研究

主要研究了高架点源的污染,即``羽流扩散'. 由于大气边界层内的湍流运动,引 起污染物扩散的主要因素是湍流扩散, 研究湍流扩散有两种基本方法：统计理论和梯度理论. 采用拉格朗日粒子扩散模型 (Lagrangian particle dispersion model, LPDM)来评价羽流扩散的平均浓度,分别计算了表征扩散程度的3个参数：沿烟轴方向的 着地浓度、水平方向的羽流扩散宽度和垂直方向的羽流扩散宽度. 从计算结果和风洞实验的 数据对比来看,用该模式可以较好地模拟大气边界层内复杂气候条件下的羽流扩散.     

分形介质中的反常扩散

本文利用适当的标度变换将欧氏空间中的经典扩散方程抗议为分形点阵上的标准扩散方程，并应用该标准扩散方程证明了分形布朗粒子的运动服从反常扩散，同时还讨论了一般分形扩散方程及其渐进解。     

冷喷涂中粒子与基体的高速冲击过程

利用非线性动力分析软件LS-DYNA,对冷喷涂中铜粒子与铜基体的碰撞过程进行了3维数值分 析。研究沉积粒子与被喷涂基体的变形行为以及影响数值模拟结果的关键因素。数值模拟结果表明:通过3 维模型获得的铜粒子与铜基体在碰撞结束后的沉积形貌与实验观察吻合很好,比采用2维模型获得的模拟结 果更精确。另外,3维模型的网格尺寸对喷涂粒子的临界速度和粒子的变形行为也有很大影响,随着网格尺 寸的减小,粒子的临界速度及扁平率逐渐减小,压缩率逐渐增大,利用插补法获得的网格尺寸为0时的铜粒子 临界速度、扁平率及压缩率与实验结果吻合很好。     

两种流型下雾滴表面气溶胶粒子的扩散沉积

用浓度边界层理论分析了势流和Hadamard流情形下雾滴或云滴前向驻点处小粒子的扩散效应。结果表明：存在静电效应时，有限直径小粒子在雾滴表面的扩散通量与两个无量纲数参数有效；对于无外场作用时的点粒子情形，上述两种流型下的扩散沉积特征与粘性流结果亦不相同。     

Effects of vortex pairing on particle dispersion in turbulent shear flows

Particle dispersion in large-scale dominated turbulent shear flow is investigated numerically with special emphasis on the effects of the vortex-pairing phenomenon. The particle dispersion is visualized numerically by following the particle trajectories in a flow consisting of large vortices which are undergoing pairing interaction. The flow field is generated by a discrete vortex method. Important global and local fiow quantities from the numerical simulation compare reasonably well with experimental measurements.

For both cases of point sources with continuous particle release and an initially distributed line source, the particle dispersion results demonstrate that the extent of particle dispersion depends strongly on the Stokes number, the ratio of the particle aerodynamic response time to the characteristic time of the vortex-pairing flow field. Particles with relatively small Stokes numbers disperse laterally at approximately the saine rate as that of the fluid particles and particles with large Stokes numbers disperse much less than the fluid particles. Particles with intermediate Stokes numbers (0.5-5) may be dispersed laterally farther than the fiuid particles and may actually be flung out of the vortex structures. Due to the strong particie entrainment power, the flow during the vortex-pairing process seems to produce higher particle lateral dispersion than the pre-pairing and post-pairing flows.     

NUMERICAL STUDY OF PARTICLE DISTRIBUTION IN WAKE OF LIQUID-PARTICLE FLOWS PAST A CIRCULAR CYLINDER USING DISCRETE VORTEX METHOD

Particle-laden water flows past a circular cylinder were numerically investigated. The discrete vortex method (DVM) was employed to evaluate the unsteady water flow fields and a Lagrangian approach was applied for tracking individual solid particles. A dispersion function was defined to represent the dispersion scale of the particle. The wake vortex patterns, the distributions and the time series of dispersion functions of particles with different Stokes numbers were obtained. Numerical results show that the particle distribution in the wake of the circular cylinder is closely related to the particle's Stokes number and the structure of wake vortices: (1) the intermediate sized particles with Stokes numbers, St, of 0.25, 1.0 and 4.0 can not enter the vortex cores and concentrate near the peripheries of the vortex structures, (2) in the circular cylinder wake, the dispersion intensity of particles decreases as St is increased from 0.25 to 4.0.     

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

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

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

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

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

COMPUTATION OF DILUTE PARTICULATE LAMINAR FLOW OVER A BACKWARD-FACING STEP

Particle-laden flows are calculated for a classical laminar backward- facing step problem. The particle tracks are calculated using a recently developed exponential Lagrangian tracking scheme. The behaviour of the particle-laden flow is considered for various inlet for Reynolds number, Stokes numbers and void fractions. Doping the flow with low-Stokes-number particles has the effect of increasing the inlet inertia of the flow and this increases the strength of the recirculation behind the step. High-Stokes- number particles are dominated by gravitational effects which affect the flow accordingly. Differences between the single-phase flow and the particle-laden flows are therefore dependent on the Stokes number and increase linearly with void fraction.     

稀疏两相湍流的惯性颗粒倾向性弥散

气相采用大涡模拟、颗粒相采用拉格朗日轨道模型的方法对后台阶突扩流、充分发展槽道流和圆湍射流3种典型的稀疏气固两相流动进行了数值模拟,研究了颗粒倾向性弥散的特征和规律. 研究表明颗粒的跟随性和倾向性相联系,颗粒惯性和大涡结构同时决定颗粒的倾向性分布特征. Stokes数量级为1(气相时间参考尺度取为宏观特征时间尺度)左右的颗粒,倾向性分布特征最强烈. 颗粒倾向分布于低涡量(或是低脉动速度)的湍流区域.      

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

Large eddy simulation of polydispersed inertial particles using two-way coupled PDF-PBE

This paper presents large eddy simulation (LES) of polydispersed turbulent recirculating flows using a two-way coupled probability density function of the population balance equation (PDF-PBE). A stochastic Monte Carlo method is adopted to solve the PDF-PBE on an ensemble of notional Lagrangian particles and the method of Stokes binning, that was recently developed by Salehi et al. (2017) is employed to explicitly treat effects of inertia. The PDF-PBE is applied to the experiment of Boŕee et al. (2001) which studied dispersion of polysized inertial particles in a bluff body configuration. The particle mass loading is 22% where the dispersed elements affect the carrier phase velocity. The simulations are performed using both the standard Smagorinsky and the Wall-Adapting Local Eddy-viscosity (WALE) subgird turbulent models. It is found that the subgrid model has a significant impact on the results, particularly on the carrier velocity. The WALE model shows a better agreement with the measurements. Different boundary conditions are tested for injection of notional particles. It is demonstrated both dispersed and carrier phase velocities are initially sensitive to the particle boundary condition whereas the difference between tested conditions becomes marginal further downstream. The best results are obtained for the particle boundary condition that accounts for effects of inertia at the inlet. Finally, the sensitivity of the PDF-PBE simulations to the number of Stokes bins is studied. It is found that eight Stokes bins are enough to accurately model the polysized particles dispersion in a bluff body configuration. The results indicate that particle dispersion is notably sensitive to the number of Stokes bins but the particle velocity predictions are much less sensitive with small variation in velocity results over the range of Stokes bins that were tested.     

Numerical investigation of gas-particle flow in the primary air pipe of a low NOx swirl burner – The DEM-CFD method

The gas-particle flow in the primary air pipe (PAP) of a low NO_x swirl burner was investigated using the computational fluid dynamics (CFD) coupled with the discrete element method (DEM). The mathematical models were validated using the measured values obtained at the outlet of the primary pipe through a phase Doppler anemometer (PDA) system. Particles of different Stokes numbers in the primary air pipe (PAP) were investigated, and the effects of the structure of the primary air pipe and the particle–particle interaction on particle dispersion were analyzed. The results indicate that particles under the combined effects of the Venturi pipe and the spindle body are concentrated into a narrow band area and that the PAP structure can more efficiently concentrate particles with large Stokes numbers. The formed fuel rich/lean jet persists for a long distance out of the burner, thereby favoring of air-staged combustion and NO_x reduction. The particle collision frequency and its fluctuation range increase as the particle Stokes number increases. The collisions among particles result in an increase of the spanwise dispersion of particles. Experimental results indicate that the models that take particle–particle collision into consideration are more able to predict particle concentration.