[发明专利]一种基于动态荧光粉浓度确定均匀扁平微通道平均流速和剪切力的方法

摘要

本发明提供了一种基于动态荧光粉浓度确定均匀扁平微通道平均流速和剪切力的方法，属于细胞生物力学实验装置技术领域。所用装置包括动态荧光粉溶液产生装置、均匀扁平微流控芯片、荧光显微镜和废液回收容器；本发明通过两组可编程控制的注射泵和注射器产生动态荧光粉溶液，动态荧光粉溶液在均匀扁平微通道中的传输过程满足Taylor‑Aris弥散方程，荧光粉溶液在微通道中的传输过程利用荧光显微镜进行实时记录并得到一系列的荧光图像。对荧光图像分析得到微通道一段距离内荧光粉溶液浓度随时间的变化，通过对荧光粉溶液在均匀扁平微通道中的Taylor‑Aris弥散进行逆向求解，计算出均匀扁平微流通道内的平均流速和底部剪切力大小。

主权项

一种基于动态荧光粉浓度确定均匀扁平微通道平均流速和剪切力的方法，其特征在于如下步骤：该方法采用的装置包括动态的荧光粉溶液产生装置、均匀扁平的微流控芯片、荧光显微镜和废液回收容器；其中动态的荧光粉溶液产生装置包括可编程控制的泵、可编程控制的注射器和三通接口，可编程控制的泵和可编程控制的注射器通过三通接口导入至均匀扁平的微流控芯片，均匀扁平的微流控芯片上的废液直接通入废液回收容器；荧光显微镜进行实时记录并得到一系列的荧光图像；从均匀扁平的微流控芯片入口处加载浓度随时间变化的荧光粉溶液，保证宽度x方向的荧光粉浓度相同；动态荧光粉溶液在微通道中传输受到流动的影响，满足对流‑扩散公式∂φ∂t+uz∂φ∂z=D(∂2φ∂y2+∂2φ∂z2)---(1)]]>其中，t为时间，x,y,z分别是宽度、高度、长度方向的坐标，φ＝φ(y,z,t)是荧光粉溶液浓度，uz＝uz(y,t)是荧光粉溶液流体速度，D是荧光粉溶液扩散系数；由于均匀扁平微通道几何尺寸小，且微通道中的荧光粉溶液流体运动为小雷诺数流动，Womersley数小，满足准定常假设条件，微通道中荧光粉溶液的流速和底部剪切力分别满足uz(y,t)==3u‾z(t)2[1-(2yH)2]---(2)]]>τw=η∂uz∂y|y=-H/2=6ηu‾z(t)H---(3)]]>其中，为高度方向的平均流速，H为待检测的均匀扁平微通道的高度，η表示粘度系数；由于微通道高度小，荧光粉溶液在高度方向上形成均匀浓度；高度方向上的平均浓度定义为φ‾=1H∫-H/2H/2φ(y,z,t)dy---(4)]]>满足Taylor‑Aris弥散公式∂φ‾∂t+u‾z∂φ‾∂z=Deff∂2φ‾∂z2---(5)]]>Deff称为有效扩散系数，满足Deff=D[1+1210(u‾zHD)2]---(6)]]>以空间步长Δz将长度沿z方向均匀离散，网格点为zi，其中i＝1,2,.....I+1，同时用时间步长Δt将时间t均匀离散，时间网格点为tk，其中k＝1,2,......K+1，则公式(5)用有限差分近似为φ‾ik-φ‾ik-1Δt+u‾z(tk)φ‾i+1k-φ‾i-1k2Δz=(D+u‾z(tk)2H2210D)φ‾i+1k-2φ‾ik+φ‾i-1k(Δz)2---(7)]]>其中，I表示空间离散网格数，K表示时间离散网格数，分别表示tk时刻zi‑1、zi、zi+1位置处的荧光粉浓度，表示tk‑1时刻zi位置处的荧光粉浓度；通过荧光显微镜测得各个时刻微通道中的荧光粉浓度分布，得到时间间隔为Δt的一系列荧光图像；把荧光图像的每个像素点看作荧光粉浓度的采样点，使相邻像素间的距离为上述的Δz；对公式(7)进一步整理为关于的公式如下：aiu‾z(tk)2+biu‾z(tk)+ci=0---(8)]]>其中，ai=φ‾i+1k-2φ‾ik+φ‾i-1k(Δz)2H2210Dbi=-φ‾i+1k-φ‾i-1k2Δzci=φ‾i+1k-2φ‾ik+φ‾i-1k(Δz)2D-φ‾ik-φ‾ik-1Δt---(9)]]>在公式(8)和(9)中，D和H为已知常数，利用tk时刻z＝zi部位及前后Δz位置的浓度以及tk‑1时刻z＝zi部位的浓度代入公式(9)计算ai，bi，ci的值；公式(8)是一个关于未知变量的一元二次公式，该公式的解即为tk时刻的平均流速根据上述数值方法，相邻时刻的荧光图像上n个相邻像素点位置z＝zi，i＝1,2,.....n的荧光粉溶液浓度一共组成n‑2个公式(8)形式的一元二次公式；由于相邻像素间荧光粉溶液浓度差异很小，为了减少误差，将n个相邻像素点位置z＝zi，i＝1,2,.....n的荧光粉溶液浓度构成的n‑2个公式叠加平均，得到如下公式au‾z(tk)2+bu‾z(tk)+c=0---(10)]]>式中，进一步得出，根据公式(10)和(11)，通过tk时刻和tk‑1时刻n个相邻像素点位置荧光粉溶液浓度求出系数a、b和c，进而对一元二次公式(10)求解即得tk时刻的平均流速根据平均流速再根据公式(3)计算微通道底部的剪切力大小。