CN114740222B - Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids - Google Patents

Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids Download PDF

Info

Publication number
CN114740222B
CN114740222B CN202210255200.6A CN202210255200A CN114740222B CN 114740222 B CN114740222 B CN 114740222B CN 202210255200 A CN202210255200 A CN 202210255200A CN 114740222 B CN114740222 B CN 114740222B
Authority
CN
China
Prior art keywords
laser
speed camera
grid
time
velocity field
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
CN202210255200.6A
Other languages
Chinese (zh)
Other versions
CN114740222A (en
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202210255200.6A priority Critical patent/CN114740222B/en
Publication of CN114740222A publication Critical patent/CN114740222A/en
Application granted granted Critical
Publication of CN114740222B publication Critical patent/CN114740222B/en
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/26Measuring 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 optical wave

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明公开了一种动静叶栅间三维速度场均匀性测量装置及测量方法,属于多相流测试领域。本发明提出了一种动静叶栅间三维速度场均匀性测量装置及测量方法,利用光的干涉设计了光路系统,在待测流场当中形成交错激光网格,实现示踪分子的标记与激发。利用同步控制器实现增强型高速相机和脉冲激光器的时序控制,可实现动静叶栅间三维速度场及其均匀性的同步测量。同时,本发明根据激光光强特性重构磷光网格图像,可将速度场精度提高到亚像素精度。

Figure 202210255200

The invention discloses a measuring device and a measuring method for the uniformity of a three-dimensional velocity field between dynamic and static blade cascades, belonging to the field of multiphase flow testing. The invention proposes a three-dimensional velocity field uniformity measurement device and measurement method between dynamic and static cascades. The optical path system is designed by using light interference, and an interlaced laser grid is formed in the flow field to be measured to realize the marking and excitation of tracer molecules. . The synchronous controller is used to realize the timing control of the enhanced high-speed camera and the pulse laser, and the synchronous measurement of the three-dimensional velocity field and its uniformity between the moving and static cascades can be realized. At the same time, the invention reconstructs the phosphorescence grid image according to the characteristics of laser light intensity, which can improve the precision of the velocity field to sub-pixel precision.

Figure 202210255200

Description

一种动静叶栅间三维速度场均匀性测量装置及测量方法A device and method for measuring the uniformity of three-dimensional velocity field between dynamic and static cascades

技术领域technical field

本发明属于多相流测试领域,尤其是一种动静叶栅间三维速度场均匀性测量装置及测量方法。The invention belongs to the field of multiphase flow testing, in particular to a measuring device and a measuring method for the uniformity of a three-dimensional velocity field between dynamic and static blade cascades.

背景技术Background technique

在现代航空航天领域,发动机中动静叶栅间流速的高精度测量以及流场稳定性测量具有重大的意义。但是动静叶栅间空间小、速度大,成为测量的难点(动静叶栅间距8~10mm,转速可达几千转到几万转)。在加工过程中由于制造公差,会造成流场的不稳定,需要对流场周期内均匀性进行高精度测量。传统的测量方法,如:毕托管、热风线速仪,对流场干扰大,响应慢,难以适应小空间中高速流场的测量。近几年来,随着光学设备和高速摄影技术发展,一些光学非接触测速技术也得到了快速的发展,如:分子标记测速技术、激光多普勒测速法(LDV)和粒子图像测速法(PIV)等。其中,激光多普勒测速法和粒子图像测速法需要在流场中加入示踪颗粒,因为示踪颗粒的加入,会影响叶轮的正常运转,且由于颗粒的惯性较大,颗粒的滞后性很大,同时对流场的均匀性测量难以实现。In the field of modern aerospace, the high-precision measurement of the flow velocity between the dynamic and static cascades in the engine and the measurement of the stability of the flow field are of great significance. However, the space between the dynamic and static cascades is small and the speed is high, which has become a difficulty in measurement (the distance between the dynamic and static cascades is 8-10mm, and the speed can reach several thousand to tens of thousands of revolutions). Due to manufacturing tolerances in the processing process, the flow field will be unstable, and high-precision measurement of the uniformity of the flow field within a cycle is required. Traditional measurement methods, such as Pitot tubes and hot air linear velocity meters, have great interference to the flow field and slow response, making it difficult to adapt to the measurement of high-speed flow fields in small spaces. In recent years, with the development of optical equipment and high-speed photography technology, some optical non-contact velocimetry technologies have also been developed rapidly, such as: molecular marker velocimetry technology, laser Doppler velocimetry (LDV) and particle image velocimetry (PIV). )wait. Among them, laser Doppler velocimetry and particle image velocimetry need to add tracer particles in the flow field, because the addition of tracer particles will affect the normal operation of the impeller, and due to the large inertia of the particles, the hysteresis of the particles is very large. It is difficult to measure the uniformity of the flow field at the same time.

发明内容Contents of the invention

本发明的目的在于克服上述现有技术的缺点,提供一种动静叶栅间三维速度场均匀性测量装置及测量方法。The object of the present invention is to overcome the above-mentioned shortcomings of the prior art, and provide a three-dimensional velocity field uniformity measurement device and measurement method between dynamic and static blade cascades.

为达到上述目的,本发明采用以下技术方案予以实现:In order to achieve the above object, the present invention adopts the following technical solutions to achieve:

一种动静叶栅间三维速度场均匀性测量装置,包括脉冲激光器和光路系统,脉冲激光器用于发射激光,所述激光进入光路系统内改变光路,之后经光路系统的出口变成两束薄片状激光,所述激光薄片进入动叶片和静叶片之间待测流场中形成交错激光网格;A device for measuring the uniformity of a three-dimensional velocity field between moving and static cascades, including a pulse laser and an optical path system. The pulse laser is used to emit laser light. The laser light enters the optical path system to change the optical path, and then passes through the exit of the optical path system to become two beams of thin sheets. Laser, the laser sheet enters the flow field to be measured between the moving blade and the stationary blade to form a staggered laser grid;

所述待测流场中采用联乙酰或丙酮作为示踪分子,待测流场为无氧环境;In the flow field to be measured, biacetyl or acetone is used as a tracer molecule, and the flow field to be measured is an oxygen-free environment;

待测流场的两侧分别设有第一增强型高速相机和第二增强型高速相机,第一增强型高速相机和第二增强型高速相机分别通过连接线连接有同步控制器,同步控制器连接有图像处理系统;A first enhanced high-speed camera and a second enhanced high-speed camera are respectively arranged on both sides of the flow field to be measured, and the first enhanced high-speed camera and the second enhanced high-speed camera are respectively connected to a synchronous controller through a connection line, and the synchronous controller Connected with an image processing system;

所述第一增强型高速相机和第二增强型高速相机用于同步的拍摄示踪分子受激发后产生的磷光图像。The first enhanced high-speed camera and the second enhanced high-speed camera are used for synchronously shooting phosphorescence images generated after the tracer molecules are excited.

进一步的,光路系统包括反光镜、分光镜和空气劈尖;Further, the optical path system includes mirrors, beam splitters and air wedges;

当一束激光通过空气劈尖后,分成多个薄片状激光,随后通过分光镜、反光镜将一束激光分成两束激光,两束激光进入待测流场中形成交错激光网格。When a laser beam passes through the air wedge, it is divided into multiple sheet-shaped lasers, and then the laser beam is divided into two laser beams by a beam splitter and a mirror. The two laser beams enter the flow field to be measured to form an interlaced laser grid.

进一步的,空气劈尖包括第一石英玻璃和第二石英玻璃,第一石英玻璃和第二石英玻璃一端相接触,另一端之间设有金属丝;Further, the air wedge includes first quartz glass and second quartz glass, one end of the first quartz glass and the second quartz glass are in contact, and a metal wire is arranged between the other ends;

两个石英玻璃之间形成一个空气薄膜,所述空气薄膜利用光的干涉将一束激光分成多个薄片状;An air film is formed between the two quartz glasses, and the air film uses light interference to divide a beam of laser light into multiple thin slices;

激光能够在第一石英玻璃的上表面激光全穿透,在下表面一部分发生发射,另一部分发生折射,在第二石英玻璃上表面发生全发射。The laser can fully penetrate the upper surface of the first quartz glass, emit on a part of the lower surface, refract the other part, and fully emit on the upper surface of the second quartz glass.

一种动静叶栅间三维速度场均匀性测量方法,基于本发明的动静叶栅间三维速度场均匀性测量装置的进行,包括以下步骤:A method for measuring the uniformity of the three-dimensional velocity field between static and dynamic cascades, based on the implementation of the device for measuring the uniformity of the three-dimensional velocity field between the static and dynamic cascades of the present invention, comprising the following steps:

步骤一、进行标定实验,计算第一增强型高速相机、第二增强型高速相机与实际空间位置关系,从而得到图像像素坐标aStep 1. Carry out a calibration experiment, and calculate the relationship between the first enhanced high-speed camera, the second enhanced high-speed camera and the actual space position, so as to obtain the image pixel coordinate a

32`和实际空间坐标之间的转换关系式:32` and the conversion relationship between the actual space coordinates:

Figure GDA0004001803650000031
Figure GDA0004001803650000031

式中:(u,v)代表像素坐标;(Xw,Yw,Zw)代表空间实际坐标;kij代表求解参数,i=1,2,3,j=1,2,3,4;In the formula: (u, v) represents the pixel coordinates; (X w , Y w , Z w ) represents the actual space coordinates; k ij represents the solution parameters, i=1,2,3, j=1,2,3,4 ;

步骤二、在示踪分子的磷光寿命内,利用第一增强型高速相机和第二增强型高速相机获得两个连续时刻t1、t2的磷光图像,分别表示为:

Figure GDA0004001803650000032
Step 2. Within the phosphorescence lifetime of the tracer molecule, use the first enhanced high-speed camera and the second enhanced high-speed camera to obtain two consecutive phosphorescence images at time t 1 and t 2 , which are expressed as:
Figure GDA0004001803650000032

步骤三、对于t1时刻第一增强型高速相机、第二增强型高速相机获得的磷光图像

Figure GDA0004001803650000033
和/>
Figure GDA0004001803650000034
根据磷光最大强度点进行网格交点一一匹配,得到实际空间中同一个网格交点在两个不同相机拍摄图片中像素点的坐标,随后将两个像素点坐标带入式,计算出实际网格空间点的坐标(Xw,Yw,Zw)t1;Step 3: Phosphorescence images obtained by the first enhanced high-speed camera and the second enhanced high-speed camera at time t1
Figure GDA0004001803650000033
and />
Figure GDA0004001803650000034
Match the grid intersection points one by one according to the phosphorescence maximum intensity point, and obtain the pixel coordinates of the same grid intersection point in the actual space in the pictures taken by two different cameras, and then bring the coordinates of the two pixel points into the formula to calculate the actual grid The coordinates of the grid space point (X w , Y w , Z w ) t1 ;

步骤四、分别对第一增强型高速相机、第二增强型高速相机获得两个时刻的磷光图像

Figure GDA0004001803650000035
和/>
Figure GDA0004001803650000036
和/>
Figure GDA0004001803650000037
进行互相关计算,得到同一个相机捕捉到两个不同时刻磷光网格图片中网格交点的对应关系;利用互相关算法计算得到t2时刻两个相机磷光图像/>
Figure GDA0004001803650000038
和/>
Figure GDA0004001803650000039
中网格交点的像素点的坐标,带入式,计算出t2时刻网格交点在实际空间中的坐标(Xw,Yw,Zw)t2;Step 4. Obtain phosphorescence images at two moments for the first enhanced high-speed camera and the second enhanced high-speed camera respectively
Figure GDA0004001803650000035
and />
Figure GDA0004001803650000036
and />
Figure GDA0004001803650000037
Carry out cross-correlation calculations to obtain the corresponding relationship between grid intersection points in two phosphorescent grid pictures captured by the same camera at different times; use cross-correlation algorithm to calculate phosphorescence images of two cameras at time t2 >
Figure GDA0004001803650000038
and />
Figure GDA0004001803650000039
The coordinates of the pixel points of the grid intersection in the middle are brought into the formula to calculate the coordinates (X w , Y w , Z w ) t2 of the grid intersection in the actual space at time t2 ;

步骤五、根据t1时刻和t2时刻网格交点在实际空间中的位置(Xw,Yw,Zw)t1、(Xw,Yw,Zw)t2与时间的关系,计算得到流场的三维速度场为

Figure GDA00040018036500000310
Step 5. According to the relationship between the position (X w , Y w , Z w ) t1 and (X w , Y w , Z w ) t2 of the grid intersection point in the actual space at time t 1 and time t 2 and the time, calculate The three-dimensional velocity field of the flow field is
Figure GDA00040018036500000310

进一步的,步骤二还包括:Further, step 2 also includes:

对所述磷光图像进行去噪处理,同时增强图像对比度;performing denoising processing on the phosphorescence image while enhancing image contrast;

对于单条磷光网格图像,在垂直于网格线的方向采用高斯函数对磷光光束进行拟合。For a single phosphorescent grid image, a Gaussian function is used to fit the phosphorescent beam in the direction perpendicular to the grid lines.

进一步的,通过对脉冲激光器、第一增强型高速相机和第二增强型高速相机的时序控制,使得对每个动叶片经过同一个静叶片时的速度场进行测量,通过比较不同动叶片之间的速度场的差异,从而实现流场的周期内均匀性的测量。Further, through the timing control of the pulsed laser, the first enhanced high-speed camera and the second enhanced high-speed camera, the velocity field of each moving blade passing the same stationary blade can be measured, and by comparing the The difference of the velocity field, so as to realize the measurement of the uniformity of the flow field in the cycle.

进一步的,测量速度场均匀性具体操作为:Further, the specific operation of measuring the uniformity of the velocity field is as follows:

当动叶片运动的转速为n,动叶片个数为k,则动叶片旋转一周的时间为

Figure GDA0004001803650000041
此时同一个静叶片经历2个相邻的叶片所需要的时间为/>
Figure GDA0004001803650000042
此时脉冲激光器的频率为f=nk,高速相机拍摄时间/>
Figure GDA0004001803650000043
When the rotating speed of the moving blade is n and the number of moving blades is k, the time for the moving blade to rotate one circle is
Figure GDA0004001803650000041
At this time, the time required for the same stationary blade to experience two adjacent blades is />
Figure GDA0004001803650000042
At this time, the frequency of the pulsed laser is f=nk, and the shooting time of the high-speed camera />
Figure GDA0004001803650000043

相邻两次脉冲激光激发时间满足:

Figure GDA0004001803650000044
The excitation time of two adjacent pulse lasers satisfies:
Figure GDA0004001803650000044

Figure GDA0004001803650000045
时,脉冲激光器的频率为/>
Figure GDA0004001803650000046
高速摄像机拍摄时间
Figure GDA0004001803650000047
when
Figure GDA0004001803650000045
, the frequency of the pulsed laser is />
Figure GDA0004001803650000046
High-speed camera shooting time
Figure GDA0004001803650000047

进一步的,在t1时刻,磷光网格图像初始分布布置在两个动叶片中间位置。Further, at time t 1 , the initial distribution of the phosphorescent grid images is arranged in the middle of the two rotor blades.

与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:

本发明提出的一种动静叶栅间三维速度场均匀性测量装置及测量方法,利用光的干涉设计了光路系统,在待测流场当中形成交错激光网格,实现示踪分子的标记与激发。利用同步控制器实现增强型高速相机和脉冲激光器的时序控制,可实现动静叶栅间三维速度场及其均匀性的同步测量。同时,本发明根据激光光强特性重构磷光网格图像,可将速度场精度提高到亚像素精度。The invention proposes a three-dimensional velocity field uniformity measurement device and measurement method between dynamic and static cascades. The optical path system is designed by using light interference, and an interlaced laser grid is formed in the flow field to be measured to realize the marking and excitation of tracer molecules. . The synchronous controller is used to realize the timing control of the enhanced high-speed camera and the pulse laser, and the synchronous measurement of the three-dimensional velocity field and its uniformity between the moving and static cascades can be realized. At the same time, the invention reconstructs the phosphorescence grid image according to the characteristics of laser light intensity, which can improve the precision of the velocity field to sub-pixel precision.

进一步,根据激光的波长设计了空气劈尖,可进一步提高标记激光的质量,提高网格的空间分辨率,且在流场中形成一个高质量的网格磷光图像。Furthermore, the air wedge is designed according to the wavelength of the laser, which can further improve the quality of the marking laser, improve the spatial resolution of the grid, and form a high-quality grid phosphorescence image in the flow field.

进一步,本发明根据动静叶栅的转速、半径、叶片数,设计了磷光拍摄时序,确定了相机频率、脉冲激光器频率、以及相机曝光时间等,可实现流场均匀性的测量,并提高的测量的精确性。Further, according to the rotating speed, radius, and blade number of the static and dynamic cascades, the present invention designs the phosphorescence shooting sequence, determines the camera frequency, pulse laser frequency, and camera exposure time, etc., so that the measurement of the uniformity of the flow field can be realized, and the measurement accuracy can be improved. the accuracy.

进一步,本发明利用分子标记测速技术测量动静叶栅间三维速度场,对流场无干扰,测量精度高。Further, the present invention uses the molecular marker velocity measurement technology to measure the three-dimensional velocity field between the dynamic and static blade cascades, without interference to the flow field, and has high measurement accuracy.

附图说明Description of drawings

图1为动静叶栅间三维速度场测量装置的结构图;Figure 1 is a structural diagram of the three-dimensional velocity field measurement device between the moving and static cascades;

图2为光路系统的结构图;Fig. 2 is the structural diagram of optical path system;

图3为空气劈尖的结构图;Fig. 3 is the structural diagram of air wedge;

图4为磷光网格图像。Figure 4 is an image of the phosphorescent grid.

其中:01-脉冲激光器;02-激光;03-第一数据传输线;04-第一增强型高速相机;05-光路系统;06-动叶片;07-第二增强型高速相机;08-第二数据传输线;09-静叶片;10-同步控制器;11-第三数据传输线;12-图像处理系统;01-脉冲激光器;13-反光镜;14-分光镜;15-空气劈尖;16-第一石英玻璃;17-金属丝;18-第二石英玻璃。Among them: 01-pulse laser; 02-laser; 03-first data transmission line; 04-first enhanced high-speed camera; 05-optical system; 06-moving blade; 07-second enhanced high-speed camera; 08-second Data transmission line; 09-stator blade; 10-synchronous controller; 11-third data transmission line; 12-image processing system; 01-pulse laser; 13-mirror; 14-beam splitter; 15-air splitter; 16- First quartz glass; 17 - wire; 18 - second quartz glass.

具体实施方式Detailed ways

为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to enable those skilled in the art to better understand the solutions of the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only It is an embodiment of a part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts shall fall within the protection scope of the present invention.

需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first" and "second" in the description and claims of the present invention and the above drawings are used to distinguish similar objects, but not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used are interchangeable under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having", as well as any variations thereof, are intended to cover a non-exclusive inclusion, for example, a process, method, system, product or device comprising a sequence of steps or elements is not necessarily limited to the expressly listed instead, may include other steps or elements not explicitly listed or inherent to the process, method, product or apparatus.

本发明采用分子标记测速技术实现动静叶栅间三维速度场和流场均匀性的同步测量。其中,分子标记测速技术是一个非侵入式、非接触式激光诊断技术。其原理是采用特定波长的激光激发流场中示踪分子发出磷光,一定时间间隔后直接捕捉磷光信号确定示踪分子的位置,进而根据示踪分子位移-时间关系求解速度。分子标记测速技术能够有效的克服激光多普勒测速法和粒子图像测速法的缺点,能够很好的应用于动静叶栅间速度场的高精度测量,同时采用时序控制,可实现流场周期内均匀性测量。本发明基于分子标记测速技术提出了一种动静叶栅间三维速度场均匀性测量装置及测量方法,可克服以往测试方法的缺点,实现动静叶栅间三维速度场和速度场均匀性的同步测量。The invention adopts the molecular marker velocity measurement technology to realize the synchronous measurement of the three-dimensional velocity field and the uniformity of the flow field between the dynamic and static blade cascades. Among them, the molecular marker velocimetry technology is a non-invasive, non-contact laser diagnostic technology. The principle is to use a laser with a specific wavelength to excite tracer molecules in the flow field to emit phosphorescence. After a certain time interval, the phosphorescence signal is directly captured to determine the position of the tracer molecules, and then the velocity is calculated according to the displacement-time relationship of the tracer molecules. Molecular marker velocimetry technology can effectively overcome the shortcomings of laser Doppler velocimetry and particle image velocimetry, and can be well applied to high-precision measurement of the velocity field between dynamic and static cascades. Uniformity measurement. The present invention proposes a three-dimensional velocity field uniformity measurement device and measurement method between dynamic and static cascades based on molecular marker velocity measurement technology, which can overcome the shortcomings of previous testing methods and realize simultaneous measurement of the three-dimensional velocity field and velocity field uniformity between dynamic and static cascades .

下面结合附图对本发明做进一步详细描述:The present invention is described in further detail below in conjunction with accompanying drawing:

参见图1,图1为动静叶栅间三维速度场测量装置的结构图,一种动静叶栅间三维速度场的高精度测量装置,由动静叶栅结构、光路系统、图像捕捉系统、图像处理四部分组成。动静叶栅结构由动叶片06和静叶片09组成,该结构可存在于发动机,燃气轮机当中。光路系统主要由脉冲激光器01、反光镜13、空气劈尖15和分光镜14构成,用于在流场待测区域中构建交错激光网格,实现对待测流场内示踪分子的激发与标记。图像捕捉系统主要由同步控制器10、第一增强型高速相机04和第二增强型高速相机07构成,可获得连续两个不同时刻示踪分子的磷光网格图像。图像处理系统12主要根据示踪分子的磷光网格图像重建磷光网格三维模型,根据磷光强度分布特性,开发互相关算法得到不同时刻网格空间位置之间关系,从而得到三维速度场。Refer to Fig. 1, Fig. 1 is a structural diagram of the three-dimensional velocity field measuring device between the dynamic and static cascades, a high-precision measurement device for the three-dimensional velocity field between the dynamic and static cascades. It consists of four parts. The moving and static cascade structure is composed of moving blades 06 and stationary blades 09, and this structure may exist in engines and gas turbines. The optical path system is mainly composed of a pulsed laser 01, a mirror 13, an air wedge 15, and a beam splitter 14, which are used to construct a staggered laser grid in the area to be measured in the flow field to realize the excitation and marking of tracer molecules in the flow field to be measured. The image capture system is mainly composed of a synchronous controller 10, a first enhanced high-speed camera 04 and a second enhanced high-speed camera 07, which can obtain two consecutive phosphorescence grid images of traced molecules at different times. The image processing system 12 mainly reconstructs the three-dimensional model of the phosphorescence grid based on the phosphorescence grid image of the tracer molecules, and develops a cross-correlation algorithm to obtain the relationship between the spatial positions of the grid at different times according to the distribution characteristics of the phosphorescence intensity, thereby obtaining the three-dimensional velocity field.

动静叶栅结构可存在发动机燃气轮机当中,采用石英玻璃对待测量区域进行可视化,同时为了避免拍摄过程其他荧光对示踪分子磷光的影响,应该对叶片进行氧化发黑处理。采用示踪分子为联乙酰、或者丙酮,此时流场中的气体应该无氧存在。氧气会造成示踪分子磷光淬灭。The dynamic and static blade cascade structure can exist in the engine gas turbine, and the area to be measured is visualized by using quartz glass. At the same time, in order to avoid the influence of other fluorescence on the phosphorescence of the tracer molecules during the shooting process, the blades should be oxidized and blackened. If the tracer molecule is biacetyl or acetone, the gas in the flow field should be free of oxygen. Oxygen quenches tracer molecule phosphorescence.

参见图2,图2为光路系统的结构图,光路系统中一束激光通过空气劈尖15分成多个激光薄片,随后通过分光镜14、反光镜13将一束激光分成两束激光并且在待测流场中形成交错激光网格。利用光的干涉搭建高质量的激光网格,其中示踪分子受激发出磷光和激光强度成正比,根据激光强度特性,开发磷光图像重构算法,可将速度场精度提高的亚像素精度。Referring to Fig. 2, Fig. 2 is a structural diagram of the optical path system. In the optical path system, a beam of laser light is divided into a plurality of laser slices by an air wedge 15, and then a beam of laser light is divided into two beams of laser light by a beam splitter 14 and a reflector 13, and then the beam of laser light is divided into two laser beams in the waiting state. A staggered laser grid is formed in the flow field. Using the interference of light to build a high-quality laser grid, in which the phosphorescence of the tracer molecules is proportional to the laser intensity. According to the characteristics of the laser intensity, the phosphorescence image reconstruction algorithm can be developed to improve the accuracy of the velocity field to sub-pixel accuracy.

参见图3,图3为空气劈尖的结构图,空气劈尖15主要由两块石英玻璃和金属丝组成,两个石英玻璃之间形成一个空气薄膜,利用光的干涉将一束激光分成多个薄片。其中,第一石英玻璃16上表面激光可全穿透,下表面一部分发生发射,一部分发生折射;第二石英玻璃18上表面发生全发射。Referring to Fig. 3, Fig. 3 is the structural diagram of the air splitter, the air splitter 15 is mainly composed of two pieces of quartz glass and a metal wire, an air film is formed between the two quartz glasses, and a beam of laser light is divided into multiple parts by interference of light. slices. Wherein, the upper surface of the first quartz glass 16 can fully penetrate the laser light, a part of the lower surface is emitted, and a part is refracted; the upper surface of the second quartz glass 18 is fully emitted.

激光在空气薄膜上下表面发生干涉,形成明暗相间的激光薄片。由于金属丝直径D很小,两相干光相遇时光程差约为

Figure GDA0004001803650000071
形成明暗条纹的条件为:The laser light interferes on the upper and lower surfaces of the air film, forming light and dark laser flakes. Since the diameter D of the metal wire is very small, the path difference when two coherent lights meet is about
Figure GDA0004001803650000071
The conditions for forming light and dark stripes are:

Figure GDA0004001803650000072
Figure GDA0004001803650000072

相邻两条明纹间空气薄膜的厚度差

Figure GDA0004001803650000073
可知相邻两条明纹之间的距离为l=Δdsinθ,由于空气薄膜的角度从θ很小,近似为/>
Figure GDA0004001803650000074
金属丝直径D和相邻两条明纹间距l之间的关系为/>
Figure GDA0004001803650000075
Thickness difference of air film between two adjacent bright grains
Figure GDA0004001803650000073
It can be seen that the distance between two adjacent bright lines is l=Δdsinθ, since the angle of the air film is very small from θ, it is approximated as />
Figure GDA0004001803650000074
The relationship between the wire diameter D and the distance l between two adjacent open grains is />
Figure GDA0004001803650000075

当动静叶栅间距离为10mm,待测流场空间小,对激光网格质量要求高。选用联乙酰作为示踪份子,选用激光波长355nm,L为20mm,两条明纹之间的距离为1mm,此时金属丝直径为3.55×10-6m。利用6条明纹激光在待测流场形成一个6×6激光标记网格,从而实现待测流场的标记和激发。When the distance between the dynamic and static cascades is 10mm, the space of the flow field to be measured is small, and the quality of the laser grid is required to be high. Diacetyl was selected as the tracer molecule, the laser wavelength was 355nm, L was 20mm, the distance between two bright lines was 1mm, and the diameter of the metal wire was 3.55×10 -6 m. A 6×6 laser marking grid is formed in the flow field to be measured by using 6 bright stripe lasers, so as to realize the marking and excitation of the flow field to be measured.

图像捕捉系统中增强型高速相机由高速摄像机、图像增强器、镜头以及滤波片等组成,可实现磷光信号的精准捕捉,测试中需要根据示踪分子磷光特性的不同选择不同尺寸的滤波片。同步控制器10通过第三数据传输线11连接有图像处理系统12,同步控制器10通过第一数据传输线03连接第一增强型高速相机04,同步控制器10通过第二数据传输线08连接第二增强型高速相机07。同时,利用同步控制器10触发控制第一增强型高速相机04、第二增强型高速相机07以及脉冲激光器01,保证两个相机拍摄的每一帧图像是同步的,从而实现脉冲激光器和增强型高速相机时序控制。图像捕捉系统中可在示踪分子的磷光寿命内,连续获得两个不同时刻t1和t2流场中的磷光网格图形。其中,连续拍摄的时间由相机频率决定

Figure GDA0004001803650000081
为了避免求解速度场时产生错误的信息,相邻图像中网格节点的最大位移应保持在网格间距的一半以下/>
Figure GDA0004001803650000082
此时,相机的频率为/>
Figure GDA0004001803650000083
且由于拍摄流场为高速流场,相机每一帧图像的曝光时间不宜太长,否则将会造成较大的拖影,对速度场的求解精度有较大的影响。为了保证测量精度,曝光时间Δt应满足/>
Figure GDA0004001803650000084
The enhanced high-speed camera in the image capture system consists of a high-speed camera, an image intensifier, a lens, and a filter, which can accurately capture phosphorescence signals. During the test, filters of different sizes need to be selected according to the phosphorescence characteristics of the tracer molecules. The synchronous controller 10 is connected with the image processing system 12 through the third data transmission line 11, the synchronous controller 10 is connected with the first enhanced high-speed camera 04 through the first data transmission line 03, and the synchronous controller 10 is connected with the second enhanced high-speed camera 04 through the second data transmission line 08. Type High Speed Camera 07. Simultaneously, use synchronous controller 10 to trigger and control the first enhanced high-speed camera 04, the second enhanced high-speed camera 07 and the pulsed laser 01 to ensure that each frame of images captured by the two cameras is synchronized, thereby realizing the pulsed laser and enhanced High-speed camera timing control. In the image capture system, within the phosphorescence lifetime of the tracer molecule, the phosphorescence grid patterns in the flow field at two different times t 1 and t 2 can be obtained continuously. Among them, the continuous shooting time is determined by the camera frequency
Figure GDA0004001803650000081
To avoid erroneous information when solving the velocity field, the maximum displacement of grid nodes in adjacent images should be kept below half the grid spacing />
Figure GDA0004001803650000082
At this point, the frequency of the camera is />
Figure GDA0004001803650000083
And because the shooting flow field is a high-speed flow field, the exposure time of each frame of the camera image should not be too long, otherwise it will cause a large smear, which will have a great impact on the solution accuracy of the velocity field. In order to ensure the measurement accuracy, the exposure time Δt should satisfy />
Figure GDA0004001803650000084

图像处理部分主要利用图像捕捉系统得到的示踪分子磷光图像,并结合理论分析与图像处理技术,得到叶栅间三维速度场。在联乙酰磷光寿命内,利用第一增强型高速相机04和第二增强型高速相机07分别获得两个连续时刻的磷光图像,得到的磷光网格图像如图4所示分别表示为:

Figure GDA0004001803650000091
获得包含以下步骤:The image processing part mainly uses the tracer molecular phosphorescence image obtained by the image capture system, and combines theoretical analysis and image processing technology to obtain the three-dimensional velocity field between the cascades. Within the biacetyl phosphorescence lifetime, use the first enhanced high-speed camera 04 and the second enhanced high-speed camera 07 to obtain phosphorescence images at two consecutive moments, and the obtained phosphorescence grid images are shown in Figure 4 as follows:
Figure GDA0004001803650000091
Obtaining involves the following steps:

步骤一、进行标定实验,计算两相机与实际空间位置关系,从而得到图像像素坐标和实际空间坐标之间的转换关系式:Step 1. Carry out a calibration experiment, and calculate the relationship between the two cameras and the actual space position, so as to obtain the conversion relationship between the image pixel coordinates and the actual space coordinates:

Figure GDA0004001803650000092
Figure GDA0004001803650000092

式中:(u,v)代表像素坐标;(Xw,Yw,Zw)代表空间实际坐标;kij代表求解参数,i=1,2,3,j=1,2,3,4);In the formula: (u, v) represents the pixel coordinates; (X w , Y w , Z w ) represents the actual space coordinates; k ij represents the solution parameters, i=1,2,3, j=1,2,3,4 );

步骤二、对得到的磷光图像进行预处理,采用高斯滤波器、均值滤波器对图像进行去噪处理,去除背景噪声,利用图像处理软件,对图像进行统一处理,增强图像对比度,提高磷光图像亮度;Step 2, preprocessing the obtained phosphorescence image, using Gaussian filter and mean value filter to denoise the image to remove background noise, using image processing software to uniformly process the image, enhancing the image contrast, and improving the brightness of the phosphorescence image ;

步骤三、对于单条磷光网格图像,在垂直于网格线的方向采用高斯函数对其磷光光束进行拟合,从而提高图像空间分辨率,将图像精度提高到亚像素;Step 3. For a single phosphorescent grid image, a Gaussian function is used to fit the phosphorescent beam in the direction perpendicular to the grid lines, thereby improving the spatial resolution of the image and improving the image accuracy to sub-pixels;

步骤四、对于t1时刻两个相机获得的磷光图像

Figure GDA0004001803650000093
根据磷光最大强度点进行网格交点一一匹配,得到实际空间中同一个网格交点在两个不同相机拍摄图片中像素点的位置,随后将两个像素点坐标带入式2,计算出实际网格空间点的坐标(Xw,Yw,Zw)t1;Step 4. Phosphorescence images obtained by two cameras at time t 1
Figure GDA0004001803650000093
Match the grid intersection points one by one according to the phosphorescence maximum intensity point, and obtain the pixel position of the same grid intersection point in the actual space in the pictures taken by two different cameras, and then bring the coordinates of the two pixel points into Equation 2 to calculate the actual Coordinates (X w , Y w , Z w ) t1 of grid space points;

步骤五、分别对相机1和相机2获得两个时刻的磷光图像

Figure GDA0004001803650000094
和/>
Figure GDA0004001803650000095
和/>
Figure GDA0004001803650000096
进行互相关计算,得到同一个相机捕捉到两个不同时刻磷光网格图片中网格交点的一一对应关系。利用互相关算法计算得到t2时刻两个相机磷光图像(/>
Figure GDA0004001803650000097
和/>
Figure GDA0004001803650000098
)中网格交点的像素点的坐标,带入式2,计算出t2时刻网格交点在实际空间中的坐标(Xw,Yw,Zw)t2;Step 5: Obtain phosphorescence images at two moments for camera 1 and camera 2 respectively
Figure GDA0004001803650000094
and />
Figure GDA0004001803650000095
and />
Figure GDA0004001803650000096
Perform cross-correlation calculations to obtain the one-to-one correspondence between grid intersections in two phosphorescent grid pictures captured by the same camera at different times. Using the cross-correlation algorithm to calculate the phosphorescence images of the two cameras at time t 2 (/>
Figure GDA0004001803650000097
and />
Figure GDA0004001803650000098
), bring the coordinates of the pixel points of the grid intersection point into formula 2, and calculate the coordinates (X w , Y w , Z w ) t2 of the grid intersection point in the actual space at time t2 ;

步骤六、根据t1时刻和t2时刻网格交点在实际空间中的位置(Xw,Yw,Zw)t1、(Xw,Yw,Zw)t2与时间的关系,从而计算得到流场的三维速度场为

Figure GDA0004001803650000101
Figure GDA0004001803650000102
Step 6. According to the relationship between the position (X w , Y w , Z w ) t1 and (X w , Y w , Z w ) t2 of the grid intersection in the actual space at time t 1 and time t 2 and the relationship between time and calculation The three-dimensional velocity field of the flow field is obtained as
Figure GDA0004001803650000101
Figure GDA0004001803650000102

进一步,由于实际制造加工过程中,因为工艺的差异,会造成动静叶栅间流场的不均匀,为了测量动静叶栅间三维速度场的大小以及流场的不均匀程度,可通过对脉冲激光器、增强型高速相机的时序控制对每个动叶片经过同一个静叶片时的速度场进行测量,通过比较不同动叶片之间的速度场的差异,从而实现流场的周期内均匀性的测量。速度场不均匀性测量可描述为:Furthermore, due to the difference in the actual manufacturing process, the flow field between the static and dynamic cascades will be uneven due to the difference in the process. In order to measure the size of the three-dimensional velocity field and the unevenness of the flow field between the dynamic and static cascades, the , The timing control of the enhanced high-speed camera measures the velocity field of each moving blade when it passes the same stationary blade, and by comparing the difference in the velocity field between different moving blades, the measurement of the uniformity of the flow field within a period is realized. Velocity field inhomogeneity measurement can be described as:

当动叶片运动的转速为n,动叶片个数为k,叶片旋转一周的时间为

Figure GDA0004001803650000103
此时同一个静叶片需经历2个相邻的叶片所需要的时间为/>
Figure GDA0004001803650000104
此时脉冲激光器的频率为f=nk,高速相机拍摄时间/>
Figure GDA0004001803650000105
由于示踪分子的磷光寿命可达1ms,为了保证同一时间待测流场中不会出现两组磷光网格图像,相邻两次脉冲激光激发时间应满足:/>
Figure GDA0004001803650000106
Figure GDA0004001803650000107
时,为了保证待测流场中不会出现两组磷光网格图像,脉冲激光器的频率为
Figure GDA0004001803650000108
高速摄像机拍摄时间/>
Figure GDA0004001803650000109
虽然脉冲激光器频率减小,但高速相机拍摄也随之变周期增大,将导致相机的内存增大。同时,为了保证测量的准确性,初始时刻,磷光网格图像初始分布应布置在两个动叶片中间位置。When the rotating speed of the moving blade is n, the number of moving blades is k, and the time for the blade to rotate one circle is
Figure GDA0004001803650000103
At this time, the time required for the same stationary blade to experience two adjacent blades is />
Figure GDA0004001803650000104
At this time, the frequency of the pulsed laser is f=nk, and the shooting time of the high-speed camera />
Figure GDA0004001803650000105
Since the phosphorescence lifetime of tracer molecules can reach 1 ms, in order to ensure that two sets of phosphorescence grid images do not appear in the flow field to be measured at the same time, the excitation time of two adjacent pulse lasers should meet: />
Figure GDA0004001803650000106
when
Figure GDA0004001803650000107
, in order to ensure that two groups of phosphorescent grid images do not appear in the flow field to be measured, the frequency of the pulsed laser is
Figure GDA0004001803650000108
High-speed camera shooting time/>
Figure GDA0004001803650000109
Although the frequency of the pulsed laser decreases, the shooting cycle of the high-speed camera also increases, which will lead to an increase in the memory of the camera. At the same time, in order to ensure the accuracy of the measurement, the initial distribution of the phosphorescent grid image should be arranged in the middle of the two moving blades at the initial moment.

一种动静叶栅间速度场均匀性测量流程如下:A measurement process of velocity field uniformity between static and dynamic cascades is as follows:

1)根据待测流场的特征,选择示踪分子种类,确定激发激光波长,空气劈尖的结构,滤波片尺寸;1) According to the characteristics of the flow field to be measured, select the type of tracer molecules, determine the wavelength of the excitation laser, the structure of the air wedge, and the size of the filter;

2)根据动叶片的转速、叶片数,确定脉冲激光器频率和拍摄周期;2) Determine the pulse laser frequency and shooting period according to the rotating speed and number of blades;

3)根据叶片的转速和叶片直径确定估算流场速度Vest,max,计算相机频率和曝光时间;3) Determine the estimated flow field velocity V est,max according to the blade speed and blade diameter, and calculate the camera frequency and exposure time;

4)设置好各个设备参数,开展相机标定实验,进行动静叶栅间速度场测量实验,同时计算流场速度Vact4) Set the parameters of each equipment, carry out the camera calibration experiment, conduct the velocity field measurement experiment between the dynamic and static cascades, and calculate the flow field velocity V act at the same time;

5)比较Vest,max和Vact,max的大小,若Vest,max>Vact,max,试验结束;若Vest,max<Vact,max,则需要根据Vact,max重新计算相机频率和曝光时间,重复实验,直到结束。5) Compare V est,max and V act,max , if V est,max >V act,max , the test is over; if V est,max <V act,max , you need to recalculate the camera according to V act,max Frequency and exposure time, repeat the experiment until the end.

以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。The above content is only to illustrate the technical ideas of the present invention, and cannot limit the protection scope of the present invention. Any changes made on the basis of the technical solutions according to the technical ideas proposed in the present invention shall fall within the scope of the claims of the present invention. within the scope of protection.

Claims (6)

1.一种动静叶栅间三维速度场均匀性测量方法,其特征在于,基于动静叶栅间三维速度场均匀性测量装置进行,所述动静叶栅间三维速度场均匀性测量装置,包括脉冲激光器(01)和光路系统(05),脉冲激光器(01)用于发射激光(02),所述激光(02)进入光路系统(05)内改变光路,之后经光路系统(05)的出口变成两束薄片状激光,所述激光薄片进入动叶片(06)和静叶片(09)之间待测流场中形成交错激光网格;1. A method for measuring the uniformity of the three-dimensional velocity field between the static and dynamic cascades, characterized in that, based on a measuring device for the uniformity of the three-dimensional velocity field between the static and dynamic cascades, the measuring device for the uniformity of the three-dimensional velocity field between the static and dynamic cascades includes a pulse The laser (01) and the optical path system (05), the pulse laser (01) is used to emit laser light (02), and the laser light (02) enters the optical path system (05) to change the optical path, and then changes the optical path through the exit of the optical path system (05). into two beams of sheet-shaped lasers, and the laser sheets enter the flow field to be measured between the moving blade (06) and the stationary blade (09) to form a staggered laser grid; 所述待测流场中采用联乙酰或丙酮作为示踪分子,待测流场为无氧环境;In the flow field to be measured, biacetyl or acetone is used as a tracer molecule, and the flow field to be measured is an oxygen-free environment; 待测流场的两侧分别设有第一增强型高速相机(04)和第二增强型高速相机(07),第一增强型高速相机(04)和第二增强型高速相机(07)分别通过连接线连接有同步控制器(10),同步控制器(10)连接有图像处理系统(12);A first enhanced high-speed camera (04) and a second enhanced high-speed camera (07) are respectively arranged on both sides of the flow field to be measured, and the first enhanced high-speed camera (04) and the second enhanced high-speed camera (07) are respectively A synchronous controller (10) is connected through a connection line, and the synchronous controller (10) is connected with an image processing system (12); 所述第一增强型高速相机(04)和第二增强型高速相机(07)用于同步的拍摄示踪分子受激发后产生的磷光图像,包括以下步骤:The first enhanced high-speed camera (04) and the second enhanced high-speed camera (07) are used for synchronous shooting of phosphorescent images produced after the tracer molecules are excited, including the following steps: 步骤一、进行标定实验,计算第一增强型高速相机(04)、第二增强型高速相机(07)与实际空间位置关系,从而得到图像像素坐标和实际空间坐标之间的转换关系式:Step 1. Carry out a calibration experiment, and calculate the relationship between the first enhanced high-speed camera (04), the second enhanced high-speed camera (07) and the actual space position, thereby obtaining the conversion relation between the image pixel coordinates and the actual space coordinates:
Figure FDA0004001803640000011
Figure FDA0004001803640000011
式中:(u,v)代表像素坐标;(Xw,Yw,Zw)代表空间实际坐标;kij代表求解参数,i=1,2,3,j=1,2,3,4;In the formula: (u, v) represents the pixel coordinates; (X w , Y w , Z w ) represents the actual space coordinates; k ij represents the solution parameters, i=1,2,3, j=1,2,3,4 ; 步骤二、在示踪分子的磷光寿命内,利用第一增强型高速相机(04)和第二增强型高速相机(07)获得两个连续时刻t1、t2的磷光图像,分别表示为:
Figure FDA0004001803640000012
Figure FDA0004001803640000013
Step 2. Within the phosphorescence lifetime of the tracer molecule, use the first enhanced high-speed camera (04) and the second enhanced high-speed camera (07) to obtain two consecutive phosphorescence images at time t 1 and t 2 , respectively expressed as:
Figure FDA0004001803640000012
Figure FDA0004001803640000013
步骤三、对于t1时刻第一增强型高速相机(04)、第二增强型高速相机(07)获得的磷光图像
Figure FDA0004001803640000021
和/>
Figure FDA0004001803640000022
根据磷光最大强度点进行网格交点一一匹配,得到实际空间中同一个网格交点在两个不同相机拍摄图片中像素点的坐标,随后将两个像素点坐标带入式(2),计算出实际网格空间点的坐标(Xw,Yw,Zw)t1
Step 3, for the phosphorescence images obtained by the first enhanced high-speed camera (04) and the second enhanced high-speed camera (07) at time t1
Figure FDA0004001803640000021
and />
Figure FDA0004001803640000022
According to the maximum intensity point of phosphorescence, the grid intersection points are matched one by one, and the pixel coordinates of the same grid intersection point in the actual space in the pictures taken by two different cameras are obtained, and then the coordinates of the two pixel points are brought into formula (2) to calculate coordinates (X w , Y w , Z w ) t1 of the actual grid space point;
步骤四、分别对第一增强型高速相机(04)、第二增强型高速相机(07)获得两个时刻的磷光图像
Figure FDA0004001803640000023
和/>
Figure FDA0004001803640000024
Figure FDA0004001803640000025
和/>
Figure FDA0004001803640000027
进行互相关计算,得到同一个相机捕捉到两个不同时刻磷光网格图片中网格交点的对应关系;利用互相关算法计算得到t2时刻两个相机磷光图像/>
Figure FDA0004001803640000028
和/>
Figure FDA0004001803640000029
中网格交点的像素点的坐标,带入式(2),计算出t2时刻网格交点在实际空间中的坐标(Xw,Yw,Zw)t2
Step 4. Obtain phosphorescence images at two moments for the first enhanced high-speed camera (04) and the second enhanced high-speed camera (07) respectively
Figure FDA0004001803640000023
and />
Figure FDA0004001803640000024
Figure FDA0004001803640000025
and />
Figure FDA0004001803640000027
Carry out cross-correlation calculations to obtain the corresponding relationship between grid intersection points in two phosphorescent grid pictures captured by the same camera at different times; use cross-correlation algorithm to calculate phosphorescence images of two cameras at time t2 >
Figure FDA0004001803640000028
and />
Figure FDA0004001803640000029
The coordinates of the pixel points of the grid intersection in the middle grid are brought into formula (2), and the coordinates (X w , Y w , Z w ) t2 of the grid intersection in the actual space at time t2 are calculated;
步骤五、根据t1时刻和t2时刻网格交点在实际空间中的位置(Xw,Yw,Zw)t1、(Xw,Yw,Zw)t2与时间的关系,计算得到流场的三维速度场为
Figure FDA00040018036400000210
Step 5. According to the relationship between the position (X w , Y w , Z w ) t1 and (X w , Y w , Z w ) t2 of the grid intersection point in the actual space at time t 1 and time t 2 and the time, calculate The three-dimensional velocity field of the flow field is
Figure FDA00040018036400000210
通过对脉冲激光器(01)、第一增强型高速相机(04)和第二增强型高速相机(07)的时序控制,使得对每个动叶片(06)经过同一个静叶片(09)时的速度场进行测量,通过比较不同动叶片(06)之间的速度场的差异,从而实现流场的周期内均匀性的测量。Through the timing control of the pulse laser (01), the first enhanced high-speed camera (04) and the second enhanced high-speed camera (07), the time when each moving blade (06) passes the same stationary blade (09) The velocity field is measured, and the measurement of the uniformity of the flow field within a cycle is realized by comparing the difference of the velocity field between different moving blades (06).
2.根据权利要求1所述的动静叶栅间三维速度场均匀性测量方法,其特征在于,光路系统(05)包括反光镜(13)、分光镜(14)和空气劈尖(15);2. The method for measuring the uniformity of the three-dimensional velocity field between moving and static cascades according to claim 1, wherein the optical path system (05) includes a mirror (13), a beam splitter (14) and an air wedge (15); 当一束激光通过空气劈尖(15)后,分成多个薄片状激光,随后通过分光镜(14)、反光镜(13)将一束激光分成两束激光,两束激光进入待测流场中形成交错激光网格。When a beam of laser light passes through the air wedge (15), it is divided into multiple lamellar laser beams, and then the laser beam is divided into two laser beams by the beam splitter (14) and mirror (13), and the two laser beams enter the flow field to be measured A staggered laser grid is formed in the 3.根据权利要求1所述的动静叶栅间三维速度场均匀性测量方法,其特征在于,空气劈尖(15)包括第一石英玻璃(16)和第二石英玻璃(18),第一石英玻璃(16)和第二石英玻璃(18)一端相接触,另一端之间设有金属丝(17);3. The method for measuring the uniformity of the three-dimensional velocity field between moving and static cascades according to claim 1, wherein the air wedge (15) comprises a first quartz glass (16) and a second quartz glass (18), the first One end of the quartz glass (16) is in contact with the second quartz glass (18), and a metal wire (17) is arranged between the other ends; 两个石英玻璃之间形成一个空气薄膜,所述空气薄膜利用光的干涉将一束激光分成多个薄片状;An air film is formed between the two quartz glasses, and the air film uses light interference to divide a beam of laser light into multiple thin slices; 当激光入射时,激光能够在第一石英玻璃(16)的上表面激光全穿透,在下表面一部分发生发射,另一部分发生折射,在第二石英玻璃(18)上表面发生全发射。When the laser is incident, the laser can fully penetrate the upper surface of the first quartz glass (16), partly emit on the lower surface, refract the other part, and fully emit on the upper surface of the second quartz glass (18). 4.根据权利要求1所述的动静叶栅间三维速度场均匀性测量方法,其特征在于,步骤二还包括:4. The method for measuring the uniformity of the three-dimensional velocity field between the static and dynamic cascades according to claim 1, wherein step 2 further comprises: 对所述磷光图像进行去噪处理,同时增强图像对比度;performing denoising processing on the phosphorescence image while enhancing image contrast; 对于单条磷光网格图像,在垂直于网格线的方向采用高斯函数对磷光光束进行拟合。For a single phosphorescent grid image, a Gaussian function is used to fit the phosphorescent beam in the direction perpendicular to the grid lines. 5.根据权利要求4所述的动静叶栅间三维速度场均匀性测量方法,其特征在于,测量速度场均匀性具体操作为:5. The method for measuring the uniformity of the three-dimensional velocity field between the dynamic and static cascades according to claim 4, wherein the specific operation of measuring the uniformity of the velocity field is as follows: 当动叶片(06)运动的转速为n,动叶片(06)个数为k,则动叶片(06)旋转一周的时间为
Figure FDA0004001803640000031
此时同一个静叶片(09)经历2个相邻的叶片所需要的时间为/>
Figure FDA0004001803640000032
此时脉冲激光器的频率为f=nk,高速相机拍摄时间/>
Figure FDA0004001803640000033
When the rotating speed of the moving blade (06) is n and the number of moving blades (06) is k, the time for the moving blade (06) to rotate one circle is
Figure FDA0004001803640000031
At this time, the time required for the same stationary blade (09) to experience two adjacent blades is />
Figure FDA0004001803640000032
At this time, the frequency of the pulsed laser is f=nk, and the shooting time of the high-speed camera />
Figure FDA0004001803640000033
相邻两次脉冲激光激发时间满足:
Figure FDA0004001803640000034
The excitation time of two adjacent pulse lasers satisfies:
Figure FDA0004001803640000034
Figure FDA0004001803640000035
时,脉冲激光器的频率为/>
Figure FDA0004001803640000036
高速摄像机拍摄时间/>
Figure FDA0004001803640000037
when
Figure FDA0004001803640000035
, the frequency of the pulsed laser is />
Figure FDA0004001803640000036
High-speed camera shooting time/>
Figure FDA0004001803640000037
6.根据权利要求4所述的动静叶栅间三维速度场均匀性测量方法,其特征在于,在t1时刻,磷光网格图像初始分布布置在两个动叶片中间位置。6. The method for measuring the uniformity of the three-dimensional velocity field between the moving and static cascades according to claim 4, characterized in that, at time t1 , the initial distribution of the phosphorescence grid image is arranged in the middle of the two moving blades.
CN202210255200.6A 2022-03-15 2022-03-15 Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids Active CN114740222B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210255200.6A CN114740222B (en) 2022-03-15 2022-03-15 Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210255200.6A CN114740222B (en) 2022-03-15 2022-03-15 Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids

Publications (2)

Publication Number Publication Date
CN114740222A CN114740222A (en) 2022-07-12
CN114740222B true CN114740222B (en) 2023-03-28

Family

ID=82276284

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210255200.6A Active CN114740222B (en) 2022-03-15 2022-03-15 Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids

Country Status (1)

Country Link
CN (1) CN114740222B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115469115B (en) * 2022-11-14 2023-01-31 中国空气动力研究与发展中心设备设计与测试技术研究所 Fluorescence detection method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03209169A (en) * 1990-01-11 1991-09-12 Mitsubishi Heavy Ind Ltd Method for measuring velocity of flow using laser induced fluorescent method
CN105242066A (en) * 2015-09-30 2016-01-13 东南大学 Synchronous measuring apparatus and method for wall jet gas-solid two-phase velocity fields
JP2019132840A (en) * 2018-01-29 2019-08-08 国立大学法人 東京大学 Three-dimensional fluid measuring method
CN111458532A (en) * 2020-04-17 2020-07-28 东北电力大学 System and method for synchronously measuring speed field and density field
CN111579248A (en) * 2020-05-14 2020-08-25 中国空气动力研究与发展中心 Optical synchronous measurement method for combustion flow field of scramjet engine

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3281755B2 (en) * 1995-04-05 2002-05-13 三菱重工業株式会社 Measurement device by laser induced fluorescence method
US20060146910A1 (en) * 2004-11-23 2006-07-06 Manoochehr Koochesfahani Method and apparatus for simultaneous velocity and temperature measurements in fluid flow
CN104807501A (en) * 2015-04-15 2015-07-29 上海交通大学 UV laser-induced phosphorescence imaging type synchronous speed and temperature measuring method and system
CN105973852B (en) * 2016-03-16 2019-01-11 上海交通大学 Fuel jet distribution of concentration test device and its implementation method
CN107478267B (en) * 2017-07-11 2020-02-21 上海交通大学 Synchronization test method of temperature field and velocity field of 3D flow field based on light field camera
CN108362906B (en) * 2018-02-12 2021-06-08 西北核技术研究所 A molecular marker two-dimensional velocity vector measurement device and method
CN108519495B (en) * 2018-02-28 2020-12-04 哈尔滨工程大学 A device and method for measuring internal flow field velocity of turbomachinery
CN109084914B (en) * 2018-09-17 2020-07-28 上海交通大学 A full-field measurement system and method for high-temperature and high-speed turbulent heat flux
CN111596092A (en) * 2020-05-18 2020-08-28 天津大学 Laser-induced phosphorescence particle-based imaging speed measurement device and method
CN113376401B (en) * 2021-04-28 2022-12-09 西安交通大学 A flow-controllable tracer molecule adding device and adding method thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03209169A (en) * 1990-01-11 1991-09-12 Mitsubishi Heavy Ind Ltd Method for measuring velocity of flow using laser induced fluorescent method
CN105242066A (en) * 2015-09-30 2016-01-13 东南大学 Synchronous measuring apparatus and method for wall jet gas-solid two-phase velocity fields
JP2019132840A (en) * 2018-01-29 2019-08-08 国立大学法人 東京大学 Three-dimensional fluid measuring method
CN111458532A (en) * 2020-04-17 2020-07-28 东北电力大学 System and method for synchronously measuring speed field and density field
CN111579248A (en) * 2020-05-14 2020-08-25 中国空气动力研究与发展中心 Optical synchronous measurement method for combustion flow field of scramjet engine

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
激光诱导荧光法用于内燃机燃烧可视化的研究进展;马骁 等;《车用发动机》;20080229(第1期);第7-13页 *

Also Published As

Publication number Publication date
CN114740222A (en) 2022-07-12

Similar Documents

Publication Publication Date Title
CN110220610B (en) A system capable of simultaneously measuring combustion temperature field and velocity field
CN102323445B (en) Method for measuring flowing velocity of recycle gas in multiple laser cavity
CN105004466B (en) A high-precision non-contact pneumatic friction resistance measurement method and measurement device
CN106644353A (en) PIV wind tunnel test method for ship air flow field measurement
JP3718686B2 (en) Plane and space time-series fluid velocity measurement system
CN111458532A (en) System and method for synchronously measuring speed field and density field
CN114740222B (en) Device and method for measuring uniformity of three-dimensional velocity field between moving blade and static blade grids
CN204988587U (en) Pneumatic frictional resistance measuring device of high accuracy non -contact
CN113252295A (en) Flow field measuring instrument and flow field measuring method based on rectangular pyramid sensor
CN107883889A (en) Vibration test 3 D deformation measurement apparatus and method based on laser speckle interferometry
CN108362906B (en) A molecular marker two-dimensional velocity vector measurement device and method
CN114354036B (en) Method and device for synchronously measuring surface pressure and three-dimensional shape of motion model
CN111473943A (en) A PIV near-wall data processing method, device and storage medium
Neuendorf et al. On large streamwise structures in a wall jet flowing over a circular cylinder
CN204269339U (en) A kind of various visual angles projection oldered array deviation imaging system
CN103048112B (en) Device and method for measuring disturbance of cross flow fan impeller of pulsed gas laser
CN114353697A (en) Ablation process structure surface topography measuring method
CN112036072A (en) Three-dimensional tracer particle matching method and velocity field measuring device
Zhang et al. Development of a digital image projection technique to measure wind-driven water film flows
Post et al. Two-color particle-imaging velocimetry using a single argon-ion laser
CN116718344A (en) Multi-parameter thrust vector jet optical calibration method
CN101210935A (en) A Measuring Method of Transient Characteristics of Secondary Flow Field
WO2023023961A1 (en) Piv image calibration apparatus and method based on laser linear array
Orloff et al. High-resolution LDA measurements of Reynolds stress in boundary layers and wakes
Adrian Engineering applications of particle image velocimeters

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