CN116558771B - Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility - Google Patents
Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility Download PDFInfo
- Publication number
- CN116558771B CN116558771B CN202310849791.4A CN202310849791A CN116558771B CN 116558771 B CN116558771 B CN 116558771B CN 202310849791 A CN202310849791 A CN 202310849791A CN 116558771 B CN116558771 B CN 116558771B
- Authority
- CN
- China
- Prior art keywords
- flow field
- pressure
- measured
- mach number
- calibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000523 sample Substances 0.000 title claims abstract description 95
- 238000005259 measurement Methods 0.000 title claims abstract description 28
- 238000011156 evaluation Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 5
- 238000011088 calibration curve Methods 0.000 claims description 51
- 230000003068 static effect Effects 0.000 claims description 46
- 238000009530 blood pressure measurement Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000003491 array Methods 0.000 claims description 6
- 101100391172 Dictyostelium discoideum forA gene Proteins 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L27/00—Testing or calibrating of apparatus for measuring fluid pressure
- G01L27/002—Calibrating, i.e. establishing true relation between transducer output value and value to be measured, zeroing, linearising or span error determination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
The invention belongs to the field of experimental aerodynamics, and discloses a five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility. The method for evaluating the measurement uncertainty of the five-hole probe spherical coordinate system with consideration of compressibility comprises the following steps of calibrating the ground of the five-hole probe; a five-hole probe measurement test of a flow field to be measured; solving parameters of a flow field to be measured; and (5) evaluating the uncertainty of the flow field parameters to be tested. According to the method for evaluating the measurement uncertainty of the five-hole probe spherical coordinate system with consideration of compressibility, provided by the invention, the uncertainty evaluation method of the measurement result in the five-hole probe spherical coordinate system based on the Monte Carlo simulation method is established, the influence of compressibility caused by different Mach numbers is taken into consideration, and the method is provided for the refinement evaluation of the measurement result of the five-hole probe based on the spherical coordinate system in high-speed flow, so that the method has engineering application value.
Description
Technical Field
The invention belongs to the field of experimental aerodynamics, and particularly relates to a five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility.
Background
The air inlet channel, the engine, the spray pipe and the like of the aircraft propulsion system have compact internal structures, narrow space, complex flow structure and difficult flow fine measurement. The five-hole probe has the advantages of high measurement precision, high result reliability, convenience in use and the like, is a main means for measuring the flow direction of the flow points in the prior aeroengine, has very wide application space and value, can simultaneously obtain the key flow field information such as the flow speed and direction of the flow points, total pressure, static pressure and the like, and provides a powerful test means for analyzing the flow structure of the complex flow field.
Because the pressure measuring holes of the five-hole probes cannot ensure ideal consistency and symmetry in the processing and manufacturing process, each five-hole probe has unique mechanical construction characteristics and pneumatic characteristics, the response characteristics of the five-hole probes cannot be described by using the completely same mathematical model, and the probe calibration test is a key link for obtaining the response function of the five-hole probes and ensuring that the five-hole probes obtain correct flow field data.
At present, the available literature data are all calibration and application methods for measuring data of the five-hole probe in a rectangular coordinate system, namely, the pitch angle and the yaw angle of the probe are changed to obtain response data of the probe in different postures, but the influence of compressibility is not considered. The calibration data of the five-hole probe provided by the American Aeroprobe company as a five-hole probe design, manufacturing, calibration and application tap enterprise is given based on a spherical coordinate system, and because the calibration and application method of the measurement data under the spherical coordinate system is not disclosed, a five-hole probe user can only use the five-hole probe calibrated by paying a program of the Aeroprobe company for high-volume rent, thereby severely restricting the technical development of the five-hole probe and causing the waste of scientific research expenses to a certain extent. Moreover, the existing temporary literature data provide a method for evaluating the uncertainty of the measurement result of the five-hole probe spherical coordinate system, so that the fine evaluation of the measurement result of the five-hole probe spherical coordinate system is severely restricted.
Currently, there is a need to develop a five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility.
Disclosure of Invention
The invention aims to provide a five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility.
The five-hole probe used in the five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility is a circular tube type probe, the number of a central pressure hole is 1, the number of a right side pressure hole is 2, the number of a left side pressure hole is 3, the number of a lower pressure hole is 4, and the number of an upper pressure hole is 5.
The five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility comprises the following steps:
s10, calibrating the ground of the five-hole probe;
s11, performing a calibration test on the five-hole probe by using ground calibration equipment, wherein the ground calibration equipment has a polar angleAnd azimuth->Mechanism (S)>，/>The method comprises the steps of carrying out a first treatment on the surface of the Obtaining a five-well probe calibration data set including Mach number->Polar angle->Azimuth angle->Total pressure->Static pressure->And five-hole probe five pressure values corresponding to five pressure holes +.>；
S12 at each Mach numberRespectively calculating the following calibration coefficients to obtain each Mach number +.>Calibration coefficient set for the next five-well probe:
polar angle calibration coefficient:，
azimuth calibration coefficient:，
total pressure calibration coefficient:，
static pressure calibration coefficient:，
wherein ,，
based on each Mach number obtainedThe calibration coefficient set of the lower five-well probe, obtain each Mach number +.>A lower five-hole probe polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total pressure characteristic calibration curve and a static pressure characteristic calibration curve;
s20, a five-hole probe measurement test of a flow field to be measured;
placing the calibrated five-hole probe in a flow field to be measured to obtain pressure measurement values of five pressure holes；
S30, solving parameters of a flow field to be measured;
s31, calculating to obtain polar angle calibration coefficientsAnd azimuth alignment coefficient->；
According to the five-hole probe pressure measurement value obtained in step S20Calculating to obtain polar angle calibration coefficient +.>Azimuth calibration coefficient->：
，
wherein ,；
s32, setting a Mach number iteration initial value and calculation accuracy of a flow field to be detected;
assuming that the flow field Mach number iteration value is，nFor iteration count variable, initial valuen=0, mach number iterative computation accuracy is +.>；
S33, obtaining Mach number iteration value through interpolationA lower calibration curve;
each Mach number obtained according to step S12Under five-hole probe polar angle characteristic calibration curve, azimuth angle characteristic calibration curve, total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain Mach number iteration value +.>A lower polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total voltage characteristic calibration curve, and a static voltage characteristic calibration curve;
s34, obtaining a polar angle iteration value through interpolationAnd azimuth iteration value->；
The Mach number iteration value obtained according to step S33The lower polar angle characteristic calibration curve and the azimuth angle characteristic calibration curve are interpolated to obtain polar angle calibration coefficients +.>And azimuth alignment coefficient->Corresponding polar angle iteration value +.>And azimuth iteration value->；
S35, interpolation is carried out to obtain a total pressure calibration coefficient iteration valueIterative value to static pressure calibration coefficient>；
The Mach number iteration value obtained according to step S33Under total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain polar angle iteration value +.>And azimuth iteration value->Corresponding total pressure calibration coefficient iteration value +.>Iterative value to static pressure calibration coefficient>；
In the above, the corner mark beltnIs the iteration value of the corresponding variable,for the total pressure calibration coefficient variable, +.>The iteration value of the total pressure calibration coefficient is obtained; />For the static pressure calibration coefficient variable, then->The static pressure calibration coefficient iteration value;
s36, calculating to obtain Mach number iteration value；
Calculating the iteration value of the total pressure calibration coefficient according to the isentropic relationAnd static pressure calibration coefficient iteration value->Corresponding Mach number iteration value->：
；
S37, carrying out iterative computation to obtain flow field parameters to be detected;
judgingAnd->Is a size relationship of (a): if->Then the calculation is finished, the Mach number of the flow field to be measuredThe method comprises the steps of carrying out a first treatment on the surface of the If->Let->Repeating the steps S33-S36 until the conditions are satisfiedThe method comprises the steps of carrying out a first treatment on the surface of the After the calculation is finished, the Mach number of the flow field to be measured is +.>Polar angle +.>Azimuth angle +.>Total pressure of flow field to be measured->Static pressure +.>；
S40, evaluating uncertainty of parameters of the flow field to be tested;
evaluating the uncertainty of the flow field parameters to be tested by using a Monte Carlo simulation method;
s41, constructing a probability density model of the pressure measurement value of the five-hole probe;
pressure measurement value of five-hole probe in flow field to be measuredUncertainty of +.>By measuring range of pressure scanning valve>And accuracy determination;
the method comprises the steps of carrying out a first treatment on the surface of the Five-well probe the five-well pressure measurements were in agreement: mean value is measured value->The standard deviation is normal distribution corresponding to 1/3 of uncertainty, namely, the standard deviation accords with the 3 sigma principle;
s42, randomly sampling a probability density model;
generating random arrays of five-hole pressure measurement values of five-hole probes conforming to normal distribution in step S41 by utilizing MATLAB software, wherein the data volume of each array is as followsNFive random arrays are respectively recorded as，/>The data amount of the random array of the pressure measurement value of each hole of the five-hole probe which accords with the normal distribution in the step S41; the method comprises the steps of carrying out a first treatment on the surface of the
S43, calculating parameters of a flow field to be measured;
to be used forAs a five-hole probe, measuring the pressure of five pressure holes in a flow field to be measured; repeating the step S30 to obtain the Mach number array +.>Polar angle array->Azimuth array->Total pressure array->Static pressure array；
S44, evaluating uncertainty of parameters of the flow field to be tested;
analyzing the statistical characteristics of the flow field parameter array to be measured obtained in the step S43 to obtain the standard deviation of the flow field parameter array to be measured:
mach number array standard deviation of flow field to be measured，
Polar angle array standard deviation，
Standard deviation of azimuth angle array，
Total pressure array standard deviation，
Static pressure array standard deviation，
According to the principle of 3 sigma, the uncertainties of parameters of the flow field to be measured are respectively:
mach number uncertainty of flow field to be measured，
Polar angle uncertainty of flow field to be measured，
Uncertainty of azimuth angle of flow field to be measured，
Uncertainty of total pressure of flow field to be measured，
Static pressure uncertainty of flow field to be measured。
According to the method for evaluating the measurement uncertainty of the five-hole probe spherical coordinate system with consideration of compressibility, provided by the invention, the uncertainty evaluation method of the measurement result in the five-hole probe spherical coordinate system based on the Monte Carlo simulation method is established, the influence of compressibility caused by different Mach numbers is taken into consideration, and the method is provided for the refinement evaluation of the measurement result of the five-hole probe based on the spherical coordinate system in high-speed flow, so that the method has engineering application value.
Drawings
FIG. 1 is a schematic diagram of the numbering of five pressure measurement holes of a five-hole probe used in the method for evaluating the measurement uncertainty of a five-hole probe ball coordinate system taking compressibility into consideration;
fig. 2 is a flowchart of a five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility according to the present invention.
Detailed Description
The invention is described in detail below with reference to the drawings and examples.
Example 1:
as shown in fig. 1, the five-hole probe used in the five-hole probe spherical coordinate data processing method considering compressibility according to the present embodiment is a round tube probe, the number of the center pressure hole is 1, the number of the right pressure hole is 2, the number of the left pressure hole is 3, the number of the lower pressure hole is 4, and the number of the upper pressure hole is 5. As shown in fig. 2, the method for processing five-hole probe spherical coordinate data considering compressibility in this embodiment includes the following steps:
s10, calibrating the ground of the five-hole probe;
s11, performing a calibration test on the five-hole probe by using ground calibration equipment, wherein the ground calibration equipment has a polar angleAnd azimuth->Mechanism (S)>，/>The method comprises the steps of carrying out a first treatment on the surface of the Obtaining a five-well probe calibration data set including Mach number->Polar angle->Azimuth angle->Total pressure->Static pressure->Five-hole probe five pressure holesCorresponding pressure value->；
S12 at each Mach numberRespectively calculating the following calibration coefficients to obtain each Mach number +.>Calibration coefficient set for the next five-well probe:
polar angle calibration coefficient:，
azimuth calibration coefficient:，
total pressure calibration coefficient:，
static pressure calibration coefficient:，
wherein ,，
based on each Mach number obtainedThe calibration coefficient set of the lower five-well probe, obtain each Mach number +.>A lower five-hole probe polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total pressure characteristic calibration curve and a static pressure characteristic calibration curve;
s20, a five-hole probe measurement test of a flow field to be measured;
placing the calibrated five-hole probeIs arranged in a flow field to be measured to obtain pressure measurement values of five pressure holes；
S30, solving parameters of a flow field to be measured;
s31, calculating to obtain polar angle calibration coefficientsAnd azimuth alignment coefficient->；
According to the five-hole probe pressure measurement value obtained in step S20Calculating to obtain polar angle calibration coefficient +.>Azimuth calibration coefficient->：
，
wherein ,；
s32, setting a Mach number iteration initial value and calculation accuracy of a flow field to be detected;
assuming that the flow field Mach number iteration value is，nFor iteration count variable, initial valuen=0, mach number iterative computation accuracy is +.>；
S33, obtaining Mach number iteration value through interpolationA lower calibration curve;
each Mach number obtained according to step S12Under five-hole probe polar angle characteristic calibration curve, azimuth angle characteristic calibration curve, total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain Mach number iteration value +.>A lower polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total voltage characteristic calibration curve, and a static voltage characteristic calibration curve;
s34, obtaining a polar angle iteration value through interpolationAnd azimuth iteration value->；
The Mach number iteration value obtained according to step S33The lower polar angle characteristic calibration curve and the azimuth angle characteristic calibration curve are interpolated to obtain polar angle calibration coefficients +.>And azimuth alignment coefficient->Corresponding polar angle iteration value +.>And azimuth iteration value->；
S35, interpolation is carried out to obtain a total pressure calibration coefficient iteration valueIteration with static pressure calibration coefficientsValue->；
The Mach number iteration value obtained according to step S33Under total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain polar angle iteration value +.>And azimuth iteration value->Corresponding total pressure calibration coefficient iteration value +.>Iterative value to static pressure calibration coefficient>；
In the above, the corner mark beltnIs the iteration value of the corresponding variable,for the total pressure calibration coefficient variable, +.>The iteration value of the total pressure calibration coefficient is obtained; />For the static pressure calibration coefficient variable, then->The static pressure calibration coefficient iteration value;
s36, calculating to obtain Mach number iteration value；
Calculating the iteration value of the total pressure calibration coefficient according to the isentropic relationAnd static pressure calibration coefficient iteration value->Corresponding Mach number iteration value->：
；
S37, carrying out iterative computation to obtain flow field parameters to be detected;
judgingAnd->Is a size relationship of (a): if->Then the calculation is finished, the Mach number of the flow field to be measuredThe method comprises the steps of carrying out a first treatment on the surface of the If->Let->Repeating the steps S33-S36 until the conditions are satisfiedThe method comprises the steps of carrying out a first treatment on the surface of the After the calculation is finished, the Mach number of the flow field to be measured is +.>Polar angle +.>Azimuth angle +.>To be treated withTotal pressure of flow field>Static pressure +.>；
S40, evaluating uncertainty of parameters of the flow field to be tested;
evaluating the uncertainty of the flow field parameters to be tested by using a Monte Carlo simulation method;
s41, constructing a probability density model of the pressure measurement value of the five-hole probe;
pressure measurement value of five-hole probe in flow field to be measuredUncertainty of +.>By measuring range of pressure scanning valve>And accuracy determination;
the method comprises the steps of carrying out a first treatment on the surface of the Five-well probe the five-well pressure measurements were in agreement: mean value is measured value->The standard deviation is normal distribution corresponding to 1/3 of uncertainty, namely, the standard deviation accords with the 3 sigma principle;
s42, randomly sampling a probability density model;
generating random arrays of five-hole pressure measurement values of five-hole probes conforming to normal distribution in step S41 by utilizing MATLAB software, wherein the data volume of each array is as followsNFive random arrays are respectively recorded as，/>To conform to step S41Normally distributed five-hole probe data volume of random array of pressure measurement value of each hole
S43, calculating parameters of a flow field to be measured;
to be used forAs a five-hole probe, measuring the pressure of five pressure holes in a flow field to be measured; repeating the step S30 to obtain the Mach number array +.>Polar angle array->Azimuth array->Total pressure array->Static pressure array；
S44, evaluating uncertainty of parameters of the flow field to be tested;
analyzing the statistical characteristics of the flow field parameter array to be measured obtained in the step S43 to obtain the standard deviation of the flow field parameter array to be measured:
mach number array standard deviation of flow field to be measured，
Polar angle array standard deviation，
Standard deviation of azimuth angle array，
Total pressure array standard deviation，
Static pressure array standard deviation，
According to the principle of 3 sigma, the uncertainties of parameters of the flow field to be measured are respectively:
mach number uncertainty of flow field to be measured，
Polar angle uncertainty of flow field to be measured，
Uncertainty of azimuth angle of flow field to be measured，
Uncertainty of total pressure of flow field to be measured，
Static pressure uncertainty of flow field to be measured。
Although the embodiments of the present invention have been disclosed above, it is not limited to the use listed in the specification and the embodiments, but it can be fully applied to various fields suitable for the present invention. It will be apparent to one skilled in the art that the invention is not limited to the specific details and illustrations shown and described herein without departing from the general concepts defined in the claims and their equivalents.
Claims (1)
1. The five-hole probe ball coordinate system measurement uncertainty evaluation method considering compressibility is characterized in that the five-hole probe ball coordinate system measurement uncertainty evaluation method considering compressibility comprises the following steps:
s10, calibrating the ground of the five-hole probe;
s11, performing a calibration test on the five-hole probe by using ground calibration equipment, wherein the ground calibration equipment has a polar angleAnd azimuth angleMechanism (S)>，/>The method comprises the steps of carrying out a first treatment on the surface of the Obtaining a five-well probe calibration data set including Mach number->Polar angle->Azimuth angle->Total pressure->Static pressure->And five-hole probe five pressure values corresponding to five pressure holes +.>；
S12 at each Mach numberRespectively calculating the following calibration coefficients to obtain each Mach number +.>Calibration coefficient set for the next five-well probe:
polar angle calibration coefficient:，
azimuth calibration coefficient:，
total pressure calibration coefficient:，
static pressure calibration coefficient:，
wherein ,，
based on each Mach number obtainedThe calibration coefficient set of the lower five-well probe, obtain each Mach number +.>A lower five-hole probe polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total pressure characteristic calibration curve and a static pressure characteristic calibration curve;
s20, a five-hole probe measurement test of a flow field to be measured;
placing the calibrated five-hole probe in a flow field to be measured to obtain pressure measurement values of five pressure holes；
S30, solving parameters of a flow field to be measured;
s31, calculating to obtain polar angle calibration coefficientsAnd azimuth alignment coefficient->；
According to the five-hole probe pressure measurement value obtained in step S20Calculating to obtain polar angle calibration coefficient +.>Azimuth calibration coefficient->：
，
wherein ,；
s32, setting a Mach number iteration initial value and calculation accuracy of a flow field to be detected;
assuming that the flow field Mach number iteration value is，nFor iteration count variable, initial valuenThe Mach number iterative calculation accuracy is =0；
S33, obtaining Mach number iteration value through interpolationA lower calibration curve;
each Mach number obtained according to step S12Under five-hole probe polar angle characteristic calibration curve, azimuth angle characteristic calibration curve, total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain Mach number iteration value +.>A lower polar angle characteristic calibration curve, an azimuth angle characteristic calibration curve, a total voltage characteristic calibration curve, and a static voltage characteristic calibration curve;
s34, obtaining a polar angle iteration value through interpolationAnd azimuth iteration value->；
The Mach number iteration value obtained according to step S33The lower polar angle characteristic calibration curve and the azimuth angle characteristic calibration curve are interpolated to obtain polar angle calibration coefficients +.>And azimuth alignment coefficient->Corresponding polar angle iteration value +.>And azimuth iteration value->；
S35, interpolation is carried out to obtain a total pressure calibration coefficient iteration valueIterative value to static pressure calibration coefficient>；
The Mach number iteration value obtained according to step S33Under total pressure characteristic calibration curve and static pressure characteristic calibration curve, interpolation is carried out to obtain polar angle iteration value +.>And azimuth iteration value->Corresponding total pressure calibration coefficient iteration value +.>Iterative value to static pressure calibration coefficient>；
In the above, the corner mark beltnIs the iteration value of the corresponding variable,for the total pressure calibration coefficient variable, +.>The iteration value of the total pressure calibration coefficient is obtained; />For the static pressure calibration coefficient variable, then->The static pressure calibration coefficient iteration value;
s36, calculating to obtain Mach number iteration value；
Calculating the total pressure calibration system according to the isentropic relationNumerical iteration valueAnd static pressure calibration coefficient iteration value->Corresponding Mach number iteration value->：
；
S37, carrying out iterative computation to obtain flow field parameters to be detected;
judgingAnd->Is a size relationship of (a): if->Then the calculation is finished, the Mach number of the flow field to be measuredThe method comprises the steps of carrying out a first treatment on the surface of the If->Let->Repeating the steps S33-S36 until the conditions are satisfiedThe method comprises the steps of carrying out a first treatment on the surface of the After the calculation is finished, the Mach number of the flow field to be measured is +.>Polar angle +.>Azimuth angle +.>Total pressure of flow field to be measured->Static pressure +.>；
S40, evaluating uncertainty of parameters of the flow field to be tested;
evaluating the uncertainty of the flow field parameters to be tested by using a Monte Carlo simulation method;
s41, constructing a probability density model of the pressure measurement value of the five-hole probe;
pressure measurement value of five-hole probe in flow field to be measuredUncertainty of +.>By measuring range of pressure scanning valve>And accuracy determination;
the method comprises the steps of carrying out a first treatment on the surface of the Five-well probe the five-well pressure measurements were in agreement: mean value is measured value->The standard deviation is normal distribution corresponding to 1/3 of uncertainty, namely, the standard deviation accords with the 3 sigma principle;
s42, randomly sampling a probability density model;
generating five-hole pressure measurement of five-hole probes conforming to normal distribution in step S41 by utilizing MATLAB softwareRandom arrays of magnitudes, each array having a data size ofNFive random arrays are respectively recorded as，/>The data amount of the random array of the pressure measurement value of each hole of the five-hole probe which accords with the normal distribution in the step S41;
s43, calculating parameters of a flow field to be measured;
to be used forAs a five-hole probe, measuring the pressure of five pressure holes in a flow field to be measured; repeating the step S30 to obtain the Mach number array +.>Polar angle array->Azimuth array->Total pressure array->Static pressure array；
S44, evaluating uncertainty of parameters of the flow field to be tested;
analyzing the statistical characteristics of the flow field parameter array to be measured obtained in the step S43 to obtain the standard deviation of the flow field parameter array to be measured:
mach number array standard deviation of flow field to be measured，
Polar angle array standard deviation，
Standard deviation of azimuth angle array，
Total pressure array standard deviation，
Static pressure array standard deviation，
According to the principle of 3 sigma, the uncertainties of parameters of the flow field to be measured are respectively:
mach number uncertainty of flow field to be measured，
Polar angle uncertainty of flow field to be measured，
Uncertainty of azimuth angle of flow field to be measured，
Uncertainty of total pressure of flow field to be measured，
Static pressure uncertainty of flow field to be measured。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310849791.4A CN116558771B (en) | 2023-07-12 | 2023-07-12 | Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310849791.4A CN116558771B (en) | 2023-07-12 | 2023-07-12 | Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116558771A CN116558771A (en) | 2023-08-08 |
| CN116558771B true CN116558771B (en) | 2023-09-01 |
Family
ID=87498693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310849791.4A Active CN116558771B (en) | 2023-07-12 | 2023-07-12 | Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116558771B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110514391A (en) * | 2019-08-12 | 2019-11-29 | 北京航空航天大学 | A kind of five pore pressure force probes measurement three-dimensional flow field uncertainty evaluation method |
| CN110514390A (en) * | 2019-08-12 | 2019-11-29 | 北京航空航天大学 | A kind of three pore pressure force probes measurement two-dimensional flow field uncertainty evaluation method |
| CN114486164A (en) * | 2022-04-18 | 2022-05-13 | 中国空气动力研究与发展中心空天技术研究所 | Three-hole probe gradient flow field test calibration method based on neural network |
| CN115358026A (en) * | 2022-09-21 | 2022-11-18 | 大连温特纳科技有限公司 | Five-hole probe data processing method based on multivariate linear regression and surface fitting |
| CN116007884A (en) * | 2023-03-21 | 2023-04-25 | 中国空气动力研究与发展中心高速空气动力研究所 | Five-hole probe spherical coordinate data processing method considering compressibility |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100329081A1 (en) * | 2009-06-26 | 2010-12-30 | Eric Sullivan | Method for non-destructively evaluating rotary earth boring drill components and determining fitness-for-use of the same |
| US20230211873A1 (en) * | 2022-01-04 | 2023-07-06 | The Florida State University Research Foundation, Inc. | Method for Algorithmic Optimization of Active Flow Control Actuator Placement and Parameters |
-
2023
- 2023-07-12 CN CN202310849791.4A patent/CN116558771B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110514391A (en) * | 2019-08-12 | 2019-11-29 | 北京航空航天大学 | A kind of five pore pressure force probes measurement three-dimensional flow field uncertainty evaluation method |
| CN110514390A (en) * | 2019-08-12 | 2019-11-29 | 北京航空航天大学 | A kind of three pore pressure force probes measurement two-dimensional flow field uncertainty evaluation method |
| CN114486164A (en) * | 2022-04-18 | 2022-05-13 | 中国空气动力研究与发展中心空天技术研究所 | Three-hole probe gradient flow field test calibration method based on neural network |
| CN115358026A (en) * | 2022-09-21 | 2022-11-18 | 大连温特纳科技有限公司 | Five-hole probe data processing method based on multivariate linear regression and surface fitting |
| CN116007884A (en) * | 2023-03-21 | 2023-04-25 | 中国空气动力研究与发展中心高速空气动力研究所 | Five-hole probe spherical coordinate data processing method considering compressibility |
Non-Patent Citations (1)
| Title |
|---|
| 基于蒙特卡罗法的五孔探针测量不确定度评定;马宏伟 等;航空动力学报;第11卷(第37期);第2588-2591页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116558771A (en) | 2023-08-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116007884B (en) | Five-hole probe spherical coordinate data processing method considering compressibility | |
| CN109521444B (en) | Self-adaptive least square fitting estimation algorithm for GPS horizontal velocity field of crustal movement | |
| CN111190060A (en) | Antenna array fault diagnosis method considering array errors in impulse noise environment | |
| US9851441B2 (en) | Method and system for generating a distance velocity azimuth display | |
| CN115358026A (en) | Five-hole probe data processing method based on multivariate linear regression and surface fitting | |
| CN116558771B (en) | Five-hole probe sphere coordinate system measurement uncertainty evaluation method considering compressibility | |
| CN116147724B (en) | Metering method suitable for ultrasonic water meter | |
| CN109900789B (en) | Imaging method for internal defects of longitudinal section of tree | |
| CN113031514B (en) | R-test calibration uncertainty evaluation method based on metrology | |
| CN103077274A (en) | High-precision curve modeling intelligent method and device | |
| CN111811827B (en) | Product performance consistency inspection method based on Rayleigh distribution | |
| CN116558770B (en) | Five-hole probe measurement uncertainty evaluation method considering compressibility | |
| CN109859323B (en) | Method for weighting forest stand space pattern based on triangular network model | |
| CN117634938A (en) | Production quality detection method and system for seamless steel pipe | |
| CN109164439A (en) | Air index calculation method on a kind of electric wave ray | |
| CN115436874A (en) | Three-dimensional passive positioning method | |
| CN108008337A (en) | The computational methods and device of electric energy meter measurement error uniformity | |
| Wu et al. | Simulation and model validation of positional uncertainty of line feature on manual digitizing a map | |
| Milanovic et al. | Numerical calibration of a conical five-hole probe for supersonic measurements | |
| CN116558769B (en) | Five-hole probe calibration data application method considering compressibility | |
| CN107063218A (en) | Data processing method and device | |
| CN108398600B (en) | Impedance spectrum measuring method and tester based on Prony state space algorithm | |
| CN113155384A (en) | Sensor arrangement method for reducing uncertainty of structural damping ratio identification | |
| CN113158501A (en) | Sensor arrangement method for reducing structural natural frequency identification uncertainty | |
| CN115310028B (en) | Calculation method for uncertainty in measurement of included angle between axis of blind hole and plane |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |