CN116358823A - High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method - Google Patents

High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method Download PDF

Info

Publication number
CN116358823A
CN116358823A CN202310604920.3A CN202310604920A CN116358823A CN 116358823 A CN116358823 A CN 116358823A CN 202310604920 A CN202310604920 A CN 202310604920A CN 116358823 A CN116358823 A CN 116358823A
Authority
CN
China
Prior art keywords
pulsation
uncertainty
wind tunnel
mass flow
speed wind
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.)
Granted
Application number
CN202310604920.3A
Other languages
Chinese (zh)
Other versions
CN116358823B (en
Inventor
杜钰锋
熊能
汪路路
吴琦
李聪健
高川
王伟仲
罗太元
余强
夏雨
郭旦平
林俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN202310604920.3A priority Critical patent/CN116358823B/en
Publication of CN116358823A publication Critical patent/CN116358823A/en
Application granted granted Critical
Publication of CN116358823B publication Critical patent/CN116358823B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/02—Wind tunnels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/16—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements
    • 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)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

The invention belongs to the field of experimental aerodynamics, and discloses a method for evaluating the uncertainty of mass flow and total temperature pulsation of free incoming flow of a high-speed wind tunnel. The invention relates to a method for evaluating the uncertainty of mass flow and total temperature pulsation of free incoming flow of a high-speed wind tunnel, which comprises the following steps: installing a hot wire probe; connecting a hot wire anemometer; a first high-speed wind tunnel measurement test; a second high-speed wind tunnel measurement test; calculating mass flow pulsation and total temperature pulsation; and calculating mass flow pulsation uncertainty and total temperature pulsation uncertainty. The uncertainty evaluation method for the free incoming flow mass flow and the total temperature pulsation of the high-speed wind tunnel can provide the uncertainty of the free incoming flow mass flow pulsation and the total temperature pulsation of the high-speed wind tunnel, so that the reliability of apparent dynamic flow field quality and disturbance modal parameter solving results is improved, the reliability of high-speed wind tunnel test results is further improved, and the method has engineering application value.

Description

High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method
Technical Field
The invention belongs to the field of experimental aerodynamics, and particularly relates to a method for evaluating the uncertainty of mass flow and total temperature pulsation of free incoming flow of a high-speed wind tunnel.
Background
It is well known that the accurate simulation of the flow phenomenon under the real flight condition by using the wind tunnel test is an initial and pursued target of experimental hydrodynamic development, especially a high-speed wind tunnel, and plays an important role in aerodynamic force evaluation and aerodynamic shape refinement design of large advanced aircrafts such as various airliners, transport planes and the like. However, the disturbance modes existing in the free incoming flow of the high-speed wind tunnel can influence the aerodynamics, acoustics, thermal property and other characteristics of the test model, and the caused wind tunnel test result errors can cause design deviation of aerodynamic parameters such as lift force, resistance coefficient and the like in the aircraft design process, so that the load capacity of the aircraft can be caused to have estimation errors, and the economy and the safety of the aircraft are severely restricted. Therefore, it is important to precisely and quantitatively evaluate the disturbance modes existing in the free incoming flow of the high-speed wind tunnel.
The hot wire anemometer is the most suitable testing means for finely measuring the disturbance modes existing in the free incoming flow of the high-speed wind tunnel due to the advantages of high frequency response, high sensitivity, good economy and the like. The mass flow pulsation and the total temperature pulsation of the free incoming flow of the high-speed wind tunnel can be obtained through the measurement of the hot wire anemometer, and the apparent dynamic flow field quality and the disturbance modal parameters of the free incoming flow can be obtained according to a series of fitting processing, so that the method is very important for the accuracy analysis of the test result of the high-speed wind tunnel. However, uncertainty evaluation on free incoming flow mass flow pulsation and total temperature pulsation of the high-speed wind tunnel is not developed at present, so that reliability of apparent dynamic flow field quality and disturbance modal parameter solving results cannot be evaluated, reliability evaluation of high-speed wind tunnel test results is limited, and safety and economy of developing an aircraft through the high-speed wind tunnel test are severely restricted.
Currently, development of a method for evaluating mass flow rate of free incoming flow and uncertainty of total temperature pulsation of high-speed wind tunnel is needed.
Disclosure of Invention
The invention aims to provide a method for evaluating the uncertainty of the mass flow rate of free incoming flow and the total temperature pulsation of a high-speed wind tunnel, which is used for overcoming the defects of the prior art.
The invention relates to a method for evaluating the uncertainty of mass flow and total temperature pulsation of free incoming flow of a high-speed wind tunnel, which comprises the following steps:
s10, installing a hot wire probe;
the one-dimensional hot wire probe is arranged on the supporting rod, is well clamped by the clamping mechanism and is arranged on the middle bracket of the high-speed wind tunnel;
s20, connecting a hot wire anemometer;
opening a constant-temperature hot wire anemometer, connecting a data acquisition card, and setting the sampling frequency and the sampling time of the hot wire anemometer;
s30, a first high-speed wind tunnel measurement test;
setting the overheat ratio of the hot wire anemometer
Figure SMS_1
Starting a high-speed wind tunnel, and setting the operation Mach number>
Figure SMS_2
Recording the total free incoming flow temperature of the high-speed wind tunnel>
Figure SMS_3
Synchronously recording output voltage of hot wire anemometer>
Figure SMS_4
And calculates the first output voltage ripple of the hot wire anemometer>
Figure SMS_5
；
S40, performing a second high-speed wind tunnel measurement test;
maintaining the operation parameters of the high-speed wind tunnel unchanged, and changing the overheat ratio of the hot wire anemometer
Figure SMS_6
Record the output voltage of the hot wire anemometer +.>
Figure SMS_7
And calculates the second output voltage ripple of the hot wire anemometer +.>
Figure SMS_8
；
S50, calculating mass flow pulsation
Figure SMS_9
And total temperature pulsation->
Figure SMS_10
；
Hot wire anemometer output voltage pulsation
Figure SMS_11
Mass flow pulsation->
Figure SMS_12
Total temperature pulsation->
Figure SMS_13
The relationship between them is as follows:
Figure SMS_14
，
in the above-mentioned method, the step of,Gthe expression is as follows, which is the sensitivity coefficient of the total temperature pulsation:
Figure SMS_15
，
in the above-mentioned method, the step of,
Figure SMS_18
for reference temperature coefficient of resistance->
Figure SMS_20
For reference resistance, ">
Figure SMS_22
Is a temperature recovery coefficient->
Figure SMS_17
Is cold resistance value, ">
Figure SMS_21
For calibration coefficient, the five parameters are free to flow in the high-speed wind tunnel at total temperature +.>
Figure SMS_23
The data obtained in S30 and S40 form the output voltage pulsation +.>
Figure SMS_24
Mass flow pulsation->
Figure SMS_16
Pulsation of total temperature
Figure SMS_19
The equation set between is as follows:
Figure SMS_25
；
in the above-described system of equations,
Figure SMS_26
for the overheat ratio of the hot-wire anemometer>
Figure SMS_27
The corresponding sensitivity coefficient of the total temperature pulsation,
Figure SMS_28
for the overheat ratio of the hot-wire anemometer>
Figure SMS_29
Solving the equation set to obtain the free incoming flow mass flow pulsation of the high-speed wind tunnel according to the corresponding total temperature pulsation sensitivity coefficient>
Figure SMS_30
Total temperature pulsation->
Figure SMS_31
The expression of (2) is as follows:
Figure SMS_32
；
s60, calculating mass flow pulsation uncertainty
Figure SMS_33
And total temperature pulsation uncertainty +.>
Figure SMS_34
；
Figure SMS_35
，
In the above-mentioned method, the step of,
Figure SMS_37
is->
Figure SMS_39
Uncertainty of->
Figure SMS_41
Is->
Figure SMS_38
Uncertainty of->
Figure SMS_40
Is that
Figure SMS_42
Uncertainty of->
Figure SMS_43
Is->
Figure SMS_36
Is not determined by the degree of uncertainty of (2).
Further, the calculated mass flow pulsation uncertainty and the total temperature pulsation uncertainty of the S60
Figure SMS_44
Uncertainty of individual items +.>
Figure SMS_45
The calculation method of (2) is as follows:
S61.
Figure SMS_46
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_47
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure SMS_48
the calculation method of (2) is as follows:
Figure SMS_49
，
in the method, in the process of the invention,
Figure SMS_51
is total temperature->
Figure SMS_53
Is equal to the total temperature transmission of the high-speed wind tunnelSensor measurement accuracy, < >>
Figure SMS_55
Is a overheat ratio->
Figure SMS_52
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure SMS_54
Corresponding coefficients
Figure SMS_56
And->
Figure SMS_57
Corresponding coefficient->
Figure SMS_50
The calculation method comprises the following steps:
Figure SMS_58
；
Figure SMS_59
corresponding coefficient->
Figure SMS_60
The calculation method comprises the following steps:
Figure SMS_61
；
S62.
Figure SMS_62
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_63
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure SMS_64
the calculation method of (2) is as follows:
Figure SMS_65
，
in the method, in the process of the invention,
Figure SMS_67
is total temperature->
Figure SMS_69
Is equal to the measurement accuracy of the total temperature sensor of the high-speed wind tunnel, +.>
Figure SMS_71
Is a overheat ratio->
Figure SMS_68
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure SMS_70
Corresponding coefficients
Figure SMS_72
And->
Figure SMS_73
Corresponding coefficient->
Figure SMS_66
The calculation method comprises the following steps:
Figure SMS_74
；
Figure SMS_75
corresponding coefficient->
Figure SMS_76
The calculation method comprises the following steps:
Figure SMS_77
；
S63.
Figure SMS_78
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_79
Is calculated by the method;
Figure SMS_80
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_81
The calculation method of (2) is as follows:
Figure SMS_82
；
S64.
Figure SMS_83
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_84
Is calculated by the method;
Figure SMS_85
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_86
The calculation method of (2) is as follows:
Figure SMS_87
；
S65.
Figure SMS_88
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure SMS_89
Is calculated by the method;
Figure SMS_90
the calculation method has been given in S61, corresponding coefficients +.>
Figure SMS_91
The calculation method of (2) is as follows:
Figure SMS_92
；
S66.
Figure SMS_93
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure SMS_94
Is calculated by the method;
Figure SMS_95
the calculation method is given in S62, corresponding coefficient +.>
Figure SMS_96
The calculation method of (2) is as follows: />
Figure SMS_97
Figure SMS_98
；
S67.
Figure SMS_99
And its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_100
Is calculated by the method;
Figure SMS_101
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_102
The calculation method of (2) is as follows:
Figure SMS_103
；
S68.
Figure SMS_104
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_105
Is calculated by the method;
Figure SMS_106
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_107
The calculation method of (2) is as follows:
Figure SMS_108
。
the uncertainty evaluation method for the free incoming flow mass flow and the total temperature pulsation of the high-speed wind tunnel can provide the uncertainty of the free incoming flow mass flow pulsation and the total temperature pulsation of the high-speed wind tunnel, so that the reliability of apparent dynamic flow field quality and disturbance modal parameter solving results is improved, the reliability of high-speed wind tunnel test results is further improved, and the method has engineering application value.
Drawings
FIG. 1 is a flow chart of a method for analyzing uncertainty of internal resistance measurement of a high-speed wind tunnel ventilation model.
Detailed Description
The invention is described in detail below with reference to the drawings and examples.
Example 1:
as shown in fig. 1, the method for evaluating the uncertainty of the mass flow rate and the total temperature pulsation of the free incoming flow of the high-speed wind tunnel in the embodiment comprises the following steps:
s10, installing a hot wire probe;
the one-dimensional hot wire probe is arranged on the supporting rod, is well clamped by the clamping mechanism and is arranged on the middle bracket of the high-speed wind tunnel;
s20, connecting a hot wire anemometer;
opening a constant-temperature hot wire anemometer, connecting a data acquisition card, and setting the sampling frequency and the sampling time of the hot wire anemometer;
s30, a first high-speed wind tunnel measurement test;
setting the overheat ratio of the hot wire anemometer
Figure SMS_109
Starting a high-speed wind tunnel, and setting the operation Mach number>
Figure SMS_110
Recording the total free incoming flow temperature of the high-speed wind tunnel>
Figure SMS_111
Synchronously recording output voltage of hot wire anemometer>
Figure SMS_112
And calculates the first output voltage ripple of the hot wire anemometer>
Figure SMS_113
；
S40, performing a second high-speed wind tunnel measurement test;
maintaining the operation parameters of the high-speed wind tunnel unchanged, and changing the overheat ratio of the hot wire anemometer
Figure SMS_114
Record the output voltage of the hot wire anemometer +.>
Figure SMS_115
And calculates the second output voltage ripple of the hot wire anemometer +.>
Figure SMS_116
；
S50, calculating mass flow pulsation
Figure SMS_117
And total temperature pulsation->
Figure SMS_118
；
Hot wire anemometer output voltage pulsation
Figure SMS_119
Mass flow pulsation->
Figure SMS_120
Total temperature pulsation->
Figure SMS_121
The relationship between them is as follows:
Figure SMS_122
，
in the above-mentioned method, the step of,Gthe expression is as follows, which is the sensitivity coefficient of the total temperature pulsation:
Figure SMS_123
，
in the above-mentioned method, the step of,
Figure SMS_125
for reference temperature coefficient of resistance->
Figure SMS_128
For reference resistance, ">
Figure SMS_130
Is a temperature recovery coefficient->
Figure SMS_124
Is cold resistance value, ">
Figure SMS_127
For calibration coefficients, the five parameters are all fixed values under the condition that the total temperature of free incoming flow of the high-speed wind tunnel is unchanged, and the data obtained in S30 and S40 form output voltage pulsation +.>
Figure SMS_129
Mass flow pulsation->
Figure SMS_131
Pulsation of total temperature
Figure SMS_126
The equation set between is as follows:
Figure SMS_132
；
in the above-described system of equations,
Figure SMS_133
for the overheat ratio of the hot-wire anemometer>
Figure SMS_134
The corresponding sensitivity coefficient of the total temperature pulsation,
Figure SMS_135
for the overheat ratio of the hot-wire anemometer>
Figure SMS_136
Solving the equation set to obtain the free incoming flow mass flow pulsation of the high-speed wind tunnel according to the corresponding total temperature pulsation sensitivity coefficient>
Figure SMS_137
Total temperature pulsation->
Figure SMS_138
The expression of (2) is as follows:
Figure SMS_139
；
s60, calculating mass flow pulsation uncertainty
Figure SMS_140
And total temperature pulsation uncertainty +.>
Figure SMS_141
；
Figure SMS_142
，
In the above-mentioned method, the step of,
Figure SMS_145
is->
Figure SMS_147
Uncertainty of->
Figure SMS_149
Is->
Figure SMS_144
Uncertainty of->
Figure SMS_146
Is->
Figure SMS_148
Uncertainty of->
Figure SMS_150
Is->
Figure SMS_143
Is not determined by the degree of uncertainty of (2).
Further, the calculated mass flow pulsation uncertainty of S60
Figure SMS_151
And total temperature pulsation uncertainty +.>
Figure SMS_152
The calculation method of each item uncertainty in (a) is as follows: />
Figure SMS_153
Figure SMS_154
S61.
Figure SMS_155
And its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_156
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure SMS_157
the calculation method of (2) is as follows:
Figure SMS_158
，/>
in the method, in the process of the invention,
Figure SMS_161
is total temperature->
Figure SMS_163
Is equal to the measurement accuracy of the total temperature sensor of the high-speed wind tunnel, +.>
Figure SMS_165
Is a overheat ratio->
Figure SMS_160
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure SMS_162
Corresponding coefficients
Figure SMS_164
And->
Figure SMS_166
Corresponding coefficient->
Figure SMS_159
The calculation method comprises the following steps:
Figure SMS_167
，
Figure SMS_168
corresponding coefficient->
Figure SMS_169
The calculation method comprises the following steps:
Figure SMS_170
；
S62.
Figure SMS_171
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_172
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure SMS_173
the calculation method of (2) is as follows:
Figure SMS_174
，
in the method, in the process of the invention,
Figure SMS_175
is total temperature->
Figure SMS_178
Is equal to the measurement accuracy of the total temperature sensor of the high-speed wind tunnel, +.>
Figure SMS_180
Is a overheat ratio->
Figure SMS_177
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure SMS_179
Corresponding coefficients
Figure SMS_181
And->
Figure SMS_182
Corresponding coefficient->
Figure SMS_176
The calculation method comprises the following steps:
Figure SMS_183
，/>
Figure SMS_184
corresponding coefficient->
Figure SMS_185
The calculation method comprises the following steps:
Figure SMS_186
；
S63.
Figure SMS_187
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_188
Is calculated by the method;
Figure SMS_189
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_190
The calculation method of (2) is as follows:
Figure SMS_191
；
S64.
Figure SMS_192
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_193
Is calculated by the method;
Figure SMS_194
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_195
The calculation method of (2) is as follows:
Figure SMS_196
；/>
S65.
Figure SMS_197
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure SMS_198
Is calculated by the method;
Figure SMS_199
the calculation method has been given in S61, corresponding coefficients +.>
Figure SMS_200
The calculation method of (2) is as follows:
Figure SMS_201
；
S66.
Figure SMS_202
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure SMS_203
Is calculated by the method;
Figure SMS_204
the calculation method is given in S62, corresponding coefficient +.>
Figure SMS_205
The calculation method of (2) is as follows:
Figure SMS_206
；
S67.
Figure SMS_207
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_208
Is calculated by the method;
Figure SMS_209
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_210
The calculation method of (2) is as follows: />
Figure SMS_211
；
S68.
Figure SMS_212
And its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure SMS_213
Is calculated by the method;
Figure SMS_214
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure SMS_215
The calculation method of (2) is as follows:
Figure SMS_216
。
although embodiments of the invention have been disclosed in the foregoing description and illustrated in the drawings, it will be understood by those skilled in the art that the present invention is not limited to the specific details and illustrations of features and steps set forth herein, and that all features of the invention disclosed, or steps of the method or process, except for mutually exclusive features and/or steps, may be combined in any manner without departing from the principles of the invention.

Claims (2)

1. The method for evaluating the uncertainty of the mass flow and the total temperature pulsation of the free incoming flow of the high-speed wind tunnel is characterized by comprising the following steps of:
s10, installing a hot wire probe;
the one-dimensional hot wire probe is arranged on the supporting rod, is well clamped by the clamping mechanism and is arranged on the middle bracket of the high-speed wind tunnel;
s20, connecting a hot wire anemometer;
opening a constant-temperature hot wire anemometer, connecting a data acquisition card, and setting the sampling frequency and the sampling time of the hot wire anemometer;
s30, a first high-speed wind tunnel measurement test;
setting the overheat ratio of the hot wire anemometer
Figure QLYQS_1
Starting a high-speed wind tunnel, and setting the operation Mach number>
Figure QLYQS_2
Recording the total free incoming flow temperature of the high-speed wind tunnel>
Figure QLYQS_3
Synchronously recording output voltage of hot wire anemometer>
Figure QLYQS_4
And calculates the first output voltage ripple of the hot wire anemometer>
Figure QLYQS_5
；
S40, performing a second high-speed wind tunnel measurement test;
maintaining the operation parameters of the high-speed wind tunnel unchanged, and changing the overheat ratio of the hot wire anemometer
Figure QLYQS_6
Record the output voltage of the hot wire anemometer +.>
Figure QLYQS_7
And calculates the second output voltage ripple of the hot wire anemometer +.>
Figure QLYQS_8
；
S50, calculating mass flow pulsation
Figure QLYQS_9
And total temperature pulsation->
Figure QLYQS_10
；
Hot wire anemometer output voltage pulsation
Figure QLYQS_11
Mass flow pulsation->
Figure QLYQS_12
Total temperature pulsation->
Figure QLYQS_13
The relationship between them is as follows:
Figure QLYQS_14
，
in the above-mentioned method, the step of,Gthe expression is as follows, which is the sensitivity coefficient of the total temperature pulsation:
Figure QLYQS_15
，
in the above-mentioned method, the step of,
Figure QLYQS_17
for reference temperature coefficient of resistance->
Figure QLYQS_21
For reference resistance, ">
Figure QLYQS_23
Is a temperature recovery coefficient->
Figure QLYQS_18
Is cold resistance value, ">
Figure QLYQS_20
For calibration coefficient, the five parameters are free to flow in the high-speed wind tunnel at total temperature +.>
Figure QLYQS_22
The data obtained in S30 and S40 form the output voltage pulsation +.>
Figure QLYQS_24
Mass flow pulsation->
Figure QLYQS_16
Pulsation of total temperature
Figure QLYQS_19
The equation set between is as follows:
Figure QLYQS_25
；
in the above-described system of equations,
Figure QLYQS_26
for the overheat ratio of the hot-wire anemometer>
Figure QLYQS_27
Corresponding toTotal temperature pulsation sensitivity coefficient,/->
Figure QLYQS_28
For the overheat ratio of the hot-wire anemometer>
Figure QLYQS_29
Solving the equation set to obtain the free incoming flow mass flow pulsation of the high-speed wind tunnel according to the corresponding total temperature pulsation sensitivity coefficient>
Figure QLYQS_30
Total temperature pulsation->
Figure QLYQS_31
The expression of (2) is as follows:
Figure QLYQS_32
；
s60, calculating mass flow pulsation uncertainty
Figure QLYQS_33
And total temperature pulsation uncertainty +.>
Figure QLYQS_34
；
Figure QLYQS_35
；
In the above-mentioned method, the step of,
Figure QLYQS_37
is->
Figure QLYQS_39
Uncertainty of->
Figure QLYQS_41
Is->
Figure QLYQS_38
Uncertainty of->
Figure QLYQS_40
Is->
Figure QLYQS_42
Uncertainty of->
Figure QLYQS_43
Is->
Figure QLYQS_36
Is not determined by the degree of uncertainty of (2).
2. The method for evaluating the uncertainty of the mass flow rate and the total temperature pulsation of the free incoming flow of the high-speed wind tunnel according to claim 1, wherein the uncertainty of the calculated mass flow rate pulsation of the S60 is
Figure QLYQS_44
And total temperature pulsation uncertainty
Figure QLYQS_45
The calculation method of each item uncertainty in (a) is as follows:
S61.
Figure QLYQS_46
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_47
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure QLYQS_48
the calculation method of (2) is as follows:
Figure QLYQS_49
，
in the method, in the process of the invention,
Figure QLYQS_52
is total temperature->
Figure QLYQS_54
Is equal to the measurement accuracy of the total temperature sensor of the high-speed wind tunnel, +.>
Figure QLYQS_56
Is a overheat ratio->
Figure QLYQS_51
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure QLYQS_53
Corresponding coefficient->
Figure QLYQS_55
And
Figure QLYQS_57
corresponding coefficient->
Figure QLYQS_50
The calculation method comprises the following steps:
Figure QLYQS_58
，
Figure QLYQS_59
corresponding coefficient->
Figure QLYQS_60
The calculation method comprises the following steps:
Figure QLYQS_61
；
S62.
Figure QLYQS_62
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_63
Is calculated by the method;
according to S50GIs represented by the formula (i),
Figure QLYQS_64
the calculation method of (2) is as follows:
Figure QLYQS_65
，
in the method, in the process of the invention,
Figure QLYQS_67
is total temperature->
Figure QLYQS_70
Is equal to the measurement accuracy of the total temperature sensor of the high-speed wind tunnel, +.>
Figure QLYQS_72
Is a overheat ratio->
Figure QLYQS_68
Uncertainty in time, given by the constant temperature hot wire anemometer parameters, < >>
Figure QLYQS_69
Corresponding coefficient->
Figure QLYQS_71
And->
Figure QLYQS_73
Corresponding coefficient->
Figure QLYQS_66
The calculation method comprises the following steps:
Figure QLYQS_74
；
Figure QLYQS_75
corresponding coefficient->
Figure QLYQS_76
The calculation method comprises the following steps:
Figure QLYQS_77
；
S63.
Figure QLYQS_78
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_79
Is calculated by the method;
Figure QLYQS_80
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure QLYQS_81
The calculation method of (2) is as follows:
Figure QLYQS_82
；
S64.
Figure QLYQS_83
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_84
Is calculated by the method;
Figure QLYQS_85
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure QLYQS_86
The calculation method of (2) is as follows:
Figure QLYQS_87
；
S65.
Figure QLYQS_88
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure QLYQS_89
Is calculated by the method;
Figure QLYQS_90
the calculation method has been given in S61, corresponding coefficients +.>
Figure QLYQS_91
The calculation method of (2) is as follows:
Figure QLYQS_92
；
S66.
Figure QLYQS_93
and its corresponding coefficient in the uncertainty of the total temperature pulsation +.>
Figure QLYQS_94
Is calculated by the method;
Figure QLYQS_95
the calculation method is given in S62, corresponding toCoefficient->
Figure QLYQS_96
The calculation method of (2) is as follows:
Figure QLYQS_97
；
S67.
Figure QLYQS_98
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_99
Is calculated by the method;
Figure QLYQS_100
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure QLYQS_101
The calculation method of (2) is as follows:
Figure QLYQS_102
；
S68.
Figure QLYQS_103
and its corresponding coefficient in mass flow pulsation uncertainty +.>
Figure QLYQS_104
Is calculated by the method;
Figure QLYQS_105
the constant temperature hot wire anemometer output voltage parameter is given, corresponding coefficient +>
Figure QLYQS_106
The calculation method of (2) is as follows:
Figure QLYQS_107
。
CN202310604920.3A 2023-05-26 2023-05-26 High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method Active CN116358823B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310604920.3A CN116358823B (en) 2023-05-26 2023-05-26 High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310604920.3A CN116358823B (en) 2023-05-26 2023-05-26 High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method

Publications (2)

Publication Number Publication Date
CN116358823A true CN116358823A (en) 2023-06-30
CN116358823B CN116358823B (en) 2023-08-18

Family

ID=86922410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310604920.3A Active CN116358823B (en) 2023-05-26 2023-05-26 High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method

Country Status (1)

Country Link
CN (1) CN116358823B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117949692A (en) * 2024-03-27 2024-04-30 江西省气象探测中心 Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09159571A (en) * 1995-12-08 1997-06-20 Mitsubishi Heavy Ind Ltd Wind direction variation device for wind tunnel
DE19933665A1 (en) * 1999-07-17 2001-01-18 Bosch Gmbh Robert Device for detecting a pulsating variable
RU2002132839A (en) * 2002-12-05 2004-06-10 Российский государственный университет нефти и газа им. И.М. Губкина METHOD FOR DETERMINING THE PULSING GAS FLOW CONSUMPTION AT THE MAIN GAS PIPELINE
JP2007046924A (en) * 2005-08-08 2007-02-22 Hitachi Ltd Wind tunnel system and pulsation flow producing device
FR2942038A1 (en) * 2009-02-09 2010-08-13 Airbus France Aerodynamic coefficient's rigid value determining method for aircraft e.g. A320 aircraft, involves determining experimental value of aerodynamic coefficient, and determining rigid value of aerodynamic coefficient
JP3181789U (en) * 2012-12-10 2013-02-21 正裕 岩永 Hot-wire anemometer verification device
CA2950032A1 (en) * 2014-05-23 2015-11-26 Centre National De La Recherche Scientifique (Cnrs) Device and method for measuring a physical magnitude of a fluid flow
KR20170112005A (en) * 2016-03-30 2017-10-12 한국과학기술원 A measurement method of dynamic stability derivatives of a flight vehicle by utilizing magnetic force
JP2018021770A (en) * 2016-08-01 2018-02-08 東洋ゴム工業株式会社 Wind-tunnel test device
RU2697918C1 (en) * 2018-12-30 2019-08-21 Федеральное государственное бюджетное образовательное учреждение высшего образования "Пермский государственный национальный исследовательский университет" Method for measuring acoustic pulsations of a gas stream
CN110516342A (en) * 2019-08-22 2019-11-29 北京理工大学 A kind of compressible cavitating flows Numerical Predicting Method of propeller based on OpenFOAM platform
CN210426948U (en) * 2019-09-25 2020-04-28 中国空气动力研究与发展中心高速空气动力研究所 Hypersonic wind tunnel model bottom pressure measuring device
CN111624363A (en) * 2020-05-27 2020-09-04 中国空气动力研究与发展中心高速空气动力研究所 Hot wire anemometer calibration method considering influence of fluid compressibility
CN111766039A (en) * 2020-07-10 2020-10-13 中国空气动力研究与发展中心高速空气动力研究所 Method for calculating measurement result of compressible fluid disturbance mode of subsonic wind tunnel
CN113418675A (en) * 2021-08-06 2021-09-21 中国空气动力研究与发展中心设备设计与测试技术研究所 Hot wire measurement wind tunnel flow field disturbance modal method
CN115575074A (en) * 2022-12-09 2023-01-06 中国空气动力研究与发展中心高速空气动力研究所 High-speed wind tunnel ventilation model internal resistance measurement uncertainty analysis method
CN115788936A (en) * 2022-12-29 2023-03-14 沈阳鼓风机集团股份有限公司 Aerodynamic force analysis method of vaneless diffuser

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09159571A (en) * 1995-12-08 1997-06-20 Mitsubishi Heavy Ind Ltd Wind direction variation device for wind tunnel
DE19933665A1 (en) * 1999-07-17 2001-01-18 Bosch Gmbh Robert Device for detecting a pulsating variable
RU2002132839A (en) * 2002-12-05 2004-06-10 Российский государственный университет нефти и газа им. И.М. Губкина METHOD FOR DETERMINING THE PULSING GAS FLOW CONSUMPTION AT THE MAIN GAS PIPELINE
JP2007046924A (en) * 2005-08-08 2007-02-22 Hitachi Ltd Wind tunnel system and pulsation flow producing device
FR2942038A1 (en) * 2009-02-09 2010-08-13 Airbus France Aerodynamic coefficient's rigid value determining method for aircraft e.g. A320 aircraft, involves determining experimental value of aerodynamic coefficient, and determining rigid value of aerodynamic coefficient
JP3181789U (en) * 2012-12-10 2013-02-21 正裕 岩永 Hot-wire anemometer verification device
CA2950032A1 (en) * 2014-05-23 2015-11-26 Centre National De La Recherche Scientifique (Cnrs) Device and method for measuring a physical magnitude of a fluid flow
KR20170112005A (en) * 2016-03-30 2017-10-12 한국과학기술원 A measurement method of dynamic stability derivatives of a flight vehicle by utilizing magnetic force
JP2018021770A (en) * 2016-08-01 2018-02-08 東洋ゴム工業株式会社 Wind-tunnel test device
RU2697918C1 (en) * 2018-12-30 2019-08-21 Федеральное государственное бюджетное образовательное учреждение высшего образования "Пермский государственный национальный исследовательский университет" Method for measuring acoustic pulsations of a gas stream
CN110516342A (en) * 2019-08-22 2019-11-29 北京理工大学 A kind of compressible cavitating flows Numerical Predicting Method of propeller based on OpenFOAM platform
CN210426948U (en) * 2019-09-25 2020-04-28 中国空气动力研究与发展中心高速空气动力研究所 Hypersonic wind tunnel model bottom pressure measuring device
CN111624363A (en) * 2020-05-27 2020-09-04 中国空气动力研究与发展中心高速空气动力研究所 Hot wire anemometer calibration method considering influence of fluid compressibility
CN111766039A (en) * 2020-07-10 2020-10-13 中国空气动力研究与发展中心高速空气动力研究所 Method for calculating measurement result of compressible fluid disturbance mode of subsonic wind tunnel
CN113418675A (en) * 2021-08-06 2021-09-21 中国空气动力研究与发展中心设备设计与测试技术研究所 Hot wire measurement wind tunnel flow field disturbance modal method
CN115575074A (en) * 2022-12-09 2023-01-06 中国空气动力研究与发展中心高速空气动力研究所 High-speed wind tunnel ventilation model internal resistance measurement uncertainty analysis method
CN115788936A (en) * 2022-12-29 2023-03-14 沈阳鼓风机集团股份有限公司 Aerodynamic force analysis method of vaneless diffuser

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
KRAETZIG B ET AL: "temperature and heat flux measurement on hot models in short-duration facilities", JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER, vol. 29, no. 1, pages 37 - 46 *
朱博 等: "定/变热线过热比跨超声速流场湍流度测量", 航空动力学报, vol. 37, no. 9 *
杜钰锋 等: "亚声速风洞可压缩流体扰动模态分析", 航空学报, vol. 42, no. 6 *
杜钰锋 等: "变热线过热比可压缩流湍流度测量方法优化", 航空学报, vol. 40, no. 12 *
黄勇, 钱丰学, 董立新: "基于符号计算的风洞试验测量不确定度评估", 流体力学实验与测量, no. 03 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117949692A (en) * 2024-03-27 2024-04-30 江西省气象探测中心 Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer
CN117949692B (en) * 2024-03-27 2024-05-31 江西省气象探测中心 Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer

Also Published As

Publication number Publication date
CN116358823B (en) 2023-08-18

Similar Documents

Publication Publication Date Title
CN112304563B (en) Wind tunnel test method for researching influence of transition on aerodynamic characteristics of hypersonic aircraft
CN116358823B (en) High-speed wind tunnel free incoming flow mass flow and total temperature pulsation uncertainty evaluation method
Krishnamoorthy et al. Temperature-dissipation measurements in a turbulent boundary layer
CN111624363B (en) Hot wire anemometer calibration method considering influence of fluid compressibility
CN113639856B (en) Beam bridge modal frequency identification method considering environmental temperature influence
CN111707439B (en) Hyperbolic fitting method for compressible fluid turbulence measurement test data
CN110362902B (en) Single-source dynamic load identification method based on interval dimension-by-dimension analysis
CN103940577B (en) A kind of appraisal procedure of wind tunnel test balance based on acceleration signal energy
Jamroz Relationship between dynamic coefficients of two temperature sensors under nonstationary flow conditions
CN116989972A (en) Method for evaluating uncertainty of heat measurement test data based on thermoelectric simulation network test
CN115575074B (en) High-speed wind tunnel ventilation model internal resistance measurement uncertainty analysis method
CN114813096B (en) Multi-layer sintering net test model selection method for hypersonic wind tunnel
CN116181425A (en) Precision compensation method for high-temperature dynamic stress test of turbine rotor blade of engine
CN110044635A (en) A kind of thermal balance determination method based on real-time data analysis
CN117271938A (en) Method, system, device and medium for calculating relative error of test performance parameters
CN115435894A (en) Blade tip timing vibration stress inversion method based on simulated annealing algorithm
RU2714612C1 (en) Method of identifying nonlinear systems
Poljak Uncertainty Analysis of the NASA Glenn 8-by 6-Foot Supersonic Wind Tunnel 2019 Characterization Test 14-Foot, 5.8 Percent Test Section Porosity
John et al. Test bed calibration by trend analysis for reliability of a turboshaft engine performance
CN118225278B (en) Batch test method for temperature sensor chips
CN110261432A (en) Non-embedded multilayer materials thermal coefficient steady state measurement method
Finke Temperature-Compensated Calibration and Measurements with Cross-Style Hotwire Anemometers
CN114528670B (en) Method for detecting tensile strength of casting
BENNETT Measurements of periodic flow in rotating machinery
US8040958B2 (en) Method for measuring correlation between frequency response functions

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