基于变截面微通道的黏弹性流体内颗粒高效分离研究

黏弹性流体广泛存在于自然界中,如人体中的血液等。实现黏弹性流体中不同尺寸微颗粒的高效分离对于生命科学和临床医学等领域有着重要的意义。本文基于对黏弹性流体中的微颗粒先富集再分离的思想,首先通过渐缩截面微通道,改变弹性升力的方向,并增强微颗粒在微通道中的受力,实现不同微颗粒的高效富集。而后,利用不同粒径微颗粒在层流状态下的运动特性差异,进一步实现对不同尺寸微颗粒的高效分离。实验结果表明,在维森伯格数Wi为17.5至34.9的范围内,聚乙烯吡咯烷酮(PVP)黏弹性流体中10 μm与4 μm两种微颗粒可实现完全分离。此外,本文还研究了不同流动参数和通道几何结构对黏弹性流体中颗粒分离的影响。与其他用于黏弹性流体中微颗粒分离的微流控技术相比,本文提出的颗粒分离方法具有分离精度和效率高、通道长度短等优点。基于此方法的微流控芯片技术在生物医学等领域有着巨大的应用潜力。     

微尺度稀薄效应下颗粒沉降的格子Boltzmann方法模拟

基于多松弛格子Boltzmann模型,对竖直细长微通道内颗粒自由沉降过程进行模拟,分析气体稀薄效应、初始位置以及颗粒间相互作用对微颗粒沉降特性的影响.研究表明：随Knudsen数增大,微通道内气体稀薄效应增强,颗粒表面气体滑移速度增大,气相流体有效粘度减小,颗粒相同运动状态下受到气体阻力相应减小,颗粒沉降平衡速度明显增大;不同初始位置颗粒沉降过程存在明显差异,初始位置偏离中心线颗粒将发生水平方向位移且呈振荡趋势,最终稳定于中心线平衡位置;在微尺度双颗粒沉降DKT现象过程中,气体稀薄效应影响颗粒运动特性,后颗粒跟随过程明显增长.     

三层材料纳米颗粒形状对消光特性的影响

通过高分辨率电子束光刻方法制备了不同形状的三层复合材料纳米颗粒,研究了这种纳米颗粒的形状变化对消光特性的影响。测试结果表明,当入射波偏振方向平行于短轴时,随着长宽比的增大,共振峰位置发生蓝移;当光源偏振方向平行于长轴时,随着长宽比的增大,共振峰位置发生红移。还用时域有限差分算法以及表面等离波子的Lorentz模型对纳米颗粒的消光特性进行数值计算,所得的消光频谱曲线、共振峰位置变化趋势与实验基本一致。此外,还研究了主体材料层厚度对消光特性的影响,发现其厚度在20～90nm变化时,共振峰发生3～115nm的蓝移。     

通道宽度和初始流量对颗粒稀疏流-密集流转变临界开口的影响

通过实验研究斜面上二维颗粒流，当出口尺寸减小到临界值Dc时，发生稀疏流到密集流的突变.发现临界尺寸Dc与初始流量和通道宽度有关，通道宽度一定的情况下，临界开口尺寸Dc近似随初始流量Q0的平方根增大.在初始流量Q0一定时，临界开口尺寸Dc随通道宽度W近似线性增大.给出了这些关系的表达式，理论计算与实验观测结果一致.同时，也讨论了通道宽度影响临界开口尺寸的原因. 关键词： 颗粒物质 颗粒流     

非牛顿流体在介质中的自然对流与传热

基于修正的Darcy定律,数值模拟了上端开口,从底部等温加热的多孔介质中黏弹性流体的自然对流与传热,其中使用了特征松弛时间λ和迟滞时间ε来表征黏弹性流体的性质。我们选择了三组不同的λ和ε.发现与牛顿流体中对应的情况相比,黏弹性流体的传热特性有很大不同。另外,上端开口多孔介质中的黏弹性流体比上端闭口多孔介质中黏弹性流体的自然对流与传热也有很大不同.在不同的Ra数范围,传热规律服从不同的标度律。     

湿颗粒堆力学特性的离散元法模拟研究

利用基于线性黏聚接触模型的离散元法对不同颗粒系统的堆积过程进行了数值模拟研究,分析了颗粒形状和湿颗粒间液桥力对颗粒堆积形态的影响机理,获得了球形和块状湿颗粒堆基底表面所受的法向力以及堆中颗粒间的法向力和切向力"中心凹陷"式的分布规律,讨论了颗粒形状和黏聚能量密度对基底表面作用力和颗粒间作用力的影响.研究结果表明,颗粒形状和液桥力对颗粒堆的堆积形态具有显著的影响.堆积角随着黏聚能量密度的增加而增大,并且相同条件下的块状颗粒堆积角大于球形颗粒.颗粒形状和黏聚能量密度对基底表面所受作用力和堆中颗粒间的作用力变化及最大幅值均有影响作用.当黏聚能量密度值逐渐增大时,颗粒堆的作用力最大幅值均逐渐增大,并且块状颗粒堆的作用力最大幅值大于球形颗粒堆.当黏聚能量密度值过大时,颗粒堆力学特性更加复杂,液桥力对颗粒堆积特性的影响作用大于颗粒形状的影响.     

滴形颗粒在通道中沉降的格子Boltzmann-虚拟区域方法模拟

采用格子Boltzmann-虚拟区域方法对滴形颗粒在垂直通道中的沉降过程进行直接数值模拟.通过格子Boltzmann方法求解N-S方程,流体与固体之间的相互作用通过虚拟区域方法描述.研究雷诺数在10^-2到100范围内颗粒形状因子对其摩擦系数和阻力系数的影响.为便于比较,给出了圆形颗粒沉降的结果.结果发现,当雷诺数小于1的时候,颗粒的摩擦系数始终保持常数,而当雷诺数大于1时摩擦系数随雷诺数的增大而增大.此外,当雷诺数小于约30时圆形颗粒的摩擦系数和阻力系数均小于滴形颗粒,而当雷诺数大于30时情况正好相反.颗粒周围的压力分布证明了这一结论.     

颗粒物质从稀疏流到密集流转变的普适规律

二维颗粒流实验研究发现，当初始稀疏颗粒流量固定，出口的尺寸减小到一临界值，或固定出口的尺寸，颗粒流量增大到一临界值时，都会发生流量的突然减小，从稀疏流转变为密集流．得到了临界流量与出口尺寸、颗粒尺寸及通道宽度之间的普适标度关系，揭示了“瓶颈效应”的物理本质．     

二维颗粒流通道宽度效应的分子动力学模拟

通过用分子动力学方法对颗粒物质流的计算机模拟，研究发现增大通道宽度可以使二维颗粒流从稀疏流转变为密集流状态.通过对不同通道宽度下，固定开口为9.5d的颗粒流和 漏斗口以上9.5d×8d区域记录的模拟结果分析，发现随通道宽度增大，密度变大、温 度降低.当“颗粒温度”T较低时(T/m<0.05 J/kg),颗粒流内部接触数开始超过1.2 ，同时出现较为牢固的横向链状颗粒团簇，是造成流量突变以及密集流的原因. 关键词： 颗粒物质 颗粒流 计算机模拟     

不同纳米流体在微通道内的流动传热特性研究

采用数值模拟的方法研究了不同工质在微通道内流动传热特性的差异。对比了去离子水、纳米流体Al₂O₃/Water、CuO/Water、TiO₂/Water、Cu/Water等工质在微通道内的流动传热特性,并研究了纳米颗粒的浓度对流动换热特性的影响。结果表明：CuO/Water作为冷却工质时的对流换热系数比水增加了9.6%,微通道底面平均温度降低了2.6 K,换热性能明显优于其他几种纳米流体。由于纳米颗粒的加入,纳米流体的粘度比水大,进出口的压降比水大。纳米颗粒的体积分数越大,对流换热系数越大,纳米流体在微通道内的换热性能越好。     

Inertial focusing and rotating characteristics of elliptical and rectangular particle pairs in channel flow

The lattice Boltzmann method is used to study the inertial focusing and rotating characteristics of two-dimensional elliptical particles and rectangular particles in channel flow. The results show that both elliptical particles and rectangular particles initially located on one side and two sides of channel centerline migrate first towards the equilibrium position. Then, the single-line particle train with an increasing spacing and the staggered particle train with stable spacing are formed. The axial spacing of the staggered particle pair increases with aspect ratio and Reynolds number increasing. The staggered elliptical or rectangular particle pairs form perpendicular orientation angles, which will be more obvious at larger aspect ratio and lower Reynolds number. The single-line particle trains with different shapes seldom form the perpendicular orientation angle.     

Analytical solutions for particle transport through an inclined channel with gravitational effect

The combined mechanisms of Brownian diffusion and gravity settling are considered to investigate the transport and deposition of particles in an inclined rectangular channel in the laminar flow regime. The exact analytical solution for particle concentration is obtained with some reasonable transformations and the conventional method of separation of variables. The effects of Peclet number, depositional parameter, incline angle, and uphill and downhill airflows on the mean concentration of particles are discussed in detail. The results indicate that the exact solution obtained in the present study can describe the transport and deposition behavior of particles due to the combined mechanisms throughout a channel at any inclination angle.     

Thermosolutal convection in a Maxwellian viscoelastic fluid in porous medium

The thermosolutal convection in a layer of Maxwellian viscoelastic fluid heated and soluted from below in porous medium is considered. The effects of uniform magnetic field and uniform rotation on the thermosolutal convection are also considered. For stationary convection, the Maxwellian viscoelastic fluid behaves like a Newtonian fluid. The sufficient conditions for the nonexistence of overstability are obtained. The critical Rayleigh number is found to increase with the increase in magnetic field, rotation and stable solute gradient.     

Stability analysis of a Maxwell fluid in a porous medium heated from below

Based on a modified-Darcy–Brinkman–Maxwell model, stability analysis of a horizontal layer of Maxwell fluid in a porous medium heated from below is performed. By solving the eigenvalue problems, the critical Rayleigh number, wave number and frequency for overstability are determined. It is found that the critical Rayleigh number for overstability decreases as the relaxation time increases and the elasticity of a Maxwell fluid has a destabilizing effect on the fluid layer in porous media. On the other hand, the critical Rayleigh number for overstability increases by increasing the porous parameter which acts to stabilize the system. In limiting cases, some previous results for viscoelastic fluids in nonporous media are recovered from our results.     

Wave propagation in fluid-conveying viscoelastic single-walled carbon nanotubes with surface and nonlocal effects

In this paper, the transverse wave propagation in fluid-conveying viscoelastic single-walled carbon nanotubes is investigated based on nonlocal elasticity theory with consideration of surface effect. The governing equation is formulated utilizing nonlocal Euler-Bernoulli beam theory and Kelvin-Voigt model. Explicit wave dispersion relation is developed and wave phase velocities and frequencies are obtained. The effect of the fluid flow velocity, structural damping, surface effect, small scale effects and tube diameter on the wave propagation properties are discussed with different wave numbers. The wave frequency increases with the increase of fluid flow velocity, but decreases with the increases of tube diameter and wave number. The effect of surface elasticity and residual surface tension is more significant for small wave number and tube diameter. For larger values of wave number and nonlocal parameters, the real part of frequency ratio raises.     

Linear stability analysis of the natural convection in inclined rotating parallel plates

The linear stability analysis of the natural convection in a rectangular tilted infinite cavity filled with a Boussinesq fluid subject to Coriolis force is presented. The bottom and top surfaces have fixed temperatures. Both unstable and stable thermal conditions are studied (heated from below and heated from above respectively). The rotation axis passes through the center and it is orthogonal to the hot and cold surfaces. The stability equations were solved using the Tau–Chebyshev spectral method. The critical Rayleigh number and critical wave number were obtained for several rotation rates and different orientation of convective oblique rolls in a range of inclination of the cavity from 0 to 120 degrees. The stability analysis show that rotation rate affects the basic velocity profile, onset of convection, wave number and critical orientation of convective rolls.     

Subcritical finite-amplitude solutions for plane Couette flow of viscoelastic fluids

Plane Couette flow of viscoelastic fluids is shown to exhibit a purely elastic subcritical instability at a very small-Reynolds number in spite of being linearly stable. The mechanism of this instability is proposed and the nonlinear stability analysis of plane Couette flow of the Upper-Convected Maxwell fluid is presented. Above a critical Weissenberg number, a small finite-size perturbation is sufficient to create a secondary flow, and the threshold value for the amplitude of the perturbation decreases as the Weissenberg number increases. The results suggest a scenario for weakly turbulent viscoelastic flow which is similar to the one for Newtonian fluids as a function of Reynolds number.     

Choice of thickness ratio of a coated beam used for investigating the complex modulus of viscoelastic materials

The coated beam method has been widely used for investigating the complex modulus of elasticity of relatively soft viscoelastic materials; e.g., materials for vibration damping. With this method a rectangular section metal beam is coated with the viscoelastic material on one side, or symmetrically on both sides, and the resonances of the bending vibration of the beam are investigated. In this paper a procedure is described for finding the ratio of coating thickness to metal thickness required to obtain accurate results for the complex modulus in such investigations. The contradiction experienced earlier—that a relatively large thickness ratio is required for precision while the nature of the resonance method is such as to demand a small one—is analyzed mathematically. The relationships derived serve as a basis for optimizing the choice of the ratio. Diagrams and a procedure for this are worked out and presented. Finally, a few experimental results obtained with suitable and unsuitable thickness ratios 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.     

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.