CN110889217B - Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine - Google Patents

Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine Download PDF

Info

Publication number
CN110889217B
CN110889217B CN201911142671.0A CN201911142671A CN110889217B CN 110889217 B CN110889217 B CN 110889217B CN 201911142671 A CN201911142671 A CN 201911142671A CN 110889217 B CN110889217 B CN 110889217B
Authority
CN
China
Prior art keywords
combustion chamber
temperature
fuel
outlet
bout
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
CN201911142671.0A
Other languages
Chinese (zh)
Other versions
CN110889217A (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.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
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 Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN201911142671.0A priority Critical patent/CN110889217B/en
Publication of CN110889217A publication Critical patent/CN110889217A/en
Application granted granted Critical
Publication of CN110889217B publication Critical patent/CN110889217B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a correction method of outlet temperature of a one-dimensional simulation model of a combustion chamber of a gas turbine. The method comprises the following steps: according to the principle of conservation of mass and conservation of energy, a pressure dynamic model and a temperature dynamic model of the combustion chamber are obtained; step two: setting the temperature of fuel and determining heat Q; the heat Q is the heat released by the combustion of the fuel minus the enthalpy difference of the outlet working medium between the outlet temperature of the combustion chamber and the ambient temperature. The invention has the advantage of enabling the calculation of the outlet temperature of the combustion chamber to be more accurate.

Description

Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine
Technical Field
The invention relates to the technical field of simulation models of combustion chambers of gas turbines, in particular to a method for correcting outlet temperature of a one-dimensional simulation model of a combustion chamber of a gas turbine.
Background
The combustor is an important component of the gas turbine, and the combustor outlet temperature determines the temperature ratio of the gas turbine, which has an important effect on the performance of the gas turbine. Definition of fuel calorific value: the unit mass or unit volume of fuel is completely burned and the combustion products cool to a pre-combustion temperature and release heat.
The outlet temperature in the combustion chamber model is calculated according to the principle of conservation of energy, air enters the combustion chamber and is mixed with fuel for combustion, chemical energy is released in the combustion process, and the temperature of residual air which does not participate in combustion and combustion products for absorbing the chemical energy rises.
The chemical energy of the combustion chamber is replaced by the heat value of the fuel in the previous model, but the calculated outlet temperature of the combustion chamber is higher than the experimental value, so that the simulation calculation of the outlet temperature of the combustion chamber is inaccurate.
Therefore, there is a need to develop a method for correcting the outlet temperature of a one-dimensional simulation model of a combustion chamber of a gas turbine that improves the accuracy of the calculation of the outlet temperature of the combustion chamber.
Disclosure of Invention
The invention aims to provide a correction method of the outlet temperature of a one-dimensional simulation model of a combustion chamber of a gas turbine, so that the calculation of the outlet temperature of the combustion chamber is more accurate.
In order to achieve the above purpose, the technical scheme of the invention is as follows: the method for correcting the outlet temperature of the one-dimensional simulation model of the combustion chamber of the gas turbine is characterized by comprising the following steps of: comprises the following steps of the method,
step one: according to the principle of conservation of mass and conservation of energy, a pressure dynamic model and a temperature dynamic model of the combustion chamber are obtained;
step two: setting the temperature of fuel and determining heat Q; the heat Q is the fuel calorific value minus the enthalpy difference of the outlet working medium between the outlet temperature of the combustion chamber and the ambient temperature.
In the above technical scheme, in the first step, according to the principles of conservation of mass and conservation of energy, a pressure dynamic model and a temperature dynamic model of a combustion chamber are obtained;
wherein, the pressure dynamic model is shown as a formula (1-1):
Figure BDA0002281376720000021
in the formula (1-1): p (P) bout Is the combustion chamber outlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is the gas constant of fuel gas, kg/s;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
v is the volume of the combustion chamber, m 3 ；
t is time, S;
the temperature dynamic model is shown in the formula (1-2):
Figure BDA0002281376720000022
in the formula (1-2): p (P) bout Is the combustion chamber outlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is a gas constant of fuel gas;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
k is an adiabatic constant;
h bin j/kg is the enthalpy value of the inlet working medium of the combustion chamber;
h bout j/kg is the enthalpy value of the working medium at the outlet of the combustion chamber;
v is the volume of the combustion chamber, m 3 ；
C pg The specific heat capacity of the fuel gas is fixed, J/(kg.K);
q is the heat released by the fuel, J;
t is time, s.
In the above technical scheme, in the second step, the fuel chemical formula is C x H y O z S v Heat value of H u ，C x H y O z S v The chemical reaction formula of the combustion reaction is shown in the formula (1-3):
Figure BDA0002281376720000031
setting the fuel temperature to T 0 Determining the heat quantity Q as represented by the formula (1-4) according to the chemical expression:
Figure BDA0002281376720000032
in the formulae (1-3) and (1-4):
Figure BDA0002281376720000033
is CO 2 Relative molecular mass;
Figure BDA0002281376720000034
is H 2 O relative molecular mass;
Figure BDA0002281376720000035
is SO 2 Relative molecular mass;
Figure BDA0002281376720000036
is C 8 H 16 Relative molecular mass;
Figure BDA0002281376720000037
is CO 2 At T 0 Enthalpy at temperature, J/kg;
Figure BDA0002281376720000038
is H 2 O is at T 0 Enthalpy at temperature, J/kg;
Figure BDA0002281376720000039
is SO 2 At T 0 Enthalpy at temperature, J/kg;
Figure BDA00022813767200000310
is CO 2 At a temperature of T bout Enthalpy value at time, J/kg;
Figure BDA0002281376720000041
is H 2 O has a temperature T bout Enthalpy value at time, J/kg;
Figure BDA0002281376720000042
is SO 2 At a temperature of T bout Enthalpy value at time, J/kg;
η b is combustion chamber efficiency;
x, y, z, v represent the molecular weight of each element in the fuel expression;
G f kg/s is the fuel flow;
H u is the heating value of the fuel, J/kg.
The invention has the following advantages:
(1) The invention considers that the outlet temperature of the combustion chamber of the gas turbine is obviously higher than the temperature of the fuel before combustion, so that the enthalpy difference of the combustion product at the outlet temperature and the temperature before combustion needs to be subtracted from the heat released by the fuel in order to calculate the outlet temperature of the combustion chamber more accurately; in order to accurately reflect the temperature distribution condition of the outlet of the combustion chamber of the gas turbine, the method combines the traditional calculation method taking the calorific value of fuel as the calorific value, considers partial enthalpy which is not released by the outlet working medium of the combustion chamber due to high temperature in the process of calculating the calorific value released by the combustion of the fuel, and removes the enthalpy difference between the outlet temperature of the combustion chamber and the environmental temperature (25 ℃) in the calculation of the outlet temperature, thereby enabling the calculation of the outlet temperature of the combustion chamber to be more accurate;
(2) According to the invention, the enthalpy difference between the outlet temperature and the inlet temperature of the combustion product is subtracted on the basis of the fuel heat value by correcting the fuel heat value, and the enthalpy difference is taken as the fuel heat release amount, so that the simulation calculation of the outlet temperature of the combustion chamber is more accurate; the problem that the temperature of the simulation outlet of the combustion chamber is higher than the measured data, which is caused by the fact that the heat value of the fuel is taken as the heat release amount of the fuel to perform simulation calculation in the prior art, is solved.
Drawings
FIG. 1 is a graph comparing turbine outlet temperature with measured experimental values in two simulation modes.
As shown in fig. 1, it can be seen that, under the same fuel injection quantity under the same working condition, the maximum temperature error after the turbine is 10 ℃ and the maximum relative error after the turbine is corrected according to the method of the invention (i.e. the simulation data curve after the improvement is shown in fig. 1) is 2 ℃ by using the fuel calorific value as a combustion chamber calculation model of the calorific value (i.e. the simulation data curve before the improvement is shown in fig. 1); the analysis shows that the temperature after the turbine is calculated by the method provided by the invention is closer to a true value, which proves that the simulation model of the combustion chamber is more accurate, so that the calculation of the outlet temperature of the combustion chamber is more accurate.
Detailed Description
The following detailed description of the invention is, therefore, not to be taken in a limiting sense, but is made merely by way of example. While making the advantages of the present invention clearer and more readily understood by way of illustration.
As can be seen with reference to the accompanying drawings: the method for correcting the outlet temperature of the one-dimensional simulation model of the combustion chamber of the gas turbine comprises the following steps,
step one: according to the principle of conservation of mass and conservation of energy, a pressure dynamic model and a temperature dynamic model of the combustion chamber are obtained;
step two: setting the temperature of fuel and determining heat Q; the heat Q is the fuel heating value minus the enthalpy difference of the outlet working medium between the outlet temperature of the combustion chamber and the ambient temperature.
Further, in the first step, a pressure dynamic model and a temperature dynamic model of the combustion chamber are obtained according to the principles of conservation of mass and conservation of energy;
wherein, the pressure dynamic model is shown as a formula (1-1):
Figure BDA0002281376720000051
in the formula (1-1): p (P) bout Is the combustion chamber outlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is a gas constant of fuel gas;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
v is the volume of the combustion chamber, m 3 ；
t is time, S;
the temperature dynamic model is shown in the formula (1-2):
Figure BDA0002281376720000061
in the formula (1-2): p (P) bout Is a combustion chamberOutlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is a gas constant of fuel gas;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
k is an adiabatic constant;
h bin j/kg is the enthalpy value of the inlet working medium of the combustion chamber;
h bout j/kg is the enthalpy value of the working medium at the outlet of the combustion chamber;
v is the volume of the combustion chamber, m 3 ；
C pg The specific heat capacity of the fuel gas is fixed, J/(kg.K);
q is the heat released by the fuel, J;
t is time, s;
when the components of the working medium at the inlet and the outlet of the combustion chamber are the same, the factors for determining the steady-state value of the temperature at the outlet of the combustion chamber only comprise the heat release quantity of the fuel; therefore, whether the fuel exotherm is accurate directly determines the accuracy of the combustion chamber outlet temperature calculation.
Further, in step two, considering that the outlet temperature of the combustion chamber of the gas turbine is obviously higher than the temperature of the fuel before combustion, in order to calculate the outlet temperature of the combustion chamber more accurately, the enthalpy difference between the outlet temperature and the temperature before combustion of the combustion product needs to be subtracted from the heat released by the fuel; the chemical formula of the fuel is C x H y O z S v Heat value of H u ，C x H y O z S v The chemical reaction formula of the combustion reaction is shown in the formula (1-3):
Figure BDA0002281376720000071
setting the fuel temperature to T 0 Determining the heat quantity Q as represented by the formula (1-4) according to the chemical expression:
Figure BDA0002281376720000072
in the formulae (1-3) and (1-4):
Figure BDA0002281376720000073
is CO 2 Relative molecular mass;
Figure BDA0002281376720000074
is H 2 O relative molecular mass;
Figure BDA0002281376720000075
is SO 2 Relative molecular mass;
Figure BDA0002281376720000076
is C 8 H 16 Relative molecular mass;
Figure BDA0002281376720000077
is CO 2 At T 0 Enthalpy at temperature, J/kg;
Figure BDA0002281376720000078
is H 2 O is at T 0 Enthalpy at temperature, J/kg;
Figure BDA0002281376720000079
is SO 2 At T 0 Enthalpy at temperature, J/kg;
Figure BDA00022813767200000710
is CO 2 At a temperature of T bout Enthalpy value at time, J/kg;
Figure BDA00022813767200000711
is H 2 O has a temperature T bout Enthalpy value at time, J/kg;
Figure BDA00022813767200000712
is SO 2 At a temperature of T bout Enthalpy value at time, J/kg;
η b is combustion chamber efficiency;
x, y, z, v represent the molecular weight of each element in the fuel expression;
G f kg/s is the fuel flow;
H u is the heating value of the fuel, J/kg.
In particular, since the outlet temperature of the combustion chamber is high, direct measurement is difficult, and the turbine post-temperature is commonly used for replacing the outlet temperature of the combustion chamber in actual calculation (industry practice) for quantitatively researching the change rule of the outlet temperature of the combustion chamber.
Other non-illustrated parts are known in the art.

Claims (2)

1. The method for correcting the outlet temperature of the one-dimensional simulation model of the combustion chamber of the gas turbine is characterized by comprising the following steps of: comprises the following steps of the method,
step one: according to the principle of conservation of mass and conservation of energy, a pressure dynamic model and a temperature dynamic model of the combustion chamber are obtained;
step two: setting the temperature of fuel and determining heat Q; the heat Q is the enthalpy difference of the fuel calorific value minus the outlet working medium between the outlet temperature of the combustion chamber and the ambient temperature;
in the second step, the fuel has a chemical formula of C x H y O z S v Heat value of H u ，C x H y O z S v The chemical reaction formula of the combustion reaction is shown in the formula (1-3):
Figure FDA0004135809180000011
setting the fuel temperature to T 0 The heat quantity Q is determined according to the chemical expression and is shown as the formula (1-4):
Figure FDA0004135809180000012
in the formulae (1-3) and (1-4):
Figure FDA0004135809180000013
is CO 2 Relative molecular mass;
Figure FDA0004135809180000014
is H 2 O relative molecular mass;
Figure FDA0004135809180000015
is SO 2 Relative molecular mass;
Figure FDA0004135809180000016
is C 8 H 16 Relative molecular mass;
Figure FDA0004135809180000017
is CO 2 At T 0 Enthalpy at temperature, J/kg;
Figure FDA0004135809180000018
is H 2 O is at T 0 Enthalpy at temperature, J/kg;
Figure FDA0004135809180000019
is SO 2 At T 0 Enthalpy at temperature, J/kg;
Figure FDA00041358091800000110
is CO 2 At a temperature of T bout Enthalpy value at time, J/kg;
Figure FDA00041358091800000111
is H 2 O has a temperature T bout Enthalpy value at time, J/kg;
Figure FDA0004135809180000021
is SO 2 At a temperature of T bout Enthalpy value at time, J/kg;
η b is combustion chamber efficiency;
x, y, z, v represent the molecular weight of each element in the fuel expression;
G f kg/s is the fuel flow;
H u is the heating value of the fuel, J/kg.
2. The method for correcting the outlet temperature of a one-dimensional simulation model of a combustion chamber of a gas turbine according to claim 1, wherein: in the first step, a pressure dynamic model and a temperature dynamic model of a combustion chamber are obtained according to the principles of conservation of mass and conservation of energy;
wherein the dynamic pressure model of the combustion chamber is shown in formula (1-1):
Figure FDA0004135809180000022
in the formula (1-1): p (P) bout Is the combustion chamber outlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is a gas constant of fuel gas;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
v is the volume of the combustion chamber, m 3 ；
t is time, S;
the temperature dynamic model is shown in the formula (1-2):
Figure FDA0004135809180000023
in the formula (1-2): p (P) bout Is the combustion chamber outlet pressure, pa;
T bout is the combustion chamber outlet temperature, T;
R g is a gas constant of fuel gas;
G f is fuel flow, kg/s;
G bin kg/s for combustor inlet flow;
G bout kg/s for combustor outlet flow;
k is an adiabatic constant;
h bin j/kg is the enthalpy value of the inlet working medium of the combustion chamber; h is a bout J/kg is the enthalpy value of the working medium at the outlet of the combustion chamber; v is the volume of the combustion chamber, m 3 ；
C pg The specific heat capacity of the fuel gas is fixed, J/(kg.K); q is the heat released by the fuel, J;
t is time, s.
CN201911142671.0A 2019-11-20 2019-11-20 Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine Active CN110889217B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911142671.0A CN110889217B (en) 2019-11-20 2019-11-20 Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911142671.0A CN110889217B (en) 2019-11-20 2019-11-20 Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine

Publications (2)

Publication Number Publication Date
CN110889217A CN110889217A (en) 2020-03-17
CN110889217B true CN110889217B (en) 2023-06-27

Family

ID=69748061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911142671.0A Active CN110889217B (en) 2019-11-20 2019-11-20 Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine

Country Status (1)

Country Link
CN (1) CN110889217B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112484075B (en) * 2020-12-01 2023-04-21 中国航发沈阳发动机研究所 Combustion chamber outlet temperature field correction method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008163746A (en) * 2006-12-26 2008-07-17 Tokyo Electric Power Co Inc:The Method and device for calculating gas turbine inlet temperature
CN104896457A (en) * 2015-06-03 2015-09-09 华北电力科学研究院有限责任公司 All-level boiler heating surface flue-gas temperature determination method and device
CN106649945A (en) * 2016-09-30 2017-05-10 赵志渊 Calculation method for determining generated power and smoke discharge parameters of single-shaft gas turbine
CN108843451A (en) * 2018-05-31 2018-11-20 中国航发沈阳发动机研究所 Gas-turbine combustion chamber outlet temperature calculation method
CN110268148A (en) * 2017-02-23 2019-09-20 三菱日立电力系统株式会社 Gas turbine control device, gas turbine complete equipment and gas turbine control method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140114619A1 (en) * 2012-10-23 2014-04-24 Tenneco Automotive Operating Company Inc. Burner Outlet Designs for Locomotive Burner Integration

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008163746A (en) * 2006-12-26 2008-07-17 Tokyo Electric Power Co Inc:The Method and device for calculating gas turbine inlet temperature
CN104896457A (en) * 2015-06-03 2015-09-09 华北电力科学研究院有限责任公司 All-level boiler heating surface flue-gas temperature determination method and device
CN106649945A (en) * 2016-09-30 2017-05-10 赵志渊 Calculation method for determining generated power and smoke discharge parameters of single-shaft gas turbine
CN110268148A (en) * 2017-02-23 2019-09-20 三菱日立电力系统株式会社 Gas turbine control device, gas turbine complete equipment and gas turbine control method
CN108843451A (en) * 2018-05-31 2018-11-20 中国航发沈阳发动机研究所 Gas-turbine combustion chamber outlet temperature calculation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钟华贵等.《燃气分析测量高温燃气温度的方法》.《燃气分析测量高温燃气温度的方法》.2005,第460-466. *

Also Published As

Publication number Publication date
CN110889217A (en) 2020-03-17

Similar Documents

Publication Publication Date Title
CA2858631C (en) Control of the gas composition in a gas turbine power plant with flue gas recirculation
CN108843451B (en) Gas-turbine combustion chamber outlet temperature calculation method
CN103080560B (en) Method for operating gas compressor, and gas turbine provided with gas compressor
JP5837111B2 (en) Method for adjusting natural gas temperature in a fuel supply line of a gas turbine engine
US11592430B2 (en) Method for estimating a combustion characteristic of a gas that may contain dihydrogen
US9249738B2 (en) Method for automatic closed-loop control of one or more combustion temperatures in a gas turbine installation and method for determination of the water content in the working fluid of a gas turbine installation
CN110889217B (en) Correction method for outlet temperature of one-dimensional simulation model of combustion chamber of gas turbine
CN102428260B (en) Method for determining the spark advance of a heat engine
CN104344959B (en) The test method and device of single cylinder engine simulation complete machine
CN110568018B (en) Method for calculating raw coal moisture in medium-speed coal mill on line
CN104903563A (en) Method and arrangement for controlling fuel supply for gas turbine
US10415484B2 (en) Engine control system
CN113255248B (en) Method for calculating fuel machine power of single-shaft combined cycle unit
CN109655489B (en) Method and device for metering standard coal consumption of coal-fired unit
JP6907207B2 (en) Gas turbine combustion temperature control using air injection system
CN113361040B (en) Combustion chamber outlet temperature evaluation method under complete engine condition
CN114061691A (en) Method and system for measuring gas consumption of aerospace gas supply system
KR102234570B1 (en) A method for calibrating the calorific value of fuel for improving the accuracy of estimation of the efficiency of the thermal power plant
CN113032914A (en) Gas turbine performance calculation and simulation method
CN101183068B (en) Method for correcting basic state temperature absorption coefficient of light
CN108051474B (en) Method for prejudging service life of impregnating varnish storage process for air-cooled mould pressing integral-immersion steam turbine generator
CN115034083B (en) Induced draft fan power consumption correction calculation method in station service power consumption calculation of thermal power generating unit
CN204101295U (en) The test unit of single cylinder engine simulation complete machine
CN110619929B (en) Soft measurement method for flue gas flow of pulverized coal boiler of power station
US8965537B2 (en) Method for ascertaining process values for a process control

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