Name: 测量的方向和旋转速率方法3D打印在湍流颗粒 (视频) | JoVE Journal
Uploaded: 2016-06-24T13:00:00.000+00:00
Duration: 12 min 34 s
Description: 了解利用3D打印技术制造各向异性颗粒并在湍流中测量其运动的方法。该技术通过立体视频成像和相机校准，实现复杂形状颗粒方向和旋转速率的精确测量，为流体力学研究提供新工具。

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

感谢Susantha Wijesinghe设计并构建了我们所使用的图像压缩系统。本研究得到了美国国家科学基金会（NSF）DMR-1208990号资助项目的支持。

1. 颗粒的制备
<ol>
	<li>使用三维计算机辅助设计软件创建粒子模型。每种模型导出为一个独立文件，文件格式需兼容所使用的3D打印机。
	<ol>
		<li>使用“圆”命令绘制直径为0.3 mm的圆，然后使用“拉伸”功能生成长度为3 mm的圆柱体。</li>
		<li>制作一个“十字”结构：两个正交的圆柱体共享同一中心；制作一个“三轴”结构：三个相互正交的圆柱体共享同一中心；制作一个“四面体”结构：四个圆柱体共享一个共同端点，彼此之间的夹角为109.5°；制作一个“平面三叉”结构：三个圆柱体位于同一平面，共享一个共同端点，彼此之间的夹角为120°。</li>
		<li>若需使圆柱体（此后称为粒子的“臂”）相对于彼此倾斜，可使用“三维旋转”命令，在圆柱一端的圆形截面直径上绘制一条直线，然后输入所需的旋转角度。</li>
		<li>使用“合并”命令将不同的臂连接成一个完整的密封整体。</li>
		<li>再次使用“三维旋转”命令倾斜整个物体，确保没有任何臂沿垂直或水平轴方向，因为沿这些轴方向的臂更容易出现缺陷、断裂或压扁。</li>
		<li>将每个物体分别导出为独立文件，文件格式需适用于3D打印机。</li>
	</ol>
	</li>
	<li>从专注于增材制造的商业供应商处订购每种类型约10,000个粒子，或在可用的设施中自行打印。粒子应使用聚合物挤出式3D打印机打印，该打印机需采用可溶解的不同材料作为支撑基质。
	<ol>
		<li>由于大量粒子的排布与打印过程较为耗时，应在计划实验前至少三周订购粒子。确保粒子以“高分辨率模式”打印，因为这些粒子接近许多3D打印机的最小特征尺寸，若以较低分辨率打印，臂的对称性会降低，且更易断裂。</li>
	</ol>
	</li>
</ol>
2. 颗粒的制备
<ol>
	<li>配制一种盐溶液，使颗粒在其中呈中性浮力，以尽量减少储存过程中颗粒臂的弯曲，并确保在分析时无需考虑重力和浮力的影响。
	<ol>
		<li>通过将颗粒浸入由水与氯化钙（CaCl₂）混合的溶液中，测试颗粒的平均密度。2) 在约 1.20 g/cm 的密度下3.
		<ol>
			<li>为测定水的密度，首先将空的100 ml容量瓶置于天平上并归零。取下容量瓶，加入含有CaCl混合的水2将烧瓶重新放回天平上，并将所得质量除以 100 ml。 
			注意：由于 1 ml = 1 cm3，1 g/ml = 1 g/cm³3.</li>
			<li>在多种不同溶液密度下测试颗粒，密度范围从1.16 g/cm³开始3 至 1.25 g/cm3，约0.01 g/cm3 逐步增加密度。在每种密度下测试多个颗粒，因为并非所有颗粒的密度都相同：在相同溶液中，部分颗粒会下沉，部分呈中性浮力，部分则会上浮。</li>
		</ol>
		</li>
		<li>记录数小时后颗粒平均达到中性浮力时的密度。 
		注意：测得的密度可能与颗粒生产商提供的体相密度有显著差异。</li>
		<li>混合约 400 kg 的 CaCl2 直至溶液达到2.1.1–2.1.2中记录的密度，总体积约为1,600 L。</li>
		<li>每种粒子类型（jacks、tetrads）各取出约 1 L 该混合溶液 等等。用于颗粒储存。将每升溶液分别置于不同的容器中，室温保存。其余溶液存放在大型储液罐中，室温保存。</li>
	</ol>
	</li>
</ol>
<img alt="微塑料显微图像；颗粒分析；尺寸标尺5 mm；环境研究。" src="/files/ftp_upload/53599/53599fig1.jpg" data-altupdatedat="1787630408000" data-alt="图1"> 
图1. 树脂去除过程中的不同阶段的样品。 a) 颗粒所附着的支撑树脂块。 b) 从其余部分分离出的一个单独树脂块。 c-e) 通过手工完成的多个树脂去除阶段。 f) 经NaOH浴处理并用罗丹明B染色后的单个样品。 <a href="https://www.jove.com/files/ftp_upload/53599/53599fig1large.jpg" target="_blank">请点击此处查看该图的高清版本。</a>
<ol start="2">
	<li>通过轻轻敲击将颗粒所包裹的支撑材料中的大块（约 5 mm × 320 mm，如图1a所示部分）破碎成小段（约 5 mm × 5 mm，图1b），手动松解颗粒周围的支撑材料，然后逐段用手揉搓，直至大部分多余的树脂脱落（图1c-e）。通过此方式去除多余树脂，以减少步骤 2.2.1 - 2.2.4 所需配制的 NaOH 溶液的用量。
	<ol>
		<li>将剩余的树脂块置于质量分数为 10% 的氢氧化钠（NaOH）溶液中，并放入超声波清洗仪中处理一小时。树脂与颗粒材质不同，因此 NaOH 可去除树脂而不会永久影响颗粒。 
		注意：该溶液具有腐蚀性，在超声波清洗过程中会发热。</li>
		<li>过滤分离颗粒。
		<ol>
			<li>为过滤颗粒，使用带有 0.1016 cm × 0.13462 cm 塑料孔的网布制作漏斗。将漏斗置于用于废弃 NaOH 溶液的容器上方，缓慢倾倒溶液使其通过滤网。按照环境健康与安全指南处置 NaOH 溶液。</li>
		</ol>
		</li>
		<li>用水轻轻冲洗颗粒，然后将其浸入新的质量分数为 10% 的 NaOH 溶液中，在超声波清洗仪中继续处理半小时。</li>
		<li>按照 2.2.2.1 的方法过滤分离颗粒，并将其储存在 2.1.4 中配制的密度匹配溶液中，直至颗粒硬化。操作时需小心，因为 NaOH 溶液会暂时软化颗粒。 
		注意：若颗粒未储存在密度匹配溶液中，部分臂状结构可能发生弯曲。将颗粒在密度匹配溶液中浸泡数小时，还可使塑料内部的部分孔隙被液体填充。</li>
	</ol>
	</li>
	<li>使用罗丹明-B（Rhodamine-B）与水的混合溶液对颗粒进行染色，使其在绿光激光照射下产生荧光。
	<ol>
		<li>配制 1 L 浓度为 0.5 g/L 的罗丹明-B 水溶液（后续简称为“染色液”）。 
		警告：有毒。</li>
		<li>将染色液加热至 50 至 80 °C，具体温度取决于颗粒材料。较硬的塑料使用较高温度；温度过高会导致臂状结构弯曲。</li>
		<li>取约 2,500 个颗粒（在密度匹配储存液中体积约为 25 ml，松散填充）放入染色液中，保持温度在 80 °C 持续 2 至 3 小时，使染料充分渗入聚合物。当颗粒变为粉红色（如图1f所示）时取出。 
		注意：加热会暂时软化颗粒。</li>
		<li>过滤分离颗粒，冲洗后再储存于 2.1.4 中配制的指定溶液中。颗粒会损失少量染料，导致溶液呈粉红色，但用自来水冲洗有助于防止染料过度流失。 
		注意：由于染色后颗粒的平均密度发生变化，需按照 2.1.1–2.1.2 重新测试，确定颗粒平均呈中性浮力时所需的新溶液密度。</li>
	</ol>
	</li>
	<li>根据需要调整 2.1.3 中的大量 CaCl2 溶液的密度。重复 2.1.4 步骤，并分离出新的密度匹配溶液。按照环境健康与安全法规处置旧的储存溶液（其中已含有少量罗丹明-B 染料）。</li>
	<li>对后续每组约 2,500 个颗粒重复 2.3.2–2.3.4 步骤，将相同形状的颗粒统一储存在 2.4 步骤中配制的相同密度匹配溶液中，并与不同形状的颗粒分开存放。 
	注意：大约重复 5 次 2.3.2–2.3.4 后，罗丹明-B 溶液的浓度将不足以有效染色。</li>
	<li>按照环境健康与安全法规处置 2.3.1 中配制的溶液，然后重复 2.3.1，重新配制 0.5 g/L 的新染色液用于后续染色。</li>
	<li>每完成 5 次 2.3.2–2.3.4 操作后，重复一次 2.6 步骤。</li>
</ol>
3. 实验与光学装置
<img alt="采用Nd:YAG激光器、反射镜和分束器进行三维物体重建的多视角摄影测量装置。" src="/files/ftp_upload/53599/53599fig2.jpg" data-altupdatedat="1787630408000" data-alt="图2"> 
图2. 实验装置示意图。 在振荡网格之间的八角形流动中，位于四个摄像机焦点位置的中心观测区域由一束绿色Nd:YAG激光照明。a) 侧视图，显示四个摄像机的布置及其与计算机的连接方式。图源自13。b) 俯视图，展示激光器、反射镜和透镜的配置，以实现中心区域的均匀照明。<a href="https://www.jove.com/files/ftp_upload/53599/53599fig2large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<ol>
	<li>准备相机。
	<ol>
		<li>使用分辨率至少为1兆像素、帧速率为450帧/秒的相机。</li>
		<li>布置相机，使每台相机均对准观测体积的中心并聚焦于该中心。可使用较少数量的相机，但粒子臂之间相互遮挡会限制取向测量的准确性，相机数量越少，实验越容易受到此效应的影响。使用四台以上的相机同样可能提高取向测量精度，因为这会降低所有相机上均出现臂被遮挡的可能性，而这是不确定性的主要来源。</li>
		<li>在满足装置约束条件的前提下，将相机以较大的角度（约90°）相互布置。按照图2所示方式放置相机，以平衡实验可操作性与各相机之间的夹角大小。通过在装置上设置垂直于每台相机观测方向的观察窗口，以尽量减少光学畸变。</li>
		<li>在每台相机上使用200 mm的微距镜头，以便在半米的工作距离下获得所需的测量体积。四台相机共同观测到的区域决定了检测体积，约为3 × 3 × 3 cm3。</li>
		<li>校准相机，以实现从测得的像素位置到三维空间坐标的转换。
		<ol>
			<li>将光圈设为f/11，并在相机上安装532 nm陷波滤光片，以滤除激光光，同时允许长波长的荧光通过并到达相机。</li>
			<li>将图像校准掩模置于水箱中，注入2.4步骤中的主体溶液，并照亮掩模。</li>
			<li>调整相机，使每台相机均能看见掩模，并且全部聚焦在掩模上的同一点。仔细对准相机，以优化检测体积的形状。</li>
			<li>从此时起，尽量减少对光学装置的任何改动。</li>
			<li>采集并存储每台相机拍摄的掩模图像。</li>
			<li>将溶液从水箱中排出，并泵回其先前储存的位置。</li>
			<li>从校准图像中提取表征每台相机位置、观测方向、放大倍数和光学畸变的参数。通过识别在四台相机上均可看见的校准掩模上的标记点，并定义这些点之间的距离来完成此过程。利用这些信息，采用标准校准方法提取相关参数。 
			注：基本校准方法见Tsai, 19877。本实验所用的具体实现方法见Oullette et al.3。希望开发相机校准软件的研究人员也可参考OpenPTV4。</li>
			<li>使用动态校准流程创建最终的校准文件。该步骤在获取示踪粒子数据后进行。采用非线性最小二乘法搜索，优化相机校准参数，使多个相机观测到的粒子位置之间的差异最小化。这些方法详见参考文献8和9。</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>使用平均功率可达50 W的Q开关绿色Nd:YAG激光器（以下简称“激光器”），照射水箱中心处流动均匀的区域，形成横截面直径约为3 cm的圆柱形照明区。8 
	注：激光器功率标定于5 kHz脉冲频率下。本实验中脉冲频率为900 Hz，此时输出功率显著降低。
	<ol>
		<li>利用分束器将激光束分成两束，并使用反射镜将其中一束导入水箱前部，另一束则垂直于第一束导入水箱侧面。</li>
		<li>在水箱外部、光束入射位置的对面放置两个额外的反射镜，以将光线反射回水箱内部，从而实现更均匀的照明，显著降低阴影效应。 
		注：反向传播光束产生的干涉效应的空间尺度极小，不会对本实验产生显著影响。</li>
	</ol>
	</li>
</ol>
4. 进行实验
<ol>
	<li>准备从每台相机录制视频。
	<ol>
		<li>编写一个实时去除不需要的图像数据的图像压缩系统。10,13
		<ol>
			<li>如果相机视野中没有颗粒，则不要保存图像。</li>
			<li>仅保存亮像素的位置和亮度，而非整幅图像。 
			注意：由于每个颗粒通常覆盖约5,000个亮像素，且视野中很少同时出现多个颗粒，因此图像压缩系统可显著减少使用高速相机记录数小时所需的存储空间。</li>
		</ol>
		</li>
		<li>准备数据采集软件。</li>
	</ol>
	</li>
	<li>制备1 x 1 x 1 m的湍流流动3 使用两个同步振荡的8 cm网格式网格的八角形水箱。8
	<ol>
		<li>泵入 CaCl2 将2.4中的溶液放入真空室中，并在室中静置过夜以脱气，从而减少实验中的气泡产生。</li>
		<li>通过真空室泵送溶液 0.2 µm 过滤至将进行实验的八角形水箱中。</li>
	</ol>
	</li>
	<li>运行实验。</li>
	<li>选择一种颗粒类型（示踪颗粒、十字形颗粒、叉形颗粒、四联体或三联体）用于第一轮实验，通过装置顶部的加样口将全部10,000个该类型颗粒加入水中。加样后关闭该加样口。
	<ol>
		<li>打开激光器。</li>
		<li>将相机和激光器设置为响应外部触发信号，并将触发频率分别设为相机450 Hz、激光器900 Hz。使用外部触发以确保所有相机同时开始采集，并在整个记录过程中保持同步。</li>
		<li>打开激光光阑。</li>
		<li>将网格设置为选定频率（1 或 3 Hz）并启动运行。在开始数据采集前，先运行网格约 1 分钟，以使湍流充分发展。</li>
		<li>记录 106 帧数以控制文件大小，并防止图像压缩系统中可能出现的错误影响过多数据。</li>
		<li>关闭激光光阑并停止相机触发。重置图像压缩系统和相机。
		<ol>
			<li>通过查看每个文件的部分内容，检查视频文件是否损坏。</li>
		</ol>
		</li>
		<li>重复 4.4.1 - 4.4.6 直至 107 已针对选定的颗粒在选定的网格频率下记录了图像。</li>
	</ol>
	</li>
	<li>将网格频率更改为在 4.4.4 中未选择的频率，并重复步骤 4.4.4 至 4.4.7</li>
	<li>排空水箱并过滤水体以去除所有颗粒。如有需要，可保留步骤 2.4 中储存水内的颗粒。</li>
	<li>对所有颗粒类型重复步骤 4.4 - 4.6。</li>
	<li>所有实验完成后，需再次校准相机，步骤同 3.1.5-3.1.5.7。</li>
</ol>
5. 数据分析
注意：本协议部分概述了获取颗粒取向和旋转速率所采用的流程。所使用的具体程序、测试图像和校准文件已作为本出版物的补充材料提供，并向所有感兴趣的读者开放使用。（详见补充文件“MATLAB_files.zip”中的“Use_Instructions.txt”文件。）
<ol>
	<li>利用相机标定参数，从多个相机拍摄的颗粒图像中获取其三维位置和取向。
	<ol>
		<li>在每一帧中，确定颗粒在四个图像上的中心位置。本实验中所有颗粒均具有足够的对称性，因此从任意视角观察时，颗粒的中心均位于图像上亮像素区域的几何中心。</li>
		<li>通过立体匹配该颗粒在四个相机上同时获取的二维位置，计算其三维位置3,8。</li>
		<li>建立一个可投影到各相机成像平面上的颗粒数值模型，以模拟该相机图像中的光强分布。
		<ol>
			<li>将颗粒建模为若干杆状结构的组合。利用3.1.5.7和3.1.5.8节中的相机标定参数，将每根杆的两个端点投影到各个相机上，然后根据软件协议，在二维空间中建模光强分布：沿杆的宽度方向采用高斯函数，沿长度方向采用费米-狄拉克函数。</li>
			<li>采用这种二维光强建模方式，以降低数据分析的计算成本。虽然对荧光颗粒建立完整的三维模型可进一步优化该方法，但会显著增加计算量。</li>
			<li>点击“运行”（Run）开始分析。</li>
		</ol>
		</li>
		<li>选择颗粒取向的初始估计值。
		<ol>
			<li>若正在分析该颗粒首次出现的帧，则初始估计值可设为一组随机的欧拉角。</li>
			<li>若该颗粒在至少前一帧中已出现，则以前一帧所确定的取向作为初始估计值。</li>
		</ol>
		</li>
		<li>进行非线性最小二乘拟合，以确定颗粒的取向。
		<ol>
			<li>优化三个三维位置坐标和三个欧拉角，使模型在四个相机上的二维投影与实测光强之间的平方差最小，具体操作依据软件协议。 
			注：定义欧拉角存在多种惯例。本文将角度 (φ, θ, ψ) 定义如下：φ 为绕 z 轴的初始旋转，生成新坐标轴 x' 和 y'；θ 为绕 x' 轴的旋转，生成新坐标轴 z' 和 y''；ψ 为绕新 z' 轴的旋转。11</li>
		</ol>
		</li>
		<li>选择相对于前一帧旋转量最小的取向。对于“jack”形颗粒，所得到的欧拉角对应24种对称取向之一；对于“tetrad”形颗粒，为12种对称取向之一；对于“cross”形颗粒，为8种对称取向之一；对于“triad”形颗粒，为6种对称取向之一。 
		注：5.1.6节中的方法假设颗粒在相邻帧之间的旋转不超过其内部角度的一半。该假设的合理性在“讨论”部分中予以说明。</li>
	</ol>
	</li>
	<li>将位置和欧拉角随时间的变化数据保存下来。</li>
	<li>利用这些数据提取刚体旋转速率及其他相关物理量。</li>
</ol>

<ol><li>Marcus, G., Parsa, S., Kramel, S., Ni, R., Voth, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Measurements+of+the+Solid-body+Rotation+of+Anisotropic+Particles+in+3D+Turbulence%5BTitle%5D)+AND+%22New+J.+Phys%22%5BJournal%5D)">Measurements of the Solid-body Rotation of Anisotropic Particles in 3D Turbulence.</a> New J. Phys. 16, 102001(2014).</li>
<li>Bretherton, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+motion+of+rigid+particles+in+a+shear+flow+at+low+Reynolds+number%5BTitle%5D)+AND+%22J.+Fluid+Mech%22%5BJournal%5D)">The motion of rigid particles in a shear flow at low Reynolds number.</a> J. Fluid Mech. 14 (02), 284-304 (1962).</li>
<li>Oullette, N., Xu, H., Bodenschatz, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+quantitative+study+of+three-dimensional+Lagrangian+particle+tracking+algorithms%5BTitle%5D)+AND+%22Exp.+in+Fluids%22%5BJournal%5D)">A quantitative study of three-dimensional Lagrangian particle tracking algorithms.</a> Exp. in Fluids. 40 (2), 301-313 (2006).</li>
<li>Open Source Particle Tracking Velocimetry. , Open PTV Consortium. Available from:
 <a target="_blank" href="www.openptv.net">www.openptv.net</a> (2014).</li>
<li>Parsa, S., Calzavarini, E., Toschi, F., Voth, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Rotation+Rate+of+Rods+in+Turbulent+Fluid%5BTitle%5D)+AND+%22Phys.+Rev.+Lett%22%5BJournal%5D)">Rotation Rate of Rods in Turbulent Fluid.</a> Phys. Rev. Lett. 109 (13), 134501(2012).</li>
<li>Parsa, S., Voth, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Inertial+Range+Scaling+in+Rotations+of+Long+Rods+in+Turbulence%5BTitle%5D)+AND+%22Phys.+Rev.+Lett%22%5BJournal%5D)">Inertial Range Scaling in Rotations of Long Rods in Turbulence.</a> Phys. Rev. Lett. 112 (2), 024501(2014).</li>
<li>Tsai, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+versatile+camera+calibration+technique+for+high-accuracy+3d+machine+vision+metrology+using+off-the-shelf+tv+cameras+and+lenses%5BTitle%5D)+AND+%22IEEE+Journal+of+Robotics+and+Automation%22%5BJournal%5D)">A versatile camera calibration technique for high-accuracy 3d machine vision metrology using off-the-shelf tv cameras and lenses.</a> IEEE Journal of Robotics and Automation. 3 (4), 323-344 (1987).</li>
<li>Blum, D., Kunwar, S., Johnson, J., Voth, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Effects+of+nonuniversal+large+scales+on+conditional+structure+functions+in+turbulence%5BTitle%5D)+AND+%22Phys.+Fluids%22%5BJournal%5D)">Effects of nonuniversal large scales on conditional structure functions in turbulence.</a> Phys. Fluids. 22 (1), 015107(2010).</li>
<li>Mann, J., Ott, S., Andersen, J. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Experimental+study+of+relative%2C+turbulent+diffusion%5BTitle%5D)+AND+%22RISO+Internal+Report%22%5BJournal%5D)">Experimental study of relative, turbulent diffusion.</a> RISO Internal Report. , R-1036 (1999).</li>
<li>Chan, K., Stich, D., Voth, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Real-time+image+compression+for+high-speed+particle+tracking%5BTitle%5D)+AND+%22Rev.+Sci.+Instrum%22%5BJournal%5D)">Real-time image compression for high-speed particle tracking.</a> Rev. Sci. Instrum. 78 (2), 023704(2007).</li>
<li>Goldstein, H., Poole, C., Safko, J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aGoldstein%2C+H+%22Classical+Mechanics%2C+3rd+Edition%22">Classical Mechanics, 3rd Edition</a>. , Addison-Wesley Publishing Company. 134-180 (2002).</li>
<li>Parsa, S. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22Rotational+dynamics+of+rod+particles+in+fluid+flows%22">Rotational dynamics of rod particles in fluid flows</a>. , Wesleyan University. Ph.D. Thesis (2013).</li>
<li>Wijesinghe, S. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22Measurement+of+the+effects+of+large+scale+anisotropy+on+the+small+scales+of+turbulence%22">Measurement of the effects of large scale anisotropy on the small scales of turbulence</a>. , Wesleyan University. Ph.D. Thesis (2012).</li>
<li>Wallace, J., Foss, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Measurement+of+Vorticity+in+Turbulent+Flows%5BTitle%5D)+AND+%22Annu.+Rev.+Fluid+Mech%22%5BJournal%5D)">The Measurement of Vorticity in Turbulent Flows.</a> Annu. Rev. Fluid Mech. 27, 469-514 (1995).</li>
<li>Su, L., Dahm, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Scalar+imaging+velocimetry+measurements+of+the+velocity+gradient+tensor+field+in+turbulent+flows.+I.+Assessment+of+errors%5BTitle%5D)+AND+%22Phys.+Fluids%22%5BJournal%5D)">Scalar imaging velocimetry measurements of the velocity gradient tensor field in turbulent flows. I. Assessment of errors.</a> Phys. Fluids. 8, 1869-1882 (1996).</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>Frish, M., Webb, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Direct+measurement+of+vorticity+by+optical+probe%5BTitle%5D)+AND+%22J.+Fluid+Mech%22%5BJournal%5D)">Direct measurement of vorticity by optical probe.</a> J. Fluid Mech. 107, 173-200 (1981).</li>
<li>Zimmerman, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Tracking+the+dynamics+of+translation+and+absolute+orientation+of+a+sphere+in+a+turbulent+flow%5BTitle%5D)+AND+%22Rev.+Sci.+Instrum%22%5BJournal%5D)">Tracking the dynamics of translation and absolute orientation of a sphere in a turbulent flow.</a> Rev. Sci. Instrum. 82 (3), 033906(2011).</li>
<li>Zimmerman, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Rotational+Intermittency+and+Turbulence+Induced+Lift+Experienced+by+Large+Particles+in+a+Turbulent+Flow%5BTitle%5D)+AND+%22Phys.+Rev.+Lett%22%5BJournal%5D)">Rotational Intermittency and Turbulence Induced Lift Experienced by Large Particles in a Turbulent Flow.</a> Phys. Rev. Lett. 106 (15), 154501(2011).</li>
<li>Klein, S., Gibert, M. athieu, Bérut, A., Bodenschatz, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Simultaneous+3D+measurement+of+the+translation+and+rotation+of+finite-size+particles+and+the+flow+field+in+a+fully+developed+turbulent+water+flow%5BTitle%5D)+AND+%22Meas.+Sci.+Technol%22%5BJournal%5D)">Simultaneous 3D measurement of the translation and rotation of finite-size particles and the flow field in a fully developed turbulent water flow.</a> Meas. Sci. Technol. 24 (2), 1-10 (2013).</li>
<li>Bellani, G., Byron, M., Collignon, A., Meyer, C., Variano, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Shape+effects+on+turbulent+modulation+by+large+nearly+neutrally+buoyant+particles%5BTitle%5D)+AND+%22J.+Fluid+Mech%22%5BJournal%5D)">Shape effects on turbulent modulation by large nearly neutrally buoyant particles.</a> J. Fluid Mech. 712, 41-60 (2012).</li>
</ol>

作者声明无竞争性财务利益。

在实验流体力学领域，长期以来人们一直认为测量湍流中粒子的涡量和旋转是重要的研究目标。小球体在湍流中的刚体旋转等于流体涡量的一半，但由于球体具有旋转对称性，直接测量其刚体旋转一直较为困难。传统上，流体涡量的测量依赖于复杂的多传感器热线探针14，但这类传感器仅能在具有较大平均速度的气流中获得单点涡量数据。此外，还发展了其他一些涡量测量方法。例如，Su 和 Dahm 采用基于标量图像的流场测速技术15，而 Lüthi、Tsinober 和 Kinzelbach 则使用了三维粒子追踪测速技术16。通过追踪单个粒子的旋转来测量涡量的方法最早由 Frish 和 Webb 开创，他们利用一种涡量光学探针测量了实心球形粒子的旋转17。该探针利用带有嵌入式平面晶体的小粒子，这些晶体充当镜子，产生随粒子旋转而改变方向的光束。近年来，已发展出通过成像标记在粒子表面的图案18,19或嵌入透明水凝胶粒子中的荧光粒子20来测量大尺寸球形粒子旋转运动的方法。为追踪各向异性粒子，Bellani 等. 采用了定制模具制备的...

在最近的一项研究中，我们提出使用由多个细长臂构成的粒子来测量湍流中粒子的旋转运动1。这些粒子可通过3D打印机制造，并且可以利用多台相机精确测量其位置、取向和旋转速率。借助细长体理论的工具可证明，这些粒子具有相应的有效椭球体2，其旋转运动与各自对应的有效椭球体完全相同。具有等长对称臂的粒子旋转方式与球体相同。其中一种粒子称为“jack”（ jack 粒子），其三个臂相互垂直并连接于中心。通过调节 jack 粒子各臂的相对长度，可形成等效于任意三轴椭球体的粒子。若将其中一个臂的长度设为零，则形成一个“cross”（十字形粒子），其等效椭球体为圆盘状。由细长臂构成的粒子所占实体体积远小于其对应的实心椭球体。因此，它们的沉降速度更慢，更易于实现密度匹配。这使得研究比使用实心椭球体粒子时大得多的粒子成为可能。此外，由于这些粒子遮挡其他粒子的光线比例较小，因此可在更高的粒子浓度下进行成像。
本文记录了3D打印颗粒的制备与追踪方法。多个研究团队已开发出通过多相机观测到的颗粒位置来追踪球形颗粒平动运动的工具3,4。Parsa et al.5 将该方法扩展至利用多相机观测到的杆状颗粒的位置和取向来实现对其运动的追踪。本文中，我们提出了制备多种形状颗粒并重建其三维取向的方法。这为将具有复杂形状颗粒的三维追踪技术拓展至广泛的新应用领域提供了可能。
由于可设计的颗粒形状种类繁多，该技术具有广阔的发展潜力。这些形状中的许多在环境流动中具有直接应用价值，例如浮游生物、种子和冰晶均呈现出极为多样的形态。颗粒旋转与湍流小尺度基本特性之间的关联6表明，对这些颗粒旋转行为的研究为探讨湍流级联过程提供了新的视角。

<table><tbody><tr><td>Condor Nd:YAG 50 W 激光器</td><td>Quantronics</td><td>532-30-M</td><td></td></tr><tr><td>高速相机</td><td>Basler</td><td>A504k</td><td></td></tr><tr><td>高速相机</td><td>Mikrotron</td><td>EoSens Mc1362</td><td></td></tr><tr><td>罗丹明-B</td><td>ScienceLab.com</td><td>SLR1465</td><td></td></tr><tr><td>氢氧化钠</td><td>Macron</td><td>7708</td><td>颗粒状。</td></tr><tr><td>500 Connex 3D 打印机</td><td>Objet</td><td></td><td>用于制造较小的颗粒。颗粒由 RP+M（快速成型加制造）公司提供。</td></tr><tr><td>VeroClear</td><td>Stratasys</td><td>RGD810</td><td>Objet 构建材料。</td></tr><tr><td>Clear Form 1 光敏树脂</td><td>Formlabs</td><td></td><td></td></tr><tr><td>柱面和球面透镜</td><td></td><td></td><td></td></tr><tr><td>200 mm 微距相机镜头</td><td></td><td></td><td>F 卡口。</td></tr><tr><td>超声波清洗机</td><td>Sonicator</td><td></td><td></td></tr><tr><td>氯化钙</td><td>Spectrum Chemical Mfg. Corp.</td><td>CAS 10043-52-2</td><td>颗粒状。</td></tr><tr><td>LabVIEW 系统设计软件</td><td>National Instruments</td><td></td><td>用于触发相机、控制网格以及触发激光器。</td></tr><tr><td>XCAP 软件</td><td>EPIX</td><td></td><td>与 LabVIEW 配合使用以触发相机。</td></tr><tr><td>MATLAB</td><td>Mathworks</td><td></td><td>用于所有图像和数据分析。从多幅图像中提取三维取向的程序随本出版物一并提供。</td></tr><tr><td>OpenPTV：开源粒子追踪测速法</td><td>OpenPTV 联盟</td><td></td><td></td></tr><tr><td>ParaView</td><td>Kitware</td><td></td><td></td></tr><tr><td>AutoCAD</td><td>AutoDesk</td><td></td><td>用于设计所有颗粒。所有颗粒设计的截图均来自 AutoCAD。</td></tr><tr><td>孔径为 0.040 × 0.053 英寸的筛网</td><td>Industrial Netting</td><td>XN5170&amp;ndash;43.5</td><td></td></tr><tr><td>相机滤光片</td><td>Schneider Optics</td><td>B+W 040M</td><td></td></tr></tbody></table>

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

Methods for Measuring the Orientation and Rotation Rate of 3D-printed Particles in Turbulence

本文介绍了测量各向异性粒子在湍流流体中旋转和平移运动的实验方法。采用3D打印技术制备具有细长臂并在共同中心连接的粒子，研究的形状包括十字形（两条相互垂直的杆）、星形（三条相互垂直的杆）、三联形（三条呈三角平面对称的杆）和四面体形（四条呈四面体对称的臂）。文中描述了制备约10,000个荧光染色粒子的方法。通过四台同步摄像机记录粒子在振荡网格间湍流（Rλ = 91）中的运动，获得其取向和刚体旋转速率的时间分辨测量结果。在这种相对低雷诺数的流动中，被携带的粒子足够小，可近似为椭球形示踪粒子。我们展示了粒子位置和取向的三维轨迹随时间变化的结果，以及对其旋转速率的测量数据。

图3a 显示了来自我们某台相机拍摄的一个四联体图像，其下方是根据该四联体轨迹某一段计算得到的欧拉角曲线图（图3c）。在图3b中，协议5 - 5.3中描述的取向识别算法的结果被叠加在四联体图像上。图3a中四联体的臂并不遵循用于构建模型的简单强度分布（协议5.1.3.1）。所有颗粒均表现出这一特性。此外，观测到的强度对臂之间的夹角、照明方向和观察方向具有非平凡的依赖关系12。尽管模型未包含上述任何因素，但仍能非常准确地测量颗粒的取向。
当通过最小二乘法拟合确定取向后，会保存粒子中心的三维坐标以及三个欧拉角（φ, θ, ψ），这三个角度用于定义其取向矩阵11。对于每个粒子在所有四台相机视野中的帧，均执行此操作。这些数据可用于重建粒子在整个观测体积内的完整轨迹，...

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

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

介绍了测量各向异性粒子在湍流流体中旋转和平移运动的实验方法。3D打印……

本实验的总体目标是利用3D打印机制造多种形状的颗粒，并测量其在湍流中的运动。该方法有助于回答流体力学中的关键问题，例如颗粒在湍流中如何旋转和取向。其主要优势在于，能够对具有复杂形状的颗粒进行高精度的时间分辨取向测量，这是以往方法所无法实现的。该方法为研究湍流中圆盘和球体的运动提供了深入见解，并可用于测量复杂形状颗粒在多种其他情况下的动力学行为。通常，初次使用该方法的研究人员在三维相机校准以及将取向识别算法适配于特定颗粒时会遇到困难。实验步骤将由 Guy Marcus 演示，他此前曾使用十字形和星形颗粒完成了原始实验。首先使用计算机辅助设计软件创建用于打印的颗粒模型。在 AutoCAD 中创建三联体时，首先使用圆命令绘制直径为 0.3 毫米的圆，然后使用拉伸命令将其拉伸成长度为 3 毫米的圆柱体。这种尺寸的圆柱体是构成颗粒的基本单元。现在使用旋转命令，将旋转轴设为沿X轴方向，旋转角度设为90度。调整视图以留出空间，便于操作其他圆柱体。接下来，复制该圆柱体以构成三联体的第二条臂。在该圆柱体末端，应用沿 Z 轴方向的三维旋转命令，将圆柱体旋转 120 度以形成三联体。将旋转后的圆柱体移至指定位置，以形成三联体的一条臂。圆柱体两端的中心应重合。以相同方式制作并定位第三条臂。然后，旋转物体，使圆柱体不沿垂直或水平轴方向。当所有臂部就位后，全部选中，并使用合并命令将它们连接成一个单一物体。将文件保存为适用于3D打印的格式。通过遵循类似的步骤制备实验中要使用的其他颗粒。除了三联体之外，还包括由两个正交圆柱体交叉组成且具有共同中心的十字形颗粒，由三个相互正交的圆柱体组成且具有共同中心的木偶形（jack）颗粒，以及四个圆柱体共用一个端点且彼此之间夹角为109.5度的四联体（tetrad）颗粒。每种类型的颗粒均需以高分辨率模式打印约10,000个。准备将打印出的微粒储存在一种使其呈中性浮力的溶液中。为此，使用氯化钙溶液，并加入约1,600升水，调节至中性浮力密度。每种微粒各取出约一升溶液用于储存。将一升溶液拿到实验台上，开始处理颗粒。本实验所用的颗粒为商业打印制备，到达时被包裹在用作支撑材料的树脂中。通过轻轻将大块材料破碎成小片来释放颗粒。对每一片段进行手动揉搓，直至大部分多余的树脂脱落。这将使后续步骤更加容易。为了进一步清洁颗粒，先在烧杯中配制质量分数为10%的氢氧化钠溶液。然后将部分颗粒及其残留的树脂包埋物放入烧杯中，以去除树脂。将盛有溶液的烧杯放入超声波清洗仪中。 首先将颗粒置于超声波清洗机中一小时，清洗后经冲洗再额外超声半小时。为回收颗粒，需准备第二个烧杯用于承接来自清洗机的氢氧化钠溶液，并准备好滤器。超声清洗步骤完成后，戴手套将烧杯从清洗机中取出，将内容物通过滤器倒入第二个烧杯中。在将颗粒转移至密度匹配溶液之前，先将其冲洗，以防止颗粒变形。这些颗粒已从树脂中分离出来，并在超声波浴中清洗干净，现可进行染色处理，使其在绿激光照射下发出荧光。染色时，使用一升罗丹明B（Rhodamine-B）水溶液，并将其置于加热板上。将染料加热至80摄氏度，以处理该颗粒材料。从储存容器中取少量颗粒松散地填充一个小容器，数量约为2,500个颗粒。将颗粒加入染料溶液中，并在浸泡期间保持温度。大约两到三个小时后，将溶液和颗粒从热源移开。使用滤网和另一个容器将颗粒从溶液中过滤分离。当颗粒因受热变软时，需小心操作，避免损坏。接着，仔细冲洗颗粒以去除多余的染料。完成后，将其储存于中性浮力溶液中，直至实验需要时再取出使用。实验将在一个装有摄像机的水槽中进行，摄像机用于记录流场中心的粒子。该示意图概述了水槽的设置，水槽具有八边形横截面，采用振荡网格产生湍流，以及测量区域和摄像机的位置。为确保方向测量的准确性，至少使用四台摄像机。每台相机的分辨率应至少达到100万像素，帧速率为450帧每秒，并连接至专用计算机。将相机以较大的夹角布置，确保它们均聚焦于目标测量区域。为尽量减少光学畸变，在水箱上安装了垂直于相机观测方向的观察窗。接下来，将图像校准掩模放置在水箱的测量区域内。这将有助于分析相机采集的实验数据。为模拟实验条件，开始向水箱中注入整体中性浮力溶液。当水箱注满后，设置光源以照亮遮罩。接着，返回操作相机。将每台相机对焦于遮罩上的同一位置，然后采集并存储用于校准的图像。排空水箱后，继续通过放置一块钕钇铝石榴石激光器来搭建实验装置。使用分光镜将部分激光束导向一面反射镜。光线从反射镜出发，经过透镜组形成柱状光束，随后进入水箱。在储液槽的另一侧放置一面镜子，将射出的激光反射回其原路径。利用分束器的另一输出端以及类似的镜片和透镜组合，产生一束与第一束光垂直的柱状光束，如顶部示意图所示。光束出射处的镜子有助于形成更均匀的照明。然后，检查装有脱气中性浮力氯化钙溶液的容器。在继续操作前，确保摄像头和数据采集软件已准备就绪。准备好将约10,000个已为本实验准备好的四分体加入容器中。将它们移至水箱顶部端口附近的位置。准备就绪后，打开端口，将颗粒加入水箱的水中。待所有颗粒加入后，关闭端口。接下来，前往激光器处。打开激光器并开启激光光阑。然后使用网格控制装置，使网格以三赫兹的频率开始振荡。在水槽中，让网格振荡约一分钟，以使湍流充分发展。此时，开始数据采集，并记录约一百万帧图像。运行结束后，停止数据采集并关闭激光光阑。所有实验完成后，移至储液罐出口处。在此处放置滤膜以从溶液中去除颗粒，并排空储液罐。将收集到的颗粒保存在中性浮力溶液中，以备后续使用。使用软件对记录的图像进行数据分析，以确定颗粒的位置和取向。以下步骤在软件中实现：从四个相机同步采集的颗粒图像开始，利用二维信息确定其三维位置。利用相机标定参数，将被照明颗粒的三维数值模型投影到每台相机上。这会在所有四台相机上生成颗粒光强度的二维模型。接下来，运行非线性最小二乘拟合程序，以确定最符合四台相机所观测图像的模型取向。结合来自四个摄像头的所有信息，以最小化投影到每个摄像头上的模型与观测到的颗粒之间的差异。这些步骤将应用于所有在四个摄像头视野内均可见颗粒的帧。这是四联体在位置和方向上的演化过程。注意四联体臂方向的不连续变化，若无保护机制则可能发生这种情况。该四联体演化的第二个版本表明，软件可确保每一帧的取向角相对于前一帧为最小旋转角度。位置和取向信息将被保存以供后续使用。这是实验中四联体的一个相机图像。在第二幅图像中，叠加的模型表示算法所确定的取向。尽管该二维强度模型较为简单，该方法仍能对颗粒取向进行精确测量。在实验中，利用取向求解算法获取每个粒子在整个轨迹过程中倾角的变化情况。这是对四面体颗粒穿过观察体积的完整轨迹进行的重建。通过结合粒子在四个相机视野中的每一帧图像中的取向信息和三维坐标，实现该重建。轨迹的长度为229帧。请注意，粒子未按比例绘制。该数据还可用于确定测得的翻滚速率的概率密度函数。红色方块表示十字交叉的数据，蓝色圆圈表示千斤顶的数据。实线代表球体数值模拟的结果，与实验数据吻合得非常好。一旦掌握，若操作得当，该技术可在三天内完成。在实施该实验过程中，需确保颗粒在不同密度的溶液中停留数小时。此外，必须保证在整个实验过程中相机正确对准、固定牢固并保持同步。完成此操作后，可采用其他方法（如示踪粒子追踪）来测量3D打印粒子在周围流场的完整速度梯度张量作用下的响应。观看本视频后，您应能够充分理解如何利用各向异性粒子，结合四组立体视频图像，从湍流中提取速度梯度张量的信息。请务必注意，操作过程中涉及的罗丹明B染料和氢氧化钠具有潜在危险性，因此执行该实验时必须采取防护措施，例如佩戴手套和护目镜。

methods-for-measuring-orientation-rotation-rate-3d-printed-particles

Research

科研

JoVE Journal

Engineering

工程

测量的方向和旋转速率方法3D打印在湍流颗粒

本文介绍了一种基于高速相机和四视角分光器的三维粒子跟踪测速系统（3D-PTV）。该技术应用于雷诺数 Re ≈ 7,000...

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

介绍了两种不同的方法，用于表征从层流到湍流条件下，单个颗粒的起始运动随沉积床几何结构的变化特性。...

基于视觉的规则基底中颗粒起动的表征：从层流到湍流条件

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

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

本方案展示了利用反应性喷墨打印技术制备可自主运动、生物相容且环境友好的微搅拌器，可用于生物医学和环境领域的应用。...

基于丝素蛋白的自驱动微搅拌器的反应性喷墨打印及推进力分析

本方案介绍了一种分析3D打印笔释放物的方法。通过测量释放颗粒的粒子浓度和粒径分布来评估其特性。进一步采用透射电子显微镜（TEM）对释放的颗粒进行分析。利用电感耦合等离子体质谱法（ICP-MS）对打印耗材中的金属含量进行定量检测。...

3D打印——评估3D打印笔的颗粒物排放

3D打印笔，颗粒物排放，气溶胶，打印温度，热塑性材料，聚乳酸（PLA），丙烯腈-丁二烯-苯乙烯共聚物（ABS），室内空气质量，暴露评估，纳米颗粒

我们采用等离子体增强化学气相沉积法，在多种材料的纳米级颗粒上沉积了厚度从几纳米到数百纳米的薄膜。随后通过刻蚀去除核心材料，形成中空纳米壳层，其渗透性由壳层厚度控制。我们表征了这些涂层对小分子溶质的渗透性，并证明这些屏障能够实现核心物质在数天内的持续释放。...

等离子体聚合中空颗粒的包封与渗透特性

本文演示了在基于表面活性剂的剪切稀化粘弹性流体中测量球形颗粒终端沉降速度的实验步骤。制备了具有广泛流变特性的流体，并在无限流体以及平行壁面间的流体中，测量了多种颗粒尺寸的沉降速度。...

无约束与受限条件下基于表面活性剂的剪切稀化粘弹性流体中球形颗粒沉降速度的实验测量

声流控器件利用微流控通道内的超声波来操控、富集和分离悬浮的微米及纳米级微粒。本方案描述了一种支持体声驻波的此类器件的制备与操作方法，可在无需鞘流液辅助的情况下将颗粒聚焦于中心流线。...

用于无鞘聚焦颗粒的体声驻波声流控器件的制备与操作

采用平面频率混合磁检测技术研制了一种用于平面样品中磁性颗粒成像的扫描仪。该技术通过双频激励，记录颗粒非线性无滞回磁化产生的磁互调响应，可用于对薄层生物样品进行二维成像。...

用于平面样品中磁性颗粒成像的频率混合磁检测扫描仪

本文介绍了一种基于阻抗的溶液蒸发速率检测装置。该装置相较于传统的失重法具有明显优势：响应速度快、检测灵敏度高、样品需求量小、可进行多样品测量，并且易于拆卸以进行清洁和重复使用。...

一种基于高性能阻抗的蒸发速率检测平台

了解利用3D打印技术制造各向异性颗粒并在湍流中测量其运动的方法。该技术通过立体视频成像和相机校准，实现复杂形状颗粒方向和旋转速率的精确测量，为流体力学研究提供新工具。

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

探索更多视频