Simulator Laboratory

Simulator Laboratory

Sim Labs for Thermodynamics and Thermal Power Plant Simulator

Serhat Beyenir & Sanja Boskovic

BCIT/BCcampus

Vancouver, B.C.

1

Dedication

To our students who taught us to be better teachers than we were many years ago.

Contents

2

About the Book

These Simulator Laboratory Worksheets are open educational resources that align with BCIT’s Power Engineering Thermodynamics and Thermal Power Plant Simulator courses.

Open educational resources (OER) are defined as “teaching, learning, and research resources that reside in the public domain or have been released under an intellectual property license that permits their free use and re-purposing by others.” (Hewlett Foundation).

Funding for this project was provided through an OER Grant from the B.C. Open Textbook Project. The B.C. Open Textbook Project began in 2012 with the goal of making post-secondary education in British Columbia more accessible by reducing student cost through the use of openly licensed textbooks and other OER. The B.C. Open Textbook Project is administered by BCcampus and funded by the British Columbia Ministry of Advanced Education, Skills and Training and the Hewlett Foundation.

The BCcampus Writing Guidelines for Style and Tone were referenced during the copy editing and proof reading phases of this project and the Open Textbooks Theme by Brad Payne was used in producing the book.

3

Credits

Cover image by Gretar Ívarsson is found in the public domain.

Images of simulator user interface courtesy of Kongsberg Digital Simulation. We wish to thank Kongsberg for their continuous support.

4

Preface

There are certain theoretical concepts that power engineering students struggle with and there are limited number of textbooks or videos that address these concepts directly linked to the field of power engineering. When we present those challenging concepts in class we also want to make sure that our students develop safe working habits and environmental sensibilities related to sustainable operation of the power engineering equipment.

The primary purpose of this project was to design and develop simulator lab worksheets (SIM LABS) that would be used as ancillary resources for Thermodynamics and Thermal Power Plant Simulator (TPPS) courses. We have limited time to teach theory and we also have to offer practical skills to our students, so having as many SIM LABS as possible can bring those challenging concepts closer to our students and enhance their hands-on skills.

We consider this project as a first step to writing a comprehensive open textbook in the future. Through this project we aimed to

A final note to instructors who teach similar courses: We will be glad to share the “Instructor Guides and Solutions” to all SIM LABS included in this booklet. Please email us at
Sanja_Boskovic@bcit.ca or Serhat_Beyenir@bcit.ca

5

Acknowledgments

We owe many thanks to Dr. David Porter (former Associate Vice President) for his leadership in open education during his short tenure at BCIT albeit very impactful. We are equally indebted to BCcampus and openBCIT for co-sponsoring our project.

This was a rare project in which we did not need to disentangle any bureaucratic red tape and that was due to rare people who were involved since day one of our work. Special thanks go to James Rout (Associate Vice President), Ken Jeffery and Lin Brander (co-chairs of the OpenEd Working Group) as well as Alison Griffin, Deirdre Grace, Holly Stanley and Rosario Passos in alphabetical order. These individuals kept us motivated and we truly appreciate their friendliness in equal measures.

Also thanks to those whom we never met in person but received support within Open Grants.

And last but not least, we would like to thank our students, the entire Class of 2017 for their hard work, cooperation and feedback: Adrian Pontellini, Amar Manhas, April Derald Burias, Asmamaw Gebremariam, Brodie Pulles, Chelsea  Chiao Shiun Lo, Cole Tewinkel, Darrick Cheung, Devin Hughes, Harjit Kang, Jesiah de Gouveia, John Aldana, Joshua Lewis, Kevin Hamann, Matthew Fermill, Nathan Kwong, Riley Walker, Ryan Ngun, Sean McGilp, Terence Ka Chun Pang, Tommy Yoo, Trevor Rigby, Tylar Hodson and Vangeli Asimakis.

1

Boiler Efficiency

Learning Objectives

Operate the Plant in the following modes and compute the Boiler Thermal Efficiency

  • 230 MW Burning HFO
  • 80% Load Burning Bio Fuel
  • 80% Load Burning Coal
  • 80% Load Burning Coal with soot built up.

Theory

First we define efficiency. There are many ways to present it, for example, efficiency is the ability to do something without wasting energy, effort or time. Put it mathematically, it is the ratio of useful output to thermal energy input. Layperson often uses the terms efficiency and effectiveness interchangeably however efficiency is to do with minimizing waste and effectiveness to do with maximizing output (more on this in Heat Exchangers Lab). Boiler efficiency is sometimes defined as combustion efficiency which is computed by the ratio of the burner’s capability to burn fuel completely to the unburnt fuel and excess air in the exhaust. Thermal efficiency on the other hand, indicates the heat exchanger’s (i.e. boiler’s) capacity to transfer heat from the combustion process to the water or steam in the boiler. In general, the maximum boiler efficiency attainable from a boiler depends on such factors as method of burning the fuel, design of the furnace and heat transfer surfaces. In addition, the type of fuel, boiler load and operational practices influence the boiler efficiency. In this lab, we focus on the fuel type, boiler load and best practices.

Boiler Thermal Efficiency

Next we express the Boiler Thermal Efficiency as follows:

\eta_{boiler} = \frac{Energy\ to\ steam}{Energy\ from\ fuel}

Where energy to steam is the heat transfer required to form steam. Let,

h2=specific enthalpy of steam formed, [kJ/kg],

h1=specific enthalpy of feed water, [kJ/kg].

Because the steam is formed at constant pressure, heat transfer required to form 1 kg of steam in the boiler is

Energy to Steam= (h2-h1) [kJ]

Energy from fuel is computed from the mass of fuel used and its heat (calorific) value. For coal this is the heating value as measured in a bomb calorimeter and it corresponds to the internal energy of combustion. If

mf= mass of fuel burned in a given time

ms= mass of steam generated in the same time

HV= heat value of fuel [kJ/kg]

Energy from fuel= mf x HV [kJ]

Then we can write:

\eta_{boiler} = \frac{m_{s}(h_{2}-h_{1})}{m_{f}HV}100\%

Lab Instructions

You will run 4 different initial conditions in this lab:

  • I10 230 MW Burning HFO
  • I15 80% Load Burning Bio Fuel
  • I14 80% Load Burning Coal
  • I14 80% Load Burning Coal and use MD250 to set up soot variables.

For each condition collect the relevant data to compute the Boiler Thermal Efficiency.

Hints & Tips

For data collection, use trends as shown below:

Trends sample

Trends sample

In addition to various pressure, temperature and flow values, you will need to log the following tags in your trends:

To calculate the enthalpy values, you may use an app or online tool such as the Superheated Steam Table: https://goo.gl/GdVM4U

For the coal operation with soot built up, use MD250 and set malfunctions as follows:

MD250 malfunction settings.

MD250 malfunction settings.

Deliverables

Your lab report is to include the following:

  • Trend plots: Supply all plots taken for each of the 4 conditions,
  • Computation: Use MATLAB or MS Excel and calculate the boiler thermal efficiency for the 4 conditions specified,
  • Conclusion: Write a summary (max. 500 words, in a text box if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Fundamentals of Classical Thermodynamics SI Version by G. J. Van Wylen and R. E. Sonntag: Evaluation of Actual Combustion Processes.
  • Thermal Engineering by H.L. Solberg, O.C. Cromer and A.R. Spalding: Capacity and Efficiency of Steam Generating Units.
  • Basic Engineering Thermodynamics in SI Units by R. Joel: Boiler calculations.

2

Turbine Efficiency

Learning Objectives

Operate the Plant at the following generating capacities to compute the isentropic change in enthalpy and thermal efficiency for the HP turbine.

  • 35% Load (I13),
  • 80% Load (I14),
  • 230 MW (I10).

Theory

Recall from the First and Second Law of Thermodynamics that the adiabatic process where entropy remains constant provides the maximum energy for work. As shown on the H-S coordinates, the difference in enthalpy, (H1-H2), is maximum when the lowest enthalpy (H2) is reached at the exit conditions. The ideal expansion is, therefore, a vertical line.

Turbine efficiency.

Turbine efficiency.

On the diagram above, T1, P1 and P2 are known process variables, for example, H1 is determined by using T1 and P1. H2 then can be found drawing a vertical line from P1 to P2 by following adiabatic isentropic expansion (expansion at constant entropy).

Non-ideal processes or real processes, however, do not present straight lines as shown on the Mollier diagram due to such factors as friction. If the expansion is not isentropic (i.e. entropy is not constant but it increases), the lowest enthalpy (H2) cannot be reached at the exit conditions, in other words, H2’ > H2. This means that ΔH for the ideal expansion is greater than ΔH for the non-ideal expansion between the same pressure boundaries. The internal turbine efficiency is therefore given by:

\eta_{Turbine} = \frac{Actual\ change\ in\ enthalpy}{Isentropic\ change\ in\ enthalpy}

\eta_{Turbine} = \frac{(H_{1}-H_{2\prime})}{(H_{1}-H_{2})}

The difference in enthalpy H2’-H2 is called the reheat factor and is the basis for multi-stage turbines. As can be seen on the Mollier diagram, the pressure curves are divergent. This means that the higher the pressure drop in a single stage turbine the greater the reheat factor and in turn the lower the turbine efficiency. However, if the steam is expanded through multiple stages and between each stage the steam is reheated, higher turbine efficiencies can be achieved. We will see this effect later in the Power Plant Efficiency lab.

Lab Instructions

You will run 3 different initial conditions in this lab:

  • 35% Load (I13),
  • 80% Load (I14),
  • 230 MW (I10).

For each condition collect the relevant data to compute the isentropic change in enthalpy for the HP turbine. Compare your results, which of the three conditions yield the most favourable results and why?

Hints & Tips

In addition to various pressure and temperature values; log the following tags in your trends:

To calculate the enthalpy values, you may use an app or online tool such as the Superheated Steam Table: https://goo.gl/GdVM4U

Deliverables

Your lab report is to include the following:

  • Trend plots: Supply all plots taken for each of the 3 conditions,
  • Computation: Use MATLAB or MS Excel and calculate the turbine efficiency for the 3 conditions specified,
  • Conclusion: Write a summary (max. 500 words, in a text box if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Thermodynamics and Heat Power by I. Granet: Vapor power cycles.

3

Power Plant Efficiency

Learning Objectives

Operate the Plant at full generating capacity and compute the Power Plant Efficiency when the plant is operating:

  • Under normal conditions,
  • With the cooling water temperature very high (lake water temperature: 35°C),
  • Without regeneration.

Theory

Excluding hydroelectric power plants, most power generating plants employ a type of boiler and steam turbine. A schematic diagram of a simple steam power plant is shown below:

Schematic diagram of a steam power plant

Schematic diagram of a steam power plant

High-pressure steam leaves the boiler and enters the turbine. The steam expands in the turbine and does work which enables the turbine to drive the electric generator. The exhaust steam leaves the turbine and enters the condenser where heat is transferred from the steam to cooling water. The pressure of the condensate leaving the condenser is increased in the pump thereby enabling the condensate to flow into the boiler. This thermodynamic cycle is known as the Rankine Cycle.

The Rankine Cycle Efficiency

As noted above, some heat is always lost from the steam to cooling water. In addition, feed pumps consume energy thus reducing the net work output. Rankine Cycle Efficiency then can be expressed as:

\eta_{Rankine} = \frac{Net\ work\ output}{Heat\ supplied\ in\ the\ boiler}

or

\eta_{Rankine} = \frac{W_{Turbine}-W_{Pump}}{Q_{boiler}}

referring to the diagram above and using the enthalpy values in the Rankine cycle, we can write:

\eta_{Rankine} = \frac{(h_{1}-h_{2})-(h_{4}-h_{3})}{(h_{1}-h_{4})}

Improvements to the Rankine Cycle Efficiency

Effect of Pressure and Temperature on the Rankine Cycle

If the exhaust pressure drops from P4 to P4‘ with the corresponding decrease in temperature at which heat is rejected in the condenser the net work is increased by area 1-4-4′-1′-2’-2-1 (see diagram below)

Effect of exhaust pressure

Effect of exhaust pressure

In a similar way, if the steam is superheated in the boiler, it is evident that the work is increased by area 3-3′-4′-4-3 (see diagram below):

Effect of superheating

Effect of superheating

Superheating the steam is done by increasing the time the steam is exposed to the flue gases. The result of superheating is that for a given power output, the plant using superheated steam will be of smaller size than that using dry saturated steam.

The Reheat Cycle

Above we noted that the efficiency of the Rankine cycle is increased by superheating the steam. If metals could be found that would allow us to reach higher temperatures, the Rankine cycle could be more efficient. To improve the efficiency, the reheat cycle has been developed which is shown schematically below:

Rankine Cycle with Reheat

Rankine Cycle with Reheat

In this cycle, the steam is expanded to some intermediate pressure in the turbine and is then reheated in the boiler, after which it expands in the low-pressure turbine to the exhaust pressure. Rankine Cycle with reheat thermal efficiency can be expressed as:

\eta_{thermal} = \frac{W_{12}+W_{67}-W_{43}}{Q_{41}+Q_{26}}

The Regenerative Cycle

Another variation from the Rankine cycle is the regenerative cycle, which involves the use of feedwater heaters. During the process between states 2 and 2′ the feedwater is heated and the average temperature is much lower during this process than during the vaporization process 2′-3. In other words, the average temperature at which heat is supplied in the Rankine cycle is lower than in the Carnot cycle 1′-2′-3-4-1′, and consequently the efficiency of the Rankine cycle is less than that of the corresponding Carnot cycle. The relationship between Carnot cycle and Rankine cycle is shown below.

Relationship between Carnot cycle and Rankine cycle

Relationship between Carnot cycle and Rankine cycle

In the regenerative cycle, feedwater enters the boiler at some point between 2 and 2′. As a result, the average temperature at which heat is supplied is increased. A schematic of practical cycle is shown below:

Regenerative cycle

Regenerative cycle

The Plant Thermal Efficiency

In order to calculate the overall plant thermal efficiency, we need to adjust the formulas above to incorporate heat added in the reheater sections of the boiler:

\eta_{thermal} = \frac{W_{Turbines}-W_{Pumps}}{Q_{boiler}+Q_{Reheat1}+Q_{Reheat2}}

Lab Instructions

Run the initial condition I10 230 MW_oil_auto:

  • Draw a T-S diagram of the Rankine cycle (not to scale) including reheat and regeneration,
  • Using Trend Group Directory, collect the relevant process values,
  • Calculate the overall thermal efficiency of the plant:
    • Under normal conditions,
    • When the cooling water temperature is very high (Set the Variable List Page 0100, tag#: T00305 to 35°C),
    • When all the steam extraction valves are closed (i.e. no regeneration and T00305 set to 10°C).

Hints & Tips

In this lab, you are essentially calculating the Rankine Cycle thermal efficiency. However, you need to take the reheat cycle into consideration and log the following tags in your trends:

For Boiler Feedwater Inlet Temperature, you may use the Startup Heat Exchanger Feedwater Outlet Temperature tag#: T02447.

For the second calculation, locate the Variable List Page 0100 as shown below:

Lake water temperature setting

Lake water temperature setting

For the third calculation, make sure you closed all steam extraction valves and set T00305 to 10°C:

No steam extraction

No steam extraction

To calculate the enthalpy values, you may use an app or online tool such as the Superheated Steam Table: https://goo.gl/GdVM4U

Deliverables

Your lab report is to include the following:

  • T-S diagram: As per instructions above,
  • Trend plots: Supply all plots taken for this lab,
  • Computation: Use MATLAB or MS Excel and calculate the overall thermal efficiency of the plant as per Lab Instructions.
  • Conclusion: Write a summary (max. 500 words, in a text box if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Applied Thermodynamics for Engineering Technologists by T. D. Eastop and A. McConkey: Steam Plant.
  • Fundamentals of Classical Thermodynamics SI Version by G. J. Van Wylen and R. E. Sonntag: Vapor power cycles.
  • Thermodynamics and Heat Power by I. Granet: Vapor power cycles.

4

Emissions of NOx

Learning Objectives

Operate the Plant at full generating capacity and compare the NOx emissions when the plant is operating:

  • Under normal conditions,
  • During over fire air damper failure,
  • During over burner air control failure,
  • Burning poor quality fuel,
  • With the DeNOx plant bypassed.

Theory

Combustion of coal generates considerable quantities of byproducts, some of which are considered pollutants. The byproducts are mostly water vapor which is what we see coming out of a power plant smokestack, carbon dioxide, and nitrogen that is readily available in the air we breathe, and they do not necessarily pose any direct health hazard. However, the emissions do carry small concentrations of pollutants into the atmosphere, which translate into large quantities of hazardous emissions due to the large amount of coal combusted. The main pollutants that can cause health problems are sulfur oxides, nitrogen oxides, particulate matter (see Combustion Analysis) and such trace elements as arsenic, lead and mercury.

During the combustion process in a coal-fired power plant, nitrogen from the coal and air is converted into nitric oxide (NO) and nitrogen dioxide (NO2); these nitrogen oxides are commonly known as NOx. NOx emissions contribute to the formation of acid rain.

NOx is primarily formed by two mechanisms: thermal NOx and fuel-bound NOx.

Thermal NOx formation takes place at high flame temperatures. Formation of thermal NOx increases exponentially with combustion temperature. Fuel-bound NOx formation is dependent upon the nitrogen content of the fuel.

The best way to minimize NOx formation is to reduce flame temperature, reduce excess oxygen, and/or to burn low nitrogen-containing fuels.

The DeNOx Plant Description

The purpose of the DeNOx plant is for removal of nitrogen oxide from the flue gases. The plant employs a selective catalytic reduction method. The medium used for the reduction is ammonia gas.

The DeNOx plant includes two selective catalytic reduction (SCR) reactors and an ash silo. Various dampers channel the flue gases either into or bypass the SCR-reactors.

DeNOxPlantOverview

DeNOx Plant Overview

Operation of the DeNOx Plant

The DeNOx plant is highly automatic, controlled by Programmable Logic Control sequences. The sequences are S701-Purge, S702/703-Start/stop Reactors, S704/705-Heating of Reactors, S706/707-Ammonia Injection, S708/709-Product Handling and S710/711-Soot Blowing.

Over Fire Air

To reduce NOx, a quantity of additional air above all burner planes is supplied. This over-fire-air (OFA) reduces NOx by enabling richer fuel mix in the high-temperature combustion zone at the burners (two-step combustion). Approximately 10 % of the combustion air is added as OFA air. OFA controller failure (MD200 malfunction 0881) is modeled in the simulator:

OFA controller failure.

OFA controller failure.

Over Burner Air

To further reduce NOx generation, a portion of the secondary air is split off the main duct and directed into a third channel just above the burner. The damper controlling this over-burner-air is abbreviated OBA. OBA controller failure (MD180 malfunction 0780) is modeled in the simulator:

OBA controller failure.

OBA controller failure.

Fuel Quality

In the simulator, the chemical composition of the coal or other fuels can be specified. The sum of the five components C, H, S, O and N should preferably add up to 100%, to avoid confusion, but it is not strictly necessary because the C, H, S, O and N setting is always recalculated to a 100% basis prior to using in other computations. Water and inert matter (ash/slag) should then be added. The water content varies much and has a great impact on the amount of preheating required by primary air. The simulator computes the lower heat value (including water/inert matter) and theoretical combustion air needed and flue gas produced. The air/flue gas values are given in ncm/kg (ncm=normal cubic meter).

In this lab, we will burn both default and lower quality coal for a comparison. Fuel data can be changed using Variable List page 0111 on MD180, for example:

Chemical composition of coal

Chemical composition of coal

Lab Instructions

Run the initial condition I14 80% Coal and setup trends for the following variables:

G02197 X17821 G17107 X17106 D17104

T17103 C08444 G08444 G08443 C08400

  1. Stable operation: After 5 minutes of running a stable operation, freeze simulator and print the two trends. This is the reference point for the rest of the lab.
  2. OFA damper failure: Switch to run mode and activate malfunction 0881 on MD200. After 5 minutes, freeze simulator and print the two trends. Before moving on to the next step deactivate the malfunction.
  3. OBA control failure: Switch to run mode and activate malfunction 0780 on MD180. After 5 minutes, freeze simulator and print the two trends. Before moving on to the next step deactivate the malfunction.
  4. Burning poor quality fuel: Switch to run mode and access Variable List page 0111 on MD180. Set the new values as shown below. After 5 minutes, freeze simulator and print the two trends.
    • X00820: 70.40
    • X00821: 5.10
    • X00822: 1.10
    • X00823: 12.50
    • X00824: 1.60
    • X00825: 9.30
  5. DeNOx plant bypassed: Switch to run mode and bypass the SCR 1 and SCR 2 on MD710 and MD720 respectively. After 5 minutes, freeze simulator and print the two trends.

Hints & Tips

Your trend windows should look like the following:

 

Trend sample 1: DeNOx plant data.

Trend sample 1: DeNOx plant data.

Trend sample 2: NOx stack data.

Trend sample 2: NOx stack data.

Make sure your trend printouts are labeled properly otherwise, data analysis will be very confusing.

Tabulate your data as shown below:

The deNOx plant data.

The deNOx plant data.

Percentage deviation can be expressed as

Percentage\ deviation = \frac{Current\ Value-Reference\ Value}{Reference\ Value}100\%

and tabulated as follows:

The deNOx plant deviation data.

The deNOx plant deviation data.

Deliverables

Your lab report is to include the following:

  • Trend plots: Supply all plots taken for this lab (make sure plots are labeled properly),
  • Computation: Use MATLAB or MS Excel to process your data. Calculate the percentage deviation for each operation and plot your results,
  • Conclusion: Write a summary (max. 500 words in a text box, if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Thermal Power Plant Simulator Course Manual by BCIT: The DeNOx plant

5

Heat Exchangers

Learning Objectives

Operate the Plant in full generating capacity and study the effect of heat exchanger surface area using the Low Pressure Feed Heater 3.

Theory

A heat exchanger is equipment in which heat exchange takes place between two working media that enter and exit at different temperatures. The main function of heat exchanger is to either remove heat from a hot working media or to add heat to the cold working media. Depending on direction of working media-fluid flow the heat exchanger is either parallel (concurrent) flow heat exchanger or counter flow heat exchanger (see Figures below).

Concurrent flow.

Concurrent flow.

Counter current flow.

Counter current flow.

The terms are related to how the fluid flows through their respective flow passages relative to each other. If fluids flow in the same direction such as in Figure 1 it is termed a parallel flow. If fluids flow in opposite directions as in Figure 2 it is termed counter flow.

Parallel flow in heat exchangers occurs when both fluids enter the heat exchanger at their largest temperature difference. The temperature difference becomes less over the length of the heat exchanger. In the counter flow heat exchanger, fluids enter at opposite ends and therefore at different ends of the temperature scale Figure 2. The temperature difference between two fluids is relatively constant over the length of the exchanger.

The heat transfer process that occurs in any heat exchanger can be described by the following equations.

Q_{hot} = m_{hot}c_{p hot}\Delta T_{hot}
Q_{cold} = m_{cold}c_{p cold}\Delta T_{cold}

Considering the surface area involved in heat transfer, Newton’s Law of cooling states that the rate of heat loss is proportional to the difference in temperatures between the body and its surroundings and given by,

Q = \alpha Area \Delta T

where α is called the heat transfer coefficient [W/m2K], area is taken in m2 and ΔT is the temperature difference.

Furthermore Q, heat transferred between the hot water and cold water can be calculated as follows:

Q = \frac{F(LMTD)}{R_{T}}

Or

Q = F(UA)(LMTD)

where F is the correction factor which equals 1 for this SIMLAB (it takes values between 0.5 and 1). RT is the overall resistance, U, overall heat transfer coefficient and LMTD is the Log Mean Temperature Difference.

The overall resistances can be calculated using:

R_{T}=R_{hf}+R_{w}+R_{cf}

R_{hf} = \frac{1}{A_{1}\alpha_{h}}

R_{w} = \frac{ln\frac{D_{2}}{D_{1}}}{2\pi L\lambda_{w}}

R_{cf} = \frac{1}{A_{2}\alpha_{c}}

Heat transfer coefficients ah and ac can be calculated using the following expression for Nusselt number for hot and cold water:

For cooling

\alpha_{h} = \frac{Nu_{h}\lambda_{h}}{D_{h}}

Nu_{h} = 0.3Re_{h}^{0.8}Pr_{h}^{0.3}

 

For heating

\alpha_{c} = \frac{Nu_{c}\lambda_{c}}{D_{c}}

Nu_{c} = 0.3Re_{c}^{0.8}Pr_{c}^{0.3}

And the LMTD is given by the following correlation where 1 and 2 presents the ends of the heat exchanger:

\Delta T_{LMTD} = \frac{\Delta T_{1}-\Delta T_{2}}{ln\frac{\Delta T_{1}}{\Delta T_{2}}}

Heat Exchanger Effectiveness

Recall from the Boiler Efficiency Lab that efficiency is to do with minimizing waste and effectiveness to do with maximizing output. Here we define heat exchanger effectiveness as the ratio of actual heat transfer rate to the maximum possible heat transfer rate for the given temperatures.

\epsilon = \frac{Q}{Q_{max}}

Q_{max} = C_{min}(T_{hi}-T_{ci})

Where Cmin is defined either by cold or hot fluid whichever is smaller and it is defined by:

C_{h} = m_{h}c_{ph}

C_{c} = m_{c}c_{pc}

 

Lab Instructions

Run the initial condition I10 230 MW_oil_auto and setup trends for

Heat Area Factor C34201
Q Q34228
Tcold1 T24214
Tcold2 T34214
mcold G34213
Thot1 T34204
Thot2 T34227
mhot G34203
  1. Heat Area Factor set to 0.5: Using MD420 and Variable List 4210 set C34201 to 0.5. Run the simulator with this setting for 15 minutes. Freeze and print the two trends.
  2. Heat Area Factor set to 1: As in step 1, set C34201 to 1. This is the default setting for the Heat Area Factor (heat transfer coefficient x area). After 15 minutes of running the simulator, freeze simulator and print the two trends.
  3. Heat Area Factor set to 1.5: This time, set C34201 to 1.5. Run the simulator with this setting for 15 minutes. Freeze and print the two trends.

Hints & Tips

As always, label your trends using descriptive names. In this lab, you are changing the Heat Area Factor (heat transfer coefficient x area) of LP Feed Heater 3 and comparing the temperature data. Your evaluation will be based on LMTD and Heat Exchanger Effectiveness values.

MD420 LP Feed Heater 3

MD420 LP Feed Heater 3

Variable List 4210.

Variable List 4210.

Deliverables

Your lab report is to include the following:

  • Trend plots: Supply all plots taken for this lab, make sure they are labelled properly.
  • Computation: Use MATLAB or MS Excel and calculate the LMTD and Heat Exchanger effectiveness values for the 3 tests.
  • Conclusion: Write a summary (max. 500 words, in a text box if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Applied Thermodynamics for Engineering Technologists by T. D. Eastop and A. McConkey: Heat Transfer.

6

Combustion Analysis

Learning Objectives

Operate the Plant at 80% capacity burning coal to

  • Perform combustion analyses for two types of coal,
  • Compare results.

Theory

In the Boiler Efficiency lab, we stated that Combustion Efficiency is defined as the ratio of the burner’s capability to burn fuel completely to the unburned fuel and excess air in the exhaust. In this lab, we will perform a combustion analysis.

Fossil fuels may be classified into solid, liquid and gaseous fuels. The vast majority of fuels are based on carbon (C), hydrogen (H2) or some combination of carbon and hydrogen called hydrocarbons.

During combustion, oxygen (O2) combines rapidly with C, H2, sulphur (S2) and their compounds in solid, liquid and gaseous fuels and results in the liberation of energy. Except for special applications such as oxyacetylene welding, in which a high-temperature flame is required, the O2 necessary for combustion is obtained from air. Air contains O2 and nitrogen (N2), plus negligible amounts of other gasses and for engineering purposes, may be considered to have the following percentage composition by mass:

O2: 23%
N2: 77%

The proportions in which the elements enter into the combustion reaction by mass are dependent upon the relative molecular weights as shown below:

Element Symbol Molecular Weight
Carbon C 12
Sulphur S2 32
Hydrogen H2 2
Oxygen O2 32
Nitrogen N2 28

Stoichiometric Combustion Theory

Complete combustion of simple hydrocarbon fuels forms carbon dioxide (C02) from the carbon and water (H20) from the hydrogen, so for a hydrocarbon fuel with the general composition CnHm, the combustion equation on a molar basis is as flows:

C_{n}H_{m}+yO_{2}\rightarrow aCO_{2}+bH_{2}O

Where the balance should be satisfied following the moles for any mathematcial  equation:

Carbon balance:

a=n kmol CO2/ kmol fuel

Hydrogen balance:

2b=m 

b = \frac{m}{2} kmol H2O/ kmol fuel

Oxygen balance:

2y = 2a+b

y = a+\frac{b}{2} kmol O2/ kmol fuel 

Considering that combustion occurs in air rather than in pure oxygen, the nitrogen in the air may react in the combustion process to produce nitrogen oxides. Beside, some fuels contain elements other than carbon, and these elements may react with oxygen during combustion. Also, combustion is not always complete, and the exhaust gases contain unburned and partially burned products in addition to C02 and H2O.

Air is composed of oxygen, nitrogen, and small amounts of carbon dioxide, argon, and other trace components. For the purposes of the further calculation it is perfectly reasonable to consider air as a mixture of 21% (mole basis) 02 and 79 % (mole basis) N2.  Nitrogen will be considered as an “inert” gas in the combustion calculations. The stoichiometric relation for complete combustion of a hydrocarbon fuel, CnHm, becomes

C_{n}H_{m}+y(O_{2}+\frac{79}{21}N_{2}\rightarrow aCO_{2}+bH_{2}O+cN_{2}

The balance equations are:

Carbon balance:

a=n kmol CO2/ kmol fuel

Hydrogen balance:

2b=m 

b = \frac{m}{2} kmol H2O/ kmol fuel

Oxygen balance:

2y = 2a+b

y = a+\frac{b}{2} kmol O2/ kmol fuel 

Nitrogen balance:

y\frac{79}{21}=c kmol N2/kmol fuel

Flue gas compositions are presented in terms of mole fractions as kmol of product per kmol of fuel.

Example 1.  Combustion of Octane in Air

Determine the stoichiometric air/ fuel mass ratio and product gas composition for combustion of octane ( C8H18)  in air.

C_{8}H_{18}+y(O_{2}+\frac{79}{21}N_{2}\rightarrow aCO_{2}+bH_{2}O+cN_{2}

Carbon balance:

a=n=8 kmol CO2/ kmol fuel

Hydrogen balance:

2b=18 

b = \frac{18}{2}=9 kmol H2O/ kmol fuel

Oxygen balance:

2y = 2(8)+b

y = 8+\frac{9}{2}=12.5 kmol O2/ kmol fuel 

Nitrogen balance:

y\frac{79}{21}=c=12.5\frac{79}{21}=47 kmol N2/kmol fuel

The combustion equation becomes:

C_{8}H_{18}+12.5(O_{2}+\frac{79}{21}N_{2}\rightarrow 8CO_{2}+9H_{2}O+47N_{2}

Air/ fuel mass based ratio considering that 1 kmol fuel is 114 kg of fuel ( 8*12 + 18 *(1) = 114):

12.5\frac{kmol\ O_{2}}{kmol\ fuel}\frac{1\ kmol\ fuel}{114\ kg\ fuel}32\frac{kg\ O_{2}}{kmol\ O_{2}}\frac{100\ kg\ air}{23\ kg\ O_{2}}=15.25\frac{kg\ air}{kg\ fuel}

Flue gas composition on molar basis is:

Total number of kmol of flue gasses = kmol CO2  + kmol H2O + kmol N2 = 8 +9+ 47 = 64 kmol flue gasses/ kmol fuel

CO2  = 8/64 = 12.5 %

H2O = 9/64 = 14 %

N2 = 47/64 = 73.5%

Other components and impurities in the fuel make the calculation process more complicated. For example if sulfur exists in fuel it is usually combust into sulfur dioxide (SO2). Ash, the noncombustible inorganic (mineral) impurities in the fuel, undergoes a number of transformations at combustion temperatures, will be neglected in the further calculation (ash will be assumed to be inert).

For most solid and liquid fuels, the chemical composition is on a mass basis, as determined in the ultimate analysis.

Mass Based Chemistry of Combustion

The combustion reactions are written following the stoichiometric rules as defined above. The quantity of matter entering into a reaction is equal to the quantity of matter in the products of the reaction.

The reaction for the complete combustion of C may be written as follows:

C+O2=CO2

or, if molecular weights are used,

12+32=44

1 kg C + 2⅔ kg O2 = 3⅔ kg CO2

The complete combustion of H2 occurs as follows:

2H2 +O2 =2H2O

4+32=36

1 kg H2 + 8 kg O2 = 9 kg H2O

Sulphur burns as follows:

S +O2 =SO2

32+32=64

1 kg S + 1 kg O2 = 2 kg SO2

In addition, 1 kg of O2 (stoichiometric mass of O2) is contained in 1/0.232=4.3 kg air which is the stoichiometric mass of air. This air will contain 4.3-1=3.3 kg N2. Therefore we can write:

1 kg S + 4.3 kg air = 2 kg SO2+3.3 kg N2

Procedure

If the analysis of fuel is given by mass, follow the steps below:

  1. Total O2 required: Determine the mass of O2 required for each constituent and find the total mass of O2 (Subtract any O2 which may be in the fuel)
  2. Stoichiometric air: Stoichiometric mass of air = O2 required/0.232
  3. Total mass of combustion products: Determine the % mass of each combustion product. For example, given C content of 84.9%, CO2=84.9/100*2⅔=3.11% and find the total mass of combustion products.
  4. Analysis of combustion products by mass: Suppose the total mass of combustion products in step 3 has been found as 12.09 kg/kg fuel, then CO2=3.11/12.09*100=25.74. This means that 25.74% of the flue gas is CO2. Repeat this calculation for each constituent.

Lab Instructions

Run the initial condition I14 80% Coal and setup trends for the following variables:

X00820 X00821 X00822 X00823 X00824 E23356

G02196 G02197 X32419 X02419 G00831 H00830

  1. Burning default coal: After 10 minutes of running the simulator, freeze simulator and print the two trends. This is the reference point for the next step.
  2. Burning poor quality coal: Switch to run mode and access Variable List page 0111 on MD180. Set the new values as shown below. After 10 minutes, freeze simulator and print the two trends.
    • X00820: 70.40
    • X00821: 5.10
    • X00822: 1.10
    • X00823: 12.50
    • X00824: 1.60
    • X00825: 9.30
  3. Computation: Compute the following values for both types of fuels:
    • Total O2 required
    • Stoichiometric air
    • Total mass of combustion products
    • Analysis of combustion products by mass [%]: CO2, H2O, SO2, N2
  4. Comparison: Compare your findings based on the following data:
    • Furnace outlet SOx flow
    • Furnace outlet NOx flow
    • CO content in flue gas
    • Oxygen content in flue gas
    • Theoretical combustion air
    • Coal Heat Value

Hints & Tips

In this lab, you are carrying out two combustion analyses. For data collection, set up your trends, a sample trend plot is shown below (Make sure your trend printouts are labeled properly otherwise, data analysis will be very confusing):

Sample data for coal

Sample data for coal

To change the fuel composition use the Variable List Page#: 0111 on MD180:

Chemical composition of coal

Chemical composition of coal

Deliverables

Your lab report is to include the following:

  • Trend plots: Supply all plots taken for each of the 2 fuels,
  • Computation: As per lab instructions above, perform combustion analyses using MATLAB or MS Excel.
  • Conclusion: Write a summary (max. 500 words, in a text box if using Excel) comparing your results and suggestions for further study.

Further Reading:

  • Basic Engineering Thermodynamics in SI Units by R. Joel: Combustion.
  • Thermal Engineering by H.L. Solberg, O.C. Cromer and A.R. Spalding: Fossil fuels and their combustion.
  • Stoichiometry.

1

Authors

Sanja Boskovic is the Chief Instructor at the Power Engineering Group of BCIT where she teaches Thermodynamics, Heat Transfer and Fluid Mechanics within the full time programs. Additionally, Sanja is the lead Thermodynamics and Applied Mechanics virtual classroom instructor for distance education and industry services courses. She was the 2011 recipient of BCIT Alumni Association Excellence in Teaching Award as well as the 2016 BCIT Teaching Excellence Award. Sanja is a member of the Institute of Power Engineers and American Society of Mechanical Engineers, and serves as a chair at Interprovincial Power Engineering Curriculum Committee and co-chair at Boiler Technology Committee at Technical Safety BC. Sanja holds a PhD in Mechanical Engineering and a PhD in Educational Technology and Learning Design.

 

Serhat Beyenir is an instructor at the Power Engineering Group of BCIT where he maintains web-based online courses and tutor distance education learners. He also teaches in the Power and Process Engineering program. Serhat is an open source software (OSS) advocate and open textbook proponent: He began to use OSS in 2001. Since then he has promoted Linux-variant operating systems, built content and learning management systems using Joomla, Drupal, WordPress and Moodle, and gave OSS related presentations. As a graduate student at Athabasca University, Serhat spent considerable time learning about open textbooks in his coursework and research assignments. Receiving a Master of Education degree in 2011, Serhat authored an open textbook utilizing the Connexions platform and he has been distributing his book free of charge to Power Engineering students since October 2011.

2

Revision History

Version 0.7 – December 6, 2017 (Additional copy editing and proof reading)
Version 0.6 – September 22, 2017 (Guest-lecturing at Kalmar Maritime Academy, Sweden)
Version 0.5 – March 28, 2017 (BCIT OpenEd Week presentation)
Version 0.4 – March 1, 2017 (Second Sprint event with our students)
Version 0.3 – February 1, 2017 (First Sprint event with our students)
Version 0.2 – December 19, 2016 (First draft)
Version 0.1 – November 3, 2016 (Initial development)