REVIEW OF SOME RESEARCHES ON SUSPENSION OF SOLID PARTICLE IN NON-NEWTONIAN FLUID
-
摘要: 非牛顿流体固粒悬浮流具有广泛的应用背景，其特殊的流动属性使其成为一些新兴技术领域的核心突破点。同时，该流动又比较复杂，即便是在低固粒浓度的情况下，非牛顿流体特性也会对整个系统的微结构产生重要的影响，从而进一步影响固粒的运动。本文给出了非牛顿流体方程、固粒运动方程和非牛顿流体固粒悬浮流的特征参数，分析了这些参数的作用；阐述了单个固粒在管道中的径向移动、多固粒的相互作用和聚集、多固粒形成的链状结构以及非圆球固粒运动等方面的研究成果、结果分析以及尚未解决的问题，并对以上问题进行了总结和展望，给出了需要深入研究的具体问题和内容，旨在为进一步的研究提供参考和依据。Abstract: Suspension of solid particle in non-Newtonian fluid has a wide range of applications, and its special flow properties make it the core breakthrough point in some new technology fields. Meanwhile, the flow is more complex. Even in the case of low particle concentration, the characteristics of non-Newtonian fluid have an important influence on the microstructure of the whole system, which further affects the movement of solid particles. In this paper, the non-Newtonian fluid equation, particle motion equation and characteristic parameter of suspension of solid particle in non-Newtonian fluid are given, and the effect of these parameters is analyzed. The research findings, result analysis and open questions of some topics including the radial motion of single solid particle in a pipe, the interaction and aggregation of multi-particles, the chain structure formed by multi-particles, and motion of non-spherical particles are related. Finally, the topics mentioned above are summarized and prospected, and the concrete problems and contents that need to be studied deeply are given, which is aimed at providing references and basis for further research.
-
Keywords:
- non-Newtonian fluid /
- solid particle /
- suspension /
- review
-
引言
牛顿流体固粒悬浮流已得到了广泛的研究[1-3]，自然界和实际应用中，除了水、空气、甘油等少数牛顿流体外，绝大部分是非牛顿流体.非牛顿流体是一种具有内部结构的复杂流体，其基本特征是应力和应变率之间不满足牛顿内摩擦定律，流体黏度依赖于剪切率，有些具有弹性或塑性，流体介质具有记忆效应即流体介质的力学属性与介质的历史状态有关.非牛顿流体固粒悬浮流具有广泛的应用背景，且这种流动因其特殊的属性而成为一些新兴技术领域的核心突破点，例如石油工业中钻井时常用泥浆和聚合物构成的非牛顿流体来输运岩石碎末；由清水、砂和添加剂构成的压裂液通过液压使岩石产生裂纹而产出石油或天然气；又如火箭推进系统中的凝胶燃料，该燃料具有高燃烧能量、稳定耐贮存的显著优势.
非牛顿流体固粒悬浮流的动力学特性非常复杂，即便是在低固粒浓度的情况下，非牛顿流体特性诸如剪切变稀、拉伸变稠、第一和第二法向应力的存在都会对整个系统的微结构产生重要影响.固粒在流体中的分布不均匀会使悬浮流的流变特性发生剧烈的变化.然而，非牛顿流体内部结构的复杂性，也可被用来使固粒按照预期的且在牛顿流体中无法奏效的方式运动，例如在微流控技术中，让固粒在装置的某些指定区域内运动，从而达到计数、诊断和分离等目的；又如通过对固粒运动和分布的控制，可以优化产品的生产过程.因此，了解和掌握非牛顿流体固粒悬浮流的动力学特性具有重要的学术价值和实际意义.
对非牛顿流体固粒悬浮流的系统研究起源于20世纪70年代，初始阶段主要是对简单流场进行实验研究，旨在揭示在相似的条件下与牛顿流体相比而呈现出的复杂特性.随着实验仪器与技术的发展，人们开始关注高浓度固粒情形和复杂流场.此外，随着计算机技术以及数值计算方法的发展，数值模拟逐渐成为对非牛顿流体固粒悬浮流研究的重要手段.由于对非牛顿流体固粒悬浮流研究的内容涉及流体性质、固粒性质和流场特性，覆盖的范围较宽，因此以下只对单个固粒在管道中的径向移动、多固粒的相互作用和集聚、多固粒形成链状结构以及非圆球固粒的运动进行阐述.
1. 基本方程和特征参数
1.1 非牛顿流体方程
假设流体不可压缩，非牛顿流体方程的连续性方程和动量方程为
$$\nabla \cdot \mathit{\boldsymbol{u}} = 0$$ (1) $${\rho _{\rm{f}}}(\frac{{\partial \mathit{\boldsymbol{u}}}}{{{\rm{d}}t}} + \mathit{\boldsymbol{u}} \cdot \nabla \mathit{\boldsymbol{u}}) = \nabla \cdot \sigma {\rm{ }}$$ (2) 式中， ${\pmb u}$ 是速度， $\rho_{\rm f}$ 是流体密度， ${\pmb \sigma }$ 是总应力张量.不同于牛顿流体，非牛顿流体的复杂性体现在 ${\pmb \sigma }$ 的多样性和复杂性，对于黏弹性流体， ${\pmb \sigma }$ 可表示为
$$\sigma = - p\mathit{\boldsymbol{I}} + 2{\eta _S}\mathit{\boldsymbol{D}} + \mathit{\boldsymbol{\tau }}$$ (3) 式中， $p$ 是压力， ${\pmb I}$ 是单位张量， $\eta_{S}$ 是黏性系数， ${\pmb\tau }$ 是黏弹性应力， ${\pmb D}$ 是应变率张量
$$\mathit{\boldsymbol{D}} = [\nabla \mathit{\boldsymbol{u}} + {(\nabla \mathit{\boldsymbol{u}})^{\rm{T}}}]/2$$ (4) 对不同的流体要采用不同的本构关系来表示黏弹性应力 ${\pmb\tau}$ ，例如
$$\lambda \mathit{\boldsymbol{\widehat \tau }} + \mathit{\boldsymbol{\tau }} = 2{\eta _p}\mathit{\boldsymbol{D}}$$ (5) 表示Upper-Convected Maxwell流体(方程(3) 中的 $\eta_{S}$ 为零)和Oldroyd-B流体，其中, λ是流体的松弛时间， $\eta_{p}$ 是非牛顿流体对黏度的贡献， $\hat{\pmb\tau}$ 为
$$\mathit{\boldsymbol{\widehat \tau }} \equiv \frac{{\partial \mathit{\boldsymbol{\tau }}}}{{\partial t}} + \mathit{\boldsymbol{u}} \cdot \nabla \mathit{\boldsymbol{\tau }} - {(\nabla \mathit{\boldsymbol{u}})^{\rm{T}}} \cdot \mathit{\boldsymbol{\tau }} - \mathit{\boldsymbol{\tau }} \cdot \nabla \mathit{\boldsymbol{u}}$$ (6) 1.2 固粒运动方程
在非牛顿流体中运动的固粒其平动满足牛顿第二定律，转动满足Euler方程[4]
$${m_i}\frac{{{\rm{d}}{\mathit{\boldsymbol{U}}_i}}}{{{\rm{d}}t}} = {\mathit{\boldsymbol{F}}_{{\rm{ext}},i}} + {\rm{ }}\int_{\partial {P_i}(t)} {\mathit{\boldsymbol{\sigma }} \cdot \mathit{\boldsymbol{n}}{\rm{d}}S} {\rm{ }}$$ (7) $$\frac{{{\rm{d}}({I_i}{\mathit{\boldsymbol{\omega }}_i})}}{{{\rm{d}}t}} = {\mathit{\boldsymbol{T}}_{{\rm{ext}},i}} + {\rm{ }}\int_{\partial {P_i}(t)} {(\mathit{\boldsymbol{x}} - {\mathit{\boldsymbol{X}}_i}) \times (\mathit{\boldsymbol{\sigma }} \cdot \mathit{\boldsymbol{n}}){\rm{ d}}S} $$ (8) 式中， $m_{i}$ 和 ${ I}_{i}$ 分别是第 $i$ 个固粒的质量和转动惯量， ${\pmb U}_{i}$ 和 ${\pmb \omega}_{i}$ 分别是第 $i$ 个固粒的平动速度和角速度， ${\pmb F}_{{\rm ext}, i}$ 和 ${\pmb T}_{{\rm ext}, i}$ 分别是作用在第 $i$ 个固粒的外力和外力矩， ${\pmb x}$ 是在边界 $\partial P_{i}(t)$ 上某个点的位置矢量， ${\pmb X}_{i}$ 是固粒的质心坐标， ${\pmb n}$ 是在边界 $\partial P_{i}(t)$ 处的单位法向矢量，d $S$ 是局部表面积.不同于牛顿流体情形，非牛顿流体总应力张量 ${\pmb \sigma }$ 的复杂性使得方程(7) 和方程(8) 的求解比较困难.
基于单向耦合模式，联立式(1) ~式(4) 和式(7) ~式(8) 以及相应的本构方程，基于一定的边界条件和初始条件，便可以求解得到非牛顿流体的速度、压力、应力以及固粒的平动速度和角速度，进而得到固粒的轨迹和分布.
1.3 特征参数
对于非牛顿流体固粒悬浮流，除了常用的雷诺数外，还有另外一些重要的特征参数. Deborah数表示流体的松弛时间和流场的特征时间之比；黏度比表示流体黏度与总的零剪切黏度之比；密度比是固粒的密度与流体的密度之比；重力数是作用在固粒的沉降力与黏性力之比；Weissenberg数是第一法向应力差和剪切应力之比.
2. 单个固粒在管道中的径向移动
单个固粒在管道中的径向移动的早期研究主要关注固粒在非牛顿流中运动时的阻力系数，研究发现固粒运动带动流体剪切运动，导致流体的黏度变化从而改变了阻力系数[5].固粒在管道中的径向移动是固粒在流体中运动的典型问题，具有明确的应用背景，例如在固粒分离微流控技术中，通过添加聚合物使流体具有黏弹性，进而影响固粒的运动特性来提高固粒的分离效率.对于牛顿流体，流体的惯性会使固粒沿径向迁移到管道中心和壁面之间的某个位置，不同大小的固粒具有不同的迁移速度和径向平衡位置.近年来利用该特性发展了“惯性聚焦”固粒分离微流控技术[6].然而，对于亚微米和纳米尺度的固粒，极小的固粒雷诺数导致固粒的迁移速度非常小而分离效果不明显.为此，通过在流体中添加聚合物，利用流体的黏弹性来加速固粒向管道中心的迁移，从而改进分离效果[7-9].因此，固粒在管流中的径向移动一直是研究的热点.
在“惯性聚焦”分离技术中，矩形管道比较普遍(图 1). Yang等[8]的实验结果表明，当流体惯性可忽略时，固粒在黏弹性流体中会迁移到管的中心和4个角落; 当流体惯性影响弱时，固粒均向中心迁移. Lim等[10]由实验发现，固粒在弱弹性流体中，即使流体惯性很强，固粒也向管道中心迁移，而流体的剪切变稀特性驱使固粒偏离中心线[11]. Villone等[12]研究了可忽略流体惯性时黏弹性和二次流对固粒迁移的影响, 发现PTT黏弹性流体的固粒会向中心线和附近壁面迁移，并猜测迁移到壁面后的固粒会进一步向角落迁移，迁移方向取决于初始位置和流体弹性大小. Li等[13]研究黏弹性和惯性对固粒迁移的影响发现, 流体弹性驱使固粒向中心迁移，而流体惯性驱使固粒向壁面迁移.
以上研究结果给出了一些圆球形固粒在管流中径向移动的结论，但仍有一些问题尚未弄清楚，例如, 固粒停留在管道角落的平衡点消失时所对应的临界雷诺数值大小以及管道截面的长宽比对固粒平衡点位置产生的影响.理论上预测, 雷诺数较高时，固粒在黏弹性流体中会存在角落平衡点，但实验研究没有发现这样的平衡点[10]，此时固粒在黏弹性流体中的平衡位置需要进一步确定.
3. 多固粒的相互作用和聚集
多个固粒在非牛顿流体中的相互作用和聚集对实际应用有直接影响，例如石油工业中的钻井液和水力压裂液操作时都需要输运固粒，掌握固粒在输运中的聚集特性有助于优化设计，提高石油开采效率.固粒的聚集还会影响固粒分布的均匀性和减阻传热效果等.实验发现，在黏弹性剪切变稀流体中，当间距小于临界距离时竖直放置的两个球形固粒会相互靠近；当间距大于临界距离时，则会相互排斥.该排斥可能和负尾流(即沉降固粒尾流场中速度向上)有关[14].在无弹性强剪切变稀流体中，两个相距较远的球形固粒也会聚集(图 2) [15]，前一固粒通过时产生的低黏度通道被认为是后一固粒加速的原因.水平放置的两个球形固粒在黏弹性剪切变稀流体中会相互靠近，两者转到竖直方向后稳定沉降[16]；在无弹性剪切变稀流体中的固粒也被观察到会径向靠近，但转到竖直位置后将分开[17].实验还发现，多个球形固粒在黏弹性剪切变稀流体[18]和无弹性剪切变稀流体[19]中会形成沉降方向上的固粒串.
以上仅是实验研究结果，但在理论和数值模拟方面的研究还存在一些问题.例如，基于黏弹性本构方程的数值模拟能够重现固粒的聚集和成串现象[20]，基于触变性剪切变稀模型的数值模拟也能重现两固粒在相距很远时的纵向集聚[21]，而流体弹性和剪切变稀同时存在时固粒的聚集模式尚不清楚.实验结果显示，固粒在近似无弹性的剪切变稀流体中存在聚集的现象[19]，但在无剪切变稀的黏弹性Boger流体中却不会聚集[17]，这些结果目前还无法从数值模拟中得到，这是由于尽管数值模拟中的本构模型符合剪切流中所测量的流变特性，但未必能确切描述聚合物在复杂流场下微结构的变化.
4. 多固粒形成的链状结构
多个固粒在非牛顿流体中形成链状结构是非牛顿多相流中重要的现象之一.所谓链状结构, 即若干固粒沿流向形成稳定的队列结构，就像链条一样，固粒间距很小且作整体运动(图 3).固粒形成链状结构后将显著改变流动的宏观特性，例如热传导性显著提高[22]，也能被用来对亚微米和纳米量级的固粒进行更有效的收集、分类与分离[23-25].
Michele等[26]首次由实验发现，黏弹性流体中的固粒在剪切作用下会出现链状结构，Giesekus[27]认为，固粒间润滑流动减弱了流体的运动速度及剪切率，法向应力差在固粒周围形成差异，从而导致固粒相互靠近并形成链状结构.
Feng等[29]则认为, 弹性引起的法向应力差改变了固粒附近的压力分布，从而导致链状结构的形成. Won等[30]由实验发现，在Boger流体中的固粒虽然会沿流向聚集但不会形成链状结构，而在剪切变稀黏弹性流体中，则存在链状结构，且该结构的特征与Weissenber数有关. Pasquino等[31-34]由实验发现，剪切率、固粒浓度、平板间距和固粒尺度是影响链状结构形成速率及链条长度的主要因素. Nie和Lin[35]由数值模拟发现了剪切流中三固粒在幂律流体中的成链过程. Van Loon等[36]由实验发现，在纯剪切变稀流体中也有微弱的链状结构，据此推断剪切变稀是形成链状结构的必要条件，而流体的弹性仅起到促进作用，并非必要条件，他们的这一推断至今没有理论解释. Hwang等[37]由数值模拟发现，在剪切变稀黏弹性流体中存在着与Weissenber数有关的链状结构，这一结论后来被实验所证实[34].他们在Oldroyd-B流体中也模拟到较弱的链状结构，并据此推断流体的弹性是形成链状结构的必要条件，而剪切变稀虽有助于形成链状结构，但其决定性的作用还不明确.这一结论与Jaensson等[38]的数值模拟结果一致，但与实验结果[30, 34, 36]和数值模拟结果[34]不一致，其原因目前没有合理的解释.
由上可见，固粒在非牛顿流体中形成链状结构的机理尚不清楚，是流体的弹性还是剪切变稀起主要作用还存在争议.目前尚未建立链状结构的长度、位置以及取向等特征与流场特性和固粒属性的定量关系.此外，以往的数值模拟大多局限于二维情形且忽略了流体和固粒的惯性，因而无法与实验结果进行比较.在对非牛顿流体中固粒形成链状结构的直接数值模拟中，需引进流体的本构方程，该方程是非线性偏微分方程，求解复杂且耗时，目前还没有兼顾精度、效率和稳定性的数值模拟方法.
5. 非圆球固粒的运动
以往的研究大多数集中在圆球固粒上，然而，许多自然现象和实际应用中的固粒是非圆球固粒.非圆球固粒在形状上为非各向同性，固粒的取向和转动与其平动高度耦合，这使得问题变得非常复杂.非圆球固粒在牛顿流体中的运动已有不少研究结果[39-41]，但在非牛顿流体中运动的研究则较少. Iso等[42-43]由实验发现，剪切流中流体的弱弹性使圆柱形固粒的取向转向涡量轴，而强弹性则使固粒转到流动-剪切平面. Gunes等[44]由实验发现，剪切流中的椭球固粒当长轴取向为涡量轴方向时，固粒则不会集聚，而当长轴转到流动-剪切平面内时，固粒则会聚集，且聚集结构与球形固粒的结构不同.椭球形固粒与流向呈约30 $^\circ$ 的角度, 组合成固粒柱. Bartram等[45]的研究表明，非圆球固粒在黏弹性流体中的旋转周期比圆球固粒长，且在较低的剪切率下将偏离Jeffery轨道. Johnson等[46]发现，对于Boger流体，由于存在非牛顿流体法向应力，椭球形的红血球细胞将偏离牛顿流体时的运动轨道.当作用在细胞上的弹性力强于布朗运动的随机力时，细胞的主轴将偏离流动方向而转到旋涡的轴向；而当随机力强于弹性力时，细胞主轴则更有可能指向流动方向. Gunes等[47]研究了黏弹性流体中由流动导致的椭球形固粒的取向变化.研究发现，当剪切率增大时，固粒的取向由随机转向Jeffery轨道；当Peclet数和Weissenberg数增加时，固粒的旋转周期变长，此时固粒的轨道逐渐呈现滚动状态，导致固粒的取向往旋涡的轴向集中，这一过程的持续时间与流场的剪切率成反比；当流体的弹性更大时，固粒又将重新指向流动方向，而在Boger流体中没有发现以上现象. Leal [48]研究了轴对称柱状颗粒通过一静止二阶黏弹性流体的运动以及在简单剪切流场中的转动.研究发现，在一级近似下，具有任意指向的前后对称的固粒, 其平动速率与在相当的牛顿流体中的情形相同，但是前后不对称的固粒则没有这样的结果. Harlen和Koch[49]分析了在高Deborah数情况下圆柱状固粒在稀聚合物溶液中的运动特性.研究发现，当Deborah数足够高时，在固粒速度扰动和平均剪切流动的共同作用下，会导致聚合物拉伸效应显著加强；非牛顿流体的应力对固粒的角速度产生扰动，使得固粒偏离Jeffery轨道而将主轴指向旋涡的轴向.与黏弹性二阶流体不同的是，这一效应不依赖于第二法向应力差. Phan-Thien等[50-51]在小Deborah情形下，通过扁球体固粒在Oldroyd-B剪切流中的运动，研究了流体的黏弹性效应.结果发现，黏弹性应力将使固粒的旋转变缓，当应力增大时，固粒将偏离Jeffery轨道. D'Avino等[52]研究了椭球固粒在黏弹性流体自由剪切流场中的运动，研究发现，在低Deborah数时，固粒由盘旋轨道转变成绕旋涡轴的滚动轨道，而在高Deborah数时，固粒的主轴则指向流动方向.
尽管非圆球固粒在非牛顿流体中的运动已有一些研究结果，但仍然有一些问题未得到研究或澄清.例如：圆柱状固粒在非牛顿流体管流中径向移动的研究还未见报道；椭球固粒为何与圆球固粒的聚集结构不同、圆柱形固粒在非牛顿流体运动中相互作用的机理与结果尚不清楚；非圆球固粒在随时间演化的非牛顿流场中的运动特性以及Stokes数、Deborah数、固粒与流体的密度比等参数对于流场以及固粒的运动与分布影响尚不清楚.
6. 总结和展望
以上对单个固粒在管道中的径向移动、多固粒的相互作用和聚集、多固粒形成的链状结构以及非圆球固粒运动等方面的研究成果进行了阐述与分析.可见，以往在这几个方面虽取得了一定的研究成果，但仍有一些问题没有得到解决.
对于单个固粒在管道中的径向移动，可以开展的研究包括：① 圆球和椭球形固粒在矩形管道内的径向迁移；② 不同雷诺数下，固粒停留在矩形管道角落的平衡点存在与消失的特征；③ 流变特性、雷诺数、矩形管道截面长宽比等因素对固粒径向迁移特性和平衡点位置的影响；④ 探讨较高雷诺数下固粒在黏弹性流体中是否有平衡点及其平衡点的位置.
在多固粒的相互作用和聚集方面，可以开展的研究包括：① 多个圆球形固粒在黏弹性Couette流中的运动规律和聚集结构；② 在两竖直圆筒间侧向剪切流中，多个圆球形固粒在纯剪切变稀、黏弹性剪切变稀、无剪切变稀黏弹性流体中的运动规律和聚集结构；③ 流体弹性和剪切变稀对固粒聚集所起的作用以及两者同时存在时固粒的聚集模式；④ 多个椭球形固粒、圆柱形固粒在黏弹性Couette流中的运动规律和聚集结构以及流变特性和颗粒长径比对聚集的影响.
在多固粒形成的链状结构方面，可以开展的研究包括：① 将用于牛顿流体的格子Boltzmann虚拟区域方法与非牛顿流体本构方程结合，探索用于非牛顿流体的格子Boltzmann虚拟区域方法；② 研究Couette流中多固粒在纯剪切变稀、黏弹性剪切变稀、无剪切变稀黏弹性流体中形成链状结构的机理；③ 了解Couette流中流变特性、惯性、颗粒浓度、管道与固粒尺度比等因素对链状结构的长度、位置以及取向的影响；④ 给出三维情况下考虑流体和固粒的惯性时，链状结构形成的条件.
-
[1] 宋晓阳, 及春宁, 许栋.明渠湍流边界层中颗粒的运动及分布.力学学报, 2015, 47(2):231-241 doi: 10.6052/0459-1879-14-164 Song Xiaoyang, Ji Chunning, Xu Dong. Distribution and motion of particles in the turbulent boundary layer of channel flow. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(2):231-241 (in Chinese) doi: 10.6052/0459-1879-14-164
[2] 丁珏, 李家骅, 邱骁等.蒙特卡洛方法数值研究大气颗粒物动力学效应和辐射传输性质.力学学报, 2016, 48(3):557-565 http://lxxb.cstam.org.cn/CN/abstract/abstract145833.shtml Ding Jue, Li Jiahua, Qiu Xiao, et al. Numerical study on dynamics effect and rediation transfer characteristics of atmospheric particle by Monte Carlo method. Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(3):557-565 (in Chinese) http://lxxb.cstam.org.cn/CN/abstract/abstract145833.shtml
[3] 王帅, 于文浩, 陈巨辉等.鼓泡流化床中流动特性的多尺度数值模拟.力学学报, 2016, 48(3):585-592 http://lxxb.cstam.org.cn/CN/abstract/abstract145781.shtml Wang Shuai, Yu Wenhao, Chen Juhui, et al. Multi-scale simulation on hydrodynamic characteristics in bubbling fluidized bed. Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(3):585-592 (in Chinese) http://lxxb.cstam.org.cn/CN/abstract/abstract145781.shtml
[4] D'Avino G, Maffettone PL. Particle dynamics in viscoelastic liquids. Journal of Non-Newtonian Fluid Mechanics, 2015, 215:80-104 doi: 10.1016/j.jnnfm.2014.09.014
[5] Chhabra RP. Bubbles, Drops and Particles in Non-Newtonian Fluids. New York:Taylor & Francis Group, 2007
[6] Di Carlo D, Irimia D, Tompkins RG, et al. Continuous inertial focusing, ordering, and separation of particles in microchannels. Proceedings of the National Academy of Sciences, 2007, 104 (48):18892-18897 doi: 10.1073/pnas.0704958104
[7] Leshansky A, Bransky A, Korin N, et al. Tunable nonlinear viscoelastic "focusing" in a microfluidic device. Phys Rev Lett, 2007, 98(23):234501 doi: 10.1103/PhysRevLett.98.234501
[8] Yang S, Kim JY, Lee SJ, et al. Sheathless elasto-inertial particle focusing and continuous separation in a straight rectangular microchannel. Lab Chip, 2011, 11(2):266-273 doi: 10.1039/C0LC00102C
[9] Karimi A, Yazdi S, Ardekani AM. Hydrodynamic mechanisms of cell and particle trapping in microfluidics. Biomicrofluids, 2013, 7(2):021501 doi: 10.1063/1.4799787
[10] Lim EJ, Ober TJ, Edd JF, et al. Inertio-elastic focusing of bioparticles in microchannels at high throughput. Nature Communications, 2014, 5:4120 http://people.seas.harvard.edu/~tober/Publications_files/TJO9.pdf
[11] Liu C, Xue C, Chen X, Shan L, et al. Size-based separation of particles and cells utilizing viscoelastic effects in straight microchannels. Analytical Chemistry, 2015, 87(12):6041-6048 doi: 10.1021/acs.analchem.5b00516
[12] Villone MM, Avino GD', Hulsen MA, et al. Particle motion in square channel flow of a viscoelastic liquid:migration vs. secondary flows. Journal of Non-Newtonian Fluid Mechanics, 2013, 195:1-8 doi: 10.1016/j.jnnfm.2012.12.006
[13] Li G, McKinley GH, Ardekani AM. Dynamics of particle migration in channel flow of viscoelastic fluids. Journal of Fluid Mechanics, 2015, 785:486-505 doi: 10.1017/jfm.2015.619
[14] Riddle MJC, Narvaez C, Bird RB. Interactions between two spheres falling along their line of centers in a viscoelastic fluid. Journal Non-Newtonian Fluid Mechanics, 1977, 2(17):23-35 http://www.academia.edu/10573503/Interactions_between_two_spheres_falling_along_their_line_of_centers_in_a_viscoelastic_fluid
[15] Daugan S, Talini L, Herzhaft B, et al. Aggregation of particles settling in shear-thinning fluids. Part 1. Two-particle aggregation. The European Physical Journal E, 2002, 7(1):55 https://www.researchgate.net/publication/237430629_Aggregation_of_particles_settling_in_shear-thinning_fluids_Part_1_Two-particle_aggregation
[16] Joseph DD, Liu YJ, Poletto M, et al. Aggregation and dispersion of spheres falling in viscoelastic liquids. Journal of Non-Newtonian Fluid Mechanics, 1994, 54(6):45-86 https://www.researchgate.net/publication/222439800_Aggregation_and_dispersion_of_spheres_falling_in_viscoelastic_liquids
[17] Gheissary G, van den Brule BHAA. Unexpected phenomena observed in particle settling in non-Newtonian media. Journal of NonNewtonian Fluid Mechanics, 1996, 67(1):1-18 https://www.researchgate.net/publication/223679532_Unexpected_phenomena_observed_in_particle_settling_in_non-Newtonian_media
[18] Bobroff S, Phillips R. Nuclear magnetic resonance imaging investigation of sedimentation of concentrated suspensions in nonNewtonian fluids. Journal of Rheology, 1998, 42(1-2):1419-1436 https://www.researchgate.net/publication/243529792_Nuclear_magnetic_imaging_investigation_of_sedimentation_of_concentrated_suspensions_in_non-Newtonian_fluids
[19] Daugan S, Talini L, Herzhaft B, et al. Sedimentation of suspensions in shear-thinning fluids. Oil & Gas Science Technology, 2004, 59(1):71-80 https://www.researchgate.net/publication/241350216_Sedimentation_of_Suspensions_in_Shear-Thinning_Fluids
[20] Yu ZS, Phan-Thien N, Fan Y, et al. Viscoelastic mobility problem of a system of particles. Journal of Non-Newtonian Fluid Mechanics, 2002, 104:87-124 doi: 10.1016/S0377-0257(02)00014-9
[21] Yu ZS, Wachs A, Peysson Y. Numerical simulation of particle sedimentation in shear-thinning fluids with a fictitious domain method. Journal of Non-Newtonian Fluid Mechanics, 2006, 136(2-3):126-139 doi: 10.1016/j.jnnfm.2006.03.015
[22] Vaia R, Giannelis EP. Polymer nanocomposites:status and opportunities. MRS Bull, 2001, 26(5):394-401 doi: 10.1557/mrs2001.93
[23] Sun X, Tabakman SM, Seo WS, et al. Separation of nanoparticles in a density gradient:FeCo@C and gold nanocrystals. Angewandte Chemie International Edition Engl, 2009, 48(5):939-942 doi: 10.1002/anie.v48:5
[24] Hao J, Pan TW, Glowinski R, et al. A fictitious domain/distributed Lagrange multiplier method for the particulate flow of Oldroyd-B fluids:a positive definiteness preserving approach. Journal of NonNewton Fluid Mechanics, 2009, 156(1):95-111 https://www.researchgate.net/publication/223099193_A_fictitious_domaindistributed_Lagrange_multiplier_method_for_the_particulate_flow_of_Oldroyd-B_fluids_A_positive_definiteness_preserving_approach
[25] Surendra Balaji Devarakonda, Han J, Ahn CH, et al. Bioparticle separation in non-Newtonian fluid using pulsed flow in micro-channels. Microfluidics and Nanofluidics, 2007, 3(4):391-401 doi: 10.1007/s10404-006-0131-6
[26] Michele J, Padzold R, Donis R. Alignment and aggregation effects in suspensions of spheres in non-Newtonian media. Rheologica Acta, 1977, 16(3):317-321 doi: 10.1007/BF01523742
[27] Giesekus H. Particle movement in flows of non-Newtonian fluids. Z Angew Math Mechanics, 1978, 58:T26-T37
[28] Van Loon S, Fransaer J, Clasen C, et al. String formation in sheared suspensions in rheologically complex media:the essential role of shear thinning. Journal of Rheology, 2014, 58(1):237-254 doi: 10.1122/1.4853455
[29] Feng J, Huang PY, Joseph DD. Dynamic simulation of sedimentation of solid particles in an Oldroyd-B fluid. Journal of Non-Newtonian Fluid Mechanics, 1996, 63(1):63-88 doi: 10.1016/0377-0257(95)01412-8
[30] Won D, Kim C. Alignment and aggregation of spherical particles in viscoelastic fluid under shear flow. Journal of Non-Newtonian Fluid Mechanics, 2004, 117(2-3):141-146 doi: 10.1016/j.jnnfm.2004.01.005
[31] Pasquino R, D'Avino G, Maffettone PL, et al. Migration and chaining of noncolloidal spheres suspended in a sheared viscoelastic medium. Experiments and numerical simulations. Journal of NonNewtonian Fluid Mechanics, 2014, 203:1-8 doi: 10.1016/j.jnnfm.2013.10.006
[32] Pasquino R, Snijkers F, Grizzuti N, et al. Directed self-assembly of spheres into a two-dimensional colloidal crystal by viscoelastic stresses. Langmuir, 2010, 26(5):3016-3019 doi: 10.1021/la904775c
[33] Pasquino R, Snijkers F, Grizzuti N, et al. The effect of particle size and migration on the formation of flow-induced structures in viscoelastic suspensions. Rheologica Acta, 2010, 49(10):993-1001 doi: 10.1007/s00397-010-0466-5
[34] Pasquino R, Panariello D, Grizzuti N. Migration and alignment of spherical particles in sheared viscoelastic suspensions. A quantitative determination of the flow-induced self-assembly kinetics. Journal of Colloid Interface Science, 2013, 394(4):49-54 https://www.researchgate.net/publication/233983504_Migration_and_alignment_of_spherical_particles_in_sheared_viscoelastic_suspensions_A_quantitative_determination_of_the_flow-induced_self-assembly_kinetics
[35] Nie DM, Lin JZ. Behavior of three circular particles in a confined power-law fluid under shear. Journal of Non-Newtonian Fluid Mechanics, 2015, 221:76-94 doi: 10.1016/j.jnnfm.2015.04.004
[36] Van Loon S, Fransaer J, Clasen C, et al. String formation in sheared suspensions in rheologically complex media:the essential role of shear thinning. Journal of Rheol, 2014, 58(1):237-254 doi: 10.1122/1.4853455
[37] Hwang W, Hulsen MA. Structure formation of non-colloidal particles in viscoelastic fluids subjected to simple shear flow. Macromolecular Materials & Engineering, 2011, 296(3-4):321-330 https://www.researchgate.net/publication/228046397_Structure_Formation_of_Non-Colloidal_Particles_in_Viscoelastic_Fluids_Subjected_to_Simple_Shear_Flow
[38] Jaensson NO, Hulsena MA, Anderson PD. Simulations of the startup of shear flow of 2D particle suspensions in viscoelastic fluids:structure formation and rheology. Journal of Non-Newtonian Fluid Mechanics, 2015, 225:70-85 doi: 10.1016/j.jnnfm.2015.09.006
[39] Kazi SN, Duffy GG, Chen XD. Validation of heat transfer and friction loss data for fibre suspensions in a circular and a coaxial pipe heat exchanger. International Journal of Thermal Sciences, 2014, 79(4):146-160 https://www.researchgate.net/profile/G_Duffy2/publication/260212861_Validation_of_heat_transfer_and_friction_loss_data_for_fibre_suspensions_in_a_circular_and_a_coaxial_pipe_heat_exchanger/links/0deec531164ee9202e000000.pdf
[40] Lin JZ, Xia Y, Ku XK. Flow and heat transfer characteristics of nanofluids containing rod-like particles in a turbulent pipe flow. International Journal of Heat and Mass Transfer, 2016, 93(1):57-66 https://www.researchgate.net/publication/283801931_Flow_and_heat_transfer_characteristics_of_nanofluids_containing_rod-like_particles_in_a_turbulent_pipe_flow
[41] Shao XM, Zhang XL, Yu ZS, et al. Numerical studies on the dynamics of an open triangle in a vertically oscillatory flow. Journal of Fluid Mechanics, 2016, 788:381-406 doi: 10.1017/jfm.2015.703
[42] Iso Y, Koch DL, Cohen C. Orientation in simple shear flow of semidilute fiber suspensions 1. Weakly elastic fluids. Journal NonNewtonian Fluid Mech, 1996, 62(2-3):115-134 doi: 10.1016/0377-0257(95)01404-7
[43] Iso Y, Koch DL, Cohen C. Orientation in simple shear flow of semidilute fiber suspensions 2. Highly elastic fluids. Journal NonNewtonian Fluid Mech, 1996, 62(2-3):135-153 doi: 10.1016/0377-0257(95)01405-5
[44] Gunes DZ, Scirocco R, Mewis J, et al. Flow-induced orientation of nonspherical particles:effect of aspect ratio and medium rheology. Journal of Non-Newtonian Fluid Mechanics, 2008, 155(1-2):39-50 doi: 10.1016/j.jnnfm.2008.05.003
[45] Bartram E, Goldsmith HL, Mason SG. Particle motions in non-newtonian media llI. Further observations in elasticoviscous fluids. Rheological Acta, 1971, 14(9): 776-782 https://www.researchgate.net/publication/227016835_Particle_motions_in_non-Newtonian_media
[46] Johnson SJ, Salem NJ, Fuller GG. Dynamics of colloidal particles in sheared, non-Newtonian fluids. Journal of Non-Newtonian Fluid Mechanics, 1990, 34(1): 89-121 doi: 10.1016/0377-0257(90)80013-P
[47] Gunes DZ, Scirocco R, Mewis J, et al. Flow-induced orientation of nonspherical particles: effect of aspect ratio and medium rheology. Journal of Non-Newtonian Fluid Mechanics, 2008, 155(1-2): 39-50 doi: 10.1016/j.jnnfm.2008.05.003
[48] Leal LG. The slow motion of slender rod-like particles in a second-order Fluid. Journal of Fluid Mechanics, 1975, 69(2): 305-337 doi: 10.1017/S0022112075001450
[49] Harlen OG, Koch DL. Simple shear-flow of a suspension of fibers in a dilute polymer-solution at high Deborah number. Journal of Fluid Mechanics, 1993, 252: 187-207 doi: 10.1017/S0022112093003726
[50] Phan-Thien N, Fan XJ. Viscoelastic mobility problem using a boundary element method. Journal of Non-Newtonian Fluid Mechanics, 2002, 105(2-3): 131-152 doi: 10.1016/S0377-0257(02)00079-4
[51] Nguyen-Hoang H, Phan-Thien N, Khoo BC. et al. Completed double layer boundary element method for periodic fibre suspension in viscoelastic fluid. Chemical Engineering Science, 2008, 63(15): 3898-3908 doi: 10.1016/j.ces.2008.04.058
[52] D'Avino G, Hulsen MA, Greco F, et al. Bistability and metabistability scenario in the dynamics of an ellipsoidal particle in a sheared viscoelastic fluid. Physical Review E, 2014, 89(4): 043006 doi: 10.1103/PhysRevE.89.043006
-
期刊类型引用(0)
其他类型引用(2)
下载:
计量
- 文章访问数: 1542
- HTML全文浏览量: 162
- PDF下载量: 839
- 被引次数: 2
