CN115292946B - 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 PDFInfo
- Publication number
- CN115292946B CN115292946B CN202210977690.0A CN202210977690A CN115292946B CN 115292946 B CN115292946 B CN 115292946B CN 202210977690 A CN202210977690 A CN 202210977690A CN 115292946 B CN115292946 B CN 115292946B
- Authority
- CN
- China
- Prior art keywords
- pressure
- total
- turbine
- pressure turbine
- outlet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/08—Thermal 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 indexes of the small turbines; 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 total outlet pressure of the last small turbine relative to the total outlet pressure of the high-pressure turbine, and if the change rate is larger than a preset value, updating the first high-pressure turbine efficiency initial guess value until the change rate is smaller than the preset value. The application realizes the calculation of the variable specific heat ratio of the high-pressure turbine efficiency, and greatly improves the calculation precision of the high-pressure turbine efficiency.
Description
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 variable specific heat calculation.
Background
In order to evaluate the efficiency of the components such as the fan, the air compressor, the turbine and the like of the aero-engine, the total temperature and the total pressure of the inlet and the outlet of the corresponding components are obtained in a measuring or calculating mode, then a constant engine adiabatic index kg (1.33 for fuel gas and 1.4 for air) is given, and the efficiency of the corresponding components is calculated according to the engine principle.
The efficiency of corresponding parts of the engine can be estimated and calculated preliminarily by adopting the method, the calculation accuracy is acceptable for the small temperature change range of the inlet and outlet of fans, compressors and the like, but the calculation accuracy is acceptable for the high-pressure turbine part, the inlet and outlet temperature change is large (from thousands of hundreds of Ks to hundreds of Ks), and the calculation accuracy is slightly high and can be calculated by adopting the average value of the inlet and outlet heat insulation indexes kg or kg of the inlet heat insulation indexes of the high-pressure turbine due to the fact that the heat insulation indexes kg of the engine are related to the gas temperature and the gas-oil ratio, and the calculation accuracy is also greatly error.
Disclosure of Invention
In order to solve one of the problems, the application provides a high-pressure turbine efficiency evaluation method and a high-pressure turbine efficiency evaluation device based on variable specific heat calculation, so that the efficiency of a high-pressure turbine can be evaluated more accurately.
The first aspect of the application provides a high-pressure turbine efficiency evaluation method based on variable specific heat calculation, which mainly comprises the following steps:
s1, calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating the oil-gas ratio according to the air flow rate at the inlet of the high-pressure turbine and the fuel flow meter;
s2, dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the inlet total temperature T of each small turbine according to the inlet total temperature and the outlet total temperature of the high-pressure turbine and the equal interval cooling principle in (i) And the total temperature T of the outlet out (i)；
S3, according to the inlet total temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k of the ith small turbine g (i)；
Step S4, determining turbine efficiency eta of the ith small turbine TH(i) And based on the inlet total temperature T of the ith small turbine in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total pressure P of the ith small turbine out (i) Wherein the total inlet pressure of the first small turbine is the inlet pressure of the high-pressure turbineTotal pressure P of mouth 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-S4 until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out ；
Step S6, determining the total outlet pressure P of the last small turbine out (num) relative to the high pressure turbine outlet total pressure P 5 If the change rate is smaller than the preset value, taking the first high-pressure turbine efficiency initial guess value given in the step S1 as the final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value in the step S1, and repeating the steps until the change rate is smaller than the preset value.
Preferably, in step S2, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the outlet of the high-pressure turbine, num is the number of small turbines divided by the high-pressure turbine.
Preferably, in step S4, the outlet total pressure P of the ith small turbine is determined out (i) The method comprises the following steps:
preferably, in step S6, updating the first high-pressure turbine efficiency preliminary guess includes:
step S61, setting a first guess value eta of the efficiency of the second high-pressure turbine THtmp Setting a multiple of the first high pressure turbine efficiency initial guess;
step S62, according to the first guess value eta of the efficiency of the second high-pressure turbine TH Performing the inner circulation of the steps S2-S5, determining the total outlet pressure of the last small turbine, and assigning the total outlet pressure to the second variable P 5outtmp ；
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing the outer loop of the steps S2-S6:
wherein eta TH(new) A new first high pressure turbine efficiency initial guess is made.
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 variable specific heat calculation, mainly comprising:
the efficiency and 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 the oil-gas ratio according to the air flow rate at the inlet of the high-pressure turbine and the fuel flow meter;
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 inlet total temperature T of each small turbine according to the inlet total temperature and the outlet total temperature of the high-pressure turbine and the equal interval cooling principle in (i) And the total temperature T of the outlet 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 of the ith small turbine g (i)；
A total pressure calculation module of each small turbine outlet for determining turbine efficiency eta of the ith small turbine TH(i) And based on the inlet total temperature T of the ith small turbine in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total 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 final small turbine outlet total pressure calculation module is used for repeatedly calling the modules until the outlet of the final small turbine is calculatedTotal pressure P of mouth out (num) and assigned to the first variable P 5out ；
Residual-based cycle control module for determining the total outlet pressure P of the last small turbine out (num) relative to the high pressure turbine outlet total pressure P 5 If the change rate is smaller than a preset value, taking the given first high-pressure turbine efficiency initial guess value as a final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value until the change rate is smaller than the preset value.
Preferably, in the high-pressure turbine dividing module, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the outlet of the high-pressure turbine, num is the number of small turbines divided by the high-pressure turbine.
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) The method comprises the following steps:
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 a multiple of the first high pressure turbine efficiency initial guess;
an internal circulation unit for estimating a first high-pressure turbine efficiency value eta based on the second high-pressure turbine efficiency value eta TH Determining the total outlet pressure of the last small turbine and assigning to a second variable P 5outtmp ；
An outer circulation unit for updating the first high pressure turbine efficiency initial guess:
wherein eta TH(new) A new first high pressure turbine efficiency initial guess is made.
Preferably, the set multiple is 1.01 times.
According to the application, the high-pressure turbine is decomposed into a plurality of small turbines for calculation according to a certain temperature drop interval, so that the calculation of the variable specific heat ratio of the high-pressure turbine efficiency is realized, and the calculation precision of the high-pressure turbine efficiency can be greatly improved.
Drawings
FIG. 1 is a flow chart of a preferred embodiment of the high pressure turbine efficiency assessment method based on variable specific heat calculation of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application become more apparent, the technical solutions in the embodiments of the present application will be described in more detail with reference to the accompanying 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 application. The embodiments described below by referring to the drawings are exemplary and intended to illustrate the present application and should not be construed as limiting the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present application without making any inventive effort, are intended to fall within the scope of the present application. Embodiments of the present application will be described in detail below with reference to the accompanying drawings.
In general, in order to obtain higher-precision high-pressure turbine efficiency, a variable specific heat calculation formula (1) is adopted to calculate an adiabatic index kg according to the gas-oil ratio and the total temperature of gas, and then the efficiency of the high-pressure turbine can be calculated according to a formula (2) according to the total inlet-outlet temperature and the total pressure of the high-pressure turbine obtained by a measurement or calculation mode.
kg=f (oil/gas ratio, total temperature) … … … … … … … … … … … … … … (1)
Wherein:
t4 and P4 are the total inlet temperature and total inlet pressure of the high-pressure turbine;
t5 and P5 are the total temperature and total pressure of the outlet of the high-pressure turbine;
kg is the adiabatic index.
To calculate the high pressure turbine efficiency more accurately, the high pressure turbine calculation process may be divided into numerous "small turbines" for calculation. And (3) dividing each small turbine according to fixed temperature drop, sequentially calculating the total temperature and total pressure of the outlets of each small turbine according to the initial value of the efficiency of the high-pressure turbine calculated by the formula (2), and finally constructing residual errors according to the outlet pressure of the high-pressure turbine, and updating and iterating the initial value of the efficiency of the high-pressure turbine until the residual errors meet the error requirement. Specifically as shown in fig. 1, the first aspect of the present application provides a high-pressure turbine efficiency evaluation method based on variable specific heat calculation, which mainly includes:
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 inlet air flow of the high-pressure turbine and the fuel flow meter.
S2, dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the inlet total temperature T of each small turbine according to the inlet total temperature and the outlet total temperature of the high-pressure turbine and the equal interval cooling principle in (i) And the total temperature T of the outlet out (i)。
In some alternative embodiments, in step S2, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the high-pressure turbine outlet, num is the high-pressure turbineThe number of small turbines divided.
S3, according to the inlet total temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k of the ith small turbine g (i)。
Step S4, determining turbine efficiency eta of the ith small turbine TH(i) And based on the inlet total temperature T of the ith small turbine in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total 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.
In some alternative embodiments, in step S4, the outlet total pressure P of the ith small turbine is determined out (i) The method comprises the following steps:
step S5, repeating the steps S2-S4 until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out The method comprises the steps of carrying out a first treatment on the surface of the Namely P 5out ＝P out (num)。
Step S6, determining the total outlet pressure P of the last small turbine out (num) relative to the high pressure turbine outlet total pressure P 5 If the change rate is smaller than the preset value, taking the first high-pressure turbine efficiency initial guess value given in the step S1 as the final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value in the step S1, and repeating the steps until the change rate is smaller than the preset value.
In some alternative embodiments, in step S6, updating the first high pressure turbine efficiency preliminary guess includes:
step S61, setting a first guess value eta of the efficiency of the second high-pressure turbine THtmp Setting a multiple of the first high pressure turbine efficiency initial guess;
step S62, according to the first guess value eta of the efficiency of the second high-pressure turbine TH PerformingThe internal circulation of the steps S2-S5, the total outlet pressure of the last small turbine is determined and assigned to the second variable P 5outtmp ；
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing the outer loop of the steps S2-S6:
wherein eta TH(new) A new first high pressure turbine efficiency initial guess is made.
In some alternative embodiments, in step S61, the set multiple is 1.01 times.
In this embodiment, compare P 5out And P 5 Errors of (e.g.)Terminating the calculation; such asSetting eta THtmp ＝η TH *1.01, the total outlet pressure P of the last small turbine is carried out again out (num) calculation and assignment to P 5outtmp The method comprises the steps of carrying out a first treatment on the surface of the Thereafter, press->Updating eta TH . Finally repeating the process until a new eta 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 variable specific heat calculation corresponding to the above method, mainly comprising:
the efficiency and 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 the oil-gas ratio according to the air flow rate at the inlet of the high-pressure turbine and the fuel flow meter;
a high-pressure turbine dividing module for dividing the high-pressure turbine from an inlet to an outlet into a plurality of small partsThe turbine, according to the inlet total temperature and outlet total temperature of the high-pressure turbine, determining the inlet total temperature T of each small turbine according to the principle of equal interval cooling in (i) And the total temperature T of the outlet 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 of the ith small turbine g (i)；
A total pressure calculation module of each small turbine outlet for determining turbine efficiency eta of the ith small turbine TH(i) And based on the inlet total temperature T of the ith small turbine in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total 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 final small turbine outlet total pressure calculation module is used for repeatedly adjusting the modules until the outlet total pressure P of the final small turbine is calculated out (num) and assigned to the first variable P 5out ；
Residual-based cycle control module for determining the total outlet pressure P of the last small turbine out (num) relative to the high pressure turbine outlet total pressure P 5 If the change rate is smaller than a preset value, taking the given first high-pressure turbine efficiency initial guess value as a final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value until the change rate is smaller than the preset value.
In some alternative embodiments, in the high-pressure turbine partition module, an inlet total temperature T of an ith small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the high-pressure turbine outlet, num is the small division of the high-pressure turbineThe number of turbines.
In some alternative embodiments, in the calculating module of the outlet total pressure of each small turbine, the outlet total pressure P of the ith small turbine is determined out (i) The method comprises the following steps:
in some alternative embodiments, 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 a multiple of the first high pressure turbine efficiency initial guess;
an internal circulation unit for estimating a first high-pressure turbine efficiency value eta based on the second high-pressure turbine efficiency value eta TH Determining the total outlet pressure of the last small turbine and assigning to a second variable P 5outtmp ；
An outer circulation unit for updating the first high pressure turbine efficiency initial guess:
wherein eta TH(new) A new first high pressure turbine efficiency initial guess is made.
In some alternative embodiments, the set factor is 1.01.
According to the application, the high-pressure turbine is decomposed into a plurality of small turbines for calculation according to a certain temperature drop interval, so that the calculation of the variable specific heat ratio of the high-pressure turbine efficiency is realized, and the calculation precision of the high-pressure turbine efficiency can be greatly improved.
While the application has been described in detail in the foregoing general description and with reference to specific embodiments thereof, it will be apparent to one skilled in the art that modifications and improvements can be made thereto. Accordingly, such modifications or improvements may be made without departing from the spirit of the application and are intended to be within the scope of the application as claimed.
Claims (6)
1. A high pressure turbine efficiency evaluation method based on variable specific heat calculation, characterized by comprising:
s1, calculating a first high-pressure turbine efficiency initial guess value eta according to a high-pressure turbine efficiency formula TH Calculating the oil-gas ratio according to the inlet air flow of the high-pressure turbine and the fuel flow meter, and obtaining the total inlet temperature T of the high-pressure turbine 4 Total temperature T of high-pressure turbine outlet 5 Total pressure P at inlet of high-pressure turbine 4 Total pressure P at high-pressure turbine outlet 5 ；
S2, dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet, and determining the inlet total temperature T of each small turbine according to the inlet total temperature and the outlet total temperature of the high-pressure turbine and the equal interval cooling principle in (i) And the total temperature T of the outlet out (i)；
S3, according to the inlet total temperature T of the ith small turbine in (i) And said gas-oil ratio determining the adiabatic index k of the ith small turbine g (i)；
Step S4, determining turbine efficiency eta of the ith small turbine TH(i) And based on the adiabatic index k of the ith small turbine g (i) Total inlet temperature T in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total 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, wherein the total outlet pressure P of the ith small turbine is determined out (i) The method comprises the following steps:
step S5, repeating the steps S2-S4 until the outlet total pressure P of the last small turbine is calculated out (num) and assigned to the first variable P 5out ；
Step S6, determining the total outlet pressure P of the last small turbine out (num) relative to high pressureTurbine outlet total pressure P 5 If the change rate is smaller than the preset value, taking the first high-pressure turbine efficiency initial guess value given in the step S1 as the final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value in the step S1, and repeating the steps until the change rate is smaller than the preset value;
wherein, in step S6, updating the first high-pressure turbine efficiency initial guess value includes:
step S61, setting a first guess value eta of the efficiency of the second high-pressure turbine THtmp Setting a multiple of the first high pressure turbine efficiency initial guess;
step S62, according to the first guess value eta of the efficiency of the second high-pressure turbine THtmp Performing the inner circulation of the steps S2-S5, determining the total outlet pressure of the last small turbine, and assigning the total outlet pressure to the second variable P 5outtmp ；
Step S63, updating the first high-pressure turbine efficiency initial guess value, and executing the outer loop of the steps S2-S6:
wherein eta TH(new) For a new first high pressure turbine efficiency preliminary guess, for returning to step S1, replacing the old first high pressure turbine efficiency preliminary guess η TH 。
2. The high-pressure turbine efficiency evaluation method based on variable specific heat calculation as recited in claim 1, wherein in step S2, the inlet total temperature T of the ith small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the outlet of the high-pressure turbine, num is the number of small turbines divided by the high-pressure turbine.
3. The high-pressure turbine efficiency evaluation method based on the variable specific heat calculation according to claim 1, wherein in step S61, the set multiple is 1.01 times.
4. A high-pressure turbine efficiency evaluation device based on variable specific heat calculation, characterized by comprising:
the efficiency and 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 the oil-gas ratio according to the inlet air flow of the high-pressure turbine and the fuel flow meter, and obtaining the total inlet temperature T of the high-pressure turbine 4 Total temperature T of high-pressure turbine outlet 5 Total pressure P at inlet of high-pressure turbine 4 Total pressure P at high-pressure turbine outlet 5 ；
A high-pressure turbine dividing module for dividing the high-pressure turbine into a plurality of small turbines from an inlet to an outlet according to the adiabatic index k of the high-pressure turbine g (i) The total inlet temperature and the total outlet temperature are determined according to the principle of equal interval cooling in (i) And the total temperature T of the outlet 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 of the ith small turbine g (i)；
A total pressure calculation module of each small turbine outlet for determining turbine efficiency eta of the ith small turbine TH(i) And based on the inlet total temperature T of the ith small turbine in (i) Total outlet temperature T out (i) Total inlet pressure P in (i) Determining the outlet total 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, wherein the total outlet pressure P of the ith small turbine is determined out (i) The method comprises the following steps:
the final small turbine outlet total pressure calculation module is used for repeatedly adjusting the modules until the outlet total pressure P of the final small turbine is calculated out (num) and assigned to the first variable P 5out ；
Residual-based cycle control module for determining the total outlet pressure P of the last small turbine out (num) relative to the high pressure turbine outlet total pressure P 5 If the change rate is smaller than a preset value, taking a given first high-pressure turbine efficiency initial guess value as a final value of the high-pressure turbine efficiency, otherwise, updating the first high-pressure turbine efficiency initial guess value until the change rate is smaller than the preset value;
wherein, 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 a multiple of the first high pressure turbine efficiency initial guess;
an internal circulation unit for estimating a first high-pressure turbine efficiency value eta based on the second high-pressure turbine efficiency value eta THtmp Determining the total outlet pressure of the last small turbine and assigning to a second variable P 5outtmp ；
An outer circulation unit for updating the first high pressure turbine efficiency initial guess:
wherein eta TH(new) For new first high pressure turbine efficiency initial guess value, for return efficiency and gas ratio calculation module, replace old first high pressure turbine efficiency initial guess value eta TH 。
5. The high-pressure turbine efficiency evaluation device based on variable specific heat calculation as claimed in claim 4, wherein in the high-pressure turbine dividing module, an inlet total temperature T of an i-th small turbine in (i) The total outlet temperature T of the ith small turbine is the total outlet temperature of the ith-1 th small turbine out (i) The method comprises the following steps:
T out (i)＝T 4 -ΔT*i；
wherein T is 4 For the total inlet temperature of the high-pressure turbine, Δt= (T 5 -T 4 ) /num, wherein T 5 For the total temperature of the outlet of the high-pressure turbine, num is the number of small turbines divided by the high-pressure turbine.
6. The high-pressure turbine efficiency evaluation device based on variable specific heat calculation according to claim 4, wherein the set multiple is 1.01 times.
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 CN115292946A (en) | 2022-11-04 |
| CN115292946B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116087613B (en) * | 2023-04-07 | 2023-06-30 | 沃德传动（天津）股份有限公司 | Reciprocating compressor energy efficiency calculation system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
| 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100749620B1 (en) * | 2005-03-02 | 2007-08-14 | 가부시키가이샤 덴소 | Controller for internal combustion engine with supercharger |
| AT513137B1 (en) * | 2012-09-13 | 2014-02-15 | Avl List Gmbh | Method for monitoring at least one exhaust gas turbocharger |
| JP6751680B2 (en) * | 2017-02-16 | 2020-09-09 | 日野自動車株式会社 | State quantity estimation device |
-
2022
- 2022-08-15 CN CN202210977690.0A patent/CN115292946B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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)
| Title |
|---|
| 影响涡轮增压器涡轮端效率的几个因素;王华;《内燃机及配件》(第5期);第9-10页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115292946A (en) | 2022-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115292946B (en) | High-pressure turbine efficiency evaluation method and device based on variable specific heat calculation | |
| EP1882215B1 (en) | Analysis method | |
| CN108108528B (en) | One-dimensional matching design method for power turbine of split-shaft type aeroderivative | |
| US20100042367A1 (en) | Analysis method | |
| JP2012031799A (en) | Thermal efficiency analyzing method, performance evaluation method, thermal efficiency analysis program, and performance evaluation program of geothermal power plant | |
| CN110489706B (en) | Simplified calculation method for energy efficiency index EEI of plate heat exchanger | |
| Chiu et al. | On the steady-state performance of the Poisson double GWMA control chart | |
| CN107870611A (en) | Workshop analogue means and workshop analogy method | |
| 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 | |
| Herrera et al. | A comparative analysis of turbine rotor inlet temperature models | |
| CN114861810A (en) | Coal gasification device process diagnosis method and system | |
| CN111275320B (en) | Performance adjustment data processing method, system and storage medium of generator set | |
| CN110889217B (en) | Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine | |
| CN113255207A (en) | Iterative multi-output-Markov chain-based multi-step prediction method for argon fraction variable of air separation system | |
| CN112670997A (en) | Electric heating energy source system time sequence probability load flow calculation method considering photovoltaic uncertainty | |
| CN107153408B (en) | A kind of filtering method of fired power generating unit DCS data | |
| CN114491809A (en) | Method for evaluating combustion efficiency of aero-engine | |
| WO2009016210A1 (en) | Method for designing turbo-machine blades, specifically gas turbine and steam turbine blades | |
| CN112949166B (en) | Gas turbine steady-state model coefficient determination method | |
| CN116882189A (en) | Correction method for heat exchange calculation model of air lubricating oil radiator of aeroengine | |
| CN114896787A (en) | Method suitable for calculating internal trajectory of engine at any temperature | |
| CN106126803A (en) | Refrigeration system analogy method and device | |
| Schneider | Transient Thermodynamic Simulation of the DLR Turbine Test Facility NG-Turb |
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 |