CZ:Cold Storage/GF Engineering Material

Move into cold storage due to editors decision. Chris Day 08:04, 2 February 2009 (UTC)

This article addresses analysis of energy conversion as it pertains to power cycles (Carnot, Brayton, Otto and Diesel), power cycle components/processes (compression, combustion and expansion) and compressible flow (nozzle, diffuser and thrust) ideal operation, T - s diagrams and performance trends.

Power cycle analysis

In this section dealing power cycles, ideal and isentropic operation of Carnot, Brayton, Otto and Diesel cycles and their preformance trends are presented.

Carnot cycle

In the presented Carnot Cycle analysis, only air is considered as the working fluid behaving as a perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a Carnot Cycle schematic layout. Air enters a compressor at point 1, and exits the compressor at point 2. Isentropic compression is considered with no entropy change. Air enters a heat exchanger -- heat addition -- at point 2, and exits the heat exchanger at point 3. At a constant temperature, heat addition takes place. Air enters a turbine at point 3, and exits the turbine at point 4. Isentropic expansion is considered with no entropy change. Air enters a heat exchanger -- heat rejection -- at point 4, and exits the heat exchanger at point 1. At a constant temperature, heat rejection takes place. It should be mentioned that air at point 1 enters the compressor and the cycle is repeated.

Figure 1 - Carnot Cycle Schematic Layout

Figure 2 presents a Carnot Cycle entropy vs temperature diagram.

Figure 1 - Carnot Cycle T - s Diagram

Figure 3 presents the Carnot Cycle efficiency as a function of the heat addition temperature. It should be noted that the compressor inlet tempearture is at standard ambient conditions of 298 [K] and 1 [atm] of absolute pressure.

Figure 1 - Carnot Cycle Efficiency

Figure 4 presents the Carnot Cycle efficiency as a function of the heat rejection temperature. It should be noted that the turbine inlet temperature is at 800 [K].

Figure 1 - Carnot Cycle Efficiency

The Carnot Cycle efficiency increases with an increase in the heat addition temperature when the heat rejection temperature does not change at all. Furthermore, the Carnot Cycle efficiency decreases with an increase in the heat rejection temperature when the heat addition temperature does not change at all.

The Carnot Cycle efficiency is not dependant on the working fluid properties.

The governing equations are as follows:

T₂/T₁ = (p₂/p₁)^(k-1)/k

T_3/T₄ = (p₃/p₄)^k-1)/k

ϰ = c_p/c_v

pv = RT

η = 1 - T₁/T₂

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

Brayton cycle

In the presented Brayton cycle analysis, only air is considered as the working fluid behaving as a perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

A gas turbine is a heat engine that uses a high temperature, high pressure gas as the working fluid. Combustion of a fuel in air is usually used to produce the needed temperatures and pressures in the gas turbine, which is why gas turbines are often referred to as combustion turbines. Expansion of the high temperature, high pressure working fluid takes place in the gas turbine. The gas turbine shaft rotation drives an electric generator and a compressor for the working fluid, air, used in the gas turbine combustor. Many gas turbines also use a heat exchanger called a recouperator to impart turbine exhaust heat into the combustor's air/fuel mixture. Gas turbines produce high quality heat that can be used to generate steam for combined heat and power and combined-cycle applications, significantly enhancing efficiency.

The initial state of the working fluid, air, is assumed to be at point 1 at the atmospheric conditions, 298 [K] temperature and 1 [atm] pressure. The working fluid is compressed, isentropically, along line 1-2 by a compressor and enters the combustor. At a constant pressure, combustion takes place, fuel is added to the combustor and the working fluid temperature raises. High temperature working fluid exits the combustor at point 3. Then the working fluid enters a gas turbine where an isentropic expansion occurs, producing power. The working fluid exits the gas turbine at point 4. It should be mentioned that the working fluid at point 1 enters the compressor and the cycle is repeated.

A schematic layout of the ideal, open Brayton cycle is shown in Figure 1.

Figure 1 - Brayton Cycle Schematic Layout

The changes in the state of the working fluid, air, are presented in a T - s diagram in Figure 2.

Figure 2 - Brayton Cycle Temperature vs Entropy Diagram

It should be pointed out that this material deals with the open Brayton cycle. The way how the T - s diagram is presented, it describes a closed Brayton cycle -- this would require a heat exchanger after point 4 where the working fluid would be cooled down to point 1 and the cycle repeats. Therefore, the T - s diagram is presented as a closed Brayton cycle to allow easier understanding and derivation of the Brayton cycle thermal efficiency -- heat addition and heat rejection.

The gas turbine and compressor are connected by shaft so the considerable amount of work done on the gas turbine is used to power the compressor.

It can be noticed from the T - s diagram that the work done on the gas turbine is greater than the work necessary to power the compressor -- constant pressure lines in the T - s diagram diverge by going to the right side (entropy wise).

The thermal cycle efficiency can be given as a function of specific external work (specific net power output) and added heat to the working fluid as follows:

η = w/q_h = (w_t - w_c)/q_h = (q_h - q_l)/q_h

or

η = 1 - q_l/q_h = 1 - (c_p(T₄ - T₁))/(c_p(T₃ - T₂)) = 1 - (T₁(T₄/T₁ - 1))/(T₂(T₃/T₂ - 1))

where

η - thermal efficiency [/]

w - specific external work (specific net power output) [kJ/kg]

w_t - expansion specific power output [kJ/kg]

w_c - compression specific power input [kJ/kg]

W_t - expansion power output [kW]

W_c - compression power input [kW]

q_h - added heat to the working fluid [kJ/kg]

q_l - heat rejected from the working fluid [kJ/kg]

c_p - specific heat at constant pressure [kJ/kg*K]

c_v - specific heat at constant volume [kJ/kg*K]

m - working fluid mass flow rate [kg/s]

For isentropic turbine and compressor

p₃/p₄ = (T₃/T₄)^ϰ/(ϰ-1)

p₂/p₁ = (T₂/T₁)^ϰ/(ϰ-1)

Knowing that

p₃/p₄ = p₂/p₁

where

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

p₁, p₂, p₃, p₄ - pressure values at points 1, 2, 3 and 4 [atm]

T₁, T₂, T₃, T₄ - temperature values at points 1, 2, 3 and 4 [K]

It follows that

T₃/T₄ = T₂/T₁

or

T₃/T₂ = T₄/T₁

Therefore, after some mathematical operations the thermal efficiency is:

η = 1 - T₁/T₂ = 1 - T₄/T₃

If the temperature ratio is substituted in terms of the pressure ratio:

η = 1 - 1/r_p^(ϰ-1)/ϰ

where

r_p - pressure ratio p₂/p₁ [/]

The derived Brayton cycle thermal efficiency depends on the pressure ratio and on the working fluid. In Figure 3, it is shown the thermal efficiency increases with the pressure ratio. In this material, air is the only working fluid analyzed. However, it can be concluded that the thermal efficiency increases with the higher value for ϰ, which is a ratio of gas specific heat values (cp/cv).

Figure 3 - Brayton Cycle Efficiency

Here, two general performance trends are considered. First, impact of the gas turbine inlet temperature and compressor pressure ratio on the Brayton cycle specific power output and second, impact of the working mass flow rate on the Brayton cycle power output.

Figure 4 presents the results of the first performance trend, while Figure 5 presents the results of the second trend.

Figure 4 - Brayton Cycle Specific Power Output

Figure 5 - Brayton Cycle Power Output

One can notice that the Brayton cycle specific power output increases with an increase in the gas turbine inlet temperature. Furthermore, the increase is greater for the higher compressor pressure ratio.

One can notice that the Brayton cycle power output increases with an increase in the working fluid mass flow rate. The increase is greater for the higher compressor pressure ratio.

The general assumptions are as follows:

Compression and expansion processes are reversible and adiabatic – isentropic. The working fluid has the same composition throughout the cycle. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

p₃/p₄ = (T₃/T₄)^ϰ/(ϰ-1)

p₂/p₁ = (T₂/T₁)^ϰ/(ϰ-1)

ϰ = c_p/c_v

pv = RT

w_c = c_p(T₂ - T₁)

W_c = c_p(T₂ - T₁)m

w_t = c_p(T₃ - T₄)

W_t = c_p(T₃ - T₄)m

η = 1 - 1/r_p^(ϰ-1)/ϰ

r_p = p₂/p₁

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

T₃ = 900, 1,200 and 1,500 [K]

p₃ = 5 and 15 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

m = 50, 100 and 150 [kg/s]

The results are as follows:

Table 1 - Brayton Cycle Efficiency vs Pressure Ratio

Pressure Ratio [/]	Efficiency [%]
5	36.92
10	48.22
15	53.87
20	57.53
25	60.16

Table 2 - Specific Power Output vs Pressure Ratio for a few Gas Turbine Inlet Values

Specific Power Output [kW/kg/s]	Gas Turbine Inlet Temperature [K]
Pressure Ratio [/]	900	1,200	1,500
5	159	270	381
15	137	300	462

Table 3 - Power Output vs Pressure Ratio for a few Mass Flow Rates

Power Output [MW]	Mass Flow Rate [kg/s]
Pressure Ratio [/]	50	100	150
5	19.1	38.1	57.2
15	23.1	46.2	69.3

The Brayton cycle thermal efficiency depends on the pressure ratio and on the working fluid. The thermal efficiency increases with the pressure ratio. Also, the thermal efficiency increases with the higher value for ϰ, which is a ratio of gas specific heat values (c_p/c_v).

The Brayton cycle specific power output increases with an increase in the gas turbine inlet temperature. Furthermore, the increase is greater for the higher compressor pressure ratio.

The Brayton cycle power output increases with an increase in the working fluid mass flow rate. The increase is greater for the higher compressor pressure ratio.

Otto cycle

In the presented Otto cycle analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains an Otto cycle pressure vs volume diagram. Air enters a cylinder at point 1 when compression starts, and compression ends at point 2. Isentropic compression is considered with no entropy change. Heat addition starts at point 2, and heat addition ends at point 3. At a constant volume, air gets heated and the working fluid temperature raises. Expansion starts at point 3, and expansion ends at point 4. Isentropic expansion is considered with no entropy change. Air heat rejection starts at point 4, and heat rejection ends at point 1. At a constant volume, air gets cooled and the working fluid temperature goes down. It should be mentioned that air at point 1 enters the compression process again and the cycle is repeated.

Figure 1 - Otto Cycle Pressure vs Volume Diagram

Figure 2 presents an Otto Cycle entropy vs temperature diagram.

Figure 2 - Otto Cycle Temperature vs Entropy Diagram

The thermal cycle efficiency can be given as a function of specific external work (specific net power output) and added heat to the working fluid as follows:

η = w/q_h = (w_e - w_c)/q_h = (q_h - q_l)/q_h

or

η = 1 - q_l/q_h = 1 - (c_v(T₄ - T₁))/(c_v(T₃ - T₂)) = 1 - (T₁(T₄/T₁ - 1))/(T₂(T₃/T₂ - 1))

where

η - thermal efficiency [/]

w - specific external work (specific net power output) [kJ/kg]

w_e - expansion specific power output [kJ/kg]

w_c - compression specific power input [kJ/kg]

W_e - expansion power output [kW]

W_c - compression power input [kW]

q_h - added heat to the working fluid [kJ/kg]

q_l - heat rejected from the working fluid [kJ/kg]

c_p - specific heat at constant pressure [kJ/kg*K]

c_v - specific heat at constant volume [kJ/kg*K]

m - working fluid mass flow rate [kg/s]

For isentropic compression and expansion

T₂/T₁ = (p₂/p₁)^(ϰ-1)/ϰ = (V_1/V₂)^(ϰ-1)

T₄T₃ = (p₄/p₃)^(ϰ-1)/ϰ = (V₃/V₄)^(ϰ-1)

Knowing that

p₃/p₄ = p₂/p₁

where

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

p₁, p₂, p₃, p₄ - pressure values at points 1, 2, 3 and 4 [atm]

T₁, T₂, T₃, T₄ - temperature values at points 1, 2, 3 and 4 [K]

It follows that

T₃/T₄ = T₂/T₁ = (V₁/V₂)^(ϰ-1) = ε^(ϰ-1)

where

ε = V₁/V₂

Therefore, after some mathematical operations the thermal efficiency is:

η = 1 - T₁/T₂ = 1 - 1/(V₁/V₂)^(ϰ-1)

If the temperature ratio is substituted in terms of the volume/compression ratio:

η = 1 - 1/_ε^(ϰ-1)

Figure 3 presents the Otto cycle efficiency as a function of the compression ratio. It should be noted that the compressor inlet temperature is at standard ambient conditions of 298 [K] and 1 [atm] of absolute pressure. Final combustion temperature is 1,500 [K].

Figure 3 - Otto Cycle Efficiency vs Compression Ratio

Figure 4 presents the Otto cycle power output as a function of combustion temperature and compression ratio. It should be noted that the number of revolutions is 60 [1/s] for a given geometry of a four cylinder and four stroke Otto engine.

Figure 4 - Otto Cycle Power Output

One can notice that the Otto cycle efficiency increases with an increase in the compression ratio. One can notice that the Otto cycle power output increases with an increase in the combustion temperature and that the Otto cycle power output is greater for higher compression ratio values.

The general assumptions are as follows:

Working fluid is air. There is no friction. Compression and expansion are isentropic -- there is no entropy change. During heat addition and heat rejection, the air temperature does change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T₂/T₁ = (V₁/V₂)^(ϰ-1)

T₃/T₄ = (V₄/V₃)^(ϰ-1)

ϰ = c_p/c_v

pv = RT

w_c = c_p(T₂ - T₁)

W_c = c_p(T₂ - T₁)m

w_e = c_p(T₃ - T₄)

W_e = c_p(T₃ - T₄)m

η = 1 - 1/ε^{(ϰ -1)}

ε = V₁/V₂

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

T₃ = 900, 1,200 and 1,500 [K]

ε = 5, 10 and 15 [/]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

The results are as follows:

Table 1 - Otto Cycle Efficiency vs Compression Ratio

Compression Ratio [/]	Otto Cycle Efficiency [%]
2.5	30.7
5	47.5
7.5	55.3
10	60.2
12.5	63.6

Table 2 - Otto Cycle Power Output

Power Output [kW]	Combustion Temperature [K]
Compression Ratio [/]	900	1,200	1,500
5	123.0	233.9	344.9
10	71.0	211.7	352.3
15	10.1	164.7	319.3

The Otto cycle efficiency increases with an increase in the compression ratio. Also, the Otto cycle power output increases with an increase in the combustion temperature and the Otto cycle power output is greater for higher compression ratio values.

Diesel cycle

In the presented Diesel cycle analysis, only air is considered as the working fluid behaving as a perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a Diesel cycle pressure vs volume diagram. Air enters a cylinder at point 1 when compression starts, and compression ends at point 2. Isentropic compression is considered with no entropy change. Heat addition starts at point 2, and heat addition ends at point 3. At a constant pressure, air gets heated and the working fluid temperature raises. Expansion starts at point 3, and expansion ends at point 4. Isentropic expansion is considered with no entropy change. Air heat rejection starts at point 4, and heat rejection ends at point 1. At a constant volume, air gets cooled and the working fluid temperature goes down. It should be mentioned that air at point 1 enters the compression process again and the cycle is repeated.

Figure 1 - Diesel Cycle Pressure vs Volume Diagram

Figure 2 presents a Diesel cycle entropy vs temperature diagram.

Figure 2 - Diesel Cycle Temperature vs Entropy Diagram

The thermal cycle efficiency can be given as a function of specific external work (specific net power output) and added heat to the working fluid as follows:

The thermal cycle efficiency can be given as a function of specific external work (specific net power output) and added heat to the working fluid as follows:

η = w/q_h = (w_e - w_c)/q_h = (q_h - q_l)/q_h

or

η = 1 - q_l/q_h = 1 - (c_v(T₄ - T₁))/(c_p(T₃ - T₂)) = 1 - (T₁(T₄/T₁ - 1))/(ϰT₂(T₃/T₂ - 1))

where

η - thermal efficiency [/]

w - specific external work (specific net power output) [kJ/kg]

w_e - expansion specific power output [kJ/kg]

w_c - compression specific power input [kJ/kg]

W_e - expansion power output [kW]

W_c - compression power input [kW]

q_h - added heat to the working fluid [kJ/kg]

q_l - heat rejected from the working fluid [kJ/kg]

c_p - specific heat at constant pressure [kJ/kg*K]

c_v - specific heat at constant volume [kJ/kg*K]

m - working fluid mass flow rate [kg/s]

For isentropic turbine and compressor

T₂/T₁ = (p₂/p₁)^(ϰ-1)/ϰ = (V₁/V₂)^(ϰ-1)

T₄/T₃ = (p₄/p₃)^(ϰ-1)/ϰ = (V₃/V₄)^(ϰ-1)

Knowing that

p₃/p₄ = p₂/p₁

where

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

p₁, p₂, p₃, p₄ - pressure values at points 1, 2, 3 and 4 [atm]

T₁, T₂, T₃, T₄ - temperature values at points 1, 2, 3 and 4 [K]

Knowing that

s₃ - s₂ = s₄ - s₁

and

s₃ - s₂ = c_pln(T₃/T₂)

s₄ - s₁ = c_vln(T₄/T₁)

It follows

(T₃/T₂)^ϰ = T₄/T₁

It follows that

T₃/T₄ = T₂/T₁ = (V₁/V₂)^(ϰ-1) = ε^(ϰ-1)

where

ε = V₁/V₂

φ = V₃/V₂

Therefore, after some mathematical operations the thermal efficiency is:

η = 1 - (T₁((T₃/T₂)^ϰ - 1)))/(ϰT₂(T₃/T₂ - 1))

If the temperature ratio is substituted in terms of the volume/compression ratio:

η = 1 - (φ^ϰ - 1)/(ϰε^(ϰ-1)(φ – 1))

Figure 3 presents the Diesel cycle efficiency as a function of the compression ratio. It should be noted that the compressor inlet temperature is at standard ambient conditions of 298 [K] and 1 [atm] of absolute pressure. The cut off ratio is 3 [/].

Figure 3 - Diesel Cycle Efficiency vs Compression Ratio

Figure 4 presents the Diesel cycle power output as a function of compression and cut off ratio values. It should be noted that the number of revolutions is 60 [1/s] for a given geometry of a four cylinder and four stroke Diesel engine.

Figure 4 - Diesel Cycle Power Output

One can notice that the Diesel cycle efficiency increases with an increase in the compression ratio. One can notice that the Diesel cycle power output increases with an increase in the compression ratio and that the Diesel cycle power output is greater for higher cut off ratio values.

The general assumptions are as follows:

Working fluid is air. There is no friction. Compression and expansion are isentropic -- there is no entropy change. During heat addition and heat rejection, the air temperature does change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T₂/T₁ = (V₁/V₂)^(ϰ-1)

T₃/T₄ = (V₄/V₃)^(ϰ-1)

ϰ = c_p/c_v

pv = RT

w_c = c_p(T₂ - T₁)

W_c = c_p(T₂ - T₁)m

w_e = c_p(T₃ - T₂) + c_p - T4)

W_e = c_p(T₃ - T₂)m + c_p(T₃ - T₄)m

η = 1 - (φ^(ϰ–1))/(ϰε^(ϰ-1)(φ - 1))

ε = V₁/V₂

φ = V₃/V₂

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

ε = 5, 7.5, 10, 12.5 and 15 [/]

φ = 3 and 4 [/]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

The results are as follows:

Table 1 - Diesel Cycle Efficiency vs Compression Ratio

Compression Ratio [/]	Diesel Cycle Efficiency [%]
5	31.4
7.5	41.7
10	48.0
12.5	52.5
15	55.8

Table 2 - Diesel Cycle Power Output

Power Output [kW]	Compression Ratio [/]
Cut Off Ratio [/]	5	7.5	10	12.5	15
3	919.0	1,230.2	1,438.7	1,701.4	1,894.4
4	1,266.8	1,733.6	2,113.9	2,440.5	2,729.9

The Diesel cycle efficiency increases with an increase in the compression ratio. Also, the Diesel cycle power output increases with an increase in the compression ratio and the Diesel cycle power output is greater for higher cut off ratio values.

Power cycle components/processes analysis

In this section dealing with power cycle components/processes, ideal compression, combustion and expansion operation and their preformance trends are presented.

Compression

In the presented compression analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a compression schematic layout. Air enters a compressor at point 1, and exits the compressor at point 2. Isentropic compression is considered with no entropy change.

Figure 1 - Compression Schematic Layout

Figure 2 presents a compression entropy vs temperature diagram.

Figure 2 - Compression Temperature vs Entropy Diagram

Figure 3 presents compression specific power input requirements for a few typical pressure ratio values. It should be noted that the air enters the compressor at standard ambient conditions of 298 [K] and 1 [atm] of absolute pressure.

Figure 3 - Compression Specific Power Input

Figure 4 presents compression power input requirements for two typical pressure ratio values and a few different air mass flow rate values.

Figure 4 - Compression Power Input

One can notice that both compression specific power input and power input increase with the pressure ratio. As the air mass flow increases, the compression power input requirements increase too.

The general assumptions are as follows:

Working fluid is air. There is no friction and heat transfer. Compression is isentropic -- there is no entropy change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T₂/T₁ = (p₂/p₁)^(ϰ-1)/ϰ

ϰ = c_p/c_v

pv = RT

w = c_p(T₂ - T₁)

W = c_p(T₂ - T₁)m

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

p₂ = 5, 10 and 15 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

m = 50, 100 and 150 [kg/s]

The results are as follows:

Table 1 - Specific Power Input vs Pressure Ratio

Pressure Ratio [/]	Specific Power Input [kW/kg/s]
5	174
10	278
15	349

Table 2 - Power Input vs Pressure Ratio for a few Mass Flow Rates

Power Input [MW]	Mass Flow Rate [kg/s]
Pressure Ratio [/]	50	100	150
5	8.68	17.37	26.06
15	17.47	34.94	52.41

Both compression specific power input and power input increase with the pressure ratio. As the air mass flow increases, the compression power input requirements increase too.

Combustion

A combustion analysis for a few typical fuel cases (carbon, hydrogen, sulfur, coal, oil and gas) when the fuel reacts with air at stoichiometric conditions is presented.

Furthermore, in the presented combustion analysis, both fuel and air are at standard inlet combustion conditions of 298 [K] and 1 [atm] of absolute pressure. Furthermore, combustion is complete and with no heat loss.

During the combustion, a large amount of reactants' chemical energy gets released in the form of thermal energy.

Enthalpy of combustion (HHV) is the difference between the reactants enthalpy value minus the combustion products enthalpy value at the standard reference temperature, which is 298 [K].

When the reactants enthalpy value is equal to the combustion products enthalpy value, one can calculate the combustion products flame temperature (adiabatic temperature).

The plot in Figure 1 depicts the reactants and combustion products enthalpy value change with an increase in the temperature.

Figure 1 - Reactants and Combustion Products Enthalpy vs Temperature

Physical properties for both reactants and combustion products are very important and need to be known in order to carry out successful combustion calculations.

The plot in Figure 2 depicts how the reactants and combustion products species enthalpy values change with the temperature. The physical properties provided in this plot come from the JANAF Thermochemical Data - Tables, 1970.

Figure 2 - Reactants and Combustion Products Species Enthalpy vs Temperature

It is interesting to note that the enthalpy value for basic combustion elements such as carbon (C), hydrogen (H₂), sulfur (S), oxygen (O₂) and nitrogen (N₂) is equal to zero at the standard combustion conditions of 298 [K] and 1 [atm].

Also, it should be mentioned that for ideal gas species, the enthalpy value is only dependent on the temperature.

In addition to knowing the reactants and combustion products physical properties, for any kind of combustion analysis and calculations, it is important to know both fuel and oxidant compositions.

For solid and liquid type fuels, the fuel composition is given on the weight basis for a unit mass amount. For the gas type fuels, the fuel composition is provided on the mole/volume basis for a unit volume amount. In this analysis, methane (CH₄) is the only gas fuel considered. In order to keep the combustion analysis simple and straightforward, the CH₄ composition is provided on the weight basis. Oxidant composition is usually given on the mole/volume basis.

Table 1 provides the fuel composition.

Table 1: Fuel Composition (Proportion by Mass)

Fuel	C	H	S	N	O	H2O	CH4
Carbon	1.000	0.000	0.000	0.000	0.000	0.000	-
Hydrogen	0.000	1.000	0.000	0.000	0.000	0.000	-
Sulfur	0.000	0.000	1.000	0.000	0.000	0.000	-
Coal	0.780	0.050	0.030	0.040	0.080	0.020	-
Oil	0.860	0.140	0.000	0.000	0.000	0.000	-
Gas	-	-	-	-	-	-	1.000

Table 2 provides the standard air composition.

Table 2: Oxidant/Air Composition

Oxidant	N [kg/kg]	O [kg/kg]	N2 [kmol/kmol]	O2 [kmol/kmol]
Air	0.767	0.233	0.790	0.210

Again, in this combustion analysis, only the stoichiometric combustion is analyzed. Results of such analysis are provided, including a combustion products composition on weight and mole/volume basis, a flame temperature, a stoichiometric ratio and a higher heating value (HHV).

Table 3 provides the combustion products composition on the weight basis.

Table 3: Combustion Products on Weight Basis (Proportion by Mass)

Fuel	CO2	H2O	SO2	N2	O2
Carbon	0.295	0.000	0.000	0.705	0.000
Hydrogen	0.000	0.255	0.000	0.745	0.000
Sulfur	0.000	0.000	0.378	0.622	0.000
Coal	0.249	0.041	0.005	0.705	0.000
Oil	0.202	0.080	0.000	0.718	0.000
Gas	0.151	0.124	0.000	0.725	0.000

Table 4 provides the combustion products composition on the volume basis.

Table 4: Combustion Products on Mole Basis (Proportion by Moles [kmol])

Fuel	CO2	H2O	SO2	N2	O2
Carbon	0.210	0.000	0.000	0.790	0.000
Hydrogen	0.000	0.347	0.000	0.653	0.000
Sulfur	0.000	0.000	0.210	0.789	0.000
Coal	0.170	0.068	0.002	0.759	0.000
Oil	0.132	0.129	0.000	0.739	0.000
Gas	0.095	0.190	0.000	0.715	0.000

When considering coal, oil and gas as the fuel, coal has the largest amount of CO₂ in the combustion products on both weight and mole basis.

Table 5 provides the combustion products flame temperature, stoichiometric ratio and the fuel higher heating value.

Table 5: Other Fuel Characteristics

Fuel	Flame Temperature [K]	Stoichiometric Ratio (See Note Below!)	HHV [kJ/kg]
Carbon	2,460	11.444	32,779.8
Hydrogen	2,525	34.333	141,866.8
Sulfur	1,972	4.292	9,261.3
Coal	2,484	10.487	32,937.9
Oil	2,484	14.649	48,055.1
Gas	2,327	17.167	50,151.2
Note: Stoichiometric ratio is the mass of air required for complete combustion of a unit mass of fuel. Thus, 1 [kg] of carbon fuel requires 11.444 [kg] of air for complete, ideal combustion.

Stoichiometric ratio is the mass of air required for complete combustion of a unit mass of fuel. Thus, 1 [kg] of carbon fuel requires 11.433 [kg] of air for complete, ideal combustion.

Today, global warming is becoming more evident and it is being said that it is primarily caused by CO₂ emissions. A detailed combustion analysis, as it is provided here, can be very useful in determining different fuel and technology scenarios that would result in the reduction of current CO₂ emissions.

The general assumptions:

Both fuel and air are at standard inlet combustion conditions of 298 [K] and 1 [atm] of absolute pressure.

Furthermore, combustion is complete and with no heat loss.

The governing equations are as follows:

Enthalpy of combustion (HHV) is the difference between the reactants enthalpy value minus the combustion products enthalpy value at the standard reference temperature, which is 298 [K].

When the reactants enthalpy value is equal to the combustion products enthalpy value, one can calculate the combustion products flame temperature (adiabatic temperature).

Input data are as follows:

Table 1: Fuel Composition (Proportion by Mass)

Fuel	C	H	S	N	O	H2O	CH4
Carbon	1.000	0.000	0.000	0.000	0.000	0.000	-
Hydrogen	0.000	1.000	0.000	0.000	0.000	0.000	-
Sulfur	0.000	0.000	1.000	0.000	0.000	0.000	-
Coal	0.780	0.050	0.030	0.040	0.080	0.020	-
Oil	0.860	0.140	0.000	0.000	0.000	0.000	-
Gas	-	-	-	-	-	-	1.000

Table 2: Oxidant/Air Composition

Oxidant	N [kg/kg]	O [kg/kg]	N2 [kmol/kmol]	O2 [kmol/kmol]
Air	0.767	0.233	0.790	0.210

The results are as follows:

Table 3: Combustion Products on Weight Basis (Proportion by Mass)

Fuel	CO2	H2O	SO2	N2	O2
Carbon	0.295	0.000	0.000	0.705	0.000
Hydrogen	0.000	0.255	0.000	0.745	0.000
Sulfur	0.000	0.000	0.378	0.622	0.000
Coal	0.249	0.041	0.005	0.705	0.000
Oil	0.202	0.080	0.000	0.718	0.000
Gas	0.151	0.124	0.000	0.725	0.000

Table 4: Combustion Products on Mole Basis (Proportion by Moles [kmol])

Fuel	CO2	H2O	SO2	N2	O2
Carbon	0.210	0.000	0.000	0.790	0.000
Hydrogen	0.000	0.347	0.000	0.653	0.000
Sulfur	0.000	0.000	0.210	0.789	0.000
Coal	0.170	0.068	0.002	0.759	0.000
Oil	0.132	0.129	0.000	0.739	0.000
Gas	0.095	0.190	0.000	0.715	0.000

Table 5: Other Fuel Characteristics

Fuel	Flame Temperature [K]	Stoichiometric Ratio (See Note Below!)	HHV [kJ/kg]
Carbon	2,460	11.444	32,779.8
Hydrogen	2,525	34.333	141,866.8
Sulfur	1,972	4.292	9,261.3
Coal	2,484	10.487	32,937.9
Oil	2,484	14.649	48055.1
Gas	2,327	17.167	50,151.2
Note: Stoichiometric ratio is the mass of air required for complete combustion of a unit mass of fuel. Thus, 1 [kg] of carbon fuel requires 11.444 [kg] of air for complete, ideal combustion.

Expansion

In the presented expansion analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains an expansion schematic layout. Air enters a turbine at point 1, and exits the turbine at point 2. Isentropic expansion is considered with no entropy change.

Figure 1 - Expansion Schematic Layout

Figure 2 presents an expansion entropy vs temperature diagram.

Figure 2 - Expansion Temperature vs Entropy Diagram

Figure 3 presents expansion specific power output values for a few typical pressure ratio values. It should be noted that the air enters the turbine at the temperature of 1,500 [K] and the turbine exhaust pressure is always equal to the standard ambient pressure -- 1 [atm] of absolute pressure.

Figure 3 - Expansion Specific Power Output

Figure 4 presents expansion power output values for two typical pressure ratio values and a few different air mass flow rate values.

Figure 4 - Expansion Power Output

One can notice that both expansion specific power output and power output increase with the pressure ratio. As the air mass flow increases, the expansion power output values increase too.

The general assumptions are as follows:

Working fluid is air. There is no friction and heat transfer. Expansion is isentropic -- there is no entropy change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T₁/T₂ = (p₁/p₂)^(ϰ-1)/ϰ

ϰ = c_p/c_v

pv = RT

w = c_p(T₁ - T₂)

W = c_p(T₁ - T₂)m

Input data are as follows:

T₁ = 1,500 [K]

p₁ = 5, 10 and 15 [atm]

p₂ = 1 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

m = 50, 100 and 150 [kg/s]

The results are as follows:

Table 1 - Specific Power Output vs Pressure Ratio

Pressure Ratio [/]	Specific Power Output [kW/kg/s]
5	555
10	726
15	811

Table 2 - Power Output vs Pressure Ratio for a few Mass Flow Rates

Power Output [MW]	Mass Flow Rate [kg/s]
Pressure Ratio [/]	50	100	150
5	27.76	55.51	83.27
15	40.56	81.13	121.69

Both expansion specific power output and power output increase with the pressure ratio. As the air mass flow increases, the expansion power output values increase too.

Compressible flow analysis

In this section dealing with compressible flow, ideal and isentropic operation of nozzle, diffuser and thrust and their preformance trends are presented.

Nozzle

In the presented nozzle analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a nozzle schematic layout. Air enters a nozzle at point 1, and exits the nozzle at point 2. Isentropic expansion is considered with no entropy change.

Figure 1: Nozzle Schematic Layout

Figure 2 presents a nozzle entropy vs temperature diagram.

Figure 2 - Nozzle Temperature vs Entropy Diagram

Figure 3 presents nozzle performance -- stagnation over static temperature and pressure values -- as a function of the Mach Number. Only subsonic nozzle operation is considered. It should be noted that air enters the nozzle at the stagnation conditions of 1,500 [K] and 10 [atm] of absolute pressure.

Figure 3 - Nozzle Performance vs Mach Number

One can notice that nozzle stagnation over static temperature and pressure ratio values increase with an increase of the Mach Number.

The general assumptions are as follows:

Working fluid is air. There is no friction and heat transfer. Expansion is isentropic -- there is no entropy change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T_t /T = (1 + M²(ϰ - 1)/2)

p_t/p = (1 + M²(ϰ - 1)/2)^ϰ/(ϰ-1)

T_t/T = (p_t/p)^(ϰ-1)/ϰ

v = (2c_p(T_t - T)^)1/2

v_s = (ϰRT)^1/2

M = v/v_s

ϰ = c_p/c_v

pv = RT

Input data are as follows:

T₁ = 1,500 [K]

p₁ = 10 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

M = 0.39, 0.67 and 1 [/]

The results are as follows:

Table 1 - Nozzle Performance vs Mach Number

Mach Number [/]	Temperature Ratio [/]	Pressure Ratio [/]
0.39	1.03	1.11
0.67	1.09	1.35
1.00	1.19	1.86

Nozzle stagnation over static temperature and pressure ratio values increase with an increase of the Mach Number.

Diffuser

In the presented diffuser analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a diffuser schematic layout. Air enters a diffuser at point 1, and exits the diffuser at point 2. Working fluid inlet velocity gets reduced to zero resulting in the stagnation temperature and pressure increase. Isentropic process is considered with no entropy change.

Figure 1 - Diffuser Schematic Layout

Figure 2 presents a diffuser entropy vs temperature diagram.

Figure 2 - Diffuser Temperature vs Entropy Diagram

Figure 3 presents diffuser performance -- stagnation over static temperature and pressure values -- as a function of the Mach Number. Only subsonic diffuser operation is considered. It should be noted that the air enters the diffuser at the static conditions of 298 [K] and 1 [atm] of absolute pressure.

Figure 3 - Diffuser Performance vs Mach Number

One can notice that diffuser stagnation over static temperature and pressure ratio values increase with an increase of the Mach Number.

The genaral assumptions are as follows:

Working fluid is air. There is no friction and heat transfer. Isentropic process -- there is no entropy change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T_t/T = (1 + M²(ϰ - 1)/2)

p_t/p = (1 + M²(ϰ - 1)/2)^ϰ/(ϰ-1)

T_t/T = (p_t/p)^(ϰ-1)/ϰ

T_t = T + v²/(2c_p)

v_s = (ϰRT)^1/2

M = v/v_s

ϰ = c_p/c_v

pv = RT

Input data are as follows:

T₁ = 298 [K]

p₁ = 1 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

M = 0.29, 0.58 and 0.95 [/]

The results are as follows:

Table 1 - Diffuser Performance vs Mach Number

Mach Number [/]	Temperature Ratio [/]	Pressure Ratio [/]
0.29	1.017	1.06
0.58	1.067	1.25
0.95	1.182	1.80

Diffuser stagnation over static temperature and pressure ratio values increase with an increase of the Mach Number.

Thrust

In the presented thrust analysis, only air is considered as the working fluid behaving as the perfect gas -- specific heat has a constant value. Ideal gas state equation is valid -- pv = RT.

Figure 1 contains a thrust schematic layout. Air enters a nozzle at point 1, and exits the nozzle at point 2. Isentropic expansion is considered with no entropy change.

Figure 1 - Thrust Schematic Layout

Figure 2 presents a thrust entropy vs temperature diagram.

Figure 2 - Thrust Temperature vs Entropy Diagram

Figure 3 presents thrust performance as a function of the nozzle inlet stagnation temperature and pressure for a few fixed values such as: working fluid mass flow rate, nozzle outlet Mach Number and ambient pressure. Only subsonic nozzle operation is considered.

Figure 3 - Thrust Performance vs Inlet Stagnation Temperature and Pressure

One can notice that thrust values increase with an increase of the inlet stagnation temperature and pressure.

The general assumptions are as follows:

Working fluid is air. There is no friction and heat transfer. Expansion is isentropic -- there is no entropy change. Ideal gas state equation is valid -- pv = RT. Air behaves as a perfect gas -- specific heat has a constant value.

The governing equations are as follows:

T_t/T = (1 + M²(ϰ - 1)/2)

p_t/p = (1 + M²(ϰ - 1)/2)^ϰ/(ϰ-1)

T_t/T = (p_t/p)^{(ϰ -1)/k}

v = (2c_p(T_t – T))^1/2

v_s = (ϰRT)^1/2

M = v/v_s

ϰ = c_p/c_v

pv = RT

Thrust = vm + (p - p_a)A

Input data are as follows:

T₁ = 900, 1,200 and 1,500 [K]

p₁ = 5, 10 and 15 [atm]

c_p = 1.004 [kJ/kg*K]

ϰ = c_p/c_v - for air ϰ = 1.4 [/]

m = 1 [kg/s]

M = 0.85 [/]

p_a = 1 [atm]

The results are as follows:

Table 1 - Thrust Performance vs Inlet Stagnation Temperature and Pressure

Thrust [N]	Inlet Stagnation Temperature [K]
Inlet Stagnation Pressure [atm]	900	1,200	1,500
5	797.7	922.3	1,031.1
10	873.5	1,009.5	1,128.6
15	898.7	1,038.7	1,161.3

Thrust values increase with an increase of the inlet stagnation temperature and pressure.

References

JANAF Thermochemical Data - Tables, 1970