13-hole probe and fluid velocity vector measuring method and system

Abstract

The invention discloses a thirteen-hole probe and a fluid velocity vector measurement method and system, and the method comprises the steps: carrying out the partitioning of the thirteen-hole probe according to the physical position of each pressure hole in the surface of a spherical probe; adjusting the angle attitude of the thirteen-hole probe, obtaining multiple groups of calibration data under a known flow field, and calculating an angle coefficient characteristic curve graph of each subarea by using the calibration data; performing dimensionless processing on the pore pressure data of each pressure pore in the known flow field to obtain a dimensionless pressure coefficient for fitting, and establishing a general calculation model; and collecting actually measured pore pressure data of the thirteen-hole probe in the flow field to be measured, positioning the subarea to which the actually measured pore pressure data belongs, and calling the angle characteristic curve of the corresponding subarea in the general calculation model to calculate the flow field parameters of the flow field to be measured. A plurality of pressure holes are formed in the spherical probe and used for receiving airflow in multiple directions, the defect that only airflow in the axis direction of the probe can be received in the prior art is overcome, and large-angle bias flow measurement is achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of probe measurement, in particular to a thirteen-hole probe and a fluid velocity vector measurement method and system. BACKGROUND

[0002] In the field of fluid velocity measurement, the fluid such as air flow or water flow is usually three-dimensional unsteady flow, and the velocity and direction of the fluid usually have a certain influence on the performance of the equipment placed in the fluid. In some more complex flow state researches, such as turbulent flow, separated flow, multiphase flow, etc., it is more necessary to analyze the flow structure and evolution law through vector data.

[0003] At present, the measurement technology of air flow is increasingly mature. The pitot tube can capture the total pressure of the air flow through the front total pressure hole, capture the ambient static pressure through the side wall static pressure hole, and then calculate the air flow velocity by using the Bernoulli equation through the difference between the total pressure and the static pressure, or use the hot-wire anemometer to use the characteristics of the air flow taking away heat, and use the change of the hot-wire resistance to deduce the air flow velocity.

[0004] However, the pitot tube can only measure one-dimensional velocity along the probe axis direction, and the hot-wire anemometer is easy to be polluted by dust, oil stains and water droplets in the air flow, resulting in a decrease in measurement accuracy or even failure. Therefore, both the pitot tube and the hot-wire anemometer have defects in the process of air flow velocity vector measurement. SUMMARY

[0005] Therefore, the present application aims to provide a thirteen-hole probe and a fluid velocity vector measurement method and system, which aims to solve the problems of limitations and easy environmental influence in the process of measuring air flow velocity by using the pitot tube or the hot-wire anemometer.

[0006] To achieve the above-mentioned purpose, the present application provides a fluid velocity vector measurement method, which is applied to a thirteen-hole probe, and the method comprises the following steps: According to the physical positions of the pressure holes on the surface of the spherical probe, the thirteen-hole probe is divided into nine partition regions, wherein each partition region is composed of at least three pressure holes; The thirteen-hole probe is subjected to a calibration experiment, the angle posture of the thirteen-hole probe is adjusted, a plurality of sets of calibration data under a known flow field are obtained, and the angle coefficient characteristic curve of each partition region is calculated by using the calibration data; The hole pressure data of each pressure hole under the known flow field are subjected to dimensionless processing, and the dimensionless pressure coefficient is obtained, the dimensionless pressure coefficient is fitted, and a general calculation model is established; Collecting the measured orifice pressure data of the thirteen-hole probe under the flow field to be measured, positioning the partition to which the measured orifice pressure data belongs, and calling the angle characteristic curve of the corresponding partition in the general calculation model to calculate the flow field parameters of the flow field to be measured.

[0007] According to an aspect of the above technical solution, the step of partitioning the thirteen-hole probe is: The first zone is composed of the first, second, third, fourth and fifth holes; the second zone is composed of the second, third and sixth holes; the third zone is composed of the second, fourth and seventh holes; the fourth zone is composed of the third, fifth and eighth holes; the fifth zone is composed of the fourth, fifth and ninth holes; the sixth zone is composed of the second, sixth, seventh and tenth holes; the seventh zone is composed of the third, sixth, eighth and eleventh holes; the eighth zone is composed of the fourth, seventh, ninth and twelfth holes; and the ninth zone is composed of the fifth, eighth, ninth and thirteenth holes.

[0008] According to an aspect of the above technical solution, the step of calibrating the thirteen-hole probe, adjusting the angle posture of the thirteen-hole probe, and obtaining multiple sets of calibration data under the known flow field is: Calibrating the thirteen-hole probe under the known flow field, adjusting the posture of the thirteen-hole probe under different pitch angles and different azimuth angles using a double-turntable, and obtaining calibration data of the thirteen-hole probe under each different posture; After obtaining the calibration data, partitioning the orifice pressure data according to a preset rule, matching the calibration data in the nine partition regions, and using the matched calibration data to calculate the angle coefficient characteristic curve of each partition.

[0009] According to an aspect of the above technical solution, in the step of non-dimensionalizing the orifice pressure data of each pressure hole under the known flow field and obtaining the non-dimensional pressure coefficient, fitting the non-dimensional pressure coefficient, and establishing a general calculation model, From the calibration data, the orifice pressure data of each pressure hole is selected, and after obtaining the orifice pressure data of each pressure hole, the orifice pressure data is non-dimensionalized and converted into a non-dimensional pressure coefficient, wherein the non-dimensional pressure coefficient includes an angle coefficient, a total pressure coefficient and a static pressure coefficient. The non-dimensional pressure coefficient is fitted using a least squares fitting method to establish a general calculation model:

[0010] wherein, represents the pitch angle θ and the azimuth angle , the total pressure coefficient K O and the static pressure coefficient K QAny one of the four dimensionless pressure coefficients, subscript i represents the i-th measuring point; K represents the calibration coefficient, the superscript A represents the calibration coefficient corresponding to all a certain sequence of K, and the subscript represents the corresponding term in the expansion; wherein any one of the dimensionless pressure coefficients in any partition region comprises 15 calibration coefficients.

[0011] According to an aspect of the above technical solution, after establishing the general calculation model, for any one partition region, the dimensionless pressure coefficients of the calibration points of the thirteen-hole probe in different angle postures in the partition region are input into the general calculation model, and a calibration matrix is obtained:

[0012] The simplified form is:

[0013]

[0014] Wherein, the A matrix contains the dimensionless pressure coefficients of n calibration points, the K matrix contains the calibration coefficients of n calibration points, K A T The transpose matrix of K A is K A , and K is a matrix composed of calibration coefficients.

[0015] The application also provides a fluid velocity vector measurement system for implementing the fluid velocity vector measurement method, and the system comprises: A region partition module is configured to partition the thirteen-hole probe according to the physical positions of the pressure holes in the surface of the spherical probe, and divide the surface of the spherical probe into nine partition regions, wherein any one of the partition regions is composed of at least three pressure holes. A curve calculation module is configured to perform a calibration experiment on the thirteen-hole probe, adjust the angle posture of the thirteen-hole probe, obtain a plurality of sets of calibration data under a known flow field, and calculate the angle coefficient characteristic curve of each partition region by using the calibration data. A model establishment module is configured to perform dimensionless processing on the hole pressure data of each pressure hole under the known flow field, obtain dimensionless pressure coefficients, and fit the dimensionless pressure coefficients to establish a general calculation model. A parameter calculation module is configured to collect the measured hole pressure data of the thirteen-hole probe under the to-be-measured flow field, locate the partition to which the measured hole pressure data belongs, and call the angle characteristic curve of the corresponding partition in the general calculation model to calculate the flow field parameters of the to-be-measured flow field.

[0016] The application further provides a thirteen-hole probe, which comprises a spherical probe head, a shell support rod fixed to the spherical probe head, and a plurality of pressure guide channels arranged in the shell support rod. The spherical probe head comprises thirteen pressure holes symmetrically arranged on the spherical probe head body, wherein a first hole is arranged at the center of the spherical probe head and is arranged opposite to the shell support rod, and the remaining pressure holes are evenly arranged on planes at 30°, 90° and 120° with respect to the central axis of the spherical probe head, wherein each layer is provided with four pressure holes, and the pressure holes in each layer are staggered. Among them, the second hole, the third hole, the fourth hole and the fifth hole are arranged on the plane at 30° with respect to the central axis of the spherical probe head, the sixth hole, the seventh hole, the eighth hole and the ninth hole are arranged on the plane at 90° with respect to the central axis of the spherical probe head, and the tenth hole, the eleventh hole, the twelfth hole and the thirteenth hole are arranged on the plane at 120° with respect to the central axis of the spherical probe head.

[0017] The application further provides a computer readable storage medium, which stores a computer program, and the program is executed by a processor to realize the fluid velocity vector measurement method.

[0018] The application further provides an electronic device, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, and the processor realizes the fluid velocity vector measurement method when executing the computer program.

[0019] In summary, according to the fluid velocity vector measurement method provided by the application, the thirteen-hole probe is placed in a known flow field, the attitude of the thirteen-hole probe at different angles is adjusted to obtain calibration data, the angle coefficient characteristic curve of each partition region is calculated by using the calibration data, the hole pressure data of each pressure hole is processed in a dimensionless manner, the correlation between the pressure difference and the flow field parameter is focused on to eliminate the influence of the absolute pressure, and a general calculation model is constructed by using the least square fitting method, so as to be used for subsequent testing of the flow field parameter of the to-be-tested flow field. The plurality of pressure holes arranged on the spherical probe head are used to receive airflow at multiple directions, which overcomes the shortcoming that the traditional technology can only receive airflow along the axis direction of the probe, realizes large-angle deflection measurement, and the general calculation model can be used to directly calculate the flow field parameter of the to-be-tested flow field by measuring the hole pressure data in the pressure hole and focusing on the pressure data of the pressure hole, so that the calculation result is more accurate and reliable.

[0020] Additional aspects and advantages of the application will be made apparent from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0021] Figure 1Structure diagram of the spherical probe in the embodiment one of the present application; Figure 2 Structure diagram of the spherical probe in the embodiment one of the present application from another perspective; Figure 3 Radial sectional view of the thirteen-hole probe in the embodiment one of the present application; Figure 4 Axial sectional view of the thirteen-hole probe in the embodiment one of the present application; Figure 5 Flow chart of the fluid velocity vector measurement method in the embodiment two of the present application; Figure 6 Interpolation calculation example diagram in the embodiment two of the present application; Figure 7 Structure diagram of the fluid velocity vector measurement system in the embodiment three of the present application; Figure 8 Structure block diagram of the electronic device in the embodiment five of the present application.

[0022] Explanation of the components in the drawings: Spherical probe D, shell support B, pressure guide channel C DETAILED DESCRIPTION In order to make the objects, features and advantages of the present application more apparent, the specific embodiments of the present application will be described in detail below with reference to the drawings. The drawings show several embodiments of the present application. However, the present application can be realized in many different forms and is not limited to the embodiments described herein. On the contrary, these embodiments are provided to make the disclosure of the present application more thorough and comprehensive.

[0023] It should be noted that when an element is referred to as being "fixedly attached" to another element, it can be directly on the other element or there can be intervening elements. When an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements can be present. As used herein the terms "vertical", "horizontal", "left", "right", "upper", "lower", and the like as well as the like terms are used for explanation only and not to indicate or imply that specifically oriented operation or construction of devices or elements is required or preferred, and therefore should not be taken to mean limiting of the present application.

[0024] In the present application, unless otherwise explicitly specified and limited, the terms "mounting", "connecting", "connecting", "fixing" and the like should be understood in a broad sense, for example, it can be fixedly connected, or it can be detachably connected, or it can be integrally connected; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium; it can be the communication inside two elements. For those skilled in the art, the specific meaning of the above-mentioned terms in the present application can be understood according to the specific circumstances. The term "and / or" used herein includes any and all combinations of one or more related listed items.

[0025] Embodiment one Please refer to Figures 1-4 , it is a structure diagram of a thirteen-hole probe provided in the embodiment one of the present application, the thirteen-hole probe includes a spherical probe head, a shell support rod fixed on the spherical probe head, and a plurality of pressure guide channels arranged in the shell support rod, wherein: The flow environment in modern engineering and scientific research often contains large-angle deflection, for example, large-angle landing of aircraft, strong vortex in blade passage, separation flow and the like. In order to test the above large-angle deflection, a plurality of pressure holes are symmetrically arranged on the spherical probe head. The present embodiment takes thirteen holes as an example for illustration.

[0026] The spherical probe head includes thirteen pressure holes symmetrically arranged on the spherical probe head body. The arrangement mode of the thirteen pressure holes is 1-4-4-4 structure, which is divided into four layers. The first layer is a first hole, which is arranged at the center of the spherical probe head and opposite to the shell support rod. The remaining pressure holes are evenly arranged on the planes with the center axis of the spherical probe head at 30°, 90° and 120°, and each layer has four pressure holes. The pressure holes in each layer are staggered, and the included angle between the connecting line of any two adjacent pressure holes in any layer and the center of the spherical probe head is 90°.

[0027] It should be noted that the second hole, the third hole, the fourth hole and the fifth hole are arranged on the plane with the center axis of the spherical probe head at 30°, the sixth hole, the seventh hole, the eighth hole and the ninth hole are arranged on the plane with the center axis of the spherical probe head at 90°, and the tenth hole, the eleventh hole, the twelfth hole and the thirteenth hole are arranged on the plane with the center axis of the spherical probe head at 120°. In addition, since a plurality of pressure guide channels are arranged in the shell support rod, in order to realize the minimum disturbance to the incoming flow, the bending angle of the internal pipeline is calculated by numerical simulation, and the curvature is represented as K is 0.29R, wherein R is the radius of the spherical head.

[0028] Embodiment two As Figure 5 shown is a flow chart of a fluid velocity vector measurement method in the embodiment two of the present application, the fluid velocity vector measurement method includes the following steps S01-S04, wherein: S01, according to the physical position of each pressure hole in the surface of the spherical probe, the thirteen-hole probe is partitioned, and the surface of the spherical probe is divided into nine partition regions, wherein any one of the partition regions is composed of at least three pressure holes.

[0029] According to the physical position of the thirteen pressure holes in embodiment one, the surface of the spherical probe is divided into nine partition regions: the first region is composed of the first hole, the second hole, the third hole, the fourth hole and the fifth hole; the second region is composed of the second hole, the third hole and the sixth hole; the third region is composed of the second hole, the fourth hole and the seventh hole; the fourth region is composed of the third hole, the fifth hole and the eighth hole; the fifth region is composed of the fourth hole, the fifth hole and the ninth hole; the sixth region is composed of the second hole, the sixth hole, the seventh hole and the tenth hole; the seventh region is composed of the third hole, the sixth hole, the eighth hole and the eleventh hole; the eighth region is composed of the fourth hole, the seventh hole, the ninth hole and the twelfth hole; and the ninth region is composed of the fifth hole, the eighth hole, the ninth hole and the thirteenth hole.

[0030] S02, a calibration experiment is performed on the thirteen-hole probe, the angle posture of the thirteen-hole probe is adjusted, a plurality of sets of calibration data under a known flow field are obtained, and an angle coefficient characteristic curve of each partition region is calculated using the calibration data.

[0031] When performing the thirteen-hole probe calibration test, in order to measure the full plane flow field (i.e. a known flow field) required by the design, a double-turntable is used to control the posture of the thirteen-hole probe at different pitch angles and different azimuth angles, and calibration data of the thirteen-hole probe at different postures are obtained, wherein the calibration data includes hole pressure data of the thirteen pressure holes, and the specific process is as follows: The pitch angle range is set to 90°, and the pitch angle interval is 5°; the azimuth angle is rotated by 360°, and the interval is 5°. The data of the first calibration point is collected, the azimuth angle of the first calibration point is 0°, and the pitch angle is 0°; the data of the next calibration point is collected, the azimuth angle turntable is rotated by 5°, the azimuth angle turntable remains unchanged, and the pitch angle turntable is rotated from 0° to 90° at an interval of 5°; the azimuth angle turntable is continuously rotated by 5°, and the data of the next calibration point is collected, the azimuth angle turntable remains unchanged, and the pitch angle turntable is rotated from 90° to 0° at an interval of 5°; the process of collecting the calibration point data is repeated until the calibration points in the range of 360° of the azimuth angle turntable are all collected; when the azimuth angle is rotated for one revolution, the last pitch angle test should be stopped at 90° position, and the zero return is performed to return to the original position, and the calibration test is completed.

[0032] Calibration data is obtained According to the subsequent, the calibration data needs to be partitioned, first, the thirteen pressure hole hole pressure data is arranged according to the pitch angle from small to large, and the inside and outside partition is judged according to the thirteen pressure hole hole pressure data, the first two pressure value size hole pressure data is found, which is the data of the first five holes, and the largest pitch angle in these data is found , all the data located at the pitch angle 0 Belong to the first area.

[0033] Then judge other test area, in the pressure value size of the hole pressure data, the hole pressure data of the first five holes is eliminated, and the smallest pitch angle In the remaining data is found. According to the practical experience, it is found that in most cases ≥ , that is, there is a part of overlap between the first area and the second area to the fifth area, that is, the position of the test area boundary, according to the concept of overlapping domain of thirteen hole probe, which does not affect the subsequent interpolation operation, which is equivalent to partition expansion, which ensures that the test point can be normally interpolated in the use process, even if < It is also true, the key is that the airflow receiving point on the spherical probe is set in the connecting line area of the four pressure holes.

[0034] Then the division of the second area to the fifth area is carried out, the selected data is sorted according to the azimuth angle from small to large, the first hole and the second hole vertical line position is 0° azimuth angle starting point, so the second area belongs to the range of azimuth angle 0° ~ 90°, after selecting this part of the data, the maximum pressure hole ladder judgment is carried out, the pressure data is sorted from large to small, the first four pressure value of the pressure hole number is recorded, and the sum is added, the data with the sum of 12 is extracted, which is the data of the second area.

[0035] Similarly, the third area is positioned according to the azimuth angle 90°~180°, and the data of the third area is screened out according to the sum of the first four hole numbers being 14; the fourth area is positioned according to the azimuth angle 180°~270°, and the data of the fourth area is screened out according to the sum of the first four hole numbers being 17; the fifth area is positioned according to the azimuth angle 270°~360°, and the data of the fifth area is screened out according to the sum of the first four hole numbers being 19. It is worth noting that the azimuth angle 0° and 360° are measured and recorded twice under each pitch angle. Theoretically, the two data should be the same, but due to environmental factors, the two data actually have some differences, so the average of the two data under each pitch angle is the data to be used. Then the data of the sixth area to the ninth area is screened out. Since the sixth area to the ninth area is staggered with the second area to the fifth area, that is, the sixth area spans the azimuth angle range of 315°~45°; the seventh area spans the azimuth angle range of 45°~135°; the eighth area spans the azimuth angle range of 225°~315°; and the ninth area spans the azimuth angle range of 135°~225°. Similar to the division method of the second area to the fifth area, first, the azimuth angle is divided, and then the data belonging to each area is determined according to the sum of the first four hole numbers: the data with the sum of the first four hole numbers being 25 belongs to the sixth area; the data with the sum of the first four hole numbers being 28 belongs to the seventh area; the data with the sum of the first four hole numbers being 32 belongs to the eighth area; and the data with the sum of the first four hole numbers being 35 belongs to the ninth area.

[0036] After the hole pressure data is divided into zones, different calibration coefficient calculation formulas are used for different zones to obtain the angle coefficient characteristic curve diagram of the corresponding zone, which is used for subsequent calculation of the flow field parameters of the measured flow field.

[0037] S03, the hole pressure data of each pressure hole in the known flow field is dimensionless processed, and the dimensionless pressure coefficient is obtained, the dimensionless pressure coefficient is fitted, and a general calculation model is established.

[0038] The hole pressure data of each pressure hole is screened out from the calibration data, the least square method is used to fit the flow characteristics, the hole pressure data of each pressure hole is dimensionless processed to convert into dimensionless pressure coefficient, and a double variable third order polynomial of pressure coefficient is obtained. When the data points match the measurement area, the corresponding pressure coefficient is calculated, and the desired flow field flow characteristics can be determined by expanding the fourth order polynomial, wherein the dimensionless pressure coefficient includes pitch angle, azimuth angle, total pressure coefficient and static pressure coefficient:

[0039] wherein, represents the pitch angle θ, the azimuth angle , the total pressure coefficient K O and the static pressure coefficient K QAny one of the four dimensionless pressure coefficients, subscript i represents the i-th measuring point; K represents the calibration coefficient, the superscript A represents the calibration coefficient corresponding to all of a certain sequence of K, and the subscript represents the corresponding term in the expansion; wherein, any one of the dimensionless pressure coefficients in any partition region contains 15 calibration coefficients, so that each measuring point, each measuring area of the thirteen-hole probe needs 60 calibration coefficients, and the thirteen-hole probe has nine partition regions, a total of 540 calibration coefficients, which can be calculated by the calibration data.

[0040] After establishing the above general calculation model, for any one partition region, input the dimensionless pressure coefficients of the calibration points of the thirteen-hole probe at different angles of the partition region into the general calculation model, and obtain the calibration matrix:

[0041] The simplified form is:

[0042]

[0043] Wherein, the A matrix contains the dimensionless pressure coefficients of n calibration points, the K matrix contains the calibration coefficients of n calibration points, K A T is the transpose matrix of K A , K A is a matrix composed of calibration coefficients.

[0044] After the calibration coefficient matrix is calculated, the calibration constant K is determined by fitting the least squares curve of the probe calibration data obtained by the calibration experiment, and the calibration process is completed.

[0045] S04, collecting the measured orifice pressure data of the thirteen-hole probe in the to-be-measured flow field, positioning the partition to which the measured orifice pressure data belongs, and calling the angle characteristic curve of the corresponding partition in the general calculation model to calculate the flow field parameters of the to-be-measured flow field.

[0046] In the flow field to be measured, the received measurement pressure values are calculated as dimensionless pressure coefficients according to the pressures of the 13 pressure holes measured by the thirteen-hole probe at a series of different incident angles. Angle coefficient characteristic curves are plotted using the angle calibration coefficients unique to each zone. After the determination of the zones is completed, the corresponding angle coefficient characteristic curves are linearly interpolated according to the calculated angle calibration coefficients. When the interpolation method is used to process the angle calibration coefficients, the following two parts of work need to be completed in sequence: first, find the smallest nearest block area that can contain the target point in the angle coefficient characteristic curve, and then perform interpolation calculation on the two angle values according to the position of the target point relative to the block area. After finding the four closest calibration points (i.e., locking the target point in the smallest block area), interpolation calculation is performed according to the distance relationship of the four edges of the quadrilateral formed by the target point and the four calibration points.

[0047] Specifically, please refer to Figure 6 , represent the proximity of the test point to this side, if =0.4, then the test point , similarly, the value of the azimuth angle of the test point is calculated, if , . .

[0048] In summary, according to the fluid velocity vector measurement method provided by the present application, the thirteen-hole probe is placed in a known flow field, the calibration data is obtained by adjusting the attitude of the thirteen-hole probe at different angles, and the angle coefficient characteristic curve of each partition region is calculated using the calibration data, and the hole pressure data of each pressure hole is dimensionless processed, focusing on the correlation between pressure difference and flow field parameters to eliminate the influence of absolute pressure size, and using the least square fitting method, a general calculation model is constructed for subsequent testing of the flow field parameters of the test flow field. The present application overcomes the shortcomings of the prior art that can only receive airflow along the axis direction of the probe by providing multiple pressure holes on the spherical probe for receiving airflow at multiple orientations, enabling large-angle deflection measurement. At the same time, the present application can directly calculate the flow field parameters of the test flow field by measuring the hole pressure data in the pressure holes and using the general calculation model, which is less affected by external environmental factors, and directly focuses on the pressure data of the pressure holes, making the calculation results more accurate and reliable.

[0049] Example Three The present application also provides a fluid velocity vector measurement system, please refer to Figure 7 , which is a structural schematic diagram of the fluid velocity vector measurement system in the second embodiment of the present application, the fluid velocity vector measurement system comprises: The regional partition module 11 is used for partitioning the thirteen-hole probe according to the physical position of each pressure hole in the spherical probe surface, and dividing the spherical probe surface into nine partition regions, wherein each partition region is composed of at least three pressure holes; The curve calculation module 12 is used for performing a calibration experiment on the thirteen-hole probe, adjusting the angle posture of the thirteen-hole probe, obtaining a plurality of sets of calibration data under a known flow field, and calculating the angle coefficient characteristic curve of each partition region by using the calibration data. The model establishment module 13 is used for performing dimensionless processing on the hole pressure data of each pressure hole under the known flow field, obtaining a dimensionless pressure coefficient, fitting the dimensionless pressure coefficient, and establishing a general calculation model. The parameter calculation module 14 is used for collecting the measured hole pressure data of the thirteen-hole probe under the to-be-measured flow field, positioning the partition to which the measured hole pressure data belongs, and calculating the flow field parameter of the to-be-measured flow field by using the angle characteristic curve of the corresponding partition in the general calculation model.

[0050] Embodiment four In another aspect, the application further provides a computer readable storage medium having one or more computer programs stored thereon, which, when executed by a processor, implement the above-mentioned fluid velocity vector measurement method.

[0051] Those skilled in the art can understand that the logic or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequence list of executable instructions for implementing the logic function, which can be specifically implemented in any computer readable storage medium for use by or in conjunction with an instruction execution system, device or equipment (such as a computer-based system, a system including a processor or other system that can fetch and execute instructions from the instruction execution system, device or equipment). For the purpose of this specification, the "computer readable storage medium" can be any device that can contain, store, communicate, propagate or transport programs for use by or in conjunction with an instruction execution system, device or equipment, or in conjunction with these instruction execution systems, devices or equipment.

[0052] More specific examples (a non-exhaustive list) of the computer-readable storage medium include the following: an electrical connection having one or more wires (electrical device), a portable computer diskette (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer readable storage medium can also be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, for example via an optical scanner, then compiled, interpreted, or otherwise processed, and stored in a computer memory in order to be executed.

[0053] Embodiment Five Figure 8 A structural block diagram of an electronic device is provided for Embodiment Four. The electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the fluid velocity vector measurement method in the above embodiments when executing the program. Figure 8 The electronic device 30 shown is merely an example and should not limit the function and scope of use of the embodiments of the present application.

[0054] As shown in Figure 8 The electronic device 30 can be in the form of a general computing device, for example, it can be a server device. The components of the electronic device 30 can include, but are not limited to, the above-mentioned at least one processor 31, the above-mentioned at least one memory 32, and a bus 33 connecting different system components, including the memory 32 and the processor 31.

[0055] The bus 33 includes a data bus, an address bus, and a control bus.

[0056] The memory 32 can include volatile memory, such as a RAM 321 (Random Access Memory), and / or cache memory 322, and can further include a ROM 323 (Read-Only Memory).

[0057] The memory 32 can further include a program tool 325 having a set of (at least one) program modules 324, such as an operating system, one or more application programs, other program modules, and program data, and each of these examples or some combination thereof can include implementation of a network environment.

[0058] The processor 31 performs various function applications and data processing by running the computer program stored in the memory 32, such as the fluid velocity vector measurement method of the present application as described above.

[0059] The electronic device 30 can also communicate with one or more external devices 34 such as a keyboard or a pointing device, by way of I / O interface 35. Furthermore, the electronic device 30 can communicate to one or more networks, such as a local area network (LAN), a wide area network (WAN), and / or the public network, such as the Internet, by way of the network adapter 36. As Figure 8 illustrated, the network adapter 36 communicates to the other components of the model-generated electronic device 30 via bus 33. It should be appreciated that although not shown, other hardware and / or software components could be used in conjunction with the model-generated electronic device 30. These include, but are not limited to, microcode, device drivers, redundant processing units, disk drive arrays, RAID systems, tape drives, and data archival storage systems, etc.

[0060] It should be noted that although several units / modules or sub-units / modules of the electronic device are mentioned in the foregoing detailed description, such a division is merely exemplary and not mandatory. Indeed, according to embodiments of the application, the features and functions of two or more units / modules described above can be embodied in one unit / module. Conversely, the features and functions of one unit / module described above can be further divided into units / modules embodied by several units / modules.

[0061] In the description of the present specification, the description referring to the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" etc. means that the particular feature, structure, material or characteristic being described in connection with the embodiment or example is included in at least one embodiment or example of the present application. Illustrative descriptions of the above terms do not necessarily mean the same embodiment or example in the present specification. Moreover, the particular feature, structure, material or characteristic being described can be combined in any one or more embodiments or examples in a suitable manner.

[0062] The above-described embodiments are merely some embodiments of the present application, which are described in a relatively specific and detailed manner. However, it should not be understood that the scope of the patent right of the present application is limited to the above-described embodiments. It should be noted that, for those skilled in the art, several modifications and improvements can be made without departing from the concept of the present application, and these all belong to the scope of protection of the present application. Therefore, the scope of protection of the patent of the present application should be subject to the appended claims.

Claims

1. A thirteen-hole probe, characterized in that, The thirteen-hole probe comprises a spherical probe head, a shell support rod fixed to the spherical probe head, and a plurality of pressure guide channels arranged in the shell support rod; The spherical probe head comprises thirteen pressure holes symmetrically arranged on the spherical probe head body, a first hole arranged at the center of the spherical probe head and arranged opposite to the shell support rod, and the remaining pressure holes evenly arranged on planes at 30°, 90° and 120° with respect to the central axis of the spherical probe head, each layer having four pressure holes, and the pressure holes of each layer being staggered. The second hole, the third hole, the fourth hole and the fifth hole are arranged on a plane at 30° with respect to the central axis of the spherical probe head, the sixth hole, the seventh hole, the eighth hole and the ninth hole are arranged on a plane at 90° with respect to the central axis of the spherical probe head, and the tenth hole, the eleventh hole, the twelfth hole and the thirteenth hole are arranged on a plane at 120° with respect to the central axis of the spherical probe head.

2. A method of fluid velocity vector measurement, characterized by, The fluid velocity vector measurement method is applied to the thirteen-hole probe of claim 1, and the method comprises: According to the physical positions of the pressure holes on the surface of the spherical probe head, the thirteen-hole probe is divided into nine partition regions, wherein each partition region is enclosed by at least three pressure holes. A calibration experiment is performed on the thirteen-hole probe, the angle posture of the thirteen-hole probe is adjusted, a plurality of sets of calibration data under a known flow field are obtained, and the angle coefficient characteristic curve of each partition region is calculated using the calibration data. The hole pressure data of each pressure hole under the known flow field is processed by dimensionless, and the dimensionless pressure coefficient is obtained. The dimensionless pressure coefficient is fitted to establish a general calculation model. The measured hole pressure data of the thirteen-hole probe under the to-be-measured flow field is collected, the partition to which the measured hole pressure data belongs is located, and the flow field parameters of the to-be-measured flow field are calculated by calling the angle characteristic curve of the corresponding partition in the general calculation model.

3. The fluid velocity vector measurement method of claim 2, wherein, The step of dividing the thirteen-hole probe comprises: The first partition region is composed of the first hole, the second hole, the third hole, the fourth hole and the fifth hole; the second partition region is composed of the second hole, the third hole and the sixth hole; the third partition region is composed of the second hole, the fourth hole and the seventh hole; the fourth partition region is composed of the third hole, the fifth hole and the eighth hole; the fifth partition region is composed of the fourth hole, the fifth hole and the ninth hole; the sixth partition region is composed of the second hole, the sixth hole, the seventh hole and the tenth hole; the seventh partition region is composed of the third hole, the sixth hole, the eighth hole and the eleventh hole; the eighth partition region is composed of the fourth hole, the seventh hole, the ninth hole and the twelfth hole; and the ninth partition region is composed of the fifth hole, the eighth hole, the ninth hole and the thirteenth hole.

4. The fluid velocity vector measurement method of claim 2, wherein, The step of performing a calibration experiment on the thirteen-hole probe and adjusting the angle posture of the thirteen-hole probe to obtain a plurality of sets of calibration data under a known flow field comprises: A calibration experiment is performed on the thirteen-hole probe under a known flow field, the posture of the thirteen-hole probe under different pitch angles and different azimuth angles is adjusted using a double-turntable, and calibration data of the thirteen-hole probe under each different posture is obtained; After obtaining the calibration data, the hole pressure data is partitioned according to a preset rule, the calibration data in the nine partition regions is matched, and the angle coefficient characteristic curve of each partition is calculated using the matched calibration data.

5. The method of fluid velocity vector measurement according to claim 2, wherein, The step of non-dimensionalizing the hole pressure data of each pressure hole under the known flow field, obtaining the non-dimensional pressure coefficient, fitting the non-dimensional pressure coefficient, and establishing a general calculation model, The hole pressure data of each pressure hole is screened from the calibration data, and after obtaining the hole pressure data of each pressure hole, the hole pressure data is non-dimensionalized and converted into a non-dimensional pressure coefficient, wherein the non-dimensional pressure coefficient includes an angle coefficient, a total pressure coefficient, and a static pressure coefficient. The non-dimensional pressure coefficient is fitted by using a least square fitting method to establish a general calculation model: wherein, denotes the pitch angle θ, the azimuth angle , the total pressure coefficient K O and the static pressure coefficient K Q any one of the four dimensionless pressure coefficients, the subscript i denotes the i-th measuring point; K represents the calibration coefficient, the superscript A indicates the calibration coefficient corresponding to all of a certain sequence of K, and the subscript represents the corresponding term in the expansion; wherein, any one of the dimensionless pressure coefficients in any sub-region comprises 15 calibration coefficients.

6. The method of fluid velocity vector measurement according to claim 5, wherein, After the general calculation model is established, for any one partition region, the non-dimensional pressure coefficients of the calibration points of the thirteen-hole probe in different angle postures in the partition region are input into the general calculation model, and a calibration matrix is obtained: The simplified form is: where A matrix contains dimensionless pressure coefficients of n calibration points, K matrix contains calibration coefficients of n calibration points, K A T is the transpose matrix of K A A is the matrix composed of calibration coefficients.​ 7. A fluid velocity vector measurement system characterized by, The fluid velocity vector measurement system is used to implement the fluid velocity vector measurement method according to any one of claims 2-6, and the system comprises: A region partition module is configured to partition the thirteen-hole probe according to the physical positions of the pressure holes on the surface of the spherical probe, and divide the surface of the spherical probe into nine partition regions, wherein any one of the partition regions is enclosed by at least three pressure holes; A curve calculation module is configured to perform a calibration experiment on the thirteen-hole probe, adjust the angle posture of the thirteen-hole probe, obtain a plurality of groups of calibration data under a known flow field, and calculate an angle coefficient characteristic curve of each partition region using the calibration data; A model establishment module is configured to non-dimensionalize the hole pressure data of each pressure hole under the known flow field, obtain a non-dimensional pressure coefficient, fit the non-dimensional pressure coefficient, and establish a general calculation model; A parameter calculation module is configured to collect measured hole pressure data of the thirteen-hole probe under a to-be-measured flow field, locate the to-be-measured flow field, and call an angle characteristic curve of a corresponding partition region in the general calculation model to calculate a flow field parameter of the to-be-measured flow field.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The program is executed by the processor to implement the fluid velocity vector measurement method according to any one of claims 2-6.

9. An electronic device comprising a memory, a processor, and a computer program stored on the memory and running on the processor, characterized in that, The processor executes the computer program to implement the fluid velocity vector measurement method according to any one of claims 2-6.

Images