Online course and simulator for engineering thermodynamics

Generating the equations of a Thermoptim gas turbine cycle model

Modèle Thermoptim

It is a very simple cycle, the synoptic view of which is given below:

Synoptic view of a simple gas turbine

The diagram and project files are given below. Please note that they require the use of Thermoptim in English, i.e. with the inth2.zip file of this language.

GT_comb.dia

GT_comb.prj

Raw Generated Equations

There are 42 of them. They are given in this file.

GT_comb.eq

Redundancies

The redundancies identified are as follows:

Redundancy for m_dot_gasoutlet:

# 1 m_dot_gasoutlet = m_dot_turbine, equation: 1

# 2 m_dot_gasoutlet = 1.01827, equation: 2

Redundancy for m_dot_fuel:

# 1 m_dot_fuel = 0.015643905, equation: 3

# 2 m_dot_fuel = abs(Q_dot_combustionchamber/DeltaHr_combustionchamber), equation: 34

Redundancy for m_dot_compressor:

# 1 m_dot_compressor = m_dot_airinlet, equation: 8

# 2 m_dot_compressor = m_dot_airinlet, equation: 9

These redundancies relate to flow rates and are easy to resolve. The equations to keep are 1, 34 and 8, with 9 being a duplicate.

Uninitialized Variables

Il n'y a pas de variable non initialisée

List of Equation Groups

Group 1

Group 1 (13 Variables) : [T_fuel, p_2, p_4, T_3, p_fuel, m_dot_fuel, etaT_turbine, a_combustionchamber, p_airinlet, T_airinlet, m_dot_gasoutlet, m_dot_airinlet, etaT_compressor]

Group 1 (Equations):

m_dot_gasoutlet = 1.01827

m_dot_fuel = 0.015643905

m_dot_airinlet = 1.0

T_airinlet = 25.0

p_airinlet = 1.0

etaT_compressor = 0.85

p_2 = 16.0

etaT_turbine = 0.85

p_4 = 1.0

T_3 = 1065.0

a_combustionchamber = 4

T_fuel = 15.0

p_fuel = 20.0

We find here all the equations providing the data for the problem, except T_3 identified as missing.

Group 2

Group 2 (5 Variables) : [DeltaHr_combustionchamber, p_3, h_airinlet, m_dot_compressor, h_fuel]

Group 2 (Equations):

h_airinlet = calcH_TP("air";T = T_airinlet ;P = p_airinlet)

m_dot_compressor = m_dot_airinlet

m_dot_compressor = m_dot_airinlet

DeltaHr_combustionchamber = (-(-74850) +(-393520)+a_combustionchamber/2*(-242000))/16

p_3 = p_2

h_fuel = calcH_TP("CH4 ` methane";T = T_fuel ;P = p_fuel)

This second group corresponds to equations that allow the calculation of new variables through simple substitution of those from the first group.

Group 3

Group 3 (2 Variables) : [s_airinlet, m_dot_combustionchamber]

Group 3 (Equations):

s_airinlet = calcS_PH("air";P = p_airinlet;H = h_airinlet)

m_dot_combustionchamber = m_dot_compressor + m_dot_fuel

This third group corresponds to equations that allow the calculation of new variables through simple substitution of those from the first and second groups.

The process is repeated in the subsequent groups.

Unresolved Equations

There are no unresolved equations.

These are the equations that either depend on the properties of the fluid or cannot be directly solved.

Conversion to EES format

EES is a solver developed by f-Chart, which requires a license. The conversion results in a file that can be processed by the solver. The equations for calculating the fluid properties converted to this format are given below, the others remaining unchanged:

//Equation: 7

h_airinlet = enthalpy(air;T = T_airinlet) // Downstream point - air inlet

s_airinlet = entropy(air;P = p_airinlet;H = h_airinlet) // Upstream point - air inlet - Downstream point - 2

//Equation: 11

hs_2 = enthalpy(air;P = p_2;S = s_airinlet) // Downstream point - 2

//Equation: 14

T_2 = temperature(air;H = h_2) // Downstream point - 2

//Equation: 18

s_3 = entropy(burnt gases;P = p_3;H = h_3) // Upstream point - 3 - Downstream point - 4

//Equation: 19

hs_4 = enthalpy(burnt gases;P = p_4;S = s_3) // Downstream point - 4

//Equation: 21

h_4 = h_3 - etaT_turbine*(h_3 - hs_4) // Upstream point - 3 - Downstream point - 4

//Equation: 22

T_4 = temperature(burnt gases;H = h_4) // Downstream point - 4

//Equation: 23

s_4 = entropy(burnt gases;P = p_4;H = h_4) // Entropy

//Equation: 28

lambda_combustionchamber = LAMBD(T_2;T_3;a_combustionchamber)// air factor lambda

//Equation: 29

h_3 = h_products(T_3;a_combustionchamber;lambda_combustionchamber)// enthalpy of the reactants

//Equation: 30

hfict_2 = h_products(T_2;a_combustionchamber;lambda_combustionchamber)// enthalpy of a fictitious inlet point for calculating the heat released

//Equation: 38

h_fuel = enthalpy(CH4 ` methane;T = T_fuel) // Fuel point – fuel

As explained in the section "Calculations Downstream of a Combustion Chamber," all equations concerning the properties of the 'burnt gases' substance must be adapted to use the two calculation functions defined at the beginning of the generated file. In equation 38, the substance must be renamed simply to CH4.

It should be noted that EES uses the decimal separator defined in the Windows regional settings, which is '.' in the United States and ',' in Europe. Similarly, in function arguments, ',' is used in the United States and ';' in Europe. It may therefore be necessary to replace these separators according to the Windows settings.

The file that can be resolved in EES is provided below.

GT_combEES.eq

Other examples not commented

Regenerative gas turbine

Synoptic view of a regenerative gas turbine

regen_GT_CH4.dia

regen_GT_CH4.prj

Raw Equations File

regen_GT_CH4.eq

Executable file in EES

regen_GT_CH4_EES.eq

High-temperature nuclear power plant with intermediate heat exchanger

This is a helium Brayton cycle envisaged for high-temperature nuclear reactors (900 °C) of about a hundred MWe, using low-enriched uranium pebbles coated in carbon (which acts as a moderator) and using helium as a heat transfer fluid. It is a variant of the cycle presentedthe in guided exploration CTRN-5.

Its synoptic view is given below:

Synoptic view of a high-temperature nuclear power plant cycle with an intermediate exchanger

PBMR_IHX_N2_En.dia

PBMR_IHX_N2_En.prj

File of the104 raw equations

PBMR_IHX_Polytr.eq

Executable file in EES (117 equations)

PBMR_IHX_Polytr_EES.eq

copyright R. Gicquel v2024.3

Created with Scenari (new window)