直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

颗粒沉降运动的虚拟区域法直接数值模拟

文采用基于四边形网格的分布式拉格朗日乘子／虚拟区域方法(DLM／FD method)对二维方槽内775个圆形颗粒在流体中的沉降过程进行了直接数值模拟。得到了颗粒流沉降过程中流体和颗粒速度和涡量分布、流场压力分布等流动细节,展示了颗粒在沉降过程中由于相间的相互作用以及颗粒间的作用,使得颗粒流在流场内形成大小不一的旋流区,颗粒回旋着沉降,同时颗粒的尾涡影响附近颗粒的运动．本文的结果说明分布式拉格朗日乘子／虚拟区域方法对模拟存在很多颗粒的悬浮体流动是可行的。     

双颗粒在溶解条件下沉降的多相流动特性

对牛顿流体内溶解双颗粒在垂直管道中的沉降运动进行了直接数值模拟. 流体运动由守恒方程计算, 密度和黏性的变化考虑流场温度变化的影响, 通过积分粘性应力和压力获得颗粒的受力跟踪颗粒运动, 溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定. 通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合. 对双颗粒在等温流体无溶解条件和非等温流体溶解条件下的沉降过程进行了计算. 结果表明, 在一定雷诺数内, 热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒质量、颗粒表面形态的变化引起了颗粒的横向摆动, 并使颗粒沉降速度发生了变化.     

热对流条件下固液密度比对颗粒沉降运动影响的直接数值模拟

本文应用任意拉格朗日-欧拉(ALE)算法对固液两相流流场中考虑热对流的非等温颗粒在竖直通道中的沉降运动进行了数值模拟.在牛顿流体中通过积分黏性应力和压力获得颗粒的受力跟踪颗粒运动,使用有限元方法数值求解流场的N-S方程和能量方程,模型不需经验假设.通过模拟来研究颗粒沉降的运动规律和热对流下固液密度比对固液两相流的影响作用.结果表明随着固液密度比的增加,颗粒经历了稳定沉降、周期性摆动,不规则摆动等过程;热对流使颗粒的摆动幅度和沉降速度发生变化;热对流对颗粒的影响作用随着固液密度比的增加而减小.     

三维Kolmogorov流动弹性湍流直接数值模拟研究

黏弹性流体弹性湍流现象是在极低雷诺数下(惯性效应可以忽略),完全由流体弹性引起的复杂湍流现象,在此条件下牛顿流体处于层流流态.该现象不仅具有与惯性湍流类似的一些特征,还具有其独特的特性,此外,其诱发机理与惯性湍流完全不同.为了更深刻地认识弹性湍流特性及其产生机理,本文基于三维Kolmogorov流动采用直接数值模拟研究...     

无网格方法中的背景积分方案及单颗粒下降过程的数值模拟

在用无网格方法对有边界移动的气固两相流动进行数值模拟时发现,如果采用背景积分网格对求解域进行积分,将会产生压力场的数值振荡,使计算结果失真.详细分析了产生这种数值振荡的原因,提出了解决方案,并计算了一个颗粒在管道内的下降过程.计算结果表明,提出的解决方案能有效地降低压力场的数值振荡,从而使无网格方法能应用于气固两相流动的直接数值模拟.     

黏弹性流体泡状流混合层流动的PIV实验研究

本文对注入气泡的黏弹性流体气液两相泡状流混合层流动进行实验研究,含气率为0.5%。与纯水气液两相流相似,气泡的注入也会削弱黏弹性流体中的涡结构,影响大涡的卷起。雷诺应力与涡量的分布规律类似,其峰值集中在混合层中心区域。注入气泡使混合层中心区域内的雷诺应力增大,峰值区域范围减小,添加聚合物的混合层中的雷诺应力峰值增大得更多。对涡量而言,注入气泡使峰值产生非常明显的减小,黏弹性流体的涡量所受的影响比纯水的要小。     

黏弹性流体基铜纳米流体流动与传热实验研究

成功建立了流体流动阻力和换热性能测试实验台,在45℃的流体温度下,对不同铜粒子体积分数和基液浓度的纳米流体在湍流状态下的对流换热特性和流动阻力进行了实验测量。实验结果表明:黏弹性流体基液中添加纳米粒子后,在降低对应基液减阻率的同时能明显增强传热性能.例如,将1.0%体积分数的铜纳米粒子添加到质量分数为6×10~(-4)的基液中所形成的黏弹性流体基纳米流体其综合性能指数K=0.47,表现了很好的传热强化和减阻性能.     

二维单液滴内颗粒沉降特性的数值模拟研究

燃煤烟气中的含细微颗粒物浆滴经烟气夹带进入大气会对电厂周边环境造成较大的危害.高效脱除含细微颗粒物的湿烟气冷凝液滴是防控细微颗粒物污染的有效手段,而细微颗粒物在液滴内部的流动与沉降行为对液滴脱除滑动特性有重要影响.基于此,本文建立了微米尺度固体颗粒在二维单液滴内的自由沉降运动模型,数值模拟研究了颗粒粒径、液滴形貌参数、...     

圆截面直流道中微粒黏弹性聚焦机理研究

微粒黏弹性聚焦技术近年来受到了广泛的研究重视, 但影响粒子聚焦特性的关键参数调控机理仍不清楚. 基于此目的, 本文量化研究了圆截面直流道中非牛顿流体诱导微粒黏弹性聚焦的行为, 给出了流速和流道长度对粒子聚焦特性的调控机理. 具体而言: 首先, 对比分析不同黏度牛顿流体(水和22 wt%甘油水溶液)和非牛顿流体(8 wt%聚乙烯吡咯烷酮水溶液)中粒子横向迁移行为, 发现非牛顿流体中粒子将在弹性力主导下聚焦至流道中心区域, 而牛顿流体中粒子则在惯性升力主导下迁移形成Segré-Silberberg圆环. 其次, 量化分析粒子尺寸和驱动流速对黏弹性聚焦效果的影响, 发现随着流速的增加, 粒子聚焦效果逐渐变好并最终趋于稳定, 且大粒子较小粒子具有更好的聚焦效果. 最后, 研究粒子沿流道长度的动态聚焦过程, 推导并验证了粒子聚焦所需安全流道长度的数学模型, 发现大粒子聚焦所需安全流道长度显著短于小粒子. 上述研究结果对于提升粒子黏弹性聚焦机理和过程的理解, 实现微粒聚焦特性的灵活控制具有非常重要的意义.     

Direct numerical simulation of elastic turbulence and its mixing-enhancement effect in a straight channel flow

In this paper,we present a direct numerical simulation(DNS) of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number(Re = 1) in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.     

Sedimentation of Brownian Particles in a Gravitational Potential

The settling of a Brownian particle in a semi-infinite fluid bounded by a bottom plane is studied on the basis of Smoluchowski's exact solution of the equation describing diffusion in the gravitational potential. Expressions are derived for the mean height and the variance of height at some time after starting at an initial height. These quantities show interesting behavior as a function of time. It is shown that for certain initial heights the Boltzmann entropy does not increase steadily. It increases at first but then decreases to its equilibrium value.     

Characteristics and generation of elastic turbulence in a three-dimensional parallel plate channel using direct numerical simulation

Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.     

矩形网格上VOF运动界面重构的流体体积分数保持法

针对处理运动界面问题的流体体积函数(VOF)法,提出一种新的流体体积分数保持界面重构算法.该方法在单个网格内用斜线段近似运动界面,并要求相邻网格内的斜线段在公共边上的交点重合,通过保持网格内的流体体积分数不变,在边界网格上建立非线性方程组,通过求解非线性方程组,确定运动界面边界线与网格线的交点位置,最终重构出运动界面的形状.应用算例模拟结果表明,该方法能够高精度,较少尖点地重构运动界面,证实了该方法的有效性和可行性.     

Correlation of Two Coupled Particles in Viscoelastic Medium

Considering the viscoelastic memory effect, we study the correlated motion of two hydrodynamically coupled colloidal particles, each of which confined in a harmonic potential well, in a Kelvin-type and Maxwell-type viscoelastic medium. We find that viscoelastic relaxation plays a significant role in modifying the correlation, particularly the cross correlation. We also find that both the real and imaginary parts of the response function are significantly different from the viscous medium case. In particular there is a phase shift between the vanishing imaginary part and the maximal real part of the response function in a viscoelastic medium. In addition imaginary part of the cross correlation response function exhibits a net energy loss （gain） behavior when the elasticity parameter of the medium is larger （smaller） than the critical value for Kelvin （Maxwell） viscoelastic fluid. Some implication of our results and their connection with previous works are discussed.     

Correlation of Two Coupled Particles in Viscoelastic Medium

Considering the viscoelastic memory effect, we study the correlated motion of two hydrodynamically coupled colloidal particles, each of which confined in a harmonic potential well, in a Kelvin-type and Maxwell-type viscoelastic medium. We find that viscoelastic relaxation plays a significant role in modifying the correlation, particularly the cross correlation. We also find that both the real and imaginary parts of the response function are significantly different from the viscous medium case. In particular there is a phase shift between the vanishing imaginary part and the maximal real part of the response function in a viscoelastic medium. In addition imaginary part of the cross correlation response function exhibits a net energy loss (gain) behavior when the elasticity parameter of the medium is larger (smaller) than the critical value for Kelvin (Maxwell) viscoelastic fluid. Some implication of our results and their connection with previous works are discussed.