本文描述了一种能够获取光学不透明流体中或通过光学不透明几何结构的二维速度场的回声粒子图像测速（EPIV）系统，并报道了在管道流动中的验证性测量结果。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本文内容
摘要
摘要
流体流动中质量、动量和能量的输运最终由流体速度场的时空分布决定。1 因此，要理解、预测和控制流体流动，前提在于能够以足够的空间和时间分辨率测量速度场。2 对于在光学不透明流体中或通过光学不透明几何结构进行的速度测量，回波粒子图像测速技术（EPIV）是一种有吸引力的诊断方法，可用于生成速度的“瞬时”二维场。3,4,5,6 本文介绍了通过将商用医用超声设备7与运行商用粒子图像测速（PIV）软件的计算机8集成而构建的EPIV系统的操作规程，并报告了在哈根-泊肃叶（i.e.，层流管流）流动中的验证测量结果。
在进行EPIV测量时，将相控阵探头连接到医用超声成像仪，通过在不同时间脉冲激励压电探头元件来生成二维超声图像。每个探头元件向流体中发射超声波脉冲，流体中的示踪粒子（自然存在或人为加入）将超声回波反射回探头并被记录下来。反射超声波的幅度及其相对于发射信号的时间延迟被用于生成所谓的B型（亮度模式）二维超声图像。具体而言，时间延迟用于确定流体中散射体的位置，而幅度则用于设定散射体的亮度强度。获取单幅B型图像所需的时间t由对相控阵探头所有元件完成一次脉冲激励所需时间决定。在采集多幅B型图像时，系统的帧率以每秒帧数表示，即(fps) = 1/δt。（有关超声成像的综述，参见文献9。）
对于典型的EPIV实验，帧率通常介于20–60 fps之间，具体取决于流动条件，并需采集100–1000帧显示示踪粒子在流场中空间分布的B模式图像。采集完成后，B模式超声图像通过以太网连接传输至运行PIV商业软件的计算机。利用PIV软件，可获得示踪粒子的位移场， D(x,y)[像素]（其中 x 和 y 分别表示超声图像中的水平和垂直空间位置）是通过对连续的超声 B 模式图像应用互相关算法而获得的。10 速度场， u(x,y)[m/s]，根据图像对之间的时间步长和位移场确定， ΔT[s]，以及图像放大倍数， M[米/像素] 即, u(x,y) = MD(x,y)/ΔT图像之间的时间步长 ΔT = 1/帧每秒 + D(x,y)/B，其中 B[pixels/s] 是超声探头扫过图像宽度所需的时间。在本研究中，M = 77[μm/pixel]。 fps = 49.5[1/s]，且 B = 25,047[像素/秒]。获取速度场后，可对其进行分析以计算感兴趣的流动参数。
方案
1. 创建可测量的流体流动
- 将在甘油水溶液（50% 甘油 - 50% 水）的管道流动中演示 EPIV 验证测量。实验装置的示意图见图 1。
- 将名义直径为 10 μm 的空心玻璃微球以约 17 重量 ppm 的浓度添加到流体中。这些空心玻璃微球作为超声造影剂，其尺寸和密度的选择使其能够被动跟随流体流动。10
- 向泵提供固定电压，以引入已知的流速。选择该流速使得 U << ΔX/δt，其中 U 为管道中的平均流速，ΔX 为 EPIV 测量体积的线性长度，δt 为图像之间的时间步长，即，与超声系统的 fps 相比，流动必须足够“慢”。3
2. 校准超声设备
- 将超声探头安装到管道外壁。在超声探头表面涂抹水性耦合凝胶，以尽量减少超声波束在探头表面与管道壁之间传播时的信号损失。
- 开启超声成像设备电源。系统完成加载后，将自动开始实时超声图像采集。
- 使用超声设备控制面板上的Depth Control（深度调节）旋钮设置图像深度。
- 使用超声设备控制面板上的2D Gain（二维增益）旋钮设置整体图像增益。
- 调节Time Gain Compensation (TGC)（时间增益补偿）滑块，以减弱来自管道壁的散射信号，并补偿超声信号随深度增加而产生的衰减。
- 使用可分配功能的控制旋钮调节图像宽度、聚焦位置、探头工作频率和帧率。这四个旋钮位于控制面板左上方，其功能随系统运行模式的不同而变化。在当前使用的二维（2D）模式下，从左至右依次对应宽度、聚焦、频率和帧率。需注意，由于超声成像的基本原理9，这四个参数本质上相互关联。因此，在特定的超声扫描条件下（即，EPIV实验中），空间分辨率与时间分辨率之间存在权衡关系。
- 参见图2，其中展示了在管道流动中加入10 μm空心玻璃微球后的典型超声图像。注意，由于横向分辨率有限，玻璃微球在图像中沿横向被展宽，呈现为椭球状。
3. 数据收集
- 按下超声控制面板上的 新建检查 按钮以开始新的实验。
- 通过在 姓氏 中输入 Pipe Flow，在 名字 中输入当天日期，在 患者编号 中输入测试编号，创建一个新的“患者”。
- 创建“患者”后，超声扫描将开始，直至达到预设的 1000–1500 张图像上限，之后将开始新一轮扫描循环。在达到预设最大图像数之前，随时可按下超声控制面板上的 冻结 按钮两次，以重新启动扫描。
- 当获取到一组质量良好的超声图像后（即 种子粒子图像清晰且种子粒子密度足够），按下超声控制面板上的冻结按钮以停止图像采集。
- 按下超声控制面板上的 动态图像循环 按钮。使用超声控制面板上的 首帧循环 旋钮选择该组图像中的第一帧，使用 末帧循环 旋钮选择该组图像中的最后一帧，以选定待分析的超声图像集。
- 按下超声控制面板上的 图像存储 按钮，以保存选定的超声图像集。
- 按下超声控制面板上的 归档 按钮，并使用鼠标光标选择 结束检查。系统将提示用户选择要保存到本地硬盘的图像或动态图像循环。选择感兴趣的动态图像循环后，退出本次检查。
- 按下超声控制面板上的 归档 按钮，并使用鼠标光标先选择 更多，然后选择 磁盘管理。磁盘管理功能将把已保存的动态图像循环传输至运行 PIV 软件的计算机。
4. 转换文件类型
- 超声图像以医学数字成像和通信（DICOM）文件类型存储在超声设备上。为了能够被PIV软件打开和读取，必须将DICOM文件转换为图片文件。目前使用运行DICOM2JPG.m的Matlab脚本将DICOM文件转换为联合图像专家小组（JPEG）文件格式。
- 随后使用LaVision公司的DaVis软件分析JPEG超声图像。
5. 使用 DaVis 计算位移场 D(x, y)
- 在电脑上双击 DaVis 图标。选择 新建项目，然后选择 PIV。
- 在工具栏中选择 导入图像，然后选择 通过编号文件导入。在下拉菜单中找到存储 JPEG 超声图像的文件夹，并双击该序列中的第一幅图像。这将导入此编号序列中的所有超声图像。
- 通常需要定义一个图像掩膜，以在待处理的超声图像中隔离出感兴趣区域（ROI）。对于管道流，掩膜用于定义管壁之间的区域（即，流体区域）。
- 进入 DaVis 的主控制面板，选择位于 当前项目 下方、包含已导入图像的标签页，然后选择标有 批量处理 的标签页。这将启用 DaVis 的矢量处理窗口，用于对导入的超声图像进行批量处理。
- 在操作列表中，通过 PIV-时间序列 树形结构选择 矢量计算参数，并选择用于矢量处理的参数。如果使用掩膜，请在 矢量计算参数 菜单中勾选 数据范围 = 使用掩膜区域。请注意，矢量计算参数 的最优选择取决于流动几何结构、流动特性、图像分辨率、示踪粒子密度以及所需的定量流动分析。10
对于管道流动测量，通常获得最佳结果的参数设置为：采用多步处理方式， interrogation 区域从 32 × 32 像素2 逐步减小至 8 × 8 像素2，重叠率为 50%。相对矢量范围限制设为 ±all(window size/2)，绝对矢量范围限制设为 ±5 像素。最后，使用 3 × 3 像素2 的中值滤波器以抑制噪声并平滑矢量场。 - 在 批量处理 界面左侧，选择要处理的图像总数，然后点击 开始处理。这将利用互相关算法计算连续超声图像之间的位移场 D(x,y)。
6. 分析矢量场
- 进行后处理和数据分析时，将DaVis中的EPIV矢量场导出为.txt文件。具体操作为：在项目界面中，选择JPEG图像分支下的矢量位移分支，然后在工具栏中选择导出选项卡，选择文件类型为ASCII .txt，选择或创建导出文件夹，最后点击导出。
- 导出的矢量场文件命名为Bxxxxx.txt，其中00001≤ xxxxx ≤ 99999，B表示缓冲区。每个文件包含四列数据：(1) 矢量在图像中的x位置，(2) 矢量在图像中的y位置，(3) 位移的x分量（即流向位移），(4) 位移的y分量（即壁面法向位移）。通过MATLAB打开并处理Bxxxxx.txt文件，首先计算速度场，需已知图像对之间的时间步长ΔT[s]和图像放大倍数M[米/像素]，即u(x,y) = MD(x,y)/ΔT。图像间的时间步长ΔT = 1/fps + D(x,y)/B，其中B[像素/秒]表示超声探头扫过图像宽度所需的时间。本研究中，M = 77[μm/像素]，fps = 49.5[1/秒]，B = 25,047[像素/秒]。随后，计算系综平均速度矢量场、平均速度的壁面法向分布以及其他感兴趣的流动参数。（参见“代表性结果”部分。）
访问受限。请登录或开始试用以查看此内容。
结果
瞬时回波粒子图像测速法（EPIV）矢量场如图3所示。矢量图显示每隔四列的矢量，背景颜色等高线图对应于速度大小。对1000幅瞬时EPIV矢量图进行集合平均后得到的平均矢量图如图4所示。与管流特性一致，速度矢量主要沿流向方向，最大速度出现在管道中心线处，而管道壁面处的速度减小至零。图5显示了速度大小的均方根（rms）脉动值。由于在哈根-泊肃叶流动中，均方根速度应严格为零，因此非零的均方根速度反映了EPIV测量中的噪声水平。靠近上壁面处较高的均方根值是由于超声波束在管壁处发生强烈的反射和折射，导致该区域图像强度较高所致（见图2）。这些靠近壁面的高强度信号会掩盖粒子信号，从而引起测量误差。通过在行方向（水平方向）对集合平均后的矢量图进行平均，得到流向平均速度的壁面法向分布曲线，如图6所示。实心黑线表示在给定实验条件下哈根-泊肃叶（层流管流）流动的理论平均流向速度分布。EPIV测量结果与理论哈根-泊肃叶分布的吻合程度在管道中心线附近最佳，在管道壁面附近最差，其中与理论偏差最大的区...
访问受限。请登录或开始试用以查看此内容。
讨论
本文介绍了一种回波粒子图像测速（EPIV）系统的操作规程，该系统能够获取在光学不透明流体中或透过光学不透明几何结构的二维速度场。EPIV 的实际应用非常适合研究工业和生物流动系统，因为在许多应用场景中都存在不透明流体的流动。本文所展示的特定系统是专门用于研究在木质纤维素乙醇生产过程中所使用的液化生物质流体的流动特性。通过管道流动中的代表性测量结果，验证了 EPIV 系统的能力。特别是，从 EPIV 矢量场中计算得到了平均速度和均方根速度分布，并成功实现了对哈根-泊肃叶（层流）管道流动的测量与量化。EPIV 的局限性在于其固有的低帧率（受限于商用超声成像系统的成像能力）以及较低的空间分辨率，这限制了可测量的速度范围和瞬态流动行为。最后，尽管我们力求使本文内容自成体系，但仍建议查阅商用超声设备的用户手册7 和 PIV 软件的用户手册8 以获得完整信息。此外，读者可参考文献9 和10，分别获取关于超声成像基础和粒子图像测速技术的全面综述。
访问受限。请登录或开始试用以查看此内容。
披露
作者无任何利益冲突需要披露。
致谢
作者们衷心感谢美国国家科学基金会（项目编号：CBET0846359）提供的支持，项目监管人为 Horst Henning Winter。
访问受限。请登录或开始试用以查看此内容。
材料
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 超声波仪器 | GE | Vivid 7 Pro | |
| 线性超声探头阵列 | GE | 10 L | |
| 直流水泵 | KNF | NF 10 KPDC | |
| 矢量处理软件 | Lavision | DaVis 7.2 | |
| 后处理软件 | Mathworks | MATLAB 7.12 | |
| 丙烯酸管路 | McMaster-Carr | 8486K531 | |
| 超声耦合凝胶 | Parker | Aquasonic 100 |
参考文献
- White, F. M. Fluid Mechanics. , McGraw Hill. New York, New York. (1994).
- Hak, M. G. ad-el Flow Control: Passive, Active, and Reactive Flow Management. , University Press. Oxford. (2000).
- Kim, B. H., Hertzberg, J. R., Shandas, R. Development and validation of echo PIV. Exp. Fluids. 36, 455-462 (2004).
- Zheng, H., Liu, L., Williams, L., Hertzberg, J. R., Lanning, C., Shandas, R. Real time multicomponent echo particle image velocimetry technique for opaque flow imaging. Appl. Phys. Lett. 88, 261915(2006).
- Beulen, B., Bijnens, N., Rutten, M., Brands, P., van de Vosse, F. Perpendicular ultrasound velocity measurement by 2D cross correlation of RF data. Part A: validation in a straight tube. Exp. Fluids. 49, 1177-1186 (2010).
- Poelma, C., Mari, J. M., Foin, N., Tang, M. -X., Krams, R., Caro, C. G., Weinberg, P. D., Westerweel, J. 3D Flow reconstruction using ultrasound PIV. Exp. Fluids. 50, 777-785 (2011).
- GE VINGMED ULTRASOUND A/A. Vivid 7/Vivid 7 PRO User’s Manual. , FC092326 edition, GE VINGMED. Horten, Norway. (1988).
- DaVis Software for Intelligent Imaging [Internet]. , LaVision. Michigan. Available from: http://www.lavision.de/en/download.php?id=3 (2013).
- Szabo, T. Diagnostic Ultrasound Imaging: Inside Out. , Elsevier Academic Press. Burlington, MA. (2004).
- Raffel, M., Willert, C., Wereley, S., Kompenhans, J. Particle Image Velocimetry: A Practical Guide. , Springer-Verlag. Berlin Heidelberg New York. (2007).
访问受限。请登录或开始试用以查看此内容。
重印与许可
