CN113297706B - Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index - Google Patents

Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index Download PDF

Info

Publication number
CN113297706B
CN113297706B CN202110679175.XA CN202110679175A CN113297706B CN 113297706 B CN113297706 B CN 113297706B CN 202110679175 A CN202110679175 A CN 202110679175A CN 113297706 B CN113297706 B CN 113297706B
Authority
CN
China
Prior art keywords
compressor
design
predicting
rotating speed
pressure ratio
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
Application number
CN202110679175.XA
Other languages
Chinese (zh)
Other versions
CN113297706A (en
Inventor
高丽敏
刘锬韬
吴锋
刘涛
冯旭栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
AECC Sichuan Gas Turbine Research Institute
Original Assignee
Northwestern Polytechnical University
AECC Sichuan Gas Turbine Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University, AECC Sichuan Gas Turbine Research Institute filed Critical Northwestern Polytechnical University
Publication of CN113297706A publication Critical patent/CN113297706A/en
Application granted granted Critical
Publication of CN113297706B publication Critical patent/CN113297706B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Control Of Positive-Displacement Air Blowers (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

The invention relates to a method for predicting efficiency points of various rotating speed peaks of a gas compressor according to design, and belongs to the field of gas compressor characteristic prediction. The data input in the prediction process is only the design index (the total pressure ratio and the efficiency of a design point) of the compressor, and the blade profile data of the compressor is not needed. According to the total pressure ratio of the design point, three coefficients are firstly constructed, and then the peak efficiency point characteristics of all equal rotating speed lines are predicted. The formula and parameters used in the invention are derived from the mining of historical compressor characteristic data. The prediction method provided by the invention can be used for predicting the peak efficiency points of the rotating speeds of the axial flow compressor and the centrifugal compressor, is particularly suitable for predicting the peak efficiency points of the rotating speeds of the axial flow compressor and the centrifugal compressor when no blade profile data exists at the initial stage of the overall design of the gas turbine and the aero-engine, and has the advantages of simple and convenient calculation and higher precision.

Description

Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index
Technical Field
The invention belongs to the field of compressor characteristic prediction, and relates to a method for predicting a compressor characteristic diagram, in particular to prediction of equal-speed-line peak efficiency points of the compressor characteristic diagram.
Background
The compressor is an important part of a gas turbine, is used for compressing incoming flow gas, and is widely applied to the fields of aviation power, ship power, thermal power generation and the like. The characteristic diagram of the gas compressor is a main representation mode of the gas performance of the gas compressor in the full operation range, and is a key basis for the actual operation of an aircraft engine and a gas turbine, the characteristic diagram describes the relationship between the efficiency and the total pressure ratio of the gas compressor and the dimensionless normalized rotating speed and the dimensionless normalized flow rate, and is generally drawn into two line graphs comprising a plurality of equal rotating speed lines, each equal rotating speed line takes the converted flow rate as an independent variable, the efficiency and the total pressure ratio in the two line graphs are respectively a dependent variable, and each equal rotating speed line generally has a peak efficiency point. When a gas turbine and an aero-engine work, the working line of a gas compressor is generally close to the connecting line of peak efficiency points, the peak efficiency points are the optimal points of the working state of the gas compressor at the rotating speed, and prediction of the peak efficiency points is of great significance to prediction of gas turbine operation and full-range characteristics.
The method generally comprises the steps of obtaining a characteristic diagram of the compressor, and measuring a test of the compressor and calculating a numerical value of a flow field of the compressor through multiple methods, wherein the two methods both need known blade profile data of the compressor. In engineering practice, in the initial design stage of the gas turbine engine, initial compressor characteristic data is needed when overall performance calculation is carried out, and at the moment, compressor blade profile data generally do not exist, so that only a performance prediction algorithm independent of the blade profile data can be used.
At present, the performance prediction of the compressor is mostly based on an interpolated value or an extrapolated rotation speed range which is not included in an existing characteristic diagram, or based on the positive and negative design of a compressor blade profile according to design indexes, and then characteristics are obtained through the calculation of flow field numbers, and the calculated amount is large.
Disclosure of Invention
Technical problem to be solved
In order to avoid the defects of the prior art, the invention provides the compressor characteristic prediction method which does not need calculation of geometric data of the blade shape and numerical values of a flow field, and only needs design indexes and simple calculation.
Technical scheme
A method for predicting the peak efficiency point of a characteristic line of a gas compressor according to design indexes comprises the following steps: total pressure ratio pi of highest efficiency of rotation speed D And efficiency η D Calculating the normalized flow m of each rotation speed peak efficiency point p Total pressure ratio of pi p And efficiency η p The calculation formula is as follows:
Figure BDA0003122172520000021
Figure BDA0003122172520000022
Figure BDA0003122172520000023
Figure BDA0003122172520000024
Figure BDA0003122172520000025
γ＝1.4
wherein n is the rotation speed, the coefficient k T 、k p And k h Are all design index total pressure ratio pi D The coefficient correlation is as follows:
k T ＝0.068exp(4.86/π D )+0.98
k P ＝0.058exp(5.65/π D )+1.157
k h ＝0.52exp(2.41/π D )+0.078。
the coefficient k T 、k p And k h A float in the range of 5% can be performed.
A computer system, comprising: one or more processors, a computer readable storage medium, for storing one or more programs, which when executed by the one or more processors, cause the one or more processors to implement the above-described method.
A computer-readable storage medium storing computer-executable instructions for performing the above-described method when executed.
A computer program comprising computer executable instructions for implementing the above method when executed.
Advantageous effects
According to the method for predicting the efficiency points of the rotating speed peak values of the gas compressor according to the design indexes, the data input in the prediction process is only the design indexes (the design point pressure ratio and the efficiency) of the gas compressor, and the blade profile data of the gas compressor is not needed. According to the design point pressure ratio, three coefficients are firstly constructed, and then the peak efficiency point characteristics of all equal rotating speed lines are predicted. The formula and parameters used in the invention are derived from the mining of historical compressor characteristic data. The prediction method provided by the invention can be used for predicting the peak efficiency points of the rotating speeds of the axial flow compressor and the centrifugal compressor, is particularly suitable for predicting the peak efficiency points of the rotating speeds of the axial flow compressor and the centrifugal compressor when no blade profile data exists at the initial stage of the overall design of the gas turbine and the aero-engine, and has the advantages of simple and convenient calculation and higher precision. The beneficial effects are as follows:
1. the method does not need the blade profile geometric parameters of the compressor and only depends on design indexes for prediction. 2. The calculated amount is small, the flow field numerical value calculation is not needed, and the adjustment range of the empirical parameters is small.
Drawings
The drawings are only for purposes of illustrating particular embodiments and are not to be construed as limiting the invention, wherein like reference numerals are used to designate like parts throughout.
FIG. 1 is a flow chart of a method for predicting the characteristics of efficiency points of different rotational speeds of a compressor depending on design indexes;
fig. 2 is a comparison between the predicted effect and the actual characteristic diagram according to the design index of a certain type of compressor, where a triangle is a predicted value, a circle is an original actual value, and a solid line is an original characteristic diagram.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
In the process of researching a large number of compressor characteristic graphs, the following results are found: the design points of different rotating speeds of the compressor present some special rules, and the prediction of each rotating speed peak efficiency point of the compressor can be realized by using the rules. The specific rule is that the specific work delta h corresponding to the peak efficiency point of each rotating speed in the working rotating speed range of the air compressor * The input torque T and the input power P are linearly related to the total pressure ratio pi, all variables are correspondingly normalized by the designed rotating speed peak value efficiency point, and the linear relation can be expressed as follows:
Δh * /Δh * D ＝k h (π/π D )+1-k h
T * /T * D ＝k T (π/π D )+1-k T
P * /P * D ＝k P (π/π D )+1-k P
further research shows that the proportionality coefficient and the total pressure ratio of the design point of the compressor are in an exponential relationship, namely, the expression in the first step is as follows:
k T ＝0.068exp(4.86/π D )+0.98
k P ＝0.058exp(5.65/π D )+1.157
k h ＝0.52exp(2.41/π D )+0.078
by utilizing the relation formula, the total pressure ratio pi is designed D After obtaining the above three proportionality coefficients, the relative torque T is used in consideration of the design point of the rotation speed n =1 and the flow rate m =1 * /T * d Relative power P * /P * d Relative ratio work Δ h * /Δh * d And the physical relation among variable rotating speed n, flow m, isentropic efficiency eta and total pressure ratio pi in a compressor characteristic diagram is as follows:
Figure BDA0003122172520000041
Figure BDA0003122172520000042
Figure BDA0003122172520000043
the flow m of each rotating speed peak efficiency point can be deduced p Total pressure ratio pi p And efficiency η p Namely, the expression in step two is as follows:
Figure BDA0003122172520000051
Figure BDA0003122172520000052
Figure BDA0003122172520000053
wherein:
Figure BDA0003122172520000054
Figure BDA0003122172520000055
γ＝1.4
in order that those skilled in the art will better understand the present invention, the following detailed description is given with reference to specific examples.
The prediction is carried out by taking a certain single-stage compressor as a true value, and the design parameter pi of the compressor is D ＝1.9195，η D =0.85123, step one specifically includes: three coefficients k T 、k p And k h Are all design index total pressure ratio pi D The relationship is as follows:
k T ＝0.068exp(4.86/π D )+0.98
k P ＝0.058exp(5.65/π D )+1.157
k h ＝0.52exp(2.41/π D )+0.078
the three expressions are provided after researching and analyzing a plurality of gas compressor characteristic diagrams with different pressure ratios, and the physical significance is as follows: at the highest efficiency point of the equal rotating speed line, the slope of the dimensionless linear relation of torque, power and specific work with the total pressure ratio at the design point, and the slope parameter has an exponential relation with the design pressure ratio. k is a radical of T 、k p And k h The three parameters represent the difference between different design configurations of the same design index, and the difference is calculated by substituting the formula:
k T ＝1.8352
k p ＝2.2579
k h ＝1.9031
the three parameters can float within about 5% according to experience and historical data, and after floating by referring to a prototype compressor, the parameter values are as follows:
k T ＝1.872
k p ＝2.383
k h ＝1.822
the second step specifically comprises: normalized flow m of each rotation speed line peak efficiency point p Total pressure ratio pi p And efficiency η p The relation among the rotating speed n, the design index and the three coefficients in the second step is as follows:
Figure BDA0003122172520000061
Figure BDA0003122172520000062
Figure BDA0003122172520000063
wherein:
Figure BDA0003122172520000064
Figure BDA0003122172520000065
γ＝1.4
FIG. 2 shows a comparison graph of a prediction result and a true value, a triangle is a predicted value, a circle is a true value, the predicted value is very close to the true value, and only a small difference exists, so that the accuracy of the prediction algorithm is high.
While the invention has been described with reference to specific embodiments, the invention is not limited thereto, and various equivalent modifications or substitutions can be easily made by those skilled in the art within the technical scope of the present disclosure.

Claims (4)

1. A method for predicting the peak efficiency point of a characteristic line of a compressor according to design indexes is characterized in that according to the design indexes: total pressure ratio pi of highest efficiency point of rotation speed D And efficiency η D Calculating the normalized flow m of each rotation speed peak efficiency point p Total pressure ratio pi p And efficiency η p The calculation formula is as follows:
Figure FDA0003922428270000011
Figure FDA0003922428270000012
Figure FDA0003922428270000013
Figure FDA0003922428270000014
Figure FDA0003922428270000015
γ＝1.4
wherein n is the rotation speed, the coefficient k T 、k p And k h Are all design indexes and total pressure ratio pi D The coefficient correlation is as follows:
k T ＝0.068exp(4.86/π D )+0.98
k P ＝0.058exp(5.65/π D )+1.157
k h ＝0.52exp(2.41/π D )+0.078。
2. the method of claim 1, wherein the coefficient k is a function of a design criterion for predicting a peak efficiency point of a compressor characteristic line T 、k p And k h A float in the range of 5% can be performed.
3. A computer system, comprising: one or more processors, a computer readable storage medium, for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the method of claim 1.
4. A computer-readable storage medium having stored thereon computer-executable instructions, which when executed, perform the method of claim 1.
CN202110679175.XA 2021-04-22 2021-06-18 Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index Active CN113297706B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110437583 2021-04-22
CN2021104375834 2021-04-22

Publications (2)

Publication Number Publication Date
CN113297706A CN113297706A (en) 2021-08-24
CN113297706B true CN113297706B (en) 2023-02-21

Family

ID=77328787

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110679175.XA Active CN113297706B (en) 2021-04-22 2021-06-18 Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index

Country Status (1)

Country Link
CN (1) CN113297706B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114135397B (en) * 2021-11-05 2023-04-25 广州珠江天然气发电有限公司 Diagnosis optimization method, device and system for gas-steam combined cycle performance

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4539841A (en) * 1984-02-13 1985-09-10 General Motors Corporation Method of determining engine cylinder compression pressure and power output
CN105930662A (en) * 2016-04-26 2016-09-07 中国科学院工程热物理研究所 Low-speed characteristic extrapolation method for gas compressor
CN108301955A (en) * 2018-01-15 2018-07-20 武汉大学 Axial-flow type PAT power generation modes optimum efficiency point parameter and performance curve predictor method
CN108953232A (en) * 2018-07-20 2018-12-07 大连海事大学 A kind of non-axisymmetric distribution stator blade axial-flow compressor
CN109684597A (en) * 2018-12-21 2019-04-26 哈尔滨工程大学 A method of it is expanded for compressor whole working condition property
CN109800519A (en) * 2019-01-25 2019-05-24 深圳友铂科技有限公司 A kind of aerofoil fan multiple spot Aerodynamic optimization design method
CN110489877A (en) * 2019-08-21 2019-11-22 中国航发沈阳发动机研究所 A kind of interpolation method suitable for Aeroengine Real Time Model
CN110701086A (en) * 2019-09-10 2020-01-17 天津大学 Compressor full-working-condition performance prediction method
CN111159909A (en) * 2019-12-31 2020-05-15 中国船舶重工集团公司第七0三研究所 Gas compressor characteristic data processing and characteristic cloud picture drawing method
CN112270041A (en) * 2020-10-27 2021-01-26 中国船舶重工集团公司第七0三研究所 Method for drawing distribution of matching points of compressor master model machine of ship gas turbine
CN112594064A (en) * 2020-11-25 2021-04-02 北京航空航天大学 S2 flow field diagnosis method based on interstage measurement parameters of axial flow compressor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203364B (en) * 2017-05-26 2020-12-22 哈尔滨工程大学 Prediction and identification method for full-working-condition characteristics of gas compressor

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4539841A (en) * 1984-02-13 1985-09-10 General Motors Corporation Method of determining engine cylinder compression pressure and power output
CN105930662A (en) * 2016-04-26 2016-09-07 中国科学院工程热物理研究所 Low-speed characteristic extrapolation method for gas compressor
CN108301955A (en) * 2018-01-15 2018-07-20 武汉大学 Axial-flow type PAT power generation modes optimum efficiency point parameter and performance curve predictor method
CN108953232A (en) * 2018-07-20 2018-12-07 大连海事大学 A kind of non-axisymmetric distribution stator blade axial-flow compressor
CN109684597A (en) * 2018-12-21 2019-04-26 哈尔滨工程大学 A method of it is expanded for compressor whole working condition property
CN109800519A (en) * 2019-01-25 2019-05-24 深圳友铂科技有限公司 A kind of aerofoil fan multiple spot Aerodynamic optimization design method
CN110489877A (en) * 2019-08-21 2019-11-22 中国航发沈阳发动机研究所 A kind of interpolation method suitable for Aeroengine Real Time Model
CN110701086A (en) * 2019-09-10 2020-01-17 天津大学 Compressor full-working-condition performance prediction method
CN111159909A (en) * 2019-12-31 2020-05-15 中国船舶重工集团公司第七0三研究所 Gas compressor characteristic data processing and characteristic cloud picture drawing method
CN112270041A (en) * 2020-10-27 2021-01-26 中国船舶重工集团公司第七0三研究所 Method for drawing distribution of matching points of compressor master model machine of ship gas turbine
CN112594064A (en) * 2020-11-25 2021-04-02 北京航空航天大学 S2 flow field diagnosis method based on interstage measurement parameters of axial flow compressor

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Joint prediction of aircraft gas turbine engine axial flow compressor off-design performance and surge line based on the expanded method of generalized functions;A. B. Mikhailova;D. A. Akhmedzyanov;Yu. M. Akhmetov;A. E. Mikhai;《Russian Aeronautics》;20141231;第57卷(第3期);第291-296页 *
几何参数变化对离心压气机性能的影响;靳军等;《科学技术与工程》;20141218(第35期);第302-308页 *
基于转静干涉效应的压气机叶片气动载荷分析;杨文军等;《工程设计学报》;20181028(第05期);第79-87页 *
组合压气机性能和气动稳定性计算分析;陈靖华等;《航空计算技术》;20141125(第06期);第45-47页 *

Also Published As

Publication number Publication date
CN113297706A (en) 2021-08-24

Similar Documents

Publication Publication Date Title
CN113297705B (en) Method for predicting characteristics of gas compressor according to design indexes
CN109684597B (en) Method for expanding all-condition characteristics of gas compressor
CN108153998B (en) Full three-dimensional pneumatic optimization design method for centrifugal blower impeller
US20170306970A1 (en) Extrapolation method of low rotational speed characteristic of compressor
CN113297706B (en) Method for predicting efficiency point of each rotating speed peak value of gas compressor according to design index
Nishida et al. Performance improvement of a return channel in a multistage centrifugal compressor using multiobjective optimization
CN112685851B (en) Multi-stage axial flow compressor pneumatic design method based on key dimensionless load control parameters
CN107203364B (en) Prediction and identification method for full-working-condition characteristics of gas compressor
CN110489877B (en) Interpolation method suitable for real-time model of aero-engine
Wang et al. Prediction method for low speed characteristics of compressor based on modified similarity theory with genetic algorithm
CN110705079B (en) Centrifugal compressor structure optimization method based on simulated annealing algorithm
CN108491572B (en) Two-dimensional flow simulation method for turbomachinery based on two-sided flow theory
CN112699503B (en) Method for designing inverse problem of S2 of axial flow compressor based on dimensionless load control parameters
CN112100862A (en) Novel turbine isentropic efficiency fast solving method based on variable specific heat capacity
CN110866312A (en) Format conversion method for turbine characteristics of aviation gas turbine engine
CN113836642B (en) Low-rotation-speed characteristic expansion method for aero-engine component
Li et al. Development and application of a throughflow method for high-loaded axial flow compressors
CN115688319A (en) Variable cycle compression system test matching characteristic modeling method
Blinov et al. Multi-criteria optimisation of subsonic axial compressor blading
Wang et al. Method of extrapolating low speed compressor curves based on improved similarity laws
CN115186443B (en) Pneumatic design method for multistage power turbine of power generation type gas turbine for inverted bathtub load ship
CN115186442B (en) Multi-stage power turbine pneumatic design method for power generation type gas turbine for decreasing load ship
CN115186441B (en) Pneumatic design method for multistage power turbine of power generation type gas turbine for incremental load ship
CN115186440B (en) Pneumatic design method for two-stage high-speed power turbine of marine power generation type gas turbine
CN114912227B (en) Non-steady-state similarity method for centripetal turbine

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