CN110465339B - 一种流固两相输运中颗粒定位的方法 - Google Patents
一种流固两相输运中颗粒定位的方法 Download PDFInfo
- Publication number
- CN110465339B CN110465339B CN201910828574.0A CN201910828574A CN110465339B CN 110465339 B CN110465339 B CN 110465339B CN 201910828574 A CN201910828574 A CN 201910828574A CN 110465339 B CN110465339 B CN 110465339B
- Authority
- CN
- China
- Prior art keywords
- solid particles
- particle
- intervals
- distributed
- solid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 336
- 239000007787 solid Substances 0.000 title claims abstract description 237
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000012530 fluid Substances 0.000 claims abstract description 53
- 230000000903 blocking effect Effects 0.000 claims abstract description 27
- 239000012071 phase Substances 0.000 claims abstract description 22
- 230000009471 action Effects 0.000 claims abstract description 9
- 239000007791 liquid phase Substances 0.000 claims abstract description 8
- 238000002347 injection Methods 0.000 claims abstract description 6
- 239000007924 injection Substances 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims description 19
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 6
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 239000007790 solid phase Substances 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000004807 localization Effects 0.000 claims 4
- 238000004364 calculation method Methods 0.000 description 25
- 238000006073 displacement reaction Methods 0.000 description 25
- 238000004088 simulation Methods 0.000 description 13
- 238000000926 separation method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/50273—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip characterised by the means or forces applied to move the fluids
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/502753—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip characterised by bulk separation arrangements on lab-on-a-chip devices, e.g. for filtration or centrifugation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/502761—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip specially adapted for handling suspended solids or molecules independently from the bulk fluid flow, e.g. for trapping or sorting beads, for physically stretching molecules
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/502769—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip characterised by multiphase flow arrangements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L2200/00—Solutions for specific problems relating to chemical or physical laboratory apparatus
- B01L2200/14—Process control and prevention of errors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Clinical Laboratory Science (AREA)
- Hematology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Fluid Mechanics (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
- Automatic Analysis And Handling Materials Therefor (AREA)
Abstract
本发明公开了一种流固两相输运中颗粒定位的方法。在微流控芯片的微通道内，用压力源从入口端将微通道注入充满液相牛顿流体，将两排固体颗粒并排放置在微通道的入口端，且位于流向中线对称的两侧，两排固体颗粒在入口段驱动压力作用下均逐颗注入到微流控芯片的微通道内，初始注入速度为0，间隙/距离相同；调节压力源的驱动压力和颗粒直径，使得两排固体颗粒分别均匀分布在微流控芯片的微通道两排内壁的壁面附近，形成颗粒间距稳定的交错分布的颗粒链。本发明根据雷诺数和固体颗粒的阻塞率优化了形成均匀稳定的颗粒间距和排布，能够在无外力条件下就形成均匀间距的颗粒链，达到精确定位和计量颗粒的目的，方法提高了颗粒定位的精度和效率。
Description
技术领域
本发明属于微流控芯片中固液两相流领域，特别涉及一种流固两相输运中颗粒定位的方法。
背景技术
颗粒惯性迁移是微流控分析系统中一种重要的分离技术。如细胞仪对血液中血细胞的检测计量，污染颗粒物在呼吸道中的运动，化工过程中颗粒的筛选与分离，微流控芯片中颗粒的高效分离输运等，具有重要的应用价值。以往的研究中多采用施加外力的方式，如磁场、电场以及复杂的微通道结构来实现颗粒在微流控芯片中会形成均匀间距的颗粒链，提高细胞仪对颗粒定位的精度。而颗粒在压力驱动的惯性作用下，能够在无其他外力条件，在简单的直管中固体颗粒就能惯性迁移到指定位置，形成均匀稳定的颗粒间距，极大的提高颗粒监测的准确度和效率。因此，研究颗粒在牛顿流体中的定位方法对提高颗粒检测精度和分离效率具有重要意义。
发明内容
本发明的目的是针对现有技术的不足，提供一种流固两相输运中颗粒定位的方法，还根据雷诺数和固体颗粒的阻塞率优化了形成均匀稳定的颗粒间距和排布，方法提高了颗粒定位的精度和效率。
本发明的目的是通过以下技术方案来实现的：
在微流控芯片的微通道内，微通道的入口端连接压力源，用压力源从入口端将微通道注入充满液相牛顿流体，将两排固体颗粒并排放置在微通道的入口端，且位于微通道的流向中线对称的两侧，两排固体颗粒平行于微通道且在入口端之外，两排固体颗粒到微通道的流向中线的径向距离相同，分别到微通道两侧内壁之间的距离相同，
两排固体颗粒在入口段驱动压力作用下均逐颗注入到微流控芯片的微通道内，固体颗粒的初始注入速度为0，两排固体颗粒中保持流向注入的相邻两颗固体颗粒的间隙/距离相同；
调节压力源的驱动压力和颗粒直径，固体颗粒受固定的驱动压力所引起的惯性和相邻固体颗粒相互靠近时的排斥力作用而向微通道下游运动，使得两排固体颗粒分别均匀分布在微流控芯片的微通道两排内壁的壁面附近，形成颗粒间距稳定的交错分布的颗粒链。
所述的微通道为直通道。
所述的颗粒间距稳定的交错分布的颗粒链，是指两排固体颗粒的各颗固体颗粒在流动方向交错交替分布，并且：
其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
或者是，两排固体颗粒中的固体颗粒均以相同间距或者规律间距地间隔分布；
或者是，两排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布。
所述的液相牛顿流体为水或甘油之类。
选择圆形固体颗粒作为固体相，且固体颗粒密度与液相牛顿流体的密度相等。具体实施中，圆形固体颗粒可以为细胞等。
所述的压力源采用压力泵或注射泵。
所述的调节压力源的驱动压力和颗粒直径，具体是：调节压力源的驱动压力使得雷诺数为14、82～120，雷诺数(Re＝ρUmaxH/m，ρ为牛顿流体的密度，Umax为驱动压力产生的速度，H为微流控芯片的直通道高度，m为牛顿流体的粘度)，调节颗粒直径使得阻塞率为0.125～0.4，阻塞率(k＝D/H，D为固体颗粒直径，H为微流控芯片的直通道高度)，使得固体颗粒间的间距能达到均匀稳定，形成固定的颗粒链。
当调节压力源的驱动压力使得雷诺数为14且调节颗粒直径使得阻塞率为0.125～0.3时，形成以下颗粒间距稳定的交错分布的颗粒链：其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布。
当调节压力源的驱动压力使得雷诺数为82且调节颗粒直径使得阻塞率为0.2～0.3时，形成以下颗粒间距稳定的交错分布的颗粒链：其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布。
当调节压力源的驱动压力使得雷诺数为120且调节颗粒直径使得阻塞率为0.125～0.4时，形成以下颗粒间距稳定的交错分布的颗粒链：两排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
并且相邻两固体颗粒的间距为正弦变化，间距变化随雷诺数的增加而增加，相邻两固体颗粒的间距随下游流向越远而越大。
本发明选择微流控芯片的直通道为牛顿流体和固体颗粒运动的场所，驱动压力的压力源通常采用压力泵或注射泵，微通道入口段用压力驱动，将牛顿流体先注满微流控芯片的直通道。
本发明比较相邻两固体颗粒的间距随固体颗粒沿流向位移的变化，获得了固体颗粒间距能达到均匀稳定时所对应驱动压力产生的情形。
基于上述技术方案，本发明实施例至少可以产生如下技术效果：
本发明仅控制驱动压力和固体颗粒直径，无需其他外力条件下，固体颗粒在简单的微流控芯片的直通道中就能运动到指定位置，得到均匀稳定的颗粒间距，达到对颗粒位置准确定位的目的，极大地提高了流固两相输运中颗粒定位的精度和效率，可以用于极大的提高了细胞仪对细胞检测的精度，提高效率和节约成本，为在微流控芯片上实现高效计数、分离提供了简单有效的方法。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明，其中：
图1是并排的颗粒在直管中注入过程示意图；图中，1、微流控芯片的微通道；2、液相牛顿流体；3、固体颗粒；4、压力源。
图2是雷诺数14，阻塞率0.125时，牛顿流体中颗粒间距随位移的变化；
图3是雷诺数14，阻塞率0.2时，牛顿流体中颗粒间距随位移的变化；
图4是雷诺数14，阻塞率0.3时，牛顿流体中颗粒间距随位移的变化；
图5是雷诺数14，阻塞率0.4时，牛顿流体中颗粒间距随位移的变化；
图6是雷诺数82，阻塞率0.125时，牛顿流体中颗粒间距随位移的变化；
图7是雷诺数82，阻塞率0.2时，牛顿流体中颗粒间距随位移的变化；
图8是雷诺数82，阻塞率0.3时，牛顿流体中颗粒间距随位移的变化；
图9是雷诺数82，阻塞率0.4时，牛顿流体中颗粒间距随位移的变化；
图10是雷诺数120，阻塞率0.125时，牛顿流体中颗粒间距随位移的变化；
图11是雷诺数120，阻塞率0.2时，牛顿流体中颗粒间距随位移的变化；
图12是雷诺数120，阻塞率0.3时，牛顿流体中颗粒间距随位移的变化；
图13是雷诺数120，阻塞率0.4时，牛顿流体中颗粒间距随位移的变化；
图14是相邻两固体颗粒的间距随雷诺数和阻塞率的变化。
具体实施方式
以下是本发明的具体实施例并结合附图，对本发明的技术方案做进一步的描述，但本发明并不限于这些实施例。
本发明保护一种流固两相输运中颗粒定位的方法，该方法以微流控芯片中的直通道内的固液两相流为研究对象，通过调整驱动压力和固体颗粒直径，圆形固体颗粒在微通道只受流场作用，无其他外力时，就能自动运动到确定位置，形成均匀颗粒间距的颗粒链，提高了流固两相输运中颗粒定位的精度和效率，有利于细胞检测仪对细胞及颗粒的检测计量。
本发明的具体实施例如下：
具体实施通过格子玻尔兹曼方法，模拟计算中通过驱动压力4将牛顿流体2注满微流控芯片的直通道1中，圆形固体颗粒3并排在微流控芯片的直通道1中线两侧相同高度的入口处，在入口段驱动压力4作用下逐排注入到微流控芯片的直通道1内，在给定的驱动压力4作用下，固体颗粒3在驱动压力4引起的惯性作用向下游运动时，固体颗粒3都均匀分布在微流控芯片的直通道1的壁面附近，固体颗粒3自动形成颗粒间距稳定、位置确定的交错分布颗粒链。
如图1并排的颗粒在直管中注入过程示意图所示。微流控芯片的直通道1的高度为H，长度为L。微流控芯片的直通道1中的流体为牛顿流体3，微流控芯片的直通道1的入口处给予驱动压力4，使牛顿流体2和固体颗粒3运动，初始时刻圆形的固体颗粒3在微流控芯片的直通道1的中线两侧，离微流控芯片的直通道1上下壁面无量纲高度(h/H)分别为0.25和0.75的位置处，从入口静止释放，相邻的固体颗粒3的间距W为两倍的固体颗粒3直径。在驱动压力4作用和固体颗粒3间相互作用下，固体颗粒3都均匀分布在微流控芯片的直通道1的壁面附近，固体颗粒3自动形成颗粒间距稳定、位置确定的交错分布的颗粒链。数值模拟计算中取固体颗粒3位置的无量纲高度y/H，固体颗粒3流向的无量纲运动长度x/H，相邻两固体颗粒3的无量纲颗粒间距dp/D，微流控芯片的直通道1中线以上的固体颗粒3从右向左依次命名为Pu1，Pu2，Pu3，Pu4，Pu5，Pu6，微流控芯片的直通道1中线以下的固体颗粒3从右向左依次命名为Pd1，Pd2，Pd3，Pd4，Pd5，Pd6，无量纲颗粒间距依次为dd1-u1,du1-d2,…,dd6-u6。
具体实施例1
通过格子玻尔兹曼方法模拟计算的结果如图2，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为14，固体颗粒3的阻塞率(k)为0.125，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图2的插图所示，相邻两固体颗粒3的间距起初会波动，但最终大部分固体颗粒3形成了均匀稳定间距的交错分布颗粒链，只有小部分固体颗粒3间距略大，但交错分布颗粒链结构清晰。
具体实施例2
通过格子玻尔兹曼方法模拟计算的结果如图3，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为14，固体颗粒3的阻塞率(k)为0.2，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图3的插图所示，相邻两固体颗粒3的间距起初会波动，但固体颗粒3很快形成了均匀稳定间距的交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.125时低。
具体实施例3
通过格子玻尔兹曼方法模拟计算的结果如图4，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为14，固体颗粒3的阻塞率(k)为0.3，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图4的插图所示，相邻两固体颗粒3的间距起初会波动，但固体颗粒3很快形成了均匀稳定间距的交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.125时低。
具体实施例4
通过格子玻尔兹曼方法模拟计算固体颗粒3的运动过程如图5，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为14，固体颗粒3的阻塞率(k)为0.4，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，相邻两固体颗粒3的间距在该雷诺数时，出现无规律性，固体颗粒3稳定的交错分布结构消失，此时颗粒间相互作用大于流场驱动压力4，其起到了主导作用，所以固体颗粒3的稳定交错结构消失。
具体实施例5
通过格子玻尔兹曼方法模拟计算的结果如图6，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为82，固体颗粒3的阻塞率(k)为0.125，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图6的插图所示，相邻两固体颗粒3的间距出现大幅度波动，但最终固体颗粒3形成了均匀稳定间距的交错分布颗粒链。
具体实施例6
通过格子玻尔兹曼方法模拟计算的结果如图7，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为82，固体颗粒3的阻塞率(k)为0.2，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图7的插图所示。相邻两固体颗粒3的间距在一定范围内出现正弦形状的波动，波动幅度随雷诺数的增加而增加，但此时固体颗粒3仍然能形成交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.125时低。
具体实施例7
通过格子玻尔兹曼方法模拟计算的结果如图8，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为82，固体颗粒3的阻塞率(k)为0.3，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图8的插图所示。相邻两固体颗粒3的间距在一定范围内出现正弦形状的波动，波动幅度随雷诺数的增加而增加，但此时固体颗粒3仍然能形成交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.2时进一步降低。
具体实施例8
通过格子玻尔兹曼方法模拟计算的结果如图9，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为82，固体颗粒3的阻塞率(k)为0.4，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图9的插图所示。相邻两固体颗粒3的间距波动幅度随雷诺数的增加而增加，但随着阻塞率的增加，颗粒间距最终趋于稳定，形成稳定交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.3时进一步降低。
具体实施例9
通过格子玻尔兹曼方法模拟计算的结果如图10，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为120，固体颗粒3的阻塞率(k)为0.125，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图10的插图所示，相邻两固体颗粒3的间距波动幅度进一步增加，此时均匀稳定间距的交错分布颗粒链不再明显。
具体实施例10
通过格子玻尔兹曼方法模拟计算的结果如图11，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为120，固体颗粒3的阻塞率(k)为0.2，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图11的插图所示，相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.125时进一步降低，但固体颗粒3形成交错分布颗粒链依然不明显。
具体实施例11
通过格子玻尔兹曼方法模拟计算的结果如图12，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为120，固体颗粒3的阻塞率(k)为0.3，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图12的插图所示，相邻两固体颗粒3的间距在一定范围内出现正弦形状的波动，波动幅度随雷诺数的增加而增加，但此时固体颗粒3仍然形成了交错分布颗粒链，而且相邻两固体颗粒3的间距波动幅度比相同条件时阻塞率为0.125和0.2时低。
具体实施例12
通过格子玻尔兹曼方法模拟计算的结果如图13，计算时调整驱动压力4和固体颗粒3的直径，使流场的雷诺数(Re)为120，固体颗粒3的阻塞率(k)为0.4，得到牛顿流体2中相邻固体颗粒3的间距随固体颗粒3位移的变化；固体颗粒3在牛顿流体2中分别往靠近固体颗粒3的微流控芯片的直通道1的壁面方向运动，颗粒分布如图13的插图所示，相邻两固体颗粒3的间距出现了正弦形状的波动，波动幅度随雷诺数的增加而增加，但此时固体颗粒3仍然形成了交错分布颗粒链。
本发明专利将固体颗粒3在不同的雷诺数和阻塞率时的间距进行比较，结果如图14所示，颗粒间距随着阻塞率的增加而降低，随雷诺数的增加而增加，而且综合比较具体实施例1-8，可以得到，在低雷诺数时，固体颗粒3的颗粒间距稳定，当雷诺数增加到82和120时，固体颗粒3间距的波动幅度大幅度增加，固体颗粒3间距不再稳定。
此时增加固体颗粒3的直径，使阻塞率增加，固体颗粒3的间距波动幅度稳定到固定范围内。为了达到稳定的颗粒间距，优选出固体颗粒的阻塞率为0.125-0.4时，颗粒间距能达到均匀稳定时所对应驱动压力4产生的雷诺数为14-82。所以可以通过调整驱动压力4和固体颗粒3的直径，得到位置和间距确定的颗粒链，有利于提高细胞仪对颗粒定位和检测的精度。
Claims (5)
1.一种流固两相输运中颗粒定位的方法，其特征在于：
方法在以下条件下进行实施：
在微流控芯片的微通道（1）内，微通道（1）的入口端连接压力源（4），用压力源（4）从入口端将微通道（1）注入充满液相牛顿流体（2），将两排固体颗粒（3）并排放置在微通道的入口端，且位于微通道的流向中线对称的两侧，两排固体颗粒（3）在入口段驱动压力作用下均逐颗注入到微流控芯片的微通道内，固体颗粒的初始注入速度为0，两排固体颗粒中保持流向注入的相邻两颗固体颗粒的间隙/距离相同；
方法处理具体为：
调节压力源的驱动压力和颗粒直径，使得两排固体颗粒分别均匀分布在微流控芯片的微通道两排内壁的壁面附近，形成颗粒间距稳定的交错分布的颗粒链；
选择圆形固体颗粒作为固体相，且固体颗粒（3）密度与液相牛顿流体（2）的密度相等；
所述的颗粒间距稳定的交错分布的颗粒链，是指两排固体颗粒的各颗固体颗粒在流动方向交错交替分布，并且：
其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
或者是，两排固体颗粒中的固体颗粒均以相同间距或者规律间距地间隔分布；
或者是，两排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
当调节压力源的驱动压力使得雷诺数为14且调节颗粒直径使得阻塞率为0.125～0.3时，形成以下颗粒间距稳定的交错分布的颗粒链：其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
当调节压力源的驱动压力使得雷诺数为82且调节颗粒直径使得阻塞率为0.2～0.3时，形成以下颗粒间距稳定的交错分布的颗粒链：其中一排固体颗粒中的固体颗粒以相同间距或者规律间距地间隔分布，另一排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；
当调节压力源的驱动压力使得雷诺数为120且调节颗粒直径使得阻塞率为0.125～0.4时，形成以下颗粒间距稳定的交错分布的颗粒链：两排固体颗粒中的固体颗粒沿流向以连续的几个相靠近聚集为一组而形成不同组，组之间间隔分布，组内的固体颗粒以相同间距或者规律间距地间隔分布；并且相邻两固体颗粒（3）的间距为正弦变化，间距变化随雷诺数的增加而增加，相邻两固体颗粒（3）的间距随下游流向越远而越大。
2.根据权利要求1所述的一种流固两相输运中颗粒定位的方法，其特征在于：所述的微通道为直通道。
3.根据权利要求1所述的一种流固两相输运中颗粒定位的方法，其特征在于：所述的液相牛顿流体为水或甘油。
4.根据权利要求1所述的一种流固两相输运中颗粒定位的方法，其特征在于：所述的压力源采用压力泵或注射泵。
5.根据权利要求1所述的一种流固两相输运中颗粒定位的方法，其特征在于：所述的调节压力源的驱动压力和颗粒直径，具体是：调节压力源的驱动压力使得雷诺数为14、82～120，调节颗粒直径使得阻塞率为0.125～0.4，使得固体颗粒间的间距能达到均匀稳定，形成固定的颗粒链。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910828574.0A CN110465339B (zh) | 2019-09-03 | 2019-09-03 | 一种流固两相输运中颗粒定位的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910828574.0A CN110465339B (zh) | 2019-09-03 | 2019-09-03 | 一种流固两相输运中颗粒定位的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110465339A CN110465339A (zh) | 2019-11-19 |
| CN110465339B true CN110465339B (zh) | 2021-02-09 |
Family
ID=68514732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910828574.0A Active CN110465339B (zh) | 2019-09-03 | 2019-09-03 | 一种流固两相输运中颗粒定位的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110465339B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113769797B (zh) * | 2021-09-02 | 2023-03-14 | 浙江理工大学 | 一种流固两相输运中的微尺度颗粒直径测定的方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6915245B1 (en) * | 2000-09-14 | 2005-07-05 | General Atomics | Method of simulating a fluid by advection of a time-weighted equilibrium distribution function |
| CN107110765A (zh) * | 2014-11-03 | 2017-08-29 | 通用医疗公司 | 在微流体装置中浓缩颗粒 |
| CN107213929A (zh) * | 2017-06-06 | 2017-09-29 | 国家纳米科学中心 | 一种基于界面效应的微纳颗粒分离系统 |
| CN107349981A (zh) * | 2016-05-10 | 2017-11-17 | 李榕生 | 一种采取新的流体驱动方式的微流控芯片装置 |
| CN108037170A (zh) * | 2017-12-13 | 2018-05-15 | 西安理工大学 | 基于惯性迁移与介电电泳的复合式细胞分离系统及方法 |
| CN109800469A (zh) * | 2018-12-25 | 2019-05-24 | 上海交通大学 | 基于ib-lb方法对多颗粒链颗粒平衡间距进行预测的模拟仿真方法 |
Family Cites Families (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0658B2 (ja) * | 1984-11-07 | 1994-01-05 | 株式会社日立製作所 | 細胞取扱い装置 |
| JP4614947B2 (ja) * | 2003-03-28 | 2011-01-19 | イングラン・リミテッド・ライアビリティ・カンパニー | 粒子を仕分けるためおよび性別で仕分けられた動物の精子を提供するための装置、および方法 |
| US20070237706A1 (en) * | 2006-04-10 | 2007-10-11 | International Business Machines Corporation | Embedded nanoparticle films and method for their formation in selective areas on a surface |
| US10052571B2 (en) * | 2007-11-07 | 2018-08-21 | Palo Alto Research Center Incorporated | Fluidic device and method for separation of neutrally buoyant particles |
| EP2664916B1 (en) * | 2007-04-02 | 2017-02-08 | Acoustic Cytometry Systems, Inc. | Method for manipulating a fluid medium within a flow cell using acoustic focusing |
| US20110163013A1 (en) * | 2008-05-30 | 2011-07-07 | Eppendorf Ag | Apparatus and Method for Moving Particles in a Fluid |
| US8356714B2 (en) * | 2009-06-02 | 2013-01-22 | Georgia Tech Research Corporation | Microfluidic device for separation of particles |
| US9353342B2 (en) * | 2010-01-21 | 2016-05-31 | Emd Millipore Corporation | Cell culture and gradient migration assay methods and devices |
| JP2013152171A (ja) * | 2012-01-26 | 2013-08-08 | Tohoku Univ | 血液中の粒子分離装置 |
| CN102746986B (zh) * | 2012-07-13 | 2014-07-16 | 中国科学院大连化学物理研究所 | 一种基于微流控芯片的肿瘤细胞迁移动力学监测方法 |
| DE102013210952B4 (de) * | 2013-06-12 | 2020-03-19 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren zur Bestimmung von ungelösten Teilchen in einem Fluid |
| CN203569079U (zh) * | 2013-10-31 | 2014-04-30 | 吉林大学 | 一种用于细胞筛选的密封腔分离器 |
| US20170262559A1 (en) * | 2016-03-11 | 2017-09-14 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and Systems for Simulating Nanoparticle Flux |
| CN108073743A (zh) * | 2016-11-14 | 2018-05-25 | 中国科学院力学研究所 | 基于非牛顿效应聚焦分离亚微米纳米颗粒的系统及方法 |
| CN109967150B (zh) * | 2019-04-24 | 2021-07-23 | 河海大学常州校区 | 一种用于操控微纳米颗粒的惯性微流控芯片 |
-
2019
- 2019-09-03 CN CN201910828574.0A patent/CN110465339B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6915245B1 (en) * | 2000-09-14 | 2005-07-05 | General Atomics | Method of simulating a fluid by advection of a time-weighted equilibrium distribution function |
| CN107110765A (zh) * | 2014-11-03 | 2017-08-29 | 通用医疗公司 | 在微流体装置中浓缩颗粒 |
| CN107349981A (zh) * | 2016-05-10 | 2017-11-17 | 李榕生 | 一种采取新的流体驱动方式的微流控芯片装置 |
| CN107213929A (zh) * | 2017-06-06 | 2017-09-29 | 国家纳米科学中心 | 一种基于界面效应的微纳颗粒分离系统 |
| CN108037170A (zh) * | 2017-12-13 | 2018-05-15 | 西安理工大学 | 基于惯性迁移与介电电泳的复合式细胞分离系统及方法 |
| CN109800469A (zh) * | 2018-12-25 | 2019-05-24 | 上海交通大学 | 基于ib-lb方法对多颗粒链颗粒平衡间距进行预测的模拟仿真方法 |
Non-Patent Citations (4)
| Title |
|---|
| Inertial migration of circular particles in Poiseuille flow of a power-law fluid;Xiao Hu et al.;《Physics of Fluids》;20190724;第31卷(第7期);全文 * |
| Inertial migration of neutrally buoyant spheres in a pressure-driven ow through square channels;Kazuma Miura et al.;《Journal of Fluid Mechanics》;20140515;第749卷;全文 * |
| 微通道内单颗粒周围微观流场的直接定量观测;张润林、潘振海、吴慧英;《工程热物理学报》;20180430;第39卷(第4期);全文 * |
| 管流中颗粒"惯性聚集" 现象的研究进展及其在微流动中的应用;王企鲲、孙仁;《力学进展》;20121125;第42卷(第6期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110465339A (zh) | 2019-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Taitel et al. | A model for slug frequency during gas-liquid flow in horizontal and near horizontal pipes | |
| Armaly et al. | Experimental and theoretical investigation of backward-facing step flow | |
| Yan et al. | Numerical simulation of junction point pressure during droplet formation in a microfluidic T-junction | |
| CN110465339B (zh) | 一种流固两相输运中颗粒定位的方法 | |
| Kumara et al. | Particle image velocimetry for characterizing the flow structure of oil–water flow in horizontal and slightly inclined pipes | |
| Gleghorn et al. | Transport and collision dynamics in periodic asymmetric obstacle arrays: Rational design of microfluidic rare-cell immunocapture devices | |
| CN110836766B (zh) | 一种双流体喷嘴的检测装置及检测方法 | |
| Wang et al. | Investigating the flow characteristics of air-lift pumps operating in gas–liquid two-phase flow | |
| CN113769797B (zh) | 一种流固两相输运中的微尺度颗粒直径测定的方法 | |
| Franklin et al. | Morphology and displacement of dunes in a closed-conduit flow | |
| Hasan et al. | Countercurrent bubble and slug flows in a vertical system | |
| Wei et al. | IB-LBM simulation on blood cell sorting with a micro-fence structure | |
| Pantzali et al. | Three-component particle velocity measurements in the bottom section of a riser | |
| Petrova et al. | Estimations of gas flow maldistribution in packed‐bed columns | |
| Hossain et al. | CFD investigation of particle deposition around bends in a turbulent flow | |
| Alekseenko et al. | Investigation of secondary waves dynamics in annular gas–liquid flow | |
| Diaper et al. | Crossflow Pressure Drop and Flow Distribution within a Tube Bundle Using Computational Fluid Dynamics | |
| Lopez et al. | Experimental investigation of falling liquid film in vertical downward two-phase pipe flow | |
| CN110261062A (zh) | 一种减阻剂筛选与评价装置及使用方法 | |
| Mudawwar | Interfacial instabilities of air-driven liquid films | |
| CN115608111A (zh) | 一种基于微通道在线检测气泡实时浓度的方法 | |
| Charru et al. | Sand transport over a barchan dune | |
| Zhang et al. | Void fraction effect on bubble deformation of turbulent bubbly flow in a large square duct | |
| Kolf et al. | The vortex chamber as an automatic flow-control device | |
| Cageao et al. | Experiments on the dynamics of subaerial two-phase debris flows |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |