CN112985758B - 测量复杂气液两相流中涡体尺度的实验装置及方法 - Google Patents

测量复杂气液两相流中涡体尺度的实验装置及方法 Download PDF

Info

Publication number
CN112985758B
CN112985758B CN202110156276.9A CN202110156276A CN112985758B CN 112985758 B CN112985758 B CN 112985758B CN 202110156276 A CN202110156276 A CN 202110156276A CN 112985758 B CN112985758 B CN 112985758B
Authority
CN
China
Prior art keywords
probes
probe
base
cross
vortex
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.)
Active
Application number
CN202110156276.9A
Other languages
English (en)
Other versions
CN112985758A (zh
Inventor
白瑞迪
王航
刘雪纯
汤荣才
牛志攀
赵维扬
刘善均
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN202110156276.9A priority Critical patent/CN112985758B/zh
Publication of CN112985758A publication Critical patent/CN112985758A/zh
Application granted granted Critical
Publication of CN112985758B publication Critical patent/CN112985758B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明提供一种测量复杂气液两相流中涡体尺度的实验装置及方法。所述实验装置包括探针、基座、模拟信号‑数字信号模块、探针供电电源和数据处理端，所述基座为“十”字交叉型，基座上沿十字型的横部和竖部分别等间距排列有若干探针，基座内部为空腔结构，探针竖立固定在基座上，针头朝向空中，尾部穿过基座***基座的空腔中，每根探针尾部在基座空腔内均连接有一根导线，每根导线从基座底部引出基座再与模拟信号‑数字信号模块连接，所述模拟信号‑数字信号模块通过数据传输线与数据处理端连接；所述探针分别与探针供电电源连接。本发明装置和计算方法可计算三维空间涡体尺度的矢量值，能够真实准确的反映复杂气液两相流中涡体尺寸的流动规律。

Description

测量复杂气液两相流中涡体尺度的实验装置及方法
技术领域
本发明属于海洋工程、化工工程、环境工程以及水利水电工程中复杂的气液两相流领域，具体涉及气液两相流中流速分布和涡体尺度等微观特性。
背景技术
在复杂的气液两相流流动中，如海洋工程中波浪破碎、化工工程的射流、环境工程中的鼓泡塔、水利工程中的水跃等，大量的气泡在各种作用力下被卷入到液体中。宏观角度，一方面由于气体的卷入，液体由单一颜色变为“乳白”色；另一方面气体的卷入改变原有液体的密度，使得液面抬升等。从微观角度，由于气泡的卷入会改变原有流体的流动产生影响；另一方面卷入液体的气泡往往以单一气泡或气泡簇的形式存在，气泡的尺寸从毫米到厘米级别都存在，而不同尺寸的气泡大小对微观空间的紊动流速、涡体尺寸的影响更为复杂。
目前关于复杂气液两相流中的研究，一般采用电压探针或光学探针对微观气泡尺寸、掺气浓度等进行测量。目前较为广泛采用的一种前后两根探针布置方式，通过气泡通过两根探针的时间差异，获得空间某一点顺水流方向的流速、掺气浓度、气泡个数及尺寸等。它无法反映三维空间涡体尺度的大小与方向。但这种实验装置仅仅只能对水流方向的流动特性和气泡特性进行测量，而对复杂气液两相流中三维流速分布和涡体尺度特性等无法进行测量。
发明内容
本发明的目的是针对现有实验技术的不足，提供一种测量复杂气液两相流中涡体尺度的实验装置及方法，对复杂的三维流动的流速分布与涡体尺度等进行分析，为复杂气液两相流的研究提供技术支撑。
本发明提供的测量复杂气液两相流中涡体尺度的实验装置，包括探针、基座、模拟信号-数字信号模块、探针供电电源和数据处理端，所述基座为“十”字交叉型，基座上沿十字型的横部(X轴方向)和竖部(Y轴方向)分别等间距排列有若干探针，所述基座内部为空腔结构，所述探针竖立固定在基座上，探针针头朝向空中，探针尾部穿过基座***基座的空腔中，每根探针尾部在基座空腔内均连接有一根导线，每根导线从基座底部引出基座，再与模拟信号-数字信号模块连接，将探针采集的模拟信号转化为数字信号，所述模拟信号-数字信号模块通过数据传输线与数据处理端连接；所述探针分别与探针供电电源连接。
上述实验装置，进一步地，所述探针由不锈钢外壳包裹银丝内芯构成，探针针尖部的银丝内芯暴露在外，探针尾部的银丝与对应的导线相连接，探针的不锈钢外壳与探针供电电源连接。银丝直径为0.1mm，不锈钢外径的直径为0.7mm。为增加探针在基座上的稳定性，探针外形设置为圆锥体形。为保证结构的稳定，防止被水流冲坏，探针焊接在十字架底座上。
上述实验装置，进一步地，基座横部和竖部上的探针数量相等，分别为4～10 根。
上述实验装置，进一步地，所述探针相对基座的高度，从原点处的探针开，沿X轴正负方向递减，沿Y轴正负方向递减。优选地，X轴和Y轴上，相邻两根探针之间的高程差相等，相邻两根探针之间的水平间距相等。
上述实验装置，进一步地，由于装置工作时探针需要逆水流方向，因此基座的横部和竖部为与探针朝向方向一致的流线型。
上述实验装置，进一步地，所述模拟信号-数字信号模块作用是将模拟信号转变为数字信号，可选用现有市售具有相应功能的模块，如美国生产的NI采集器、成都泰斯特高速数据采集器等。这些采集器有多个信号接口，可同时对多个探针信号进行采集。
上述实验装置，进一步地，数据处理端为配备有相应的数据处理软件的电脑。所述处理软件可选用Signal Express。
上述实验装置，进一步地，所述探针供电电源为5V或10V直流电压输出的电源转换模块。
本发明提供的上述测量复杂气液两相流中涡体尺度的实验装置工作原理是：通过外部一10V直流电源与探针相接对探针进行供电，由于液体具有导电性，气体无导电性，根据探针针尖在气液两相流中液体与气体中电压信号的不同，得到相应的电压随时间的变化的模拟信号，在数据传输线传输到数据处理端，得到电压信号V随时间t的变化曲线，根据曲线再进行后续的分析和处理，计算出三维空间涡体尺度的矢量值，实现对复杂的三维流动的流速分布与涡体尺度的分析。
本发明提供基于上述实验装置的测量复杂气液两相流中涡体尺度的方法，包括以下步骤：
(1)数据采集
对探针进行编号，十字型基座的X轴的正负两个方向分别定义为Xm、Xn方向， Y轴的正负两个方向分别定为Yi、Yj方向；X轴中间(原点处)探针标记为1#，1# 探针一侧的探针逐次标记为2#、3#、……、m#，1#探针另一侧的探针逐次标记为 2#、3#、……、n#；在基座的Y轴方向上，中间探针(原点处)探标记为1#，1# 探针一侧的探针逐次标记为2#、3#、……、i#，1#探针另一侧的探针逐次标记为2#、 3#、……、j#；X轴方向上相邻探针的间距记为ΔX，Y轴方向上相邻探针的间距记为ΔY，优选ΔX＝ΔY；X轴方向上相邻探针的高差记为Δz，Y轴方向上相邻探针的高差记为Δz；
将整个装置放在水里中，并将探针的针尖迎着水流方向；运行电脑上的数据处理软件，如Signal Express软件，在数据处理软件里面设置测量的时间为t和采样的频率为F，得到电压信号V随时间t的变化曲线；
(2)对Xm轴方向探针进行信号处理
将1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RX12)max，以及对应的时间tX12；然后依次分别对探针1#与3#、1#与4#……1#与m#探针所测的电压信号进行互相关分析，并获得最大的相关系数值(RX1m)max与对应的时间tX1m；
根据每个互相关曲线最大相关系数对应的时间tX12、tX13、……、tX1m，求得相应的速度大小VX12、VX13、……、VX1m，由于两个探针之间有一定固定的距离和固定的角度(根据ΔX、ΔY与ΔZ的距离确定)，因此可获得速度的矢量值：
Figure BDA0002933614220000031
Figure BDA0002933614220000032
Figure BDA0002933614220000033
其中，关于ΔLX的计算公如下：
Figure BDA0002933614220000034
Figure BDA0002933614220000041
根据互相关分析方法，得到互相关曲线，分别对每个互相关曲线进行积分，积分的时间区域t>tX1m得到相应的相关时间尺度大小TX12、TX13、……、TX1m值：
Figure BDA0002933614220000042
Figure BDA0002933614220000043
Figure BDA0002933614220000044
同样，根据互相关分析方法，得到1#探针与X方向上其他探针的互相关曲线，通过对互相关图进行积分，可获得涡体时间尺度TXX12、TXX13、……、TXX1m值：
Figure BDA0002933614220000045
Figure BDA0002933614220000046
Figure BDA0002933614220000047
分别将两探针之间的速度矢量与与涡体时间尺度进行乘积，可获得相应的涡体尺度的矢量值
Figure BDA0002933614220000048
Figure BDA0002933614220000049
Figure BDA00029336142200000410
Figure BDA00029336142200000411
将获得的涡体尺度的矢量值进行矢量相加，即可得到Xm方向上的涡体尺度矢量值；
Figure BDA0002933614220000051
3)对Yi轴方向布置的探针进行信号处理
将Y轴方向上的1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RY12)max，以及对应的时间tY12；然后依次分别对探针1#与3#、1#与 Y4#……1#与i#等探针所测的电压信号进行互相关分析，并获得最大的相关系数值 (RY1i)max与对应的时间tY1i；
根据每个互相关曲线最大相关系数对应的时间tY12、tY13、……、tY1i，可求得相应的速度大小VY12、VY13、……、VY1i，由于两个针之间有一定固定的距离和相应的夹角，因此可获得速度的矢量值；
Figure BDA0002933614220000052
Figure BDA0002933614220000053
Figure BDA0002933614220000054
其中，关于ΔLY的计算公如下：
Figure BDA0002933614220000055
Figure BDA0002933614220000056
分别对每个互相关曲线进行积分，积分的时间区域t>tY1m，得到相应的时间尺度大小TY12、TY13、……、TY1i值：
Figure BDA0002933614220000057
Figure BDA0002933614220000058
Figure BDA0002933614220000059
根据1#探针与Y方向上其他探针的互相关图，通过积分求解，可获得涡体时间尺度TYY12、TYY13、……、TYY1m值：
Figure BDA0002933614220000061
Figure BDA0002933614220000062
Figure BDA0002933614220000063
分别将两探针之间的速度矢量与与时间尺度进行乘积，可获得相应的涡体尺度的矢量值
Figure BDA0002933614220000064
Figure BDA0002933614220000065
Figure BDA0002933614220000066
Figure BDA0002933614220000067
将获得的涡体尺度的矢量值进行矢量相加，即得到Yi方向上的涡体尺度矢量值；
Figure BDA0002933614220000068
(4)采用同样方法获得Xn与Yj方向上涡体尺度的矢量值
Figure BDA0002933614220000069
和
Figure BDA00029336142200000610
最后将 X和Y轴四个方向Xm、Xn与Yi、Yj的涡体尺度的矢量值进行相加，可获得该测点处的涡体尺度的三维矢量值:
Figure BDA00029336142200000611
与现有技术相比本发明具有以下有益效果：
现有的测量技术，是在一平面内布置两根探针，通过每次改变二者的间距大小，通过多次测量获得二维平面内涡体尺度的标量值(有大小无方向)。本发明实验装置与相应的算法，可计算三维空间涡体尺度的矢量值，能够真实准确的反映复杂气液两相流中涡体尺寸的流动规律。
附图说明
图1为本发明测量装置的三维示意图；
图2本本发明测量装置连接关系图；
图3为本发明装置X方向探针布置图(a)和设备Y方向探针布置图(b)；
图4为探针所采集电压信号随时间变化曲线；
图5为1#探针与X方向各探针信号的互相关图；
图6为1#探针与Y方向各探针信号的互相关图；
图7本发明方法的计算流程图。
图中，1—探针；2—基座；3—导线；4—数据传输线；5—模拟信号-数字信号模块；6—数据处理端；7—探针供电电源。
具体实施方式
下面通过具体实施方式对本发明做进一步说明。以下所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书内容所作的等效变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
实施例1
测量复杂气液两相流中涡体尺度的实验装置，包括探针1、基座2、模拟信号- 数字信号模块5、探针供电电源和数据处理端6，所述基座为“十”字交叉型，由于装置工作时探针需要逆水流方向，因此基座的横部和竖部为与探针朝向方向一致的流线型。基座上沿十字型的横部(X轴方向)和竖部(Y轴方向)分别等间距排列有若干探针，所述基座内部为空腔结构，所述探针竖立固定在基座上，探针针头朝向空中，探针尾部穿过基座***基座的空腔中，每根探针尾部在基座空腔内均连接有一根导线3，每根导线从基座底部引出基座，再与模拟信号-数字信号模块连接，将探针采集的模拟信号转化为数字信号，所述模拟信号-数字信号模块通过数据传输线与数据处理端连接；所述探针分别与探针供电电源连接。由于装置工作时探针需要逆水流方向，因此基座的横部和竖部为与探针朝向方向一致的流线型。
所述探针由不锈钢外壳包裹银丝内芯构成，探针针尖部的银丝内芯暴露在外，探针尾部的银丝与对应的导线相连接，探针的不锈钢外壳与探针供电电源连接。银丝直径为0.1mm，不锈钢外壳外径为0.7mm。为增加探针在基座上的稳定性，探针外形设置为圆锥体形。为保证结构的稳定，防止被水流冲坏，探针焊接在十字架底座上。
基座横部和竖部上的探针数量相等，分别为4根。所述探针相对基座的高度，从原点处的探针开，沿X轴正负方向递减，沿Y轴正负方向递减。X轴和Y轴上，相邻两根探针之间的高程差相等，相邻两根探针之间的水平间距相等。
所述模拟信号-数字信号模块作用是将模拟信号转变为数字信号，选用美国生产的NI采集器。采集器有多个信号接口，可同时对多个探针信号进行采集。数据处理端为配备有相应的数据处理软件的电脑。所述处理软件可选用Signal Express。
实施例2
基于实施例1实验装置的测量复杂气液两相流中涡体尺度的方法，包括以下步骤：
(1)数据采集
对探针进行编号，十字型基座的X轴的正负两个方向分别定义为Xm、Xn方向， Y轴的正负两个方向分别定为Yi、Yj方向；X轴中间(原点处)探针标记为1#，1# 探针一侧的探针逐次标记为2#、3#、……、m#，1#探针另一侧的探针逐次标记为 2#、3#、……、n#；在基座的Y轴方向上，中间探针(原点处)探标记为1#，1# 探针一侧的探针逐次标记为2#、3#、……、i#，1#探针另一侧的探针逐次标记为2#、 3#、……、j#；X轴方向上相邻探针的间距记为ΔX，Y轴方向上相邻探针的间距记为ΔY，优选ΔX＝ΔY；X轴方向上相邻探针的高差记为Δz，Y轴方向上相邻探针的高差记为Δz；
将整个装置放在水里中，并将探针的针尖迎着水流方向；运行电脑上的数据处理软件Signal Express软件，在数据处理软件里面设置测量的时间为t和采样的频率为F，得到电压信号V随时间t的变化曲线；
(2)对Xm轴方向探针进行信号处理
将1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RX12)max，以及对应的时间tX12；然后依次分别对探针1#与3#、1#与4#……1#与m#探针所测的电压信号进行互相关分析，并获得最大的相关系数值(RX1m)max与对应的时间tX1m；
根据每个互相关曲线最大相关系数对应的时间tX12、tX13、……、tX1m，求得相应的速度大小VX12、VX13、……、VX1m，由于两个探针之间有一定固定的距离和固定的角度(根据ΔX、ΔY与ΔZ的距离确定)，因此可获得速度的矢量值：
Figure BDA0002933614220000091
Figure BDA0002933614220000092
Figure BDA0002933614220000093
其中，关于ΔLX的计算公如下：
Figure BDA0002933614220000094
Figure BDA0002933614220000095
根据本互相关分析方法，得到互相关曲线。分别对每个互相关曲线进行积分，积分的时间区域t>tX1m(即图中阴影区域)，如图4，得到相应的相关时间尺度大小TX12、TX13、……、TX1m值：
Figure BDA0002933614220000096
Figure BDA0002933614220000097
Figure BDA0002933614220000098
同样，根据互相关分析方法，得到1#探针与X方向上其他探针的互相关曲线，通过对互相关图进行积分，可获得涡体时间尺度TXX12、TXX13、……、TXX1m值：
Figure BDA0002933614220000099
Figure BDA0002933614220000101
Figure BDA0002933614220000102
分别将两探针之间的速度矢量与与涡体时间尺度进行乘积，可获得相应的涡体尺度的矢量值
Figure BDA0002933614220000103
Figure BDA0002933614220000104
Figure BDA0002933614220000105
Figure BDA0002933614220000106
将获得的涡体尺度的矢量值进行矢量相加，即可得到Xm方向上的涡体尺度矢量值；
Figure BDA0002933614220000107
3)对Yi轴方向布置的探针进行信号处理
将Y轴方向上的1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RY12)max，以及对应的时间tY12；然后依次分别对探针1#与3#、1#与 Y4#……1#与i#等探针所测的电压信号进行互相关分析，并获得最大的相关系数值 (RY1i)max与对应的时间tY1i；
根据每个互相关曲线最大相关系数对应的时间tY12、tY13、……、tY1i，可求得相应的速度大小VY12、VY13、……、VY1i，由于两个针之间有一定固定的距离和相应的夹角，因此可获得速度的矢量值；
Figure BDA0002933614220000108
Figure BDA0002933614220000109
Figure BDA0002933614220000111
其中，关于ΔLY的计算公如下：
Figure BDA0002933614220000112
Figure BDA0002933614220000113
分别对每个互相关曲线进行积分，积分的时间区域t>tY1m(即图中阴影区域)，如图4所示，得到相应的时间尺度大小TY12、TY13、……、TY1i值：
Figure BDA0002933614220000114
Figure BDA0002933614220000115
Figure BDA0002933614220000116
根据1#探针与Y方向上其他探针的互相关图，通过积分求解，可获得涡体时间尺度TYY12、TYY13、……、TYY1m值：
Figure BDA0002933614220000117
Figure BDA0002933614220000118
Figure BDA0002933614220000119
分别将两探针之间的速度矢量与与时间尺度进行乘积，可获得相应的涡体尺度的矢量值
Figure BDA00029336142200001110
Figure BDA00029336142200001111
Figure BDA00029336142200001112
Figure BDA0002933614220000121
将获得的涡体尺度的矢量值进行矢量相加，即得到Yi方向上的涡体尺度矢量值；
Figure BDA0002933614220000122
(4)采用同样方法获得Xn与Yj方向上涡体尺度的矢量值
Figure BDA0002933614220000123
和
Figure BDA0002933614220000124
最后将 X和Y轴四个方向Xm、Xn与Yi、Yj的涡体尺度的矢量值进行相加，可获得该测点处的涡体尺度的三维矢量值:
Figure BDA0002933614220000125
。

Claims (4)

1.测量复杂气液两相流中涡体尺度的方法，其特征在于，基于测量复杂气液两相流中涡体尺度的实验装置，该装置包括探针、基座、模拟信号-数字信号模块、探针供电电源和数据处理端，所述基座为“十”字交叉型，基座上沿十字型的横部和竖部分别等间距排列有若干探针，所述基座内部为空腔结构，所述探针竖立固定在基座上，探针针头朝向空中，探针尾部穿过基座***基座的空腔中，每根探针尾部在基座空腔内均连接有一根导线，每根导线从基座底部引出基座，再与模拟信号-数字信号模块连接，将探针采集的模拟信号转化为数字信号，所述模拟信号-数字信号模块通过数据传输线与数据处理端连接；所述探针分别与探针供电电源连接；所述探针由不锈钢外壳包裹银丝内芯构成，探针针尖部的银丝内芯暴露在外，探针尾部的银丝与对应的导线相连接，探针的不锈钢外壳与探针供电电源连接；基座十字型的横部为X轴方向，竖部为Y轴方向，所述探针相对基座的高度，从基座的十字交叉点上的探针开始，沿X轴正负方向递减，沿Y轴正负方向递减；X轴和Y轴上，相邻两根探针之间的高程差相等，相邻两根探针之间的水平间距相等；
所述方法包括以下步骤：
(1)数据采集
对探针进行编号，十字型基座的X轴的正负两个方向分别定义为Xm、Xn方向，Y轴的正负两个方向分别定为Yi、Yj方向；X轴中间探针标记为1#，1#探针一侧的探针逐次标记为2#、3#、……、m#，1#探针另一侧的探针逐次标记为2#、3#、……、n#；在基座的Y轴方向上，中间探针探标记为1#，1#探针一侧的探针逐次标记为2#、3#、……、i#，1#探针另一侧的探针逐次标记为2#、3#、……、j#；X轴方向上相邻探针的间距记为ΔX，Y轴方向上相邻探针的间距记为ΔY，ΔX＝ΔY；X轴方向上相邻探针的高程差记为Δz，Y轴方向上相邻探针的高程差记为Δz；
将整个装置放在水里中，并将探针的针尖迎着水流方向；运行电脑上的数据处理软件Signal Express软件，在数据处理软件里面设置测量的时间为t和采样的频率为F，得到电压信号V随时间t的变化曲线；
(2)对Xm轴方向探针进行信号处理
将1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RX12)max，以及对应的时间tX12；然后依次分别对探针1#与3#、1#与4#……1#与m#所测的电压信号进行互相关分析，并获得最大的相关系数值(RX1m)max与对应的时间tX1m；
根据每个互相关曲线最大相关系数对应的时间tX12、tX13、……、tX1m，求得相应的速度大小VX12、VX13、……、VX1m，由于两个探针之间有一定固定的距离和固定的角度，因此可获得速度的矢量值：
Figure FDA0003831705940000021
Figure FDA0003831705940000022
Figure FDA0003831705940000023
其中，关于ΔLX的计算公如下：
Figure FDA0003831705940000024
Figure FDA0003831705940000025
根据互相关分析方法，得到互相关曲线，分别对每个互相关曲线进行积分，积分的时间区域t>tX1m得到相应的相关时间尺度大小TX12、TX13、……、TX1m值：
Figure FDA0003831705940000026
Figure FDA0003831705940000027
Figure FDA0003831705940000028
同样，根据互相关分析方法，得到1#探针与X方向上其他探针的互相关曲线，通过对互相关图进行积分，获得涡体时间尺度TXX12、TXX13、……、TXX1m值：
Figure FDA0003831705940000029
Figure FDA0003831705940000031
Figure FDA0003831705940000032
分别将两探针之间的速度矢量与涡体时间尺度进行乘积，获得相应的涡体尺度的矢量值
Figure FDA0003831705940000033
Figure FDA0003831705940000034
Figure FDA0003831705940000035
Figure FDA0003831705940000036
将获得的涡体尺度的矢量值进行矢量相加，即得到Xm方向上的涡体尺度矢量值；
Figure FDA0003831705940000037
3)对Yi轴方向布置的探针进行信号处理
将Y轴方向上的1#与2#探针所获得电压信号进行互相关分析，获得最大的相关系数值(RY12)max，以及对应的时间tY12；然后依次分别对探针1#与3#、1#与Y4#……1#与i#等探针所测的电压信号进行互相关分析，并获得最大的相关系数值(RY1i)max与对应的时间tY1i；
根据每个互相关曲线最大相关系数对应的时间tY12、tY13、……、tY1i，得相应的速度大小VY12、VY13、……、VY1i，由于两个探针之间有一定固定的距离和相应的夹角，因此获得速度的矢量值；
Figure FDA0003831705940000038
Figure FDA0003831705940000039
Figure FDA0003831705940000041
其中，关于ΔLY的计算公如下：
Figure FDA0003831705940000042
Figure FDA0003831705940000043
分别对每个互相关曲线进行积分，积分的时间区域t>tY1i，得到相应的时间尺度大小TY12、TY13、……、TY1i值：
Figure FDA0003831705940000044
Figure FDA0003831705940000045
Figure FDA0003831705940000046
根据1#探针与Y方向上其他探针的互相关图，通过积分求解，获得涡体时间尺度TYY12、TYY13、……、TYY1i值：
Figure FDA0003831705940000047
Figure FDA0003831705940000048
Figure FDA0003831705940000049
分别将两探针之间的速度矢量与时间尺度进行乘积，获得相应的涡体尺度的矢量值
Figure FDA00038317059400000410
Figure FDA00038317059400000411
Figure FDA00038317059400000412
Figure FDA0003831705940000051
将获得的涡体尺度的矢量值进行矢量相加，即得到Yi方向上的涡体尺度矢量值；
Figure FDA0003831705940000052
(4)采用同样方法获得Xn与Yj方向上涡体尺度的矢量值
Figure FDA0003831705940000053
和
Figure FDA0003831705940000054
最后将X和Y轴四个方向Xm、Xn与Yi、Yj的涡体尺度的矢量值进行相加，获得该测点处的涡体尺度的三维矢量值:
Figure FDA0003831705940000055
2.根据权利要求1所述测量复杂气液两相流中涡体尺度的方法，其特征在于，基座横部和竖部上的探针数量相等，为4～10根。
3.根据权利要求1所述测量复杂气液两相流中涡体尺度的方法，其特征在于，基座的横部和竖部为与探针朝向方向一致的流线型。
4.根据权利要求1所述测量复杂气液两相流中涡体尺度的方法，其特征在于，所述探针供电电源为5V或10V电压输出的电源转换模块。
CN202110156276.9A 2021-02-04 2021-02-04 测量复杂气液两相流中涡体尺度的实验装置及方法 Active CN112985758B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110156276.9A CN112985758B (zh) 2021-02-04 2021-02-04 测量复杂气液两相流中涡体尺度的实验装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110156276.9A CN112985758B (zh) 2021-02-04 2021-02-04 测量复杂气液两相流中涡体尺度的实验装置及方法

Publications (2)

Publication Number Publication Date
CN112985758A CN112985758A (zh) 2021-06-18
CN112985758B true CN112985758B (zh) 2022-12-23

Family

ID=76347067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110156276.9A Active CN112985758B (zh) 2021-02-04 2021-02-04 测量复杂气液两相流中涡体尺度的实验装置及方法

Country Status (1)

Country Link
CN (1) CN112985758B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI20001685A0 (fi) * 2000-07-19 2000-07-19 Tr Tech Int Oy Mittausjärjestelmä ja menetelmä hiukkasnopeuden ja/tai hiukkasnopeusjakauman ja/tai hiukkaskoon ja/tai hiukkaskokojakauman mittaamiseksi
CN102998343B (zh) * 2012-12-05 2014-10-15 中国科学技术大学 一种基于阵列式单极电导探针的两相流层析成像***
CN107389742A (zh) * 2017-06-19 2017-11-24 天津大学 用于水包油乳状液测量的截面分布式阵列传感器及方法
CN108534835B (zh) * 2018-05-07 2020-05-19 中国核动力研究设计院 两相流界面参数测量方法
CN109596526A (zh) * 2018-11-29 2019-04-09 天津大学 双光纤阵列传感器气液两相流气相分布参数测量方法
CN109916947A (zh) * 2019-03-28 2019-06-21 中国农业大学 基于热脉冲原理测量饱和土壤水流向和流速的方法和装置
CN110794164A (zh) * 2019-12-13 2020-02-14 中国科学院大学 强磁场下液态金属速度场高时空精度的测量***及方法

Also Published As

Publication number Publication date
CN112985758A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112343774B (zh) 一种漂浮式海上风力发电装置的大尺度模型试验***及制作方法
Taddei et al. Characterisation of drag and wake properties of canopy patches immersed in turbulent boundary layers
CN202599852U (zh) 气液两相泡状流中气泡的识别及测量装置
Zhang et al. Vector high-resolution marine turbulence sensor based on a MEMS bionic cilium-shaped structure
Armandei et al. Hydroelastic buffeting assessment over a vertically hinged flat plate
CN105092684A (zh) 一种气液两相流中微气泡体积浓度测量装置
CN103335813A (zh) 一种气流反相对冲集沙仪
CN112985758B (zh) 测量复杂气液两相流中涡体尺度的实验装置及方法
JP2018084537A (ja) 風計測装置
CN106767721B (zh) 一种峡谷狭深河段水深测量方法和装置
García-Baena et al. Drag reduction of a blunt body through reconfiguration of rear flexible plates
CN101762436A (zh) 液体密度测量装置和测量方法
CN107621532A (zh) 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法
CN113715982B (zh) 一种外悬式喷水推进船模自航试验测试***
Steiner et al. Large-scale vortex structures in turbulent wakes behind bluff bodies. Part 2. Far-wake structures
Li et al. Numerical investigation of flow motion and performance of a horizontal axis tidal turbine subjected to a steady current
CN110333047B (zh) 一种正交悬臂梁结构的二维湍流传感器
CN202075303U (zh) 一种微差压式高精度测风装置
CN202928632U (zh) V型内锥流量计
CN200976007Y (zh) 气象传感器风速探头装置
CN104568656A (zh) 液体密度测量装置和测量方法
CN101226132B (zh) 快速测量细胞尺寸的微芯片及测量方法
CN210442358U (zh) 一种测量含沙量的传感器装置
CN203772890U (zh) 一种手持中空流速仪
Myers et al. Design, commissioning and performance of a device to vary the turbulence in a recirculating flume

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
GR01 Patent grant
GR01 Patent grant