Simulation of Power Plants by Excel
Simulation of Power Plants by Excel
Simulation of Power Plants by Excel
University of Tripoli
g.fellah@uot.edu.ly
Abstract
The aim of this work is to demonstrate the power of Excel software in simulating
the performance of thermal systems. A steam power plant is taken as an example for the
demonstration. The idea is to show how one can analyze and simulate the performance of
thermal systems within the Excel environment.
The thermodynamic properties are obtained by the Excel tools developed for
thermodynamics. The platform of these tools is the Microsoft Excel. The obtained
properties using these tools are tested and found in good agreements with other sources.
Energy balance is performed on each heater to find the extracted mass flow rates.
The obtained equations are solved simultaneously by using Excel software. The functions
MMULT and MINVERSE are used to multiply the matrices and to find the inverse and
hence the mass flow rates. The solver facility is used to optimize the thermal efficiency
and/or the plant irreversibility.
2
Nomenclature
cP specific heat [kJ/kg.K] i inlet
h enthalpy [kJ/kg] k component
İ irreversibility rate [kW] o ambient
𝑚̇ mass flow rate rate [kg/s] P product
P pressure [kPa]
𝑄̇ heat transfer [kW] Abbr.
s entropy [kJ/kg.K] BFP boiler feed pump
T temperature [o C] CEP condensate pump
w work done kJ/kg Cond condenser
𝑊̇ power [kW] CW cooling water
x mass fraction [-] DCA drain cooler approach [o C]
DEA deaerator
Greek HPH high pressure heater
first law efficiency [-] HPT high pressure turbine
effectiveness [-] IPT intermediate pressure turbine
Ψ̇ exergy rate [kW[ LHV lower heating value [kJ/kg]
specific exergy [kJ/kg] LPH low pressure heater
LPT low pressure turbine
Subscript RH reheater
d destroyed RHS right hand side
e exit SG steam generator
F fuel TTD terminal temperature difference [o C]
1. Introduction
Spreadsheets programs are used to analyze many problems in different engineering areas.
They offer an attractive alternative to conventional programming that allows ready
experimentation with numerical algorithms [1].
The design and development of a Microsoft Excel-based Power System Load Flow
Analysis tool and its application for system planning and operation were demonstrated by
Musti and Ramkhelawan [4]. They developed a simple desktop tool which provides an
interactive and simplified interface for users to store different systems with different
operating conditions and then to observe the response of the system.
3
Musti [5] presented the design and development of a Microsoft Excel-based tool for Power
System Static State Estimation. The tool can be effectively used to understand the process
of state estimation and its real-time application. It contains Newton-Raphson load flow that
provides system measurements, which are used as inputs to the state estimator that uses the
popular weighted least square (WLS) algorithm.
El-Awad presented [6] an Add-Ins that determines the thermodynamic properties for
various fluids by using Microsoft Excel. The Add-Ins provides property functions for ideal
gases, saturated and superheated water, saturated and superheated refrigerants for vapor-
compression (VC) systems, binary solutions of ammonia-water and water-lithium bromide
for vapor-absorption (VA) refrigeration systems, reacting mixtures, and atmospheric air.
Excel-Thermax platform for performing energy and exergy analyses of the evaporative
regenerative gas-turbine (ERGT) cycle was introduced by El-Awad [7].
A spreadsheet was developed to present a model that uses the effectiveness-NTU method
to explicitly take into consideration the design particulars of the regenerator, such as its
size and overall heat-transfer coefficient was presented by El-Awad [8].
A paper deals with the use of Microsoft Excel as an educational tool for conducting basic
engineering analyses related to thermal-fluid systems was presented by El-Awad [9]. The
paper focuses on using Excel and its Goal-Seek command for solving thermal-fluid
problems that require iterative solutions by presenting three related examples from the
subjects of heat-transfer, fluid dynamics, and thermodynamics.
The thermodynamic properties of water are not included in Microsoft Excel. Those
properties which are needed to perform the thermodynamic analysis can be inserted by
using Add-Ins files. The following steps should be taken to implement the Add-Ins file of
thermodynamic properties of water:
4
To make sure that it works, for example, find the enthalpy for water at given
pressure and temperature, just write “= h”, choose h_PT_H2O and set the pressure
and temperature in order, and you will have the enthalpy, see Figures 1a and 1b.
All other properties can be generated for any arbitrary two independent properties.
The steam power which will be simulated within the Excel software environment
is shown in Figure 2. All properties will be generated by the Add-Ins fluid package
[10]. The seven linear equation which is obtained by taking energy balance on the
seven heaters will be solved by using Excel facilities to obtain the mass fraction
which is extracted from the steam turbine.
5
Figure 2: Schematic of the steam power plant
6
3.1. Thermodynamic modeling
For the analysis, steady-state, steady flow processes are assumed. Pressure drop due
to friction, heat exchange with surroundings, the change in kinetic and potential energies
are neglected.
𝑊̇ (𝑀𝑊)
𝑚̇𝑐𝑦𝑐𝑙𝑒 = (2𝑏)
𝑘𝐽
𝑤( )
𝑘𝑔
𝐸𝑛𝑒𝑟𝑔𝑦 (𝑠𝑜𝑢𝑔ℎ𝑡)
𝜂= (3)
𝐸𝑛𝑒𝑟𝑔𝑦 (𝑐𝑜𝑠𝑡)
Exergy is defined as the maximum useful reversible work that can be obtained from a given
mass in a given state when the mass brought reversibly into thermodynamic equilibrium
with the environment.
Ψ̇ = 𝑚̇[(ℎ − ℎ0 ) − 𝑇0 (𝑠 − 𝑠0 )] (4)
And the exergy balance for a given component can be written as;
𝑁 𝑁 𝑁
𝑇0
∑ (1 − ) 𝑄̇𝑘 + ∑ Ψ̇𝑖,𝑘 = ∑ Ψ̇𝑒,𝑘 + 𝑊̇𝑘 + 𝐼𝑘̇ (5)
𝑇
𝑖 𝑖 𝑒
By using the definitions of Fuel-Product-Loss (F-P-L). Fuel and Product are expressed by
exergy flow. Exergy balance for a single component (k) is given as:
7
𝛹̇̇̇ 𝐹 = 𝛹̇̇̇ 𝑃 + 𝛹̇̇̇ 𝐷 (6)
Where𝛹̇̇̇ 𝐹 , 𝛹̇̇̇ 𝑃 𝑎𝑛𝑑 𝛹̇̇̇ 𝐷 are the input exergy (fuel), exergy rate of the desired product, and
the exergy destroyed (irreversibility) during the process, respectively.
Ψ̇𝑃 Ψ̇𝐷
𝜀𝑘 = =1− (7)
Ψ̇𝐹 Ψ̇𝐹
The definitions of F-P (Fuel-Product) for the current power unit are given in Table 2. The
effectiveness of the power cycle is given as:
𝑊̇𝑛𝑒𝑡
𝜀= (8)
𝑚̇𝑓𝑢𝑒𝑙 × 𝜓𝑓𝑢𝑒𝑙
𝐿𝐻𝑉
≈ 1.00565 (9)
𝜓𝑓𝑢𝑒𝑙
HPT Ψ̇11 − Ψ̇12 − Ψ̇13 𝑊̇𝐻𝑃𝑇 LPH1 Ψ̇18 − Ψ̇27 Ψ̇6 − Ψ̇5
IPT Ψ̇14 − Ψ̇15 − Ψ̇16 − Ψ̇17 𝑊̇𝐼𝑃𝑇 LPH2 Ψ̇19 + Ψ̇28 − Ψ̇29 Ψ̇5 − Ψ̇4
LPT Ψ̇17 − Ψ̇18 − Ψ̇19 − Ψ̇20 − Ψ̇21 − Ψ̇22 𝑊̇𝐿𝑃𝑇 LPH3 Ψ̇20 + Ψ̇30 − Ψ̇31 Ψ̇4 − Ψ̇3
Cond Ψ̇22 − Ψ̇34 − Ψ̇1 Ψ̇𝑐𝑤2 − Ψ̇𝑐𝑤1 LPH4 Ψ̇21 + Ψ̇32 − Ψ̇33 Ψ̇3 − Ψ̇2
8
3.2. Steps of the analysis
The main steps of the analysis by using the Excel spreadsheet is summarized in the
following steps:
(i) Step 1: Assign the properties values for the given states.
(ii) Step 2: Use your knowledge in thermodynamic to generate the properties,
such as the enthalpy, entropy, and specific volume when necessary as shown
in Table 3.
(iii) Step 3: To find the mass fraction for each heater, solve the set of the linear equations
which are developed by taking energy balance on the feed water heaters (Eq. (2a)).
Use MINVERSE to create the inverse matrix, and MMULT to multiply the two
matrices. Here, the generated inverse matrix is multiplied by the original one to
obtain the unknowns (the mass fractions).
(iv) Step 4: Calculate the exergy at the all state points by using Eq. (4).
(v) Step 5: Calculate the irreversibility (Eq. (6)), the first law efficiency (Eq. (3))
and the effectiveness (Eq. (8)).
The results could be optimized by selecting the objective function, which might be the
efficiency (maximizing) or the plant irreversibility (minimizing). Also, the constraints
must be selected, which might be the bleeding pressures by letting them vary in a certain
range, for example, +/- 10 %.
9
(i) Step 1: Select data from the menu.
(ii) Step 2: select solver, see Fig. 3. Solver Add-Ins must be activated through:
𝑂𝑝𝑡𝑖𝑜𝑛𝑠 → 𝑎𝑑𝑑 𝑖𝑛𝑠 → 𝑔𝑜 (𝑒𝑥𝑐𝑒𝑙 𝑎𝑑𝑑 𝑖𝑛𝑠) → 𝑐ℎ𝑒𝑐𝑘 𝑡ℎ𝑒 𝑠𝑜𝑙𝑣𝑒𝑟 𝑎𝑑𝑑 𝑖𝑛𝑠 𝑏𝑜𝑥
(iii) Select maximize (for thermal efficiency), or minimize for the irreversibility.
(iv) Select the objective function cell (efficiency or irreversibility).
(v) Select the constraints cells and set the limits.
(vi) Select the optimization method.
(vii) Push on the solve button to implement the optimization.
10
4. Results and discussions
STATE P(kPa) T(oC) h(kJ/kg) s(kJ/kg.K) STATE P(kPa) T(oC) h(kJ/kg) s(kJ/kg.K)
1 7 39 163.3655 0.5591 19 298 209 2883.997 7.3540
2 977 39 164.6685 0.5603 20 161 153 2777.243 7.3987
3 977 90 375.7926 1.1868 21 69 90 2650.336 7.4605
4 977 113 476.7254 1.4563 22 7 39 2367.544 7.6204
5 977 133 560.9334 1.6687 23 4350 216 925.5933 2.4761
6 977 155 655.8304 1.8960 24 2030 213 925.5933 2.4819
7 977 179 758.2299 2.1286 25 2030 188 799.9385 2.2175
8 17490 183 783.0186 2.1426 26 977 179 799.9385 2.2209
9 17490 210 905.063 2.4024 27 549 139 584.5979 1.7278
10 17490 253 1098.458 2.7857 28 298 133 584.5979 1.7287
11 17490 538 3389.603 6.3851 29 298 119 499.9804 1.5180
12 4350 338 3053.422 6.4841 30 161 113 499.9804 1.5187
13 4350 338 3053.422 6.4841 31 161 95 398.7221 1.2519
14 4350 538 3529.322 7.1597 32 69 90 398.7221 1.2526
15 2030 429 3312.014 7.2150 33 69 45 187.3257 0.6351
16 977 337 3130.352 7.2682 34 7 39 187.3257 0.6358
17 977 337 3130.352 7.2682 cw-in 101.325 15 63.07903 0.2245
18 549 271 3003.61 7.3097 cw_out 101.325 32 134.1932 0.4642
Solving the seven linear equations which are developed by taking mass and energy
balance on the feed water heaters (use Eqs. (1 &2a)) will result in the mass fractions (x12,
x15, x16, x18, x19, x20 and x21) which are extracted for each heater. The linear equations are
set in a matrix form. The functions MMULT and MINVERSE are used to multiply the
matrices and to find the inverse, see Tables. 5a and 5b.
11
Table 5b: The inverse matrix
The heat transfer and work done per kilogram of the steam are calculated by using the
first law of thermodynamics (Eq. (2a)). Those are, heat transfer to the steam generator
(qh), heat transfer in the condenser (qc), work done by the high pressure turbine (whp),
the work done by the intermediate pressure turbine (wip), work done by the low pressure
turbine (wlp), work done by the condensate pump (wcp) and the work done by the feed
pump (wfp). The net work done (wnet) is then calculated and tabulated in Table 7.
qh 2723.791 kJ/kg
qc -1557.86 kJ/kg
whp 336.181 kJ/kg
wip 354.7077 kJ/kg
wlp 567.7104 kJ/kg
wCP -1.08353 kJ/kg
wfp -24.7887 kJ/kg
wnet 1232.727 kJ/kg
12
By knowing the net power output (350 MW), Eq (2b) is applied to calculate the cycle
mass flow rate (283.9234 kg/s).
All other mass flow rates are calculated by multiplying the mass fractions by the cycle
mass flow rate, the results are tabulated in Table 8.
The rate of heat transfer and power can be calculated and tabulated by multiplying the
specific quantities from Table 7 by the cycle mass flow rate, the results in MW are shown
in Table 9.
Qh 773.3480 [MW]
Qc -442.3136 [MW]
Whp 95.4496 [MW]
WIP 100.7098 [MW]
WLp 161.1863 [MW]
WCP -0.3076 [MW]
WFP -7.0380 [MW]
Wnet 350 [MW]
The thermal efficiency (Eq. (3)), fuel exergy (Eq. (9)), effectiveness (Eq. (8)) and all other
heat transfer (Eq. (2a)) is now calculated and tabulated in Table 10.
13
Table 10: Selected results
Exergy at each state is calculated (Eq. (4)) and tabulated in Table 11.
The total irreversibility of the steam cycle is found (Eq. (6)) equals to 554.709 MW as
seen in Table 12.
14
Table 12: Irreversibilities of the cycle components.
Optimization:
The next step of the analysis is an attempt to reduce the cycle irreversibility. For
optimization:
(i) Select data from the menu and then go to the solver.
(ii) Select the objective function (cell contains the total irreversibility).
(iii) Select minimize (for the irreversibility).
(iv) Select the constraints cells as shown in Table 13, and set the limits as shown in
Table 14.
State P [kPa]
11 17490
12 4350
18 549
19 298
20 161
21 69
15
Table 14: Setting the limits
The new property values are now generated as shown in Table 15.
16
Table 15: Property values after optimization
Selected results are shown in Table 16. As it can be seen both the heat flow rate (and
hence the rate of the fuel consumption) and the rate of the input exergy are reduced to
902.8341 MW and 897.752 MW, respectively. A minor increase in the thermal efficiency
is noticed.
17
Table 17, shows the exergy at each state after optimization, and Table 18, shows the
irreversibility for each component and hence the total irreversibility is calculated.
As it can be seen, the total irreversibility is reduced to 547.7519 with a small increase in
the effectiveness
5. Conclusions
REFERENCES
[10] http://www.me.ua.edu/ExcelinME/index.htm.
19