CN108194412A - The uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 and method for designing impeller - Google Patents

Abstract

The invention discloses a kind of uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 and method for designing impeller, main purpose is to improve the efficiency of CO2 compressor group, the power consumption of CO2 compressor group is reduced, ensures that compressor has wider condition range, while can improve the stability of rotor.The Stage includes：Impeller, vane diffuser, bend and return channel, wherein, the impeller is the entry position for being located at Stage with reference to low frictional resistance wedge shape impeller, the impeller of rule design according to particular vane thickness；The vane diffuser is equipped in the outlet of the impeller, the return channel is located at the outlet port of Stage, is connected between the vane diffuser and the return channel by the bend；The machine Mach number Mu of the Stage₂=0.6~0.9, discharge coefficient Φ₁=0.0056, design point energy head coefficient τ=0.61, polytropic efficiency η under each Mach number_pol=0.64~0.645, the range of flow of application for design point 75% 146%.

Description

The uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 and method for designing impeller

Technical field

The present invention relates to Design of Stage field, more particularly to a kind of uniaxial CO2 compressor end of discharge coefficient 0.0056 Segment model grade and method for designing impeller.

Background technology

The exploitation of modern centrifugal compressor is mostly according to user demand, and based on existing Stage database, utilization is similar Theory carries out conceptual design.Therefore, the quality of Design of Stage directly affects the performance of entire product.It is filled as urea synthesizing The important component put, CO2 compressor group have become major compression because having the characteristics that power pressure is high, density is big Most there are technology content, the most one of product of emphasis in machine producer.

At present, CO2 compressor what discharge coefficient behind high pressure cylinder is all very small, some discharge coefficient Φ₁Even It can be to 0.004 or lower.However, the polytropic efficiency of this small discharge coefficient Stage is generally lower than 60%, it far can not Meets the needs of market.Therefore, the meaning for researching and developing the small discharge coefficient Stage of efficient CO2 units latter end is very great, is The technical issues of design CO2 compressor field is urgently to be resolved hurrily at present.

Invention content

In view of this, the present invention provides a kind of uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 and impeller is set Meter method, main purpose are to improve the efficiency of CO2 compressor group, reduce the power consumption of CO2 compressor sets, ensure compression equipment There is wider condition range, while the stability of rotor can be improved.

According to the present invention on one side, a kind of uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056, packet are provided It includes：

The Stage is located in compressor, including impeller (1), vane diffuser (2), bend (3) and return channel (4), Wherein, the impeller (1) is according to low frictional resistance wedge shape impeller of the particular vane thickness distribution with reference to rule design, the impeller (1) positioned at the entry position of Stage；The vane diffuser (2), the return channel are equipped in the outlet of the impeller (1) (4) pass through the bend (3) phase between the outlet port of Stage, the vane diffuser (2) and the return channel (4) Connection；The machine Mach number Mu of the Stage₂=0.6~0.9, discharge coefficient Φ₁=0.0056, design point energy head coefficient τ= 0.61, polytropic efficiency η under each Mach number_pol=0.64~0.645, the range of flow of application is the 75%-146% of design point.

According to the present invention on the other hand, a kind of uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 is provided The design method of middle impeller：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, ensures that the Yplus values of wall surface first layer grid are small In 1, while the length-width ratio of grid is less than 3000, extends than being less than 3, the minimal orthogonality of grid is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

The uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 provided by the invention is single shaft CO2 compressor latter ends Special purpose model grade, efficient, design point energy head coefficient height, hub ratio are big, and vane thickness and height are larger, using the present invention Stage can cause CO2 compressor to have higher operational efficiency and wider condition range, reduce the work(of CO2 compressor group Consumption, while the stability of rotor can be improved.In addition, the uniaxial CO2 compressor latter end of discharge coefficient 0.0056 provided by the invention The design method of impeller in Stage can achieve the purpose that reduce wetted perimeter, effectively reduce the friction damage between fluid and wall surface It loses, so as to design low frictional resistance impeller.

Above description is only the general introduction of technical solution of the present invention, in order to better understand the technological means of the present invention, And it can be implemented in accordance with the contents of the specification, and in order to allow above and other objects of the present invention, feature and advantage can It is clearer and more comprehensible, below the special specific embodiment for lifting the present invention.

Description of the drawings

By reading the detailed description of hereafter preferred embodiment, it is various other the advantages of and benefit it is common for this field Technical staff will become clear.Attached drawing is only used for showing the purpose of preferred embodiment, and is not considered as to the present invention Limitation.And throughout the drawings, the same reference numbers will be used to refer to the same parts.In the accompanying drawings：

Fig. 1 is the schematic diagram of CO2 compressor Stage that discharge coefficient provided in an embodiment of the present invention is 0.0056；

Fig. 2 is the CO2 compressor Stage impeller disengaging bicker that discharge coefficient provided in an embodiment of the present invention is 0.0056 Schematic diagram；

Fig. 3 be discharge coefficient provided in an embodiment of the present invention be 0.0056 CO2 compressor Stage vane diffuser into Angle of outlet schematic diagram；

Fig. 4 is the CO2 compressor Stage return channel inlet and outlet that discharge coefficient provided in an embodiment of the present invention is 0.0056 Angle schematic diagram；

Fig. 5 flows to impeller blade load distribution schematic diagram for edge provided in an embodiment of the present invention；

Fig. 6 is the high section relative velocity distribution map of 10% leaf of impeller provided in an embodiment of the present invention；

Fig. 7 is the high section relative velocity distribution map of 50% leaf of impeller provided in an embodiment of the present invention；

Fig. 8 is the high section relative velocity distribution map of 90% leaf of impeller provided in an embodiment of the present invention；

Fig. 9 is meridional channel relative velocity distribution map provided in an embodiment of the present invention；

Figure 10 is the impeller outlet air-flow angular distribution provided in an embodiment of the present invention extended to relative position；

Figure 11 is the impeller outlet static pressure distribution figure provided in an embodiment of the present invention extended to relative position；

Figure 12 flows to distribution schematic diagram for wedge-shaped impeller blade angle provided in an embodiment of the present invention with vane thickness edge；

Figure 13 is wedge-shaped impeller channel profile provided in an embodiment of the present invention and conventional impellers runner silhouette contrast schematic diagram；

Figure 14 is that provided in an embodiment of the present invention extend to the impeller outlet flow angle of relative position is distributed basis for estimation figure；

Figure 15 is the impeller outlet static pressure distribution basis for estimation figure provided in an embodiment of the present invention extended to relative position；

Figure 16 be edge provided in an embodiment of the present invention flow to impeller blade suction surface and pressure face load distribution judge according to According to figure；

Figure 17 is the high section relative velocity distribution map of 10% leaf of vane diffuser provided in an embodiment of the present invention；

Figure 18 is the high section relative velocity distribution map of 50% leaf of vane diffuser provided in an embodiment of the present invention；

Figure 19 is the high section relative velocity distribution map of 90% leaf of vane diffuser provided in an embodiment of the present invention；

Figure 20 is the high section relative velocity distribution map of 10% leaf of return channel provided in an embodiment of the present invention；

Figure 21 is the high section relative velocity distribution map of 50% leaf of return channel provided in an embodiment of the present invention；

Figure 22 is the high section relative velocity distribution map of 90% leaf of return channel provided in an embodiment of the present invention；

Figure 23 is pressure ratio performance curve schematic diagram under different Mach number provided in an embodiment of the present invention；

Figure 24 is polytropic efficiency performance curve schematic diagram under different Mach number provided in an embodiment of the present invention；

Figure 25 is can head coefficient performance curve schematic diagram under different Mach number provided in an embodiment of the present invention.

Specific embodiment

The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in attached drawing Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here It is limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure Completely it is communicated to those skilled in the art.

To solve the above-mentioned problems, an embodiment of the present invention provides a kind of uniaxial CO2 compressor ends of discharge coefficient 0.0056 Segment model grade, suitable for the Modulated Design of uniaxial CO2 compressor product, as shown in Fig. 1, the Stage is located in compressor, Including impeller 1, vane diffuser 2, bend 3 and return channel 4, wherein, the impeller 1 is with reference to rule according to particular vane thickness Low frictional resistance wedge shape impeller, the impeller 1 of design are located at the entry position of Stage；It is equipped in the outlet of the impeller 1 described Vane diffuser 2, the return channel 4 is located at the outlet port of Stage, between the vane diffuser 2 and the return channel 4 It is connected by the bend 3；As in Figure 2-4, respectively illustrate impeller 1, vane diffuser 2, return channel 4 blade into Angle of outlet schematic diagram, the machine Mach number Mu of the Stage₂=0.6~0.9, discharge coefficient Φ₁=0.0056, design point energy Head coefficient τ=0.61, polytropic efficiency η under each Mach number_pol=0.64~0.645, the range of flow of application is design point 75%-146%, specific performance curve is as shown in Figure 23~25.

The hub ratio for the CO2 compressor Stage that discharge coefficient in the embodiment of the present invention is 0.0056 is very big, than general Stage hub ratio ds/D₂=0.34 improves about 33%, this can greatly improve the stabilization of the rotor using the Stage Property.

Specifically, the impeller 1 is the two-dimensional impeller of enclosed, and the basic parameter of two-dimensional impeller is as follows：Impeller outlet diameter D₂=450m, the m number of blade Z₁=11, impeller opposite outlet width b₂/D₂=0.0214, b₂Impeller outlet for the impeller (1) Width；The two-dimensional impeller is close to wheel cap and wheel disc lateral lobe piece stagger angle β_1AsAnd β_1AhIt it is 24 °, the two-dimensional impeller leans on Nearly wheel cap and the exit vane angle beta of wheel disc side_2AsAnd β_2AhIt it is 22 °, impeller blade import thickness is 4mm, and impeller blade exports Thickness is 14mm, and thickness is first increases and then decreases from the inlet to the outlet, maximum gauge 18mm, much larger than general vane thickness, I.e. vane thickness is to increase to 18mm from the 4mm of entrance, then be reduced to 14mm.

The wheel cap side of the impeller 1 and the meridional channel of reel side are by one section of circular arc and two straightways tangent with it It forms, the pressure face of the impeller (1) and the blade loading maximum position of suction surface are located at 30% length of blade.Such as Fig. 5 institutes Show, this load is distributed so that the Stage can be reduced effectively in runner with very high efficiency and wider condition range Friction loss.Under each Mach number design conditions, the flowing of the impeller has following features：

1st, impeller is on 10%, 50%, the 90% high section of leaf, the no any flow separation of flowing among impeller channel, As shown in Figure 6 to 8.

2nd, unsteady three-dimensional viscous flows CFD analysis results show on the meridional channel of the impeller also without any stream Dynamic separation, as shown in Figure 9.

3rd, 1 outlet flow angle of impeller is evenly distributed, and average outlet flow angle is 8~12 °, and fluctuation range is less than 5 degree, such as schemes Shown in 10.

4th, 1 exit static pressure distributing homogeneity of impeller is good, and fluctuation range is within the 2% of outlet average static pressure, such as Figure 11 institutes Show.

5th, the polytropic efficiency of impeller is very high, and the polytropic efficiency at impeller outlet is up to 89.5%.

Specifically, the blade of the vane diffuser 2 uses overall height short blade, it is therefore an objective to while efficiency is improved, also Take into account wider condition range.The vane diffuser 2 is close to wheel cap and wheel disc lateral lobe piece stagger angle β_3AsAnd β_3Ah It it is 8 °, the vane diffuser 2 is close to the exit vane angle beta of wheel cap and wheel disc side_4AsAnd β_4AhIt is 19 °, diffusion degree is 5.18, Flow losses at this time in vane diffuser are smaller.

The meridian molded line of the reel side of the vane diffuser 2 is made of one section of straight line, and the meridian molded line of shrouding disc side is It is made of one section of circular arc and other three sections of straight lines；Wherein, the meridian molded line of reel side is perpendicular to axial direction.

The entrance width of the vane diffuser (2) and the ratio b of the impeller outlet width of the impeller (1)₃/b₂= 0.793, the entrance relative position D of the vane diffuser (2)₃/D₂=1.067, the outlet of the vane diffuser (2) is opposite Position D₄/D₂=1.45, D₂For impeller outlet diameter, D₃For the vane diffuser (2) inlet diameter, D₄For the blade diffusion The outlet diameter of device (2).

Above-mentioned design ensure that 1 exit flow of impeller can smoothly enter vane diffuser 2, and will not expand because of runner It opens and generates apparent reflux, be particularly advantageous to the flowing of low flow rate condition；2 exit width of vane diffuser and blade expand simultaneously The ratio b of 2 entrance width of depressor₄/b₃≈0.719.The design can be very good the flow field at matching impeller outlet, and blade is inhibited to expand The flow separation of depressor shrouding disc side reduces flow losses of the vane diffuser without leaf section.As shown in Figure 17~Figure 19, which expands Almost without flow separation on depressor 2.

Optionally, the inlet and outlet width of the bend 3 compares b₅/b₄=1.298, the b₅Go out mouth width for the bend 3 Degree, the b₄For the exit width of the vane diffuser 2, the unsteady three-dimensional viscous flows CFD analysis knots of the Stage Fruit shows that at each Mach number design point 3 internal flow of bend is good, is not clearly separated, as shown in Figure 9.

Optionally, the blade of the return channel 4 be overall height banana airfoil fan, number of blade Z₂=18, the entrance phase of blade To position D₅/D₂=1.431, entrance established angle is 27 °, blade exit relative position D₆/D₂=0.751, exit installation angle is 105 °, the D₅For the blade inlet diameter of the return channel 4, the D₆Blade exit diameter for the return channel 4.Reflux The definition of device import and export established angle is as shown in Figure 4；The design can be very good the flow field of matching bend outlet, and cause reflux Device exit flow levels off to axial direction, reduces the adverse effect between grade.Meanwhile return channel blade is under each Mach number design conditions There is no apparent flow separation, as shown in Figure 20~Figure 22.The wheel cap side meridian molded line of the return channel 4 is an angled straight lines Section, the entrance of the return channel 4 and the wheel cap side circular arc of the bend 3 are tangent；Reel side meridian molded line is a vertical straight line Section, entrance and the 3 reel side circular arc of bend are tangent；The outlet section shrouding disc side and reel side of the return channel 4 are by one section Circular arc and the straightway tangent with it composition.

The uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 provided in an embodiment of the present invention is compressed for single shaft CO2 Machine latter end special purpose model grade, efficient, design point energy head coefficient height, hub ratio are big, and vane thickness and height are larger, using this The Stage of invention can cause CO2 compressor to have higher operational efficiency and wider condition range, while can improve The stability of rotor.

Further, the embodiment of the present invention also provides a kind of core for the CO2 compressor Stage that discharge coefficient is 0.0056 The design method of center portion part impeller, the method includes：

201st, according to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds.

For example, default impeller diameter D₂It can be 450mm, preset hub ratio ds/D₂Can be 0.453.

202nd, according to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller outlet established angle β_2A, impeller outlet width b₂, impeller inlet width D₀。

For the embodiment of the present invention, the step 202 can specifically include：Utilize impeller eye absolute velocity c₁, impeller Inlet circumference speed u₁And β_1A=tan^-1c₁/u₁+ i determines vane inlet established angle β_1A；

Utilize impeller outlet radial direction discharge coefficientPeripheral speed coefficientWithDetermine that impeller goes out Mouth established angle β_2A；

Utilize impeller eye mass flow G, impeller outlet blockage factor τ₂, impeller outlet specific volume compare k_v2, impeller outlet circle Circular velocity u₂, the D₂, impeller eye fluid density ρ_i, it is describedWith default impeller outlet width calculation formula, impeller is calculated Exit width b₂；

Utilize the ds, the D₂, the τ₂, it is describedThe k_v2, velocity coeffficient k_c, impeller inlet specific volume compare k_v0、 Diameter compares K_DWith default impeller inlet width calculation formula, impeller inlet width D is calculated₀；

Wherein, the τ₂Calculation formula be：

The default impeller outlet width calculation formula is：

The default impeller inlet width calculation formula is：

For example, default impeller blade import thickness d₁It can be 4mm, preset impeller blade exit thickness d₂Can be 14mm. Conventional impellers vane thickness is generally constant 4mm or so, and vane thickness is greatly improved in the present invention, and gives suitably The regularity of distribution.

203rd, according to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridian of impeller is determined The preliminary moulding of runner.

Wherein, the predefined impeller cover side meridian molded line can be customized before step 201, specifically, The impeller cover side meridian molded line can be defined as an angled straight lines.Design can be reduced by being defined as angled straight lines Degree of freedom facilitates subsequent modification.

204th, according to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine the blade type of impeller The preliminary moulding of line.

For the embodiment of the present invention, the step 204 can specifically include：According to blade angle along meridian flow direction in multinomial Formula change profile, vane thickness are along meridian flow direction in multinomial change profile, the β_1A, the β_2A, the d₂With default impeller Vane inlet thickness d₁, determine the preliminary moulding of vane type line of impeller, wherein, the maximum value position of the blade angle is 40% At length of blade, the maximum value position of vane thickness is at 80% length of blade.Such as Figure 12, determine that vane type line is preliminary with this Moulding.Figure 13 is wedge-shaped impeller channel profile and conventional impellers runner silhouette contrast schematic diagram provided in an embodiment of the present invention, can To see the runner of wedge-shaped impeller due to improving vane thickness and blade height so that flow channel cross-section length-width ratio a/b more connects 1 is bordering on, reduces wetted perimeter, so as to effectively reduce the friction loss of impeller internal flowing.

205th, according to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined.

206th, obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length.

207th, three-dimentional structured mesh division is carried out to obtained impeller pattern.

In embodiments of the present invention, ensure that the Yplus values of wall surface first layer grid are less than 1, while net by step 207 The length-width ratio of lattice is less than 3000, extends than being less than 3, the minimal orthogonality of grid is more than 15 °.

208th, Three Dimensional Viscous numerical simulation is carried out to the impeller pattern after mesh generation using CFD analysis softwares, analyzes leaf Take turns efficiency, condition range, blade flow field.

It should be noted that it is that the grid generated is imported into CFD analysis softwares using the process of CFD analysis softwares In, Three Dimensional Viscous numerical simulation is carried out to the impeller；Using Spalart-Allmaras turbulence models, working medium is using preferable empty Gas, inlet boundary condition give total temperature, stagnation pressure, axial admission, and export boundary condition gives mass flow；Then what is given is above-mentioned Parameter analysis calculates impeller adiabatic efficiency, condition range, blade flow field.

Since the variable considered in design is more, degree of freedom is larger, needs to carry out assisting setting using Multipurpose Optimal Method Meter.Three-dimensional parameterized fitting is carried out to obtained scheme, including lid side molded line, vane thickness distribution and blade angle point Cloth, give each parameter solution interval and each Mach number under the desired value of polytropic efficiency and pressure ratio, using CFD analysis softwares Three Dimensional Viscous optimization is carried out, to seek obtaining optimal case.

209th, whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range More than or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement.

Wherein, the default efficiency threshold and default condition range threshold value can carry out setting this hair according to the design needs Bright embodiment does not limit.For example, default efficiency threshold can be 80%, 85% etc., default condition range threshold value can be 70%-85%, 80%-85% etc..

210th, if it is not, then adjusting the matching of meridional channel and the vane type line, the distribution of modification blade beta angles, blade are thick Degree distribution, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal in advance until obtaining the impeller adiabatic efficiency If efficiency threshold, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default leaf Take turns the impeller of flow field design requirement.

It should be noted that if the impeller adiabatic efficiency is less than default efficiency threshold, the condition range is less than default operating mode Range threshold or the blade flow field do not meet default blade flow field design requirement, then perform redesign impeller, until It is more than or equal to default efficiency threshold to the impeller adiabatic efficiency, the condition range is more than or equal to default condition range threshold Value and the blade flow field meet the impeller of default blade flow field design requirement.

Wherein, the impeller is low frictional resistance wedge shape impeller, and the default blade flow field design requirement includes impeller in difference Flowing on the high section of leaf is without being clearly separated；Without flow separation on the meridional channel of impeller；The outlet flow angle distribution of impeller is equal It is even；Exit static pressure is evenly distributed；Load changes in distribution on blade pressure surface and suction surface is uniform.Specifically, meet impeller stream Design requirement is：I impellers 1 are on 10%, 50%, the 90% high section of leaf, the no any stream of flowing among 1 runner of impeller Dynamic separation, as shown in Figure 6 to 8；Also without any flow separation on the meridional channel of II impellers 1, as shown in Figure 9；III The outlet flow angle of impeller 1 is evenly distributed, and outlet flow angle fluctuation range is less than 5 °, as shown in figure 14；IV exit static pressures are distributed Uniformity is good, and fluctuation range is within the 2% of outlet average static pressure, as shown in figure 15；The static pressure of V blade pressure surfaces and suction surface Changes in distribution is uniform, as shown in figure 16.

In embodiments of the present invention, the step 210 is specially：To the impeller adiabatic efficiency and operating mode model of CFD analysis gained Carry out integrated judgment is enclosed, and stream field result is analyzed, when blade flow field does not meet design requirement, then adjust impeller meridian And the matching of vane type line, it is distributed by changing blade beta angles, vane thickness is distributed, meridian molded line and CFD are analyzed repeatedly Iteration, it is final so that the impeller adiabatic efficiency of gained is more than or equal to default efficiency threshold, the condition range is more than or Meet default blade flow field design requirement equal to default condition range threshold value and the blade flow field；At this point, complete design.

The design side of impeller in the uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 provided in an embodiment of the present invention Method can achieve the purpose that reduce wetted perimeter, the friction loss between fluid and wall surface effectively be reduced, so as to design low rub Hinder impeller.

Based on above-mentioned shown method, correspondingly, the embodiment of the present invention additionally provides a kind of computer readable storage medium, On be stored with computer program, which realizes following steps when being executed by processor：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, to ensure the Yplus values of wall surface first layer grid More than or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimum of grid Orthogonality is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

Based on above-mentioned shown method, the discharge coefficient 0.0056 that the embodiment of the present invention is additionally provided described in a kind of design is single Impeller unit in axis CO2 compressor latter end Stage on a memory and can handled including memory, processor and storage The computer program run on device, the processor realize following steps when performing described program：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, to ensure the Yplus values of wall surface first layer grid More than or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimum of grid Orthogonality is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

The embodiment of the present invention also provides following technical solution：

The uniaxial CO2 compressor latter end Stage of A1, a kind of discharge coefficient 0.0056, including：The Stage is located at compression In machine, including impeller (1), vane diffuser (2), bend (3) and return channel (4), wherein, the impeller (1) is according to specific Vane thickness is with reference to the regular low frictional resistance wedge shape impeller designed, the impeller (1) positioned at the entry position of Stage；In the leaf The outlet for taking turns (1) is equipped with the vane diffuser (2), and the return channel (4) is expanded positioned at the outlet port of Stage, the blade It is connected between depressor (2) and the return channel (4) by the bend (3)；The machine Mach number Mu of the Stage₂= 0.6~0.9, discharge coefficient Φ₁=0.0056, design point energy head coefficient τ=0.61, polytropic efficiency η under each Mach number_pol= 0.64~0.645, the range of flow of application is the 75%-146% of design point.

A2, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1, the leaf of the latter end Stage Take turns hub ratio ds/D₂=0.453, ds are hub diameter, D₂Impeller outlet diameter for the impeller (1).

A3, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1, the impeller (1) is enclosed Two-dimensional impeller, the basic parameter of two-dimensional impeller is as follows：

Impeller outlet diameter D₂=450mm, number of blade Z₁=11, impeller opposite outlet width b₂/D₂=0.0214, for institute State the impeller outlet width of impeller (1)；The two-dimensional impeller is close to wheel cap and wheel disc lateral lobe piece stagger angle β_1AsAnd β_1Ah It it is 24 °, the two-dimensional impeller is close to the exit vane angle beta of wheel cap and wheel disc side_2AsAnd β_2AhIt it is 22 °, impeller blade import is thick It spends for 4mm, impeller blade exit thickness is 14mm.

The uniaxial CO2 compressor latter end Stage of A4, such as A1 or A3 any one of them discharge coefficient 0.0056,

The wheel cap side of the impeller (1) and the meridional channel of reel side are by one section of circular arc and two straight lines tangent with it Section is formed, and the pressure face of the impeller (1) and the blade loading maximum position of suction surface are located at 30% length of blade.

A5, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1,

Vane diffuser (2) blade uses overall height short blade, and the vane diffuser (2) is close to wheel cap and wheel disc side Vane inlet established angle β_3AsAnd β_3AhIt it is 8 °, the vane diffuser (2) is close to the exit vane angle of wheel cap and wheel disc side β_4AsAnd β_4AhIt it is 19 °, diffusion degree is 5.18.

The uniaxial CO2 compressor latter end Stage of A6, such as A1 or A5 any one of them discharge coefficient 0.0056,

The entrance width of the vane diffuser (2) and the ratio b of the impeller outlet width of the impeller (1)₃/b₂= 0.793, the exit width of the vane diffuser (2) and the ratio b of entrance width₄/b₃≈ 0.719, the vane diffuser (2) Entrance relative position D₃/D₂=1.067, the outlet relative position D of the vane diffuser (2)₄/D₂=1.45, D₂For impeller Outlet diameter, D₃For the vane diffuser (2) inlet diameter, D₄Outlet diameter for the vane diffuser (2).

The uniaxial CO2 compressor latter end Stage of A7, such as A1 or A5 any one of them discharge coefficient 0.0056,

The meridian molded line of the reel side of the vane diffuser (2) is made of one section of straight line, and the meridian molded line of shrouding disc side It is to be made of one section of circular arc and other three sections of straight lines；Wherein, the meridian molded line of reel side is perpendicular to axial direction.

A8, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1,

The inlet and outlet width of the bend (3) compares b₅/b₄=1.298, the b₅For the exit width of the bend (3), The b₄Exit width for the vane diffuser (2).

A9, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1,

The blade of the return channel (4) be overall height banana airfoil fan, number of blade Z₂=18, the entrance relative position of blade D₅/D₂=1.431, entrance established angle is 27 °, blade exit relative position D₆/D₂=0.751, exit installation angle is 105 °, institute State D₅For the blade inlet diameter of the return channel (4), the D₆Blade exit diameter for the return channel (4).

A10, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in A1,

The wheel cap side meridian molded line of the return channel (4) be an angled straight lines section, the entrance of the return channel (4) with The wheel cap side circular arc of the bend (3) is tangent；Reel side meridian molded line is a vertical straight line section, entrance and the bend (3) Reel side circular arc is tangent；The outlet section shrouding disc side and reel side of the return channel (4) are by one section of circular arc and tangent with it straight Line segment forms.

B11, a kind of design such as uniaxial CO2 compressor latter end Stage of A1-A10 any one of them discharge coefficient 0.0056 In impeller method,

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, to ensure the Yplus values of wall surface first layer grid More than or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimum of grid Orthogonality is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

Impeller method in B12, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in B11, it is described According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller outlet peace Fill angle beta_2A, impeller outlet width b₂, impeller inlet width D₀, specifically include：

Utilize impeller eye absolute velocity c₁, impeller eye peripheral speed u₁And β_1A=tan^-1c₁/u₁+ i, determine blade into Mouth established angle β_1A；

Utilize impeller outlet radial direction discharge coefficientPeripheral speed coefficientWithDetermine that impeller goes out Mouth established angle β_2A；

Utilize impeller eye mass flow G, impeller outlet blockage factor τ₂, impeller outlet specific volume compare k_v2, impeller outlet circle Circular velocity u₂, the D₂, impeller eye fluid density ρ_i, it is describedWith default impeller outlet width calculation formula, impeller is calculated Exit width b₂；

Utilize the ds, the D₂, the τ₂, it is describedThe k_v2, velocity coeffficient k_c, impeller inlet specific volume compare k_v0、 Diameter compares K_DWith default impeller inlet width calculation formula, impeller inlet width D is calculated₀；

Wherein, the τ₂Calculation formula be：

The default impeller outlet width calculation formula is：

The default impeller inlet width calculation formula is：

Impeller method in B13, the uniaxial CO2 compressor latter end Stage of the discharge coefficient 0.0056 as described in B11, it is described According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine the preliminary moulding of vane type line of impeller, have Body includes：

Changed according to blade angle along meridian flow direction in multinomial change profile, vane thickness along meridian flow direction in multinomial Distribution, the β_1A, the β_2A, the d₂With default impeller blade import thickness d₁, determine the preliminary moulding of vane type line of impeller, Wherein, the maximum value position of the blade angle at 40% length of blade, the maximum value position of vane thickness is in 80% blade At length.

Leaf in the uniaxial CO2 compressor latter end Stage of B14, such as B11-B13 any one of them discharge coefficient 0.0056 Wheel method, the predefined impeller cover side meridian molded line are an angled straight lines.

C15, a kind of computer readable storage medium, are stored thereon with computer program, when which is executed by processor Realize following steps：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, to ensure the Yplus values of wall surface first layer grid More than or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimum of grid Orthogonality is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

D16, a kind of design such as uniaxial CO2 compressor latter end Stage of A1-A10 any one of them discharge coefficient 0.0056 In impeller unit, including memory, processor and storage on a memory and the computer program that can run on a processor, The processor realizes following steps when performing described program：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, leaf Take turns exit installation angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, the meridional channel of impeller is determined Preliminary moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, at the beginning of determining the vane type line of impeller Walk moulding；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, to ensure the Yplus values of wall surface first layer grid More than or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimum of grid Orthogonality is more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller effect are carried out to the impeller pattern after mesh generation using CFD analysis softwares Rate, condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, and whether the condition range is more than Or equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, the matching of meridional channel and the vane type line is then adjusted, the distribution of modification blade beta angles, vane thickness point Cloth, meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default effect until obtaining the impeller adiabatic efficiency Rate threshold value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default impeller stream The impeller of field design requirement, the impeller are low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without apparent point From；Without flow separation on the meridional channel of impeller；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade Load changes in distribution on pressure face and suction surface is uniform.

In the above-described embodiments, it all emphasizes particularly on different fields to the description of each embodiment, there is no the portion being described in detail in some embodiment Point, it may refer to the associated description of other embodiment.

It is understood that the correlated characteristic in the above method and device can be referred to mutually.In addition, in above-described embodiment " first ", " second " etc. be for distinguishing each embodiment, and do not represent the quality of each embodiment.

It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description, The specific work process of device and unit can refer to the corresponding process in preceding method embodiment, and details are not described herein.

Algorithm and display be not inherently related to any certain computer, virtual system or miscellaneous equipment provided herein. Various general-purpose systems can also be used together with example based on this.As described above, required by constructing this kind of system Structure be obvious.In addition, the present invention is not also directed to any certain programmed language.It should be understood that it can utilize various Programming language realizes the content of invention described herein, and the description done above to language-specific is to disclose this hair Bright preferred forms.

In the specification provided in this place, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention Example can be put into practice without these specific details.In some instances, well known method, structure is not been shown in detail And technology, so as not to obscure the understanding of this description.

Similarly, it should be understood that in order to simplify the disclosure and help to understand one or more of each inventive aspect, Above in the description of exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention：I.e. required guarantor Shield the present invention claims the more features of feature than being expressly recited in each claim.More precisely, as following Claims reflect as, inventive aspect is all features less than single embodiment disclosed above.Therefore, Thus the claims for following specific embodiment are expressly incorporated in the specific embodiment, wherein each claim is in itself Separate embodiments all as the present invention.

Those skilled in the art, which are appreciated that, to carry out adaptively the module in the equipment in embodiment Change and they are arranged in one or more equipment different from the embodiment.It can be the module or list in embodiment Member or component be combined into a module or unit or component and can be divided into addition multiple submodule or subelement or Sub-component.Other than such feature and/or at least some of process or unit exclude each other, it may be used any Combination is disclosed to all features disclosed in this specification (including adjoint claim, abstract and attached drawing) and so to appoint Where all processes or unit of method or equipment are combined.Unless expressly stated otherwise, this specification is (including adjoint power Profit requirement, abstract and attached drawing) disclosed in each feature can be by providing the alternative features of identical, equivalent or similar purpose come generation It replaces.

In addition, it will be appreciated by those of skill in the art that although some embodiments described herein include other embodiments In included certain features rather than other feature, but the combination of the feature of different embodiments means in of the invention Within the scope of and form different embodiments.For example, in the following claims, embodiment claimed is appointed One of meaning mode can use in any combination.

The all parts embodiment of the present invention can be with hardware realization or to be run on one or more processor Software module realize or realized with combination thereof.It will be understood by those of skill in the art that it can use in practice Microprocessor or digital signal processor (DSP) are realized in centrifugal compressor design device according to embodiments of the present invention The some or all functions of some or all components.The present invention is also implemented as performing method as described herein Some or all equipment or program of device (for example, computer program and computer program product).Such reality The program of the existing present invention can may be stored on the computer-readable medium or can have the form of one or more signal. Such signal can be downloaded from internet website to be obtained either providing or in the form of any other on carrier signal It provides.

It should be noted that the present invention will be described rather than limits the invention, and ability for above-described embodiment Field technique personnel can design alternative embodiment without departing from the scope of the appended claims.In the claims, Any reference mark between bracket should not be configured to limitations on claims.Word "comprising" does not exclude the presence of not Element or step listed in the claims.Word "a" or "an" before element does not exclude the presence of multiple such Element.The present invention can be by means of including the hardware of several different elements and being come by means of properly programmed computer real It is existing.If in the unit claim for listing equipment for drying, several in these devices can be by same hardware branch To embody.The use of word first, second, and third does not indicate that any sequence, can explain these words and run after fame Claim.

Claims (10)

1. a kind of uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056, which is characterized in that including：The model level In in compressor, including impeller (1), vane diffuser (2), bend (3) and return channel (4), wherein, the impeller (1) is root According to particular vane thickness with reference to the regular low frictional resistance wedge shape impeller designed, the impeller (1) positioned at the entry position of Stage； The outlet of the impeller (1) is equipped with the vane diffuser (2), and the return channel (4) is described positioned at the outlet port of Stage It is connected between vane diffuser (2) and the return channel (4) by the bend (3)；The machine Mach number of the Stage Mu₂=0.6~0.9, discharge coefficient Φ₁=0.0056, design point energy head coefficient τ=0.61, polytropic efficiency η under each Mach number_pol =0.64~0.645, the range of flow of application is the 75%-146% of design point.

2. the uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 according to claim 1, which is characterized in that

The hub ratio ds/D of the latter end Stage₂=0.453, ds are hub diameter, D₂Impeller for the impeller (1) Outlet diameter.

3. the uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 according to claim 1, which is characterized in that institute The two-dimensional impeller that impeller (1) is enclosed is stated, the basic parameter of two-dimensional impeller is as follows：

Impeller outlet diameter D₂=450mm, number of blade Z₁=11, impeller opposite outlet width b₂/D₂=0.0214, it is the impeller (1) impeller outlet width；The two-dimensional impeller is close to wheel cap and wheel disc lateral lobe piece stagger angle β_1AsAnd β_1AhIt is 24 °, The two-dimensional impeller is close to the exit vane angle beta of wheel cap and wheel disc side_2AsAnd β_2AhIt it is 22 °, impeller blade import thickness is 4mm, impeller blade exit thickness are 14mm.

It is 4. special according to the uniaxial CO2 compressor latter end Stage of 1 or 3 any one of them discharge coefficient of claim 0.0056 Sign is,

The wheel cap side of the impeller (1) and the meridional channel of reel side are by one section of circular arc and the two straightway structures tangent with it Into the pressure face of the impeller (1) and the blade loading maximum position of suction surface are located at 30% length of blade.

5. in a kind of 0.0056 uniaxial CO2 compressor latter end Stage of design claim 1-4 any one of them discharge coefficient Impeller method, which is characterized in that

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller goes out Mouth established angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, determine that the meridional channel of impeller is preliminary Moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine that the vane type line of impeller is tentatively made Type；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, is more than with the Yplus values for ensureing wall surface first layer grid Or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimal orthogonality of grid Property be more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller adiabatic efficiency, work are carried out to the impeller pattern after mesh generation using CFD analysis softwares Condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, the condition range whether be more than or Equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, then adjusting the matching of meridional channel and the vane type line, modification blade beta angles are distributed, vane thickness is distributed, Meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default efficiency threshold until obtaining the impeller adiabatic efficiency Value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default blade flow field and sets The impeller of requirement is counted, the impeller is low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without being clearly separated；Leaf Without flow separation on the meridional channel of wheel；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade pressure surface It is uniform with the load changes in distribution on suction surface.

6. the impeller method in the uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 according to claim 5, It is characterized in that, the foundation presets impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller outlet established angle β_2A, impeller outlet width b₂, impeller inlet width D₀, specifically include：

Utilize impeller eye absolute velocity c₁, impeller eye peripheral speed u₁And β_1A=tan^-1c₁/u₁+ i determines that vane inlet is pacified Fill angle beta_1A；

Utilize impeller outlet radial direction discharge coefficientPeripheral speed coefficientWithDetermine that impeller outlet is pacified Fill angle beta_2A；

Utilize impeller eye mass flow G, impeller outlet blockage factor τ₂, impeller outlet specific volume compare k_v2, impeller outlet circumference speed Spend u₂, the D₂, impeller eye fluid density ρ_i, it is describedWith default impeller outlet width calculation formula, impeller outlet is calculated Width b₂；

Utilize the ds, the D₂, the τ₂, it is describedThe k_v2, velocity coeffficient k_c, impeller inlet specific volume compare k_v0, diameter ratio K_DWith default impeller inlet width calculation formula, impeller inlet width D is calculated₀；

Wherein, the τ₂Calculation formula be：

The default impeller outlet width calculation formula is：

The default impeller inlet width calculation formula is：

7. the impeller method in the uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 according to claim 5, It is characterized in that, it is described according to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine the blade type of impeller The preliminary moulding of line, specifically includes：

It is flowed to according to blade angle along meridian flow direction in multinomial change profile, vane thickness along meridian in multinomial variation point Cloth, the β_1A, the β_2A, the d₂With default impeller blade import thickness d₁, determine the preliminary moulding of vane type line of impeller, In, the maximum value position of the blade angle is at 40% length of blade, the maximum value position of vane thickness is long in 80% blade At degree.

8. according to the leaf in the uniaxial CO2 compressor latter end Stage of claim 5-7 any one of them discharge coefficient 0.0056 Wheel method, which is characterized in that the predefined impeller cover side meridian molded line is an angled straight lines.

9. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is held by processor Following steps are realized during row：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller goes out Mouth established angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, determine that the meridional channel of impeller is preliminary Moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine that the vane type line of impeller is tentatively made Type；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, is more than with the Yplus values for ensureing wall surface first layer grid Or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimal orthogonality of grid Property be more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller adiabatic efficiency, work are carried out to the impeller pattern after mesh generation using CFD analysis softwares Condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, the condition range whether be more than or Equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, then adjusting the matching of meridional channel and the vane type line, modification blade beta angles are distributed, vane thickness is distributed, Meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default efficiency threshold until obtaining the impeller adiabatic efficiency Value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default blade flow field and sets The impeller of requirement is counted, the impeller is low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without being clearly separated；Leaf Without flow separation on the meridional channel of wheel；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade pressure surface It is uniform with the load changes in distribution on suction surface.

10. in a kind of 0.0056 uniaxial CO2 compressor latter end Stage of design claim 1-4 any one of them discharge coefficient Impeller unit, including memory, processor and storage on a memory and the computer program that can run on a processor, It is characterized in that, the processor realizes following steps when performing described program：

According to default impeller diameter D₂With default hub ratio ds/D₂, calculate hub diameter ds；

According to default impeller blade exit thickness d₂With one-dimensional thermodynamic computing as a result, determining vane inlet established angle β_1A, impeller goes out Mouth established angle β_2A, impeller outlet width b₂, impeller inlet width D₀；

According to the ds, the b₂, the D₀With predefined impeller cover side meridian molded line, determine that the meridional channel of impeller is preliminary Moulding；

According to the β_1A, the β_2A, the d₂, default impeller blade import thickness d₁, determine that the vane type line of impeller is tentatively made Type；

According to the preliminary moulding of the meridional channel and the preliminary moulding of the vane type line, impeller is primarily determined；

Obtained impeller is modeled, and extension processing is carried out to impeller inlet and outlet length；

Three-dimentional structured mesh division is carried out to obtained impeller pattern, is more than with the Yplus values for ensureing wall surface first layer grid Or equal to 1, while the length-width ratio of grid is more than or equal to 3000, extends than being more than or equal to 3, the minimal orthogonality of grid Property be more than 15 °；

Three Dimensional Viscous numerical simulation, analysis impeller adiabatic efficiency, work are carried out to the impeller pattern after mesh generation using CFD analysis softwares Condition range, blade flow field；

Whether impeller adiabatic efficiency described in comprehensive descision is more than or equal to default efficiency threshold, the condition range whether be more than or Equal to default condition range threshold value, and whether the blade flow field meets default blade flow field design requirement；

If it is not, then adjusting the matching of meridional channel and the vane type line, modification blade beta angles are distributed, vane thickness is distributed, Meridian molded line, re-using CFD are analyzed and are iterated, and are more than or equal to default efficiency threshold until obtaining the impeller adiabatic efficiency Value, the condition range are more than or equal to default condition range threshold value and the blade flow field meets default blade flow field and sets The impeller of requirement is counted, the impeller is low frictional resistance wedge shape impeller；

Wherein, the default blade flow field design requirement includes flowing of the impeller on the high section of different leaves without being clearly separated；Leaf Without flow separation on the meridional channel of wheel；The outlet flow angle of impeller is evenly distributed；Exit static pressure is evenly distributed；Blade pressure surface It is uniform with the load changes in distribution on suction surface.