1. Introduction
Energy is a crucial item in our daily life for almost everything, and its production, conversion, and uses closely affect the environment and sustainable development. Humanity has long depended on energy in the form of power and heat, among others, and has devised means of supplying these using power generating systems and heat producing systems which convert a primary energy source into these outputs. Cogeneration or combined heat and power (CHP) plants, which produce more than the output from a single energy source, provide an attractive approach due to their higher thermal efficiency than individual systems and lower environmental impact. This has enabled the widespread use of these systems all around the world.
Energy analysis, which is based on the first law of thermodynamics, provides a means to analyze the performance of thermal systems according to certain design parameters and operating conditions. However, it does not provide a clear picture of thermodynamic efficiency and losses. Exergy analysis overcomes these deficiencies and can help identify pathways to sustainable development [
1]. Exergy is a useful tool for determining the location, type, and true magnitude of exergy losses, which appear in the form of either exergy destruction or waste exergy emission. Therefore, exergy can assist in developing strategies and guidelines for more effective use of energy resources and technologies [
2]. Exergoenvironmental analysis likewise is an analysis method coupling the principles of thermodynamics and environics to assess the environmental impacts of the system in the view of thermodynamics [
3].
Numerous studies have utilized exergy and exergoenvironmental analysis to examine the inefficiencies and the impact of energy systems on the environment. Karklina et al. [
4] performed an exergy analysis in a biomass combined heat and powerplant in Jelgava, finding about 64% of exergy destruction just in the steam generator, followed by the heat exchangers of the district heating system. Ahmadi et al. [
5] performed an energy and exergoenvironmental analysis on a system designed to produce multiple outputs. The results showed that the system, compared to individual systems, exhibited higher exergy efficiencies and lower environmental impacts, suggesting that multigeneration can help mitigate greenhouse gas emissions. Eduardo et al. [
6] also performed an exergy and exergoenvironmental analysis of a solar hybrid trigeneration system. Their results led to the conclusion that the solution pump and high-pressure generator are the highest sources of inefficiency in the system, and the exergoenvironmental performance of the system can be improved by enhancing the exergy efficiency of these components. Ratlamwala et al. [
7] performed a comparative environmental compact assessment of two separate hydrogen and cooling cogeneration systems. Indices such as the environmental impact factor, environmental impact coefficient, environmental impact index, environmental impact improvement, exergetic stability factor, and exergetic sustainability factor were used in a bid to determine which system has a more favourable impact on the environment. The electrolyzer was found to score better than the steam methane reforming hydrogen production method.
Several models of powerplants with different energy sources, both conventional and renewable sources, have been designed and developed using software for powerplant simulation that can simulate the systems in both steady and dynamic states. Many researchers make use of commercial software such as Aspen Plus Dynamics, APROS (Advanced Process Simulation Software), ThermoSysPro (Modelica-based) library developed by EDF, SimECS developed by Delft University of Technology and used for dynamic simulation of energy conversion system [
8], Dymola, and VMGSim.
Different models of powerplants and their components have been designed and simulated by researchers with their steady and transient performance observed and analyzed. Alobaid et al. [
9] made use of APROS software in designing a simulation model of static and dynamic characteristics. The designed model is applied to the study of the start-up process in power plants. Starkloff et al. [
10] developed a full-scale dynamic model of a large-scale coal-fired powerplant. It investigated operation flexibility by making use of APROS. A simulation model of the dynamic operation of a regenerative gas turbine was developed by [
11], which analyzed the dynamic characteristics of a small engine with a static recuperator. Other researchers performed the numerical analysis of combined cycle powerplants and validated the models against measurements at different steady states [
12,
13].
Several other studies have been conducted on the dynamic modelling of powerplants. Maffezzoni [
14] performed an analysis which identified those basic concepts that are usually accepted to form the primary state of the knowledge in powerplant modelling and simulation. Lu [
15] explained the modelling techniques and carried out an analytical study on the fundamental steady state and dynamic models for powerplant components. In his study, a 677 MW coal and gas powerplant was modelled using dynamic characteristics with MATLAB and SIMULINK. Colonna and Putten [
16] analyzed and validated SimECS, which is dynamic modelling software for energy systems. SimECS then underwent development at the Delft University of Technology and was continually utilized in the modelling of small steam cycle powerplants [
17]. Due to the transient operation of dynamic models, it became important to utilize dynamic simulation tools for modelling such systems (powerplants). Celis et al. [
18] carried out a dynamic model simulation study on a steam turbine for a full-scope powerplant simulator. Oko and Wang [
19] also made use of Gproms in designing a dynamic model of a 500 MW coal-fired subcritical powerplant. A 5% correction error was recorded in the model validation, in comparison with the actual plant measurement. Wang et al. [
20] developed dynamic models of a 660 MW supercritical unit with the GSE software to obtain the energy consumption characteristics of power units during transient cycling processes. The study considered the heat storage, control systems unit, and energy utilization of the powerplant. The different outputs of the steady-state condition and at different cycling rates were measured and compared. The system showed an increment in the standard coal consumption rate (SCCR) by 3.57 g/kWh in the cycling rate as compared to the steady-state model. In their research, Powell and Edgar [
21] developed a dynamic simulation of a thermal energy storage system incorporated into a parabolic trough concentrated solar power (CSP) system. The thermal storage system (TSS) utilized the two-tank direct method. Xu et al. [
22] considered all the working parameters of TSS of a 1 MW solar tower powerplant by using modular modelling methods. Chen et al. [
23] developed an integrated powerplant model. Steady-state data of a subcritical powerplant with reheat and regenerative cycle was used in validating the model designed. The nominal power and efficiency of the system were 605 MW and 38.3%, respectively. Liu et al. [
24] worked on concentrated solar power plants by providing a comprehensive summary of both in operation and under construction.
There have been numerous studies in the literature on combined heat and power systems; however, this paper contributes to the literature by expanding on the effects of different variables such as the pressure ratio for improving the performance of these systems. The significance includes a better understanding of the practical implementation of multigeneration systems. Furthermore, trigeneration can be understood through exergoenvironmental analysis, although these analyses have not yet been published, particularly for systems based on organic Rankine cycles. The novelty of this study also stems from the utilization of the Thermolib software for the analysis, as few studies have utilized its robust computing capacity for such analyses. The present paper analyzes a combined heat and powerplant using the Thermolib extension library in MATLAB/Simulink. The behavior of the system analyzed consists of a Brayton cycle using methane gas as fuel, Rankine cycle, and ammonia heat pump, with the system performance and efficiency observed for varied component parameters. The effect of varying the pressure ratio of the Brayton cycle on the performance of the Brayton cycle, the combined gas and steam cycles, and the overall system is investigated. Similarly, the effect of varying the Rankine cycle turbine inlet pressure on the performance of the Rankine cycle and the overall system is investigated. The heat pump compressor pressure ratio is also varied to analyze its impact on the performance of the heat pump and the total system. The effect of varying the inlet energy from the combustion chamber of the Brayton cycle on the performance of each of the components is studied. Furthermore, its effect on the total system is analyzed. Lastly, an exergoenvironmental analysis is performed using the engineering equation solver (EES) to examine its environmental usefulness, appropriateness, and stability.
2. Materials and Methods
The system consists of a gas turbine (Brayton) cycle as the topping cycle, a steam turbine (Rankine) cycle next below it, and then an ammonia heat pump that can be used for space heating, as displayed in
Figure 1. Air at ambient temperature and pressure (state 1) is fed into the compressor of the gas turbine cycle. The air at high pressure (state 2) is then transferred into the mixing chamber where it is mixed with the fuel (state 3) needed for the gas power plant. Methane gas (CH
4) is used from the Thermolib setup block, and its data are loaded into the MATLAB workspace. In the presence of air and at constant pressure, the fuel is burnt in the combustion chamber. This gives rise to a high-temperature gas (state 4) which is transferred into the turbine. Inside the turbine, the high-temperature gas expands and produces power or electricity. The exhaust gas (state 5) which leaves the turbine possesses considerably high thermal energy, i.e., it leaves at a very high temperature (about 650 °C); therefore, it is used as the heat source for the steam cycle.
The exhaust gas from the gas cycle is transferred into the steam cycle by entering the heat recovery steam generator (HRSG). The water from the steam cycle (state 8) which enters the HRSG comes out as steam with high pressure and temperature (state 9). This steam enters the steam turbine and expands, thereby generating additional electricity. The steam leaving the turbine after expansion (state 10) is then passed through a condenser where it is cooled down into a saturated liquid, using a stream of water/air (states 16 and 17). The saturated liquid (state 11) leaving the condenser is then pumped into the HRSG to recycle the process. The exhaust gas after leaving the heat recovery steam generator (state 6) still possesses enough energy which is fed into the ammonia heat pump. Here, the ammonia flowing into the evaporator (state 8) receives heat from the exhaust gas and leaves at a higher temperature (state 13). The ammonia is then compressed into the high-pressure level of the heat pump, thereby causing its temperature to be further increased. The hot fluid is then passed through the condenser where it exchanges heat with the process air at ambient temperature (state 18), which leaves it at a much higher temperature (state 19). The ammonia leaving the condenser (state 14) is then passed through an expansion valve to reduce its pressure (state 12) to that of the HRSG.
2.1. The Model and Components
The information regarding the components of the proposed system is presented in this section.
2.1.1. Methane Gas
The fossil fuel is combusted in the combustion chamber part of the gas cycle section of the system. This, therefore, entails the addition of heat through a chemical reaction. The combustion process of natural gas is less complex compared to that of coal and biomass. The model of the combustion chamber has the oxidizer inlet where the compressed wet air enters the chamber. It also has a fuel inlet where the natural gas used (methane) enters and the flue gas outlet.
The combustion chemical reaction in the combustion chamber assumes the following form:
In a real practical reactor or combustion chamber, other gaseous byproducts are given out such as nitrogen, as air contains other elements other than oxygen.
Equation (2) is an expression of the mass balance.
Equation (3) is an expression of the energy balance.
where the subscript
refers to flue gas, and LHV represents the lower heating value of the fuel.
Upon ascertaining the fuel and mass flow rates entering the combustion chamber, the air-to-fuel ratio can be computed. According to [
25], the reaction assumes the following form:
where the
term ensures the availability of excess oxygen in the exhaust gases released at the outlet of the combustion chamber.
2.1.2. Brayton Cycle
A compressor is a turbo machine. The compressor increases the pressure of the compressible fluids. The ratio of the outlet pressure to the inlet pressure is known as the compressor pressure ratio, which shows the increase in pressure in the system. The mass balance of the compressor is defined as
The energy balance is given as
In order to calculate the required performance,
where
stands for the specific heat capacity of air and at the average temperature, and
For isentropic efficiency of the compressor,
A combined cycle powerplant comprises the gas turbine as a section of the system. The enthalpy of the flue gas is transformed into mechanical energy by the gas turbine, which in turn drives the generator that produces the electrical energy. The other output of the gas turbine allows the waste heat from the turbine to pass out. This waste heat can be passed out to a heat recovery steam generator (HRSG) which is used to produce steam, as well as drive a steam turbine. Various fuels can be used to operate a gas turbine. These fuels include natural gas (which is used for the system in this paper), crude oil, and biogas. The fuel range of a gas turbine can be extended to include coal and also biomass by applying the integrated gasification combined cycle (IGCC) [
26].
Considering that heavy-duty gas turbines have a much higher inlet temperature, in comparison with the steam turbine, they need cooling of their first row of blades. Due to this cooling, the mathematical model of the gas turbine is different from that of the steam turbine. Hence, the whole gas turbine has to be separated into different sections due to the presence of air cooling. The combustion gas mixes with the air cooling from a section, and this is computed in the mass balance equation of the turbine section. Founded on actual data, Wang and Leithner [
27] proposed a formula for the calculation of the mass flow of the cooling air needed at nominal load:
where
is the total air mass flow rate, while
is the inlet temperature in °C.
At partial loads, Palmer and Erbes [
28] developed the following equation:
where the subscript
nom denotes the nominal state.
2.1.3. Rankine Cycle
The transformation of the steam’s enthalpy into mechanical energy is performed by the steam turbine. Hence, the enthalpy before and after the turbine is used in computing the efficiency and the nominal values. The steam turbine converts the steam’s enthalpy into mechanical energy. Therefore, the enthalpy is determined by calculating the enthalpy value before the turbine unit, as well as the output and the marginal values after the turbine unit. The mechanical power generated at one end of the turbine transfers to the electrical generator. In the momentum and energy computations, the pressure and enthalpy drops are added as source terms. The pressure drop is a function of the nominal mass flow rate
at the turbine inlet nominal pressure
and the turbine outlet’s nominal pressure
.
The specific enthalpy drop in a turbine is computed using Equation (14) in the situation when the fluid is steam.
The specific enthalpy drop is computed using Equation (15) in the situation when there are water droplets in the fluid.
2.1.4. Ammonia Heat Pump
The use of naturally occurring and environmentally friendly substances as working fluids in heat pump systems helps keep the environment safe. Ammonia and hydrocarbons remain one of the most important natural working fluids. Ammonia heat pumps, therefore, provide an energy-efficient alternative for both heating and cooling for residential applications such as space heating, as well as heating nonresidential buildings. As shown in the system model of
Figure 1, the ammonia heat pump consists of four main elements:
Evaporator: Heat is received from the heat source and, through the regulation of the evaporator pressure with that of the expansion valve, ensures that, at the desired temperature, the cooling agent begins to boil.
Compressor: This increases the pressure in the heat pump system, thereby raising the temperature of the liquid vapor.
Condenser: In this element of the heat pump system, the cooling agent provides heat to the internal heating system, which is generally a water circulation system of heat supply to the radiators, in addition to the system supplying hot water.
Expansion valve or throttle valve or pressure regulator: This lowers the pressure and removes the cooling agent back to the evaporator, totally changing it into a liquid state so that it can pick up the heat once more [
29].
These are the four main parts of the ammonia heat pump, which are designed to be connected through a closed piping system. Ammonia, which is the cooling agent, does not harm the environment.
The COP for the heat pump system is expressed as
where
is the amount of heat is transferred to or from a cold reservoir, and
is the work supplied in moving heat from the cold reservoir to the hot reservoir.
Therefore, for the proposed system, the COP for the heat pump cycle is calculated as
where
represents the energy output of the condenser of the heat pump, and
represents the work performed by the compressor of the heat pump.
2.2. Exergy and Exergoenvironmental Analysis
The condition when the maximum useful work of a system is obtained when a system reversibly reaches equilibrium with its environment is termed exergy. Exergy is widely preferred in the analysis of the performance of thermodynamic systems because it deals with the quality rather than the quantity of energy, which enables losses and inefficiencies to be evaluated realistically [
30].
The chemical exergy of the methane fuel can be determined using the exergy ratio relation,
where
= 1.06 [
31], and
= 50,023 kJ/kg.
The exergy balance of the combustion chamber is given as
where
is the heat output of the combustion chamber,
is the mass flow rate of the fuel, and
is the exergy destruction in the combustion chamber. The exergy efficiency of the combustion chamber is calculated using
The exergy destruction rate of the other components in the multigeneration system can be obtained using the energy balance equations presented in
Table 1.
An exergoenvironmental analysis of a system is used to measure the impact of the system on the environment. The six parameters used are described below.
Exergoenvironmental impact factor (
): This shows how a system’s environmental impact can be reduced by decreasing its irreversibility. It is determined using
where
and
represent total exegetic destruction in the system and the overall exegetic destruction supplied to the system.
Exergoenvironmental impact coefficient (
): The value of
for an ideal case of a system is one to show that there is no exergy destruction. This is defined as
where
is the exergetic efficiency of the system.
Exergoenvironmental impact index (
): This is used to indicate how much damage a system exerts on the environment as a result of its exergy losses and exergy destruction. It is calculated as
Exergoenvironmental impact improvement (
): Reducing the exergoenvironmental impact index of a system will increase how appropriate the system is to the environment. Therefore, the exergoenvironmental improvement of a system is derived from its exergoenvironment impact index using
Exergy stability factor (
): A system is considered a stable one if its exergy stability factor is close to one. It is defined as
where
represents the desired output exergy.
Exergetic sustainability index (
): A system’s exergetic sustainability index is a product of its exergy stability factor and its exergoenvironmental impact improvement. It is desired to be as high as possible to illustrate its usefulness to the environment. Its value is derived using
3. Simulation Scenarios and Results
The parameters for the model inputs used in the simulation are given in
Table 2 and the results obtained from the simulation using the input parameters given in
Table 2 are presented in
Table 3. The total net output of the system is 24,748 kW and the heat pump supplies 1385 kW of heat. The efficiency of the Brayton, Rankine, and overall cycles is 33%, 36%, and 62%, respectively. Furthermore, parametric studies are carried out to demonstrate the performance of the system through the variation of some design parameters for better performance analysis and the results are presented in this section.
The impact of using different pressure ratios in the Brayton cycle on the compressor work, output power of the Brayton cycle, the combined power output of the Brayton and Rankine cycles, and overall power output is observed in
Figure 2. An increase in the pressure ratio from 3 to 9 causes the compressor work to increase from 4.7 MW to 11.1 MW; this is because the compressor will have to perform more work to increase the exit pressure of the air after passing through it. The net power output from the Brayton cycle increases from 7.2 MW to 12.3 MW, while the power output from the combined Brayton and Rankine cycles increases from 14.6 MW to 19.6 MW. The net power output of the overall system is derived by subtracting the amount of work needed to drive the compressor of the heat pump compressor from the net power output of the combined power cycles.
Figure 2 shows that increasing the Brayton cycle pressure ratio increases the overall net power output from 13.9 MW to 18.8 MW; this is because of the immense increase in the power output from the Brayton cycle which nearly doubles in value.
Figure 3 explains the impact of the variation of the pressure ratio on the efficiencies of the Brayton cycle, the combined cycle, and the overall system. The effects observed in
Figure 3 mirror those presented in
Figure 2. For an increase in the pressure ratio from 3 to 9, the Brayton cycle efficiency increases from 24% to 42%. Moreover, the efficiency of the combined power cycles increases from 49% to 66%. It should be noted that combining a simple Rankine cycle to recover the exhaust heat from the Brayton cycle increases the combined cycle efficiency by about 25% compared to the Brayton cycle alone. The efficiency of the overall system also increases from 56% to 73%.
Figure 4 examines the effect of varying the inlet pressure of the steam turbine on the steam turbine performance and the overall combined system performance. It can be observed that increasing the steam turbine inlet pressure from 4 MPa to 16 MPa leads to an increment in the net power output of the Rankine cycle from 6.5 MW to 7.7 MW. The net power output of the total system increases from 15.4 MW to 16 MW. This increase in the steam turbine inlet pressure also increases both the efficiency of the Rankine cycle and the efficiency of the overall system. The steam cycle efficiency increases from 31% to 38% while the overall efficiency increases from 62% to 66%.
The impact of varying the pressure ratio of the heat pump on its compressor work, the heat pump output, and the COP of the heat pump is illustrated in
Figure 5. An increase in the pressure ratio from 1.8 to 3 causes an increase in all the parameters considered. The heat pump compressor work increases from 0.13 MW to 0.29 MW, and the heat pump output increases from 1.37 MW to 1.53 MW. The COP of the heat pump, however, decreases from 11 to 5.3. The resultant changes in the Brayton cycle power output, the combined power output, the net power output of the overall system, and the heat pump output energy are shown in
Figure 6 when the combustion inlet energy rate is varied. It can be observed that increasing the combustion inlet energy rate from 25 MW to 35 MW causes an increase in the output energy rate of all the systems. The heat pump output energy rate increases from 1.38 MW to 1.39 MW, while the Brayton cycle experiences an increase in power output from 8.2 MW to 11.5 MW. The overall power output from the multigeneration system increases from 15.5 MW to 19 MW.
The effect of varying the inlet energy rate from the combustion chamber on the
and
is presented in
Figure 7. While the
decreases from about 0.51 to 0.3, the
can be seen to decrease from 2.05 to 0.6. The decreasing nature of both factors is a good indicator of the system’s reduced environmental impact, and this also shows that the system’s irreversibility is decreasing.
Figure 8 shows the effect of varying the inlet energy rate from the combustion chamber on
and
of the system. The
decreases from about 1.0 to 0.1, while the
increases from 1.05 to 6.1. The value of
for a system should be as low as possible, while a high
indicates that the system is beneficial to the environment. The effect of increasing inlet energy rate on the exergy stability factor and exergetic sustainability index of the system can be seen in
Figure 9. The
increases from 0.41 to 0.47, while the
also increases from about 0.6 to 0.64. The factors increasing point to the conclusion about the stability and sustainability of the system upon increasing the inlet combustion energy rate.
The exergy destruction rate of each component of the system as a percentage of the overall exergy destruction rate is displayed in
Table 4. It is seen that over 50% of the exergy destruction in the system occurs in the combustion chamber; this is reflective of systems which operate on fuel combustion of any kind. The large range of temperature difference of the chemical and heat transfer between the working fluid and the methane fuel causes high levels of exergy destruction in the combustion chamber [
32]; improving this multigeneration system significantly will require the fuel combustion process to be thoroughly improved.
4. Conclusions
In this paper, a combined cycle power plant was analyzed to observe the performance of a combined gas and steam power plant under varied component parameters. The proposed system was modelled using the Thermolib extension library in MATLAB. The results obtained from the system gave the net power from the Brayton cycle, Rankine cycle and overall system as 7.41 MW, 29.91 MW, and 24.75 MW respectively. The efficiencies were found to be 33%, 36%, and 62% respectively. The overall exergy efficiency of the multigeneration system was calculated as 78%. The ammonia heat pump had an output of 1.39 MW with a COP of 9.1.
The outcome of the study showed that using different pressure ratios ultimately affected the gas cycle and the overall system performance, as it increased the gas turbine compressor work, gas turbine output power, the net power output, and the overall efficiency of the system. Increasing the pressure ratio of the Brayton cycle from 3 to 9 increased the Brayton cycle output power from 7.2 MW to 12.3 MW, while the net power output increased from 13.9 MW to 18.8 MW. Varying the pressure ratio also increased the efficiency of the Brayton cycle from 24% to 42%, while the combined efficiency of both the Brayton cycle and the Rankine cycle increased from 49% to 66%. Increasing the steam turbine inlet pressure from 4 MPa to 16 MPa increased the Rankine cycle power output from 6.5 MW to 7.7 MW, while the overall efficiency of the system increased from 62% to 66%. The heat flow rate into the combustion chamber also varied from 25 MW to 35 MW, and this caused the net power output of the system to increase from 15.5 MW to 19 MW. The exergoenvironmental analysis of the system showed the appropriateness, usefulness, and stability of the multigeneration system in the environment.
The effect of varying the combustion inlet energy rate on the exergoenvironmental impact of the overall system was also analyzed. The results showed favorable exergoenvironmental effects, indicating a reduction in the system’s overall irreversibility and environmental impact, made possible by the improved combustion efficiency in the combustion chamber.