本文介绍了一种实验上可实现的模拟方法,用于研究稠密流体中的分子流体动力学过程。该技术利用振动下高恢复系数颗粒堆的粒子图像测速法,能够直接在宏观尺度上观测到强相互作用、高密度气体和液体中已知或预测存在的动力学过程。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本文介绍了一种实验上可实现的模拟方法,用于研究稠密流体中的分子流体动力学过程。该技术利用振动下高恢复系数颗粒堆的粒子图像测速法,能够直接在宏观尺度上观测到强相互作用、高密度气体和液体中已知或预测存在的动力学过程。
本文描述了一种用于研究稠密气体和液体中分子尺度流体动力学过程的类比宏观方法。该技术应用标准的流体动力学诊断手段——粒子图像测速法(PIV),用于测量:i)在短时间颗粒碰撞时间尺度上单个粒子(颗粒)的速度;ii)在短碰撞时间和长连续流动时间尺度上粒子系统的速度;iii)在稠密分子流体中已知存在的集体流体动力学模式;以及 iv)短时间和长时间尺度的速度自相关函数,后者对于理解强相互作用稠密流体系统中粒子尺度的动力学至关重要。该基本系统由成像系统、光源、振动传感器、具有已知介质的振动系统,以及PIV和分析软件组成。文中重点介绍了采用该类比技术研究分子尺度流体动力学过程所需的实验测量步骤,以及相关的理论工具概要。该方法为传统上用于分子流体动力学研究的光子散射和中子束散射技术提供了一种相对简便的替代方案。
分子流体动力学研究流体中单个分子及分子集合的动态行为与统计力学。在已开发的众多用于研究分子流体动力学系统的实验技术中1,2,光散射1,2,3、分子动力学模拟4,5,6,7以及在较小程度上的非弹性中子散射8应用最为广泛。然而,后两种技术存在显著局限性。例如,分子动力学(MD)模拟:i)受限于较小的空间和时间尺度
,仅能包含相对较少的分子数量
;ii)需使用近似的粒子间势能函数;iii)通常引入周期性边界条件,而该条件在非平衡体相流动条件下不成立;iv)目前尚无法回答一个基本科学问题:涉及单个分子或分子集合的分子尺度动力学如何受到体相非平衡流体流动的影响,并反过来与之耦合。中子散射的主要局限性则源于可利用的中子束源数量极为有限。
为了给本文介绍的模拟实验技术提供背景,我们重点介绍应用于简单稠密气体和液态流体的光散射技术。在典型的光散射实验中,一束偏振的激光被引导至包含静止流体样品的小探测体积中。然后,在相对于入射光束的某一固定角度处,检测样品内分子散射的光。根据所关注的分子动力学区域,对散射光信号的检测与分析采用光滤波或光混频检测方法。如 Berne 和 Pecora1 所述,滤波技术用于探测时间尺度短于
秒的流体状态分子动力学,其在散射后引入干涉仪或衍射光栅,从而实现对散射光谱密度的扫描。而光学混频技术则用于慢时间尺度的动力学过程,即
秒的情况,该方法在散射后采用自相关器或频谱分析仪,从测得的散射光强度中提取散射信号的频谱信息。
通常情况下,激光探针(至少工作在可见光谱范围内的探针)的波长远大于液态分子间的特征间距。在此条件下,探测光束会激发五种集体性的、慢时间尺度、长波长的流体动力学模式2,9,10(相对于特征碰撞频率而言较慢):两个黏性衰减的、反向传播的声波模式,两个未耦合的纯扩散性涡旋模式,以及一个单一的扩散性热(熵)模式。声波模式沿入射光束方向(纵向)被激发,而涡旋模式则在横向方向被激发。
仅考虑实验散射技术时,有两个基本问题超出了光散射和中子散射测量的能力范围,这两个问题位于分子液体体系的平衡态与非平衡态统计力学的核心。
1) 严谨的论证9,11 表明在经典牛顿动力学或量子动力学作用下的单个液态分子的随机碰撞及次碰撞时间尺度的动力学行为,可以重新表述为广义朗之万方程(GLE)的形式。GLE 进而构成了研究稠密气体和液体中分子非平衡态统计力学的核心理论工具。然而,由于单个(非大分子)分子的动力学行为无法通过任何散射技术分辨,目前除了分子动力学(MD)模拟之外,尚无直接方法可用于检验 GLE 的有效性。
2) 宏观连续流体动力学及微观分子流体动力学的核心均基于一个基本假设:在长度和时间尺度上远大于分子直径和碰撞时间,但远小于连续介质长度和时间尺度的条件下,局部热力学平衡(LTE)成立。在诸如纳维-斯托克斯(NS)方程之类的连续流体流动与传热模型中,必须采用LTE假设9 为了将本质上处于非平衡态的连续尺度流动和能量输运特征——如黏性剪切应力和热传导——与严格处于平衡态的热力学性质(如温度和内能)相耦合。同样,尽管微观尺度上的动量和能量输运本质上是非平衡过程,反映了耦合的微观质量、动量和能量流的出现,但这些微观过程的模型仍假设这些流是相对于局部热力学平衡(LTE)的微小扰动9据我们所知,目前尚未有针对局部热力学平衡(LTE)假设的直接实验验证。特别是,在稠密、流动且处于非平衡态的流体中,似乎尚未开展任何分子流体动力学散射实验。
本文介绍了一种模拟实验技术,该技术利用标准的粒子成像测速法(Particle Imaging Velocimetry, PIV)测量振动颗粒堆的宏观、单粒子及集体粒子动力学行为,从而间接预测、解释并揭示稠密气体和液体中单分子与多分子的流体动力学特性。本技术所依赖的物理与理论基础已在我们课题组近期发表的一篇论文中阐明12。在实验方面,宏观系统必须满足以下条件:(i)持续趋向局部宏观尺度的统计力学平衡态;(ii)存在微小且线性的偏离平衡态的行为,以模拟分子流体动力学系统中观察到的(弱)非平衡涨落。在理论方面,需满足:(i)描述稠密相互作用的N粒子系统在平衡态及弱非平衡态下的经典微观尺度统计力学模型,必须被重新表述为宏观尺度形式;(ii)所得的宏观尺度模型必须能够可靠地预测从短时间尺度(粒子碰撞时间尺度)到长时间尺度(连续流动时间尺度)范围内的单粒子与多粒子动力学行为。
本文介绍了一种新方法的详细实验方案以及由此技术获得的代表性结果。与分子动力学模拟以及光散射和中子散射方法相比,该新技术首次实现了对流动状态下强非平衡态稠密气体和液体中分子流体动力学过程的深入研究。
访问受限。请登录或开始试用以查看此内容。
1. 振动系统的准备
,其中
为介质的材料装填体积(本实验设置中,共添加 13,750 mL 水)。所添加水的体积取决于所用介质的类型。
为介质的装填密度,
为介质的质量(本实验中,计算得到的密度为 1649
)。2. 高速成像
注意:对于通过成像流动颗粒堆表面部分区域获得的颗粒速度场测量,成像区域
对应于下文步骤 2.2.4 中确定的视场(FOV)。通过在
内选择一个微小且固定的子区域
,可获得随时间变化的单个颗粒速度(在堆体表面)的测量值;其中,如下文所述,
的量级约为单个颗粒投影面积的大小。
,必须至少比施加的振动频率大一个数量级,
(此处,
Hz3. 数据采集
注意:如果实验方案曾被暂停,则需要重新启动相机。请执行步骤 3.1。如果实验方案未被暂停,则跳至步骤 3.1.2。

必须足够长,以使颗粒流动达到统计上的稳态条件;ii)
也不应过长,以免产生大量冗余数据。稳态条件出现的时间尺度需通过试错法确定,可采用多种严格程度不同的方法,例如:i) 确保在固定点或多个固定点处的颗粒平均速度达到名义上恒定的幅值;或 ii) 除了均值达到稳态外,相应的方差也达到名义上恒定的幅值。本实验中,数据采集持续了10.12秒,共获取5060帧图像。颗粒流在约1秒后进入稳定状态。
进行图像记录。相机会在指定的
时间内持续采集图像,并将这些图像保存至其内部存储器中。图2a 是所采集的5060幅图像中的一幅示例。4. 使用PIV处理视频数据
,最终通过将
除以帧间时间增量
,得到区域平均速度
,其中
指第
个 interrogation area。在当前实验中,每个 interrogation area 由 n x n = 32 x 32 像素组成;每个 210 mm x 160 mm 视场(FOV)被划分为 47 x 35 个 interrogation area,对应 1504 x 1128 像素。5. 处理振动数据
注意:如果使用不同的计算机系统或分析软件,步骤 5 可与步骤 4 同时进行。
访问受限。请登录或开始试用以查看此内容。
在展示代表性结果时,我们参考了 连续时间尺度 过程与在时间尺度上观察和预测的过程
相对于特征性颗粒碰撞时间尺度较长的,
和 粒子-时间尺...
访问受限。请登录或开始试用以查看此内容。
为了将振动颗粒柱作为研究分子流体动力学过程的宏观类比,实验人员一方面必须学习并使用四种基本测量方法,另一方面需掌握平衡态与非平衡态统计力学的一些基本要素。首先关注实验测量,这些方法包括:i)通过测量单颗粒速度自相关函数来研究单个颗粒的动力学行为;ii)测量时间平均/长时间尺度下的表面颗粒速度场;iii)测量颗粒介质的有效黏度;iv)测量介质容器在空载和装载介质状态下的振动频谱。
单颗粒速度自相关函数的测量
通过测量单个粒子速度的自相关函数,研究微尺度系统中分子或本方法中受振动颗粒的个体随机动力学行为,
2. 对于小型的, 例如双原子和三原子分子
访问受限。请登录或开始试用以查看此内容。
作者无任何利益冲突需要披露。
本工作由美国海军研究办公室(ONR N00014-15-1-0020)[Tkacik 和 Keanini] 资助,并在北卡罗来纳大学夏洛特分校赛车运动研究实验室完成。抛光介质由 Rosler 公司捐赠。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 振动抛光碗 | Raytech | AV-75 | |
| 流量计 | Peristaltic Pumps | 913 Mity Flex | |
| 电子天平 | Pelouze | 4040 | |
| 三轴加速度计 | PCB Piezotronics | PCB 356B11 | 带传感器信号调理器的加速度计 |
| 数据采集计算机 | IBM | Thinkpad | 与高速相机联用 |
| 高速相机 | Redlake | Motionxtra HG-XR | |
| 变焦镜头 | Tamron | 型号 A18 | 18-250mm F/3.5-6.3 |
| 高强度光源 | ARRI | EB 400/575 D | |
| 数据处理计算机 | Dell | Dell Precision Tower 7910 | |
| PIV 软件 | Dantec Dynamics | Dynamic Studio 2013 | 版本 3.41.38 |
| 数据采集硬件 | National Instruments | SCXI | SCXI-1000 机箱,配备 SCXI 1100 卡和 SCXI 1303 适配器 |
| 数据采集软件 | National Instruments | LabVIEW 2012 | |
| 数据处理软件 | MATHWORKS | MATLAB | |
| 抛光介质 | Rosler | RSG 10/10S | 使用多种介质类型(混合型、球形、三角形) |
| 抛光液 | Rosler | FC KFL (3%) | 含3%肥皂溶液的水溶液 |
| 刻度尺 | Swiss Precision Instruments | 13-911-3 | |
| 量筒 | Global Scientific | 601082 |
访问受限。请登录或开始试用以查看此内容。