CN113125103B - Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter - Google Patents
Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter Download PDFInfo
- Publication number
- CN113125103B CN113125103B CN202110314172.6A CN202110314172A CN113125103B CN 113125103 B CN113125103 B CN 113125103B CN 202110314172 A CN202110314172 A CN 202110314172A CN 113125103 B CN113125103 B CN 113125103B
- Authority
- CN
- China
- Prior art keywords
- section
- pressure
- pitot
- measuring point
- elliptical
- 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
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/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/065—Measuring arrangements specially adapted for aerodynamic testing dealing with flow
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a data processing method of an elliptic section flowmeter with 41 measuring points equidistantly distributed, which comprises the steps of measuring the pitot pressure of each measuring point of an elliptic section of an outlet of an air inlet channel and the static pressure of the section position through a flowmeter positioned at the downstream of the outlet of the elliptic section of the air inlet channel, obtaining the static pressure of each measuring point by utilizing a linear interpolation method, obtaining the specific flow density, the airflow Mach number and the total pressure recovery coefficient of each measuring point through calculation, respectively integrating and superposing the specific flow density, the airflow Mach number and the total pressure recovery coefficient of each measuring point in the elliptic section to obtain the flow, the superposed Mach number and the superposed total pressure recovery coefficient in the elliptic section, and respectively dividing the superposed Mach number and the superposed total pressure recovery coefficient in the elliptic section by the area of the elliptic section to obtain the average Mach number and the total pressure recovery coefficient in the elliptic section. The invention fills the blank of the post-processing method when the measuring points of the elliptical cross-section flow meter are distributed at equal intervals, and can more comprehensively obtain the performance parameters of the air inlet channel.
Description
Technical Field
The invention belongs to the technical field of experimental aerodynamics, and particularly relates to an elliptic cross-section flowmeter with 41 measuring points distributed equidistantly and a data processing method thereof.
Background
The air inlet channel is positioned at the most front end of the engine, captures incoming flow, decelerates and pressurizes the incoming flow to provide airflow for a downstream combustion chamber, and the working capacity of the air inlet channel is directly related to the working efficiency of the whole air-breathing propulsion system.
Parameters such as outlet flow of the air inlet, mach number, total pressure recovery coefficient and the like are key parameters influencing the design of the air inlet, the Mach number is the ratio of local airflow speed to sound velocity, and the total pressure recovery coefficient is the ratio of local airflow total pressure to incoming flow total pressure. The total pressure of the airflow is measured by installing a pressure measuring rake on the outlet section, the positions of measuring points on the leather supporting rake are limited by the size of a sensor and are usually distributed at equal intervals, pressure measuring holes are formed in the wall surface where the measuring section is located to measure the static pressure of the airflow, and parameters of the outlet section are obtained through data processing.
However, the current cross-sectional parameter measurement is directed at the air inlet with a rectangular cross section or a circular cross section, and with the development of the air inlet design technology, the air inlet with an elliptical outlet also becomes another common air inlet. For a rectangular cross-section air inlet channel, a measurement method of grid uniform distribution is often adopted; for a circular-section air inlet channel, the symmetry of a circle is utilized, and the average flow of an outlet is usually obtained by simply averaging pressure data of each measuring point. In the traditional method, by averaging the pressure data of each measuring point, the data of the outlet section cannot be represented, and the parameter distribution rule on the section cannot be obtained.
Disclosure of Invention
The invention aims to provide an elliptic section flowmeter with 41 measuring points distributed equidistantly and a data processing method thereof, so as to solve the technical problem that the flow, the Mach number and the total pressure recovery coefficient of an elliptic outlet of an air inlet channel are difficult to measure.
In order to achieve the purpose, the invention provides the following technical scheme:
the elliptical cross-section flow meter with the 41 measuring points equidistantly distributed comprises an elliptical shell, wherein the elliptical shell is connected with an elliptical cross-section outlet at the downstream of an air inlet, the elliptical shell is used for mounting a pitot harrow and a static pressure measuring device, 8 harrow positions are circumferentially arranged on the pitot harrow, the cross-section area of the elliptical outlet is divided into 8 parts by the 8 harrow positions, 1 pitot pressure measuring point is arranged at the central intersection point of the 8 harrow positions, 5 pitot pressure measuring points are uniformly arranged on each harrow position along the radial direction, pitot pressure sensors are mounted at the pitot pressure measuring points, the static pressure measuring device comprises 8 wall surface pressure sensors uniformly arranged on the elliptical shell along the circumferential direction, the position distribution of the wall surface pressure sensors corresponds to the harrow positions, the pitot pressure sensors are used for obtaining pitot pressure of the corresponding pitot pressure measuring points, and the wall surface pressure sensors are used for measuring static pressure of a flow field at the corresponding positions.
Another object of the present invention is to provide a data processing method for an oval cross-section flowmeter with equidistantly distributed 41 measuring points, which is applied to the oval cross-section flowmeter with equidistantly distributed 41 measuring points, and includes:
step 1: averaging the measurement results of the 8 wall surface pressure sensors to obtain the static pressure of the center of the elliptical section, and obtaining the static pressure of each pitot pressure measurement point by using a linear interpolation method along the radial direction according to the static pressure of the center of the elliptical section and the wall surface static pressure corresponding to the corresponding rake position;
and 2, step: calculating the Mach number of the airflow of each pitot pressure measuring point according to the pitot pressure and the static pressure of each pitot pressure measuring point;
and step 3: calculating to obtain the total airflow pressure of each pitot voltage measuring point according to the static pressure and the airflow Mach number of each pitot voltage measuring point;
and 4, step 4: calculating the specific flow density of each pitot pressure measuring point according to the airflow Mach number and the airflow total pressure of each pitot pressure measuring point;
and 5: calculating to obtain total pressure recovery coefficients of the pitot pressure measuring points according to the total airflow pressure and the total incoming flow pressure of the pitot pressure measuring points;
and 6: respectively integrating and superposing the specific flow density, the airflow Mach number and the total pressure recovery coefficient of each pitot pressure measuring point in the elliptical section to obtain the flow, the superposed Mach number and the superposed total pressure recovery coefficient in the elliptical section;
and 7: and respectively dividing the superposed Mach number and the superposed total pressure recovery coefficient in the elliptical section by the area of the elliptical section to obtain the average Mach number and the total pressure recovery coefficient in the elliptical section.
Preferably, the integrating and superposing the specific flow density, the airflow mach number and the total pressure recovery coefficient of each pitot pressure measuring point in the elliptical cross section to obtain the flow, the superposed mach number and the superposed total pressure recovery coefficient in the elliptical cross section includes:
step 61, a sector area is formed by the center of the elliptical section and two adjacent pitot pressure measurement points on the innermost layer, the specific flow density of each pitot pressure measurement point is interpolated by utilizing a linear interpolation method to obtain a specific flow density interpolation function of each sector area, and the interpolation function of the compared flow density is integrated in the sector area to obtain the flow in each sector area; the calculation formula of the interpolation function of the specific flow density in the sector area is as follows:
the calculation formula of the flow in the sector area is as follows:
and 62, interpolating the specific flow density of each pitot pressure measuring point by using a linear interpolation method in a first arc-shaped area consisting of two pitot pressure measuring points adjacent to the innermost layer and two pitot pressure measuring points adjacent to the corresponding outermost layer to obtain a specific flow density interpolation function of each first arc-shaped area, and integrating the interpolation function of the specific flow density in the first arc-shaped area to obtain the flow in each first arc-shaped area. Wherein the calculation formula of the interpolation function of the specific flow density in the first circular arc-shaped area is as follows:
the calculation formula of the flow in the first arc-shaped area is as follows:
and 63, approximating a second arc-shaped area consisting of two adjacent pitot pressure measuring points on the outermost layer and the corresponding elliptical shell by adopting the measuring points on the outermost layer along the radial direction, interpolating the specific flow density of each pitot pressure measuring point by adopting linear interpolation values along the circumferential direction to obtain a specific flow density interpolation function of each second arc-shaped area, and integrating the interpolation function of the specific flow density in the second arc-shaped area to obtain the flow in each second arc-shaped area. Wherein the calculation formula of the interpolation function of the specific flow density in the second circular arc-shaped area is as follows:
the calculation formula of the flow in the second arc-shaped area is as follows:
and step 64, superposing the flow in each sector area, the flow in each first circular arc area and the flow in each second circular arc area to obtain the flow in the elliptical section. Wherein, the calculation formula of the flow in the elliptical section is as follows:
step 65, replacing the specific flow density in the formulas (1), (3) and (5) with the Mach number of each pitot pressure measuring point, and repeating the steps 61-64 to obtain the superposed Mach number in the elliptical section;
and 66, replacing the specific flow density in the formulas (1), (3) and (5) with the total pressure recovery coefficient of each pitot pressure measuring point, and repeating the steps 61 to 64 to obtain the total pressure recovery coefficient superposed in the elliptical section.
By utilizing the technical scheme, the invention has the following beneficial technical effects:
the flow, the Mach number and the total pressure recovery coefficient of the outlet of the elliptic section of the air inlet channel can be obtained through the elliptic section flow meter and the data processing method thereof, the blank of the post-processing method when the measuring points of the elliptic section flow meter are equidistantly distributed is filled, so that the performance parameters of the air inlet channel can be more comprehensively obtained, and scientific basis can be provided for flow calculation of the elliptic section and further post-processing of the distribution of the Mach number and the total pressure recovery coefficient of the elliptic section, the average Mach number and the average total pressure recovery coefficient.
Drawings
FIG. 1 is a schematic structural diagram of an elliptic cross-section flowmeter with equally-spaced 41 measuring points according to an embodiment of the present invention;
fig. 2 is a flowchart of a data processing method of an elliptic cross section flowmeter with equidistantly distributed 41 measuring points according to an embodiment of the present invention.
Detailed Description
The invention is described in further detail below with reference to the following figures and embodiments:
as shown in figure 1, the flow meter with the elliptic cross sections and the equidistantly distributed 41 measuring points comprises an elliptic shell, wherein the elliptic shell is connected to an elliptic cross section outlet at the downstream of an air inlet, and the elliptic shell is used for mounting a pitot harrow and a static pressure measuring device.
For convenient calculation, a coordinate system is established by taking the center of the ellipse as the origin, the major axis of the ellipse as the x axis, the length of the ellipse as the 2a, the minor axis of the ellipse as the y axis and the length of the ellipse as the 2b, and the elliptic equation is
Let D ij For pitot pressure points (as shown in FIG. 1 for the solidThe center circle points are shown), i is the position of the sensor measuring point on each rake position i =1,2, \8230 \82305, which is increased from inside to outside along the radial direction, j is the rake position number j =1,2, \8230; \ 82308, which is increased from the x axis along the circumferential direction and is increased from the x axis along the counterclockwise direction, and D is the length of the sensor measuring point 00 Is a pitot pressure measuring point at the center of the ellipse. Let E j J =1,2, \8230 @, \ 8230:8, which are wall static pressure measuring points (shown by triangular dots in fig. 1), are increased counterclockwise from the x-axis in the circumferential direction and correspond to the numbers of the leather-supporting rakes.
The leather-supporting rake comprises 8 rake positions, the 8 rake positions divide the cross section area of the elliptic outlet into 8 equal parts, wherein 4 forward rake positions are distributed in a cross shape and coincide with the coordinate axis, and the included angle between the other 4 oblique rake positions and the elliptic major axis is theta 0 = arctan (b/a); 1 pitot pressure measuring point is arranged at the center intersection point of 8 rake positions, 5 pitot pressure measuring points are uniformly arranged on each rake position along the radial direction, and the total number of the pitot pressure measuring points is 41; and each pitot pressure measuring point is provided with a pitot pressure sensor, the static pressure measuring device comprises 8 wall surface pressure sensors which are uniformly arranged on the elliptical shell along the circumferential direction, the position distribution of the wall surface pressure sensors corresponds to the rake positions, the pitot pressure sensors are used for acquiring the pitot pressure of the corresponding pitot pressure measuring point, and the wall surface pressure sensors are used for measuring the static pressure of the flow field at the corresponding position.
Referring to fig. 2, a data processing method for an elliptic cross section flowmeter with equidistantly distributed 41 measuring points is applied to the elliptic cross section flowmeter with equidistantly distributed 41 measuring points, and includes:
step 1: averaging the measurement results of the 8 wall surface pressure sensors to obtain the static pressure of the center of the elliptical section, and obtaining the static pressure of each pitot pressure measurement point by utilizing a linear interpolation method along the radial direction according to the static pressure of the center of the elliptical section and the wall surface static pressure corresponding to the corresponding rake position;
in particular, through 8 measuring points E j The static pressure of (2) is averaged to calculate the center D of the ellipse 00 Static pressure p of 00 (ii) a For the measuring point Dij on the jth rake position, use D 00 At static pressure p 00 And E j At static pressure p E j Obtaining D by linear interpolation in the radial direction ij Static pressure p at (i =1,2, \ 8230; 5) ij As shown in the formula (8),
step 2: calculating the Mach number of the airflow of each pitot pressure measuring point according to the pitot pressure and the static pressure of each pitot pressure measuring point;
in particular, D is obtained by means of a Pitot-rake measurement ij Pressure p for skin t ij Calculating the ratio p of the skin pressure to the static pressure t ij /p ij And further calculating D by the formula (2) ij Mach number M of treated air flow ij ，
Where γ is the gas specific heat ratio and γ =1.4 for air.
And step 3: calculating to obtain the total airflow pressure of each pitot voltage measuring point according to the static pressure and the airflow Mach number of each pitot voltage measuring point;
in particular, by means of D ij At static pressure p ij And D ij Mach number M of treated air flow ij Calculating D by equation (3) ij Total pressure p of exhaust gas flow 0 ij ，
And 4, step 4: calculating the specific flow density of each pitot pressure measuring point according to the airflow Mach number and the airflow total pressure of each pitot pressure measuring point;
specifically, the specific flow density m is calculated A ij = ρ u, i.e. the product of the gas density ρ and the velocity u,
in the formula T 0 ij The total temperature of the gas flow can be approximated or measured by a thermocouple or the like.
And 5: calculating to obtain total pressure recovery coefficients of the pitot pressure measuring points according to the total airflow pressure and the total incoming flow pressure of the pitot pressure measuring points;
specifically, the total pressure recovery coefficient σ ij Is the total pressure p 0 ij And total pressure p of incoming flow 0 ∞ Ratio of σ ij ＝p 0 ij /p0 ∞ 。
Step 6: respectively integrating and superposing the specific flow density, the airflow Mach number and the airflow total pressure of each pitot pressure measuring point in the elliptical section to obtain the flow, the superposed Mach number and the superposed total pressure recovery coefficient in the elliptical section;
preferably, the step of integrating and superposing the specific flow density, the airflow mach number and the airflow total pressure of each pitot pressure measuring point in the elliptical cross section to obtain the flow, the superposed mach number and the superposed total pressure recovery coefficient in the elliptical cross section includes:
step 61, interpolating the specific flow density of each pitot pressure measuring point by utilizing a linear interpolation method in a sector area consisting of the center of the elliptical section and two adjacent pitot pressure measuring points on the innermost layer to obtain a specific flow density interpolation function of each sector area, and integrating the interpolation function of the compared flow density in the sector area to obtain the flow in each sector area;
step 62, interpolating the specific flow density of each pitot pressure measuring point by using a linear interpolation method in a first arc-shaped area consisting of two pitot pressure measuring points adjacent to the innermost layer and two pitot pressure measuring points adjacent to the corresponding outermost layer to obtain a specific flow density interpolation function of each first arc-shaped area, and integrating the interpolation function of the specific flow density in the first arc-shaped area to obtain the flow in each first arc-shaped area;
step 63, approximating the outermost measuring points in the radial direction by using the two adjacent pitot pressure measuring points on the outermost layer and a second arc-shaped area consisting of the corresponding elliptical shell, interpolating the specific flow density of each pitot pressure measuring point by using linear interpolation values in the circumferential direction to obtain a specific flow density interpolation function of each second arc-shaped area, and integrating the interpolation function of the specific flow density in the second arc-shaped area to obtain the flow in each second arc-shaped area;
step 64, superposing the flow in each sector area, the flow in each first arc-shaped area and the flow in each second arc-shaped area to obtain the flow in the oval cross section;
step 65, replacing the specific flow density with the Mach number of each pitot pressure measuring point, and repeating the step 61 to the step 64 to obtain the superposed Mach number in the elliptical section;
and 66, replacing the specific flow density with the total pressure recovery coefficient of each pitot pressure measuring point, and repeating the steps 61-64 to obtain the total pressure recovery coefficient superposed in the elliptical section.
Specifically, for ease of calculation, standard equations are aimed at
Can be represented by the following parametric equation (12)
In the formula, t is a radial parameter, t is more than or equal to 0 and less than or equal to 1, alpha is a circumferential angle parameter, and alpha is more than or equal to 0 and less than or equal to 2 pi.
In the parametric equation, measure point D ij Has a coordinate of (t) i ,α j ) In the invention, the measuring points are distributed at equal intervals, the coordinates of the measuring points can be represented by intervals delta t and delta alpha, and t i ＝iΔt，Δt＝1/6，α j = j Δ α, Δ α = pi/4, the elliptical area S can be calculated according to equation (6),
S＝∫∫abtdtdα (13)
however, specific flow density m A ij Are discrete points, in order to achieve a specific flow density mA ij Integration in the elliptical cross section requires interpolation using discrete points to obtain a continuous interpolation function, and then integration is performed on the interpolation function. The interpolation function is divided into three parts according to the discrete measurement point distribution positions in the ellipse.
First part, ellipseA sector area consisting of the circle center and the adjacent measuring points at the innermost layer (t is more than or equal to 0 and less than or equal to t 1 ) Within each sector, to the ellipse centre D 00 And a first layer measuring point D 1j 、D 1,j+1 The data is subjected to linear interpolation to obtain an interpolation function m of the specific flow density A (t,α)
The flow q in each sector can be obtained by integrating the interpolation function of the specific flow density m 0j
Second, an annular region (t) surrounded by the innermost measuring point and the outermost measuring point 1 ≤t≤t 5 ) In a first circular arc area consisting of 4 adjacent measuring points, pair D ij 、D i,j+1 、D i+1,j 、D i+1,j+1 Linear interpolation is carried out on the data to obtain an interpolation function m A (t,α)
The flow q in each sector can be obtained by integrating the interpolation function of the specific flow density m ij
Third part, part t between the outermost measuring point and the elliptical shell 5 ≤t≤t 6 I.e. a second circular arc-shaped zone, which is only the zone inner side D 5j If the data of the measuring points are interpolated outwards, larger deviation can occur, so that the measuring points at the outermost layer are adopted to approximate along the radial direction, linear interpolation is adopted along the circumferential direction, and the data are interpolated according to the D 5j 、D 5,j+1 、E j 、E j+1 Fractional flow density interpolation function m in composed region A (t,α)
The flow q in each sector can be obtained by integrating the interpolation function of the specific flow density m 5j
Finally, the flow rates of all parts are superposed, and the flow rate q in the elliptic section is calculated m
According to the periodicity of the ellipse along the circumferential direction, when j =8, m A i,j+1 ＝m A i1 。
The specific flow density in the interpolation formulas (1), (3) and (5) is replaced by the Mach number or the total pressure recovery coefficient of each pitot pressure measuring point, and the distribution of the Mach number or the total pressure recovery coefficient in the oval section can be obtained; and (3) replacing the specific flow density in the flow formula (7) with a Mach number and a total pressure recovery coefficient to obtain a superposed Mach number and a superposed total pressure recovery coefficient.
And 7: and respectively dividing the Mach number superposed in the elliptical section and the total pressure recovery coefficient superposed in the elliptical section by the area of the elliptical section to obtain the average Mach number and the total pressure recovery coefficient in the elliptical section.
Specifically, the average mach number and the total pressure recovery coefficient in the cross section can be obtained by dividing the superposed mach number and the total pressure recovery coefficient by the elliptical area pi ab.
In the examples, the description of the relevant symbols in the formula is shown in table 1:
TABLE 1
The foregoing is merely an example of the present invention and common general knowledge in the art of designing and/or characterizing particular aspects and/or features is not described in any greater detail herein. It should be noted that, for those skilled in the art, without departing from the technical solution of the present invention, several variations and modifications can be made, which should also be regarded as the protection scope of the present invention, and these will not affect the effect of the implementation of the present invention and the practicability of the patent. The scope of the claims of the present application shall be determined by the contents of the claims, and the description of the embodiments and the like in the specification shall be used to explain the contents of the claims.
Claims (3)
1. The utility model provides a 41 measurement stations equidistance distribution's oval cross-section flowmeter which characterized in that: the device comprises an elliptical shell, the elliptical shell is connected with an elliptical cross-section outlet at the downstream of an air inlet, the elliptical shell is used for mounting a leather-supporting rake and a static pressure measuring device, 8 rake positions are arranged on the leather-supporting rake along the circumferential direction, the cross-section area of the elliptical outlet is divided into 8 equal parts by the 8 rake positions, 1 leather-supporting pressure measuring point is arranged at the central intersection point of the 8 rake positions, 5 leather-supporting pressure measuring points are uniformly arranged on each rake position along the radial direction, a leather-supporting pressure sensor is mounted at each leather-supporting pressure measuring point, the static pressure measuring device comprises 8 wall surface pressure sensors uniformly arranged on the elliptical shell along the circumferential direction, the position distribution of the wall surface pressure sensors corresponds to the rake positions, the leather-supporting pressure sensors are used for acquiring the leather-supporting pressure of the corresponding leather-supporting pressure measuring points, and the wall surface pressure sensors are used for measuring the static pressure of a flow field at the corresponding positions.
2. A data processing method of an elliptic cross section flowmeter with equidistantly distributed 41 measuring points is applied to the elliptic cross section flowmeter with equidistantly distributed 41 measuring points, which is characterized by comprising the following steps:
step 1: averaging the measurement results of the 8 wall surface pressure sensors to obtain the static pressure of the center of the elliptical section, and obtaining the static pressure of each pitot pressure measurement point by utilizing a linear interpolation method along the radial direction according to the static pressure of the center of the elliptical section and the wall surface static pressure corresponding to the corresponding rake position;
and 2, step: calculating the Mach number of the airflow of each pitot pressure measuring point according to the pitot pressure and the static pressure of each pitot pressure measuring point;
and step 3: calculating to obtain the total airflow pressure of each pitot voltage measuring point according to the static pressure and the airflow Mach number of each pitot voltage measuring point;
and 4, step 4: calculating the specific flow density of each pitot pressure measuring point according to the airflow Mach number and the airflow total pressure of each pitot pressure measuring point;
and 5: calculating to obtain a total pressure recovery coefficient of each pitot pressure measuring point according to the total airflow pressure and the total incoming flow pressure of each pitot pressure measuring point;
step 6: respectively integrating and superposing the specific flow density, the airflow Mach number and the total pressure recovery coefficient of each pitot pressure measuring point in the elliptical section to obtain the flow, the superposed Mach number and the superposed total pressure recovery coefficient in the elliptical section;
and 7: and respectively dividing the superposed Mach number and the superposed total pressure recovery coefficient in the elliptical section by the area of the elliptical section to obtain the average Mach number and the total pressure recovery coefficient in the elliptical section.
3. The data processing method of the elliptic cross section flowmeter with the 41 measuring points equidistantly distributed as claimed in claim 2, wherein the step of integrating and superposing the specific flow density, the airflow mach number and the total pressure recovery coefficient of each pitot pressure measuring point in the elliptic cross section to obtain the flow, the superposed mach number and the superposed total pressure recovery coefficient in the elliptic cross section comprises the following steps:
step 61, interpolating the specific flow density of each pitot pressure measuring point by utilizing a linear interpolation method in a sector area consisting of the center of the elliptical section and two adjacent pitot pressure measuring points on the innermost layer to obtain a specific flow density interpolation function of each sector area, and integrating the interpolation function of the compared flow density in the sector area to obtain the flow in each sector area; the calculation formula of the interpolation function of the specific flow density in the sector area is as follows:
the calculation formula of the flow in the sector area is as follows:
step 62, interpolating the specific flow density of each pitot pressure measurement point by using a linear interpolation method in a first arc-shaped area consisting of two pitot pressure measurement points adjacent to the innermost layer and two pitot pressure measurement points adjacent to the outermost layer corresponding to the two adjacent pitot pressure measurement points to obtain a specific flow density interpolation function of each first arc-shaped area, and integrating the interpolation function of the specific flow density in the first arc-shaped area to obtain the flow in each first arc-shaped area, wherein the calculation formula of the specific flow density interpolation function in the first arc-shaped area is as follows:
the calculation formula of the flow in the first circular arc-shaped area is as follows:
step 63, approximating a second arc-shaped area composed of two adjacent pitot pressure measuring points on the outermost layer and a corresponding elliptical shell by adopting the measuring points on the outermost layer along the radial direction, interpolating the specific flow density of each pitot pressure measuring point by adopting linear interpolation values along the circumferential direction to obtain a specific flow density interpolation function of each second arc-shaped area, and integrating the interpolation function of the specific flow density in the second arc-shaped area to obtain the flow in each second arc-shaped area; wherein the calculation formula of the interpolation function of the specific flow density in the second circular arc-shaped area is as follows:
the calculation formula of the flow in the second arc-shaped area is as follows:
step 64, superposing the flow in each sector area, the flow in each first arc-shaped area and the flow in each second arc-shaped area to obtain the flow in the oval cross section; wherein, the calculation formula of the flow in the elliptical section is as follows:
step 65, replacing the specific flow density in the formulas (1), (3) and (5) with the Mach number of each pitot pressure measuring point, and repeating the steps 61-64 to obtain the superposed Mach number in the elliptical section;
step 66, replacing the specific flow density in the formulas (1), (3) and (5) with the total pressure recovery coefficient of each pitot pressure measuring point, and repeating the steps 61 to 64 to obtain the total pressure recovery coefficient superposed in the elliptical section;
wherein a is a semi-major axis of the elliptical section; b is a semi-short axis with an elliptical cross section; i is the position number of a measuring point on each rake position, i =1,2, \8230, and \82305; j is the rake position number, j =1,2, \8230, 8; m is A 00 Is the specific flow density at the center of the elliptical cross section, m A ij For each measuring point D ij Specific flow density; m is A (t, α) is a specific flow density interpolation function; t is a radial parameter, t is more than or equal to 0 and less than or equal to 1, and delta t is a radial parameter interval; alpha is a circumferential angle parameter, alpha is more than or equal to 0 and less than or equal to 2 pi, and deltaAlpha is the circumferential angle parameter interval; q. q.s m Is the flow rate in the elliptical section; q. q.s m 0j Flow rate q in a sector area formed by two pitot pressure measuring points adjacent to the center and the innermost layer of the elliptical cross section m ij Traffic in different cells.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110314172.6A CN113125103B (en) | 2021-03-24 | 2021-03-24 | Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110314172.6A CN113125103B (en) | 2021-03-24 | 2021-03-24 | Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113125103A CN113125103A (en) | 2021-07-16 |
| CN113125103B true CN113125103B (en) | 2023-02-17 |
Family
ID=76774178
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110314172.6A Active CN113125103B (en) | 2021-03-24 | 2021-03-24 | Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113125103B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113670565B (en) * | 2021-08-12 | 2022-06-07 | 同济大学 | Wind field measuring device and measuring method for wind power generation high tower model test |
| CN114061896B (en) * | 2021-11-17 | 2023-04-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method for processing internal resistance measurement test data of lifting body ventilation model |
| CN114528778B (en) * | 2022-01-24 | 2022-11-18 | 中国航空工业集团公司北京长城计量测试技术研究所 | Air flow measuring method based on bell mouth flow tube speed distribution model |
| CN114200091B (en) * | 2022-02-18 | 2022-05-27 | 中国矿业大学(北京) | Method for measuring underground roadway harmful substance concentration spatial distribution data |
| CN118376381B (en) * | 2024-06-24 | 2024-08-30 | 中国空气动力研究与发展中心低速空气动力研究所 | Pressure testing device and method for pressure testing rake and air inlet channel |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010285871A (en) * | 2009-06-09 | 2010-12-24 | Tigers Polymer Corp | Ventilation duct |
| CN103063399A (en) * | 2013-01-09 | 2013-04-24 | 浙江大学 | Framed bent device for wind tunnel flow field |
| CN106482923A (en) * | 2016-10-28 | 2017-03-08 | 北京航天长征飞行器研究所 | Test device is demarcated in a kind of flow field being applied under hot environment |
| CN107560822A (en) * | 2017-10-31 | 2018-01-09 | 中国商用飞机有限责任公司 | Anti-icing wind tunnel test equipment for front edge of engine air inlet and test method thereof |
| CN109098995A (en) * | 2017-06-21 | 2018-12-28 | 中车青岛四方机车车辆股份有限公司 | A kind of traction electric machine cooling system air duct and draught fan performance test method and device |
| CN208333564U (en) * | 2018-06-08 | 2019-01-04 | 华北电力大学 | A kind of axial flow blower flow measurement device based on cross probe |
| CN111426480A (en) * | 2020-04-17 | 2020-07-17 | 西安翔迅科技有限责任公司 | Aeroengine exhaust measuring device |
| CN111551339A (en) * | 2020-04-28 | 2020-08-18 | 中国空气动力研究与发展中心高速空气动力研究所 | Design method of special pressure measuring rake for deformed flow channel of wing-body fusion aircraft ventilation model |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE602004024194D1 (en) * | 2003-05-19 | 2009-12-31 | Michael Trachtenberg | METHOD AND DEVICE FOR GAS SEPARATION |
| FR2860293B1 (en) * | 2003-09-26 | 2006-01-20 | Airbus France | DEVICE AND METHOD FOR MEASURING PARAMETERS IN FLIGHT TESTING OF AN AIRCRAFT |
| US7401453B2 (en) * | 2004-09-02 | 2008-07-22 | Graphic Packaging International, Inc. | Packaging system having loading carousel |
| US7810391B2 (en) * | 2008-04-04 | 2010-10-12 | Honda Motor Co., Ltd. | Aerodynamic hood lift and deflection measurement device |
| CN202101706U (en) * | 2011-04-22 | 2012-01-04 | 燕山大学 | Electromagnetic flowmeter sensor with elliptical inner pipeline |
| US9027392B2 (en) * | 2013-05-07 | 2015-05-12 | Honeywell International Inc. | Multi-hole probes and methods for manufacturing multi-hole probes |
| CN204740120U (en) * | 2015-06-12 | 2015-11-04 | 中国能源建设集团广东省电力设计研究院有限公司 | A make wind rectification system for floating fan model test |
| CN106500949B (en) * | 2016-10-31 | 2019-01-29 | 武汉一冶钢结构有限责任公司 | A kind of wind-tunnel round corner section elliptical ring shell forming method |
-
2021
- 2021-03-24 CN CN202110314172.6A patent/CN113125103B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010285871A (en) * | 2009-06-09 | 2010-12-24 | Tigers Polymer Corp | Ventilation duct |
| CN103063399A (en) * | 2013-01-09 | 2013-04-24 | 浙江大学 | Framed bent device for wind tunnel flow field |
| CN106482923A (en) * | 2016-10-28 | 2017-03-08 | 北京航天长征飞行器研究所 | Test device is demarcated in a kind of flow field being applied under hot environment |
| CN109098995A (en) * | 2017-06-21 | 2018-12-28 | 中车青岛四方机车车辆股份有限公司 | A kind of traction electric machine cooling system air duct and draught fan performance test method and device |
| CN107560822A (en) * | 2017-10-31 | 2018-01-09 | 中国商用飞机有限责任公司 | Anti-icing wind tunnel test equipment for front edge of engine air inlet and test method thereof |
| CN208333564U (en) * | 2018-06-08 | 2019-01-04 | 华北电力大学 | A kind of axial flow blower flow measurement device based on cross probe |
| CN111426480A (en) * | 2020-04-17 | 2020-07-17 | 西安翔迅科技有限责任公司 | Aeroengine exhaust measuring device |
| CN111551339A (en) * | 2020-04-28 | 2020-08-18 | 中国空气动力研究与发展中心高速空气动力研究所 | Design method of special pressure measuring rake for deformed flow channel of wing-body fusion aircraft ventilation model |
Non-Patent Citations (1)
| Title |
|---|
| 矩形截面高超声速进气道高焓风洞实验;袁化成 等;《航空动力学报》;20101231;2748-2755 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113125103A (en) | 2021-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113125103B (en) | Data processing method for 41 measuring point equidistantly distributed oval cross-section flow meter | |
| CN113405763B (en) | Measuring point arrangement method of '3 +4+ 5' rice-shaped harrow flow meter | |
| Burton et al. | A literature review of low pressure steam turbine exhaust hood and diffuser studies | |
| Zhang et al. | Criteria for the matching of inlet and outlet distortions in centrifugal compressors | |
| Gaetani et al. | Investigation of the flow field in a high-pressure turbine stage for two Stator-Rotor axial gaps—part II: unsteady flow field | |
| CN106545411B (en) | Simulate the design method of the direct-connected experimental rig of Distorted Flow Field | |
| CN110990955A (en) | Hypersonic speed Bump air inlet channel design method and system | |
| CN111780949B (en) | CFD analysis-based total pressure correction method for high-speed air inlet channel precursor wind tunnel experiment | |
| Levine | Two-dimensional inflow conditions for a supersonic compressor with curved blades | |
| CN106919730B (en) | Wind power plant wake flow determination method adopting wind speed attenuation factor | |
| Xue et al. | Experimental and numerical analysis of different unsteady modes in a centrifugal compressor with variable vaned diffuser | |
| Oh et al. | A numerical investigation of the dual performance characteristics of a small propeller fan using viscous flow calculations | |
| Wolter et al. | Development of a rotating rake array for boundary-layer-ingesting fan-stage measurements | |
| Asghar et al. | Investigation of a passive flow control device in an s-duct inlet of a propulsion system with high subsonic flow | |
| CN113361028B (en) | Two-dimensional design method of volute | |
| An et al. | Effects of inclination angle, orientation angle, and hole length on film cooling effectiveness of rectangular diffusion holes | |
| Luan et al. | Numerical study on aerodynamic noise performances of axial spacing in a contra-rotating axial fan | |
| CN113503262B (en) | Simulation method for high-low pressure turbine of two-stage supercharging system | |
| CN105738067B (en) | The fast determination method of parameter after a kind of twice homonymy oblique shock wave is intersecting | |
| Xue et al. | Experimental and numerical analysis of unsteady flow structure in a centrifugal compressor with variable vaned diffuser | |
| Atassi et al. | Implementation of nonreflecting boundary conditions for the nonlinear Euler equations | |
| Yang et al. | Surrogate-based optimization of nacelle intake with fan-intake interaction: Mitigating flow separation under crosswind | |
| Zhou et al. | Multi-Objective optimization of an IGV for a large axial fan based on NSGA-II | |
| Saha et al. | Influence of prehistory and leading edge contouring on aero performance of a three-dimensional nozzle guide vane | |
| CN106780765A (en) | Stream flow field 3-D quantitative reconstructing method in a kind of supersonic speed |
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 |