CN116539916A - 一种管道流动传质特性的测量方法和测量*** - Google Patents
一种管道流动传质特性的测量方法和测量*** Download PDFInfo
- Publication number
- CN116539916A CN116539916A CN202310542171.6A CN202310542171A CN116539916A CN 116539916 A CN116539916 A CN 116539916A CN 202310542171 A CN202310542171 A CN 202310542171A CN 116539916 A CN116539916 A CN 116539916A
- Authority
- CN
- China
- Prior art keywords
- trace
- flow
- particle
- flow field
- particles
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000012546 transfer Methods 0.000 title claims abstract description 30
- 239000002245 particle Substances 0.000 claims abstract description 181
- 238000010586 diagram Methods 0.000 claims abstract description 58
- 238000000037 particle-tracking velocimetry Methods 0.000 claims abstract description 37
- 239000012530 fluid Substances 0.000 claims abstract description 30
- 238000002156 mixing Methods 0.000 claims abstract description 30
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims description 68
- 230000010349 pulsation Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 8
- 239000013598 vector Substances 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 102100020760 Ferritin heavy chain Human genes 0.000 description 4
- 101001002987 Homo sapiens Ferritin heavy chain Proteins 0.000 description 4
- 230000002035 prolonged effect Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000007850 fluorescent dye Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 239000004952 Polyamide Substances 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011010 flushing procedure Methods 0.000 description 2
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 229920002647 polyamide Polymers 0.000 description 2
- 238000001454 recorded image Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007730 finishing process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000001499 laser induced fluorescence spectroscopy Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
- G01P5/20—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using particles entrained by a fluid stream
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明提供了一种管道流动传质特性的测量方法和测量***,该测量方法包括如下步骤:利用粒子跟踪测速方法PTV获得示踪粒子轨迹信息,获得短迹线流场图;步骤S2,采用流程图编译算法,获得涵盖目标时间步长的全流程长迹线流场图;根据延长后的示踪粒子轨迹信息包含的时间和位置信息对示踪粒子进行分类;利用全流程长迹线流场图,计算每条迹线示踪粒子的停留时间PRT,根据轨迹分类里每一类轨迹路径的示踪粒子停留时间,可以计算出每一时刻示踪粒子PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。采用本发明的技术方案,能够测量出单一主流管道的传质特性,使得到的粒子停留时间信息和耗尽比等参数更加准确。
Description
技术领域
本发明属于流体流动测量技术领域,尤其涉及一种管道流动传质特性的测量方法和测量***。
背景技术
在工业过程中,流动混合是一个十分重要的过程,从简单的单相混合、到复杂的多相混合反应***,其中的反应速率、产率都高度依赖混合效果。混合不当的话,将浪费不可逆的加工条件,降低产品的质量,导致下游精加工工艺成本和废物处理成本的增加,还有可能管道因冷热不均匀造成管壁应力疲劳破裂导致事故发生。为了提高产出、降低成本和事故发生率,在反应过程中对流动传质特性的监控、测量流体的滞留时间、加快流体的输运和混合就显得尤为重要。使用可视化技术不仅可实时显示流场中各位置混合效果的差异,还可以直观反映混合过程随时间的变化规律。
依赖路径信息的动力学特征在广泛的流体尺度上发挥着重要作用,特别是在生物医学、环境和气象流动方面。粒子跟踪测速方法测量产生描述流体微团在空间和时间中的运动的迹线使其成为研究传输现象的理想方法。路径线包含了直接计算随体导数用于提取压力场或拉格朗日加速度)。长的迹线包含了示踪粒子经过感兴趣区域的整个流动过程情况,对于研究夹带和测量与路径相关的量(如粒子停留时间)是不可或缺的。
目前,传统的流动传质特性测量多采用激光诱导荧光技术(PLIF)实验方法,其原理是在流体中溶解特定分子结构的荧光染色剂,利用一定波长的激光照射流场测量区,被照射到的示踪剂激发出荧光信号,通过相机等设备接收到荧光信号后,从荧光的分布,可以探测样品粒子的种类,从荧光的强弱可以得知粒子的浓度以及温度,利用其空间分辨率可以测量粒子的空间浓度和温度分布。
PLIF方法适用于旁路射流与主流混合的情形,例如涡环的夹带、反应器及引擎燃烧的传质测量。然而管道流动为单一主流流动时,此时管道没有旁路,无法在旁路流动中加入荧光染色剂,故PLIF传质测量技术并不适用于单一主流管道流动。作为一个定量的可视化技术,PLIF测量***复杂,不适合对大尺度范围内的混合过程进行全局测量。此外还有激光器发射波长单一,不同位置需和不同波长的激光器相配合,维护费用高等缺点。
发明内容
针对以上技术问题,本发明公开了一种管道流动传质特性的测量方法和测量***,能够测量出单一主流管道的传质特性。
对此,本发明采用的技术方案为:
一种管道流动传质特性的测量方法,包括如下步骤:
步骤S1,利用粒子跟踪测速方法PTV(particle tracking velocimetry,粒子追踪测速)获得示踪粒子轨迹信息,获得短迹线流场图;
步骤S2,采用流程图编译算法对所述短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图;也就是使用迹线延伸算法将获得的原始拉格朗日数据进行扩展,获得延长后的示踪粒子轨迹信息;
步骤S3,根据延长后的示踪粒子轨迹信息包含的时间和位置信息对示踪粒子进行分类;
步骤S4,利用全流程长迹线流场图,根据式(1)计算每条迹线示踪粒子的停留时间PRT,
PRT=tn-t0 (1),
其中,tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间;
根据轨迹分类里每一类轨迹路径的示踪粒子停留时间,计算出每一时刻示踪粒子PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。
粒子滞留时间(PRT)是表征流动停滞与循环区域的重要指标,其定义为流体在测量区域内的滞留时间。采用上述计算方式,不需要空间或时间平均,可以瞬间测量域内被跟踪的每一个粒子的PRT。通过跟踪单个颗粒,可以映射到高PRT颗粒它们的来源,并用于量化它们受到复杂的流动条件和/或几何形状的相互作用。
由于实际流动状况是三维流动,通过粒子跟踪测速仪PTV跟踪粒子的过程中,粒子位置信息有可能突然丢失,不能完整地贯穿整个感兴趣的区域,因而二维测量跟踪粒子的平均时间步长太短,无法直接开展有效的拉格朗日流场分析。而采用本发明技术方案,通过流程图编译算法对PTV测量获得的短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图。
迹线延伸后得到的粒子路径将会覆盖整个感兴趣的区域或者整个测量时间,大幅增加了每一时刻的PRT信息。
作为本发明的进一步改进,根据粒子的运动轨迹以及源端和终端位置对粒子进行分类,将粒子轨迹分为六类:粒子通过狭窄后在射流中停留整个长度的射流路径、粒子通过狭窄后与剪切层相互作用的夹带路径、起源于再循环区并在记录过程中被冲向下游的冲刷路径、从下游进入测量区域改变方向并从同一方向再次流出的再循环路径、从下游进入测量区域并留在测量区域内的逆流路径、从记录开始就存在并保持在测量区域内的停滞路径。此分类只是一个具体的举例,可以按照自己的意愿任意分类。
作为本发明的进一步改进,步骤S1采用现有技术的粒子跟踪测速方法PTV获得粒子轨迹信息。
进一步优选的,步骤S1包括:
设计实验模型,在管路中安装实验模型,确定测量区域,且模型安装时测量流场部分对应充满水的透明矩形水箱一侧;将高速相机放置在实验段的一侧,激光片光源从水箱一侧的透明窗射入,照亮实验模型测量截面,同时高速相机拍摄测量区域时观测流场图像。向水中加入一定量的跟随运动的示踪粒子,将整个管道充满水,开动精密齿轮泵,使得水在管道内循环流动,从观察窗以设定采集频率和曝光时间同步测量视窗的粒子图像。最后对粒子图像通过数据采集分析***上的PTV分析软件对粒子进行检测、匹配和跟踪。使用高斯拟合对记录的图像进行峰值匹配,在每一帧中识别出数千个粒子,得到短迹线流场图。
进一步的,所述激光片光源包括532nm绿光固体激光器和PTV光学镜组。
进一步的,管路上设有可编程控制精密齿轮泵,采用此技术方案,管路中的流体可以满足正弦型脉冲流,流速满足函数u=um[1+λsin(2πft)],喉部平均速度um和振幅比λ=u0/um,f为脉动频率。
进一步的,在可编程控制精密齿轮泵产生脉冲流时,释放示踪粒子如60μm球形聚酰胺颗粒。
进一步的,该实验模型为轴对称狭窄模型。
进一步的,所述实验模型的轮廓线满足正弦曲线函数: 式中,x为沿管道方向的长度(-0.5≤x≤0.5),D为管道直径,d为喉部直径,L为狭窄总长度的一半(L=D)。
进一步的,所述实验模型放置在充满水的矩形光学水箱中,使得实验模型的管道的弯曲壁变得光学透明,以减少光学折射率的不同造成的扭曲管道壁。
进一步的,所述实验模型的材质为有机玻璃,保证模型的透明度。
作为本发明的进一步改进,步骤S2中,所述流程图编译算法包括:对PTV测量获得的示踪粒子轨迹信息,向前(ti-ti+1)、向后(ti-ti-1)双时间步流场图使用局部加权回归获得流场全流程的趋势函数,然后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长,得到延长后的各示踪粒子的路径线。
进一步的,使用高斯拟合对相机记录的图像进行峰值匹配,在每一帧中识别出数千个粒子。峰值匹配将检测到的粒子在允许的速度范围内逐帧配对,计算单个速度矢量。对每个矢量的空间相干性与其相邻矢量进行验证,并将有效的矢量连接起来形成轨迹。然后根据路径长度对粒子轨迹进行过滤,如果一条轨迹延伸到少于5张连续图像上,或者任何一个单独的步长与前一步长的差异超过20%,则丢弃该轨迹。将最小轨迹的长度由5个时间步长改为2个时间步长对有效轨迹数没有显著影响,并且短的轨迹容易引入噪声误差。增加可接受的轨迹长度反而会降低粒子密度,并对粒子扩展产生负面影响。
作为本发明的进一步改进,步骤S2还包括:使用跨越五帧的三阶Savitzky-Golay滤波器对延长后的各示踪粒子的路径线进行平滑处理。
作为本发明的进一步改进,所述的管道流动传质特性的测量方法还包括:
步骤S5,根据全流程长迹线流场图,根据公式(2)计算一个周期后测量区域的流体滞留百分比R(T),
R(T)=(nT/n0)×100% (2)
其中,n0为脉动起始时图谱中的示踪粒子个数,nT为脉动周期结束时n0个示踪粒子中仍留在测量区域内的示踪粒子个数;
然后根据公式(3)计算耗尽比,
η=(1-R(T))/T* (3)
其中,R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。
采用现有技术单独通过粒子跟踪测速仪得到的粒子轨迹的粒子位置信息有可能在任意时刻突然丢失,或者任意时刻在感兴趣区域中间区域突然出现使得粒子轨迹信息数量增加,粒子轨迹信息数量增加或者减少都将影响该时刻感兴趣区域内有效粒子个数。通过迹线延伸技术,能够将不完整的轨迹信息补充到整个拍摄周期或者覆盖整个感兴趣区域,得到的耗尽系数更加准确。
作为本发明的进一步改进,所述的管道流动传质特性的测量方法还包括:
步骤S6,利用延长后的示踪粒子轨迹信息,采用公式(4)计算混合系数mc,
mc=1-∑(Ah,i/Af) (4)
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。具体而言,采用单个示踪粒子的面积与第i个判定区域内的示踪粒子的个数得到所占面积Ah,i。
本发明还公开了一种管道流动传质特性的测量***,其包括:
短迹线流场图获取模块,用于获取粒子跟踪测速方法PTV得到的示踪粒子轨迹信息,获得短迹线流场图;
全流程长迹线流场图获取模块,用于采用流程图编译算法对所述短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图;
轨迹分类模块,用于根据延长后的示踪粒子轨迹信息包含的时间和位置信息对示踪粒子进行分类;
示踪粒子停留时间计算模块,用于利用全流程长迹线流场图,根据式(1)计算每条迹线示踪粒子的停留时间PRT,
PRT=tn-t0 (1),
其中,tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间;
根据轨迹分类里每一类轨迹路径的示踪粒子停留时间,可以计算出每一时刻示踪粒子PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。
作为本发明的进一步改进,所述流程图编译算法包括:对PTV测量获得的示踪粒子轨迹信息,向前(ti-ti+1)/向后(ti-ti-1)双时间步流场图使用局部加权回归获得流场全流程的趋势函数,然后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长,得到延长后的各示踪粒子的路径线。
作为本发明的进一步改进,对延长后的各示踪粒子的路径线使用跨越五帧的三阶Savitzky-Golay滤波器进行平滑处理。
作为本发明的进一步改进,所述的管道流动传质特性的测量***还包括:
耗尽比计算模块,根据全流程长迹线流场图,根据公式(2)计算一个周期后测量区域的流体滞留百分比R(T),
R(T)=(nT/n0)×100% (2)
其中,n0为脉动起始时图谱中的示踪粒子个数,nT为脉动周期结束时n0个示踪粒子中仍留在测量区域内的示踪粒子个数;
然后根据公式(3)计算耗尽比,
η=(1-R(T))/T* (3)
其中,R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。
作为本发明的进一步改进,所述的管道流动传质特性的测量***还包括:
混合系数计算模块,利用延长后的示踪粒子轨迹信息,采用公式(4)计算混合系数mc,
mc=1-∑(Ah,i/Af) (4)
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。
与现有技术相比,本发明的有益效果为:
采用本发明的技术方案,通过跟踪流动粒子得到相关的时间和位置信息,在不增加旁路的情况下,能够测量出单一主流管道的传质特性;使用迹线延伸方法得到的粒子路径将会覆盖整个感兴趣的区域或者整个测量时间,大幅增加了每一时刻的PRT信息,解决了现有技术采用PTV直接测量获得的流场图迹线时间步长较短,无法直接开展有效的拉格朗日流场分析的问题;通过迹线延伸技术,能够将不完整的轨迹信息补充到整个拍摄周期或者覆盖整个感兴趣区域,得到的耗尽比等参数更加准确。
附图说明
图1是本发明实施例的实验装置图。
图2是本发明实施例的迹线延伸原理图。
图3是本发明实施例的轨迹分类示意图。
图4是本发明实施例的粒子耗尽示意图;其中,(a)为初始时刻,(b)为结束时刻。
图5是本发明实施例的粒子混合示意图。
附图标记包括:
1-可编程控制精密齿轮泵,2-轴对称狭窄模型,3-光学水箱,4-激光器,5-PTV光学镜组,6-流场测量图窗,7-CMOS高速相机,8-信号发出***,9-数据采集分析***;
3-1射流路径,3-2夹带路径,3-3冲刷路径,3-4再循环路径,3-5逆流路径,3-6停滞路径。
具体实施方式
下面对本发明的较优的实施例作进一步的详细说明。
一种管道流动传质特性的测量方法,包括如下步骤:
(1)通过粒子跟踪测速方法PTV获得粒子轨迹信息
搭建实验平台,如图1所示,包括在管路上设置可编程控制精密齿轮泵1和光学水箱3,将轴对称狭窄模型2设置在管路的光学水箱3中,所述轴对称狭窄模型2光学水箱3内设有流场测量图窗6,将532nm绿光固体连续激光器4和PTV光学镜组5组成的片光源设置在流场测量图窗6的一侧,并朝着流场测量图窗6,将CMOS高速相机7位于流场测量图窗6的一侧,并拍摄该流场测量图窗6的图像。532nm绿光固体连续激光器4发出3mm直径的激光束通过PTV光学镜组的扩散成3mm厚度的片光源,能够满足低速PTV实验对光源的要求。数据采集***9向信号发生***8发出采集指令后,信号发生***8向CMOS高速相机7同步发出触发信号,在532nm绿光固体连续激光器4照亮下,CMOS高速相机7以设定采集频率和曝光时间同步流场测量图窗6的粒子图像,并使得示踪粒子在高速射流区域位移不超过8像素。最后粒子图像通过数据采集分析***上的PTV分析软件处理,得到流场的高频速度信息,即短迹线流场图。
本实施例中,轴对称狭窄模型采用有机玻璃加工保证了模型的透明度,且狭窄模型轮廓线满足正弦曲线函数:式中,x为沿管道方向的长度(-0.5≤x≤0.5),D为管道直径,d为喉部直径,L为狭窄总长度的一半。轴对称狭窄模型放置在充满水的光学水箱中,使得管道的弯曲壁变得光学透明,以减少光学折射率的不同造成的扭曲管道壁。
本实施例中设计满足正弦型脉冲流,流速满足函数u=um[1+λsin(2πft)],喉部平均速度um和振幅比λ=u0/um,f为脉动频率。启动可编程控制精密齿轮泵循环抽取工作流体,并释放60nm球形聚酰胺颗粒。
数据采集分析***上的PTV分析软件对粒子进行检测、匹配和跟踪。
(2)通过路径扩展方法获得全流程流场图
利用流场图编译算法(FMC)获得全流程流场图谱,进而直接使用流场图进行流体输运分析,得到掺混系数、耗尽比、流体滞留时间等流体输运特征量。
使用高斯拟合对记录的图像进行峰值匹配,在每一帧中识别出数千个粒子。峰值匹配将检测到的粒子在允许的速度范围内逐帧配对,计算单个速度矢量。对每个矢量的空间相干性与其相邻矢量进行验证,并将有效的矢量连接起来形成轨迹。然后根据路径长度对粒子轨迹进行过滤,如果一条轨迹延伸到少于5张连续图像上,或者任何一个单独的步长与前一步长的差异超过20%,则丢弃该轨迹。将最小轨迹的长度由5个时间步长改为2个时间步长对有效轨迹数没有显著影响,并且短的轨迹容易引入噪声误差。增加可接受的轨迹长度反而会降低粒子密度,并对粒子扩展产生负面影响。
该方法针对PTV直接测量获得的流场图迹线时间步长较短,无法直接开展有效的拉格朗日流场分析的问题,通过流程图编译算法对PTV测量获得的短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图。其核心算法为流场图编译算法,即对PTV测量获得的向前(ti-ti+1)/向后(ti-ti-1)双时间步流场图使用局部加权回归(Lowess)获得流场全流程的趋势函数,而后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长。如图2所示给出了ti时刻的正向和反向时间拟合示意图。
鉴于示踪粒子浓度太低时来构建拟合函数会产生较大的嘈杂结果,而示踪粒子浓度太高时将过度平滑实际流动物理过程,可以使用跨越五帧的三阶Savitzky-Golay滤波器对生成的路径线进行平滑处理。
(3)轨迹分类
如图3所示,根据粒子的运动轨迹以及源端和终端位置对粒子进行分类,将粒子轨迹分为六类:粒子通过狭窄后在射流中停留整个长度的射流路径3-1、粒子通过狭窄后与剪切层相互作用的夹带路径3-2、起源于再循环区并在记录过程中被冲向下游的冲刷路径3-3、从下游进入测量区域改变方向并从同一方向再次流出的再循环路径3-4、从下游进入测量区域并留在测量区域内的逆流路径3-5、从记录开始就存在并保持在测量区域内的停滞路径3-6。
(4)利用全流程流场图计算粒子停留时间
粒子滞留时间(PRT)是表征流动停滞与循环区域的重要指标,其定义为流体在测量区域内的滞留时间。全流程图谱中每条迹线的滞留时间可简单计算为:
PRT=tn-t0 (1);
其中tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间。不需要空间或时间平均,可以瞬间测量域内被跟踪的每一个粒子的PRT。通过跟踪单个颗粒,可以映射到高PRT颗粒它们的来源,并用于量化它们受到复杂的流动条件和/或几何形状的相互作用。
计算上述轨迹分类里每一类轨迹路径粒子的粒子停留时间,可以计算出每一时刻粒子PRT的大小和位置分布,并且能够知道该粒子的来源,能够间接地反映流体混合效果。
由于实际流动状况是三维流动,通过粒子跟踪测速仪PTV跟踪粒子的过程中,粒子位置信息有可能突然丢失,不能完整地贯穿整个感兴趣的区域,因而二维测量跟踪粒子的平均时间步长太短,而使用迹线延伸方法得到的粒子路径将会覆盖整个感兴趣的区域或者整个测量时间,大幅增加了每一时刻的PRT信息。
(5)利用延长后的粒子轨迹信息计算耗尽比
如图4所示为粒子耗尽示意图。耗尽比为表征流体输运、混合的重要参量。其定义为单位时间内流出测量区域的流体质量与测量区域初始流体质量的比值,即
η=(1-R(T))/T* (3)
先计算R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。查询全流程图谱R(T)=(nT/n0)×100% (2)
其中,n0为脉动起始时图谱中的示踪粒子个数,nT为脉动周期结束时,n0个示踪粒子中仍留在测量区域内的示踪粒子个数。
单独通过粒子跟踪测速仪得到的粒子轨迹的粒子位置信息有可能在任意时刻突然丢失,或者任意时刻在感兴趣区域中间区域突然出现使得粒子轨迹信息数量增加,粒子轨迹信息数量增加或者减少都将影响该时刻感兴趣区域内有效粒子个数。通过迹线延伸技术,能够将不完整的轨迹信息补充到整个拍摄周期或者覆盖整个感兴趣区域,得到的耗尽系数更加准确。
(6)利用延长后的粒子轨迹信息计算混合系数
如图5所示为粒子混合示意图。掺混系数即测量区域内的流体混合的比率,计算公式为:
mc=1-∑(Ah,i/Af) (4),
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。Ah,i为采用单个示踪粒子的面积与第i个判定区域内的示踪粒子的个数的乘积得到。
本发明还公开了一种管道流动传质特性的测量***,其包括:
短迹线流场图获取模块,用于获取粒子跟踪测速方法PTV得到的示踪粒子轨迹信息,获得短迹线流场图;
全流程长迹线流场图获取模块,用于采用流程图编译算法对所述短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图;所述流程图编译算法包括:对PTV测量获得的示踪粒子轨迹信息,向前(ti-ti+1)/向后(ti-ti-1)双时间步流场图使用局部加权回归获得流场全流程的趋势函数,然后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长,得到延长后的各示踪粒子的路径线。并对延长后的各示踪粒子的路径线使用跨越五帧的三阶Savitzky-Golay滤波器进行平滑处理。
轨迹分类模块,用于根据延长后的示踪粒子轨迹信息包含的时间和位置信息对示踪粒子进行分类;
示踪粒子停留时间计算模块,用于利用全流程长迹线流场图,根据式(1)计算每条迹线示踪粒子的停留时间PRT,
PRT=tn-t0 (1),
其中,tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间;
根据轨迹分类里每一类轨迹路径的示踪粒子停留时间,可以计算出每一时刻示踪粒子PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。
耗尽比计算模块,根据全流程长迹线流场图,根据公式(2)计算一个周期后测量区域的流体滞留百分比R(T),
R(T)=(nT/n0)×100% (2)
其中,n0为脉动起始时图谱中的示踪粒子个数,nT为脉动周期结束时n0个示踪粒子中仍留在测量区域内的示踪粒子个数;
然后根据公式(3)计算耗尽比,
η=(1-R(T))/T* (3)
其中,R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。
混合系数计算模块,利用延长后的示踪粒子轨迹信息计算混合系数mc,
mc=1-∑(Ah,i/Af) (4)
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1.一种管道流动传质特性的测量方法,其特征在于:包括如下步骤:
步骤S1,利用粒子跟踪测速方法PTV获得示踪粒子轨迹信息,获得短迹线流场图;
步骤S2,采用流程图编译算法对所述短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图;
步骤S3,根据全流程长迹线流场图包含的时间和位置信息对示踪粒子进行分类;
步骤S4,利用全流程长迹线流场图,根据式(1)计算每条迹线示踪粒子的停留时间PRT,
PRT = tn - t0 (1),
其中,tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间;
根据轨迹分类里每一类轨迹路径的示踪粒子的停留时间,计算出每一时刻示踪粒子的PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。
2.根据权利要求1所述的管道流动传质特性的测量方法,其特征在于:步骤S2中,所述流程图编译算法包括:对PTV测量获得的示踪粒子的轨迹信息,向前(ti-ti+1)、向后(ti-ti-1)双时间步流场图使用局部加权回归获得流场全流程的趋势函数,然后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长,得到延长后的各示踪粒子的路径线。
3.根据权利要求2所述的管道流动传质特性的测量方法,其特征在于:步骤S2还包括:使用跨越五帧的三阶Savitzky-Golay滤波器对延长后的各示踪粒子的路径线进行平滑处理。
4.根据权利要求3所述的管道流动传质特性的测量方法,其特征在于,还包括:
步骤S5,根据全流程长迹线流场图,根据公式(2)计算一个周期后测量区域的流体滞留百分比R(T),
R(T)=(nT/n0)×100% (2)
其中,n0为脉动起始时图谱中的示踪粒子个数,nT为脉动周期结束时n0个示踪粒子中仍留在测量区域内的示踪粒子的个数;
然后根据公式(3)计算耗尽比,
η=(1-R(T))/T* (3)
其中,R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。
5.根据权利要求3所述的管道流动传质特性的测量方法,其特征在于,还包括:
步骤S6,利用全流程长迹线流场图包含的示踪粒子的轨迹信息,采用公式(4)计算混合系数mc,
mc=1-∑(Ah,i/Af) (4)
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。
6.一种管道流动传质特性的测量***,其特征在于:其包括:
短迹线流场图获取模块,用于获取粒子跟踪测速方法PTV得到的示踪粒子的轨迹信息,获得短迹线流场图;
全流程长迹线流场图获取模块,用于采用流程图编译算法对所述短迹线流场图进行编码,进而获得涵盖目标时间步长的全流程长迹线流场图;
轨迹分类模块,用于根据全流程长迹线流场图包含的时间和位置信息对示踪粒子进行分类;
示踪粒子停留时间计算模块,用于利用全流程长迹线流场图,根据式(1)计算每条迹线示踪粒子的停留时间PRT,
PRT = tn - t0 (1),
其中,tn为示踪粒子离开测量区域的时间,t0为进入测量区域的时间;
根据轨迹分类里每一类轨迹路径的示踪粒子的停留时间,可以计算出每一时刻示踪粒子PRT的大小和位置分布,获得该示踪粒子的来源信息,反映流体混合效果。
7.根据权利要求6所述的管道流动传质特性的测量***,其特征在于:所述流程图编译算法包括:对PTV测量获得的示踪粒子轨迹信息,向前(ti-ti+1)、向后(ti-ti-1)双时间步流场图使用局部加权回归获得流场全流程的趋势函数,然后使用趋势函数预测流场中迹线的运动轨迹,将短迹线延伸至需要的时间步长,得到延长后的各示踪粒子的路径线。
8.根据权利要求7所述的管道流动传质特性的测量***,其特征在于:对延长后的各示踪粒子的路径线使用跨越五帧的三阶Savitzky-Golay滤波器进行平滑处理。
9.根据权利要求7所述的管道流动传质特性的测量***,其特征在于,还包括:
耗尽比计算模块,根据全流程长迹线流场图,根据公式(2)计算一个周期后测量区域的流体滞留百分比R(T),
R(T) = (nT /n0 ) × 100% (2)
其中,n0为脉动起始时图谱中的示踪粒子的个数,nT为脉动周期结束时n0个示踪粒子中仍留在测量区域内的示踪粒子的个数;
然后根据公式(3)计算耗尽比,
η = (1-R(T)) /T* (3)
其中,R(T)为一个周期后测量区域的流体滞留百分比,T*为无量纲脉动周期。
10.根据权利要求7所述的管道流动传质特性的测量***,其特征在于,还包括:
混合系数计算模块,利用全流程长迹线流场图的示踪粒子轨迹信息,采用公式(4)计算混合系数mc,
mc=1-∑(Ah,i/Af) (4)
其中,Ah,i为第i个判定区域内流体来自底层示踪粒子所占面积,Af为判定区域面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310542171.6A CN116539916B (zh) | 2023-05-15 | 2023-05-15 | 一种管道流动传质特性的测量方法和测量*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310542171.6A CN116539916B (zh) | 2023-05-15 | 2023-05-15 | 一种管道流动传质特性的测量方法和测量*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116539916A true CN116539916A (zh) | 2023-08-04 |
CN116539916B CN116539916B (zh) | 2024-07-16 |
Family
ID=87452069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310542171.6A Active CN116539916B (zh) | 2023-05-15 | 2023-05-15 | 一种管道流动传质特性的测量方法和测量*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116539916B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221357A (ja) * | 1997-02-12 | 1998-08-21 | Ricoh Co Ltd | 流体計測装置および流体計測方法 |
US5905568A (en) * | 1997-12-15 | 1999-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stereo imaging velocimetry |
JP2005265486A (ja) * | 2004-03-16 | 2005-09-29 | Univ Nihon | 燃料の燃焼伝播解析及び火炎燃え広がりの測定装置、及び該装置を用いた燃料の火炎伝播モードの試験方法 |
JP2010218460A (ja) * | 2009-03-18 | 2010-09-30 | Tokyo Electric Power Co Inc:The | Ptv試験装置 |
US20110299738A1 (en) * | 2009-02-19 | 2011-12-08 | La Vision Gmbh | Method for Determining Flow Conditions |
US20140002617A1 (en) * | 2012-06-27 | 2014-01-02 | The Board Of Trustees Of The University Of Illinois | Particle tracking system and method |
CN203705468U (zh) * | 2013-12-26 | 2014-07-09 | 天津大学 | 基于红外摄像的新型表面流场测速*** |
CN105488760A (zh) * | 2015-12-08 | 2016-04-13 | 电子科技大学 | 基于流场的虚拟图像拼接方法 |
CN108760234A (zh) * | 2018-06-05 | 2018-11-06 | 哈尔滨工程大学 | 一种基于piv和ptv技术的流体流动与固体运动信息同步测试方法和装置 |
CN111721652A (zh) * | 2020-06-08 | 2020-09-29 | 重庆交通大学 | 一种基于piv和ptv技术的水力侵蚀细观机理观测装置及试验方法 |
US20210166404A1 (en) * | 2018-06-12 | 2021-06-03 | King Abdullah University Of Science And Technology | Single-camera particle tracking system and method |
CN114140490A (zh) * | 2021-10-21 | 2022-03-04 | 华南理工大学 | 一种运动粒子的二维及三维示踪轨迹续连方法 |
CN114612860A (zh) * | 2022-02-28 | 2022-06-10 | 北京交通大学 | 基于计算机视觉的轨道交通站内客流识别与预测方法 |
-
2023
- 2023-05-15 CN CN202310542171.6A patent/CN116539916B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221357A (ja) * | 1997-02-12 | 1998-08-21 | Ricoh Co Ltd | 流体計測装置および流体計測方法 |
US5905568A (en) * | 1997-12-15 | 1999-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stereo imaging velocimetry |
JP2005265486A (ja) * | 2004-03-16 | 2005-09-29 | Univ Nihon | 燃料の燃焼伝播解析及び火炎燃え広がりの測定装置、及び該装置を用いた燃料の火炎伝播モードの試験方法 |
US20110299738A1 (en) * | 2009-02-19 | 2011-12-08 | La Vision Gmbh | Method for Determining Flow Conditions |
JP2010218460A (ja) * | 2009-03-18 | 2010-09-30 | Tokyo Electric Power Co Inc:The | Ptv試験装置 |
US20140002617A1 (en) * | 2012-06-27 | 2014-01-02 | The Board Of Trustees Of The University Of Illinois | Particle tracking system and method |
CN203705468U (zh) * | 2013-12-26 | 2014-07-09 | 天津大学 | 基于红外摄像的新型表面流场测速*** |
CN105488760A (zh) * | 2015-12-08 | 2016-04-13 | 电子科技大学 | 基于流场的虚拟图像拼接方法 |
CN108760234A (zh) * | 2018-06-05 | 2018-11-06 | 哈尔滨工程大学 | 一种基于piv和ptv技术的流体流动与固体运动信息同步测试方法和装置 |
US20210166404A1 (en) * | 2018-06-12 | 2021-06-03 | King Abdullah University Of Science And Technology | Single-camera particle tracking system and method |
CN111721652A (zh) * | 2020-06-08 | 2020-09-29 | 重庆交通大学 | 一种基于piv和ptv技术的水力侵蚀细观机理观测装置及试验方法 |
CN114140490A (zh) * | 2021-10-21 | 2022-03-04 | 华南理工大学 | 一种运动粒子的二维及三维示踪轨迹续连方法 |
CN114972409A (zh) * | 2021-10-21 | 2022-08-30 | 华南理工大学 | 一种运动粒子的二维及三维示踪轨迹续连方法 |
CN114612860A (zh) * | 2022-02-28 | 2022-06-10 | 北京交通大学 | 基于计算机视觉的轨道交通站内客流识别与预测方法 |
Non-Patent Citations (3)
Title |
---|
DENIS KLIMENTJEW: "3D hypergraph-oriented air flow analysis based on PTV", THE 2010 IEEE INTERNATIONAL CONFERENCE ON INFORMATION AND AUTOMATION, 23 June 2010 (2010-06-23) * |
LIU, C: "Two-Dimensional Digital Particle Tracking Velocimetry Algorithm Based on the Image of Particle Trace", JOURNAL OF COASTAL RESEARCH, 31 December 2007 (2007-12-31) * |
刘文倩: "流道形状和供水压力对曝气滴灌迷宫流道内气泡运动特征的影响", 中国优秀硕士学位论文全文数据库 (农业科技辑), 15 February 2023 (2023-02-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN116539916B (zh) | 2024-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Westerweel et al. | The turbulent/non-turbulent interface at the outer boundary of a self-similar turbulent jet | |
Law et al. | Measurement of mixing processes with combined digital particle image velocimetry and planar laser induced fluorescence | |
Wang et al. | Second-order integral model for a round turbulent buoyant jet | |
Cherdantsev et al. | Simultaneous application of two laser-induced fluorescence approaches for film thickness measurements in annular gas-liquid flows | |
Mistry et al. | Entrainment at multi-scales across the turbulent/non-turbulent interface in an axisymmetric jet | |
Koochesfahani et al. | Simultaneous whole-field measurements of velocity and concentration fields using a combination of MTV and LIF | |
Alekseenko et al. | Investigation of waves interaction in annular gas–liquid flow using high-speed fluorescent visualization technique | |
Costantini et al. | Feasibility of skin-friction field measurements in a transonic wind tunnel using a global luminescent oil film | |
Cheng et al. | Time Resolved Scanning PIV measurements at fine scales in a turbulent jet | |
CN116539916B (zh) | 一种管道流动传质特性的测量方法和测量*** | |
RU121082U1 (ru) | Устройство бесконтактной оптико-лазерной диагностики нестационарных режимов вихревых течений | |
Ramaprabhu et al. | Simultaneous measurements of velocity and density in buoyancy-driven mixing | |
Di Cristo | Particle imaging velocimetry and its applications in hydraulics: a state-of-the-art review | |
Knopp et al. | Experimental and numerical investigation of turbulent boundary layers with strong pressure gradients | |
Walker et al. | Measurement of gas transfer across wind-forced wavy air–water interfaces using laser-induced fluorescence | |
Kräuter et al. | High resolution 2-D fluorescence imaging of the mass boundary layer thickness at free water surfaces | |
Ferrari | Image analysis techniques for the study of turbulent flows | |
Song et al. | An optimization scheme for the measurement of liquid jet parameters with rainbow refractometry based on Debye theory | |
Saga et al. | A comparative study of the PIV and LDV measurements on a self-induced sloshing flow | |
CN117214050B (zh) | 一种液固两相流速度分布的计量装置及计量方法 | |
Petrov et al. | High resolution experiments of velocity and concentration fluctuations in a jet flow | |
CN104535462A (zh) | 一种原位实时测量悬移质浓度和级配的装置及方法 | |
Seitz et al. | Thermal Characterization of a Turbulent Free Jet With Planar Laser Induced Fluorescence (PLIF) | |
Bresson et al. | Simultaneous temperature, concentration, and velocity fields measurements in a turbulent heated jet using combined laser-induced fluorescence and PIV | |
Hu et al. | Simultaneous velocity and concentration measurements in a turbulent jet flow by using PIV-PLIF combined system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 518000 Taoyuan Street, Nanshan District, Shenzhen City, Guangdong Province Applicant after: Harbin Institute of Technology,Shenzhen(Shenzhen Institute of science and technology innovation Harbin Institute of Technology) Address before: 518000 Taoyuan Street, Nanshan District, Shenzhen City, Guangdong Province Applicant before: HARBIN INSTITUTE OF TECHNOLOGY (SHENZHEN) |
|
GR01 | Patent grant |