CN109100533B - 一种利用超声波测量叶栅通道气流速度场方法 - Google Patents

一种利用超声波测量叶栅通道气流速度场方法 Download PDF

Info

Publication number
CN109100533B
CN109100533B CN201811092621.1A CN201811092621A CN109100533B CN 109100533 B CN109100533 B CN 109100533B CN 201811092621 A CN201811092621 A CN 201811092621A CN 109100533 B CN109100533 B CN 109100533B
Authority
CN
China
Prior art keywords
blade
ultrasonic
sensors
measuring
airflow
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
CN201811092621.1A
Other languages
English (en)
Other versions
CN109100533A (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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201811092621.1A priority Critical patent/CN109100533B/zh
Publication of CN109100533A publication Critical patent/CN109100533A/zh
Application granted granted Critical
Publication of CN109100533B publication Critical patent/CN109100533B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明涉及一种利用超声波测量叶栅通道气流速度场方法,采用超声波声时法原理,利用超声波层析成像技术来实现。通过在叶栅通道内布置超声波传感器阵列,结合超声波信号同步扫描和接收装置,构建时差法超声波层析成像系统。当叶栅通道有气流流过时,超声波层析成像系统自动获取顺流、逆流状态下的超声波传播时间信号,并利用传播时间信号构建信号矩阵,最后利用重建算法反演得到叶栅内部叶片间的气流速度场。测量系统结构简单,方法可实现无损测量、在线测量,相比于其他探针式测量方式,测量过程中无需在流场插入探针,测量方法非接触,且不会对流场本身产生干扰。

Description

一种利用超声波测量叶栅通道气流速度场方法
技术领域
本发明涉及一种流体测试技术,特别涉及一种利用超声波测量叶栅通道气流速度场方法。
背景技术
流速测量是流体流动研究的基本要求之一。初期的流体力学研究采用管道平均流速法,通过计量一定时间内通过某个断面积的流体流量,换算流动速度,这无疑是不准确的,某些情况下(例如气体流动)甚至是不可能的。第一个准确的测量点流速的方法是“皮托管”法,根据伯努利定理,利用压力和速度的数量关系,通过测量某点的总压和静压,计算得到该点的流速。这种方法简单易行,但是空间分辨率低,同时皮托管的大小对流动有一定程度的“污染”,而且由于测压管路阻尼,只能测量稳态流动的流速,对非稳态流动是无能为力的。
叶栅通道内的气流流动是极其复杂的。除了高速和旋转这两个重要基本特征外,还存在各种复杂的流动现象。比如压气机内流动由于增压流动,会出现失速流动、并发展成喘振,危及机组的安全。燃气轮机燃烧室流动受热声震荡作用,可能发生爆震,破坏燃烧室的高温防护内衬。透平流动中,冷却工质喷射并与主流掺混,可能造成叶片局部热疲劳断裂。多级压气机和透平的动静叶片排之间的高速相对运动,可能发生气流诱导的叶片颤振,危及叶片安全等。这些问题的解决,都依赖于我们对相关流动现象的认识。只有在充分认识了流动规律后,才有可能主动控制相关流动,实现相应机组的高效安全可靠的运行。
目前,应用于叶栅通道气流速度场测量的非接触式方法多为基于示踪粒子的光学成像方法,尚未出现利用超声波层析成像技术测量叶栅通道气流速度场的报道。基于此,本发明提出一种利用超声波层析成像技术,基于超声波传播声时法原理测量叶栅通道内的气流速度,从而发明一种利用超声波测量叶栅通道内气流速度场的方法。
发明内容
本发明是针对现在叶栅通道内气流速度场测量方法存在的问题,提出了一种利用超声波测量叶栅通道气流速度场方法,利用超声波层析成像技术,通过测量声波在气流顺流、逆流不同情况下的声时信号,并构建信号矩阵,利用重建算法反演得到叶栅内部叶片间的气流速度场,实现叶栅通道气流速度场无损测量方法。
本发明的技术方案为:一种利用超声波测量叶栅通道气流速度场方法,具体包括如下步骤:
1)在叶栅通道内相邻叶片A和B上分别布置一排微型收发一体式超声波传感器,形成传感器阵列,传感器耦合层和叶片表面平齐;
2)计算机控制信号发生器发出触发脉冲,并由前置放大器进行触发脉冲功率放大,功率放大后的触发脉冲变成超声波激励脉冲,由多路选择开关控制接通传感器阵列中叶片A的某一路待激励超声波传感器,完成一次声波发射;传感器阵列中叶片B上的接收传感器也在多路选择开关控制下连接信号调理电路,进行信号放大、调理和甄别,并经采集卡送入到计算机中,完成传感器阵列中任一传感器的有效超声波渡越时间提取;
3)以扫描方式,叶片A上N个传感器依次发射声波,叶片B上N个传感器接收声波,完成传感器阵列中叶片A和B上的所有传感器有效超声波渡越时间提取,形成信号向量,实现N×N路径的顺流、逆流时间差值信号获取;
4)利用信号向量进行速度场重建,构建层析成像系统:在测量区内构建网格,叶片A和B之间的叶栅通道气流区域为测量区,其中每个网格交点作为节点D1,D2…DN的速度分解为相互垂直方向,即x,y方向;静态流体中声速为V0,流动引起的速度改变为ΔV,则对应总时差改变量ΔT为此声路径上每一网格时滞的积分量;在扫描过程中,叶栅叶片上布置的传感器对数为N,叶片A上N个传感器阵列依次发射声波发射为一次轮循,一次轮循可获得NxN条声路径上有效信号,并且每条路径的总时差该变量为路径上每个网格时滞的积分量,逐次完成N个传感器发射声波,获得NxN条路径上的有效信号采集;利用层析成像重建算法反演求解可得到气流速度场,气流速度场为反演得到的每个网格区域因气流流动带来的时差改变值组成。
所述每个网格区域因气流流动带来的时差改变值
Figure BDA0001804743710000031
其大小与气流流速和传感器之间距离相关,可表示为气流速度V与测量装置结构系数的乘积。
本发明的有益效果在于:本发明利用超声波测量叶栅通道气流速度场方法,基于超声声时法原理测量叶栅通道内气流速度场的方法,测量系统结构简单,方法可实现无损测量、在线测量,相比于其他探针式测量方式,测量过程中无需在流场插入探针,测量方法非接触,且不会对流场本身产生干扰。
附图说明
图1为本发明叶栅通道内相邻叶片传感器阵列布置示意图;
图2为本发明利用超声波测量叶栅通道气流速度场系统原理示意图。
具体实施方式
利用超声波测量叶栅通道气流速度场方法实现具体步骤如下:
1、如图1所示,在叶栅通道内相邻叶片A和B上分别布置一排微型收发一体式超声波传感器,形成传感器阵列,传感器耦合层和叶片表面平齐。
2、如图2所示,计算机1控制信号发生器2发出触发脉冲,并由前置放大器3进行触发脉冲功率放大,功率放大后的触发脉冲变成超声波激励脉冲,由多路选择开关4控制接通传感器阵列7中叶片A的某一路待激励超声波传感器,完成一次声波发射。传感器阵列7中叶片B上的接收传感器也在多路选择开关4控制下连接信号调理电路5,进行信号放大、调理和甄别,并经采集卡6送入到计算机1中,完成传感器阵列7中任一传感器的有效超声波渡越时间提取。
3、以扫描方式,叶片A上N个传感器依次发射声波,叶片B上N个传感器接收声波,完成传感器阵列中叶片A和B上的所有传感器有效超声波渡越时间提取,形成信号向量,实现N×N路径的顺流、逆流时间差值信号获取;
4、利用信号向量进行速度场重建,构建层析成像系统,在测量区Q(叶片A和B之间的叶栅通道气流区域)内构建一定数量网格,其中每个节点D1,D2…DN(网格交点为节点)的速度可以分解为相互垂直方向,即x,y方向。按照时差法原理,两个方向的气流速度均会对超声传播时间产生影响,静态流体中声速为V0,流动引起的速度改变为ΔV,则对应总时差改变量ΔT为此声路径上每一个网格时滞的积分量。在扫描过程中,一次轮循(叶片A上一个传感器一次声波发射,叶片B上N个传感器依次接收,那就有N条路径,叶片A上N个传感器依次发射,最终得到NxN条路径为一次轮循,从而实现N×N路径的顺流、逆流时间差值信号获取。)可获得NxN条声路径上有效信号,每条路径的总时差该变量都为路径上每个网格时滞的积分量,逐次完成所有路径上的有效信号采集。最后,利用层析成像重建算法反演求解可得到气流速度场,气流速度场为反演得到的每个网格因气流流动带来的时差改变值组成。
一个传感器一次声波发射,依次N个传感器接收,那就有N条路径,其中某一路径方向上可以写出声波在介质中时差一般形式:
Figure BDA0001804743710000041
其中,ΔT为某一路径方向上的总时差改变量(通过收发传感器、发射接收时间等可知),
Figure BDA0001804743710000042
Figure BDA0001804743710000043
分别为此路径对应的发射器和接收器的位置(已知),
Figure BDA0001804743710000044
为此声路径上某一网格区域因气流流动带来的时差改变值(所求),l为此声路径总长(已知)。
考虑一次声波发射所有声路径方向有效信号,将积分方程离散化,可得到时差信号的线性方程组表达式:
Figure BDA0001804743710000045
其中,N为叶栅叶片上布置的传感器对数,其大小受叶片宽度、传感器尺寸、以及成像精度约束。M为叶栅相邻叶片之间划分网格的数量。其大小受叶栅相邻叶片之间距离的约束,
Figure BDA0001804743710000046
实质是气流在网格处流动而引起的声渡越时间差,其大小与气流流速和传感器之间距离相关,可以表示为气流速度V与测量装置结构系数的乘积。因此,可构造出一个未知数为网格节点处速度参数的线性方程组:
TN×1=SN×M·VM×1
其中,S为与测量装置结构和气流对象有关的敏感系数矩阵(通过正演模型设计),V为测量区内网格节点处的速度参数,T为测量时得到的声时信号向量。方程的典型求解算法可以是Landweber迭代算法、Tikhonov正则化算法、牛顿-拉夫逊算法及共轭梯度算法等。
通过奇异值分解的算法能够获得稳定且可靠的解,同时具有收敛速度快的优点。此外,在实际问题时,亦可以通过截断奇异值分解实现降低矩阵维数,减小运算量、增大求解速度。
本发明基于超声声时法原理测量叶栅通道内气流速度场的方法,测量系统结构简单,方法可实现无损测量、在线测量,相比于其他探针式测量方式,测量过程中无需在流场插入探针,测量方法非接触,且不会对流场本身产生干扰。

Claims (2)

1.一种利用超声波测量叶栅通道气流速度场方法,其特征在于,具体包括如下步骤:
1)在叶栅通道内相邻叶片A和B上分别布置一排微型收发一体式超声波传感器,形成传感器阵列,传感器耦合层和叶片表面平齐;
2)计算机控制信号发生器发出触发脉冲,并由前置放大器进行触发脉冲功率放大,功率放大后的触发脉冲变成超声波激励脉冲,由多路选择开关控制接通传感器阵列中叶片A的某一路待激励超声波传感器,完成一次声波发射;传感器阵列中叶片B上的接收传感器也在多路选择开关控制下连接信号调理电路,进行信号放大、调理和甄别,并经采集卡送入到计算机中,完成传感器阵列中任一传感器的有效超声波渡越时间提取;
3)以扫描方式,叶片A上N个传感器依次发射声波,叶片B上N个传感器接收声波,完成传感器阵列中叶片A和B上的所有传感器有效超声波渡越时间提取,形成信号向量,实现N×N路径的顺流、逆流时间差值信号获取;
4)利用信号向量进行速度场重建,构建层析成像系统:在测量区内构建网格,叶片A和B之间的叶栅通道气流区域为测量区,其中每个网格交点作为节点,记为D1,D2…DN,每个节点的速度分解为相互垂直方向,即x,y方向;静态流体中声速为V0,流动引起的速度改变为ΔV,则某一路径方向上的总时差改变量为ΔT,ΔT为对应此声路径上每一网格时滞的积分量;在扫描过程中,叶栅叶片上布置的传感器对数为N,叶片A上一个传感器一次声波发射,叶片B上N个传感器依次接收,叶片A上N个传感器依次发射声波,叶片A上N个传感器依次发射完为一次轮循,一次轮循可获得N× N条声路径上有效信号,并且每条声路径的总时差改变量为此声路径上每个网格时滞的积分量,逐次完成N个传感器发射声波,获得N× N条路径上的有效信号采集;利用层析成像重建算法反演求解可得到气流速度场,气流速度场为反演得到的每个网格区域因气流流动带来的时差改变值组成。
2.根据权利要求1所述利用超声波测量叶栅通道气流速度场方法,其特征在于,每个网格区域因气流流动带来的时差改变值
Figure FDA0002504104790000021
其大小与气流流速和传感器之间距离相关,可表示为气流速度V与测量装置结构系数的乘积。
CN201811092621.1A 2018-09-19 2018-09-19 一种利用超声波测量叶栅通道气流速度场方法 Active CN109100533B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811092621.1A CN109100533B (zh) 2018-09-19 2018-09-19 一种利用超声波测量叶栅通道气流速度场方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811092621.1A CN109100533B (zh) 2018-09-19 2018-09-19 一种利用超声波测量叶栅通道气流速度场方法

Publications (2)

Publication Number Publication Date
CN109100533A CN109100533A (zh) 2018-12-28
CN109100533B true CN109100533B (zh) 2020-11-24

Family

ID=64866713

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811092621.1A Active CN109100533B (zh) 2018-09-19 2018-09-19 一种利用超声波测量叶栅通道气流速度场方法

Country Status (1)

Country Link
CN (1) CN109100533B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114814284B (zh) * 2022-04-21 2023-11-07 上海理工大学 一种反射式安装的超声阵列测量叶栅流场的方法与装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1774617A (zh) * 2003-01-13 2006-05-17 塞德拉公司 使用超声波传感器阵列确定管道内的流体速度的设备和方法
CN101852629A (zh) * 2010-07-14 2010-10-06 唐山汇中威顿仪表有限公司 一种渠道或河道用超声波流量传感器阵列及其布置方法
CN102607653A (zh) * 2012-04-06 2012-07-25 清华大学 一种基于超声波的管道流速成像和流量测量方法
CN107102166A (zh) * 2017-03-28 2017-08-29 天津大学 超声多普勒多相流流速分布检测设备

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7237440B2 (en) * 2003-10-10 2007-07-03 Cidra Corporation Flow measurement apparatus having strain-based sensors and ultrasonic sensors
US7673526B2 (en) * 2006-11-01 2010-03-09 Expro Meters, Inc. Apparatus and method of lensing an ultrasonic beam for an ultrasonic flow meter
US9989387B2 (en) * 2014-04-01 2018-06-05 Saudi Arabian Oil Company Flow data acquisition and telemetry processing systems

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1774617A (zh) * 2003-01-13 2006-05-17 塞德拉公司 使用超声波传感器阵列确定管道内的流体速度的设备和方法
CN101852629A (zh) * 2010-07-14 2010-10-06 唐山汇中威顿仪表有限公司 一种渠道或河道用超声波流量传感器阵列及其布置方法
CN102607653A (zh) * 2012-04-06 2012-07-25 清华大学 一种基于超声波的管道流速成像和流量测量方法
CN107102166A (zh) * 2017-03-28 2017-08-29 天津大学 超声多普勒多相流流速分布检测设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于GMRES的超声过程层析成像算法研究;刘金全 等;《传感技术学报》;20131031;第26卷(第10期);第1395-1400页 *
气固两相流超声过程层析成像系统研究;田昌 等;《仪器仪表学报》;20160331;第37卷(第3期);第586-592页 *

Also Published As

Publication number Publication date
CN109100533A (zh) 2018-12-28

Similar Documents

Publication Publication Date Title
Ostermann et al. Properties of a sweeping jet emitted from a fluidic oscillator
US9746360B2 (en) Nonintrusive performance measurement of a gas turbine engine in real time
Dennis et al. On the limitations of Taylor's hypothesis in constructing long structures in a turbulent boundary layer
CN103969022A (zh) 一种高超声速风洞湍流度间接测量方法
JP2010117360A (ja) 補正した空気流量の直接的非貫入型測定のためのシステム及び方法
CN107085123A (zh) 一种基于声学技术的多路径管道流速测量装置及其方法
JPH11505020A (ja) 空気データ装置
Gil Synthetic jet Reynolds number based on reaction force measurement
CN109100533B (zh) 一种利用超声波测量叶栅通道气流速度场方法
CN102288779B (zh) 一种高精度抗干扰超声波风速风向测量方法
CN102914589B (zh) 利用超声波检测甲烷浓度的检测方法
Marusic et al. A wall-shear stress predictive model
CN106500951A (zh) 测量高超声速气流参数的测量探头、测量系统和方法
Chen et al. Effect of transducer shadowing of ultrasonic anemometers on wind velocity measurement
Gu et al. Pressure prediction model based on artificial neural network optimized by genetic algorithm and its application in quasi-static calibration of piezoelectric high-pressure sensor
CN105300654B (zh) 一种低速风洞内的超声精确测距系统
JP4688253B2 (ja) 超音波流量計
Menna et al. The Mean Flow Structure Around and Within a Turbulent Junction or Horseshoe Vortex—Part I: The Upstream and Surrounding Three-Dimensional Boundary Layer
Schmirler et al. Double probe recovery temperature Anemometry
Wang et al. Structure optimal design and performance test of airfoil shear probes
US9200982B2 (en) Phased array turbomachine monitoring system
Bakhtar et al. Calibration characteristics of a three-hole probe and a static tube in wet steam
Schlienger et al. Measuring Unsteady 3D Flow with a single pressure transducer
Pettinari et al. Time-of-flight estimation in acoustic pyrometry: sensitivity to pulse characteristics
Hobeck et al. Design and analysis of dual pressure probes for predicting turbulence-induced vibration in low velocity flow

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