Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
Click here for the English version

Engineering

湍流应用三维粒子跟踪测速:一喷流病例

Published: February 27, 2016 doi: 10.3791/53745

Summary

这里描述了基于高速摄像头,四视图分路器上的三维粒子跟踪测速(3D-PTV)系统。该技术是下游在雷诺数Re≈7000施加到从圆形管在十个直径附近的喷流。

Abstract

3D-PTV是一个定量的流量测量技术,该技术的目的是跟踪使用图像序列的立体记录三个维度的一组颗粒的拉格朗日路径。其基本组件,功能,约束和3D-PTV拓扑结构由高速摄像头,支持四视图分路器的优化技巧进行了描述和本文中讨论。该技术是在重新 ≈7000施加到圆形射流的中间流场(5 <X / D <25)。在欧拉帧拉格朗日流动特征和湍流量围绕射流原点的下游10的直径,并在从所述喷射芯各种径向距离估计。拉格朗日属性包括轨迹,速度和选定的颗粒的加速度,以及流路,其从所述的Frenet-Serret等式获得的曲率。在位于十的横平面围绕喷射心轴的三维速度和湍流场的推定直径射流的下游与文献相比较,并在从喷射芯各种径向距离得到大规模流向速度运动的功率谱。

Introduction

湍流喷流是在工程应用中普遍存在。这种流动的详细描述是在从大型环保排放系统涵盖电子微尺度器件的实际问题广泛至关重要。由于其对一些广阔的应用前景影响,喷气流量已经深入研究了1 - 4。几个实验技术,包括热线测速仪4 - 8,激光多普勒测速仪(LDV)4,9 - 12,和粒子图像测速(PIV)12 - 16,都被用来描述飞机在大范围的雷诺数和边界流向条件。近日,采用3D-PTV研究湍流/非湍流界面的喷气流17,18已做了一些研究。 3D-PTV是特别合适的描述复杂湍流音响的技术从不同的角度的视场。它允许粒子轨迹在使用多视点立体参考拉格朗日框架体积内重建。该技术最早是由19张和介绍了由RACCA和杜威20进一步发展。此后,许多方面做了改进的3D-PTV算法和实验设置21 - 24。与这些成绩和以往的作品,该系统已被成功地用于研究各种流体的现象,如在4米x 2米×2米25,室内气流场26的一个结构域的大型流体运动,脉动流27和主动脉血流28

一个三维PTV测量的工作原理包括数据采集系统设置,记录/前处理,校准,3D的对应,时间跟踪和后处理。精确的校准允许的粒子位置的精确检测秒。在三个以上的图像视图检测到的粒子的对应允许基于对极几何形状的三维粒子位置的重建。从连续的图像帧的连锁导致定义粒子轨迹S(t)一个时间跟踪。所述3D-PTV系统的优化是必需最大化多粒子可追溯性的可能性。

优化的第一步是获取合适的数据采集系统,包括高速摄像机,照明光源和播种颗粒的功能。与询问体积的大小沿所述照相机分辨率定义像素尺寸,因此,所需要的籽晶颗粒大小,这应该是比单一像素大。检测的颗粒的质心通过取由亮度21加权粒子像素的平均位置具有子像素精度估计。摄像机的帧率是密切associat编与雷诺数和检测到的颗粒链接的能力。较高的帧速率可以用于解决更快流量或粒子的数量较多,因为当图像之间的平均位移超过所述颗粒的平均分离的跟踪变得更加困难。

快门速度,光圈和灵敏度有三个因素在图像采集来考虑。快门速度应该足够快,以最小化颗粒,从而降低颗粒质心位置的不确定性围绕模糊。相机光圈应调整到询问体积的景深,以减少的容积之外探测粒子的概率。由于摄像机的最大感光度是固定的,作为帧速率的增加,所需要的必要的光以照亮颗粒应相应增加。不像PIV,复合光学设置以及高功率激光器没有严格要求在3D-PTV,只要光源足够粪便从示踪粒子到相机羊羔。连续LED或卤素灯属于跳过同步21,需要良好的成本效益的方案。

在3D-PTV,像其他光流的测量技术,示踪粒子速度被假设为在本地瞬时流体速度29。然而,这仅仅是对空的直径和惯性的理想示踪剂的情况下示踪颗粒应足够大,以由照相机捕获。有限微粒的保真度可以通过斯托克斯数 T来确定, 颗粒的弛豫时间尺度的比率和感兴趣紊流结构的时间尺度。一般来说,S T应当比1小很多对于 T≤0.1流跟踪误差在1%以下30。在深入讨论可以在梅等人发现29 - 31 例如 50-200微米)32,而较小的颗粒( 例如 ,1-50微米)33,34可以用高功率激光可使用( 例如 80-100瓦CW激光)。具有高反射率对于一个给定波长的光粒子,就像在卤素灯银涂层,可以放大自己的印记成图像。接种密度为成功的3D-PTV测量的另一个重要参数。少数颗粒导致低轨迹数,而颗粒的数量过多会导致在建立对应和跟踪歧义。在建立对应歧义包括重叠粒子和检测沿定义极线多名候选人。在跟踪的过程中,模糊性由于高seedin是因为颗粒的相对短平均间隔的发生g密度。

第二步骤是在记录/预处理最佳设置,以提高图像质量。照相的设置,如增益和黑电平(G B)的,在优化图像质量方面发挥了重要的作用。黑电平限定在图像的最暗部分的亮度水平,而增益放大的图像的亮度。对G&B水平的轻微变化可以显著影响可追溯性的可能性。事实上,高G&B可能过亮的图像,并最终损坏相机传感器。为了说明这一点,在流程再造G&B水平的影响也研究过这个文章。在预处理步骤中,图像被过滤,用一个高通滤波器,以从颗粒强调光散射。像素大小和灰度调整以最大化询问体积内的颗粒探测。

t“的>优化的第三步是在立体成像,这是基于对极几何,照相机参数(焦距,主点,并且失真系数),以及折射率的变化的精确校准。此过程是必不可少的,以尽量减少在3D基准目标点的重构误差。极几何使用相对距离和倾斜角度(摄像机和询问体积之间)从目标图像。沿着通过询问体积相机视图折射率变化可以基于质谱的方法来考虑 21上 。在该实验中,以规则分布的目标点的三维阶梯状结构被用作靶。

在3D-PTV实验中,虽然只需要两个图像,以确定一个三维粒子位置,典型地更多的摄像机被用来减少模糊度21。昂贵的设置有多个高速摄像机另一种是viEW分离器,由Hoyer 等人提出。35,用于使用3D-PTV的和由Gulean 等人 28的生物医学应用最近应用。该视图分路器由一个金字塔形的镜子(关于此主镜)和四个调节后视镜(关于此副镜)组成。在这项工作中,一个四视图分离器和一个摄像机被用来从四个摄像机模仿立体成像。该系统是用于表征一个管喷射的中间流场与直径d H = 1cm到重新 ≈7,000拉格朗日和欧拉在约14.5-18.5直径下游射流原点帧。

Subscription Required. Please recommend JoVE to your librarian.

Protocol

1.实验室安全

  1. 查看所选照明光源的安全指导方针( 激光,LED产业,卤素)。
    注意:在这个实验中,一组五个250瓦卤素聚光灯被用作照明。该光源的基本安全性和推荐方面描述如下。
    1. 避免与卤灯,它在高温下(〜3000色温)下操作的直接接触。
    2. 采集数据时,为了避免加热的条件下考虑流动保持光ON而已。
    3. 远离可燃所有材料靠近光源,包括任何种类的纸张。

2.实验装置

  1. 选择合适的镜头
    1. 选择低像差镜头,以避免校准问题。推荐镜头类型长焦或微透镜。
    2. 确保镜头涵盖的obj中所需的字段的视图(FOV)等的距离,O,通过估计所需的放大倍率,
      注意:倍率是相机芯片到FOV的长度的比例,和对象距离可以计算为O = F(1 / M + 1),其中f是透镜的焦距。在该实验中,摄像头芯片的长度为20.34毫米和相应的视场,或主镜,为50毫米的有限物距ø≤250毫米。 (物距被抑制由于在摄像机和视图分离器被安装在滑动件的有限长度。)的倍率为M = 20.34 / 50 = 0.41,和与对象的距离的给定范围近似焦距为f ≤72.7毫米。因此,微透镜具有60mm的焦距用于以f / 2.8D的焦率。
  2. 安装和调整与分流查看相机。
    1. 通过滑动MIR水平的审讯卷的主镜的中心沿着垂直安装后ROR和修正柱支架镜像。注意,此步骤是在安装次级镜之前进行。
    2. 安装摄像机,并设置图像重合的中心与主镜的中心。
    3. 调整摄像机视图的宽度和高度刚好盖住通过控制感兴趣区域(ROI)的区域主镜的录音软件设置。这个过程是为了减少图像尺寸和图像噪声。注意:在这个实验中,主反射镜的尺寸为5×5厘米2(1728点¯x1728像素)。
    4. 设计一个(可定制)3D校准目标。它应包围整个调查音量。确保分路器的每个视图捕获所有目标标记,以便统一校准。
      注:在本演示中,目标是3D印刷用Vero细胞回塑料。它有一个台阶状与尺寸35×35×30mm的3,用直径1毫米的白色目标点分隔符编2,5 mm开始,在5mm和10mm的垂直,流向和翼展方向。目标的精确几何形状是关键的,因为它伸入一个校准模型和摄像机位置相对于所述水槽。
  3. 将校准目标进审讯音量。
    1. 对相机面临的一个可调节高度的平台上标定目标位置。
    2. 通过调整目标平台的高度相匹配与询问体积的中心校准目标的中心的高度。
      注意:在本实施例中,校准目标的中心标记与射流喷嘴的中心20厘米的高度调平。气泡水平仪可以用于水平目标。
  4. 安装和调整四视图分路器的次镜。
    1. 找到从确保充分捕捉审讯的距离主镜。这为0.2μm此演示(FIGURË1)。
    2. 在它的大致位置时,在来自主镜的各侧摄像机视图大致与副镜对准贴装副镜。通过将其固定到辅助镜的垂直安装柱固定副镜。
    3. 重复此过程,对于其他三面镜子。检查所有次级镜的几何对称性相对于所述主​​镜。
    4. 使通过调整次级反射镜的反射镜安装,以确保每个的四个视图包围整个校准目标最终的调整。检查镜'的位置和角度的一个有效的方法是使用激光指示器来可视化各视图的图像的路径。
    5. 检查用于通过移动一个反射镜子图像的重叠。
      注意:如果只有一个视图改变,那么重叠区域可以忽略不计。否则,重复步骤2.4.2至2.4.5,直到重叠区域最小化。
  5. 将光源(S)直接面对审讯音量。确保将相机调光,以避免相机传感器损坏时覆盖。
    1. 检查光源是在整个调查体积均匀分布的。
    2. 增强的光强度,如果需要的话,由下光源直接放置一个放大镜。注意:在这个实验中,使用焦距f 0 = 450nm 10mm的平凸放大镜加紧的照明。

3.建立优化

  1. 打开并调整相机设置以提高图像质量。
    1. 调整镜头倍率,直到通过主镜的反射聚焦同样在所有四个视图次要镜子。
    2. 检查从视图分离器的图像是对称的,并且通过从四个视图观察校准目标图像的对称性捕获询问体积。
    3. 调整f数来捕获从摄像机最接近的和最远的校准目标点。
      注意:这使得相机捕捉示踪粒子只有在审讯量的深度。在这个例子中,F数是11。
    4. 设置所需的帧速率550赫兹(这取决于具体的应用,见简介)和录音软件相应最大限度的光敏感度。
  2. 通过一个直播相机视图观察分离器的每个视图的颗粒密度差检查主镜的每个视图照明。
    注意:如果多个光源用来照亮所述询问体积,它是可能具有在分离器的每个视图的差异。在该实验中,顶部的两个次级镜,因为照明来自顶部接收较少的光。在水槽底部的使用平面镜可以帮助减少光变跨越的意见。
  3. 采用3D-PTV光源前,请关闭房间内的灯光背景。
  4. 调整相机的G&B的水平,以更好地捕获粒子光散射。记录与各种G&B水平几个短序列,发现通过观察粒子轨迹的分布和密度的最优的一个。
    注意:在这个实验中,G&B水平的范围是0-500,和黑色(B)的电平设定为500,以亮调光器的光散射,而,增益(G)设定介质,300,以适度放大图像信号和避免过度使图像变亮。

4.校准

  1. 放置在侦查卷中的校准目标添加示踪粒子之前,采取一些校正图像。使用调光器光源( LED闪光灯光)照射目标。
  2. 除以校准图像分成四个独立的子图像,使一个文本文件包含指ENCE协调目标标记的位置。 OpenPTV软件(http://www.openptv.net)在这里被用于此目的。
    注意:以后的处理等同于采用多个摄像机设置的用户。
  3. 点击“创建校准”选项卡中保存图像,并在步骤4.2在软件的“CAL”文件夹中获得的文本文件后启动校准过程。
  4. 点击“修改校准参数”选项卡,然后选择“校准定向参数”选项卡,每个拆分视图中心和校准目标的原点之间定义放大,旋转角度和距离。
    注意:第一个行是从原点校准目标在X,Y,Z方向上的相机传感器的距离。第二行表示角度,单位为弧度,绕x,Y,Z轴。接着,一个3×3的数据表示旋转矩阵。然后,两个下列各行是x的针孔距离和X,Y,Z方向上的位置。
  5. 点击“检测”和“显示初始猜测”调查的“猜测”点与检测到的目标点相匹配。
  6. 所有四个视图,直到'猜'点重复步骤4.4与所设定的校正图像的对齐。
  7. 点击“取向”重构询问体积的方向。
    注意:校准可以通过调整透镜畸变和仿射变换得到改善。现在,调查体积校准并准备处理数据。见笔者的附加 ​​说明论文36校准过程。

5.流量设置/数据采集

  1. 来自帧速率的估计在每个摄像机视图捕获颗粒的最大数量二维最大流速。在这个演示中,给定速度为U≈0.4米/秒,帧速率为550赫兹和4〜×4×4 立方厘米体积审讯。这导致〜每帧1000个粒子。
  2. 打开相机,在步骤3中得到的最优化设置。
  3. 添加播种颗粒和等待数平均停留时间,从而使流动达到稳定状态。如果需要加入更多的粒子,但避免高接种密度,在步骤5.1,这可能会导致模糊度估计。
    注意:在该实例中,〜1.6克为1.1g / cm 3的密度为100微米银包覆空心陶瓷球体被用作种子用于流体介质(2×0.4×0.4米3)。
  4. 记录流图像的期望数目。
    注意:在这个实验中,在550赫兹9000图像用记录软件被抓获。如果相机和/或视图分路器移动(即使是轻微的运动可以严重影响结果)重复步骤2.4至5.3。

6.数据处理(通过OpenPTV软件)

  1. 除以在步骤5.4获得分成四个独立的子图像的原始图像。
  2. 点击“开始”选项卡下的“初始化/重新启动”,从四个视图加载初始图像。
  3. 右键单击“运行”目录,然后单击“主参数的摄像头,折射率,粒子识别,图像序列,观察量和标准对应的个数来控制数量。
    1. 定义的用于“常规”选项卡下的实验相机(意见)的数量。在这个实验中,设置的摄像机的数量为4。
    2. 定义沿着“折射率”选项卡下的摄影机视图折射率。
    3. 定义最小和像素检测以及灰度值的阈值MAX数量,优化了“粒子识别”选项卡下的粒子检测的所有四个视图数。德像素的最小和最大号tection和灰度阈值确定用于颗粒探测像素大小和亮度级别。它消除了噪音和颗粒失焦。
    4. 定义'的序列处理参数“下图像处理的数量。
    5. 定义“观察音量”选项卡下的观察量。
    6. 定义下的“标准对应”,包括立体匹配总频带参数(毫米)对应的相关性。
  4. 点击“预处理”标签下的“高通滤波器”。这加剧在所有四个观点粒子光散射。
  5. 点击“粒子探测”,以确定检测到的颗粒在所有四个视图亚像素级的重心。重复步骤6.2和6.3,直到类似于在步骤5.1计算颗粒的预期数检测颗粒的数量。
  6. 点击“对应关系”,建立立体correspondenCES在每个视图。
    注意:要重建检测到粒子的三维位置,对应至少应该从三个意见确定。
  7. 点击“3D位置”来获得基于所述校准检测粒子三维位置。
  8. 点击“不带显示屏的序列'重复步骤,从6.4到6.7的所有图像序列的过程。
    注意:这将创建包含检测粒子与制表符分隔文本文件格式框架总结每个图像设定'rt_is'文件。
  9. 右键单击“运行”目录,然后单击“跟踪参数”来定义半径球体,( dvxmin和dvymin毫米/帧)的参数,以搜索候选粒子进行跟踪。
  10. 点击“跟踪不显示'来定义步骤6.7获得重建的颗粒的标识(ID)。
    注意:用于使用一个四帧predic跟踪相关的相邻帧的序列器预测-校正方案24。这个过程创造了包含在帧中检测到粒子的跟踪信息的每个图像集ptv_is文件;前两列分别示出了粒子的ID中的前一帧和下一帧。
  11. 点击“显示轨迹”可视化轨迹中的每个摄影机视图。

7.后处理(祈愿)

注:范围和后处理的类​​型取决于个人的需求和它,因此,可定制。这里,指向基底计算简单说明作为例子。

  1. 获得拉格朗日帧数据(通过Matlab的)。
    1. 提取每个粒子及其相关的ID从ptv_is文件的三维位置。它允许连接的图像序列中检测粒子重建轨迹。
    2. 计算从给定帧速率的每个轨迹的速度和颗粒的加速度。在此示范,粒子的速度和加速度由低通用移动三次样条34,37滤波该位置信号来计算。
    3. 使含三维位置,三维速度,加速度3D和时间戳以及每个轨迹的轨迹ID字段的结构阵列格式。在该数据格式中, 结构体数组的长度代表轨迹号。
  2. 获得欧拉帧数据(通过Matlab的)。
    1. 结构体阵列(步骤7.1.3)转换到使用每个颗粒的时间戳的时间之一。这将在步骤7.1.3得到同样的结构排列的结构,但结构体数组的长度现在代表帧号,这是在这个实验9000元。
    2. 内插的时间结构数组为每个时间段三维网格,以获得欧拉坐标瞬时速度场。在这个演示中,有趣的GridDataction在Matlab用于执行内插。

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

的照片,并将设置的示意图示于 12。校准目标,在基准标记上反射的视图分离器和3D校准重建在 3中示出。公认的校准目标的RMS为7.3微米,5.7微米和141.7微米的流向X,翼展方向y和深度z方向。在z -协调的相对较高的RMS是由于相对于那些在其它方向和相比x和y坐标与z轴四个视图相对小角度的减小目标点。在每个在任何给定时刻的四个视图所检测的粒子为10 3的数量级,其中所检测的微粒,成功的三维重建的数量减少到大约一半,由于这一事实,即仅在交叉区域AR颗粒É捕获。 视频 1示出了从四个视点分离器的喷流的一个高速视频样本。

在中间场区域周围并穿过x / D H = 16以径向距离平面的四种代表性粒子轨迹的样品R / D H = V0,1.5,3从喷射芯在 4中示出。正如所料,在给定的时间间隔(Δt的≈1秒)更长的轨迹围绕喷射芯观察。在射流(转/天小时 ≥2)的边缘,示踪颗粒具有短和更复杂的轨迹。 5显示了所有成功重建粒子轨迹穿过x /天小时 = 16平面。在所选域的粒子速度呈现分布广泛,从近0〜0.6 单位 j,其中ūĴ6示出了一个颗粒的穿越的x /天小时 = 16围绕喷射芯平面上的情况下, 6 B,6 C和 6 d显示粒子的轨迹,速度和加速度的规格化时间的函数的3个分量。值得强调的是​​,当地的粒子加速可以是几十倍的标准重力。粒子轨迹允许通过所谓星期五获得所述粒子轨迹的特定功能ENET-Serret框架。它描述了沿S(T)的正交向量(切线,法线,副法线)的变化。特别相关的是曲率,κ,这是曲率半径的倒数,ρ,并定义为:

式(1)

哪里公式3 =博士/ ds是轨迹的切线单位矢量,r是位置矢量的粒子的(欧几里德空间)作为时间的函数,它可以写为一个函数, ,R(S)= R(T (多个))。曲率,κ,计算为跨越x/d H = 16和x/d H = 17架飞机中的所有粒子。平均曲率, 公式3作为从喷射芯r中距离的函数的计算公式为: 公式(2)

其中,ΔR= 0.2D h的在这里使用。 7示出了公式3 = F(r)的由天小时归一化。它显示了相对低的,几乎恒定公式3通过管道的圆形横截面,转/天小时 ≤0.5限定的区域内。在从在x / D H = 16平面射流核心较大的距离, 公式3单调增加。类似的趋势是在在x / D H = 17飞机获得的,但减少公式3射流核心(R /天小时 ≥0.5)之外。值得强调的是​​,这种流动FEATURE只能与3D-PTV技术来推断。根据各级的G&B设置数据的质量在联颗粒,以表1所示三维重建颗粒的其余部分的比计进行评估。在最高链路比率在300 500 G&B设置被观察到。

欧拉的流动性可以通过网格内插,它模仿的3D粒子图像测速技术(3D-PIV)来实现。重要的是要注意,由于每次都在追踪的相对低的颗粒,需要真正模拟PIV质量为欧拉描述帧的一个显著更大的数字。这是高阶的统计数据( 例如 ,湍流强度和雷诺应力)的估计更为关键。在关于各种G&B水平射流核心流向速度在 8中示出。测量是与theoret相比iCal的行为:

公式3

其中U 0(x)是在喷射芯的流向速度,B≈6是常数,并且x 0是虚拟原点38。该图显示设置G&B水平的相关性。 9显示在x / D H = 16平面射流的平均速度分布。

最后,在位置流向速度的大规模运动的光谱分布φ(F)R / D H = 0,0.6%,及1在x / D H = 10平面在 10中示出。巴特沃斯低通滤波器施加到速度的时间序列与截止频率f C = 200赫兹。

图1
1: 在实验装置示意图 请点击这里查看该图的放大版本。

图2
2:实验装置这说明相机的各种意见和四图像视图分离器,水槽和审讯量:(左上)俯视图,(左下)相机和视图分路器系统的后视图,(顶部中间,底部中间)的整体实验装置的侧视图,(右)放大的视图在喷射流播种颗粒。 <A HREF =“htt​​ps://www.jove.com/files/ftp_upload/53745/53745fig2large.jpg”目标=“_空白”>点击此处查看该图的放大版本。

图3
3: 校准:( )校准目标,(B)图像集从视图分路器校准目标的;( )从校准目标的基准标记的3D识别请点击这里查看这个大版本数字。

图4
4: 选择 粒子轨迹在r/d H = 0, 1。5,3。 请点击此处查看该图的放大版本。

图5
5: 粒子轨迹穿过 x / D H = 16 平面,其中速度被表示为一个颜色电平在该图所示的询问体积载(x)的/dħ之间图1 (14.5,18.5),y/d ^ h 图1 (-2,2),并且z/dħ.JPG“/>( - 2,2),其中(x,Y,Z)=(0,0,0)位于射流原点的中心沿着个别轨迹的速度,由散速度归一化。 ü0,被示为色彩等级。 请点击此处查看该图的放大版本。

图6
6:( )粒子的轨迹,(B)置换,(C)速度和任意粒子( )加速,请点击此处查看该图的放大版本。

图7 /> 7: 颗粒的曲率:图表显示x/d意味着粒子的曲率从飞机喷气核心的径向距离的函数h = 16和x/d H = 17, 请点击这里查看一个更大的版本这个数字。

图8
8:在(X)/d 小时 内射流核心流向速度 图1 (15,18)关于各种G&B水平。三G&B水平都包含(300 500(最佳),300及250,100及250)。745fig8large.jpg“目标=”_空白“>点击此处查看该图的放大版本。

图9
9: x / D H = 16 的流向速度分量的非维分布 请点击此处查看该图的放大版本。

图10
10: 在位于 R / D h 的点的流向速度分量的 功率谱 φ(F) = 0 x /天 小时 (喷射芯),0。6, 1 = 16 飞机。 请点击此处查看该图的放大版本。

视频1
视频1:从四视图分路器的喷射流的视频样本,比在550 fps的获得实际速度慢10倍右键点击下载 )。

表格1

ŧ能够1:三个相互联系粒子对 3 D-重建颗粒 的其余部分 在不同的G&B水平 G&B级别包括:(100 250 300 250 300 500)。

Subscription Required. Please recommend JoVE to your librarian.

Discussion

3D-PTV具有很大的潜力,以解开各种湍流,如在低层大气中25,室内空气分布26,或脉动流的大规模动荡运动的复杂的物理主动脉拓扑28在许多其他之中。但是,它的优点和局限性,以及经验的理解是必不可少最大限度地发挥其潜力。试错的初步测试和优化设置,包括帧速率,照明光源,G&B级和图像滤波方法详尽的迭代,直接与重建一套( 例如示踪)粒子的拉格朗日路径的能力相关。值得注意的是,关键的协议的步骤,如这里证明的,是G&B水平的调整和视场的照明(的卤素灯的灯光,放大透镜组合,并从水槽的底部反射镜)是重要的。

ENT“>这些调整有助于调查的四个视图内优化光散射。查明高保真测量实验设置后,彻底修改和故障排除提出,要计算基于帧速率精确轨迹的最大数量,照相机的分辨率以及调查体积的大小。虽然捕获颗粒的数目可以与更高的帧速率可以增加,这是值得注意的是,相对于PIV在3D-PTV跟踪的粒子数低得多。3D的最大的潜在-PTV是其独特的描述多粒子的拉格朗日路径的能力。在这个演示中,视图分路器设置被实施,以避免使用多个摄像头膨胀,但是,需要注意的是这种设置需要更高的摄像头是非常重要的分辨率和限制样品体积的大小。

在这项研究中,一圆形射流的中间场特征是肛门yzed与3D-PTV技术。该方法允许获得从欧拉和拉格朗日帧流动的重要特征。特别是,颗粒作为径向距离的函数的平均曲率使用粒子轨迹的拉格朗日特征表征首次两个截面平面。公认的校准目标的RMS 7.3微米之间的范围内,以141.7微米的流向和展向的方向。虽然在跨度方向由于在z方向上的次小角度,高相对误差可能不被完全克服,它可以进一步通过在z方向加入更多的目标点,如使用2D校准目标在各个位置降低(多平面校准)。

总体而言,3D-PTV是,可以在其他的一些问题,包括时间相关的流量或活性标量的动态施加一个有用的技术。例如,它可以是给STU非常有用DY动荡和种水生环境之间的相互作用。

Subscription Required. Please recommend JoVE to your librarian.

Disclosures

作者宣称,他们没有竞争的经济利益。

Acknowledgments

这项工作是由机械科学与工程学院,伊利诺伊大学厄巴纳 - 香槟分校,部作为莱昂纳多P.查莫洛的开办包的一部分支持。

Materials

Name Company Catalog Number Comments
Mikrotron 4CXP 4 lanes CXP-6 CoaXPress ImageOps CAMMC4082 High-speed camera
Active Silicon FireBird CoaX Frame Grabber ImageOps FBD-4XCXP6 Frame Grabber
100 μm silver-coated hollow ceramic spheres Potters Industries LLC AG-SL150-30-TRD Seeding Paritcles
StreamPix6 Upstate Technical Equipment CO.,INC MISNOR-STP-6-S-CL Camera appliation
Four-view splitter Photrack AG Customized part and necessary if performing 3D-PTV with one camera
250 Watts Spotlight Halogen General Electrics 23719 Light source
OpenPTV (Software) OpenPTV (http://www.openptv.net) Open source particle tracking software (Note: available as a service for anyone who wants to use it without all the installation mess or computer power availability problems).

DOWNLOAD MATERIALS LIST

References

  1. Wygnanski, I., Fiedler, H. Some measurements in the self preserving jet. , Cambridge university press. (1968).
  2. Rajaratnam, N. Turbulent jets. , Elsevier. (1976).
  3. Panchapakesan, N., Lumley, J. Turbulence measurements in axisymmetric jets of air and helium. Part 1. Air jet. J Fluid Mech. 246, 197-223 (1993).
  4. Hussein, H. J., Capp, S. P., George, W. K. Velocity measurements in a high-Reynolds-number, momentum-conserving, axisymmetric, turbulent jet. J Fluid Mech. 258, 31-75 (1994).
  5. Yule, A. Large-scale structure in the mixing layer of a round jet. J Fluid Mech. 89, 413-432 (1978).
  6. Yule, A., Chigier, N., Ralph, S., Boulderstone, R., Venturag, J. Combustion-transition interaction in a jet flame. AIAA Journal. 19, 752-760 (1981).
  7. Quinn, W. Upstream nozzle shaping effects on near field flow in round turbulent free jets. Eur J Mech B-Fluid. 25, 279-301 (2006).
  8. Mi, J., Nathan, G. J., Luxton, R. E. Centreline mixing characteristics of jets from nine differently shaped nozzles. Exp Fluids. 28, 93-94 (2000).
  9. Karlsson, R. I., Eriksson, J., Persson, J. LDV measurements in a plane wall jet in a large enclosure. DTIC [Internet]. , Available from: http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADP008905 (1992).
  10. Liepmann, D., Gharib, M. The role of streamwise vorticity in the near-field entrainment of round jets. J Fluid Mech. 245, 643-668 (1992).
  11. Oh, S. K., Shin, H. D. A visualization study on the effect of forcing amplitude on tone-excited isothermal jets and jet diffusion flames. Int J Energ Res. 22, 343-354 (1998).
  12. Cenedese, A., Doglia, G., Romano, G., De Michele, G., Tanzini, G. LDA and PIV velocity measurements in free jets. Exp Therm Fluid Sci. 9, 125-134 (1994).
  13. Wang, H., Peng, X., Lin, W., Pan, C., Wang, B. Bubble-top jet flow on microwires. Int J Heat Mass Tran. 47, 2891-2900 (2004).
  14. Shestakov, M. V., Tokarev, M. P., Markovich, D. M. 3D Flow Dynamics in a Turbulent Slot Jet: Time-resolved Tomographic PIV Measurements. 17th Int Symp on Applications of Laser Techniques to Fluid Mechanics. , (2014).
  15. Bridges, J., Wernet, M. P. Measurements of the aeroacoustic sound source in hot jets. AIAA [Internet]. , Available from: http://arc.aiaa.org/doi/abs/10.2514/6.2003-3130 (2003).
  16. Scarano, F., Bryon, K., Violato, D. Time-resolved analysis of circular and chevron jets transition by tomo-PIV. 15th Int Symp on Applications of Laser Techniques to Fluid Mechanics. , (2010).
  17. Holzner, M., Liberzon, A., Nikitin, N., Kinzelbach, W., Tsinober, A. Small-scale aspects of flows in proximity of the turbulent/nonturbulent interface. Phys Fluids. 19, 071702 (2007).
  18. Holzner, M., et al. A Lagrangian investigation of the small-scale features of turbulent entrainment through particle tracking and direct numerical simulation. J Fluid Mech. 598, 465-475 (2008).
  19. Chang, T. P., Wilcox, N. A., Tatterson, G. B. Application of image processing to the analysis of three-dimensional flow fields. Opt Eng. 23, 283-287 (1984).
  20. Racca, R., Dewey, J. A method for automatic particle tracking in a three-dimensional flow field. Exp Fluids. 6, 25-32 (1988).
  21. Maas, H. G., Gruen, D., Papantoniou, D. Particle tracking velocimetry in three-dimensional flows. Exp Fluids. 15, 133-146 (1993).
  22. Kasagi, N., Matsunaga, A. Three-dimensional particle tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow. Int J Heat Fluid Fl. 16, 477-485 (1995).
  23. Virant, M., Dracos, T. 3D PTV and its application on Lagrangian motion. Meas Sci Technol. 8, 1539 (1997).
  24. Willneff, J. A spatio-temporal matching algorithm for 3 D particle tracking velocimetry. , Mitteilungen- Institut fur Geodasie und Photogrammetrie an der Eidgenossischen Technischen Hochschule Zurich. Zurich. (2003).
  25. Rosi, G. A., Sherry, M., Kinzel, M., Rival, D. E. Characterizing the lower log region of the atmospheric surface layer via large-scale particle tracking velocimetry. Exp Fluid. 55, 1-10 (2014).
  26. Fu, S., Biwole, P. H., Mathis, C. Particle Tracking Velocimetry for indoor airflow field: A review. Build Environ. 87, 34-44 (2015).
  27. Kolaas, J., Jensen, A., Mielnik, M. Visualization and measurements of flows in micro silicon Y-channels. Eur Phys J E. 36, 1-11 (2013).
  28. Gülan, U., et al. Experimental study of aortic flow in the ascending aortavia Particle Tracking Velocimetry. Exp Fluids. 53, 1469-1485 (2012).
  29. Mei, R. Velocity fidelity of flow tracer particles. Exp Fluids. 22, 1-13 (1996).
  30. Tropea, C., Yarin, A. L., Foss, J. F. Springer handbook of experimental fluid mechanics. 1, Springer Science & Business Media. (2007).
  31. Melling, A. Tracer particles and seeding for particle image velocimetry. Meas Sci Technol. 8, 1406 (1997).
  32. Hering, F., Leue, C., Wierzimok, D., Jähne, B. Particle tracking velocimetry beneath water waves. Part I: visualization and tracking algorithms. Exp Fluids. 23, 472-482 (1997).
  33. Biferale, L., et al. Lagrangian structure functions in turbulence: A quantitative comparison between experiment and direct numerical simulation. Phys Fluids. 20, 065103 (2008).
  34. Lüthi, B., Tsinober, A., Kinzelbach, W. Lagrangian measurement of vorticity dynamics in turbulent flow. J Fluid mech. 528, 87-118 (2005).
  35. Hoyer, K., et al. 3d scanning particle tracking velocimetry. Exp Fluids. 39, 923-934 (2005).
  36. Kim, J. -T. Three-dimensional particle tracking velocimetry for turbulence applications. , UIUC. http://chamorro.mechse.illinois.edu/3d.htm (2015).
  37. Lüthi, B. Some aspects of strain, vorticity and material element dynamics as measured with 3D particle tracking velocimetry in a turbulent flow. ETH Zürich. , Nr. 14893 (2002).
  38. Pope, S. B. Turbulent flows. , Cambridge university press. (2000).

Tags

工程,108期,机械工程,流体力学,流量的可视化,射流,拉格朗日框架,粒子跟踪测速,轨迹,动荡。
湍流应用三维粒子跟踪测速:一喷流病例
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Kim, J. T., Kim, D., Liberzon, A.,More

Kim, J. T., Kim, D., Liberzon, A., Chamorro, L. P. Three-dimensional Particle Tracking Velocimetry for Turbulence Applications: Case of a Jet Flow. J. Vis. Exp. (108), e53745, doi:10.3791/53745 (2016).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter