CN115292946A - High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation - Google Patents

High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation Download PDF

Info

Publication number
CN115292946A
CN115292946A CN202210977690.0A CN202210977690A CN115292946A CN 115292946 A CN115292946 A CN 115292946A CN 202210977690 A CN202210977690 A CN 202210977690A CN 115292946 A CN115292946 A CN 115292946A
Authority
CN
China
Prior art keywords
pressure
turbine
total
pressure turbine
efficiency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210977690.0A
Other languages
Chinese (zh)
Other versions
CN115292946B (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.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine 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 AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN202210977690.0A priority Critical patent/CN115292946B/en
Publication of CN115292946A publication Critical patent/CN115292946A/en
Application granted granted Critical
Publication of CN115292946B publication Critical patent/CN115292946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Control Of Turbines (AREA)

Abstract

The application belongs to the technical field of engine tests, and particularly relates to a high-pressure turbine efficiency evaluation method and device based on variable specific heat calculation. The method comprises the steps of S1, calculating a first high-pressure turbine efficiency initial guess value according to a high-pressure turbine efficiency formula; s2, dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the total inlet temperature and the total outlet temperature of each small turbine; s3, calculating the adiabatic index of each small turbine; s4, determining the total outlet pressure of each small turbine; s5, repeating the steps until the total outlet pressure of the last small turbine is calculated; and S6, determining the change rate of the outlet total pressure of the last small turbine relative to the outlet total pressure of the high-pressure turbine, and if the change rate is greater than a preset value, updating the initial guess value of the efficiency of the first high-pressure turbine until the change rate is smaller than the preset value. The method and the device realize the calculation of the ratio of heat to heat of the high-pressure turbine efficiency, and greatly improve the calculation precision of the high-pressure turbine efficiency.

Description

High-pressure turbine efficiency evaluation method and device based on heat-variable ratio calculation
Technical Field
The application belongs to the technical field of engine tests, and particularly relates to a high-pressure turbine efficiency evaluation method and device based on heat-variable ratio calculation.
Background
In order to evaluate the efficiency of components such as a fan, a compressor, a turbine and the like of an aircraft engine, the total temperature and the total pressure of an inlet and an outlet of the corresponding component are obtained in a measuring or calculating mode, then a constant engine heat insulation index kg (1.33 for gas and 1.4 for air) is given, and the efficiency of the corresponding component is calculated according to the engine principle.
By adopting the method, the efficiency of corresponding parts of the engine can be preliminarily evaluated and calculated, the variation range of the inlet and outlet temperatures of a fan, a gas compressor and the like is small, the calculation accuracy is acceptable, but the variation of the inlet and outlet temperatures of high-pressure turbine parts is large (reduced from thousands of hundreds of K to hundreds of K), the difference between the adiabatic index at the inlet of the high-pressure turbine and the adiabatic index at the outlet is large because the engine adiabatic index kg is related to the gas temperature and the oil-gas ratio, the original technical scheme can only adopt the adiabatic index kg at the inlet or the adiabatic index kg at the outlet, the calculation accuracy is relatively low, the average value of the two is adopted for calculation when the accuracy is slightly high, but the calculation accuracy also has large errors.
Disclosure of Invention
In order to solve one of the problems, the application provides a method and a device for evaluating the efficiency of a high-pressure turbine based on the calculation of the heat-variation ratio, so that the efficiency of the high-pressure turbine can be evaluated more accurately.
The application provides a high-pressure turbine efficiency evaluation method based on heat-to-variable ratio calculation in a first aspect, which mainly comprises the following steps:
step S1, calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating an oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow;
s2, enabling the high-pressure turbine to flow from an inlet to an outletDividing the inlet into a plurality of small turbines, and determining the total inlet temperature T of each small turbine according to the total inlet temperature and the total outlet temperature of the high-pressure turbine and the equal-interval cooling principle in (i) And total outlet temperature T out (i);
S3, according to the total inlet temperature T of the ith small turbine in (i) And the oil-gas ratio determines the adiabatic index k of the ith small turbine g (i);
Step S4, determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And total pressure at the inlet P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the total inlet pressure of the first small turbine is the total inlet pressure P of the high-pressure turbine 4 The total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine;
step S5, repeating the steps S2 to S4 until the total outlet pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out
S6, determining the total outlet pressure P of the last small turbine out (num) in relation to the total pressure P at the outlet of the high-pressure turbine 5 If the change rate is smaller than the preset value, the first high-pressure turbine efficiency initial guess value given in the step S1 is used as the final value of the high-pressure turbine efficiency, otherwise, the first high-pressure turbine efficiency initial guess value in the step S1 is updated, and the steps are repeated until the change rate is smaller than the preset value.
Preferably, in step S2, the inlet total temperature T of the ith small turbine is in (i) The total outlet temperature of the ith-1 th small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and the number of the small turbines divided by the high-pressure turbine is num.
Preferably, in step S4, the total outlet pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure BDA0003799008780000021
preferably, in step S6, updating the first high-pressure turbine efficiency initial guess value includes:
step S61, setting a first guess value eta of the second high-pressure turbine efficiency THtmp Setting multiple of the first high-pressure turbine efficiency initial guess value;
step S62, according to the first guess value eta of the second high-pressure turbine efficiency TH Performing an inner loop of steps S2-S5, determining the total outlet pressure of the last small turbine, and assigning a second variable P 5outtmp
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing an outer loop of the steps S2 to S6:
Figure BDA0003799008780000022
wherein eta is TH(new) The value is guessed for the new first high pressure turbine efficiency.
Preferably, in step S61, the set multiple is 1.01 times.
The second aspect of the present application provides a high-pressure turbine efficiency evaluation device based on heat-to-variable ratio calculation, mainly including:
the efficiency and oil-gas ratio calculation module is used for calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating an oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow;
the high-pressure turbine dividing module is used for dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the total inlet temperature T of each small turbine according to the total inlet temperature and the total outlet temperature of the high-pressure turbine and the equal-interval cooling principle in (i) And total outlet temperature T out (i);
A small turbine adiabatic index calculation module for calculating the total inlet temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k for the ith small turbine g (i);
Each small turbine outlet total pressure calculation module is used for determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And inlet total pressure P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the total inlet pressure of the first small turbine is the total inlet pressure P of the high-pressure turbine 4 The total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine;
the last small turbine outlet total pressure calculation module is used for repeatedly calling the modules until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out
A residual-based loop control module for determining the total outlet pressure P of the last small turbine out (num) in relation to the total pressure P at the outlet of the high-pressure turbine 5 If the change rate is smaller than the preset value, the given first high-pressure turbine efficiency initial guess value is used as the final value of the high-pressure turbine efficiency, otherwise, the first high-pressure turbine efficiency initial guess value is updated until the change rate is smaller than the preset value.
Preferably, in the high-pressure turbine division module, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature of the ith-1 th small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and the number of the small turbines divided by the high-pressure turbine is num.
Preferably, in the total outlet pressure calculation module of each small turbine, the total outlet pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure BDA0003799008780000031
preferably, the residual-based loop control module includes:
a second high-pressure turbine efficiency initial guess value setting unit for setting a second high-pressure turbine efficiency initial guess value eta THtmp Setting multiple of the first high-pressure turbine efficiency initial guess value;
an internal circulation unit for initially guessing the value eta according to the efficiency of the second high-pressure turbine TH Determining the total outlet pressure of the last small turbine and assigning a second variable P 5outtmp
An external circulation unit for updating the first high pressure turbine efficiency initial guess value:
Figure BDA0003799008780000041
wherein eta is TH(new) And guessing the new first high-pressure turbine efficiency.
Preferably, the set multiple is 1.01 times.
This application decomposes high-pressure turbine into a plurality of "little turbine" according to certain temperature drop interval and calculates, realizes the variable specific heat ratio of high-pressure turbine efficiency and calculates, the computational accuracy of improvement high-pressure turbine efficiency that can be very big.
Drawings
FIG. 1 is a flow chart of a preferred embodiment of the method for estimating the efficiency of a high pressure turbine based on the heat-transfer ratio calculation according to the present application.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present application clearer, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the drawings in the embodiments of the present application. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are some, but not all embodiments of the present application. The embodiments described below with reference to the accompanying drawings are exemplary and intended to be used for explaining the present application and should not be construed as limiting the present application. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments in the present application without making any creative effort belong to the protection scope of the present application. Embodiments of the present application will be described in detail below with reference to the drawings.
Generally, in order to obtain high-pressure turbine efficiency with higher precision, a heat-variable specific heat calculation formula (1) is adopted to calculate an adiabatic index kg according to an oil-gas ratio and a total temperature of gas, and then the efficiency of the high-pressure turbine can be calculated according to a formula (2) according to total temperature and total pressure of an inlet and an outlet of the high-pressure turbine obtained by measurement or calculation.
kg = f (oil-gas ratio, total temperature of 823060, 82308230, 823082308230823082308230823082308260 82308230; \8230; \ 8230; \ 8230;' 8230; (1)
Figure BDA0003799008780000042
Wherein:
t4 and P4 are total temperature and total pressure of an inlet of the high-pressure turbine;
t5 and P5 are total temperature and total pressure of an outlet of the high-pressure turbine;
kg is the adiabatic index.
In order to calculate the efficiency of the high-pressure turbine more accurately, the calculation process of the high-pressure turbine can be divided into an infinite number of small turbines for calculation. And each small turbine is divided according to a fixed temperature drop, the total outlet temperature and the total pressure of each small turbine are sequentially calculated according to the high-pressure turbine efficiency initial value calculated by the formula (2), finally, a residual error is constructed according to the high-pressure turbine outlet pressure, and the high-pressure turbine efficiency initial value is updated and iterated until the residual error meets the error requirement. Specifically, as shown in fig. 1, a first aspect of the present application provides a method for estimating efficiency of a high-pressure turbine based on heat-variation ratio calculation, which mainly includes:
step S1, calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH And calculating the oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow.
Step S2, applying the high pressureThe turbine is divided into a plurality of small turbines from an inlet to an outlet, and the total inlet temperature T of each small turbine is determined according to the total inlet temperature and the total outlet temperature of the high-pressure turbine and the equal interval cooling principle in (i) And total outlet temperature T out (i)。
In some alternative embodiments, in step S2, the total inlet temperature T of the ith small turbine in (i) The total outlet temperature of the ith-1 small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and num is the number of small turbines divided by the high-pressure turbine.
S3, according to the total inlet temperature T of the ith small turbine in (i) And the oil-gas ratio determines the adiabatic index k of the ith small turbine g (i)。
Step S4, determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And inlet total pressure P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the inlet total pressure of the first small turbine is the inlet total pressure P of the high-pressure turbine 4 And the total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine.
In some alternative embodiments, in step S4, the total outlet pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure BDA0003799008780000051
step S5, repeating the steps S2-S4 until the total outlet pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out (ii) a I.e. P 5out =P out (num)。
S6, determining the total outlet pressure P of the last small turbine out Relative height of (num)Total pressure P at outlet of pressure turbine 5 If the change rate is smaller than the preset value, the first high-pressure turbine efficiency initial guess value given in the step S1 is used as the final value of the high-pressure turbine efficiency, otherwise, the first high-pressure turbine efficiency initial guess value in the step S1 is updated, and the steps are repeated until the change rate is smaller than the preset value.
In some alternative embodiments, updating the first high-pressure turbine efficiency preliminary guess value in step S6 includes:
step S61, setting the initial guess value eta of the efficiency of the second high-pressure turbine THtmp Setting multiple of the first high-pressure turbine efficiency initial guess value;
step S62, according to the first guess value eta of the second high-pressure turbine efficiency TH Performing an inner loop of steps S2-S5, determining the total outlet pressure of the last small turbine, and assigning a second variable P 5outtmp
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing an outer loop of the steps S2 to S6:
Figure BDA0003799008780000061
wherein eta is TH(new) The value is guessed for the new first high pressure turbine efficiency.
In some alternative embodiments, in step S61, the set multiple is 1.01 times.
In this example, P is compared 5out And P 5 Error of (2), e.g.
Figure BDA0003799008780000062
The calculation is terminated; such as
Figure BDA0003799008780000063
Then set η THtmp =η TH *1.01, the total outlet pressure P of the last small turbine is resumed out (num) is calculated and assigned to P 5outtmp (ii) a Then press against
Figure BDA0003799008780000064
Updating eta TH . Finally, the process is repeated until the new eta is obtained TH The calculated turbine outlet pressure residual meets the error requirement of step S6.
The second aspect of the present application provides a high-pressure turbine efficiency evaluation device based on the heat-variation ratio calculation corresponding to the above method, which mainly includes:
an efficiency and gas-oil ratio calculation module for calculating the first high-pressure turbine efficiency initial guess value eta according to the high-pressure turbine efficiency formula TH Calculating the oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow;
the high-pressure turbine dividing module is used for dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the total inlet temperature T of each small turbine according to the total inlet temperature and the total outlet temperature of the high-pressure turbine and the equal-interval cooling principle in (i) And total outlet temperature T out (i);
A small turbine adiabatic index calculation module for calculating the total inlet temperature T of the ith small turbine in (i) And the oil-gas ratio determines the adiabatic index k of the ith small turbine g (i);
Each small turbine outlet total pressure calculation module is used for determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And inlet total pressure P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the total inlet pressure of the first small turbine is the total inlet pressure P of the high-pressure turbine 4 The total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine;
the last small turbine outlet total pressure calculation module is used for repeatedly calling the modules until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out
A residual-based loop control module for determining the total outlet pressure P of the last small turbine out (num) in relation to the total pressure P at the outlet of the high-pressure turbine 5 If the change rate is less than the predetermined valueAnd if not, updating the first high-pressure turbine efficiency initial guess value until the change rate is smaller than a preset value.
In some alternative embodiments, in the high pressure turbine division module, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature of the ith-1 th small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and num is the number of small turbines divided by the high-pressure turbine.
In some alternative embodiments, in the total outlet pressure calculation module of each small turbine, the total outlet pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure BDA0003799008780000071
in some optional embodiments, the residual based loop control module comprises:
a second high-pressure turbine efficiency initial guess value setting unit for setting a second high-pressure turbine efficiency initial guess value eta THtmp Setting multiple of the first high-pressure turbine efficiency initial guess value;
an internal circulation unit for initially guessing the value eta according to the efficiency of the second high-pressure turbine TH Determining the total outlet pressure of the last small turbine and assigning a second variable P 5outtmp
An outer loop unit for updating the first high pressure turbine efficiency initial guess value:
Figure BDA0003799008780000072
wherein eta is TH(new) Initial guess value for first high pressure turbine efficiency。
In some alternative embodiments, the set multiple is 1.01 times.
This application decomposes high-pressure turbine into a plurality of "little turbine" according to certain temperature drop interval and calculates, realizes the variable specific heat ratio of high-pressure turbine efficiency and calculates, the computational accuracy of improvement high-pressure turbine efficiency that can be very big.
Although the present application has been described in detail with respect to the general description and specific embodiments, it will be apparent to those skilled in the art that certain modifications or improvements may be made based on the present application. Accordingly, such modifications and improvements are intended to be within the scope of this invention as claimed.

Claims (10)

1. A high-pressure turbine efficiency evaluation method based on heat of variation calculation is characterized by comprising the following steps:
step S1, calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating an oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow;
s2, dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the total inlet temperature T of each small turbine according to the equal interval cooling principle according to the total inlet temperature and the total outlet temperature of the high-pressure turbine in (i) And total outlet temperature T out (i);
S3, according to the total inlet temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k for the ith small turbine g (i);
Step S4, determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And total pressure at the inlet P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the inlet total pressure of the first small turbine is the inlet total pressure P of the high-pressure turbine 4 The total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine;
step S5, repeating the stepS2-S4, until the total outlet pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out
S6, determining the total outlet pressure P of the last small turbine out (num) relative to the total pressure at the outlet of the high-pressure turbine P 5 If the change rate is smaller than the preset value, the first high-pressure turbine efficiency initial guess value given in the step S1 is used as the final value of the high-pressure turbine efficiency, otherwise, the first high-pressure turbine efficiency initial guess value in the step S1 is updated, and the steps are repeated until the change rate is smaller than the preset value.
2. The method for evaluating the efficiency of a high-pressure turbine based on the calculation of the specific heat capacity according to claim 1, wherein in step S2, the total inlet temperature T of the ith small turbine is in (i) The total outlet temperature of the ith-1 small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and num is the number of small turbines divided by the high-pressure turbine.
3. The method for estimating the efficiency of a high-pressure turbine based on a variable specific heat calculation as claimed in claim 1, wherein in step S4, the total outlet pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure FDA0003799008770000011
4. a method for estimating the efficiency of a high pressure turbine based on the calculation of the specific heat variation according to claim 1, wherein the updating the first initial guess value of the efficiency of the high pressure turbine in step S6 includes:
step S61, setting a first guess value eta of the second high-pressure turbine efficiency THtmp Is the first high pressure vortexSetting multiple of the initial guess value of the wheel efficiency;
step S62, according to the first guess value eta of the second high-pressure turbine efficiency TH Performing an inner loop of steps S2-S5, determining the total outlet pressure of the last small turbine, and assigning a second variable P 5outtmp
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing an outer loop of the steps S2 to S6:
Figure FDA0003799008770000021
wherein eta is TH(new) The value is guessed for the new first high pressure turbine efficiency.
5. The method for evaluating a high-pressure turbine efficiency based on a calculation of a specific heat capacity according to claim 1, wherein in step S61, the set multiple is 1.01 times.
6. A high-pressure turbine efficiency evaluation device based on a variable specific heat calculation, characterized by comprising:
the efficiency and oil-gas ratio calculation module is used for calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating an oil-gas ratio according to the air flow at the inlet of the high-pressure turbine and the fuel flow;
the high-pressure turbine dividing module is used for dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the total inlet temperature T of each small turbine according to the total inlet temperature and the total outlet temperature of the high-pressure turbine and the equal-interval cooling principle in (i) And total outlet temperature T out (i);
A small turbine adiabatic index calculation module for calculating the total inlet temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k for the ith small turbine g (i);
Each small turbine outlet total pressure calculation module is used for determining the turbine efficiency eta of the ith small turbine TH(i) And based on the total inlet temperature T of the ith small turbine in (i) Total outlet temperature T out (i) And total pressure at the inlet P in (i) Determining the total outlet pressure P of the ith small turbine out (i) Wherein the total inlet pressure of the first small turbine is the total inlet pressure P of the high-pressure turbine 4 The total outlet pressure of the ith small turbine is the total inlet pressure of the (i + 1) th small turbine;
the last small turbine outlet total pressure calculation module is used for repeatedly calling the modules until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out
A residual-based loop control module for determining the total outlet pressure P of the last small turbine out (num) relative to the total pressure at the outlet of the high-pressure turbine P 5 If the change rate is smaller than the preset value, the given first high-pressure turbine efficiency initial guess value is used as the final value of the high-pressure turbine efficiency, otherwise, the first high-pressure turbine efficiency initial guess value is updated until the change rate is smaller than the preset value.
7. A high-pressure turbine efficiency evaluation device based on specific heat variation calculation as claimed in claim 6, wherein in the high-pressure turbine division module, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature of the ith-1 th small turbine and the total outlet temperature T of the ith small turbine out (i) Comprises the following steps:
T out (i)=T 4 -ΔT*i;
wherein, T 4 Δ T = (T) for total high pressure turbine inlet temperature 5 -T 4 ) ,/num, wherein, T 5 And num is the total outlet temperature of the high-pressure turbine, and num is the number of small turbines divided by the high-pressure turbine.
8. The variable specific heat calculation-based high-pressure turbine efficiency evaluation device according to claim 6, wherein in each small turbine outlet total pressure calculation module, the outlet total pressure P of the ith small turbine is determined out (i) Comprises the following steps:
Figure FDA0003799008770000031
9. a high-pressure turbine efficiency estimation apparatus based on heat-of-variation calculation as described in claim 6 wherein said residual based loop control module comprises:
a second high-pressure turbine efficiency initial guess value setting unit for setting a second high-pressure turbine efficiency initial guess value eta THtmp Setting multiple of the first high-pressure turbine efficiency initial guess value;
an internal circulation unit for initially guessing the value eta according to the efficiency of the second high-pressure turbine TH Determining the total outlet pressure of the last small turbine and assigning a second variable P 5outtmp
An outer loop unit for updating the first high pressure turbine efficiency initial guess value:
Figure FDA0003799008770000032
wherein eta is TH(new) The value is guessed for the new first high pressure turbine efficiency.
10. A high-pressure turbine efficiency evaluation device based on a specific heat change calculation as set forth in claim 9, wherein said set multiple is 1.01 times.
CN202210977690.0A 2022-08-15 2022-08-15 High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation Active CN115292946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210977690.0A CN115292946B (en) 2022-08-15 2022-08-15 High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210977690.0A CN115292946B (en) 2022-08-15 2022-08-15 High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation

Publications (2)

Publication Number Publication Date
CN115292946A true CN115292946A (en) 2022-11-04
CN115292946B CN115292946B (en) 2023-09-05

Family

ID=83830278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210977690.0A Active CN115292946B (en) 2022-08-15 2022-08-15 High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation

Country Status (1)

Country Link
CN (1) CN115292946B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116087613A (en) * 2023-04-07 2023-05-09 沃德传动(天津)股份有限公司 Reciprocating compressor energy efficiency calculation system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060196182A1 (en) * 2005-03-02 2006-09-07 Denso Corporation Controller for internal combustion engine with supercharger
CN102356222A (en) * 2009-01-22 2012-02-15 雷诺股份公司 Method and device for determining the pressure upstream from the turbine of supercharging turbocharger of thermal engine
US20150275753A1 (en) * 2012-09-13 2015-10-01 Avl List Gmbh Method for Monitoring at Least One Exhaust Gas Turbocharger
CN108931379A (en) * 2018-05-24 2018-12-04 北京理工大学 A kind of turbine efficiency measurement method
JP2019094854A (en) * 2017-11-24 2019-06-20 トヨタ自動車株式会社 Control device of internal combustion engine
CN109948231A (en) * 2019-03-14 2019-06-28 中国航发湖南动力机械研究所 Engine cycle parameters analysis method and device
US20200033229A1 (en) * 2017-02-16 2020-01-30 Hino Motors, Ltd. State quantity estimating device
CN112100862A (en) * 2020-09-24 2020-12-18 南京航空航天大学 Novel turbine isentropic efficiency fast solving method based on variable specific heat capacity
CN112507586A (en) * 2020-12-02 2021-03-16 中国航发沈阳发动机研究所 Rapid assessment method for two-dimensional temperature and strength of turbine air cooling blade

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060196182A1 (en) * 2005-03-02 2006-09-07 Denso Corporation Controller for internal combustion engine with supercharger
CN102356222A (en) * 2009-01-22 2012-02-15 雷诺股份公司 Method and device for determining the pressure upstream from the turbine of supercharging turbocharger of thermal engine
US20150275753A1 (en) * 2012-09-13 2015-10-01 Avl List Gmbh Method for Monitoring at Least One Exhaust Gas Turbocharger
US20200033229A1 (en) * 2017-02-16 2020-01-30 Hino Motors, Ltd. State quantity estimating device
JP2019094854A (en) * 2017-11-24 2019-06-20 トヨタ自動車株式会社 Control device of internal combustion engine
CN108931379A (en) * 2018-05-24 2018-12-04 北京理工大学 A kind of turbine efficiency measurement method
CN109948231A (en) * 2019-03-14 2019-06-28 中国航发湖南动力机械研究所 Engine cycle parameters analysis method and device
CN112100862A (en) * 2020-09-24 2020-12-18 南京航空航天大学 Novel turbine isentropic efficiency fast solving method based on variable specific heat capacity
CN112507586A (en) * 2020-12-02 2021-03-16 中国航发沈阳发动机研究所 Rapid assessment method for two-dimensional temperature and strength of turbine air cooling blade

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王华: "影响涡轮增压器涡轮端效率的几个因素", 《内燃机及配件》, no. 5, pages 9 - 10 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116087613A (en) * 2023-04-07 2023-05-09 沃德传动(天津)股份有限公司 Reciprocating compressor energy efficiency calculation system
CN116087613B (en) * 2023-04-07 2023-06-30 沃德传动(天津)股份有限公司 Reciprocating compressor energy efficiency calculation system

Also Published As

Publication number Publication date
CN115292946B (en) 2023-09-05

Similar Documents

Publication Publication Date Title
CN107783415B (en) Thermal power generating unit constant-sliding-pressure operation curve calculation method based on DCS data mining
US20170306970A1 (en) Extrapolation method of low rotational speed characteristic of compressor
CN115292946A (en) High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation
CN108108528B (en) One-dimensional matching design method for power turbine of split-shaft type aeroderivative
CN109373778B (en) Method for evaluating and checking cooling capacity of counter-flow natural ventilation cooling tower
CN111125931B (en) Automatic optimization type searching method and device for cooling tower
CN116186947A (en) Impeller simulation design method, device, equipment, medium and product
CN107203364B (en) Prediction and identification method for full-working-condition characteristics of gas compressor
CN112304623B (en) Effective thermal efficiency prediction method of marine diesel engine based on fuel components
JP5899685B2 (en) Intake air volume calculation method
CN112800588B (en) Simulation calculation method for air intake of cabin heat exchanger under transient driving working condition
CN112459984B (en) Performance test calculation method for isothermal compressor
CN111914389B (en) Cold end system operating point determining method based on curve fitting
CN107153408B (en) A kind of filtering method of fired power generating unit DCS data
CN116933425B (en) Performance prediction method and device for centrifugal compressor
CN116882189A (en) Correction method for heat exchange calculation model of air lubricating oil radiator of aeroengine
WO2009016210A1 (en) Method for designing turbo-machine blades, specifically gas turbine and steam turbine blades
CN215374069U (en) Flow resistance test platform of closed pure water cooling equipment
CN115688554A (en) Low-rotating-speed characteristic extrapolation and correction method for rotating part of turboshaft engine
CN114021354B (en) Construction method and construction device of fan coil heat exchange quantity calculation model
CN112949166B (en) Gas turbine steady-state model coefficient determination method
CN114036875A (en) Method and system for acquiring thermophysical parameters of supercritical fluid working medium
CN107103159B (en) Cooling water flow soft measurement method of steam turbine device
CN117077518A (en) Hydrogen mixed gas turbine parameter prediction method and device, electronic equipment and storage medium
CN106126803A (en) Refrigeration system analogy method and device

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