利用示踪粒子的关联扩散判断液体表面漂移运动方法技术

本发明专利技术属于软物质及生物物理技术领域，具体为利用示踪粒子的扩散关联计算液体表面漂移运动方法。本发明专利技术根据相距一定距离s的颗粒对间的x，y轴方向上关联扩散位移来计算样品中的表面流速。主要是选定较大间距s的颗粒对，计算其关联位移，利用

【技术实现步骤摘要】

本专利技术属于主要在软物质、生物物理
，具体及涉及各类液面上存在的表 面流的测量计算方法。
技术介绍
各类液面体系在科研工业领域中普遍存在。表面流的性质往往影响了整个 体系的特征表现。通过示踪粒子的位置观察可以计算表面流的特征参数。一般都是以胶体 粒子作为示踪粒子。在胶体体系中，在没有外场力的情况下，胶体颗粒在白噪声的驱动下做 布朗运动(如图1所示的二维胶体系统所对应的图片），但在很多体系中由于表面流的存 在，示踪颗粒存在布朗运动和漂移运动。因此，要计算出样品体系中存在的表面流，去除体 系中表面流的漂移运动。传统方法是：（1)计算所有示踪颗粒的位移的平均值 Δ?_ :(Δ?_ = ?:(τ) 为t = 0时刻颗粒的初始位置，:?(f)为t = T时颗粒的位置。示踪粒子在 做布朗运动和漂移运动的情况下，在一定扩散时间τ内，其位移Δ? = + (Δ&:布朗 运动引起的位移，Δ?β:漂移运动引起的位移），理想状态下，噪声的平均值为零，即$ = 0， 胶体颗粒的平均位移及=&。但这种方法需要对胶体颗粒位移做大量的统计，才能使胶 体颗粒布朗运动位移@ =0〇这种方法计算精度误差与统计样本数目Ν平方根的倒数 (1/VI)成正比，若样品中的统计胶体颗粒数较少，误差较大，则^ 土 0，涵=. 因此，这种方法计算表面流速度不是很准确；（2)根据单个示踪粒子的平均平方位移 Δ尹2( Ti_).(MSD:mean square displacement)，计算胶体颗粒的漂移速率。j根据统计理论及 实验数据，当实验样品中颗粒密度(area fraction)较小且经过足够的扩散时间τ时，若样 品中不存在漂移运动时，颗粒的MSD即兩（τ )随τ线性变化:其中η:空间维度，D:颗粒扩散系数，τ :扩散时间。当样品系统中存在漂移运动(漂 移流恒稳)时， 漂移速度。根据τ )拟合的二次项系数，可计算出漂移运动的速度。但是，当 实验样品中颗粒密度较大时，此方法不适用。因为此时颗粒由于近邻粒子的囚笼效应，单个 颗粒做布朗运动的平均平方位移与扩散时间τ不再呈线性关系，则从(2)式中 τ )的拟合参数中，不能从二次项的拟合参数计算出漂移运动的速度。
技术实现思路
本专利技术的目的在于提供一种普适性更好而且计算精度更高的利用示踪粒子的扩 散关联计算液体表面漂移运动方法。 本专利技术提供的利用示踪粒子的扩散关联计算液体表面漂移运动方法，是通过相距 一定距离s的颗粒间的关联扩散位移随扩散时间T的关系，计算颗粒的漂移运动速度。 在示踪粒子系统中（如图1所示），颗粒间的关联扩散系数?ΦΦ、〇θθ(Φ :平行于颗粒 质心连线方向，Θ:垂直于颗粒质心连线方向）随颗粒间距s的关系为： rl、r2表征了 ΟφφΑθθ随颗粒间距s衰减快慢程度。故随着s的增大，0ΦΦ(〇θθ)衰减， 当s较大时，?ΦΦ(〇θθ)衰减到很弱的程度。?ΦΦ(〇θθ)表征了胶体中颗粒间流体力学相互作用 。为了计算样品中示踪粒子的漂移运动位移，本专利技术计算图片的x、y轴方向上的颗粒 若示踪粒子系统中存在漂移运动，则颗粒的位移（Δ?)由无规运动产生的位移 (Δ&』与漂移运动产生的位移(Δ&)两部分组成，即Δ? = Δ&+Δ^，当两颗粒间距s-定时， 在扩散时间τ内，X、y轴方向上两颗粒间的关联扩散位移、Δ/?；；为： 其中，x、y是图片的x、y轴，颗粒i、j，颗粒i布朗运动位移在X轴方向的位移， :颗粒i布朗运动位移在y轴方向的位移，颗粒i漂移运动位移在X轴方向的位移， Δγ&:颗粒i漂移运动位移在y轴方向的位移。对上式两边取平均有： f(T)为流体的流体力学引起的颗粒间的关联位移，F(t)为流体漂移引起的颗粒间 的关联位移。因为对于给定的恒定表面流，颗粒i、j的漂移运动速度相同，则有F(〇=C 2 ·(表面流速率），f (τ) = B · τ (B是颗粒间i、j相距s时两颗粒间的流 体力学关联强度，C是漂移运动在x、y轴方向上对应的速率），F(t)不随颗粒间距而变 化。Μ(τ)= ￡(τ)+Ρ(τ)、<〇)= f(T)+:F(T)，随着颗粒间距s的增大，f⑴ 减弱，F(t)不随s的变化而变化。当颗粒间距s较大时（s>>d，d:颗粒直径），f(T)减到很弱 的程度，Ζ^(τ), Ζ^(τ)主要受F(t)影响(F⑴〉〉f(T))，如图2、4所示，此时可用多 项式y = B · τ+C2 · τ2拟合计算的数据τ )、<( τ )，根据二次项系数可计算颗粒漂 移速度在x、y方向上的分量vx、vy(如图6、8)。 根据拟合的漂移运动速率，确定表面流速度：因为拟合出的漂移运动速率没有标 识速度方向，因此若要确定表面流的方向，就需要计算出所有的运动方向的可能性，然后从 中选择出需要的一组。x，y两个方向上所有的漂移速率的组合方式有四种： (Vx，Vy)、（_Vx，Vy)、（Vx，_Vy)、（ _Vx，_Vy) 从这四组漂移速度分别从颗粒it时刻原始位置(li) , 减去漂移运动Vx · τ、 Vy * T ：：分别在这四组新的颗粒位置(rUt),ri.⑴）中，选定间距S的颗粒对，分别计算所对 应的关联扩散位移，&S：T.)、<〇)，从中选出Ζ??τ)、Ζ^(τ)与τ线性最好的一 组，其所对应的 Vx、Vy即是漂移运动的速率分量，Vx、Vy前面的符号即是漂移运动速度在x、y 轴上的方向。本专利技术与传统方法相比，其优势在与：( - )传统方法（1)需要对颗粒做大量的统计，若样品中统计量不足，则计算出的漂 移运动位移不准确；(二）传统方法（2 )根据拟合颗粒的平均平方位移（MSD : mean s quar e displacement)，计算胶体颗粒的漂移速率方法。在一定的扩散时间τ内，该方法中颗粒由布 朗运动引起的MSD(2 · η · D · τ)强于两颗粒间的关联扩散位移^( τ )、_<( τ 吏得 由漂移运动速度奋引起的关联位移(?!·τ:)2项相对较弱，相当于在一个较强的背景中提取一 个较弱的信号，从而计算出的漂移运动速率不准确。而本专利技术的方法中，利用颗粒间关联扩 散位移Δ/??( τ ) = f( τ ) +丨''(〇、=: f( τ ) +F( τ ),在相同的扩散时间τ内，f(T) 远小于布朗运动引起颗粒的MSD(2 · η · D · τ)项，相当于在一个很弱的背景中提取一个很 强信号，所以计算出的漂移运动速率较准确；(三）当样品中颗粒密度较大时，利用拟合颗粒的平均平方位移方法中，颗粒的 MSDfMSD(〇与扩散时间τ非线性，τ )=/ΜΛ3_+(^· ，因此不能从二次项中拟合出漂 移速度；(四）当样品中的颗粒数为N时，根据颗粒位移平均值方法和计算颗粒MSD方法中的 统计量都为N，而利用本专利技术计算颗粒间的关联扩散位移方法的统计量为N · (N-1)，因此， 本专利技术方法能够更有效地利用实验数据，增大计算结果的统计量，使计算出的漂移速度更 准确。【附图说明】 图1将直径为2.Ομπι silica颗粒置于油水界面附近，图为采集到的二维运动图片， 其面积分数m = 2 %，实本文档来自技高网...

【技术保护点】
一种利用示踪粒子的关联扩散位移计算液体表面漂移运动方法，其特征在于：利用相距一定间距s的两颗粒间x、y轴方向上的关联扩散位移与扩散时间τ的关系ΔRxxij‾(τ)=]]>f(τ)+F(τ),]]>ΔRyyij‾(τ)=f(τ)+F(τ),]]>判断样品中颗粒是否存在漂移运动，若F(τ)＝C2·τ2，则样品中颗粒存在漂移运动；若F(τ)＝0，则样品中颗粒不存在漂移运动；其中，为x、y轴方向上的颗粒i、j间的关联位移，F(τ)为液体漂移动运动引起的颗粒i、j间的关联位移，f(τ)为流体的流体力学相互作用引起的颗粒i、j间的关联位移；τ：扩散时间，为表面流速度关联。

【技术特征摘要】

【专利技术属性】
技术研发人员：陈唯，李娜，
申请(专利权)人：复旦大学，
类型：发明
国别省市：上海;31