A thirteen-hole probe and fluid velocity vector measurement method and system

Abstract

The application discloses a thirteen-hole probe and a fluid velocity vector measurement method and system, and the method comprises the following steps: according to the physical positions of each pressure hole on the surface of a spherical probe, the thirteen-hole probe is partitioned; 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 by using the calibration data; dimensionless processing is performed on the hole pressure data of each pressure hole under the known flow field, fitting is performed on the dimensionless pressure coefficient obtained in the dimensionless processing, and a general calculation model is established; measured hole pressure data of the thirteen-hole probe under a 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.

Description

Technical Field

[0001] This invention relates to the field of probe measurement technology, specifically to a thirteen-well probe and a method and system for measuring fluid velocity vectors. Background Technology

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

[0003] Currently, airflow measurement technology is becoming increasingly mature. Pitot tubes can be used to capture the total pressure of the airflow through the front total pressure orifice and the static pressure of the environment through the side wall static pressure orifice. Then, using Bernoulli's equation, the airflow velocity can be calculated by the difference between the total pressure and the static pressure. Alternatively, a hot-wire anemometer can be used to take advantage of the characteristic that the airflow carries away heat, and the airflow velocity can be inferred by using the change in the resistance of the hot wire.

[0004] However, Pitot tubes can only measure one-dimensional velocity along the probe axis, and hot-wire anemometers are easily contaminated by dust, oil, and water droplets in the airflow, leading to decreased measurement accuracy or even failure. Therefore, both Pitot tubes and hot-wire anemometers have defects in the process of measuring airflow velocity vectors. Summary of the Invention

[0005] Based on this, the purpose of the present invention is to provide a thirteen-well probe and a fluid velocity vector measurement method and system, which aims to solve the problems of limitations and susceptibility to environmental factors in the current use of Pitot tubes or hot-wire anemometers to measure airflow velocity.

[0006] To achieve the above objectives, this invention proposes a fluid velocity vector measurement method, which is applied to a thirteen-well probe. The method includes:

[0007] Based on the physical location of each pressure hole on the surface of the spherical probe, the thirteen-hole probe is divided into nine partitioned areas on the surface of the spherical probe. Each partitioned area is surrounded by at least three pressure holes.

[0008] A calibration experiment was conducted on the thirteen-hole probe. The angle and attitude of the thirteen-hole probe were adjusted, and multiple sets of calibration data were obtained under a known flow field. The calibration data were then used to calculate the angle coefficient characteristic curve of each partition region.

[0009] The pore pressure data of each pressure orifice under a known flow field are processed to be dimensionless, and the dimensionless pressure coefficient is obtained. The dimensionless pressure coefficient is then fitted to establish a general calculation model.

[0010] Collecting the measured orifice pressure data of the thirteen-hole probe under 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.

[0011] According to an aspect of the above technical solution, the step of partitioning the thirteen-hole probe is:

[0012] The first zone is composed of the first hole, the second hole, the third hole, the fourth hole and the fifth hole; the second zone is composed of the second hole, the third hole and the sixth hole; the third zone is composed of the second hole, the fourth hole and the seventh hole; the fourth zone is composed of the third hole, the fifth hole and the eighth hole; the fifth zone is composed of the fourth hole, the fifth hole and the ninth hole; the sixth zone is composed of the second hole, the sixth hole, the seventh hole and the tenth hole; the seventh zone is composed of the third hole, the sixth hole, the eighth hole and the eleventh hole; the eighth zone is composed of the fourth hole, the seventh hole, the ninth hole and the twelfth hole; and the ninth zone is composed of the fifth hole, the eighth hole, the ninth hole and the thirteenth hole.

[0013] 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:

[0014] The thirteen-hole probe is calibrated under the 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;

[0015] After obtaining the calibration data, the orifice 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.

[0016] According to an aspect of the above technical solution, in the step of performing non-dimensional processing on the orifice pressure data of each pressure hole under the known flow field and obtaining a non-dimensional pressure coefficient, fitting the non-dimensional pressure coefficient, and establishing a general calculation model,

[0017] The orifice pressure data of each pressure hole is selected from the calibration data, after obtaining the orifice pressure data of each pressure hole, the orifice pressure data is processed in a non-dimensional manner and is 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 in a least square method fitting manner, and a general calculation model is established:

[0018]

[0019] wherein, represents the pitch angle θ and the azimuth angle , the total pressure coefficient K Oand static pressure coefficient K Q Any 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 sub-region comprises 15 calibration coefficients.

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

[0021]

[0022] The simplified form is:

[0023]

[0024]

[0025] 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.

[0026] The application also provides a fluid velocity vector measurement system for implementing the fluid velocity vector measurement method, and the system comprises:

[0027] 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 sub-regions, wherein any one of the sub-regions is composed of at least three pressure holes.

[0028] 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 an angle coefficient characteristic curve of each sub-region using the calibration data.

[0029] 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.

[0030] A parameter calculation module is configured to collect measured orifice pressure data of the thirteen-hole probe under a to-be-tested flow field, locate a corresponding partition of the measured orifice pressure data, and call an angle characteristic curve of the corresponding partition in a general calculation model to calculate a flow field parameter of the to-be-tested flow field.

[0031] 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.

[0032] The spherical probe head comprises thirteen pressure holes symmetrically arranged on a body of the spherical probe head, wherein a first hole is arranged at a 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 a central axis of the spherical probe head, wherein each of the planes has four pressure holes, and the pressure holes of each of the planes are arranged alternately.

[0033] The first hole, 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.

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

[0035] 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 implements the fluid velocity vector measurement method when executing the computer program.

[0036] In summary, according to the fluid velocity vector measurement method, 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, the angle coefficient characteristic curve of each partition area is calculated by using the calibration data, the orifice pressure data of each pressure hole is processed by using the dimensionless method, the correlation between the pressure difference and the flow field parameter is focused on, the influence of the absolute pressure is eliminated, and the general calculation model is constructed by using the least square fitting method, so as to be used for testing 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, the shortcomings of the prior art that the airflow is only received along the axis direction of the probe are overcome, the large-angle deflection measurement is realized, the orifice pressure data in the pressure holes are measured, the flow field parameter of the to-be-tested flow field is directly calculated by using the general calculation model, the influence of external environmental factors is small, the pressure data of the pressure holes are directly focused on, and the calculation result is more accurate and reliable.

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

[0038] Figure 1 Structure diagram of the spherical probe in the embodiment one of the present application;

[0039] Figure 2 Structure diagram of the spherical probe in the embodiment one of the present application from another perspective;

[0040] Figure 3 Radial sectional view of the thirteen-hole probe in the embodiment one of the present application;

[0041] Figure 4 Axial sectional view of the thirteen-hole probe in the embodiment one of the present application;

[0042] Figure 5 Flow chart of the fluid velocity vector measurement method in the embodiment two of the present application;

[0043] Figure 6 Interpolation calculation example diagram in the embodiment two of the present application;

[0044] Figure 7 Structure diagram of the fluid velocity vector measurement system in the embodiment three of the present application;

[0045] Figure 8 Structure block diagram of the electronic device in the embodiment five of the present application.

[0046] Explanation of the drawing component symbols:

[0047] Spherical probe D, shell support B, pressure guide channel C DETAILED DESCRIPTION

[0048] In order to make the objects, features and advantages of the present application more apparent, the specific embodiments of the present application are 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, the purpose of providing these embodiments is to make the disclosure of the present application more thorough and comprehensive.

[0049] It is to be understood that where an element such as a layer, region or substrate is described as being "on" another element, it can be directly on the other element or intervening elements can also be present. Where an element such as a layer, region or substrate is described as being "connected" or "coupled" to another element, it can be directly connected to the other element or intervening elements can also be present. As used herein, the term "and / or" includes any and all combinations of one or more of the associated listed items. Relative terms such as "lower," "upper," "horizontal," "vertical," "above," "below," "up," "down," "top," "bottom," and the like as can be used herein for the purposes of description only and are not intended to be limiting of an orientation or an actual position of an described apparatus or element unless specific installation is made herein. Terms such as "first," "second," "third," etc. as can be used herein are also not intended to be limiting, unless specifically stated otherwise.

[0050] In this application, terms "mounting", "connecting", "connecting", "fixing" and the like should be understood broadly, for example, can be fixed connection, can also be detachable connection, or integral connection; can be mechanical connection, can also be electrical connection; can be directly connected, can also be indirectly connected through intermediate medium, can be internal communication of two elements. For those skilled in the art, the specific meaning of the above 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 of the associated listed items.

[0051] Embodiment one

[0052] Please refer to Figures 1-4 , it is a structure schematic 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:

[0053] 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 embodiment takes thirteen holes as an example for illustration.

[0054] 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 arranged at the center of the spherical probe head and arranged 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°.

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

[0056] Embodiment two

[0057] As Figure 5 shown is a flow chart of a fluid velocity vector measurement method in embodiment two of the application, the fluid velocity vector measurement method comprises the following steps S01-S04, wherein:

[0058] S01, according to the physical position of each pressure hole in the surface of the spherical probe, the thirteen-hole probe is divided into zones, and nine zone areas are divided on the surface of the spherical probe, wherein any one of the zone areas is composed of at least three pressure holes.

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

[0060] 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 the angle coefficient characteristic curve of each zone area is calculated using the calibration data.

[0061] In the calibration test of the thirteen-hole probe, in order to measure the full plane flow field (i.e. the 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:

[0062] Set the pitch angle range to 90° with a 5° interval; rotate the azimuth angle 360° with a 5° interval. Collect data for the first calibration point, with both azimuth and pitch angles at 0°. Collect data for the next calibration point, rotating the azimuth turntable 5° while keeping it stationary, and rotating the pitch turntable from 0° to 90° in 5° increments. Continue rotating the azimuth turntable 5° and collect data for the next calibration point, again keeping it stationary, and rotating the pitch turntable back to 0° in 5° increments. Repeat this process until all calibration points within the 360° range of the azimuth turntable have been collected. After the azimuth angle has rotated one full circle, the final pitch angle test should be completed, stopping at 90° to return to the original position and the calibration test is finished.

[0063] Get calibration numbers Following this, the calibration data needs to be partitioned. First, the orifice pressure data of the thirteen pressure holes are arranged in ascending order of pitch angle. Then, based on the orifice pressure data of the thirteen pressure holes, the inner and outer partitions are determined, and the orifice pressure data with the two largest pressure values ​​are found. These data are the data for holes one through five. Finally, the largest pitch angle among these data is found. Then all those located at pitch angle 0~ All the data belongs to Zone 1.

[0064] Then, for other measurement areas, after removing the pore pressure data from holes 1 to 5 from the two pore pressure data with the highest pressure values, find the smallest pitch angle among the remaining data. Based on practical experience, it has been found that in most cases... ≥ In other words, there is some overlap between Zone 1 and Zones 2 through 5, specifically at the boundary of the measurement areas. According to the concept of the overlapping region for the thirteen-well probe, this does not affect subsequent interpolation calculations. It's equivalent to a zone expansion, ensuring that all measurement points can be interpolated normally during use, even if... < The same applies to the case where the airflow receiving point on the spherical probe is located within the area connected by the four pressure holes.

[0065] Next, the data was divided into zones two through five. The selected data was sorted from smallest to largest according to the azimuth angle. The positions of the vertical lines of holes one and two were taken as the starting points of the 0° azimuth angle. Therefore, the range of zone two was azimuth angle 0°~90°. After selecting this part of the data, the maximum pressure hole step judgment was performed. The pressure data was sorted from largest to smallest. The four pressure hole numbers with the largest pressure values ​​were recorded and summed. The data with a sum of 12 were extracted. These are the data of zone two.

[0066] Similarly, No. 3 area is located according to the azimuth angle 90°~180°, and the data of No. 3 area is screened out according to the sum of the first four hole numbers being 14; No. 4 area is located according to the azimuth angle 180°~270°, and the data of No. 4 area is screened out according to the sum of the first four hole numbers being 17; No. 5 area is located according to the azimuth angle 270°~360°, and the data of No. 5 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 actual data will have some differences, so the average of the two data under each pitch angle is the data to be used. Then the data of No. 6 area-No. 9 area is screened out. Since No. 6 area-No. 9 area is staggered with No. 2 area-No. 5 area, that is, No. 6 area spans the azimuth angle range of 315°~45°; No. 7 area spans the azimuth angle range of 45°~135°; No. 8 area spans the azimuth angle range of 225°~315°; No. 9 area spans the azimuth angle range of 135°~225°. Similar to the division method of No. 2 area-No. 5 area, first divide according to the azimuth angle, and then determine the data belonging to each area 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 No. 6 area; the data with the sum of the first four hole numbers being 28 belongs to No. 7 area; the data with the sum of the first four hole numbers being 32 belongs to No. 8 area; the data with the sum of the first four hole numbers being 35 belongs to No. 9 area.

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

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

[0069] 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 four-order polynomial can be expanded to determine the desired flow characteristics of the flow field, wherein the dimensionless pressure coefficient includes pitch angle, azimuth angle, total pressure coefficient and static pressure coefficient:

[0070]

[0071] wherein, represents the pitch angle θ, the azimuth angle , the total pressure coefficient K O and the static pressure coefficient KQ Any 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.

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

[0073]

[0074] The simplified form is:

[0075]

[0076]

[0077] 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.

[0078] 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.

[0079] 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.

[0080] 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 drawn 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.

[0081] 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 , , .

[0082] 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 area 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.

[0083] Example Three

[0084] 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:

[0085] The region partition module 11 is configured to partition the thirteen-hole probe according to the physical positions of the pressure holes in the spherical probe surface, and divide the spherical probe surface into nine partition regions, wherein each partition region is enclosed by at least three pressure holes.

[0086] The curve calculation module 12 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 an angle coefficient characteristic curve of each partition region by using the calibration data.

[0087] The model establishment module 13 is configured to perform dimensionless processing on the hole pressure data of each pressure hole under the known flow field, obtain a dimensionless pressure coefficient, and perform fitting on the dimensionless pressure coefficient to establish a general calculation model.

[0088] The parameter calculation module 14 is configured to collect measured hole pressure data of the thirteen-hole probe under a to-be-measured flow field, locate the partition to which the measured hole pressure data belongs, and calculate 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.

[0089] Embodiment four

[0090] 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.

[0091] 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.

[0092] 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.

[0093] Embodiment Five

[0094] Figure 8 A structural block diagram of an electronic device provided for Embodiment Four is shown. 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 bring any limitation to the function and use range of the embodiments of the present application.

[0095] 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.

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

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

[0098] 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.

[0099] 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.

[0100] 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 by way of the bus 33. It should be appreciated that although not shown, other hardware and / or software modules could be used in conjunction with the model-generated electronic device 30. Such as, but not limited to, microcode, device drivers, redundant processing units, disk drive arrays, RAID systems, tape drives, and data archival storage systems, etc.

[0101] 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.

[0102] 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 expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the particular features, structures, materials, or characteristics can be combined in any suitable manner in one or more embodiments or examples.

[0103] The above-described embodiments are merely illustrative of several embodiments of the present application and are not intended to limit the scope of the present application. 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 should all fall within the scope of the present application. Therefore, the scope of the present application should be determined by 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, any one layer being provided with four pressure holes, and the thirteen pressure holes being arranged in a 1-4-4-4 structure along the length direction of the shell support rod, 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, and the curvature of the bending angle of the plurality of pressure guide channels is 0.29R, R being the radius 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 zones, and nine zone areas are divided on the surface of the spherical probe head, wherein any one of the zone areas 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 an angle coefficient characteristic curve of each zone area is calculated using the calibration data; The hole pressure data of each pressure hole under the known flow field are processed by dimensionless, and a dimensionless pressure coefficient is obtained, the dimensionless pressure coefficient is fitted, and a general calculation model is established; The measured hole pressure data of the thirteen-hole probe under a to-be-measured flow field are collected, the zone area to which the measured hole pressure data belong 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 zone area in the general calculation model.

3. The fluid velocity vector measurement method of claim 2, wherein, The step of dividing the thirteen-hole probe into zones comprises: The first zone is composed of the first hole, the second hole, the third hole, the fourth hole and the fifth hole; the second zone is composed of the second hole, the third hole and the sixth hole; the third zone is composed of the second hole, the fourth hole and the seventh hole; the fourth zone is composed of the third hole, the fifth hole and the eighth hole; the fifth zone is composed of the fourth hole, the fifth hole and the ninth hole; the sixth zone is composed of the second hole, the sixth hole, the seventh hole and the tenth hole; the seventh zone is composed of the third hole, the sixth hole, the eighth hole and the eleventh hole; the eighth zone is composed of the fourth hole, the seventh hole, the ninth hole and the twelfth hole; and the ninth zone 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 are obtained. After the calibration data is obtained, 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, In the step of performing non-dimensionalization processing on the hole pressure data of each pressure hole under the known flow field and obtaining the non-dimensional pressure coefficient, the non-dimensional pressure coefficient is fitted to establish a general calculation model, The hole pressure data of each pressure hole is screened from the calibration data, the hole pressure data of each pressure hole is obtained, the hole pressure data is processed by non-dimensionalization and converted into a non-dimensional pressure coefficient, 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 method, and a general calculation model is established: 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 partition region, the non-dimensional pressure coefficients of the calibration points of the thirteen-hole probe at 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: The region partition module is used 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 composed of at least three pressure holes; The curve calculation module is used 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 the angle coefficient characteristic curve of each partition region using the calibration data; The model establishment module is used to perform non-dimensionalization processing on the hole pressure data of each pressure hole under the known flow field, obtain the non-dimensional pressure coefficient, fit the non-dimensional pressure coefficient, and establish a general calculation model; The parameter calculation module is used to collect the measured hole pressure data of the thirteen-hole probe under the to-be-measured flow field, locate the belonging partition of the measured hole pressure data, and call the angle characteristic curve of the corresponding partition in the general calculation model to calculate the 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