Name: 湍流应用三维粒子跟踪测速:一喷流病例 (视频) | JoVE Journal
Uploaded: 2016-02-27T15:00:00.000+00:00
Duration: 13 min 2 s
Description: 了解三维粒子跟踪测速技术用于湍流研究，通过高速相机与四视角分光器分析圆形管道喷流特性。该方法可追踪大量粒子时空轨迹，为流体力学提供拉格朗日视角的实验数据。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

playerTrialDesktop

本工作由伊利诺伊大学厄巴纳-香槟分校机械科学与工程系资助，属于Leonardo P. Chamorro启动计划项目的一部分。

1. 实验室安全
<ol>
	<li>查阅所选照明光源的安全使用指南（例如 激光、工业LED、卤素灯 
	注意：本实验使用一组五个250瓦卤素聚光灯作为光源。该光源的基本安全事项和使用建议如下所述。
	<ol>
		<li>避免直接接触卤素灯，因其工作温度较高（约3,000 K色温）。</li>
		<li>仅在采集数据时开启光源，以避免加热待测流体。</li>
		<li>远离光源处放置所有易燃材料，包括各种纸张。</li>
	</ol>
	</li>
</ol>
2. 实验设置
<ol>
	<li>选择合适的镜头
	<ol>
		<li>选择像差较小的镜头，以避免校准问题。推荐使用远摄镜头或微距镜头。</li>
		<li>通过估算所需放大倍数 M，确保镜头在物距 O 范围内覆盖所需的视场（FOV）。 
		注：放大倍数为相机芯片长度与视场之比，物距可由公式 O = f(1/M + 1) 计算得出，其中 f 为镜头焦距。本实验中，相机芯片长度为 20.34 mm，对应的视场（即主反射镜）为 50 mm，且物距受限于 O ≤ 250 mm（由于安装相机和视图分束器的滑轨长度有限）。因此，放大倍数 M = 20.34/50 = 0.41，在给定物距范围内对应的近似焦距为 f ≤ 72.7 mm。故本实验采用焦距为 60 mm、光圈比为 f/2.8D 的微距镜头。</li>
	</ol>
	</li>
	<li>安装并调节带视图分束器的相机
	<ol>
		<li>通过沿垂直支架滑动主反射镜并用支架固定件将其固定，使主反射镜中心与测量区域中心对齐。注意：此步骤应在安装次级反射镜之前完成。</li>
		<li>安装相机，并使图像中心与主反射镜中心重合。</li>
		<li>通过控制采集软件中的感兴趣区域（ROI）设置，调节相机视图的宽度和高度，使其刚好覆盖主反射镜。此过程可减小图像尺寸和图像噪声。注：本实验中，主反射镜尺寸为 5 × 5 cm2（对应 1,728 × 1,728 像素）。</li>
		<li>设计一个（定制的）三维校准靶标。该靶标应完全包含整个测量体积。确保分束器的每个视图都能捕捉到所有靶标标记，以便实现均匀校准。 
		注：在本演示中，靶标使用黑色 Vero 塑料 3D 打印而成，呈阶梯状，尺寸为 35 × 35 × 30 mm3，白色靶点直径为 1 mm，在垂直、流向和展向方向上的间距分别为 2.5 mm、5 mm 和 10 mm。靶标的精确几何结构至关重要，因其将用于构建校准模型，并确定相机相对于水槽的位置。</li>
	</ol>
	</li>
	<li>将校准靶标放置于测量区域内
	<ol>
		<li>将校准靶标置于一个可调高度的平台上，并使其面向相机。</li>
		<li>通过调节靶标平台的高度，使校准靶标中心与测量区域中心对齐。 
		注：在本示例中，校准靶标中心与射流喷嘴中心对齐，高度为 20 cm。可使用气泡水平仪来校准靶标水平。</li>
	</ol>
	</li>
	<li>安装并调节四视图分束器的次级反射镜
	<ol>
		<li>将主反射镜置于距离测量区域一定距离处，以确保其完整成像。本演示中该距离为 0.2 m（图 1）。</li>
		<li>将次级反射镜大致安装在合适位置，使来自主反射镜两侧的相机视图大致对准相应的次级反射镜，并通过将其固定在次级反射镜的垂直支架上来固定位置。</li>
		<li>对另外三个反射镜重复此过程。检查所有次级反射镜相对于主反射镜的几何对称性。</li>
		<li>通过调节次级反射镜的镜座进行最终调整，确保四个视图均完整包含整个校准靶标。一种高效检查反射镜位置和角度的方法是使用激光笔可视化每个视图的成像光路。</li>
		<li>通过移动一个反射镜来检查子图像之间是否存在重叠。 
		注：如果仅有一个视图发生变化，则重叠区域可忽略不计；否则，重复步骤 2.4.2 至 2.4.5，直至重叠区域最小化。</li>
	</ol>
	</li>
	<li>将光源直接面向测量区域放置。调节光源时务必遮盖相机，以防损坏相机传感器。
	<ol>
		<li>检查光源是否在整个测量体积内均匀分布。</li>
		<li>如有必要，可在光源正下方放置一个放大镜以增强光照强度。注：本实验中使用焦距 f0 = 450 mm 的平凸放大镜以增强照明。</li>
	</ol>
	</li>
</ol>
3. 设置优化
<ol>
	<li>打开相机并调整设置以提高图像质量。
	<ol>
		<li>调节镜头放大倍数，直至通过主反射镜观察时，四个次级反射镜视图中的反射像均能同时清晰聚焦。</li>
		<li>通过观察分光器四个视图中标定目标图像的对称性，检查各视图图像是否对称，并确认是否覆盖了 interrogation volume（ interrogation 体积）。</li>
		<li>调节f-number，使相机能够同时清晰捕捉到距离相机最近和最远的标定目标点。 
		注意：这可确保相机仅在 interrogation volume 的深度范围内捕捉示踪粒子。本示例中，f-number 设置为 11。</li>
		<li>根据实验需求设置所需的帧率为 550 Hz（具体值取决于实际应用，参见引言部分），并在记录软件中相应地最大化光敏感度。</li>
	</ol>
	</li>
	<li>通过实时相机画面观察分光器各视图中粒子密度的差异，检查主反射镜每个视图的照明情况。 
	注意：若使用多个光源照射 interrogation volume，则分光器各视图之间可能出现照明差异。本实验中，由于光源位于上方，上方两个次级反射镜接收到的光线较少。在水槽底部使用平面反射镜有助于减少各视图间的光照差异。</li>
	<li>在使用 3D-PTV 光源前，关闭实验室内背景光源。</li>
	<li>调节相机的 G&amp;#38;B 水平，以更好地捕捉粒子散射光。录制若干段不同 G&amp;#38;B 水平下的短序列，通过观察粒子轨迹的分布与密度，确定最佳设置。 
	注意：本实验中 G&amp;#38;B 水平范围为 0–500，黑电平（B）设为 500，以增强较弱的散射光亮度；增益（G）设为中等值 300，以适度放大图像信号并避免图像过亮。</li>
</ol>
4. 校准
<ol>
	<li>在添加示踪粒子之前，将标定目标放置于检测区域内，并拍摄几张标定图像。使用较弱的光源（例如 LED手电筒）照亮标定目标。</li>
	<li>将标定图像划分为四个独立的子图像，并创建一个文本文件，记录标定目标标记点的参考坐标位置。此处使用OpenPTV软件（http://www.openptv.net）完成此步骤。 
	注意：后续处理过程与使用多相机系统的用户完全相同。</li>
	<li>将第4.2步中获得的图像和文本文件保存至软件的“Cal”文件夹后，点击“Create calibration”选项卡以启动标定过程。</li>
	<li>点击“Edit calibration parameters”选项卡，并选择“Calibration Orientation Parameters”选项卡，以定义放大倍数、旋转角度以及每个分视图中心与标定目标原点之间的距离。 
	注意：第一行为标定目标原点到相机传感器在x,y,z方向上的距离；第二行为绕x,y,z轴的旋转角度（单位：弧度）；接下来是一个3×3的数据矩阵，表示旋转矩阵；随后两行为针孔在x和y方向上相对于图像中心的距离（单位：mm）以及焦距；最后一行为水槽玻璃相对于标定目标原点在x,y,z方向上的位置。</li>
	<li>点击“Detection”和“Show initial guess”，检查“初始估计”点是否与检测到的目标点匹配。</li>
	<li>对全部四个视图重复第4.4步，直至“初始估计”点与各标定图像中的对应点对齐。</li>
	<li>点击“Orientation”以重建检测区域的空间取向。 
	注意：可通过调整镜头畸变和仿射变换进一步优化标定结果。至此，检测区域已完成标定，可进行数据处理。有关标定过程的更多详细说明，请参见作者的学位论文36。</li>
</ol>
5. 流速设置/数据采集
<ol>
	<li>根据帧率和最大流速估算每个相机视场中捕获的粒子最大数量。在本示例中，参考流速为 U ≈ 0.4 m/s，帧率为 550 Hz， interrogation volume 约为 4 × 4 × 4 cm3，每帧图像中约包含 1,000 个粒子。</li>
	<li>使用步骤 3 中获得的优化设置开启相机。</li>
	<li>加入示踪粒子，并等待数个平均停留时间，使流动达到稳态。如需可补充添加粒子，但应避免过高的粒子密度（在步骤 5.1 中已估算），以免引起歧义。 
	注：本示例中，使用约 1.6 g 的银涂层空心陶瓷微球（粒径 100 µm，密度 1.1 g/cm3）作为流体介质（2 × 0.4 × 0.4 m3）的示踪粒子。</li>
	<li>记录所需数量的流动图像。 
	注：本实验中，使用采集软件以 550 Hz 的帧率采集了 9,000 张图像。若相机和/或视场分光器发生移动（即使轻微移动也可能严重影响结果），需重复步骤 2.4 至 5.3。</li>
</ol>
6. 数据处理（通过 OpenPTV 软件）
<ol>
	<li>将步骤 5.4 中获得的原始图像分割为四个独立的子图像。</li>
	<li>在“开始”选项卡下点击“初始化/重启”，以从四个视角加载初始图像。</li>
	<li>右键单击“Run”目录，然后点击“主参数”，以控制相机数量、折射率、粒子识别、序列图像数量、观测体积以及对应关系的判定标准。
	<ol>
		<li>在“常规”选项卡下定义实验所用的相机（视角）数量。本实验中，将相机数量设为 4。</li>
		<li>在“折射率”选项卡下定义沿相机视角方向的折射率。</li>
		<li>在“粒子识别”选项卡下设定像素检测的最小值和最大值，以及灰度值阈值，以优化四个视角中的粒子检测数量。像素检测的最小值、最大值及灰度阈值决定了用于粒子检测的像素尺寸和亮度水平，可消除噪声及失焦粒子。</li>
		<li>在“序列处理参数”选项卡下定义需处理的图像数量。</li>
		<li>在“观测体积”选项卡下定义观测体积。</li>
		<li>在“对应关系标准”选项卡下定义对应关系的相关性，包括用于立体匹配的总带宽参数（mm）。</li>
	</ol>
	</li>
	<li>在“预处理”选项卡下点击“高通滤波”。该操作可增强四个视角中粒子的光散射信号。</li>
	<li>点击“粒子检测”，以在亚像素级别确定四个视角中被检测粒子的质心。重复步骤 6.2 和 6.3，直至检测到的粒子数量接近步骤 5.1 中计算出的预期粒子数量。</li>
	<li>点击“对应关系”，以在各视角中建立立体对应关系。 
	注意：为重建检测粒子的三维位置，至少需从三个视角确定对应关系。</li>
	<li>点击“3D 位置”，基于标定结果获取检测粒子的三维坐标。</li>
	<li>点击“无显示序列处理”，对所有图像序列重复执行步骤 6.4 至 6.7。 
	注意：此过程将为每组图像生成一个“rt_is”文件，其中包含帧内检测粒子的汇总信息，采用制表符分隔的文本文件格式。</li>
	<li>右键单击“Run”目录，然后点击“跟踪参数”，以定义用于搜索跟踪候选粒子的球体半径参数（例如 dvxmin 和 dvymin，单位为 mm/帧）。</li>
	<li>点击“无显示跟踪”，为步骤 6.7 中重建的粒子定义粒子识别编号（ID）。 
	注意：该过程采用四帧预测-校正方案24，关联相邻帧序列以实现粒子跟踪。此步骤将为每组图像生成一个 ptv_is 文件，其中包含帧内检测粒子的跟踪信息；前两列分别显示粒子在前一帧和后一帧中的 ID。</li>
	<li>点击“显示轨迹”，以在每个相机视角中可视化粒子运动轨迹。</li>
</ol>
7. 后处理（可选）
注意：后处理的范围和类型取决于个人需求，因此可进行定制。此处以点基计算为例进行简要说明。
<ol>
	<li>通过拉格朗日坐标系获取数据（使用 Matlab）。
	<ol>
		<li>从 ptv_is 文件中提取每个粒子的三维位置及其对应的 ID。该步骤可实现图像序列中检测到的粒子之间的关联，从而重建粒子轨迹。</li>
		<li>根据给定的帧率，计算每条轨迹中粒子的速度和加速度。在本示例中，通过对位置信号应用移动三次样条低通滤波来计算粒子的速度和加速度34,37。</li>
		<li>构建一个结构体数组，其字段包含每条轨迹的三维位置、三维速度、三维加速度、时间戳以及轨迹 ID。在此数据格式中，struct 数组的长度表示轨迹的数量。</li>
	</ol>
	</li>
	<li>通过欧拉坐标系获取数据（使用 Matlab）。
	<ol>
		<li>利用每个粒子的时间戳，将结构体数组（步骤 7.1.3）转换为时间序列形式。这将生成与步骤 7.1.3 中类似的结构体数组结构，但此时结构体数组的长度表示帧数，在本实验中为 9,000。</li>
		<li>将每个时间帧的时序结构体数组插值到三维网格上，以获得欧拉坐标系下的瞬时速度场。在本示例中，使用 Matlab 中的 griddata 函数进行插值。</li>
	</ol>
	</li>
</ol>

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

作者声明不存在任何竞争性经济利益。

3D-PTV 在揭示多种湍流的复杂物理机制方面具有巨大潜力，例如近地面大气中的大尺度湍流运动25，室内空气分布26或主动脉拓扑结构中的脉动流28 以及其他众多因素。然而，充分理解其优势与局限性并积累相关经验，对于最大化发挥其潜力至关重要。为确定最佳设置（包括帧率、照明光源、G）而进行的试错性初步测试和 exhaustive iterations（详尽迭代）是必不可少的。&amp;#38;B 水平和图像滤波方法与重建一组（的拉格朗日路径的能力直接相关例如，示踪）颗粒。需要注意的是，关键的实验步骤（如本实验所示）是调整G&amp;#38;B 水平以及视场的照明情况（包括卤素点光源、放大透镜和水槽底部反射镜的组合）。
这些调整有助于将光散射优化至四个视图的检测范围内。在确定了高保真测量的实验参数后，应进行全面的调整和故障排查，以根据帧率、相机分辨率和检测区域体积的大小，计算出可准确追踪的粒子轨迹的最大数量。尽管提高帧率可以增加捕获的粒子数量，但需要注意的是，3D-PTV 中可追踪的粒子数量远低于PIV。3...

湍流射流在工程应用中普遍存在。对这类流动进行详细表征，在从大规模环境排放系统到电子微尺度器件的广泛实际问题中至关重要。由于射流流动对众多重要应用具有显著影响，相关研究已十分深入1-4。多种实验技术已被用于在不同雷诺数和边界条件下表征射流流动，包括热线风速仪4-8、激光多普勒测速仪（LDV）4,9-12以及粒子图像测速技术（PIV）12-16。近年来，已有少数研究采用三维粒子追踪测速（3D-PTV）来研究射流流动中湍流与非湍流的交界面17,18。3D-PTV 是一种能够从不同视角描述复杂湍流场的特殊技术，它利用多视角立体成像，在拉格朗日参考系下重建流场体积内的粒子运动轨迹。该技术最初由 Chang19 提出，随后由 Racca 和 Dewey20 进一步发展。此后，3D-PTV 的算法和实验装置已得到多项改进21-24。基于这些进展和前期工作，该系统已成功应用于研究多种流体现象，例如在 4 m × 2 m × 2 m 域内的大尺度流体运动25、室内气流场26、脉动流27以及主动脉血流28。
三维粒子追踪测速（3D-PTV）测量的工作原理包括数据采集系统的搭建、记录/预处理、标定、三维匹配、时间序列追踪以及后处理。精确的标定能够实现对粒子位置的精确定位。通过三个以上视角图像中检测到的粒子对应关系，可基于对极几何重建出粒子的三维位置。连续图像帧之间的关联形成时间序列追踪，从而确定粒子的运动轨迹 s(t)。优化3D-PTV系统对于最大限度提高多粒子可追踪性至关重要。
优化的第一步是获取合适的数据采集系统，包括高速相机、照明光源以及示踪粒子的特性。相机分辨率与 interrogation volume 的尺寸共同决定了像素大小，因此也决定了所需示踪粒子的尺寸，其尺寸应大于单个像素。通过以亮度的平方作为权重对粒子像素的位置取平均值，可实现对检测到的粒子质心的亚像素精度估计21。相机的帧率与雷诺数及粒子追踪的关联能力密切相关。更高的帧率有助于解析更快的流动或更多数量的粒子，因为当图像间粒子的平均位移超过粒子间的平均间距时，追踪将变得更加困难。
在图像采集过程中，快门速度、光圈和感光度是需要考虑的三个因素。快门速度应足够快，以尽量减少粒子周围的模糊，从而降低粒子质心位置的不确定性。相机光圈应调整至与测量区域的景深相匹配，以降低检测到该区域外粒子的概率。由于相机的最大感光度是固定的，因此随着帧率的增加，用于照亮粒子所需的光照强度也应相应增加。与粒子图像测速（PIV）不同，三维粒子跟踪测速（3D-PTV）并不严格要求复杂的光学设置和高功率激光器，只要光源能从示踪粒子充分散射到相机即可。连续发光的LED或卤素灯是具有成本效益的良好选择，且无需同步装置21。
在三维粒子追踪测速（3D-PTV）中，与其他光学流场测量技术类似，示踪粒子的速度被假定为局部瞬时流体速度29然而，这种情况仅适用于直径和惯性均为零的理想示踪粒子；实际的示踪粒子应具有足够大的尺寸，以便能被相机捕捉。有限尺寸粒子的示踪保真度可通过斯托克斯数（Stokes number）来确定。 St, 即 粒子的弛豫时间尺度与所关注的湍流结构时间尺度之比。通常情况下， St 应明显小于1。对于 St&#160; ≤0.1 的流动追踪误差低于 1%30详细讨论可参见 Mei 等.29-313D-PTV 实验推荐的示踪粒子尺寸取决于光源和相机的灵敏度。当使用卤素灯或LED灯作为照明光源时，通常采用相对较大的粒子（例如 50-200 µm32，而较小的颗粒（例如 1-50 µm33,34 可与高功率激光器一起使用（例如 80-100 瓦连续波激光器）。对于特定波长的光具有高反射率的粒子，例如在卤素光下带有银涂层的粒子，可将其标记放大为图像。示踪粒子的 seeding density（播种密度）是实现成功三维PTV测量的另一个重要参数。粒子数量过少会导致轨迹数量不足，而粒子密度过高则会引起对应关系建立和追踪过程中的歧义。对应关系建立中的歧义包括粒子重叠以及在定义的极线（epipolar line）上检测到多个候选点。在追踪过程中，由于粒子平均间距较短，高播种密度会导致追踪歧义的发生。
第二步是优化记录/预处理过程中的设置以提高图像质量。摄影参数设置，例如增益 &amp;#38; 黑电平（G&amp;#38;B）在优化图像质量方面起着重要作用。黑电平定义了图像最暗部分的亮度水平，而增益则放大图像的亮度。G 的微小变化&amp;#38;B 水平会显著影响可追溯性的可能性。事实上，较高的 G&amp;#38;B 可能使图像过亮，并最终损坏相机传感器。为说明这一点，以下阐述 G 的影响&amp;本文还研究了#38;B水平对流场重建的影响。在预处理步骤中，图像通过高通滤波器进行滤波，以增强来自粒子的光散射信号。调整像素尺寸和灰度级，以在 interrogation volume 内最大化粒子检测效果。
优化的第三步是立体成像的精确校准，该过程基于对极几何、相机参数（焦距、主点和畸变系数）以及折射率变化。此步骤对于最小化基准目标点的三维重建误差至关重要。对极几何利用相机与测量区域之间的相对距离以及目标图像的倾斜角度。沿相机视场穿过测量区域时的折射率变化可根据Mass提出的方法予以考虑。 等21在本实验中，使用具有规则分布靶点的三维阶梯状结构作为靶标。
在三维粒子追踪测速（3D-PTV）实验中，尽管仅需两幅图像即可确定粒子的三维位置，但通常会使用更多的相机以减少不确定性21。一种替代昂贵多高速相机配置的方案是视场分光器，该装置由 Hoyer 等35 提出用于 3D-PTV，最近由 Gulean 等28 应用于生物医学领域。视场分光器由一个金字塔形反射镜（以下简称主反射镜）和四个可调节反射镜（以下简称次级反射镜）组成。本研究采用一个四视场分光器和单台相机，模拟四台相机的立体成像效果。该系统用于表征管射流在射流起始位置下游约 14.5–18.5 倍管径处的中间流场，射流直径，dh = 1 cm，雷诺数 Re ≈ 7,000，分别从拉格朗日和欧拉参考系进行分析。

<table><tbody><tr><td>Mikrotron 4CXP 4 通道 CXP-6 CoaXPress</td><td>ImageOps</td><td>CAMMC4082</td><td>高速相机</td></tr><tr><td>Active Silicon FireBird CoaX 帧采集卡</td><td>ImageOps</td><td>FBD-4XCXP6</td><td>帧采集卡</td></tr><tr><td>100 &amp;mu;m 镀银中空陶瓷微球</td><td>Potters Industries LLC</td><td>AG-SL150-30-TRD</td><td>示踪粒子</td></tr><tr><td>StreamPix6</td><td>Upstate Technical Equipment CO.,INC</td><td>MISNOR-STP-6-S-CL</td><td>相机应用软件</td></tr><tr><td>四视角分光器</td><td>Photrack AG</td><td></td><td>定制部件，若使用单台相机进行三维粒子追踪测速（3D-PTV）则为必需</td></tr><tr><td>250 瓦卤素聚光灯</td><td>General Electrics</td><td>23719</td><td>光源</td></tr><tr><td>OpenPTV（软件）</td><td>OpenPTV（http://www.openptv.net）</td><td></td><td>开源粒子追踪软件（注：可为希望免于繁琐安装过程或缺乏足够计算资源的用户提供服务形式的使用支持）</td></tr></tbody></table>

用于湍流研究的三维粒子跟踪测速技术：以射流流动为例

Three-dimensional Particle Tracking Velocimetry for Turbulence Applications: Case of a Jet Flow

3D-PTV 是一种定量流动测量技术，旨在通过立体成像序列记录，追踪一组粒子在三维空间中的拉格朗日轨迹。本文描述并讨论了由高速相机与四视角分光器构成的 3D-PTV 系统拓扑结构的基本组件、特征、约束条件及优化建议。该技术应用于雷诺数 Re ≈ 7,000 的圆形射流的中间流场区域（5 &#60;x/d &#60;25）。在射流出口下游约十个直径处及距射流中心轴不同径向距离的位置，估算了拉格朗日流动特征和欧拉框架下的湍流参数。拉格朗日参数包括选定粒子的轨迹、速度和加速度，以及流线曲率，这些参数由弗勒内-塞雷方程（Frenet-Serret equation）求得。将射流中心轴周围、位于射流出口下游十个直径处横截面上的三维速度场和湍流场的估算结果与文献数据进行了对比，并在距射流中心不同径向位置获得了大尺度流向速度运动的功率谱。

装置的照片和示意图如图 1和图 2所示。标定靶、视场分光器上反射的基准标记以及三维标定重建结果如图 3所示。在流向 x、展向 y 和深度 z 方向上，识别出的标定靶点的均方根误差（RMS）分别为 7.3 µm、5.7 µm 和 141.7 µm。z 坐标方向上相对较高的 RMS 值是由于该方向上的标定靶点数量较少，且四个视角相对于 z 轴的夹角较之 x 和 y 坐标方向更小所致。在任一时刻，四个视角中检测到的粒子数量均在 103 量级。在所有检测到的粒子中，由于仅有位于交叠区域内的粒子被成功捕捉，因此成功实现三维重建的粒子数量减少至大约一半。视频 1 展示了四视角分光器拍摄的射流流动的高速视频样本。
中间场区域内四个代表性粒子轨迹的样本，这些轨迹位于 x/d 附近并穿过该区域<s...

本文介绍了一种基于高速相机和四视角分光器的三维粒子跟踪测速系统（3D-PTV）。该技术应用于雷诺数 Re ≈ 7,000 时，圆形管道下游约十个管径距离范围内的射流流动研究。

3D-PTV 是一种定量流动测量技术，旨在利用立体成像在三维空间中追踪一组粒子的拉格朗日轨迹。

本实验的目标是建立并优化用于湍流研究的粒子追踪测速系统，并利用该系统研究射流流动。该方法可通过流体的拉格朗日描述，帮助解答湍流领域中的关键科学问题。该技术的主要优势在于能够同时在时间和空间上追踪大量粒子。尽管该方法可用于研究湍流及拉格朗日参考系，但它也可应用于其他领域，例如运动捕捉或运动追踪。本实验的装置用于研究管射流的中间流场。该水箱内设有一个水槽，用于产生直径为一厘米的射流。该示意图展示了在添加三维粒子追踪测速设备之前的实验装置概览。为准备粒子追踪，需确定研究区域，即 interrogation volume（ interrogation 体积）。在本实验中，该区域为一个立方体，其各面与较大体积的对应面平行，且其最近点距离水槽末端 15 厘米。选定测量区域后，首先调整主反射镜的位置。该反射镜呈金字塔形，安装在实验水箱侧面的垂直支架上。沿支架上下滑动反射镜，使其与测量区域中心对齐，然后固定其位置。接下来，将相机安装在主镜后方。调整相机成像的中心，使其与主镜的中心重合。然后，前往计算机检查其与相机的连接情况。接下来设置相机与记录软件之间的接口。首先获取主镜及其周围区域的相机视图。使用软件中的感兴趣区域设置，调整相机视图的宽度和高度，使其恰好覆盖主镜范围。下一步是在水箱中放置一个校准靶标，以设置次级反射镜。该靶标经过定制设计，可覆盖整个测量区域。靶标上具有直径为1毫米的标定点用于校准。首先，将一个可调节高度的平台放入水箱中，用于支撑标定靶。将标定靶放置在平台上，并调整其方向使其正对相机。调节平台高度，使标定靶的中心与标定空间的中心对齐。继续完成四视图分光器的组装。调整主镜的位置，使其能够捕获完整的检测区域。次镜安装在主镜两侧的垂直安装柱上。将每个副镜安装在大致与主镜最近侧的相机视角对齐的位置，并固定副镜。最终的布局应相对于主镜呈几何对称。使用激光笔依次将光束照射到镜面上，以可视化成像光路。根据光路引导，对次镜进行最后调整，使视场包含整个校准目标。重复可视化成像光路的过程，并对剩余的每个次镜进行相应调整。调整完所有反射镜后，监测相机视图中的子图像。在计算机上观察图像的相机视图。在视图分光器处，移动其中一个次级反射镜以检查子图像的重叠情况。如果仅有一个视图随镜子移动而变化，则表明重叠可忽略不计。安装四向分光器后，关闭相机并转而调整照明。将光源从上方直接对准检测区域放置。在本实验中，将放大镜直接置于光源下方以增强光强。打开光源，目视检查整个检测区域的光照是否均匀分布。通过返回相机进行操作，开始优化实验装置。调节镜头焦距，直至通过主镜反射的光线在次镜的四个视图中均达到均匀聚焦。确保视场分光器所成图像对称，并显示检测体积。继续调节光圈值，以便能够同时捕捉校准靶标上最近点与最远点的清晰图像。相机仅会捕获检测体积内的示踪粒子。此时，切换至记录软件，设置所需的帧率。本实验使用

three-dimensional-particle-tracking-velocimetry-for-turbulence

Research

科研

JoVE Journal

Engineering

工程

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

本文介绍了一种用于多种流体流动的定量三维（3D）成像技术。基于光场成像领域的原理，我们从图像阵列中重建三维体积。我们的三维结果涵盖广泛，包括速度场和多相气泡尺寸分布。...

通过多相机光场成像确定三维流场

微粒图像测速技术（μPIV）用于可视化血流中示踪微粒的成对图像，通过互相关分析获得精确的速度分布。剪切速率、最大流速、速度分布形态及流量等参数均可由这些测量结果推导得出，这些参数均具有临床应用价值。...

微粒子图像测速法用于微小血流速度分布测量

本文介绍了一种利用高分辨率、高速粒子图像测速技术（PIV）研究边界附近瞬态流动的方法。PIV是一种非侵入式测量技术，通过优化图像与记录特性、激光片光源特性以及分析算法等多个参数约束，可应用于任何光学可访问的流动研究。...

近表面高速粒子图像测速技术

我们利用3D打印技术制备出具有各向异性的颗粒，其形状包括钩状、十字形、四联体和三联体，这些颗粒在湍流流体中的取向和旋转可通过多幅同步视频图像进行测量。...

测量湍流中3D打印颗粒取向与旋转速率的方法

本文介绍了一种利用改进的烟雾流场可视化技术观测三角翼上前缘非定常涡流的实验方案，并研究导致前缘涡破裂位置发生振荡的机制。...

通过流动可视化方法对三角翼上方流场结构的实验研究

本文所述技术提供了一种低成本且相对简单的方法，可同时测量低颗粒浓度流动中的颗粒运动学和湍流。湍流测量采用粒子图像测速技术（PIV），而颗粒运动学则通过高速相机在重叠视场内获取的图像进行计算。...

利用流动成像技术同步测量湍流与颗粒运动学

本文描述了一种能够获取光学不透明流体中或通过光学不透明几何结构的二维速度场的回声粒子图像测速（EPIV）系统，并报道了在管道流动中的验证性测量结果。...

回声粒子图像测速法

牺牲组分汽化法（VaSC）用于制造微血管结构。该方法利用牺牲性聚乳酸纤维，在激光微加工导引板提供的精确三维几何定位下，形成中空的微通道。...

利用牺牲组分汽化法制备三维微结构的过程

介绍了两种用于研究液体射流冲击高速运动表面的实验装置：气炮装置和旋转圆盘装置。这些装置用于确定将液体摩擦调节剂（LFM）施加到铁轨顶部以控制摩擦力的最佳方法。...

高速液流冲击运动表面的可视化

恢复系数是一个用于描述碰撞过程中动能损失的参数。本文中，在真空条件下建立了一种自由落体实验装置，可用于测定微米级颗粒在高撞击速度下的恢复系数参数。...

开发用于真空条件下测量恢复系数的实验装置

了解三维粒子跟踪测速技术用于湍流研究，通过高速相机与四视角分光器分析圆形管道喷流特性。该方法可追踪大量粒子时空轨迹，为流体力学提供拉格朗日视角的实验数据。

用于湍流研究的三维粒子跟踪测速技术：以射流流动为例

探索更多视频