Valery Ochkov - Thermal Studies With Excel and Mathcad
Valery Ochkov - Thermal Studies With Excel and Mathcad
Volodymyr Voloshchuk
Thermal
Engineering
Studies with
Excel, Mathcad
and Internet
By general edition of Nikolay Rogalev
Thermal Engineering Studies with Excel,
Mathcad and Internet
Valery Ochkov Konstantin Orlov
•
Volodymyr Voloshchuk
Thermal Engineering
Studies with Excel,
Mathcad and Internet
By general edition of Nikolay Rogalev
123
Authors Editor
Valery Ochkov Nikolay Rogalev
Moscow Moscow
Russia Russia
Konstantin Orlov
Moscow
Russia
Volodymyr Voloshchuk
Rivne
Ukraine
1
It is often said: “reduced to the level of the Internet”. However, the Internet has also many decent
titles with proper content and design.
v
vi Preface
2
Therefore, by the way, many authors of study books do this.
Preface vii
ix
x Contents
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 303
Introduction
3
It is often also named engineering calculator (super calculator) or engineering office.
xi
xii Introduction
those calculations become confusing after a while for the author, not to mention
those, who want to understand these calculations and to supplement them. This also
applies to the above-mentioned “programs-monsters”, in which you can look at the
source code and decompile it. However, it is almost impossible for an external to
understand how the program works.
The language of MATLAB, which people often try to compare with the
Mathcad, is not a mathematical one. It is a programming language, which is too
difficult to master.
In addition, MATLAB is not good for itself, but its highly specialized appli-
cations, which are also some “monsters” with their inherent disadvantages noted
above, are good. We should also not forget that Excel and MATLAB unlike
Mathcad cannot work with physical quantities and units of measurement, which is
very uncomfortable and prone to error [8].
We can also say something about the math programs Maple and Mathematica,
which principally can also be used for engineering calculations. Primarily, Maple
and Mathematica are programs for symbolic math and computer analytic trans-
formations. In thermal engineering calculations, you can basically use numerical
mathematics, with some elements of symbolic transformations.
So Mathcad! Why is it so good for the thermal calculations? First, of course the
fact that one of the authors, who educates thermal power engineering and works in
the field of power engineering, is well aware of this software and even wrote a few
books about it [8–20].
But there are, of course, also objective conveniences for the work with Mathcad.
Here, they are as follows:
1. Good documentation of calculations
You can print out the calculation, which is made in the environment of Mathcad,
and give it to an examination or review to someone, who has never worked on a
computer. Mathcad calculations are similar to calculations made on paper by the
system WISIWIG (What You See Is What You Get). What you see (on screen), you
will get (on paper). Printouts of calculations made in Mathcad can be left in the
archive in order to read them in 50–100 years. The combination of documentation
and calculations performed in natural math notation will allows you to understand
what has been written in the past and reproduce it without much effort in a new
software environment, which will be available by then. This is very important.
Nowadays, we can see a crisis in information technology (IT), which people name
“nightmare of old software.” Imagine a company, corporation, or university that
developed and accumulated in 30–40 years of intensive use of computers a large
number of programs of varying complexity to compute and simulate a variety of
processes and also a lot of equipment and technology in various fields of science.
New computers, which replace old outdated computers, come with new hardware
and also new operating systems. Computers are combined in local area networks,
which do not remain aloof from the process of “globalization” and integrate into the
Internet—wired or wireless. The development of wireless computers to servers was
the cause of the appearance of “clouds” in information technology, which will be
Introduction xiii
discussed in the book. Modernization of the computer park often leads to the fact
that older applications will refuse to run on new or upgraded computers, work-
stations, and servers. Sometimes, you just cannot read the program from the media
(punched cards, punched tapes, floppy disks of different diameter, obsolete “flash
cards,” etc.), because new computers do not have the appropriate reading devices
(“UI”). Either you have just to abandon such programs and to create new ones from
“zero”, or you have to spend time and money to create or acquire a certain utility
(emulator) to run old programs on new or upgraded computers. But it is not so bad.
The real trouble begins when the professionals, who establish and maintain these
programs leave the firms, corporations or universities, and the younger people, who
replace them, cannot upgrade the programs to new requirements. This stems from
the lack of documentation of the code and the lack of appropriate software tools and
a simple inability to correctly read a program written in the old style—(“dead”)
programming languages. If you have managed to recreate or upgrade an old pro-
gram using the old or new programming languages, the “nightmare of old software”
will pop up again after a while.
2. Working with units
We have already mentioned that this is a very useful tool of Mathcad [8]. We
add only that Excel, MATLAB, and other programming languages weaned us to
work with measurement units or rather have taught us to work with dimensionless
quantities and their units (SI units—pascal, kelvin, joule, watt, etc.). Keep in mind
that, again, it is very uncomfortable and prone to error in the calculation. Even the
basic SI units are inconvenient: The basic unit of pressure (Pa) is very small and
always requires multipliers such as kilo or mega, and the temperature in kelvin is
not easy to understand for normal people and requires translation into Celsius or
Fahrenheit scales. The unit of measurement is completely dedicated to the second
study of this book.
3. A flexible system of variable names
The variables and functions in Mathcad with a few exceptions have the same
names that were defined for them in various scientific and technical disciplines long
before the appearance of computers. For example, the Greek letter η with different
indices represent the thermodynamics efficiency (thermal, internal, relative, etc.—
see the pictures of the book). This along with the use of traditional writing of
mathematical operators and functions makes the “language” of Mathcad to be
available to all people (see item 1 above) without any further comment. The
remaining four features are common to other so-called mathematical programs
(MATLAB, Maple, Mathematica, SMath, Derive, etc.), but, nevertheless, we tell
about these.
4. Numerical und symbolic mathematics
Mathcad allows us to use an extensive library of numerical methods in order to
solve mathematical problems. There we can anticipate or complement the analytical
solutions (successful or unsuccessful) of the problem. Mathcad originally designed
xiv Introduction
which the reference has been made. These and other Mathcad documents can be
downloaded to your computer or local network and used as templates.
Third way. Mathcad tools allow the user to quickly write and debug functions
that return, for example, the properties of working fluids, based on the formulas,
tables, or graphics, taken from an external source—from paper or electronic books,
as well as from the Internet.
These three features can and should be used in the calculation/document, for
example, for functions with parameters as arguments to a specific point of a ther-
modynamic cycle and return the desired thermal properties of the working fluid at
this point: specific enthalpy, specific entropy, density, specific heat (isobar or iso-
choric), thermal conductivity, viscosity, etc.
Without this, it is impossible to calculate thermal engineering processes. After
this description of Mathcad, we start our book with Chap. 1.
Chapter 1
Properties of Working Fluids, Coolants
and Structural Materials for Thermal
Engineering Calculations
Abstract You will learn how to create a function that returns thermophysical
properties of working fluids, heat transfer and energy materials using interpolation
and smoothing. You will also learn what a cloud function for the thermal calcu-
lations is.
A stranger from a far-away India, which, after our conversation on the forum
became the “a far-away Indian friend”, asked visitors of the forum to create a
function that would use as arguments the pressure ranging from 1 to 300 bar (see
Table 1.1, first column) and temperature ranging from –50 to +400 °C (see the first
row of the table) and then return the density of a substance in kg/m3 (content of the
table).
The author had a function for solving this problem at his fingertips, edited it a bit
and then posted the solution on the forum—see Fig. 1.1. In Mathcad, the first
row/column designated with the index 0. By using the built-in variable ORIGIN,
the user has the opportunity to change how to designate the first row/column.
Therefore, for example, if the variable ORIGIN is set to 1, you can call the first
row/column with the index 1.
In the program shown in Fig. 1.1, you can see how double-spline interpolation is
used. What happens in the program? First, the arguments of the function ρ lose their
units x ← p/bar and y ← t/K—273.15. Without this the built-in function cspline in
Mathcad 15 would not work. In the newest version, Mathcad Prime, you can also
use cspline with units—see Fig. 1.2. The last statement of the program (function
shown in Fig. 1.1) is, because of the desired unit density. This is a common
technique for empirical formulas—see Chap. 2.
The built in function submatrix creates a matrix out of the matrix M—the outset
(vector X) and “cap” (vector Y). The line with the return statement controls the
arguments of the function ρ, so they remain within the specified pressure- and
temperature range (see also Sect. 1.4 “extrapolation” in this chapter below).
Hereinafter the same function generates a submatrix Z—it filters Table 1.1, so there
will be a submatrix created without the “cap” and the outset. Next, out of the
columns of the density matrix Z an additional line (vector Zv) for a given pressure
will be generated by spline interpolation. Our “intermediate” in this example will
have for p = 250 bar, following elements 403.703, 306.196, 249.182, 211.977,
185.394, 165.393, 136.665 and 116.848 kg/m3 for the temperatures −50, 0, 50,
100, 150, 200, 300 and 400 °C, respectively (see “cap” of Table 1.1). Further, in the
program shown in Fig. 1.1 again spline interpolation is required for the density
dependence on temperature, absent of the “cap” of the table. Spline interpolation of
two vectors is shown in Fig. 3.3.
In Fig. 1.1 the last two operators show how to call the generated function ρ(p, t)
for different values of pressure (p) and temperature (t) and with various European and
American units. Function ρ(p, t) returns the default density (ρ) with the basic SI units
(kilograms and meters), but the user also may choose other units—pounds and feet, for
example. Working with units (not just with the values and the physical quantities)—
it is, we repeat, a unique feature of Mathcad, greatly facilitates calculations and pre-
vents many errors [8], which we describe in more detail in Chap. 2.
Figure 1.2 shows a program with double spline interpolation for Mathcad Prime,
where it becomes possible to use a matrix with different dimension. The function
shown in Fig. 1.2 does not depend on the value of the system variable ORIGIN.
Again the “cap” and outset of the table in Fig. 1.2 have the units pressure and
temperature. In order to save space the units of the table content (density) not
shown.
For a function with two arguments, it is easy to construct the surface in Mathcad.
Figure 1.3 shows that the function ρ(p, t) is smooth, increases with increasing
pressure and decreases with increasing temperature.
It Fig. 1.3 shows typical “Thermophysical behavior” of many substances (gases)
in the single-phase region. The graph of the function is not only useful for visual
analysis of the behavior of matter, but also to fix possible errors and typos. Thus, if
in the matrix shown in Fig. 1.1, intentionally or accidentally the density
107.1 kg/m3 (Table 1.1 allocated in the frame) is overwritten by 170.1 kg/m3
(a very common error when publishing numbers in magazines and books, and/or
when manually entering numbers on a computer), the surface shown in Fig. 1.3,
change its form—see Fig. 1.4. However, this is not so simple. A splash on the chart
may be not just a mistake, but also maybe some anomalous behavior of substances.
Therefore, when seeing such a “mistake”, you need first to check before correcting
it—it is either a typo, or a scientific discovery. The graph is not created by inter-
polation (conducting the surface through the points)—it is just created by
approximation (smoothing—conducting surface near points) of the table.
Of course, you would have to immediately ask the “distant Indian friend” either
what this substance (gas) is or if the formula is known. This rather ordinary episode
of Mathcad problems raises a very important issue, which will be outlined in this
chapter.
and certified server with parameters of the needed substance and the server should
return to the user’s computer the needed properties.
How is it realized? A specific example: the Mathcad document shown in Fig. 1.1
placed on the server http://twt.mpei.ru/TTHB/Ro-pt.xmcdz. If the property you
have to calculate is the density of the substance is a function of pressure and
temperature, you only have to insert a reference (Fig. 1.5) to the desired
Mathcad-document (shown in Fig. 1.1) and the function ρ(p, t) is visible in your
calculation.
After the data entry of t and p and entering a reference to the file named Ro-pt,
the “cloud” function ρ(p, t) becomes available. This function is here called with the
parameters p = 20 MPa and t = 350 °C. This function further visualized by some
isobars. When entering t = 500 °C (this temperature is outside of the range covered
by Table 1.1), the function ρ(r, t) returns an error message.
On the current server, which is a collaboration of specialists from the National
Research University “Moscow Power Engineering Institute” (www.mpei.ru), Joint
Institute for High Temperatures (www.jiht.ru) RAS and Ltd “Trieru” (see www.trie.ru),
collected a large number of such “cloud” Mathcad-functions. In addition, their number
is growing.
These functions created for different reasons and in different ways.
1 Properties of Working Fluids, Coolants and Structural Materials … 7
First, if the author sees in a book or online a table similar to Table 1.1, he
instructs his students to “revive” it by using, for example, by the program shown in
Fig. 1.1. Sometimes it is enough to insert it into a new matrix and make some other
minor changes. If in the book or on the Internet the formulas showed, which used to
calculate the properties of substances, it is significantly easier. We have also
developed a technology to transfer schedules to the “cloud” function. All of these
techniques described in Ref. [21], which is released by Russian Foundation of Basic
Research (www.rfbr.ru) grant fifth additional volume reference series [22].
Second, these “cloud” functions customized for Mathcad for working in the field
of thermal and nuclear power industry [23]. Thus, for example, have been created
and placed on the server functions on thermophysical properties of ethanol for
calculation of steam turbine cycles with organic working fluids (ORC—organic
Rankine cycle), for the calculation of refrigeration units [24] and other cycles. ORC
posted at http://twt.mpei.ac.ru/MCS/Worksheets/PTU/Rankin-Ethanole.xmcd.
The main part of the “cloud” functions hosted on MPEI-JIHT-Trieru, linked, of
course, with properties of working fluids and coolants. These functions described in
Ref. [23], and the technology they use in the calculations in the articles [24, 27].
Functions related to water and steam (which is the main working body of domestic
and foreign thermal and nuclear power), based on the formulations approved by
the International Association for Properties of Water and Steam (IAPWS—
www.iapws.org; authors of this book are members of this Association). These
formulations not based on table data, but formulas. For example, the pressure and
the saturation temperature of water in the range from the triple point to a critical
point are associated with quadratic equations. However, in practice very often
calculation formulas replaced by interpolation from individual points. That is why
we pay special attention to this issue. This does not only apply for calculating the
properties of saturated (a function of one argument), but also for individual
single-phase regions and areas of subcooled liquid or superheated steam (a function
of two arguments). Waiver of calculation formulas and the transition to working
with tables can be considered a step backwards in the historical process of creating
databases on the properties of substances. At the dawn of this process, such data-
bases published in the form of tables and/or graphs (visualized tables). Reading a
table or graph without a computer or calculator is much faster than using the
formulas. Loss of accuracy here was not so important. With the advent of electronic
computing people slowly began to move away from tables, graphs and formulas
and moved to computer programs, where these formulas programmed. Formulation
process of not yet fully investigated or newly synthesized compounds carried out as
follows: create a table with the so-called skeletal material properties through
experimental measurements with various physical methods. Then out of these
discrete tabular data one common function or a set of functions is generated by
various mathematical methods for different areas or applications, and predicting
possible errors for the primarily use on computer programs—this is called the
formulation. However, often there is only a table published without the mathe-
matical processing. Having at hand powerful and convenient means of interpolation
—such, for example, which built in Mathcad (e.g., Fig. 1.1), you can refuse to
8 V. Ochkov et al.
The program shown in Fig. 1.1 is easy because the original rectangular matrix is full.
If the matrix is square shaped, it will simplify the work. Spline interpolation on a
square matrix in Mathcad shown in Fig. 1.1. If you want to cover a wide range of
parameters of a substance by interpolation, including solid, liquid and gaseous
phases, the corresponding table of material properties (density, for example) will be
1 Properties of Working Fluids, Coolants and Structural Materials … 9
Fig. 1.6 Not completely filled (sparse) matrix with fictitious numbers
Fig. 1.7 Not completely filled matrix with NaN values and linear interpolation
Source data for tables cannot only taken from directories and the Internet, but
also from other programs. Data on the properties of ethanol to calculate the
above-mentioned steam turbine cycle with organic working fluid (see separate
study) and some refrigerants (see separate study), were taken from the program
RefProp [29]. This program can connected directly to Mathcad through the
mechanism of DLL (Dynamic Link Library). However, this technology does not
work with all versions of Mathcad, and only few can do it. Therefore, (alternatively)
the technology described in Figs. 1.1 and 1.6 is used. RefProp generates the desired
vector or matrix, which will be inserted into a Mathcad function, which is similar to
that which is shown in Fig. 1.1. You can obtain a function, which returns the
desired property of a substance with the right dimensions and has dimensional
arguments, but does not work with equations of state, which hidden in the RefProp,
and can be viewed and edited if necessary. Figure 1.9 displays the process of
generating a matrix in the RefProp software with two columns that hold the tem-
perature and pressure of ethanol on the saturation line. In the table (she left in
Fig. 1.9) it is sufficient to select the desired track and copy and paste it in
Mathcad-document (right part of Fig. 1.9), in which the number will be converted
to a string (character string enclosed in quotation marks, where we need to remove
the commas and replace them by points, and then Mathcad-function str2num to
translate it into numbers). Another way to obtain direct numbers in
Mathcad-document from a program like RefProp—is to transfer data to a spread-
sheet in Excel, which will be transferred to a Mathcad-document. Data source tables
similar to that, which shown in the left half of Fig. 1.9, can often be found on the
1 Properties of Working Fluids, Coolants and Structural Materials … 11
Internet. They are easy to move in Mathcad-document with the methods described
above. Here the main “fuss” is about the problem of “point—comma as decimal
separator”. Create a function on the table with two columns is not difficult, as
shown in Fig. 3.3.
The interval between two temperatures can be increased in order to save com-
puter memory or decreased in order to increase the accuracy. Speaking of computer
memory—previously, it was one of the main limiting factors of calculations.
Therefore, interpolation was often refused, because of the required storage of large
data sets on the computer memory, and the work on equations of state proceeded.
Now this restriction eliminated.
12 V. Ochkov et al.
mind if you claim 1 in Fig. 1.9, while the second—the guess value to the desired
parameter—the pressure (see Section 2 in Fig. 1.10).
Our function, the creation of this is shown in Fig. 1.1, is continuous and
monotone dependent on each of the parameters in the selected range of pressure and
temperature. When you work with more complex regions of state you should be
noted that for such functions or there is no backward functions somewhat (which all
have the same domain, but different value). In this case, the given parameters need
to know exactly what the backward function we will seek and will need to be
created in the backward function to provide an additional argument—the initial
approximation to the desired solution. If it theoretically proved that in the defined
area the backward function is uniquely defined, we will get the same result for any
approximation. For the existence and uniqueness of constructing the backward
function of one variable its strict monotony is sufficient. In general, the search for
the backward function is a complex mathematical problem.
A specific example: the water pressure and/or steam temperature, and the
specific enthalpy can then be either determined in the 1-phase region (pressurized
water) or in the two-phase region (wet steam). This situation described in Chap. 5,
wherein the compression curve is constructed in the T, h-diagram.
Sometimes it is useful as a first approximation to use a random numbers in a
specified range as values of the backward function. Figure 1.11 illustrates the
change of the pressure dependent on the temperature for different χ (χ—the heat
capacity ratio), for χ 1.26, 1.27… 1.4. The wspPTK function based on function
wspKPT, which is included in the WaterSteamPro software (see www.wsp.ru). The
wspPTK function returns different values depending on the first approximation.
Random number generator helped us to construct the required family of curves,
although in fact they are not the curves, and the new family of points with
14 V. Ochkov et al.
discontinuities in those places where the Mathcad could not find solutions function
root.
In Fig. 1.11, as well as in Fig. 1.10, we use to create the backward function the
built-in Mathcad function root. However, as a first approximation of the pressure
we use a random value in the range of 0–30 MPa. This done with the built-in
Mathcad function rnd, which returns a random number in the range from zero to
the value of its argument. In Fig. 1.11 you can also see a Give-Find block. There
pressure and temperature is calculated at which the partial derivatives of wspKPT
are equal to zero. In mathematics, this is called a saddle point.
About Fig. 1.11 an interesting story can be told. One of MPEI’s teachers pro-
mised his students, that he will give to the one, who can construct this family of
curves using Mathcad, a six in the exam (not five—actually the best mark in
Russia). The instructor was hundred percent sure, that this can be done only with
1 Properties of Working Fluids, Coolants and Structural Materials … 15
the help of more powerful programs. Students were not able to do it, and asked the
author of this book for some help. The problem solved, but none of students
received a six in the exam, because they asked for help from outside.
In Mathcad, with “direct” functions on the properties of materials and tools for
solving equations and systems, you can understand the essence of this problem; it is
easy to create backward functions. We will show this in some studies of this book.
The issue of creating backward functions will continue in Chap. 11.
1.4 Extrapolation
the matrix M are “fictitious” number prescribed white on white. This technique
works with non-standard data and described in the following [21] (in Fig. 1.6 such
data registered “black and white”, but in brackets.)
The function ρ(NaCl) checks for wrong input data and returns an error message,
if the third argument of the function (Check) has a value of 1 (not zero). This shown
in Fig. 1.14.
Figure 1.14 shows three calls of the function ρ(NaCl): “normal” call, when the
arguments of the function are within the acceptable ranges of temperature and
concentration of the solution NaCl, “false”, when the third argument (Check) is 0,
and the function output of ρ(NaCl) is wrong, extrapolation value and clipping the
wrong case, false answer.
Figure 1.15 also shows how to create an backward function on the basis of the
function ρ(NaCl).
The new function ω(NaCl) returns the NaCl solution concentration depending
on the temperature and the density (for example: in a NaCl solution a thermometer
and hydrometer was inserted and now we want to know its concentration). If the
function is enabled ρ(NaCl) check the source data, but the backward function does
not work.
The created backward function can be supplemented by operators in order to
avoid wrong answers—see Fig. 1.16.
The function shown in Fig. 1.16, determines the concentration of NaCl in the
solution depending on its temperature and density, without any limitation, since
1 Properties of Working Fluids, Coolants and Structural Materials … 17
third argument of ρ(NaCl) is zero. But before the output response it is monitored, if
the values of the function arguments of ω(NaCl) and the values of the function are
in a reasonable range. The density value controlled by checking whether the value
does fall within the range of the minimum density (1075 kg/m3) and maximum
(1204 kg/m3) density—see the table in Fig. 1.12 and the matrix in Fig. 1.13.
Remark. Two pitfalls can be noted in the function shown in Fig. 1.16. First, the
root feature seeking zero value of functions not near the point 20 % but within the
range 10–25 %—see Fig. 1.16. Secondly, the function, the creation of which shown
in Fig. 1.16, when erroneous data returns, not text message detailing the user
response.
Speaking of extrapolation among Mathcad—functions with the root name spline
(see Fig. 1.1) can have three preposition (prefix)—l (el), p and c: lspline, psline and
cspline. These prefixes represent linear (l), square (parabolic—p) and cubic (c) extrap-
olation. Inside the area of tabular data interpolation is conducted by cubic extrapolation.
One of the authors of this book is often criticized that he “in his life has never spent
a single measurement of properties of substances, and is only concerned with the
fact that someone else takes the published experimental data, processes them, and
1 Properties of Working Fluids, Coolants and Structural Materials … 19
or the authors themselves. However, If a user, for example, need to clarify, either
the argument of sine is measured in degrees or radians, he will not check the
documentation, and just simply write x := sin(30) and see what result he gets for the
variable x. Such experiments give users an option, i.e., checking the documentation
only in severe cases (or not) and usually not finding the answer there. This appeals
to the creator as problematic. On the hotlines do not sit the creators of the program
—usually just some mildly skilled employees, and as the reader knows, it is not the
same thing. Consultant firms are likely to ask you to call back in a couple of days
during which they will perform some experiments with the program and try to find
the answer. Yes, and according to the author, you often do nothing in the meantime,
forget about the problem and will be fully absorbed to a new project. Even if it is
not this case, then the author cannot remember all the details of his program after he
calls back.
Because of this when working with programs the user often forgets that it is a
product of the mind and human hands, i.e., which is the fruit of the work of
anonymous persons, that cannot be reached via a hotline. Here, apparently, lies the
philosophical explanation (but not in any way an excuse) for the widespread use of
illegal copies of programs. Here we are talking about relatively honest people that
put the program on your computer in order to get to know more about nature and
pass on their knowledge, such as students.
We can assume that the Creator, throwing an apple on top and please them on
the head of Newton, we opened one of the secrets of his Divine Plan. The merit of
genius (Newton)—the ability to be in the right place at the right time. Program (and
not only brilliant) would also like to fall on top of us, and there is no obligation for
them to pay crazy money by Russian standards. Hence the idea of vitality and
shareware, which oppose primarily sellers and not the authors of programs.
We can assume that God created from his own image and likeness, not only
humans, but also the computer. Writing for a program—the inhalation of soul into a
lifeless pile of iron, plastic and sand. Trading bodies (organ transplants, blood
transfusions or, finally, prostitution)—is a reality of our days. Buying and selling of
souls is found only in fairy tales and fiction (the story of Faust, for example). By
purchasing the software, we buy only the “body”—disks, documentations, infor-
mation and discounts on new versions, and most importantly—after-sales service.
Hence our conclusion: experiments with programs like RefProp [29],
ThermoData Engine [30], WaterSteamPro [31] or Thermoflow (see Introduction)
have the same right to life as classic experiments.
Often you can hear the statement (another rebuke to the author) that, Mathcad is a
“frivolous” program, designed only for pupils or students. Engineers and scientists
need to work with more complex programs, like MATLAB for example.
What can I say?
1 Properties of Working Fluids, Coolants and Structural Materials … 21
In the table, placed at the beginning of the Chapter (Table 1.1 and Fig. 1.3) the
Celsius temperature scale used. This is not a typical case. Usually these tables based
on the absolute temperature—on the Kelvin scale. Working with different temper-
ature scales (Kelvin, Rankine, Celsius, Fahrenheit, and sometimes Reaumur) often
leads to confusion and errors, when we work with tables. Mathcad greatly facilitates
the work with its ability to work with units. In addition, we must remember that old
tables, graphs and formulas, which display the properties of substances, can rely on
22 V. Ochkov et al.
the temperature scale of 1968 (ITS-68). Currently a new scale since 1990—ITS-90,
is used. To set the “cloud” features on the properties of substances applied a pair of
functions, providing for a new recount temperature (ITS-90) and old (ITS-68) scales,
which eliminates possible errors in the use of old, but, nonetheless, demand data.
The site http://twt.mpei.ac.ru/PVHB/pvhbeng.html—a web version of the book
“Physical quantity” (http://twt.mpei.ac.ru/PVHB/pvhb.html) in the “thermometer”
the reader will find on-line recalculations on the temperature scale and corre-
sponding functions for embedding them in Mathcad.
Figure 1.17 shows the difference of the both temperature scales (ITS-90 and
ITS-68). On the site you can also find links to functions for recalculations of
metrological standards and description of these.
In Old Russian (Soviet) references, along with the old temperature scale can be
(see above) to meet and outdated symbol of the unit of temperature, the temperature
Fig. 1.17 Difference between old (ITS-68) and new (ITS-90) temperature scales
1 Properties of Working Fluids, Coolants and Structural Materials … 23
1.8 Uncertainly
Fig. 1.18 Calculation of the specific isobaric heat capacity of water, indicating the error
24 V. Ochkov et al.
1.9 IT-Security
1.10 Templates
Figure 1.19 shows how Word offers its users not to start writing on a blank sheet
of paper (which, incidentally, is also a template, however, a primitive one), and to
select the appropriate template.
The ready-made Word document (official letter, for example), we repeat, consist
of two parts—a constant (name of organization, its details and so on.), which was in
the template and a variable part—what we enter in the template (main part of the
letter). All this significantly accelerates and simplifies the work, and reduces the risk
of errors therein.
Mathcad also comes with a few templates. Usually they are rarely used when
working with Mathcad due to the fact that after the launch of this program not a list
of patterns (as in Word—see Fig. 1.19), and a blank document (the most primitive
pattern) appear. If we press the small button next to the icon with the image of a
new blank document in Mathcad and give the command “New” from the File menu
(Mathcad 15), then the screen will display a list of templates, which are built into
Mathcad—see Fig. 1.18.
Not only built in templates can be used (see Figs. 1.20 and 1.21), but also user
created ones. Creating a template is easy—you need to open a built-in template,
change something in it and save it in the folder “template” with a new name and the
ending *.hmct. The same operation can be also done with a finished document.
26 V. Ochkov et al.
Figure 1.19 after adding a new template the user can use his own created template
just like the built-in templates.
If you open a custom template, as shown in Fig. 1.21, Russian units can be used
in the Mathcad document (if Russian units have been introduced in this template),
which are not built into Mathcad—see Fig. 1.22. (You can also refer to the “cloud”
Mathcad document with full name http://twt.mpei.ac.ru/ochkov/units.xmcdz, and
some additional Russian and English units of measure become available. They also
used in the “cloud” document named http://twt.mpei.ru/tthb/H2O.xmcdz, to which
we often refer in this book.)
Custom pattern shown in Fig. 1.22, a built-in Mathcad standard template
(Normal), which inserted in an area named Additional units. In this area, the
operators stored kJ: = 1000 * J, g: = gm, etc., i.e.
Among Mathcad Prime 3.0 the opportunities to work with templates are
expanded. In particular, you cannot only use your own templates (My Templates—
see Fig. 1.23), but also templates from all users of the local computer network
(Shared Templates) and standard templates from the firm PTC—Developer
Mathcad.
With the ability to share templates among all users in the network, it is possible
to replace Mathcad servers in those organizations, which for various technical
reasons and/or for security reasons do not have access to the Internet.
Mathcad templates primarily contain no text (except heading, etc.) and calcu-
lation elements. Thus, Fig. 1.24 shows a template for the use with aqueous NaCl
(this will be described in detail in Chap. 3). In this pattern are following functions
integrated: a function that returns the density of NaCl in aqueous solution
depending on the concentration, expressed in different ways—the weight percent
(%), molality (mole/kg), molarity (mole/L) and the titre (g/mL), the backward
function (concentration of density), a scaling function of concentration depending
on the function of the freezing point of the solution, and other useful functions. All
this considerably simplifies and accelerates the calculation where this solution
appears—for example, the calculation of air conditioning systems with this solution
as a coolant or refrigerant.
Looking ahead, we can, say that a template for Mathcad 15 may contain links to
cloud features on the properties of water and steam—the main working body in
engineering sciences.
Figure 1.25 shows the calculation of the density of water for a given temperature
and pressure with a function named wspDPT, which is stored in the “cloud” file
named H2O.xmcdz. More details about this calculation described in Chap. 6
“Maximum density of water”.
Working with templates can greatly simplify the calculations. Starting a new
project does not require opening a new blank document and not be lazy and search
the World Wide Web right template with the right functions. The used template can
be edited and expanded—add to it, for example, new functions and return it back to
the Internet, providing new functions for other users. It is one possible way of
bulk solution, to the other—the ratio of the mass of solute to the mass of solvent,
and the third—the ratio of the amount of solute to solvent, etc. So use this data with
great care.
Another modern way to publish data about properties of substances is to create a
website, which contains a code for specific programs. Thus in Fig. 1.28 you can see
a website, which generates a code for engineering calculations with MATLAB.
In conclusion we can say that the International Association for the Properties of
Water and Steam (www.iapws.org) uses for interactive website design calculations
basic properties of this working fluid and material engineering (e.g., Fig. 1.28).
There’s a link to a white paper with instructions how to calculation the thermal
conductivity of water and steam (PDF of document), then a short description of this
1 Properties of Working Fluids, Coolants and Structural Materials … 31
Fig. 1.28 Internet generation user functions on the properties of substances for MATLAB
Fig. 1.29 Link to the on line calculation website from the site IAPWS
Fig. 1.30 Page of the site with an interactive calculation the thermal conductivity of water and
steam
but also all intermediate data. When debugging a program comparing intermediate
calculations will quickly find and resolve the error, if it took place. In addition, the
“live” calculations, links to which are shown in Fig. 1.29, contain diagrams with
area where operating point is marked, and the value of the calculation error is
shown (Fig. 1.30).
Chapter 2
Working with Physical Quantities:
Problems and Solutions
Abstract This chapter deals with the use of dimensional quantities in thermal
calculations. The main aspects of work with dimensional variables in thermal
calculations are given. The concepts of physical and empirical formulas are
extended by pseudo-empirical formula concept. Few samples are provided.
the density. However, the density was used with wrong units and therefore the tanks
have been just half-full.
A similar mistake has led to the loss of the American satellite Mars Climate
Orbiter (MCO). MCO launched on 11 December 1998 with a PH-2 Delta rocket.
The unit arrived at Mars after 9 months, at the 23 September 1999. The engine
should break the unit and go into the orbit of Mars. After 5 min, when the satellite
went behind Mars, no more signals have been detected. From the analysis of data
the scientists found out, that the device has passed the surface at a height of 57 km
instead of 140 km and simply burned in the atmosphere (unit cost $125 million).
This large deviation was caused by the fault of experts, who prepared the mission:
when calculating the braking pulse one of the groups used Imperial-units
(pound-force) and the other SI-units (newton).
You can mention another engineering error associated not with units themselves,
but with the scale of measurement. In the past was a bridge across the Rhine built.
On one side the Germans started and on the other side the Swiss started to build it.
When both side met, they recognized that the height difference is about a half meter.
The reason for the error was that the German standard of the construction zero
height is the average level of the North Sea, and in Switzerland the average level of
the Mediterranean Sea.
Any expert, who works in the field of science and technology, can tell numerous
examples of such errors and incidents associated with the incorrect handling of
units. The transition from calculations by hand to “dimensionless” programming
languages, could not either solve this problem.
Mathcad—it’s not just physical/mathematical program. Often people do not use
units and just mention the unit of the used constant/variable just in the comments.
So you can, for example, often see that the pressure is just written down in
Mathcad as: p := 120, instead of using the more clear and correct way: p := 120 atm.
What are the reasons for the under-utilization of Mathcad? Firstly, some users
are just unaware of this useful tool. It is possible to use units in Mathcad and then
nevertheless transfer the calculation to another programming language or spread-
sheet program, in which, I repeat, the variables are stored as numerical values and
their units noted in the comments.
There is a second group of Mathcad-users, which do not use physical quantities
in their calculations, explaining or justifying it by the fact that all quantities pre-
scribed in only one fundamental measurement system (e.g., international SI-Units)
and they do not have any problems with the translation of the units. This motivation
often buttressed by the fact that Mathcad-document without physical quantities is
much easier to prepare for transferring to other programming languages. In addi-
tion, a Mathcad-document with disabled mechanism to work with units of mea-
surements works faster.
The third and the main cause of the failure of many users to work physical
quantities in the calculations lies deeper. It is associated with some features and
shortcomings associated with the work with physical quantities, which even the
2 Working with Physical Quantities: Problems and Solutions 35
most experienced users sometimes force to translate them into the category of
comments.
However, using comments represents extra work. Often the user says to himself
that he will add the necessary comments later. Often this “later” doesn’t come and
the comments are never added. Nevertheless, when someone wants to use this in the
past-created document he cannot understand which values he has to use for the
variables, because he does not know which units are used there. Even for the creator
of the document it can be difficult to remember which units are need to be used,
because the document was created too long ago in the past.
Using units in Mathcad is quite simple. When you enter a numeric value you can
connect it with units through typing in a built-in/user-unit or also a functional
dependence. Additionally you can also choose a unit through the dialog box
(Fig. 2.1). Usually the multiplication sign between the numerical value and unit is
invisible. When you enter a dimensional value the option can be given to the user to
select the unit of measurement, for example, through the use of switches: the
temperature/pressure can be shown at different temperature or pressure scales with
different units of temperature/pressure.
Among Mathcad (Ver. 13, 14, 15) it became possible not only to work with the
built-in units (SI, MKS (meter, kilogram, second), CGS (centimeter, gram, second)
and the U.S. Units), but also with custom units. This in particular means that the
user can modify one of the predefined units or create even new units. Further
calculations with custom units are in the same way monitored like calculations with
old units: meters, for example, cannot be added with kilograms; which helps to
avoid errors and typos.
When you use the “ = ” operator Mathcad will display the following: ∎ = ∎ ∎.
The first place holder contains variable name, the second one the numerical value,
and the third one the unit. By clicking on the third placeholder you can change the
used unit and this convert to another system of measurement (e.g., bar → Psi). In
some cases it is useful to duplicate the result value and show this way the result
with different units; so the reader can choose, which one is the most convenient one
for him.
What forces even experienced users to avoid units in purely physical calculation
or to withdraw them from an almost finished document?
These “pitfalls”:
Some Mathcad tools are not designed to work with dimensional values. They
interrupt the work or return an error message such as “There should be no
dimension value” or, much even worse, give the wrong answer. In such cases it is
necessary to temporarily withdraw the units by dividing the variables by the unit,
and then after the calculation add the unit by multiplying the variable with the base
36 V. Ochkov and K. Orlov
unit. This method will be described in more detail when we talk in more detail
about working with empirical formulas.
If the user does not plan to use physical quantities, then the use of units should
be disabled. Therefor you have to make the following commands: press Tools,
Options, and then a new window will pop up. There you press on the button Unit
System and set the switch to No (see Fig. 2.1). This way it is possible to avoid
errors. For example, when the user wants to create a function like f(x) = x2, and
accidently uses instead of x an unit like: A, S, V, etc. Mathcad will not return an
2 Working with Physical Quantities: Problems and Solutions 37
error as expected. When using Quickplot we will see just one point instead of the
functional correlation. When turning of units we can use these former unit-names
also as variables (like x). This way mistakes can be avoided.
Symbolic math units operate in Mathcad as simple variables, they do not take
into account that 1 m is equal to 100 cm or an hour has 60 min, etc. The symbolic
math tool is a foreign element in Mathcad. It is taken out of the program Maple, in
which the units have appeared only in the eighth version (Mathcad 5–13), or out of
the program MuPAD, in which are either no units of measurement used (Mathcad
14/15 and Mathcad Prime). Symbolic Mathematics is an auxiliary tool in Mathcad,
which is rarely used in calculations, and only plays a secondary role. If you want to
use units in symbolic mathematics, than more should be done: it is necessary to use
the substitute command in Mathcad, and tell Mathcad that 1 m = 100 cm,
1 h = 60 min, etc.
In Mathcad 15 you can store only dimensionless quantities or quantities of the
same dimension in vectors/matrices, i.e., the same physical quantity like time,
strength, weight, etc. There is known only one exception to this rule. The function
Find, for example, can return a vector with mixed elements. If you need the results
of several different function united in a vector or matrix, it is possible to deprive the
units first and then bring them back. This restriction has been removed in Mathcad
Prime.
Often so-called empirical formulas are used—formulas that relate not just on
physical interrelations. Example of an empirical formula showed on Fig. 2.2 pro-
vides the information how to calculate with Mathcad the (estimated) used car value
depended on age and mileage.
The empirical formula, with which the price of an old car can be calculated, was
created with the help of statistical processing. This example shows how to sup-
plement these formulas without interference with physical quantities. When
empirical formulas are published it is always clearly marked, what units should be
used as the initial values and in what units the answer will be. In this case, the
vehicle age must be entered in years (yr) and mileage in Miles (mi). The formula
returns the car price in U.S. dollars ($). These variables must lose their units in the
formula. Figure 2.2 shows the formula of the CarPrice where variables are divided
by their assigned units, the formula itself is multiplied by the currency. (Figure 2.2
was created at a time when the U.S. dollar was worth 50 rubles.)
The list of usable units in Mathcad does not contain all units. There are missing,
for example, often used units in information technology, like bits, bytes, etc. In
Mathcad 14 are also missing currency-units (dollar, ruble, euro, etc.). Mathcad 14
deals only with “seven” SI-Units (time, length, mass, current, temperature, light
intensity and amount of substance) and their derivatives—units, which composed of
these basic units (force, energy, power, etc.). And how to deal with older Mathcad
versions, if you want to perform the so-called economical calculations, in which are
dollars, rubles, euros, etc. are present? One solution is shown in the upper part of
Fig. 2.2. The Dollar ($US) is assigned with a random unit of measurement, a
physical quantity which does not appear in this calculation—candela (unit of
luminous intensivity): $ := cd (see the first operator in Fig. 2.2). Next, you can
38 V. Ochkov and K. Orlov
connect the dollar with other currencies used in the calculation, for example,
Rub := $/50. After that the price will be displayed by default, candelas (Mathcad 14
and 15), or in base currency units, which has to be replaced by the appropriate unit—
dollars or rubles (see the last row in Fig. 2.2).
Another example shows Fig. 2.3, where we have to solve the following problem
of thermal engineering: there is a power plant efficiency η given and you need to
calculate the specific fuel consumption for electricity generation. All information
for this calculation is given. The formula for the needed fuel is: b := 12300/η. It is
said that the efficiency (η) must be expressed as a percentage (or in relative units)
and the result of the needed fuel will tell us how many gram (gm) of fuel we need
for the production of 1 kWh of electricity. For example: 12300/32 = 384.4 or
2 Working with Physical Quantities: Problems and Solutions 39
Fig. 2.4 Formula for the conversion of concentrations from the Internet
is the mistake. Although we use in this case even the right units we still get the
wrong result. The Multiplier 1000 is integrated in the formula because of the
relation between the units g/kg and cm3/L. If we use the pseudo-empirical formula
as a real empirical formula in Mathcad, we will get the correct result (see Fig. 2.6).
If we want to use the formula like a real physical formula we need to think about
the multiplier 1000. In this case we need to delete it. After correcting the formula
we will get the correct result (see Fig. 2.7).
Units support for physical quantities was till the 13th Mathcad version only
possible on an absolute scale. This means in particular that if the value is zero itself,
then it is not necessary to attribute any unit, although nevertheless it should be done:
for example, l := 0 m, in order not to disrupt the possibilities to control the units. But
often we also need a relative measurement scale for some tasks. For example, we
usually measure the temperature in degrees Celsius (relative scale), and not in Kelvin
(thermodynamic temperature—absolute scale). In Mathcad (also in new versions) we
do not have a multiplication sign between the numeric value (25) and unit (°C), when
we use a relative scale (t := 25 °C). We have more functional correlation between the
unit and the numeric value, connected by a postfix-operator.
Figure 2.8 shows one possible solution of the problem of using relative mea-
surement scales: given is the temperature t1 (heater entrance) and the temperature
difference Δt (difference between the entrance and exit). It is necessary to determine
the temperature t2 at the outlet of the heater. There you can see three objects called
°C (each time another User has to be selected in Mathcad): two constants named °C
(the first is called °C := K and the second °C := 1) and a function called °C.
Additionally you can see another function with an invisible name—(T) :=
T/K-273.15. Invisible function name can be entered by space sign (using
Ctrl-Shift-K) or using text style with white color of characters.
The names of all three objects coincide (°C). However they are different objects,
because they are used by different users in Mathcad (selectable via the button left
next to the font-button). When working with relative temperature relative scales are
three situations, in which the above introduced objects called °C are useful.
Situation 1
In the calculation it is necessary to enter the temperature in Celsius. To do this, we
use the function called °C with the postfix operator. The input of the function will
be temperature in °C. Internally Mathcad will convert the Celsius temperature into
kelvin. Therefore all calculations can be carried in Mathcad, since Mathcad works
with the absolute temperature scale.
Situation 2
You must enter the temperature difference Δt. It is very important that differences
entered with another functions. Usuall (and very bad) mistake is to use °C. By using
°C you will increase temperature difference value for 273.15 degrees. Since Δ°C is
equal to K you can use the object °C = K.
Situation 3
You want to convert the temperature from Kelvin to Celsius scale. To do this you
use the prefix operator and the invisible operator. The result will be first displayed
without any units. So u have to use additionally the °C = 1 in order to add some
“virtual” units. Otherwise you can also just write the units in a comment.
With the described situations and the three similar objects named °C, it is fully
possible to realize the work with each temperature scale: entering the temperature in
one scale and converting it, or calculating the input and output values of a
and on the other hand, they opened the possibility of making mistakes such as
addition of radian. This problem of lack of control of the dimensions you may
encounter when working with other “dimensionless” quantities.
Units referred to the number of pieces (e.g., price per piece) can cause some
additional problems. Thus, in the calculation, which is shown in Fig. 2.12, you have
to calculate the total surface of 2200 tubes with a diameter of 15 mm. First we have
to calculate the surface of one tube. The result we obtain is the surface per piece.
This value has afterwards to be multiplied with the total number of pieces to get the
final result. The surface per piece is different to the surface of one piece (different
units). For example the surface per piece cannot be added to the surface of one
piece. To create this new unit of surface per piece we introduce the new unit “piece”
(see Fig. 2.12).
2 Working with Physical Quantities: Problems and Solutions 45
In some calculations you need two variables with two different units for storing
one physical value, for example, the temperature in thermodynamic calculations,
where the variable t (or θ) stores the numerical value in degrees Celsius and T in
Kelvin. There are tools built into Mathcad in order to eliminate this undue dichotomy.
46 V. Ochkov and K. Orlov
Without this integrated tool, the calculation of the human heart (Fig. 2.10) would
be more complicated and it would be necessary to introduce operators, which
translate the units (mL and torr) in a more popular unit system, for example, in the
international SI-system: milliliter (mL) → cubic meter (m3), millimeter of mercury
(torr) → Pascal (Pa), etc.
There is another source of possible difficulties when using units. There are, if I
may say so, psychological and linguistic reasons related to the unclear definition of
physical quantities and units of measurement in everyday speech and even in
scientific papers. If two people may have a different understanding of one
expression, similar discrepancies can also occur in human-computer “conversa-
tions”. We say “light year”, where the name intends the unit of length—not time,
etc.
We continue with a problem about a pump. Thus, Fig. 2.13 shows an emergency
stop mechanism, if the source data contains an error: here was to variable p1
mistakenly assigned the unit of mass, not pressure.
But in order to avoid the emergency stop of the calculation, it is not only
possible to fix the error (p1 := 1 kg), but also to introduce a new error; (p2 := 20 kg;
see Fig. 2.14).
Of course there will also no emergency stop occur, when using no units at all
(see Fig. 2.15).
In principle, all dimensions of the control mechanism should be understood, so
that in the formulas the units occur only with certain dimensions—the examples
shown in Figs. 2.14 and 2.15 must be completed, so that they stop the calculation,
as shown in Fig. 2.13. Another possible problem is the possibility of mixing up the
input parameters p1 and p2.
The essence of this problem and its possible partial and complete solution, we
illustrate in a more common and understandable task—the creation of a function
with the name of Vc, which calculates the volume of a cylinder. In Fig. 2.16 shows
the function in Mathcad.
The figure shows that if the variables d and h have the dimensions of length, the
function will return the value of volume of the cylinder in m3. But the same
trouble-free function takes also the value of its arguments with any units or even
without any units (Fig. 2.17).
To find out how to get out of this situation, remember the wonderful principle
formulated, which Fourier presented in his classic work “The Analytical Theory of
Heat”, which was published in 1822. This principle is now called the “principle of
dimensional homogeneity” and it states that if the dimensions of each term on both
the sides of equation are same, then the physical quantity will be correct. To
illustrate the difference between the usual algebraic and “physical” equations, we
show an example given in the seminal work of Bridgman (“Regular and Chaotic
Dynamics”, 2001). There we analyze the problem of a body falling, Bridgman notes
that there are at least two equations linking the distance (s), speed (v), time (t) and
the acceleration of gravity (g):
v ¼ gt
s ¼ gt2 =2:
If we now add both as purely algebraic equation, we obtain the following equation:
v þ s ¼ gt þ gt2 =2
The figure shows an error when calling a function that returns the volume of a
cylinder, which may occur if the function arguments are reversed: instead of the
diameter, we enter the height of the cylinder. This error can be avoided by using
further monitoring tools. Here the diameter of the cylinder and its height are two
different physical quantities with the same length dimension. First, in the most
complete form this idea, how to control the inputs, has been effectively developed
in the book Huntley (Huntley G. Dimensional Analysis. Mir, Moscow, 1970),
which, in particular, suggested the use of a “vector of dimension”.
Thus, the variables have to be connected to the dimensions (length–meters, feet,
miles, etc., weight–pounds, pounds, etc.), and to the physical quantities. It should
be assumed that the two variables (diameter of the cylinder and its height) have two
different units of measurement, and their “difference”, in particular, should be
reflected in the fact that these parameters cannot, for example, be added.
In Mathcad are eight dimensions integrated: length, time, mass, current, tem-
perature, light intensity, the amount of substance, and currencies (only in the 14th
and 15th versions and Mathcad Prime). In order to avoid the possibility of mixing
up the input parameters we introduce two different lengths. The diameter is con-
nected to the unit m-d (m = meter, d = diameter) and the height is connected to the
unit m-h (m = meter, h = height). Basically the new unit m-h has not the unit meter.
Therefore we divide in the function Vc the input parameter by his unit m-h.
Afterwards the function is multiplied with m3 in order to obtain the correct result
with correct units (Fig. 2.20).
When you now swap wittingly or unwittingly the arguments of the function the
calculation is interrupted by an error message, which can be described as follows:
“The inputs are swapped”.
There are also real thermodynamic problems, which can be solved with
Mathcad, whereby we need or, at least, it would be desirable to handle different
physical variables, which have the same dimension.
You need to calculate the thermal efficiency of a combined-cycle (see Chap. 14).
To solve this problem, you need functions, which return the thermodynamic prop-
erties of the working fluids. These functions can be integrated to Mathcad by using
the already mentioned WaterSteamPro software (www.wsp.ru)—see Fig. 2.21.
50 V. Ochkov and K. Orlov
Figure 2.21 shows how you can calculate the boiling point of water for a given
pressure by using one of the functions from WaterSteamPro—function
wspTSP. Moreover, all values are non-dimensional, but it is meant that the pressure
is expressed in pascal, and the temperature is measured in kelvin. All functions
from WaterSteamPro by default (without referencing to “watersteampro.xmcd”) use
dimensionless arguments but expressed in SI units or their mixture.
But Mathcad software supports units, and this mechanism is necessary to solve
our problem and indeed in many other heat and power calculations.
In order to be also able to use units for WaterSteamPro it is necessary to make a
reference link to watersteampro.xmcd, included in the WaterSteamPro (Fig. 2.22).
What contains this file? There are firstly some additional units definitions, as
well as various constants needed for calculations in heat and power engineering (the
universal gas constant, critical points of water and steam, etc.). Secondly, and most
importantly, it makes units available for WaterSteamPro.
In the programs Mathcad Prime 1 and 2, WaterSteamPro cannot be integrated
through the mechanism of DLL. Such an opportunity came only in the third version of
Mathcad Prime. Figure 2.23 shows how to call a function named wspPST, which
returns the saturation pressure of water and steam as a function of temperature.
Calling this function in Mathcad 15 is shown in Fig. 2.21 (“dimensionless” call) and
Fig. 2.22 (“dimensional” call through the supporting document watersteamepro.
mcd). The mechanism of function redefining shown on Fig. 2.24 on the example
made im Mathcad 15 won’t work in Mathcad Prime 3 environment. WaterSteamPro
functions, that ypu can see on Fig. 2.23 must be “dimensionless” or supplemented by
the corresponding basic units (such as Pa for pressure). We will examine this methods
more closely in the Chap. 9 on the example of nuclear power plant calculation.
However in Mathcad 15 it isn’t necessary to insert reference on watersteampro.
xmcd (which redefines all the functions of WaterSteamPro program), instead, it is
possible to manually redefine only those functions that are used in the users’ project
taking into account, that different physical values might have same units. For
example in Chap. 14 (calculation of the thermal efficiency of the CCGT unit) reader
will encounter with different physical values—mass of the gas and mass of the
water/steam, which measures in one and the same units of mass.
This is one of the key aspects of working with units.
In thermophysical calculations reader may face values that have identical unit,
but different physical sense. If this fact doesn’t taken into account, it may lead to the
errors during calculations in the Mathcad environment, which we mentioned before.
52 V. Ochkov and K. Orlov
Example In the binary thermodynamic cycle (CCGT unit for example) among the
other values appear values such as specific enthalpy of the first working medium
and specific enthalpy of the second working medium. Both of these values are
typically measured in kJ/kg, but they have different kilograms. To avoid mistakes of
mixing different physical values in such calculations, it is recommended to use
built-in kilograms for one of the values and for the second assign basic SI unit, that
isn’t used in users’ calculation (for example candela—cd). In that case value m, that
usually records ratio of the first working body flow rate to the flow rate of the
second in binary thermodynamic cycles, will be dimensional value, and it will help
to avoid some typical mistakes.
The technology of entering different physical values with same units into cal-
culation shown on Fig. 2.24.
As you can see at the top of Fig. 2.24 introduced two units of mass: for the
measurement of water vapor and gas—the two working fluids of the
combined-cycle plant. Kilograms of water/steam equal to the built-in Mathcad unit
of mass (kg). Kilograms of gas equal a specific unit, which in this calculation is not
used—the unit of luminous intensity cd (candela).
The next stage in the calculation in Fig. 2.24 is the input of initial data. Notice
that in values of mass flow rate are used kilograms of water vapour and kilograms
of gas for different physical quantities that have the same mass dimension.
2 Working with Physical Quantities: Problems and Solutions 53
This allows user to have in his calculation the value of m (the ratio of the flow
rate of one of the working medium to the flow rate of another working body) not
dimensionless, but dimensional. This helps to avoid some errors in the calculation
of binary power cycles, one of which is shown in Fig. 2.24: if we had not intro-
duced two units of mass for the two working fluids, the operator Δhst + Δhgt would
not have been interrupted by the error message, and would have given an incorrect
result. The correct result (specific work of the combined cycle) is calculated from
the last two operators shown in Fig. 2.24.
54 V. Ochkov and K. Orlov
Fig. 2.25 Demonstration of the toolkit for creating different quantities with the same unit
One solution of Mathcad problems with different physical quantities of the same
dimension can be offered through an additional tab “Quantity” in the dialog box
“Properties”, which occurs if we press the right mouse button on the name of the
newly created function or variable. In the text box of this tab can be written the
desired dimension/index—the base (weight, volume, etc.) or component (air, gas,
etc.). Thus the name can be supplemented by the numbers 1, 2, 3, etc., dividing one
from the other physical value. This technique is displayed in Fig. 2.25.
In the calculation shown in Fig. 2.24 dimensionless variable m (the mass flow in a
gas turbine divided by the mass flow of steam in a steam turbine) is transformed into a
(“pseudo”-) dimensional variable. And this is, as we have seen, good: the less
dimensionless quantities, the more monitoring tools can be used. But in terms of
dimensionless quantities we also have other ones—the degree of dryness of steam (dry
weight of steam referred wet steam) and the three values of thermal efficiency (η).
2 Working with Physical Quantities: Problems and Solutions 55
If the numerators and denominators in these fractions have the same dimension,
these ratios must be dimensionless. If the numerators and denominators of these
fractions are regarded to different physical quantities, these values (m, η, etc.)
should be “pseudo”-dimensional. It is possible to use, for example, the technique
shown in Fig. 2.25.
Chapter 3
Concentration of Solutions
Valery Ochkov
Abstract In the chapter you will learn how to use the value of the concentration of
a solute in aqueous solutions in computer calculations, how to convert the species
concentration using modern computer tools.
1
Aqueous solution of this salt is often used as a heat transfer, or rather, a cooling medium in
conditioning systems (see Chap. 17). Adding a salt into the water results in decrease of the freezing
temperature—see Fig. 3.6.
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
2
Usually “measured in moles” is added here. Why?! An expression like this: “Amount of sub-
stance taken—200 mg” is often found in the literature. Certainly, here should be written “mass”
rather than “amount”. These terms are often misused in everyday life and in science—hence it is
necessary to clarify.
3
It was used in the first Mathcad document of the book, see Fig. 1.1.
4
Mass fraction is mass of dissolved solids divided by mass of the total mixture (ω), which equals
to number of dissolved solids measured in grams in 100 g of solvent (P). In this case it is important
to distinguish between this type of concentration and another one—mass of dissolved solids
divided by solvent mass.
3 Concentration of Solutions 59
Fig. 3.3 Direct and inverse functions of density and concentration (mass percentage) of NaCl
aqueous solution
calculation shown on Fig. 3.4 are differ from those containing in various reference
books, online and paper. The fact is that formulas in the reference books sustain
some coefficients (1000, 0.01, 10, etc.) connecting calculations with suitable units.
This peculiarity of Mathcad calculations with units is described in Chap. 2.
Calculation of mass fraction of the second solution is more difficult. It is required
first, access to the function returning density of solution by its mass fraction. This
was done in the calculation shown on Fig. 3.1 by reference to the cloud file
containing such a function. Second, one can calculate mass fraction of the second
3 Concentration of Solutions 61
solution only step-by-step because this value (ω2) contains in both, left and right,
parts of the equation connecting mass fraction and molarity. It was done in a
program using while cycle, see Fig. 3.4. Although, such the method of sequential
approximations can be performed without programming, in a manual mode, as
shown on Fig. 3.5.
Figure 3.5 fixes such sequence of steps: first approximation to ω2 (10 %) is
assigned and then it is refined and transferred into the first (previous) approxima-
tion. The procedure is performed in such a way until the latest couple of values
becomes closely equal (also see Sect. 6 of Introduction).
Having defined masses and mass percents of initial solutions it is easy to cal-
culate required mixture parameters. In the calculation shown on Fig. 3.1 it was
comprised and solved system of two equations. The most interesting is that the
mixture volume is not equal to sum of volumes of initial solutions. We know it
from “vodka phenomenon”: 50 ml of ethanol plus 60 ml of water give volume little
less than 100 ml. However, this phenomenon is explained by dependence of
solution density from its concentration and by other reasons.
As was noted above aqueous solution of NaCl is used as a heat transfer (cooling
medium) because temperature of such solution can be less than zero degrees
centigrade. Figure 3.6 shows the Internet site for calculation of freezing temperature
for some aqueous solutions (including NaCl).
A cold radiator in winter is awful. Although, it can be very useful in summer
heat because it substitutes an air conditioner in a room. Such conditioning system
have already been implemented. Certainly, there have been used special heat
exchangers with blowers rather than radiators. In winter it is fed by warm NaCl
solution, in summer—by NaCl solution at sub-zero temperature.
62 V. Ochkov
The calculation shown on Fig. 3.1 can be drastically simplified using function
shown on Fig. 3.7. As the function shown on Fig. 3.3 it also returns density of NaCl
solution depending on its concentration. However, argument of the new function
3 Concentration of Solutions 63
can be not only mass percent of sodium chloride but also other types of concen-
tration, namely, molality, molarity (see Fig. 3.1), and titre of the solution. Titre is a
mass concentration expressed in g/mL or g/L. The key element of the function
shown on Fig. 3.7 is Mathcad built-in function SIUnitsOf returning dimension of
its argument. It allows us to define type of the concentration inserted by a user in
program as an argument of function ωNaCl(c) and make corresponding
recalculations.
Function shown on Fig. 3.7 can be simplified—it is possible to omit cycle
defining of mass percent by molarity or titer adding three strings to matrix M that
will hold discrete values of NaCl solution concentration expressed by molality,
molarity, and titer.
Fig. 3.7 Calculation of density of NaCl solution for any type of concentration
64 V. Ochkov
Chemists in Russia and in Europe have left off measuring hardness of water in
degrees5 and substituted meq/l. Using SI system is not just recommended but
legitimized almost all over the world.
For example, according to SI system demands we cannot say “calcium con-
centration in solution equals to 1 eq/L”. There is no more such a unit—eq. Also, we
cannot say “calcium concentration in solution equals to 1 mol/L”. We should
specify type of concentration—molar or normal (to be more exact, molar concen-
tration of equivalents), that is, note what is operational unit in a given solution,
calcium ions or charges of calcium ions (cations).6 Omission may result here in
significant mistake: charges of calcium ions in volume unit twice as much as
calcium ions7: for example, if molar concentration of calcium in solution equals to
1 mol/L, normal (equivalent)—2 mol/L. Old but holding ground method of dis-
closure of this omission is to use two groups of amount of a substance: moles (to be
more exact, g-mole, mg-mole, gram-molecular, etc.) and equivalents (g-eq, mg-eq).
They are similar by physical and chemical essence but, as a rule, differ in values.
Although, sometimes is specified what type of concentration is meant. But it does
not belong to do this if it is clear from concentration unit, g/L or mole/L. Sometimes
disagreements transfers into denominator, too. It is considered that we cannot put
liters (capacity unit) into denominator, we should change them for dm3 (volume
unit, length unit in a cube). Indeed, concentration is something divided into volume,
rather than capacity. In thermal engineering, where water (plus steam) is the main
working medium, water hardness is often related to solution mass (meq/kg) rather
than volume (meq/L). Water volume varies greatly with heating (particularly with
evaporation), but mass doesn’t.
SI system eliminates units ended with equivalents and reserves only mole and
another way to note operational unit considering in concentration definition. Today,
as was mentioned above, we cannot just say “calcium concentration equals to …”
specifying operational unit of concentration—mol/L or eq/L (meq/L). It is required
to note additionally type of concentration: molar (ionic) or equivalent (normal,
“charge”). In other words, SI demands, allowing only mole and excluding equiv-
alents transfer specification of operational unit from dimension to name of physical
(more accurately chemical, physical and chemical) value.
5
Many of properties were measured in degrees earlier: hardness of water, strength of spirits,
hardness of metals, liquid viscosity, etc. SI system almost completely excludes degrees from
metrology, which were clear only to focused specialists, and change them to more “physical” units.
Now speaking of temperature we even should say just kelvin rather than degree Kelvin.
6
This thesis can be illustrated for the readers far from chemistry with a geometric analogy. For
example, we cannot say “circle size equals to 20 cm”. It is required to specify what we mean,
diameter of radius of a circle.
7
Diameter of a circle twice as much as its radius (see previous footnote).
Chapter 4
My First Power Engineering
Mathcad-Calculation
Abstract This chapter introduces how to use cloud functions and the
WaterSteamPro software for calculation of thermophysical properties of water and
steam. The thermal efficiency of steam-turbine Rankine cycle calculation is taken as
example. A comparison between spreadsheet software and mathematical package is
provided.
If we have the functional properties of the basic working fluids at hand and
functions for coolant thermal engineering problems—(see Chap. 1) and are able to
work with some dimensional quantities (see Chap. 2), we will be able to solve
easily some tasks—to calculate thermal efficiency (ηt) of a simple steam turbine
power generation cycle.
There is a very useful tool in Mathcad environment—to make links to the other
documents, we already discussed it in Chap. 1 (see Fig. 1.5, for example). It is
possible to make such a link to a document, saved on your own computer, local
network or the Internet in Mathcad 15. The command is shown in Fig. 4.1 for
Mathcad 15; by doing that you can make a link to a very interesting and useful
Mathcad-sheet named H2O, with extensions xmcdz, xmcd or mcd saved on the
Internet at http://twt.mpei.ru/tthb.
After such a reference in the Mathcad-worksheet functions will become available
(visible, as programmers say), which calculate the properties of water/steam,
depending on the parameters of the working fluid (see upper part of Fig. 4.2).
When calculating the thermal efficiency of a steam turbine cycle we need the
following functions with the prefix wsp (see in Fig. 2.21):
• wspHPT—specific enthalpy (H) of water/steam as a function of pressure
(P) and temperature (T);
• wspSPT—specific entropy (S) of water/steam as a function of pressure (P) and
temperature (T);
• wspTSP (or wspTSatP)—temperature (T) of water/steam at saturation line
(S/Sat—here saturation) as a function of pressure (P);
• wspXTS—dryness factor (X) of wet steam as a function of the temperature
(T) and the specific entropy (S);
• wspHSWT (or wspHSatWT)—specific enthalpy (H) at the water/steam satu-
ration line (S/Sat) of water (W) as a function of temperature (T);
• wspSSWT (or wspSSatWT)—specific entropy (S) at the water/steam saturation
line (S/Sat) of water (W) as a function of temperature (T);
• wspTPS—temperature (T) water/steam as a function of pressure (P) and
specific entropy (S).
Having at hand the above functions, it is easy to calculate the specific enthalpy
of water/steam at various points in the steam turbine (Rankine) cycle, the specific
work of the turbine (lt − t, turbine), the specific work to drive the feed pump (lp − p,
pump), the specific heat needed in the steam boiler (qb − b, boiler), and, finally, the
desired value of the thermal efficiency (see end of Fig. 4.2).
4 My First Power Engineering Mathcad-Calculation 67
Fig. 4.2 The calculation of the thermal efficiency of a steam turbine cycle in Mathcad 15 by using
the “cloud” functions
The calculation, shown in Fig. 4.2, can only work if the user’s computer is
connected to the Internet. If this is not the case, the link shown in Fig. 4.2 will turn
red and we get an error message (see Fig. 4.3).
There are two ways for correcting this emergency situation, which was recorded
in Fig. 4.3.
Firstly, you can upload the computer program WaterSteamPro, which we have
already mentioned in the Chap. 2 (see Figs. 2.21 and 2.22). File “H2O.xmcdz” was
designed to provide the subset of WaterSteamPro functions. And functions from
“H2O.xmcdz” can be used without installation of the WaterSteamPro. However
functions from WaterSteamPro can be used not only in Mathcad. For using
WaterSteamPro functions with units in Mathcad you only need to insert a calcu-
lation link to the file watersteampro.xmcd, supplied by the program WaterSteamPro
(In Fig. 4.3 you can find a pair of operators: assignment of units of specific enthalpy
(kJ/kg) and entropy [kJ/(kg K)]. Using such a complex variable name is done
through <Shift> + <Ctrl> + <k>).
If you click twice on a link to the “cloud” of the function (in Fig. 4.3 it does not
work, because there is no connection to the Internet), it will be downloaded and
opened the Mathcad-sheet, the beginning and the middle part of which is shown in
Fig. 4.4.
68 V. Ochkov and K. Orlov
Fig. 4.3 The calculation of the thermal efficiency of steam turbine cycle with using the functions
of the WaterSteamPro
The following figure shows the call of the function wspHPT, calculating the
specific enthalpy of water and steam as a function of pressure and temperature. The
function is hidden and password protected in a separate area. This and other
functions are created on the recommendations of the International Association for
the properties of water and steam (www.iapws.org), in which the author of this
book works.
A downloaded file (named H2O.xmcdz) can be saved on your computer or on
your local network, and you can make a link to it, as shown in Fig. 4.5.
With the functions, shown in Fig. 4.2, it is easy to create a function, that returns
the thermal efficiency ηt of a simple steam turbine power generation cycle as
function of the pressure p0, temperature T0, and condenser pressure pcond (see
Fig. 4.6).
The intermediate results are presented in Fig. 4.2, which can be used for
debugging the function in Fig. 4.6.
4 My First Power Engineering Mathcad-Calculation 69
If you put the Mathcad-document with its functions shown in Fig. 4.6 for
example, on your desktop, then the file can also be referenced (see Fig. 4.7).
Figure 4.7 does not provide a link to the file H2O.xmcdz—it is saved in the file
named Eff-Ideal-Simple-PTU.xmcdz.
The technology of creating user functions and save them on your computer, local
network or on the Internet is an effective technique for solving engineering prob-
lems including thermal engineering as we have already mentioned in Chap. 1.
Having at hand such functions it is also easy to build diagrams for thermal engi-
neering cycles [34]. In the book the reader can find a calculation of the thermal
efficiency of a simple cycle steam turbine plant (STP) with superheated steam,
made in the environment SMath.
Many engineers, scientists and technical workers perform calculations in the
spreadsheet-program MS Excel™, which was mentioned in the introduction.
Figure 4.8 shows a table in MS Excel, to which is attached
WaterSteamPro functions by DLL technology. The list of the functions can be seen
in the dialog box, shown in the same figure (see also Fig. 4.3). It is possible to
calculate the efficiency of the steam turbine cycle using the formulas shown in
Fig. 4.2. The initial data is entered in cell B2, B3 and B4. The remaining cells in the
column are filled with formulas, which carry out calculations. These formulas are
shown in Fig. 4.9. This spreadsheet is made visible after entering the command
Show Formula, tab Formula.
Comparing Figs. 4.2 (Mathcad), 4.8 and 4.9 (MS Excel), you can once again see
that the work with a spreadsheet-program is much more difficult and complicated
than the work with scientific calculator Mathcad.
There is a tool for simplifying formulas in MS Excel, fixing each entry names to
its values, which are stored in the cells. Figure 4.10 shows the Excel-calculation
displayed in Figs. 4.2 (Mathcad), 4.8 and 4.9 (Excel), with the difference that in the
formulas of column B instead of the coordinates of the cells (B5, D2, and so on),
the names of these cells (h0, p0, and so on) appear, which have been pre-assigned
by the command name of the range.
So the result of the table cell B5 in the calculation shown in Fig. 4.10, has a
name h0 (subscript here invalid) and stores the specific enthalpy of steam entering
the turbine. This parameter is calculated using the wspHPT(p0, T0), the arguments
of it are stored in the cells D2 and D3 named P0 and T0. This association between
cells is further displayed through the command Trace Precedents, tab Formula.
In these calculations (and they can be quite complex) are not only formulas used
in the cells, based on numeric values stored in other cells, as shown in our MS
Excel-tables, but also software modules written in the programming language VBA
(Visual Basic for Application). These calculations are practically impossible to
understand. But inserting them into your Mathcad-calculation is possible and
necessary. How it’s done—see below.
Some of the names of table cells, shown in Fig. 4.10, have an index: h1, s2 etc. It
is done due to the fact that the cell names h1, s2, and so on, coincide with the
“coordinates” H1, S2, and so on (H1—cell in column H and row 1).
4 My First Power Engineering Mathcad-Calculation 71
Fig. 4.7 Reference to the calculation of the efficiency of a steam turbine cycle
72 V. Ochkov and K. Orlov
Fig. 4.8 The calculation of the efficiency of steam turbine cycle with MS Excel
Fig. 4.9 Calculation of the thermal efficiency of a steam turbine cycle with MS Excel: showing
the hidden formulas
4 My First Power Engineering Mathcad-Calculation 73
Fig. 4.10 The calculation of the thermal efficiency of steam turbine cycle with MS Excel: working
with named cells
There is a very useful tool for inserting MS Excel sheets in Mathcad. This tool is
used in two cases: firstly, when there is a ready-made calculation or a fragment of a
calculation in MS Excel which you do not want to rewrite for Mathcad, and
secondly, when you want to use in your Mathcad-document any specific tool from
MS Excel—when, for example, you need to build a pie chart, which is not possible
in Mathcad.
Figure 4.11 shows the calculation of the thermal efficiency of our simple steam
turbine cycle in Mathcad with embedded MS Excel-components, which is dis-
played in Fig. 4.8.
Figure 4.12 shows a diagram of a simple steam turbine plant, for which is
performed the calculation. Such schemes are given a lot in this book. They are
usually drawn with the help of a graphical editor (using Paint, for example, a
component of Windows), and inserted into a Mathcad document. The knowledge of
schemes of thermal devices and their drawing “by hand” or a by computer—is a
form of training and testing of students and energy specialists to raise their quali-
fication. Alexey Ochkov developed a specialized program for the automation of the
work [35]. Please visit http://twt.mpei.ac.ru/CALCULON [35–37].
74 V. Ochkov and K. Orlov
Figure 4.12 shows the window for testing the knowledge loaded with a blueprint
(lower part) and a diagram drawn by the user (upper part). There is a direct link to
this demo on the Internet at http://twt.mpei.ac.ru/CALCULON/TWTSchemeQuest_
BOOK1/TWTSchemeQuest_BOOK1.htm. To run the demo, you must install a
special multimedia shell TVT Shell, a reference to the installation of which (http://
twt.mpei.ac.ru/twtshell) will automatically be offered on the site.
After drawing the user schema and clicking Ready the program compares with a
special “grading” algorithm the scheme drawn by the user with the standard. In the
case of error detection, the program will automatically indicate what is missing or
what is messed up. For example, Fig. 4.13 shows the result of the program in the
case of forgetting to draw a feed pump.
4 My First Power Engineering Mathcad-Calculation 75
Fig. 4.13 Comparison: user-drawed steam turbine circuit compared with original version
76 V. Ochkov and K. Orlov
Fig. 4.14 One of the heat engineering schemes from the “Electronic Encyclopedia of Energy”
The comparison of the two schemes is the most important part of this program.
A detailed discussion is beyond the scope of this book, but for readers, who are
interested in it, the following information may be useful.
As a basis for the comparison algorithm the original Gemini II algorithm [36] is
taken. The family of Gemini (The Twins) algorithms is based on the iterative
calculation and comparison of the codes of individual elements and “networks” of
the graph. It is widely used in the electronics industry to solve the LVS
(Layout-Versus-Schematic, comparison of schematic diagrams and final PCB
design or chip) problem.
In the process of developing an algorithm several changes in the original
algorithm have been made, in particular, the problem of the impossibility of using
“isopotential networks” has been solved, because sometimes the comparison of
technological schemes requires an exact match topology (piping, pumps, valves,
and so on).
The program for testing the knowledge of schemes is a part of the author’s
“Electronic Encyclopedia of Energy” [37], (www.trie.ru), where you can find a lot
of heat engineering schemes. One of such scheme is shown in Fig. 4.14.
Pressing the round button, shown in the lower part of Fig. 4.13, you can select
individual parts of the heat power station. In addition, from this multimedia scheme
you can quickly jump to the corresponding description in the encyclopedia.
The authors express their deep gratitude to Vsevolod Znamensky who wrote the
main part of the H2O Mathcad-sheet.
Chapter 5
“Cloud” Thermal Engineering Algorithms
Abstract This chapter explains how cloud algorithms can be used in the heat
power processes calculations for the thermal power plants, for example, the pro-
cesses of expansion or throttling of the steam, compression water in the
pump. Information technology for the Internet libraries of equations is provided.
For thermodynamic calculations one shall know not only properties of working
fluids (water, steam, air, flue gases etc.) [23, 38, 39, 40] and of thermal engineering
related materials, but also of calculation formulas.
These formulas or sets of formulas with instructions for their application (for-
mulations are taken from various reference literature sources, for example, from
[41] or from Internet and are used for relevant calculations. Such work method is
not quite convenient as it hinders execution of design calculations, and enhances
occurrence of error risk in them. In view of the above, and also considering the fact
that at present PC of almost all specialists making thermodynamic calculations and
other engineering calculations are provided with continuous and high-speed
Internet, we have developed new technology for work with formulas. It will be
demonstrated with the use of simple standard examples typical for the engineering
office (calculator) Mathcad [42], which is quite popular among the students,
engineers and research engineers.
“Cloud” technology allowing work with formulas is a development of approach
described in Chap. 4 [40] and related to use of function for calculation of energy
working fluid properties contained in Internet. We’ll illustrate work with such
functions by simple, but at the same time, quite typical examples.
When designing any facility, after introduction of initial data one shall put the
link reference to Mathcad-document (contained in Internet (in the “clouds”) at the
address http://twt.mpei.ru/TTHB/H2O.xmcdz) into the calculation (see Fig. 4.1 in
Chap. 4). After making this link reference all functions defined in
Mathcad-document to which the link reference is made will become available in the
working document. These are the functions converting heat-transfer properties of
water and water steam to the following: specific enthalpy (h), specific entropy (s),
temperature (T), dryness factor (x) etc. [23] as well as functions converting
parameters (calculations) of some thermal processes (throttling in control valve or
steam expansion (work) in steam turbine) into values of specific enthalpy and steam
dryness factor at the end of the above processes. Names of these two functions have
radicals in the name EXPANSIONPTPEFF1: EXPANSION-expansion,
PT—steam pressure (P) and temperature (T) in the initial point, P—steam pressure
in the end point and EFF (efficiency)—internal relative steam expansion process
efficiency.
This approach is an alternative to application of dedicated software tools for
calculation of working fluid properties [43], which can be downloaded and installed
on PC for use. There is no need in downloading and installation of such software
tools on the user’s part if one uses “cloud” functions. As a clear benefit of “cloud”
technology we can mention capability of on-line function updating: as a result the
user always uses the most current calculation version.
Among the drawbacks of “cloud” functions we can mention the need for Internet
access when making calculations. In this case as an alternative we can recommend
usage of standard software tools. Thus, in the absence of Internet access one can
install the software tool [43], and it will work for the examples provided in the
present book. It is provided by similar names and arguments of functions used for
calculations of working fluid properties. As another option we can recommend
downloading of Mathcad-file under the name H2O.xmcdz, with its further instal-
lation on PC (for example, on the desktop, with link reference to it from the
working calculation as it is shown in Fig. 5.1).
The use of “cloud” technologies for formulas and formula sets gives the user
access to completed, ready calculations which can be directly apply in design efforts
and which represent solutions of specific tasks.
Task 5.1. It is necessary to calculate water steam throttling in a control valve and
its further expansion in a steam turbine. These processes shall be shown in h, s- and
T, s-diagrams of water and water steam (Fig. 5.1).
1
Let us remind you that the names of all functions are prefixed with wsp. It is an abbreviation of
the WaterSteamPro software which can be downloaded from www.wsp.ru. The “cloud” functions
described in this book are a small part of WaterSteamPro software.
5 “Cloud” Thermal Engineering Algorithms 79
Fig. 5.1 Calculation process throttling and expansion of steam in the turbine
80 V. Ochkov and K. Orlov
2
We set very large pressure difference (24–10 MPa) in order to show that process on diagrams of
Fig. 5.1. At real Thermal Power Plants pressure loss in a steam pipeline from the boiler to the
turbine accounts for near 5 %. If р2 is to be entered less than р1 for those 5 %, then points 1 and 2
will coincide.
5 “Cloud” Thermal Engineering Algorithms 81
x—dryness factor. It equals to 0 for water and 1—for saturated steam. And second.
In Fig. 5.1 diagram points 1, 2 and 3 were jointed with straight lines. In Fig. 5.2
lines standing for the processes of pressure building-up in the pump and water
heating in the economizer are shown in the diagram the other way: for that purpose
vectors vp (pressure), vh (specific enthalpy), vT (temperature) and vs (specific
entropy) were generated for processes 1–2 and 2–3, which are drawn on the dia-
gram of Fig. 5.2. That’s the process in cases when the line in a diagram is not
straight or close to straight. In principle the similar way shall be taken for diagrams
of Fig. 5.1, but it’s difficult to define signatures 1–2 (throttling with efficiency
equaled to 5 %) and 2–3 (steam expansion in the turbine with internal relative
efficiency of the process equaled to 75 %), as the values of local efficiency might
differ for different areas of the process. Line 2–3 in Fig. 5.2 is not straight. One can
make sure by expanding the diagram area. That is why vectors were used for
building this line in the diagram. The same was done with line 1–2.
In Fig. 5.2 T, h-diagram was built rather than hs- or T, s-diagram. T, h-diagrams
are rarely used for demonstrating steam-turbine cycles, however they are used here
for the following reason. At present students specializing in heat and power engi-
neering and keening on computers have started building not only conventional
“planar” diagrams of steam-turbine cycles (hs- and T, s-diagrams—see Fig. 5.1),
but and 3D, dimensional diagrams. e.g., h, T, s-diagrams [34]. Two planes of such a
dimensional diagram are shown in Fig. 5.1, while the third one (Th)—in Fig. 5.2.
Looking at Fig. 5.2, one could mistakenly assume that the process of water pressure
build-up “gets into” two-phase area to the right of the saturated water line. That’s
not the case. The area to the right of the water saturation line shows two areas of
thermodynamic surface of water in h-, T- and S-coordinates. At point 2 water has
the following parameters: p = 13 MPa, T = 34.4 °C and h = 156.6 kJ/kg. But the
same parameters by T and h are applied for the “twin” of point 2—two-phase
liquid: water with steam content of 0.481 %, temperature of 34.4 °C and pressure of
5.34 kPa.
A Mathcad user can broaden the list of “cloud” Mathcad-functions contained in
the document H2O.xmcdz located at http://twt.mpei.ru/TTHB/H2O, using the tools
of equations and system solutions or zero function searches. So, in the bottom of
Fig. 5.2 creation and calling for wspPTH function is shown which converts
water/steam pressure depending on temperature and specific enthalpy (one of two
backward functions wspHPT of “cloud” functions package located at http://twt.
mpei.ru/TTHB/H2O.xmcdz). The created function wspPTH uses the
Mathcad-inbuilt function root and converts the value of pressure by temperature
and specific enthalpy close to the point of first approximation by pressure (1 MPa).
We do not include the complete list of functions as it is being continuously
enlarged.
And a little bit more about the pump.
In Chap. 2 we gave the calculation of human heart power (see Fig. 2.10) as a
pump transporting incompressible fluid. In that calculation we multiplied the value
of volume flow rate of blood by pressure difference in the heart. Sometimes this
simplified formula is applied by mistake to calculate apparatus transporting
5 “Cloud” Thermal Engineering Algorithms 83
Fig. 5.4 Calculation of water mass flow in MS Excel (demonstration of calculated values and
formulas used; work with connected WaterSteamPro)
Fig. 5.6 Calculation of water mass flow rate in SMath: operation with a function in hidden area
Fig. 5.7 Network calculation of water flow rate (a task which is inverse to the ones shown in
Figs. 5.4, 5.5 and 5.6)
And one more afterpiece. The name of the previous chapter in the book was
“My first thermo-technical Mathcad-calculation”. A book extract given further on
can be called as “My first thermo-technical calculation at Elsevier/Knovel
web-site”.
86 V. Ochkov and K. Orlov
3
A tool of Web of Science scientists publishing activity evaluation is one of developments of
Elsevier.
5 “Cloud” Thermal Engineering Algorithms 87
Fig. 5.10 Home page of e-publishing house Elsevier with a list of online tools
the pump (Figs. 5.13, 5.14, 5.15 and 5.16) relying on the function converting water
density as a function of pressure and temperature, made the start to that work.
The calculation given in Fig. 5.14 was made in SMath which can be freely
downloaded from www.smath.info. And that’s the main point. Knovel did not place
a stake on Mathcad due to the fact that it is a pretty expensive commercial product
of another company—PTC.
In interactive SMath-calculation shown in Fig. 5.14 there is a hidden area with a
function which converts water density as a function of pressure and temperature.
This area is concealed and is partially shown in Fig. 5.15. This function and lots of
5 “Cloud” Thermal Engineering Algorithms 89
other similar functions can be downloaded from the authoring web-site http://twt.
mpei.ac.ru/rbtpp/Region. It has the name wspD1PT there. If 1 is to replace with 2
in the name of this function, it will convert the density of water steam rather than
water and the properties of water steam will be defined by the second area of
formulation IAPWS-IF97 rather than the first one. If there is 3 in the functions
name, it will refer to near-critical region. If there is 5 there, it will refer to high
temperature water steam region. The fourth region of formulation IAPWS-IF97
90 V. Ochkov and K. Orlov
Fig. 5.15 Function of Knovel’s Interactive Equations converting water density (a part of the
function is shown)
defines water and water steam properties at the saturation line between water triple
point to critical point. Mathcad-file with H2O name and extension of mcd, mcdx or
mcdxz being stored in the Internet at http://twt.mpei.ru and widely used by us in
this book (see example, Fig. 5.5) includes functions with names wspD1PT,
wspD2PT, wspD3PT and wspD5PT. Besides, there is a function there with name
wspDPT, which defines the formulation region IAPWS-IF97 with parameters
(pressure and temperature) and calls the required numbered function.
One can also see a list of menu commands DOWNLD in Fig. 5.14 allowing
downloading files of calculating software SMath (Knovel Worksheet), Mathcad and
MS Excel. Figure 5.16 shows a respective Mathcad-calculation of water pump
power.
Its efficiency appears in the formula of pump power calculation. It’s easy to
prove without any measurements that this value is always less one (one hundred per
cent). It remains only to come closer to any running pump (compressor, turbine,
etc.) to hear its noise and to feel its vibration. A part of pump drive energy is used
for that noise and vibration. A perfect pump must run silently and without the least
vibration. Noise is transmitted by air, metal and water after all. Noise level is
measured in decibel which we mentioned in Chap. 2 (see Figs. 2.10 and 2.11).
Noise problem at power plants and other industrial facilities is a big problem. The
author’s server includes a section with examples of calculations related to noise
control at power plants (see http://twt.mpei.ac.ru/TTHB/1/Tupov.html) [44].
5 “Cloud” Thermal Engineering Algorithms 91
Fig. 5.16 Mathcad-calculation of pump power downloaded from Knovel’s Interactive Equations
By the way, about acoustic speed in water. Water and water steam properties are
split to thermodynamic (enthalpy, entropy, internal power, etc.) and physical
(viscosity, thermal conductivity, etc.). The term “thermal and physical properties of
water” unites these parameters. So, acoustic speed in water is neither physical nor
transport parameter as one might think at first sight, but a thermodynamic param-
eter. One can say the same about water density.
There are future plans to place functions of working media, heat carriers and
materials of heat, nuclear and industrial power engineering in Knovel’s Interactive
Equations. In particular, the data base of IVTANTERMO on thermal and physical
properties of individual substances [45]. As it is, e-version of this guide with online
calculations and a downloading option of Mathcad 15-, Mathcad Prime-, Mathcad
Express-, SMath- and Excel-files can be found at the authoring web-site http://twt.
mpei.ac.ru/TTHB/2/OIVT/IVTANThermo/Rus/index.htm (Fig. 5.17).
Properties of individual substances will be exemplified with their use similar to
what have been done to water density during calculation of pump power.
Technology of inverse calculations will be used when, e.g., water density and its
temperature are known and pressure is to be found [46].
Figure 5.18 shows a page of a web-site with calculation of heat capacity of
gaseous ozone and image of a working point at the diagram. Scaled Gibbs energy,
entropy and entropy difference for that substance are also calculated at the same
web-site. Formulas used in calculation can both be reviewed (that can be done by
referring to paper guide-book, of course), copied and pasted to your own calcula-
tion. The only thing to do is to download MS Excel table (see link to it in Fig. 5.17)
92 V. Ochkov and K. Orlov
Fig. 5.18 Interactive network calculation of specific heat capacity of gaseous ozone
with that formula in text format ready to carry over to programs created by pro-
gramming languages FORTRAN, C, Pascal, BASIC, etc. This is a specific example
of “open interactive network calculation” principle realization which is described in
this book.
Properties of individual substances at Elsevier/Knovel web-site (see Figs. 5.10,
5.11, 5.12, 5.13, 5.14, 5.15 and 5.16) shall be shown with reference to their use the
same way as it was done to water density in calculation of the pump power (see
Fig. 5.16). Technology of inverse calculations will be used when, e.g., water
density and its temperature are known and pressure is to be found [46].
Chapter 6
Symbolic Mathematics and 3D-Graphics
in Thermal Engineering
Valery Ochkov
In Chap. 4 we calculated the thermal efficiency of the simplest ideal steam turbine
cycle. Chapter 5 demonstrates how to calculate the “non-ideality” (irreversibility)
steam expansion process in the turbine and in raising water pressure in the feed
water pump. There we saw at the same time how these processes can be displayed
in 2D graphics. Now we complicate the steam turbine plant by adding two steam
extraction points for regeneration to increase the thermal efficiency of the power
plant.
Figure 6.1 shows a Mathcad calculation of the thermal efficiency of a steam
turbine cycle with two mixing regenerative heaters H-1 and H-2.
After the reference to the cloud Mathcad-sheet with the name H2O and entry of
input data (see the cycle diagram and its parameters), we introduce into the cal-
culation a system of two equations of material enthalpy balance in heaters. By
means of Mathcad’s symbolic mathematics, this is solved for the unknowns α1
(which is the fraction of steam extracted from the turbine to the first heater) and α2
(the fraction of steam extracted from the turbine to the second heater). After this
analytical solution by the help of the solve operator it is easy to calculate the cycle’s
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
Fig. 6.1 Calculation of cycle thermal efficiency of steam turbine with two steam extractions
6 Symbolic Mathematics and 3D-Graphics in Thermal Engineering 95
Fig. 6.2 Function of cycle thermal efficiency of steam turbine with two steam extractions
thermal efficiency (see Fig. 6.1) by determining the parameters of steam at its key
points and performing simple balance calculations.
You can manually change the turbine extraction pressure (p1ext and p2ext) and
observe how the thermal efficiency of the cycle is changing; although you could
entrust the work to a computer. This requires you to create a user function that
returns the thermal cycle efficiency, depending on its seven initial parameters listed
in the upper part of Fig. 6.1. This function is shown in Fig. 6.2.
96 V. Ochkov
Fig. 6.4 Optimization of the steam turbine cycle with two extractions points
Figure 6.6 shows a complementary function, created using the function shown in
Fig. 6.2. With these two functions: the main function—Fig. 6.2 and complementary
function—Fig. 6.6, we can construct the desired graph and with a tracing tool mark
on it the point of maximum thermal efficiency.
The y-axis of the curve in Fig. 6.5 records the output of the function for two
steam extraction points from the turbine. However with p1ext = p2ext this function is
also suitable for the case of one steam extraction point. Indeed the curve itself is
parametric: an independent variable in calculations and optimization will be based
on the pressure extraction pext, which changes depending on the pressure in the
condenser pcond and the live steam pressure p1, but the curve is based on two
functions of pext.
From the calculations shown in this chapter, it is clear that the optimized single
steam extraction from turbine to regeneration increases the power plant thermal
efficiency from 39.11 % (see Fig. 6.1) up to 40.125 % (see Fig. 6.6). The transition
to the two regeneration steam extractions with pressure optimization improves
thermal efficiency up to 41.372 % (see Fig. 6.4).
On the site of the book you may find calculations for steam turbine units with
three or four regeneration extractions performed by the method described above.
Figures 6.6 and 6.7 show the operators of the symbolic solution of fluid and heat
balance equations of three and four mixing heaters of a steam turbine unit, whose
diagram is shown in Fig. 6.1. This shows only part of the solution since it is too
long and written in the form of a matrix with one row and three (Fig. 6.7) or four
(Fig. 6.8) elements.
98 V. Ochkov
Fig. 6.6 Optimization of steam turbine cycle with one steam extraction point
Fig. 6.7 Solution of the equations of material and heat balances of three mixing heaters of steam
turbine unit
100 V. Ochkov
Fig. 6.8 Solution of the equations of material and heat balances of four mixing heaters of steam
turbine unit
Valery Ochkov
In the previous Studies, when calculating steam turbine cycles, we did not consider
impurities that may be contained in water and steam. Concentration of these
impurities is very small and cannot significantly affect the thermal parameters of the
working medium. But the impurities can significantly affect the performance of
power units.
The initial natural water at thermal and nuclear power plants is treated (desali-
nated, etc.) in order to avoid build-up of deposits in the steam turbine cycle.
Impurities may be transferred from boiler water to steam, and deposited in the
turbine steam path, which decreases the turbine’s efficiency and increases the load
on the thrust bearing,1 and the turbine is fraught with “under-scale” corrosion. The
main sources of boiler water impurities are:
1
To the extent that bearings cannot withstand such above-the-limit loads and fall apart. In times of
old (now we can already say this about the beginning of the era of “big” power engineering), an
operator on patrol in a turbine hall of a power plant used a stethoscope to listen to the bearings and
other “critical” parts of turbines, pumps and fans. Nowadays such “preventive medicine” control is
automated.
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
• make-up water, which is actually carefully treated, and it can bring impurities
into the boiler only due to an oversight of the staff of the chemical department of
the power plant who are involved in servicing the water treatment plants;
• leaks into a turbine condenser, which may bring a “dirty” water from a cooling
tower, a cooling pond (service water supply recirculation system with spray
ponds, which, for example, will be considered in Chap. 18) or from a river (flow
systems) into a “clean” condensate, and then into the boiler and the steam
turbine;
• corrosion processes in steam turbine circuits “enriching” water and steam with
oxides of iron, copper or other metals included in constructional materials
(alloys) of the mechanical equipment.
To reduce the accumulation of impurities in the boiler water, build-up of
deposits in the boiler, and their entrainment into the steam turbine, one should
take-up various measures: purify circuit water in a polishing plant, dose various
antiscale additives (e.g. phosphates) into the boiler, etc. Here, heat engineers come
to help to chemists conditioning make-up water, servicing a polishing plant and
preparing reagents for dosing into a boiler—see the title of this chapter. What is the
essence of this “gift” [47, 48]?
The impurities, entering with feed water in a drum boiler, are concentrated in the
boiler and removed from it with a blowdown, which makes 1–3 % steam flow (the
Pr variable). But a small part of the impurities goes into the steam, and this leads to
the above complications. The ratio of impurities concentration in the steam to
impurities concentration in the boiler water is conventionally called a cumulative
entrainment rate which consists of two different processes—entrainment of
impurities with drops of moisture into saturated steam and dissolution of impurities
in the steam. But blowdown of water from a boiler cannot be smooth, the process
should be staged, i.e. “thermal engineers’ gift to bad2 chemists” is required to be
implemented.
As you know, a gift horse is not to be looked in the mouth, but we shall not only
accept the “gift of the heat engineers” but also try to “look in its mouth” and
optimize it. Figure 7.1 shows the website, by logging on to which you can not only
simulate, but also optimize the process of three-stage evaporation of water in a
drum boiler: Water from the first (pure) chamber is partially evaporated, and par-
tially blown into the second (“dirty” salt) chamber, and from it in turn, again is
partially evaporated and partially blown into the third chamber, which is normally
designed as a remote cyclone separator.3
2
We say “bad” not in the sense of no good, but in the sense that chemists in power engineering at
some point can not, by reasons beyond their control, ensure an acceptable quality of a boiler feed
water and the quality of steam leaving a boiler.
3
The boiler drum is equipped with not only partitions separating the evaporation surface into
chambers, but also with other boiler internals (steam washing screens, moisture separators, etc.),
which are not shown in Fig. 7.1.
7 Thermal Engineers’ Gift to Water Chemists 103
On the site, shown in Fig. 7.1, you can adjust the input data, press the
Recalculate button and get an answer—the concentration of impurities in the steam
at the set proportions of steam generated in the first (pure) chamber (optimization
variable x) and in the second (“dirty”) chamber (y), and the optimum values with
which the concentration of impurities in the steam will be minimal.
This is a classic optimization problem, which includes three elements:
• Objective function—concentration of impurities in the steam to be minimized;
Fig. 7.2 shows a fragment of a Mathcad worksheet, where using the symbolic
operator solve one could solve algebraic equations for balance of impurities in
each of the three chambers and form the objective function of two arguments Cs
with forming the intermediate functions C with subscripts w1, w2, w3 (boiler
water in the three chambers) and s1, s2, s3 (steam generated in the three
chambers);
• Optimization variables x and y—steam fractions generated in the first and
second chambers of the boiler drum; Fig. 7.1 shows that a site visitor can change
the values in the same-name text fields and see that the objective function Cs
returns; with x := 1 and y = 0 we have a smooth evaporation; with x + y = 1 we
have a two-stage evaporation (x < 1), and with x + y < 1—a three-stage
evaporation; a visitor can “play” with the x and y values and watch how Cs
value varies approaching or moving away from its minimum value;
• Constraints; one of them we have already mentioned, this is x + y < 1; the two
other constraints are x ≥ 0 and y ≥ 0; these constraints lead to the fact that the
area of visualization of the objective function has the form of an equilateral
triangle (see Fig. 7.1). More precisely. This area has the form of an isosceles
7 Thermal Engineers’ Gift to Water Chemists 105
4
Figure 7.3 shows a numerical (approximated) search of the maximum of function η(x, y)—the
efficiency of the staged evaporation process.
5
The fact, that among other things, makes this problem interesting is that it can be any number of
optimization variables—from 1 to infinity. But we have to repeat that in practice they are limited to
2–3 evaporation stages. This problem is common with another problem considered in Chap. 6
here.
106 V. Ochkov
If we talk about water treatment for boilers, we can mention here that the author
has developed and placed on the server some calculations relating to the process
(see http://twt.mpei.ac.ru/ochkov/VPU_Book_New/mas/watertreatment.html) and
even opened a special website with a short name www.vpu.ru (vpu means water
treatment plant). On this site in particular, one can find calculations of the properties
of Dow Chemical ion exchange resins (www.dow.com). Information about the
properties of ion exchangers is supplied to consumers in two ways: on CD-ROM as
a “closed” tool to calculate water treatment process diagrams (ROSA, IxCalc etc.),
and as “open” technical documentation on paper or on the Internet sites.
Advantages and disadvantages of these methods for publishing information have
been already discussed here. The first method (“closed” tools on disks) can be
hardly called a publication. After starting the program, located on the disk or
downloaded from the Internet, the user is usually only given the opportunity to
enter the input data and get the final answer: no intermediate results, and much less
the formula used to obtain the answers, are given. Such “closed” information
technology is inconvenient for the manufacturer of ion exchangers: if the developer
of the program will leave the company, then it becomes very difficult, and some-
times impossible to make changes and additions in the program (“nightmare” of
legacy software, which are described in the introduction). It is worth also men-
tioning the fact that these programs are very difficult to master. Even an experienced
expert, well-versed in computers and water treatment, has to spend a lot of effort
and time (several months) to learn to work with programs such as ROSA, IxCalc
etc. After such “mastering”, an expert is often “hooked” on the program and can not
switch to other more up-to-date and flexible design tools. The weak point of the
second method for publishing technical information is that it is “dead”—for it is
impossible to calculate at once, after entering the required input data.
Moscow Power Engineering Institute (www.mpei.ru) and Ltd Triery (www.trie.ru)
developed a process of publishing technical information with the possibility of
making immediate calculations using such information. A visitor to the author’s
websites (http://twt.mpei.ac.ru/TTHB/1/Dow/index.html and http://twt.mpei.ac.ru/
TTHB/1/RH/index.html) is given an opportunity to automatically calculate using
such important process parameters of ion exchange resins as the value of slip of a
removed ion and the working exchange capacity [49].
Small Divertissement after Studies 6 and 7. One of the author’s interests,
directly related to his main job, is collection and animation of numerical methods for
solving mathematical problems. This “collection” concerning optimization is posted
on the website https://www.ptcusercommunity.com/groups/optimisation-with-
mathcad. A special place in the “collection” is given to the examples which can
be described as “obvious and incredible.” One such example, adapted to the subject
of the book is as follows. There are two mines where coal is mined, and two thermal
power plants, where the coal is burned on “just-in-time” basis, i.e. without storage.
The first mine “gives to the surface” 50 tons of coal per day, while the second—70.
The first power plant burns 40 tons of coal per day, while the second—80.
Transportation of coal from the first mine to the first power plant (the variable
m1pp1) costs 1200 rubles per tonne, from the first mine to the second power plant
7 Thermal Engineers’ Gift to Water Chemists 107
(the variable m1pp2)—1600, from the second mine to the first power plant (the
variable m2pp1)—800 and from the second mine to the second power plant (the
variable m2pp2)—1000. It is necessary to organize the transport of coal in such a
way that its cost was minimal. Figure 7.4 shows the condition of the transportation
problem and its “manual” solution: the cheapest route (the second mine—the first
thermal power plant: 800 rubles per tonne of coal transported) is fully loaded (see the
first term in the expression in the bottom in Fig. 7.4), and the remaining coal is
“scattered” on other routes.
Figure 7.5 shows the solution of the optimization problem using Mathcad 15
standard features—form the objective function (ObFunc—the cost of transporta-
tion), provide a first approximation (we took the values found through the manual
solution of the problem, as shown in Fig. 7.4), after the keyword Given the con-
straints are formed, and function Minimize finds the unknowns to minimize the
objective function when the restrictions are applied. Automatic solution of the
problem turned out to be paradoxical—nothing should be transported at the
cheapest route.
The problem of coal transportation—probably the only one in the book which is
more appropriate to be solved in MS Excel environment, rather than in Mathcad’s
environment. After all, this is purely an accounting program where the money
appear. MS Excel spreadsheet are designed specifically for economic calculations,
as noted in the introduction here.
Figure 7.6 shows calculations in MS Excel environment for the problem of
transportation of coal, which repeats the suboptimal “manual” calculations in
Mathcad’s environment, as shown in Fig. 7.4. The objective function (cost of
transportation) is stored in cell B15 and displayed in the formula bar.
The user can change the transport plan (cell B11, B12, B13 and B14) to meet the
constraints for the exit of coal from the mines and bringing it to the thermal power
plant with minimum cost of transportation. But the job could be charged to a
computer by calling a command (dialog box) Solver from the Data tab (Fig. 7.7).
108 V. Ochkov
this he should be guided with some objective function—an average specific fuel
consumption for electricity generation. One can select here either to turn off the
whole power units, or reduce the load on some of them. And on what units? On the
least or the most economical? One might find that reducing the load on the most
uneconomic unit will not be the optimal solution. At the end of the day, transporting
coal at the cheapest route leads to the cumulative rise in price of the coal trans-
portation (see Fig. 7.4)… Another small paradox in this problem, associated with a
number of variables, is mentioned in Chap. 17.
Grand Divertissement after Chaps. 6 and 7. The author showed the typescript
of the book to many people for whom the book could be of interest. One of them
expressed the opinion that the book lacks optimization problems. And our whole
life is nothing but optimization. This is an old problem, which the author put in a
ready book to fulfill the wishes of the reviewer.
The paint problem. Russian writer Mikhail Zhvanetsky is often asked questions
where he takes themes for his sketches. “Look out the window and listen to the talk
on the street”—this is the answer of the great satirist. “And how do you memorize
all this?”—It is another question. “But I can not forget!”
Everyday stories are also worth collecting for writing computer studies, which is
a hobby of the author of this book.
By profession the author is Professor in University (Moscow Power Engineering
Institute—www.mpei.ru), where he gives a course of lectures on information
technology (see the course outline at: http://twt.mpei.ac.ru/ochkov/Potoki.htm) and
related disciplines (http://twt.mpei.ac.ru/ochkov/Potoki_MOpt.htm), and manages a
group of process engineers and programmers, developing training programs and
computer simulators for thermal and nuclear power plants (www.trie.ru). Power
plants and power systems need our programs, but the notorious payments crisis
prevents them from buying that. That’s what a computer study appeared in March
1997 [50].
Joint-stock company “Tambovenergo” without free money, however, has
expressed a desire to buy our software. Kotovsk paint factory (KPF, Tambov
region) needs electricity for its operations. Moscow Power Engineering Institute
needs paint to repair its classrooms. The author’s scientific group needs a new
computer hardware, tools and, of course, the salary. To solve such problems the
humanity since the dawn of civilization has come up with the money. The transition
of Russia from socialism to the market has revived barter. In the above commodity
chain one link was missing to make it closed. Fortunately, MPEI received a batch of
computers, a couple of which, as we have agreed, could be exchanged for paint.
This combination was only part of the described computer composition, if you
remember the interpretation of the expression “chess composition”—solution to the
puzzle by drawing up a chain of moves.
The second part of the computer composition took place in Tambov and
Kotovsk—on KPF factory. In “Tambovenergo”, after take-over of the software, I
(the author proceeds to the story in the first person) was given a power of attorney
7 Thermal Engineers’ Gift to Water Chemists 111
to receive the paint products for the amount of 14 million, of course, old rubles6
against the factory’s debt for electricity and I left for Kotovsk. KPF Sales depart-
ment first flatly refused to let the paint for some incomprehensible debt out there,
and not for real money, but after the threat of blackouts and heat disconnection,
though with difficulty, they agreed. Paint, which suits me, or rather not me, but
MPEI purchasing department, costs 14,600 rubles per liter and is bottled in con-
tainers (jars, if you follow chandlery jargon, which I picked up in Kotovsk) each 15
and 55 l in volume. The cost of an empty jar was 24 and 30 thousand rubles
respectively. An employee of the KPF sales department (her name was Olga), when
typing an invoice on a computer, asked me a question what paint containers I would
take. Instinct of a longtime collector of computer studies immediately prompted that
here lies a typical and, most importantly, the real problem of linear programming,
where the objective function to be maximized—the total amount of paint (or the
cost of the paint), the variables—the number of jars filled with paint, 15 and 55 l in
volume, to be collected, and three constraints:
• the cost of paint should not exceed the 14 million rubles agreed with
“Tambovenergo” one can not take an incomplete jar (constraint for the integer
variables);
• number of jars of different sizes should not be a negative number.
Olga offered to help in solving the optimization problem and then using a
calculator figured that I needed to take 16 large and 2 small jars with a capacity of
910 l of paint for the amount of 13 million 814 thousand rubles. Remembering how
desperately I traded at “Tambovenergo” and yet increased the price of programs
from 12 to 14 million rubles, I asked Olga, if it would be possible not to loose
186,00 thousand—do not leave it in favor of “Tambovenergo”. She said no,
because she had to solve these problems almost every day, optimizing not only the
cost of paint, but its canning in containers of different sizes, and that she was a dab
hand at solving such problems.
Watching the “dance” of Olga’s fingers on the keys of the calculator and the
numbers on its display, I realized that Olga used a so-called “workers’ and peasants’
“algorithm to solve the problem: first she selected paint in a large container, and
then the rest of the money (or the volume of the container) was filled with paint in a
small container. It goes something like this when we pack a suitcase when going on
a trip—first put the big things inside, and then cram every little thing into the empty
space. I asked Olga, why she did not use to solve these problems a computer and
MS Excel spreadsheet, which was put on her computer screen as if on purpose.
I immediately offered to show how it’s done. MS Excel environment has a solver
procedure which dialog box is called by Solver command from the Data tab. In this
window, the user specifies the cell storing the objective function, which is necessary
to maximize, the cells with solver variables (at the beginning of optimization they
6
In the Mathcad-worksheet and MS Excel spreadsheet solving this problem we cut-off three zeros.
112 V. Ochkov
are either empty, or store the values of the first approximation to the maximum) and
the constraints (Fig. 7.8).
The optimization algorithm using the MS Excel solver procedure could be called
“lazy”: the user creates a calculation table and says: “With a wave of a wand do so
that… the objective function took maximum (minimum defined) value, but with
that all the constraints should be met.” To do this, the user only has to click a
7 Thermal Engineers’ Gift to Water Chemists 113
Solver. The MS Excel solver procedure gave us the old result—16 large and 2
small jars. But I did not want to give up.
A good rule—to verify the solution of a problem using not only other methods,
but also other software tools. Also, do not forget about the KISS-programming
principle. This has nothing to do with kisses, although its good attitude to the
problem and the computer is visible. KISS—an abbreviation of the English phrase
“Keep It Simple, Stupid!”. It urges to solve the problem in the easiest way and only
use sophisticated algorithms and techniques when simple methods are not suitable
because of the length of time spent on computations or due to irrational use of other
resources, human and/or computer.
The simplest way to solve the problem using a computer—it is to search through
all the options and stop at the optimum. Fortunately, the options are not so many—
1088: for the allotted 14 million rubles one could not take more than 63 small paint
jars, or not more than 16 large jars. The brute force search could be called “workers’
and peasants’” computer solver method. But among other things it can give a
hundred percent certainty that the found solution is not only correct, but also
unique, or demonstrate that there are several such solutions. A similar situation is
not uncommon in the integer linear programming problems.
So, the brute force method. By following the above rule, a new method for
solving the problem should be combined with a new software tool for its imple-
mentation. This, of course, could be also done in MS Excel environment by
drawing up a table of all solutions and/or writing a brute program using built-in
Visual Basic for Applications (VBA) language. But Olga also had Mathcad on her
computer (the phenomenon of a rabbit pulled out of a hat). Mathcad quite suc-
cessfully solves problems of various nature (including economic ones) without
turning to pure programming (BASIC, C, Pascal, etc.). Moreover, that time I was
working on one of my books on Mathcad. Example with paint could only give
more colors to the book (unintentional pun).
Protocol of paint “check weighing” in Mathcad’s environment is shown in
Fig. 7.9. Comments explain what is happening in the formulas. In the
Mathcad-worksheet a matrix is formed with the name Price, whose elements store
the paint cost values depending on the combination of packaging. Then, some
elements of the matrix Price are set to zero, if the packaging combinations are not
in the value. The rest—sleight of hand and no mathematics: use function match to
determine the line number (37) and column number (6) of the matrix Price, the
element with the maximum value is located at their intersection. Olga was
unpleasantly surprised with the answer (37 small and 6 large cans). She unwittingly
deceived me for 175 thousand rubles.
So Mathcad saved me 175 thousand, though old, but rubles. Money is not so
great, but if we add to it the new computer study placed in the book, a new lecture
topic and a new laboratory work in computer science, the game was worth the
candle.
When I came back home to Moscow from Tambov, in a relaxed atmosphere at
my home computer I one more time analyzed the problem. And that are the results.
114 V. Ochkov
Firstly, one could make the process of finding solutions in MS Excel to correctly
“explain” the paint problem through changing the initial settings. And for this it was
necessary to not be lazy and click the Settings key in the Solver dialog box. It was
enough in the new dialog box Solver Parameters to reduce the permissible
deviation from 1 to 0.01 %. After that, the right solution would be found. Frankly,
in MS Excel it is not a solver procedure that is bad but its initial settings. Very few
MS Excel users resorting to its services, click the Parameters key. Normally, those
who understand the essence of optimization do not work with MS Excel. Hence the
confusion. Table 7.1 provides a complete analysis of the optimal solutions.
Secondly, when I showed this table in MPEI procurement department, they told me
that the best option for me (for me important is money) and for MPEI (for them
important is paint) will be the forth option: The “Tambovenergo” would spend
almost all the money, and the 885 l of paint, as it may seem strange, is more than
910, and 915. The matter is that with the large packaging a lot of paint is lost due to
7 Thermal Engineers’ Gift to Water Chemists 115
Fig. 7.10 The calculation for the group of filters used for treatment of make-up water of a steam
turbine unit
overflows into a smaller container. One could take a 15-l jar to a room under repair
and make full use of it.
The wrong solution of the problem is caused not only by the bad practices or
defective software, but also by the fact that the user himself does not really know
what he wants. All programs for solving the linear programming problems require a
clear formulation of a single objective function. In solving the educational problems
the goal is clear. What is the objective in real life? However this is not mathematics,
but philosophy…
Computers in the optimization of thermal power facilities allow wide use of
brute force method, which was applied in the paint problem. For example, when
designing a make-up water treatment plant for steam turbine units, one can take
many filters of small diameter or a few large diameter filters. The best option will be
found through an exhaustive search [51].
Figure 7.10 shows a page of the website where you can set the water flow
through the filters, high-speed characteristics (range of possible changes in the
filtration rate) of one filter, and diameters of commercial standard filters. After
clicking on Recalculate button, all the options will be searched through by the
number of installed filters and their diameters to choose an optimum option whose
real rate of filtration is as close to a specified recommended rate as possible.
116 V. Ochkov
Abstract This chapter shows how to solve thermal engineering problem and to
involve for this the audience of the specialized forums in Internet. Discussed the
problem of the density of water depending on its parameters.
Water has the unique property—its density has a maximum value at a temperature
slightly above zero Celsius. Therefore, and due to the fact that solid phase of water
(ice) has a density lower than water, in winter basins generally do not freeze to the
bottom—the temperature at the bottom of a basin has a value of about +4 °C.
The properties of water, as a practically incompressible fluid, are largely
dependent on the temperature and little depends on the pressure. Nevertheless, we
will try to develop the maximum density water temperature/pressure function.
Figure 8.1 shows a corresponding calculation made in Mathcad 15 environment.
After reference to file H2O.xmcdz, one can see in the calculation the function
wspDPT, returning water and steam density D as function of pressure P and
temperature T. Then using this the “cloud” function create another function tmax(p)
with one argument p. To find the functional relationship maximum, Mathcad’s
environment provides for built-in function Maximize, but it did not work in our
problem (see the error message shown in Fig. 8.1). Probably because of the fact that
our function wspDPT in function Maximize is called without listing the arguments
in brackets.1 We had to use another Mathcad’s built-in function—root, which
1
It was the author’s assumption, which turned out to be false (see below).
Fig. 8.1 Diagram of maximum density water temperature as a function of its pressure
returns the values of its second argument t, with which the first argument (partial
derivative function wspDPT for t at fixed p) is equal to zero. The third (1 °C) and
the forth (7 °C) arguments of the function root—are the boundaries within which
the required temperature values are to be found.
Figure 8.1 shows the four parametric curves: three water density isobars as a
function of temperature at p = 0.1, 0.5 and 10 MPa (thin line, dotted line and dots)
and a solid line linking the points where the density is at its maximum. The diagram
shows that the peak shifts to the left with an increase in water pressure.
However, the author is not satisfied with the solution shown in Fig. 8.1. First of
all because creating the diagram took a very long time and because the diagram was
not created at the most interesting point—at a pressure close to atmospheric.
Before and after the diagram in Fig. 8.1 there are operators available which
include a Mathcad’s built-in function time. It has a formal argument, whose value
does not affect the result of calling the function (here we used zero), and returns
some computer operating time in seconds. Of course this time is conveying not
much,2 but using the difference between the time values, you can calculate the
execution time of individual operators of a Mathcad-document. The diagram shown
2
The thermodynamic analogy: working medium’s enthalpy and entropy values don’t add-up too
much: It is important to know the difference of these parameters at various points in the ther-
modynamic cycle.
8 Maximum Water Density 119
in Fig. 8.1 took a long time—1 to 3 min, depending on the type and load of the
computer on which these calculations are carried out. And this is understandable:
the numerical definition of the derivative—is a very complex and not entirely
correct operation in terms of “pure” mathematics. Because of this fact (the duration
of calculation), such “slow” diagrams at the time of program debugging have to be
disabled by the appropriate command which appears in the list of menu commands
called by right-clicking, and then again included in the calculation.
The second thing that the author did not like in the solution shown in Fig. 8.1—
is that the line of water maximum density as function of pressure was “torn” near
the atmospheric pressure point: function tmax cannot calculate the temperature of
maximum density of water at p = 0.1 MPa. Therefore, an attempt was made to
return to the function Maximize.
All computer programs are supported by the developers, including through the
users’ formula. If the user has some question, he can post it in the appropriate forum
and look forward to hearing other forum visitors’ responses. The Mathcad also has
its forum and not one. The main Mathcad forum (PTC Community) is located at
https://www.ptcusercommunity.com/community/mathcad. One could also come to
this forum directly from the Mathcad’s environment by choosing the appropriate
command from the Help menu (Mathcad 15, in Mathcad Prime environment this
command can be found on tab Getting started). In this forum without registration
one can only view open resources. But after the registration it will be possible to ask
questions and respond to other people.
In Mathcad forum the author asked a question whether it is possible for the
problem shown in Fig. 8.1 to use not “capricious” and slow function root with a
derivative as an argument, but function Maximize. One man from Austria Werner
Exinger responded to this “call for help”. He showed that it is enough just to swap
around pressure and temperature in the function, and it will work (Maximize
function change only last argument of input function). Our dialog recorded in
Fig. 8.2 [52].
Figure 8.3 shows the corrected calculation. It began to work much faster, despite
the fact that the number of points on the diagram became not 30 (as in Fig. 8.1), but
300. And most importantly, the diagram has ceased to be torn—function tmin began
to return a response throughout the specified range of pressure—from 0.01 to
10 MPa.
Conclusion: If, when you are working in Mathcad, difficulties arise, you should
only ask a question on a forum (see Fig. 8.2). A positive response is likely to be
obtained. Moreover, in this forum, one can create a subforum (Group), open to all,
or only to invited users. In this “cyberspace” we can work together to solve
problems, including heating engineering problems. The author have opened on a
Mathcad forum his own Group for this book, where it is possible to discuss the
book, add new or adjust old examples. It is an English-language book forum.
Russian forum, as noted in the Preface, is located at www.thermal.ru. On these
forums one could download “cloud” patterns and functions necessary for
calculations.
120 V. Ochkov and K. Orlov
Fig. 8.3 Adjusted calculation of maximum density water temperature/its pressure function
On the site of the reference book [23, 38] one can find the water/steam properties
“cloud” functions not only for Mathcad program, but also for its Russian clone—
SMath (www.smath.info). Figure 8.4 shows a page of the reference book with water
property functions (the first area of formulation IAPWS-IF97).
If you click on the link SMath in the water specific volume row after you have
installed on your computer from the site www.smath.info the SMath program itself,
a hidden area would appear in its environment with the relevant function wspV1PT
(Fig. 8.5). An inverse function named wspD1PT (water density as function of
pressure and temperature) has been used in Chap. 4 here in the description of the
modern Internet technology for publishing calculations.
Using function wspV1PT one can create function wspD1PT (water density) and
analyze it near the temperature of maximum density of water.
The problem of the maximum density of water at a fixed pressure is easy to solve
in MS Excel as well, if in addition to the WaterSteamPro attach to these spread-
sheets another add-in—Solver (Fig. 8.6).
Figure 8.7 shows MS Excel table, where in cell B2 is entered a first approxi-
mation value—the water temperature on the Celsius scale. In cell B3, this tem-
perature is shown in Kelvin (=B2 + 273.15). Cell B4 has water pressure value
entered in megapascals, which in cell B5 is converted into Pascals (=B4 * 10 ^ 6).
Cell B6 stores WaterSteamPro function wspDPT (B5; B3), that returns the density
of water as function of pressure and temperature. Figure 8.7 also shows Solver
Parameters dialog box, which called from the Data tab, and stating that the
computer has to automatically change the value in cell B2, so that the value in cell
B6 took its maximum.
122 V. Ochkov and K. Orlov
After pressing the Solver key in the Solver Parameters box (this key is not
shown in Fig. 8.7, but one could see it in Fig. 7.8 in Chap. 7) in cell B6 appears the
value 999.9748 (maximum density of water in kilograms per cubic meter at pres-
sure 0.1 MPa), and in cell B2—3.96335 (maximum water density Celsius
temperature).
8 Maximum Water Density 123
Fig. 8.8 Water density as a function of pressure and temperature: one level lines
Abstract This chapter is about calculation of steam turbine cycle of nuclear power
plant operating on wet steam. You can learn about authors Internet sites that help to
solve thermal engineering tasks.
Figure 9.1 shows a principle scheme of a steam turbine unit of a double circuit
power reactor of a nuclear power plant (NPP) [53] (so called Pressure Water
Reactor—PWR). Pressurized water flows from the nuclear reactor (NPP primary
circuit) to the steam generator (SG) and transfers heat to the boiling water of the
second circuit. Saturated steam is supplied from the steam generator to the high
pressure cylinder (HPC) of the steam turbine. Part of the live (fresh) steam from the
steam generator enters the superheater. Before being superheated, the exhaust steam
from the HPC is partially dried in the separator where some water is separated from
the wet steam, and then returned to the circuit with bypass of the low pressure
cylinder (LPC) of the steam turbine.
This complication of the thermal circuit of a nuclear steam turbine cycle com-
pared with a thermal power plant on fossil fuel (see, e.g., Fig. 4.2 in Chap. 4) is
associated with the following NPP feature: if saturated steam is supplied to a
turbine, whose outlet pressure is 5–4 kPa, moisture in the last stages of the steam
Fig. 9.1 Calculation of the thermal efficiency of the nuclear power unit with wet steam: Input data
and the unit scheme
turbine will be unacceptably high which may lead to the steam turbine failure.
Steam moisture at the outlet of a nuclear turbine without separator and superheater
can be estimated by referring to the two Internet sites [23], as shown in Figs. 9.2
and 9.3.
First, steam pressure upstream turbine (6.2 MPa—see input data in Fig. 9.1) can
be used to determine steam temperature, taking into account that this is saturated
water/steam.
Then the site, shown in Fig. 9.3, allows you to calculate (in online mode) the
turbine steam expansion process by setting the values of initial pressure, initial
temperature, final pressure and relative internal efficiency of the process.
Note. The view of the site (URL: http://twt.mpei.ru/MCS/Worksheets/PTU/h-s-
ExpEng.xmcd) slightly differs from what is shown in Fig. 9.3. The point is that the
original site was copied to the clipboard by pressing <Print Screen> key, moved to
the graphics editor and streamlined.
The calculation, shown in Fig. 9.3, demonstrates that the wet water steam, when
expanding in the turbine from initial pressure of 6.2 MPa, at the condenser pressure
of 4 kPa will consist of nearly one third water, which, as noted above, is unac-
ceptable for the last stages of the turbine, since they cannot operate with such steam
moisture.
9 Nuclear Power Plant Steam Turbine Cycle 127
Fig. 9.2 Internet site to determine the water boiling point as function of pressure
Because of this feature of the steam turbine plant, the nuclear steam turbine is
divided into two cylinders—high pressure cylinder (HPC) and low pressure
cylinder (LPC). Steam, before it enters the LP cylinder, is dried in the separator and
heated in steam-steam heater (superheater) using heat of live steam supplied from
the steam generator (see the diagram in Fig. 9.1).
Figure 9.4 shows the site we have already referred to for calculation of steam
expansion process in the turbine with another (not 4 kPa, but 0.5 MPa) pressure at
the turbine outlet.
Comparing Figs. 9.2 and 9.3, you can see that with increasing pressure at HP
cylinder outlet from 4 kPa to 0.5 MPa, moisture of steam leaving the turbine has
dropped from 27.8 to 16.8 %. This percentage of steam moisture in the last stages of
HP cylinder (16.8 %) is quite acceptable for trouble-free operation.
You can “by hand” calculate thermal efficiency of the steam turbine cycle shown
schematically in Fig. 9.1, referring to the web-sites supporting the Reference Book
[23] and shown in Figs. 9.2, 9.3 and 9.4. Thus, Fig. 9.5 represents an h, s-diagram
of steam expansion process in LP cylinder after separating and superheating the
steam. But we will proceed with our calculations on a computer with Mathcad, the
beginning of computations is shown in Fig. 9.1.
As we have noted in Chap. 2, there are water and steam property functions
Water-SteamPro (www.wsp.ru) that can be attached through DLL to Mathcad 15
(see Fig. 2.21) and Mathcad Prime 3 (see Fig. 2.23). In Mathcad Prime 3 (namely in
this software we will calculate the NPP cycle, this was already started in Fig. 9.1),
the work process has slightly changed compared to Mathcad 15. In Mathcad 15, the
function arguments with suffix wsp could not be automatically resized for given
dimension, so the arguments had to be left with no dimensions (see Fig. 2.24 in
Chap. 2). In Mathcad Prime 3, this operation (de-dimensioning of the arguments) is
128 V. Ochkov and K. Orlov
Fig. 9.3 Internet site for calculating the turbine steam expansion process
done automatically, but the WaterSteamPro functions continue (as in Mathcad 15)
to give dimensionless quantities corresponding to the base SI units: pressure in
pascals, and not in megapascals, temperature in degree Kelvin, and not by centi-
grade, etc.
To function properly in Mathcad Prime 3 with dimensional values, the
WaterSteamPro functions should be overwritten as shown in Fig. 9.6. As in Chap. 2
we highly recommend to reference file watersteampro.xmcd for Mathcad 15 or
redefine functions like it showed on Fig. 9.6. In this case all functions arguments
and results will be automatically translated to desired dimensions of
WaterSteamPro functions inputs and back with results.
In a Mathcad document the function overwriting operators are grouped into a
separate area that can be minimized by clicking on the icon “minus”, shown in the
upper left part of Fig. 9.6. Figure 9.10 shows this area minimized. In addition, this
area contains user’s entered “one-line” (without division operators) units of specific
9 Nuclear Power Plant Steam Turbine Cycle 129
Fig. 9.4 Internet diagram of steam expansion process in HP cylinder of a NPP steam turbine
entropy (kJ/(kg_K) and J/(kg_K)) and specific enthalpy (kJ/kg and J/kg).
Additionally, you can enter control of dimension of WaterSteamPro functions as
shown in Fig. 9.7.
You can clearly see in Fig. 9.7, that call of the function XEXPANSIONPTPEFF
with repositioned first two arguments (the original pressure and temperature) ends
with an error message. (Function XEXPANSIONPTPEFF was not on the list, as
shown in Fig. 9.6, since this function returns dimensionless quantity—steam
moisture at the end of expansion process in the turbine, and it is not necessary to
overwrite the function into dimensional form.)
After the WaterSteamPro functions became visible (through DLL) and dimen-
sional (through overwriting functions), it is easy in Mathcad Prime 3 to calculate
thermal efficiency of a nuclear steam turbine plant (see Fig. 9.7), that has been
started in Fig. 9.1.
130 V. Ochkov and K. Orlov
Fig. 9.5 Internet diagram of steam expansion process in LP cylinder of a NPP steam turbine
The calculation shown in Fig. 9.8 does not require special comments, and that is
a very useful consumer quality of Mathcad, which we mentioned in the introduction
and to which we will come back in Chap. 20.
9 Nuclear Power Plant Steam Turbine Cycle 131
Fig. 9.6 Overwriting of WaterSteamPro functions into dimensional functions in Mathcad Prime 3
Fig. 9.8 Calculation of thermal efficiency of steam turbine cycle of PWR NPP
Figs. 9.1 and 9.8, and checking if the formulas are entered correctly and correctness
of references to the table cells. But without these tips it would be quite difficult to
create even such a simple calculation table in MS Excel.
The author created a WaterSteamPro, which is often referred to in this book, in
the C++ programming language. But before that, all of the individual functions of
the package were written and debugged in Mathcad. It’s allows to save time for
development. Without this preliminary work, it would have been very difficult to
write a program in C++. In conclusion: When getting to the creation of a larger
project requiring the use of a programming language or an MS Excel spreadsheet, it
9 Nuclear Power Plant Steam Turbine Cycle 133
is useful to write and debug individual modules in Mathcad. It is much easier and
expedites the work by eliminating many errors. One can, of course, execute whole
project in Mathcad, but programming languages allow you to generate executable
code that runs on a computer without the programming languages themselves. As
for MS Excel (which, as a matter of fact, is installed on almost every computer),
these tables are used by many end users who want that a ready program should be
executed exactly in this software. The author of this book once developed a pro-
gram to optimize the water chemistry of a condenser cooling system of a thermal
power plant [3]. The program itself, minimizing the cost of prevention of scale
deposits in the cooling system, was written in Mathcad. Then it was manually
transcribed for MS Excel, because the staff of the power plant used to work with
these spreadsheets and did not want to learn any other programs.
Chapter 10
Isobar, Isotherm, Isochor…
Abstract This chapter describes how you can build contour lines (isobars, iso-
therms and isochores) and lines of phase transitions on the state diagrams of water
and steam.
Different types of diagrams could illustrate the thermal processes: T-, s-, h-, s-, p-, v-
and others. There are many of these diagrams in this book. Points on such diagrams
are often fixed as crossing points of isolines, e.g. isobars and isotherms (see, e.g.
Fig. 9.3 in Chap. 9). “The third time is a charm”—an isochor is missing here. Let us
use this small chapter to see how we can build this “third time” for isolines on
h, s-diagram [34].
Problem Superheated steam pressure and temperature values are given. It is nec-
essary to build in Mathcad’s environment on h, s-diagram, firstly, the water and
steam saturation line, secondly, the crosshair of isobar, isotherm and isochor
showing the specified pressure and temperature values. Figure 10.1 shows the
solution of this problem in Mathcad Prime 3.0.
The diagram in Fig. 10.1 built in Mathcad Prime 3.0 with WaterSteamPro
package connected. Without this package, you can build such a diagram in Mathcad
15 environment, making reference to the “cloud” file H2O.xmcdz, stored at
http://twt.mpei.ru/tthb, we have done this more than once, when solving some
problems in this book.
Figure 10.1 shows a collapsed Mathcad-area, which is further shown expanded
in Fig. 10.2.
Operators shown in Fig. 10.2 used for the following.
Fig. 10.2 Auxiliary operators in the problem of building the water/steam isobar, isotherm and
isochor
10 Isobar, Isotherm, Isochor… 137
To build an isochor, one needs to have functions that return values of specific
enthalpy and entropy as a function of density and temperature or pressure.
However, there are no such functions on the functions list in WaterSteamPro. But
one could create them, relying on the built-in Mathcad’s function root. This
operation described in Chap. 1 (see Fig. 1.10). This auxiliary problem could be
solved by using the third group of operators, as shown in Fig. 10.2. First, create a
function wspPDT, returning water/steam pressure as a function of density and
temperature, and second, create two functions—wspHDT and wspSDT—using
wsp-functions wspHPT and wspSPT, where the first argument will be not pressure,
but the previously created function wspPDT.
This chapter is small in size but very important to learn how to use the tool
functions when solving problems on a computer. Knowing this tool, one could
easily build other types of diagrams and other isolines on them, e.g. lines where
h = const (line of constant enthalpy, an isenthalpa) or s = const (line of constant
entropy, an adiabat). It would be very useful to learn how to use the function tool
from the point of view of making calculations without schedules. An expert (not
only thermal engineer), creating a set of his own user-defined functions, linking
different values, could easily solve complex problems by combining and inserting
into each other various previously created and debugged functions.
Chapter 11
Construction of Forward and Backward
Functions on Properties of Working Fluids
on Tabular Data Base
Volodymyr Voloshchuk
Abstract This chapter shows how to create functions that return the values of the
thermophysical properties of working fluids of power engineering, on the basis of
tabular values in different areas of the state: single-phase region, two phase region,
the saturation line.
V. Voloshchuk (&)
National University of Water Management and Natural Resources, Rivne, Ukraine
e-mail: [email protected]
Fig. 11.1 Fragment of the program for creating a function to calculate the specific entropy of the
refrigerant R410a in single-phase regions of subcooled liquid with two arguments—pressure (first
row) and temperature (first column)
values (in this case specific entropy) which are functions of the determined argu-
ments. The function arguments, pressure p and temperature T, are located in the
first row (“hat”) and first column (outset) respectively as published on paper media.
If values of pressure and temperature overflow beyond the determined boundaries
the values are not calculated. Thus, the user clearly sees the initial data applied for
the calculation. In Fig. 11.1 the input matrix corresponds to a state of subcooled
liquid. Below the critical pressure and temperature such matrixes are only partially
filled. But to realize interpolation tool built in Mathcad it is required the matrix to be
filled completely. That’s why initial matrix should be added with fictitious values.
Extrapolation can be used for this purpose. In this case it is also possible to
calculate values beyond the scope of the input matrix and to use iterative methods.
In Fig. 11.1 the region of fictitious values is separated. In the same manner the table
is created for the superheated vapor region of the refrigerant (see Fig. 11.2). For
convenience matrixes on properties of the working substance in single-phase
regions without phase transition (at supercritical regions) can be inserted addi-
tionally. In this case the matrix will be filled completely and no fictitious regions
will be required. As a result it is possible to insert several matrixes corresponding
different argument ranges.
The next step of the program is choosing the matrix which includes the specified
values of the arguments. As shown in Fig. 11.3 there are three different matrixes for
three ranges of arguments. The first matrix M1 corresponds to the single phase
region of subcooled liquid below critical pressure—so, for this matrix values of
pressure which fit this region should be applied. Moreover the matrix M1 can be
used for temperature lower than the saturation temperature at the given pressure.
For this purpose a special function for calculating the saturation temperature
depending on pressure (R410aTSLP(p)) is applied. Here the first five letters
indicate the type of working substance (R410a). The sixth letter indicates the
parameter, which is calculated (temperature T). The next two letters indicate that
11 Construction of Forward and Backward Functions on Properties … 143
Fig. 11.2 Fragment of the matrix on values of entropy of the refrigerant R410a in a superheated
vapor region as a function of two arguments—pressure (the first row) and temperature (the first
column)
the function calculates values of the saturated liquid (SL). The last letter indicates
argument on which the saturation temperature is calculated (pressure p). This
parameter (argument of the function) is duplicated in the brackets after the name of
the function. The smallest permitted value of temperature (250 K) and the range of
permitted pressure have been specified above (see Fig. 11.1). It should be noted that
refrigerant R410a is a zeotropic mixture. It is a non-ozone depleting blend of two
HFC refrigerants (R32, R125) which have different saturation temperatures at the
same pressure. It means that evaporation starts at the saturated liquid temperature of
the most volatile component and the temperature progressively rises during evap-
oration until the saturated vapor temperature of the less volatile component is
reached. That’s why saturated liquid and vapor temperatures of the R410a are
different. As a consequence functions on properties of this refrigerant at saturated
liquid and saturated vapor lines will be different too. For example, for calculation
of saturated vapor temperature for a given pressure a function R410aTSVP was
created which differs from function R410aTSLP in denoting the eighth letter by V
instead of L. In this case collocation SV points out that saturated vapor parameter is
calculated. The matrix M2 corresponds to the single phase region of superheated
vapor below critical pressure and temperature ranging from 250 to 340 K. In the
given case this matrix can be used for temperatures lower than 340 K and higher
than saturated vapor temperature for a given pressure. The left region of the matrix
can be expanded to the critical point. But in the critical region the values of the
function are rapidly changing and it is recommended to insert additional matrix
with smaller step. The matrix M3 corresponds to the single phase region of
superheated vapor below critical pressure and temperature ranging from 340 to
450 K. There is no saturation line in this matrix. Matrixes M2 and M3 border each
other along the isotherm T = 340 K. The border should be included into one of
these matrixes. If the input value of the temperature T is not included in the scope
of permissible values it can be located in the two-phase region or on the saturated
line. In this case the program returns an error message “The point is in the double
phase or saturated line” (see Fig. 11.3) using built in Mathcad operator error(S),
where S—the text of the error message.
144 V. Voloshchuk
Fig. 11.3 Fragment of the program for creating a function to calculate the specific entropy of the
refrigerant R410a in the single-phase region
The next fragment of the program for creating a function on the specific entropy
of the refrigerant R410a in the single-phase region is shown in Fig. 11.4. With the
help of built in Mathcad operator submatrix additional matrixes out of the chosen
matrix M are extracted—the first row (vector of pressures p′) and the first column
(vector of temperatures T′). The vector of pressures is transposed into a column
with the special operator built in Mathcad. Then the operator submatrix forms
matrix s′—values of specific entropy of working fluid, which is also transposed
swapping rows and columns. In Mathcad the built-in variable ORIGIN is set to 0.
It means that the first row/column is designated with the index 0.
In the program shown in Fig. 11.4 a method with double spline interpolation is
used (see also Fig. 1.1). At first out of the columns of the entropy matrix s′ an
additional line (vector s″) for a given pressure is generated by spline interpolation.
Values of vector s″ correspond to the values of temperatures from the vector T′.
Further again spline interpolation is used for calculating the specific entropy as a
function of the temperature absent in the vector of temperatures T′. The Mathcad has
its own operator of spline interpolation with the help of which it is possible to create a
function of two arguments. But in this case there is one important limitation factor—
the matrix should be square shaped. The proposed method of double spline inter-
polation can be used for rectangular matrix too. The last operator of the programming
block shown in Fig. 11.4 returns the necessary unit of the entropy—kJ/(kg K).
Fig. 11.4 Fragment of the program for creating a function on the specific entropy of the
refrigerant R410a in the single-phase region
11 Construction of Forward and Backward Functions on Properties … 145
Fig. 11.5 Graphical illustration of the created function on specific entropy of the refrigerant
R410a in sp- and sT-phase diagrams
146 V. Voloshchuk
Fig. 11.6 Fragment of the program for creating the backward function for calculating the
temperature of the refrigerant R410a in the single-phase region depending on pressure and specific
entropy
11 Construction of Forward and Backward Functions on Properties … 147
In the given case there are three matrixes regions which are graphically shown in
Fig. 11.5. For convenience matrixes M2, M3 are merged because they belong to one
region (see Fig. 11.5).
The next part of the programming block determines the matrix which contains the
input arguments. If this is a single-phase region the specific entropy which is an
argument should be less than specific entropy of saturated liquid at the specified
pressure. For calculation of specific entropy on saturation line two additionally
created functions are used. At first saturated liquid temperature is calculated with the
function R410aTSLP(p), and then this value of the temperature is used to calculate
saturated liquid specific entropy with the additional function R410aSSLT(T).
The singe phase region of the subcooled liquid of the refrigerant is specified by
matrix M1 (see Fig. 11.5). Using the built-in operator submatrix similarly as in the
case of creation of the function R410aSPT(p, T) three matrixes are generated: a
matrix of pressures p′ which is transposed into a column, a matrix of temperatures
T′ and a matrix of specific entropy which is also transposed swapping rows and
columns (see Fig. 11.7).
Next, out of the columns of the specific entropy s′ an additional line (vector s″)
for the given pressure is generated by spline interpolation. Values of vector s″
correspond to the values of temperatures from the vector T′.
In this case, unlike the previous example (see Fig. 11.4), the vector of the
specific entropy s″ is the arguments vector to perform the second spline interpo-
lation. The value of the searched temperature is generated by interpolating the
values of the matrixes s″, T′ for the given value of specific entropy s.
One control is applied in Fig. 11.7. The value of the second argument, the
specific entropy s, for the given first argument—pressure p, should remain within
the values of generated vector s″. It should be higher than the first value of the
Fig. 11.7 Fragment of the program for creating an in backward function for calculating the
temperature of the refrigerant R410a in the single-phase region of subcooled liquid
148 V. Voloshchuk
Fig. 11.8 Fragment of the program for creating backward function for calculating the temperature
of the refrigerant R410a in the single phase region of superheated vapor
Fig. 11.9 Fragment of the program for creating the backward function for calculating the
temperature of the refrigerant R410a in the double-phase region
matrix s″. If this control is not performed the created function will return a value of
temperature by extrapolation, not interpolation.
The last operator in the programming block shown in Fig. 11.7 returns the
desired unit of the temperature—K.
The similar fragment for creation of the backward function for calculation of
temperature of the refrigerant R410a in the single phase region of superheated
vapor is demonstrated (Fig. 11.8).
In the proposed backward function R410aTPS(p, s) calculation of the temper-
ature in the double-phase region is also provided. Linear interpolation is used for
this purpose (see Fig. 11.9).
A copy of the web-page for online calculation and graphical illustration of the
created backward function R410aTPS(p, s) for computing the temperature of the
refrigerant R410a depending on pressure and specific entropy is shown in
Fig. 11.10.
Figure 11.11 shows an example of using created functions on thermodynamic
properties of refrigerant R410a for calculating in Mathcad documents the amount of
heat q1_2, which is transferred to the working fluid in the isobaric process 1–2. In
this case a technology of “cloud” reference, that gives a possibility to refer to
“cloud” functions R410aHPT(p, T) and R410aTPS(p, s), which are located on the
internet-server http://twt.mpei.ac.ru/TTHB/2/R410aEng.html (see Fig. 11.12), is
proposed.
11 Construction of Forward and Backward Functions on Properties … 149
Fig. 11.10 A copy of the web-page for online calculation and graphical illustration of the created
backward function R410aTPS(p, s)
R s2 For computing the heat energy q1_2 a well-known formula can be used:
s1 Tðp; sÞds, where T(p, s)—is the backward function R410aTPS(p, s) for com-
puting temperature of the refrigerant R410a depending on pressure and specific
entropy. The proposed function can be directly used in integration and differenti-
ation operations. It substantially simplifies solving tasks of computing processes in
thermal engineering.
In the presented example it is possible to use the first law of thermodynamics—
rate of heat transferred to the refrigerant in the isobaric process 1–2 is equal
enthalpy change of the refrigerant.
150 V. Voloshchuk
Fig. 11.11 Using of the created functions on thermodynamic properties of the refrigerant R410a
for computing heat transferred to the refrigerant in isobaric process 1–2 with the help of reference
to “clouds”
Fig. 11.12 The web-page of “cloud” functions on thermophysical properties of refrigerant R410a
on saturated lines and single-phase regions
Valery Ochkov
Abstract This chapter shows how to create animated clips illustrating various
thermal processes.
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
Fig. 12.1 Interactive open network calculation process of steam expansion in the turbine
Figure 12.2 shows the calculation process of steam expansion in the turbine, which
is open not via the Internet (see Fig. 12.1), but directly in Mathcad 15 environment. In
this calculation process, the operator intended for assigning pressure values in the end
point pk is replaced with a formular in which appears the system variable FRAME. If
this variable is equal to zero, the steam pressure at the end point is equal to the
condenser pressure (100 kPa). If FRAME variable takes the value 999, the steam
pressure is equal to 24.3 MPa. Smooth change of the value of FRAME variable from
999 to 0 (in Fig. 12.2 it is equal to 99) makes it possible to create an animation and see
how the steam is expanded in a turbine. The FRAME variable can be changed not
only by Mathcad users—it can be automatically changed Mathcad itself. To do this it
is necessary and sufficient to call a dialog box Record Animation of the eponymous
command in Tools menu—see in the bottom of Fig. 12.2.
Before that, one should delete or disable the operator of manual setting of the
variable FRAME := 99, which was used to manually debug the program—form
individual animation frames. In the dialog box Record Animation specify the
12 Animation of Thermal Processes 153
initial (999) and the end (0) values of the variable FRAME (we count the frames in
a reverse order from 999 to 0) and the animation display rate (24 frames/s). In
addition, you need to broach the computer mouse to select an animation area. It is
limited by a dotted rectangle, as shown in Fig. 12.2. After that you should click
Animate. Now wait a little bit (2–3 min) until Mathcad itself starts to change the
value of FRAME variable, memorize the animation files and “put” them in the avi.
“container”. At the end of this process, the animation window will be displayed,
two frames of which are shown in Fig. 12.3.
Animations created in this way could be saved in an avi. file (see the key Save
As … in the dialog box in Fig. 12.2) and run without Mathcad. Animation could
also be placed in various forums, such as Mathcad’s user forum—PTC Community
(Fig. 12.4).
The author also collects other animations of thermal processes on the forum PTC
Community. Thus, in Fig. 12.5, you can see an animation of the processes
occurring in the steam turbine cycle with steam boiler (economizer, a drum with
downcomers and risers), steam turbine, condenser and feed pump.
The lower right corner of Fig. 12.5 includes references (More Like This) to the
animation of the other thermal cycles: Otto cycle (Fig. 12.6) and Diesel cycle
(Fig. 12.7). When creating this animation, a Mathcad opportunity was used to
change the image (image of phase of the internal combustion engine cycle) in the
calculation.
Animation in Mathcad environment can be used not only for visualization of
various processes, but also for simple calculations. One disadvantage of Mathcad is
that it cannot generate so-called executable files, for example, exe. files. Because of
this, calculations created in Mathcad can only be run when it is installed on the
computer. One solution to this problem is the Mathcad Calculation Server, which
we have repeatedly mentioned. Another partial solution to this problem may be …
animation. If our calculation depends only on one variable, we can link this variable
to the system variable FRAME, create an appropriate animation, and then run it
outside Mathcad using some video player, move the animation slider, set the value
of the argument and read the value of the function. Figure 12.8 shows a frame of
animation linking the temperature to the saturation pressure of water and steam.
One can open the animation using the slider shown in the bottom of the figure, set
the temperature with increment of 0.5 °C (K) and calculate the saturation pressure.
The other pairs of steam and water thermal properties at saturation line can be
added into this animation, such as water and steam density, specific enthalpy and
entropy, thermal conductivity, viscosity, etc.
Another aspect of the use of animation: Fig. 12.9 shows the frame of the animation
of R407c isotherm motion in coordinates “specific entropy—pressure” at a fixed
pressure and varying temperature. Considering the behavior of the moving curve one
can to some extent judge on the quality of the function created using the principle
“Nature does not like sharp corners and… oscillation” (see Fig. 1.8 in Chap. 1).
Figures 12.10 and 12.11 shows a pair of frames of the animation for scanning
isolines using one argument of a function that returns the value of the coefficient of
friction f of the fluid in a circular tube, depending on the Reynolds number Re and
154 V. Ochkov
the liquid flow and the relative roughness of the tube inner surface Δ. This problem
will be discussed in detail in Chaps. 16 (see Fig. 16.7) and 17. The function f (Re,
Δ) was created using the double interpolation method described in Chap. 1 (see
Figs. 1.1 and 1.2). So, if you use two spline interpolations, the function f (Re, Δ) is
obtained with oscillation (see Figs. 12.10 and 1.8). If a spline interpolation is
replaced with a linear interpolation, the function f (Re, Δ) gets pretty smooth (see
Fig. 12.11). Without animation this error couldn’t have been noticed.
Divertisment English blacksmith Thomas Newcomen, in 1705, built a steam
engine, which was used to pump water from mines—see the author’s
Mathcad-animation of this machine at https://www.ptcusercommunity.com/videos/
2134. A legend is connected with this machine. A boy Humphrey Potter in 1712 got
the job of operator of such a steam-driven pump. To implement the working cycle
of the machine, the boy had to open the steam supply valve when the piston was in
its lower position, and the water supply valve—at fully lifted piston. But this genius
boy figured out how to rid yourself of the repetitive work. Using sticks and rope he
coupled the handles of the valves with the stem of the steam engine.
12 Animation of Thermal Processes 157
Fig. 12.9 Scanning with animation of a function created using the properties of substances
Fig. 12.12 Frame of the animation of the steam engine with oscillating cylinder
It turned out some sort of rocker mechanism. The valves were controlled by a
piston itself, and the boy was able to play with his friends. Thus was invented the
slide valve of steam engine—an automatic steam distribution device.
Now one could only see a steam engine hardware in a museum. On the Internet,
one could find a lot of animations of steam engine alone or mounted on a steam
train or boat. It is quite realistic to create an animation of the steam engine using a
mathematical program Mathcad [57], and many heat engineering students take
interest in this.
Figure 12.12 shows the animation frame of a rather unusual steam engine with
oscillating cylinder. The animation itself is placed on the Planet PTC website
https://www.ptcusercommunity.com/videos/3033. There one could download a
Mathcad-calculation.
Chapter 13
Calculation of Gas Turbine Engine Cycle
Volodymyr Voloshchuk
Abstract This chapter describes in detail the process of calculating energy char-
acteristics of the gas turbine cycle (Brayton cycle) with modeling of the compo-
sition and thermophysical properties of air, natural gas and products of combustion
along the flow in gas turbine.
For actual thermal power plants the most comprehensive measure of their perfor-
mance is the overall efficiency that considers all losses.
It can be shown that the thermodynamic cycle, which in the given conditions has
the largest internal absolute efficiency of the cycle gi , provides the highest overall
efficiency of the plant. It gives a possibility for actual, internally irreversible, cycles
of thermal power plants that produce electricity to accept as the main measure of
their thermodynamic efficiency a parameter which is expressed by following rela-
tionship [81, 82]
lgte act
gi ¼ ; ð13:1Þ
q1 act
where lgte act —net specific work output of the actual gas turbine engine cycle;
q1 act —actual specific heat addition from outside to the cycle.
Let’s analyze influence of working fluid parameters (temperature and pressure)
on the internal absolute efficiency of the cycle of simple gas turbine engine with
continuous combustion of fuel at constant pressure sketch of which is shown in
Fig. 13.1.
V. Voloshchuk (&)
National University of Water Management and Natural Resources, Rivne, Ukraine
e-mail: [email protected]
Fig. 13.1 Sketch of the simple gas turbine engine with continuous combustion of fuel at constant
pressure: C—compressor; F_Comb—fuel combustor; GT—gas turbine; GEN—electric generator
Turbine efficiency:
Compressor efficiency:
lc h2 h1
gis ¼ ¼ :
c
lc act h2 act h1
Then
where
ð h2 h1 Þ
h2 act ¼ h1 þ :
gis c
In the cycle of gas turbine engine work input required for compression of air is
quite large and takes a significant share of the work output from the gas turbine.
Therefore, in analysis of the internal absolute cycle efficiency the ratio of these
values (/ ¼ lc =lgt ) plays an important role.
Using this ratio the net specific work output of the actual gas turbine engine
cycle can be expressed as
13 Calculation of Gas Turbine Engine Cycle 163
/
lgte act ¼ lgt gis :
gt
gis c
As a result it can be seen that gas turbine engine has network output
(lgte act [ 0) only if ðgis gt /=gis c Þ [ 0 or gis gt gis c [ /. It means that a
quality of both compressors and gas turbines which is characterized by their effi-
ciencies plays an important role for such engines.
For further analysis it may be expressed as
q1 q2 lgt lc
gt ¼ ¼ ;
q1 q1
Each of the terms in Eq. (13.5) has its own dependence on the parameters of the
cycle. Therefore, internal relative efficiency of the actual cycle gi gte is greatly
complicated by many conditions.
Calculations of different types of gas turbine engines are located on the web-site
http://twt.mpei.ru/ochkov/VPU_Book_New/mas/index.html. In the given study
approaches which were used in creation of programs for computing thermodynamic
cycles of gas turbine engines will be analyzed.
The calculated scheme of simple open gas turbine engine consisting on the
compressor, fuel combustor and gas turbine is shown in Fig. 13.1.
For determining thermodynamic parameters of working fluids WaterSteamPro
(www.wsp.ru) was used.
The calculation begins from interactive procedure of data input. The number of
entered parameters can be different. It depends on how detailed calculation should
be. An example of input data block, which is prepared on the technology of
Mathcad Calculation Server, is shown in Fig. 13.2. Such set of data enables to
conduct rather detailed thermodynamic calculation of the given gas turbine engine
cycle. In the lower part of the Fig. 13.2 operators of interaction with a user during
creation of Internet MathCad documents are presented.
For convenience acceptable units are chosen (for example, kilojoule per kilogram
of gas—kJ/(kg g)) (see Fig. 13.3). Additional procedures during calculation (for
example, procedure of unit selection, reference to the WaterSteamPro and so on)
164 V. Voloshchuk
Fig. 13.2 Input data block of the Internet document prepared on technology of Mathcad
Calculation Server
metallurgical limits on the gas turbine operating temperature go through the fuel
combustor and the gas turbine.
The compressor is the first element of the calculation which operates with moist
air. At first humidity should be determined. Relative humidity at the specified
temperature is the input data for this. Marginal pressure of water vapor of inlet air is
calculated. At temperatures lower the triple point (0.01 °C) this pressure is com-
puted as sublimation pressure, at temperatures higher 0.01 °C—as saturation
pressure (see Fig. 13.4). These values of pressure are calculated with the help of
WaterSteamPro functions—wspPSUBT(T) and wspPST(T) respectively.
Then mass and molecular humidity ratio are calculated on the base of thermo-
dynamic equations (see Fig. 13.4). The calculation mixture of the compressor inlet
Fig. 13.4 Creation of the calculation mixture of the working fluid in the compressor inlet
166 V. Voloshchuk
temperature can be derived from knowing humidity ratio of the air. For this purpose
a buil-in MathCad function concat(s1, s2,…) (returning the string formed by con-
catenating strings s1, s2, and so on) and num2str(x) (returning the number x to a
string) are used. The result string (in WaterSteamPro this is a gas specification) is
used for calculating thermodynamic properties of working fluid in the compressor.
In Fig. 13.4 with the help of a buil-in MathCad function wspgHGST(gs, T), where
gs—gas specification, and T—temperature, specific enthalpy of the moist inlet air is
calculated.
Functions in WaterSteamPro for calculating properties of ideals gases and ideal
mixtures have the prefix wspg, then calculated parameter (in the given case H—
enthalpy) and a list of arguments (gas specification GS and temperature T) are
indicated.
Both chemical denotation of single gases (CO2, H2, O2 and so on) and their
mixtures can be used as gas specification indication. Mixtures specifying is set by
rather simple method. Components of mixture are separated with semicolon and the
amount of the component in the mixture is pointed out after the component indi-
cation followed by colons. The number of moles or mass but not fractions should be
used in the amount specifying. In case of default the value means mole but if the
denotation “M” is used mass units are applied. So, the indication “AirMix:1; H2O:
0.0102062435658935” means that one mole of dry air “AirMix” and
0.0102062435658935 moles of water vapor are selected.
After creation of the mixture of the working fluid calculation of the compression
is performed. For this purpose, using WaterSteamPro functions, thermodynamic
parameters of the air at inlet and outlet of the compressor are computed. Isentropic
and actual works for the air compression are also determined (see Fig. 13.5).
For calculation of thermodynamic properties of moist air the following functions
of WaterSteamPro are used:
• wspgTGSPS(gs, p, s)—temperature as function of gas specification gs, pressure
p, specific entropy s;
• wspgHGST(gs, T)—specific enthalpy as function of gas specification gs,
pressure p, temperature T;
component gs2 in the mixture gs1. Oxygen O2 is a component and moist air is a
mixture. So, the amount of moist air required for burning 1 m3 of fuel is evaluated.
Then, using coefficients of stoichiometric combustion reaction theoretical volumes
of gases (components of combustion products) are calculated. The MathCad
functions concat(s1, s2,…) and num2str(x) create the combustion products mixture
(see Fig. 13.9).
The required temperature of the working fluid at gas turbine inlet should be
provided in the fuel combustor. Combustion with excess air, i.e. more air than is
required for stoichiometric reaction and calculated above (see Fig. 13.8), is applied
for this purpose.
13 Calculation of Gas Turbine Engine Cycle 169
Fig. 13.10 Comparing enthalpies of working fluids with temperature of determining lower
heating value of fuel
Knowing values of enthalpies of all flows entering and leaving the fuel com-
bustor it is possible to apply a heat balance in evaluating air excess for providing
the required temperature at gas turbine inlet. It should be taken into account that the
temperature of determining lower heating value of fuel can be different. That’s why
values of enthalpies of all flows entering and leaving the fuel combustor are to be
compared with this temperature (see Fig. 13.10).
The total heat balance of the fuel combustor consists of heat input with mass of
air, mass of fuel, energy release with combustion reaction and heat output with mass
of combustion products and excessive air. The heat balance equation of the fuel
combustor is presented in Fig. 13.11. Evaluation of the amount of excessive air
Fig. 13.11 Evaluation of the amount of excessive air and the contents of the turbine working fluid
170 V. Voloshchuk
made with the help of MathCad calculation block Given—Find is also shown in
Fig. 13.11. The MathCad functions concat(s1, s2,…) and num2str(x) are used for
creation of the turbine working fluid contents in the format of the WaterSteamPro
gas specification (see Fig. 13.11). The working fluid in the gas turbine is a mixture
of combustion products and excessive air for cooling. The components of the
mixture are placed in curly braces using gas specification denoting. Then after
colons the amounts of the components in kg are indicated. In the calculation shown
in Fig. 13.11 a mixture of 38.529 kg of moist excessive air and combustion pro-
ducts consisting of 1 kg of fuel and of 17.316 kg of moist air as the oxidizer is used
for creation of WaterSteamPro gas specification.
During calculation of the gas turbine inlet and outlet pressures should be
determined. The gas turbine inlet pressure is calculated using the value of the
compressor outlet pressure taking into account pressure drops in the path
“compressor-fuel combustor-gas turbine inlet”. The gas turbine outlet pressure is set
to be equal to the atmospheric pressure minus pressure drops in output path of the
engine. The next evaluation of parameters of the gas turbine is performed using
WaterSteamPro functions (see Fig. 13.12) similarly as for the compressor but
taking into account that the working fluid is expanded.
Knowing the amount of the air needed for stoichiometric burning of 1 kg of fuel,
the amount of the excessive air for 1 kg of fuel it is possible to calculate mass flow
rate of the fuel and the working fluid through the gas turbine (see Fig. 13.13).
The total amount of the heat transfer to the gas turbine engine from outside
consists of heat transferred with mass of inlet air, inlet fuel and from burning fuel.
The results of gas turbine engine calculation may be: the compressor power, the
fuel compressor power, the gas turbine power, the output power, the overall effi-
ciency (see Fig. 13.13).
For visualizing the gas turbine engine cycle it is necessary to prepare appropriate
data—vectors elements of which are parameters of working fluid in different points
of the Brayton cycle. These vectors are prepared with Mathcad-operators according
to the procedure shown in Fig. 13.14. The vectors should be created for all four
processes of the cycle (compression, heat addition, expansion and heat rejection).
A cycle is formed by a sequence of processes such that the final state of the final
process is the same as the initial state of the initial process. The stack function links
together vectors belonging to four processes. It enables to represent the cycle in
different types of property diagrams (Ts-diagram, hs-diagram and so on). The rate
of the vectors is set by a user. In the presented example there are 300 elements in
one initial vector numbered from 0 to 299. In the Mathcad environment it is
possible to create vectors in two ways: through element method (based on the
vectors with index i: pci ; pf ci and so on) or by means of the vectorize operator
(vector indicator over the entire expression) (see Fig. 13.14).
In the example shown in the Fig. 13.14 the following order of vector forming is
introduced. At first vectors of pressures values are calculated in the compressor, fuel
combustor, gas turbine and engine outlet (respectively pci ; pf ci ; pt ai ; pcooli ).
These values changes in a uniform manner from the beginning figure to the end one
within one element of the engine (for example, for the compressor the beginning
pressure is p1—inlet pressure and the final one p2—outlet pressure compressor).
172 V. Voloshchuk
Fig. 13.14 The procedure for creating vectors used in sketching the Brayton cycle
Change of the mass (it is about 2 %) and contents of the working fluid will not
be taken into account when plotting the Brayton cycle. The cycle will be sketched
for 1 kg of air. So, recalculation of the elements in the gas turbine engine where air
is not a working fluid (heat addition, expansion, heat rejection) should be done.
Besides it should be noted that Mathcad calculates entropy in J/(kg K). But if we
represent the entropy in the diagram in kJ/(kg K) we need to divide the values of y
axis by kJ/(kg K). In Mathcad temperature is calculated automatically in Kelvins. If
we want to represent the temperature in units of degrees Celsius we need to subtract
from the values of x axis 273.15 K (see Fig. 13.15).
Isobars can be plotted on the Ts-diagram which enables to evaluate visually
pressures in typical points of the Brayton cycle (compressor outlet, gas turbine inlet
13 Calculation of Gas Turbine Engine Cycle 173
and so on) and pressure drops on the components of the cycle (see Fig. 13.15). The
necessary vectors have to be formed for this purpose (see Fig. 13.16).
Let’s consider an algorithm of creation of Mathcad document for calculation of
open Brayton cycle with three-stage compression and tree-stage expansion of
working fluid.
In this case block of input data is formed similarly as in the previous example
(simple open Brayton cycle). But additional parameters characterizing three-stage
compression with intermediate cooling—pressure increase during compression
174 V. Voloshchuk
Fig. 13.17 Block of initial data and sketch of the open Brayton cycle with three-stage
compression and tree-stage expansion of working fluid prepared on technology of Mathcad
Calculation Server
before cooling and pressure drops in the cooler; and parameters characterizing
tree-stage expansion—degree of expansion of working fluid before intermediate
heat addition (see Fig. 13.17).
In the given case it is supposed that after partial compression air is cooled to the
temperature of outside. But in general this temperature can also be specified as
initial data. It is also supposed that after partial expansion the working fluid is
heated in the fuel combustors F_Comb2 and F_Comb3 to the same temperature t3
as in the F_Comb1. But again these temperatures can also be specified in initial
data.
The total work of the compressor is equal to the sum of works during com-
pression before and after air cooling (see Fig. 13.18).
Pressure drops in intermediate coolers are taken into account (see Fig. 13.19).
Calculation of the first fuel combustor (F_Comb1) and gas turbine (GT1) is the
same as in the simple gas turbine engine and described above (see Figs. 13.7, 13.8,
13.9, 13.10, 13.11 and 13.12). It is only necessary to indicate correctly the GT1
inlet and outlet pressures (see Fig. 13.20).
13 Calculation of Gas Turbine Engine Cycle 175
During calculation of the second and third fuel combustors (F_Comb2 and
F_Comb3) it should be remembered that fuel and combustion products after GT1
and GT2 respectively are supplied into these fuel combustors. That’s why calcu-
lation of F_Comb2 and F_Comb3 is performed taking into account that oxidant for
fuel burning in these combustors is a mixture of excessive air and combustion
products leaving respectively GT1 and GT2 (see Figs. 13.21 and 13.22).
13 Calculation of Gas Turbine Engine Cycle 177
Fig. 13.23 Evaluation of the amount of excessive exhausted gases after GT1 for providing the
required temperature at GT2 inlet
Fig. 13.24 Evaluation of the amount of excessive exhausted gases after GT2 for providing the
required temperature at GT3 inlet
Fig. 13.25 Result characteristics of the open Brayton cycle with three-stage compression and
tree-stage expansion of working fluid
balance. The heat with mass of combustion products and excessive amount of
exhausted gases after the GT1 are outlet terms of the F_Comb2 heat balance.
Evaluation of the amount of excessive exhausted gases leaving GT2 per 1 kg of
fuel is performed in the similar manner (see Fig. 13.24).
13 Calculation of Gas Turbine Engine Cycle 179
Fig. 13.26 Open Brayton cycle with three-stage compression and tree-stage expansion of working
fluid on the Ts-diagram
Figures 14.1, 14.2 and 14.3 shows Mathcad calculation of thermal efficiency for the
simplest ideal combined (binary) cycle, i.e. steam turbine cycle with use of
superheated steam. At this unit the boiler burner is replaced with gas turbine
(Brayton cycle) equipped with air compressor (c), combustion chamber (cc), gas
turbine (gt) and electric generator [26].
The initial part of calculation in Fig. 14.1 contains a scheme of the combined
cycle plant. There is a T, s-diagrams of thermodynamic cycles (gas turbine and water
steam cycles) on the site of the book used for the calculation. This diagram is “dead”
in the sense that changing of initial plant parameters (see assignment operators to the
right and to the left from the cycle scheme in Fig. 14.1) does not result in corre-
sponding changes in the diagram.
Fig. 14.3 Calculation of thermal efficiency for the simplest ideal combined cycle
It was done this way in the calculation shown in Fig. 14.1. Link reference was
made to the file H2O.xmcdz stored in the author’s site before calculation of specific
enthalpy h6 and specific entropy s6 of the superheated water steam. We have
discussed it in Chap. 4, where similar calculation was created. After such link
reference the functions relevant to water and water steam properties became
viewable. As it was mentioned earlier, files and function names related to water and
water steam properties are provided with wsp prefix.
Similar way “cloud” link reference was made to the file GAS.xmcdz which stores
functions with wspg prefix computing thermodynamic properties of the second
working medium of our binary cycle (see Fig. 14.2). In this case calculation is more
complicated. As a matter of fact, the first working medium (water and steam) as it
circulates in the cycle changes only its parameters and phase state. The working
medium of the second cycle at the binary plant also changes its composition: air—air
mixture with fuel—flue gases. The author’s site contains calculations of gas turbine
cycles with various degrees of complexity taking into account changes in the
working medium composition, initial air humidity degree and other aspects. One can
find one of such calculations in the site at the following address http://twtmas.mpei.
ac.ru/mas/Worksheets/orlov/gases/simple_gtu.mcd and in the Chap. 13. The cal-
culation shown in the Fig. 14.2 has as the working medium air with standard
composition (see the first assignment operator in the Figure: wm: = “Air”). One can
put symbols of other gases (H2, N2, O2 etc.) and gas mixtures into the given string
constant.
Figures 14.1 and 14.2 show calculations of thermal efficiencies for the steam
turbine cycle and gas turbine cycle (Brayton cycle) for the case when they are
separately operated with initial parameters shown in the initial part of Fig. 14.1.
Figure 14.3 provides for calculation of binary cycle thermal efficiency (see its
scheme in the Fig. 14.1), where part of gas turbine exhaust heat is used in the heat
recovery steam generator for generation of main superheated steam (point 6) out of
feed water (point 9). The control over theoretical possibility of this process is
executed by if operator inserted into the gas turbine calculation (see the fifth string
in the bottom part of Fig. 14.2): temperature of exhaust gases from the gas turbine
T4 shall be higher than the temperature of main water steam T6, fed to the steam
turbine. Otherwise, variable h4 (specific enthalpy of gas leaving the gas turbine) is
assigned string value (“Error”) instead of numeric value thus interrupting further
calculation by error message.
The calculations shown in Figs. 14.1, 14.2 and 14.3 introduce measuring units of
water and/or water/steam kilogram (kgws) and gas kilogram (kgg). This type of
calculation with separated physical quantities was discussed in the Chap. 2—See
Fig. 2.24.
The book site contains calculations of more complicated binary cycles which
diagrams are presented in Figs. 14.4, 14.5, 14.6, 14.7, 14.8, 14.9, 14.10, 14.11,
14.12 and 14.13. Internet addresses of the calculations are given in the address
strings and on the site of the book.
14 Calculation of Combined (Binary) Cycle 185
Fig. 14.5 CCPP with single pressure heat recovery steam generator
Fig. 14.6 CCPP with single pressure heat recovery steam generator with reheat and multistage air
compression
14 Calculation of Combined (Binary) Cycle 187
Fig. 14.7 CCPP with double pressure heat recovery steam generator
unit will not provide for any special cost saving. Despite of the fact that fuel
consumption for kwhr production is reduced due to the higher combined cycle
efficiency, all financial profit will be used for paying for after-sales servicing of gas
turbine supplied by foreign company (and the above servicing may be interrupted
any time due to the sanctions). Besides the turbine itself is quite expensive as it is a
customized product unlike aero gas turbines representing mass production articles.
Customized products are always expensive and prankish. It’s no coincidence that
more and more often voices are raised against application of high capacity gas
turbines at power plants.
188 V. Ochkov et al.
Fig. 14.8 CCPP with double pressure heat recovery steam generator and application of two-stage
feed water heating in the economizer
Fig. 14.11 CCPP with multistage air compression and double pressure reheat heat recovery steam
generator
190 V. Ochkov et al.
Fig. 14.12 CCPP with triple-pressure reheat heat recovery steam generator
Fig. 14.13 Comparison of gas turbine power plant and steam turbine power plant with CCPP
Chapter 15
Otto Cycle or What Is Behind
the Simplified Formula
Abstract This chapter explains how the usual simplified formula to calculate
thermal efficiency of an internal combustion engine is appears. Described Internet
sites with the “reverse” design challenge—you must not to enter raw data and get a
response, but you must enter formulas for calculations and verify their correctness
(educational task).
Otto cycle is one of the thermodynamic cycles used in the internal combustion
engines. It comprises (in its ideal version) two isochors (combustion of air-fuel
mixture in the engine cylinder and gas exhaust) and two isentropes (compression of
fuel-air mixture in the cylinder and expansion of gases in a working stroke of a
piston). In reference books and textbooks on thermodynamics one can find a for-
mula used to calculate the Otto cycle thermal efficiency. Figure 15.1 displays a
fragment of a web page with this formula for its interactive use: one can change the
compression ratio r (ratio of the cylinder maximum volume to its minimum vol-
ume), Cp and Cv (isobaric and isochoric) heat capacity of the working fluid in Otto
cycle linked to the universal gas constant R), and obtain a new ηt value to fix the
points on the efficiency/compression rate function diagram [58].
be taken outside the integral sign and we take it without a constant, and then
multiply by this constant. Figures 15.3 and 15.4 show how Mathcad symbolic
processor did that.
The specific entropy calculation is of interest in terms of measuring units used
for temperature and pressure (see Chap. 2). Mathcad symbolic processor gave a
response containing a difference between the logarithms. To use this formula for
numerical calculations, it is necessary to recall that the difference between the
logarithms is equal to the logarithm of the fraction (see Fig. 15.3), and adjust the
formula.
194 V. Ochkov and K. Orlov
Fig. 15.5 The site with the calculation and graphical display of the Otto cycle
Fig. 15.6 Start of the test on knowledge of the calculation of the Otto cycle
Fig. 15.7 End of the test on knowledge of the calculation of the Otto cycle
15 Otto Cycle or What Is Behind the Simplified Formula 197
If, calculation in Fig. 15.2 demonstrates all the formulas, then the website in
Fig. 15.7, instead of formulas offers empty text fields to be filled with the relevant
formulas to be used to calculate the thermal efficiency of the ideal Otto cycle. If you
fill out all of the formulas and press Ready key, at the end of the online test (see
Fig. 15.7) you will be given an answer (thermal efficiency of the Otto cycle), which
can be compared with the response issued by the simplified formula. Formulas
entered in the fields shown in Figs. 15.6 and 15.7, are written not by the operator
(via operator of the definite integral, for example, as shown in Figs. 15.2, 15.3 and
15.4), but in the form of text functions—int, / (division), * (multiplication), etc.
Chapter 16
Calculation of Pressure Losses in the Tube
Abstract This chapter explains how to calculate the pressure loss in the pipeline
using cloud functions for thermodynamic and transport properties of the fluid.
A description is given of finite element methods for the solving of important
engineering problems.
This chapter is designed to show how one can use “cloud” functions to solve a
simple but typical hydrotechnical problem i.e. the problem of the loss of water head
in the horizontal tubing. Such problem occurs in the design of a (cooling) water
supply loop from the cooling tower to the steam turbine condenser with a heated
water return to the cooling tower.
Loss of head (pressure) takes place also in the steam supply tubing from the
boiler to the turbine. Chapter 5 (Fig. 5.1) calculates and shows the process of steam
“delivery” from the boiler to the turbine in hs- and Ts-diagram on assumption that
the initial pressure is 24 MPa and the final pressure is 10 MPa.
In the author’s website there is a web page making possible to calculate the
process of steam expansion (throttling) and to show it in the hs-diagram (Fig. 16.1).
If we set the value of relative internal efficiency of the steam expansion process
equal to 100 %, then we obtain the process of ideal steam operation in the turbine
while steam entropy remains constant (“vertical” expansion of steam). But if we set
the value of relative internal efficiency of the steam expansion process equal to 0,
then we obtain the process of ideal throttling of steam while steam enthalpy remains
constant (“horizontal” expansion of steam). In the actual steam operation in the
turbine the value of internal efficiency of the process is close to 100 % while in the
actual throttling it is close to zero. A visitor of the web site shown in Fig. 16.1 can,
at his(her) discretion, change the value of the relative internal efficiency of the
process and obtain the required steam expansion curves.
Figure 9.9 in Chap. 9 displays a diagram of wet steam throttling making possible
to evaluate the rate of steam dryness in the initial point if one knows the initial
pressure, the pressure and the temperature in the final point.
As we have noted before, a great number of computer programs have been
created to calculate tubes of various degree of complexity and details. Such pro-
grams operate on the principle of a “black box” into which one “puts” initial data
array, “closes the cover of the box” (by pressing the key “Recalculate”), “opens the
16 Calculation of Pressure Losses in the Tube 201
cover” and “takes out” of the box the solution i.e. the parameters of the designed or
operated tubing. But as it has been already noted in the introduction, one always
wants to know what it is inside of such “black box” and whether the calculations are
being correctly made. In addition, it is useful to slightly open the cover of such
“box” for the purposes of education and self-education—for the study of mathe-
matical models built in these or those calculation programs. Here we will now
consider such particular mathematical model!
In order to solve the problem of the loss of head in a tube an engineer will recall
or find in the reference books (in the printed reference books or in the Web-versions
of reference books), a set of the corresponding formulae and rules of their appli-
cation—formulations. In addition, the engineer must know some properties of water,
in particular for this problem—its kinematic viscosity ν. Values of such water
parameters are given in many various reference books in the form of tables, graphs
or empirical formulae related water viscosity to its temperature. As a rule pressure is
assumed to not affect this property. However, data taken from the tables and for-
mulae of rather fat reference books can considerably differ. In addition, as we have
already noted many times, when transferring figures from the tables to the particular
calculation an engineer can make a mistake connected with incorrect interpolation
and even with incorrect entry of digits on the keyboard and/or incorrect interpre-
tation of measurement units of viscosity and their multipliers. All these things can
complicate and hinder calculations and increase the possibility of errors.
Dynamic viscosity of water is calculated by the formulation developed and
approved by International Association of Properties of Water and Steam (IAPWS—
www.iapws.org) with which the author of this book closely cooperates. Based on
the IAPWS formulations reference books are issued all over the world with the
tables keeping discrete values of water and steam parameters as functions of
temperature and pressure. Such tables are published in our country as well. The
most recent reference book [23] is added with the Web-version (http://twt.mpei.ac.
ru/rbtpp) to make the use of the reference book easier and faster. One of the pages
from the reference book at the URL address http://twt.mpei.ru/MCS/Worksheets/
rbtpp/tab9.xmcd is shown in Fig. 16.2.
Fig. 16.2 Web page with the on-line calculation of dynamic viscosity of water
202 V. Ochkov and K. Orlov
By getting into the Web site shown in Fig. 16.2, it is possible to change the
temperature and/or pressure values selecting the required measurement units of
these initial values, set the number of significant figures in the response, press the
key Recalculate and obtain not only the required value of water or steam dynamic
viscosity but as well the value of relative error of this figure that is different in
different pressure and temperature ranges. It is possible to slightly open such “black
box” for calculations and to learn from this Web-version of the reference book the
formulae used for calculations of water and/or steam dynamic viscosity.
As we have already noted, in Moscow Power Engineering Institute at the Power
Plant department they have created the new Internet technology to operate with
water and steam properties excluding the manual transfer of data from the printed
[23] or Web-version of the reference (Fig. 16.2) book.
The reference book [23] compiles not only on-line calculation of water and/or
steam properties but also the corresponding programmed functions that one can
directly use in the calculations by setting e.g. not water dynamic viscosity but water
temperature and water. The beginning of a list of such functions (“cloud” functions—
Web cloud-based functions) is shown in Fig. 16.3.
Such list contains the water density and dynamic viscosity values required to
calculate the loss of water head in the tubing. In order to make the water density
value versus temperature and pressure values accessible (or visible as programming
specialists say) in the Mathcad-calculation, it is necessary to know the URL of this
“cloud” function http://twt.mpei.ac.ru/rbtpp/MC-WSP/M15/wspDYNVISCTD.
xmcdz. For this purpose it is enough to point the cursor to the necessary func-
tion (see the underlined names in Fig. 16.3), right-click and choose the command
Properties. Such action will open the following window with the URL of the
required function (Fig. 16.4).
Fig. 16.4 URL of the “cloud” function “Water and steam dynamic viscosity”
It is necessary to copy this URL address and insert it in the dialog box Insert
Reference in the Mathcad program which calculates the loss of water head in the
pipeline or makes any other calculation requiring the knowledge of the water
density value ρ (Fig. 16.5).
The same way it is possible to insert in the calculation the function returning the
water dynamic viscosity μ as a function of temperature and pressure. The URL of this
“cloud” function is http://twt.mpei.ac.ru/rbtpp/MC-WSP/M15/wspDYNVISCTD.
xmcdz.
The value of kinematic viscosity ν is determined by the values of water density
and its dynamic viscosity in the calculation shown in Fig. 16.6. In the previous
calculation the internal tube diameter value d was entered and the value of water
velocity in the pipe v was calculated. These three parameters (velocity, diameter
and viscosity) are used to calculate very important dimensionless criterion—the
criterion of liquid flow in the tube i.e. Reynolds’ number Re, which determinates
not only the tube inner surface roughness but as well the desired value of the loss of
water head due to wall friction.
Reference books show the relationship between the relative dimensionless
resistance factor and the mode of liquid flow in a round tube λ either in the graphs
or by a set of formulae. Figure 16.7 shows a fragment of the Web page at URL
http://twt.mpei.ac.ru/MCS/Worksheets/Hydro/La-De-Re.xmcd displaying the
“live” nomograph (Nikuradse curve) relating the resistance factor λ under our
consideration with the Reynolds’ number Re and the relative roughness of surface
Δ (relationship between medium height of asperities (roughness) of the tube inner
surface ant its internal diameter).
16 Calculation of Pressure Losses in the Tube 205
A visitor of the web site (Fig. 16.7) can change the value Δ and Re, press the key
Recalculate and obtain (read) not only the desired figure but also can see the
character of this relationship in the most “interesting” area—in the transition area
from laminar flow to turbulent flow—the area where it is not recommended to
operate tubing. The areas to the left (liquid laminar flow) and to the right (enhanced
turbulence) of the curves shown in Fig. 16.7 are calculated by the following for-
mulae in Fig. 16.8. At the calculation portal where the web sites described in the
books are located, there is a fluid and gas dynamics reference book.
In this Web-reference book the “cloud” function is stored that returns the value
of the required friction coefficient as a function of the Reynolds’ number in the wide
range and relative roughness of the internal surface of the round tube. This function
is shown in Fig. 16.8 and entered in the calculation shown in Fig. 16.6 by means of
the above tool of the reference to the “cloud” function.
In the literature one can find various formulae for calculation of the friction
coefficient as a function of Re and Δ (Fig. 16.8 and table in Fig. 17.10 in Chap. 17).
There are other formulae for these ranges and one can take, in principle, any
formula from this set of formulae without thinking about their accuracy. The matter
is that the initial parameter i.e. roughness of the tube is very roughly estimated.
One can only approximately measure roughness of the tube. Usually it is done in
the process of hydraulic tests of tubes. Reference literature gives this important
parameter within some limits depending on the material of the tubes (metal, plastic,
ceramics, copper etc.), the method of their manufacturing (a welded tube, a
seamless drawn tube etc.) and other manufacturing features. In addition, roughness
of the tube inner surface can change in the process of its operation or at a standstill
due to corrosion and erosion or due to deposits in the tubes. In principle, when
calculating the loss of head one shall not enter absolute values of the relative
roughness Δ as initial data but certain linguistic constants (“a seamless drawn tube”,
“a plastic tube” etc.) which can be used by the method of fuzzy sets for evaluation
of the value Δ in a certain range determining along with other uncertainties (i.e.
uncertainty of the range of Reynolds’ number values and the formulae for calcu-
lation of λ in Fig. 16.8) the expected range of possible change of the calculated
value of the loss of head.
In the calculation shown in Fig. 16.6 water parameters along the length of the
tube are determined by the given temperature of 16 °C and the pressure equal to
5 atm. This is rather rough approximation which we will try to solve now.
Water pressure along the length of the tube varies, of course. As a rule, it goes
down along the liquid flow but it can also rise due to hydrostatic pressure if the
liquid in the tube is fed from bottom to top. But properties of liquids, water in
particular, weakly depend on pressure. They greatly depend on temperature. And
how will the liquid temperature and pressure change along the water flow? If water
at a temperature higher than ambient air temperature is pumped by the tube, then it
is necessary to consider in the calculation the loss of heat through the tube wall. The
attempt to solve the problem is shown in Figs. 16.9 and 16.10.
The initial data will be the following values i.e. internal diameter of the tube d,
length of the tube L, tube wall thickness ΔT, relative roughness of the inner surface
of the tube Δ, water temperature t and pressure p at the inlet of the tube, water mass
flow Q, tube wall material heat conduction, λ, heat-transfer coefficient from water to
the inner surface of the tube α1 and from the outer surface of the tube to ambient
environment α2 as well as ambient temperature tout.
For calculation the tube is divided into n sections (finite elements). Calculation
of the initial (zero) section is shown in Fig. 16.9.
This figure shows the successive calculation of the following values after the
input of the initial data: area of the inner section of the tube F, length of the tube
finite element Δl, its inner r and outer R radius. During the program debugging n
(the number of finite elements of the tube sections) can be set as small (10–100) and
after the program is ready it can be increased to 1000–10,000.
16 Calculation of Pressure Losses in the Tube 207
Fig. 16.9 Calculation of the initial section of the tubing (Var = 1: horizontal tubing)
After the value n is set, the length of the tube finite element Δl is calculated.
Then it is necessary to enter the initial (zero) elements of vectors storing the
following values on this finite elements: T—water temperature, P—water pressure,
N—specific loss of heat (enthalpy), ρ—water density, v—water velocity, h—water
specific enthalpy, μ—water dynamic viscosity, Re—Reynolds number, ΔH—loss
of head and ΔP—loss of pressure.
The formula determining the loss of head ΔH contains the summand with if
function and additional parameter Var i.e. variant. There are three variants here:
• Var = 1: horizontal tubing;
• Var = 2: vertical tubing, pump at the top;
• Var = 3: vertical tubing, pump at the bottom.
Here we consider some extreme cases. Actual tubing can contain various sec-
tions hugging the terrain and/or the tubing geometry at different process plants.
The most complicated formula is the formula for determination of water specific
loss of heat (enthalpy) N. This is a standard formula developed for solving the
problem of steady heat conduction of the coiled cylindrical tubing if there are no
internal heat sources. The boundary conditions of this problem will be the values of
the heat-transfer coefficient on the inner α1 and outer α2 surfaces of the tube.
Determination of α1- and α2-values is an individual and not least complicated
problem. For the time being we will set these values as constants. However, in the
208 V. Ochkov and K. Orlov
process of calculation we can enter in it the formulae for determination of α1- and
α2-values depending on the mode of water flow in the tube and other factors.
Figure 16.10 shows Mathcad operators entering the rest elements of the above
mentioned vectors.
Figure 16.10 shows the enumeration cycle of all tubing sections starting from the
first one. Zero (initial) tubing section has been calculated before (Fig. 16.9) outside
of the cycle including for the purpose of debugging of the calculation program.
The cycle in Fig. 16.10 shows the calculation of the successive i section of the
tubing with the use of the data for the previous (i – 1) section. The matrix M is
added with three columns i.e. data for change of water pressure P and water
temperature T as well as vector L storing the values of distances from the beginning
of the tube to the point i. During the derivation (ref. to the last operator in the
program in Fig. 16.10 with the function augment1) these data are to be firstly
1
This function built in Mathcad joins individual vectors in a matrix where these vectors become
columns.
16 Calculation of Pressure Losses in the Tube 209
Fig. 16.11 Water temperature and pressure curves in the horizontal tubing
deprived of dimensions of pressure, temperature and length and then one has to
return them in new vectors P, T and L generated outside of the program (ref. to the
last three operators in Fig. 16.10). As we have already noted many times, this is
connected with the main disadvantage of Mathcad 15 (as the calculation shown in
Figs. 16.9 and 16.10 has been created in this program) i.e. impossibility of storing
data of different dimensions in the arrays.
Figures 16.11, 16.12 and 16.13 show temperature and pressure curves in our
tubing located horizontally (Fig. 16.11), vertically with water feed from bottom to
top (Fig. 16.12) and from top to bottom (Fig. 16.13).
The curves also show the values of water saturation at the initial and final
pressure in the tubing calculated by means of the function wspTSatP. And this is
being done for some reason. The fact is that the calculation shown in Figs. 16.9,
Fig. 16.12 Water temperature and pressure curves in the vertical tubing (pump at the top)
210 V. Ochkov and K. Orlov
Fig. 16.13 Water temperature and pressure curves in the vertical tubing (pump at the bottom)
16.10, 16.11, 16.12 and 16.13 is the attempt to simulate the process of pressure hot
water feed by a rubber hose for fire extinguishing. At the end of the hose there is a
nozzle through which the water that has passed into wet steam cloud extinguishes
fire.
When feeding such water it is important to avoid its early boiling inside of the
rubber tube. The most risky case here is water feed from top to bottom. If in the two
first cases (horizontal tubing and vertical tubing with water feed from bottom to top)
the loss of pressure can be not so considerable, then in the third case the loss of
head of water can cause sharp loss of pressure and flashing of water in the rubber
tube.
On-line calculation of water parameters in the tube can be made at http://twt.
mpei.ac.ru/MCS/Worksheets/Hydro/TubeMChS.xmcd.
It is also possible to make on-line evaluation of the change of temperature inside
the tube wall. For this purpose it is enough to get in the Web site http://twt.mpei.ac.
ru/MCS/Worksheets/Thermal/Heat-Flow-1radius.xmcd, enter the required initial
data and obtain the result (Fig. 16.14).
One would always like to simplify the key formula for calculation of the loss of
head (this formula is in the box in Fig. 16.6) by transferring e.g. constants 2 and g
(gravity acceleration) from the denominator to the coefficient of resistance, λ. But
we cannot do this for several reasons. Firstly, the gravity acceleration g
(9.81 m/s2) is not a constant, strictly speaking, but a function of other values and,
mainly, a function of the altitude above sea level and terrestrial latitude. If one
makes reference in Mathcad-document to the “cloud” function at http://twt.mpei.ac.
ru/TTHB/g-h-psi.xmcdz, then in the Mathcad environment the constant g will turn
into the function g(h, φ), where h is the altitude above sea level and φ is the
terrestrial latitude (0°—the Equator, 90°—Earth poles). Besides this, one should
remember that some time in future tubing can be designed and built e.g. in the …
Moon or Mars where the value of g is absolutely not 9.81 m/s2! And this is the
16 Calculation of Pressure Losses in the Tube 211
interesting question: is our formula applicable for determination of the loss of head
in a water supply line or fuel supply line at the space station where g = 0?
Incidentally, in many calculations g is taken equal to 10 m/s2 and water density is
taken as 1000 kg/m3 to simplify the counting and with due regard to the fact that
calculations are rather rough. We also can use the loss of head we found in Fig. 16.2
for rather rough estimation of pressure drop in the environments required for
pumping water by the tube. For this purpose one shall divide the head value by ten.
So, the constant g cannot be removed from the formula of calculation of the loss
of head in the tubing and transferred to the coefficient λ. And what can we do with
“two” in the denominator of this formula? It also cannot be transferred to the
coefficient of resistance, λ, but for other reasons. Looking at two and velocity
squared in the formula of the loss of head, one can recall an old joke. A student was
in a hurry to the examination in physics, he became thoughtful and… suddenly run
into a post. “That’s good that it’s divided in half!” exclaimed the student. “What’s
in half?” his comrade asked him. “mv squared is divided in half!” the student
answered.
212 V. Ochkov and K. Orlov
Fig. 16.15 Solving the problem of the temperature profile in the cylindrical wall (beginning)
2
The function odesolve can solve some boundary value problems but by no means all. There are
special tools for boundary value problems in the Mathcad environment.
16 Calculation of Pressure Losses in the Tube 213
Fig. 16.16 Difference equation of heat conduction in the cylindrical wall and its solving by means
of Mathcad symbolic mathematics
Fig. 16.17 Solving the problem of the temperature profile by the shooting method
step-by-step approach (shooting method) i.e. build the temperature profile inside of
the cylindrical wall. This operation is shown in Fig. 16.17.
Figure 16.17 shows the given variable n i.e. a number of cylindrical sections into
which our wall is divided. The variable h is thickness of a section. Then we create the
vector R storing the radii of wall sections. We have the value of the first element of the
vector T1 equal to the wall temperature on its left (inner) edge Tc1 but we have no
value of the second element of this vector T2 for the recurrence formula. So we create
it artificially by means of the operator T2: = T1 − 0.01 K e.g. and see what the value of
the last element of the vector T (Tn) is. If it is not equal (approximately) to the wall
temperature value on the right (outer) edge Tc2, then we change the value T2 and
make the following shoot i.e. filling the elements of the vector T: Ti+1: = f(Ti−1, Ti).
In addition we keep control of the diagram in Fig. 16.18 showing three moments of
trial with different values of T2. This figure also shows the function linterp built in
Mathcad making possible to create the function T(r) by two vectors R and T by the
214 V. Ochkov and K. Orlov
Fig. 16.18 Result of solving the problem of the temperature profile by the shooting method
16 Calculation of Pressure Losses in the Tube 215
piecewise linear interpolation method, the above function is used to build the desired
temperature profile in the cylindrical wall.
Note. Figure 16.17 does not show the formulae for Ti+1 in full. This is not only
to save the space in the book but due to the fact that these formulae may be of no
interest for study and/or manual entry in the calculation. One of the two formulae to
be transferred to the calculation can be determined by the trial-and-error method. It
is also possible to attempt to simplify the initial difference equation manually by
substituting e.g. Ri+1 + Ri for 2Ri + h, and Ri + Ri−1 for 2Ri − h etc. but it’s no
good. Or rather put it that way: Simplifications and substitutions are worth trying to
do if the Mathcad fails to solve the initial equation. But it is necessary to remember
that manual substitutions and simplifications are fraught with mistakes e.g. loss of
brackets etc. True, the mechanism of measurement units used in the calculation
allows discovering many of such mistakes preventing from e.g. adding temperature
to length.
The shooting method (method of successive approximations) can be automated
by means of Mathcad programming by implementing e.g. the half-interval method
(shooting method). Such operation is implemented e.g. in Fig. 18.11 in Chap. 18.
But the manual problem solving is good so much that it can be done in the
environment of Mathcad Prime free version i.e. Mathcad Express.
If we increase the values R1 and R2 leaving the cylindrical wall thin, then it will
tend to a flat wall in shape. Figure 16.19 shows the solving of the problem of
calculation of the temperature profile in the flat endless wall with internal sources of
heat for the case when wall material heat conduction does not depend on temper-
ature. In this case the constant λ is taken out of the differential sign and the
differential equation takes on the standard for Mathcad form. Here the function
odesolve built in Mathcad is successfully used for the problem with initial condi-
tions and the above described shooting method but when establishing not the
second value of the temperature vector (Fig. 16.17) but temperature gradient on the
left edge of the wall. Figure 16.19 shows the shooting method very well when at
the right end of the diagram one can clearly see “undershoot” or “overshoot” if the
temperature gradient T’ is established incorrectly from the “artillery position” from
the left wall edge.
But let’s go back to the name of the study. The formula of the loss of head in the
tubing displays the fundamental physical law connecting the change of kinetic
energy of liquid flow due to friction on the walls and conversion of this energy to
less valuable i.e. to heat energy. In the horizontal tubing (our simplified problem)
the potential energy does not change but the kinetic energy “mv squared divided in
half” changes. So the “two” in our formula of the loss of head is rather appropriate:
it is difficult to understand the “physics” of the problem without it.
In connection with our above reasoning about gravity and “two” in the formula
of the loss of head one is put in mind of an interesting theoretical issue regarding
rivers i.e. “water-supply lines” providing the mankind, among other thing, with
electric energy. For which purpose is hydroelectric dams built? The overwhelming
majority of rather competent engineers and specialists can answer approximately as
follows. A hydraulic dam is built on the river to increase the water level in the river
216 V. Ochkov and K. Orlov
and create water head before the hydraulic turbine generating electric power. But a
dam is a purely passive structure and it is not a pump boosting the head! Our simple
problem under consideration in this book will help us to answer this question
correctly. A hydroelectric dam is built in order… to reduce the average water
velocity in the river due to the increase of the river bed cross-section. On plains
rivers low but long dams are built while on mountains rivers dams are high and
narrow. Due to this the average water velocity in the river reduces sharply and
therefore, losses due to water friction on the banks and bottom reduce. And this is
the source of energy for the turbines of a hydro power plant generating electric
power.
There is another way of generating electric power from the river which is now
used in increasing frequency, particularly on mountains rivers and brooks without
construction of any expensive and dangerous dams in seismic zones. A tube is laid
parallel to the river bed with a hydro turbine installed at the lower end of the tube.
16 Calculation of Pressure Losses in the Tube 217
Water or a part of it goes down not the rough river bed but by the smooth tube. Due
to this the potential energy of water (it should be noted once more that this energy
does not change after the damming) determined by elevation difference does not
dissipate so intensively and converts to electric power sufficient to provide a small
settlement with lighting. Such mini- and micro hydro power plants are widely
spread e.g. in Norway. Incidentally, such “power plants” are sometimes “built” at
the ends of different tubes as well. The specific example: a town is supplied with
natural gas by a high pressure gas pipeline. While household users need the gas
under very low overpressure otherwise problems of gas leakages indoors will occur.
Excess pressure can be released by gas throttling by means of pressure reducing
valves with loss of such gas flow energy. And there is another energy efficient
solution. It is possible to install a turbine generator at the end of such gas pipeline to
convert gas overpressure to electric power. Some gas turbines at large power plants
are equipped with such mini-power plants. The author has long ago seen an
interesting hydrotechnical facility in Lithuania. In one ravine water went down with
the elevation difference 2–3 m. These elevations were connected with a tube having
a piston pump with a piston drive at the lower end of the tube. That pump delivered
water to a house located much higher than the upper end of the tube. Such
hydrotechnical facilities (i.e. hydraulically activated pumps) equipped with two
pistons different in diameters located on the beams with different shoulders have
been widely used in due time for water supply in Switzerland rich in rivers and
brooks with large elevation difference.
Mention may be made of one more exotic but quite implementable method of
energy generation from a river i.e. a stream of fresh water flowing into salt seas and
oceans. There are many methods of production of fresh water from salt water e.g.
distillation, freezing, ion exchange, reverse osmosis, electrodialysis… Sea (salt)
water and electric power come to electrodialysis plants producing fresh water
(demineralized water) and concentrate. But such plants are reversible. If salt water
and fresh water is fed to such plants separately, such plant will generate electric
power due to mixing of waters of different salinity. Such electrodialysis power
stations are rather appropriate in the beds of large rivers flowing into seas and
oceans. At present the similar pilot and semi-commercial plants exist and the issue
of construction of such actual power plants is addressed.
A reverse osmosis plant is also a system to which salt water and electric power
are fed and from which two water steams go out: demineralized water (permeate)
and concentrate. But the principles of action of electrodialysis and reverse osmosis
plants are different. In the electrodialysis plant electric current is directly fed to
special electrodes forcing ions (cations and anions of salt admixtures in water) to
move, separate by means of special membranes and finally to be removed from
water. And in the reverse osmosis plants electric power is used for driving the
pumps to boost the pressure of salt water to the value (osmotic pressure and higher)
sufficient for water molecules to go (squeeze) through the finest pores of the special
membranes. As this takes place, aquated cations and anions remain in salt water and
are removed from the plant with blowdown concentrate. In principle, a reverse
osmosis plant can also be forced to generate power by feeding two water streams of
218 V. Ochkov and K. Orlov
different salinity into it. But here there will be no direct power generation as in the
electrolysis plant. In the reverse osmosis power plant i.e. in simply osmotic power
plant one has to convert the energy of water head (osmotic pressure) to electric
power by means e.g. of the above mentioned hydro turbines. Due to this no con-
struction of such power plants is now under consideration.
Now water treatment systems for heat and nuclear power stations in Russia are
subject to upgrading. Particularly, clarification units for calcification and coagula-
tion and ion exchange units are being replaced for ultrafiltration and reverse
osmosis membrane units. This results in sharp reduction of consumption of
chemicals (lime, coagulating agents, acids, alkali) for water treatment and amount
of waste water. But it also results in increase of power consumption for driving the
pumps pushing water through the membranes. In this connection in order to save
power it is reasonable to consider the issue regarding the use of a small steam
turbine instead of electric motor to drive such pumps. This is the principle of
operation of the pump delivering feed water to the boilers of the heat power plant.
Finally, one newspaper published the information that scientists proved that a
square section tube has less hydraulic resistance than a round tube. Such square
sections can under otherwise equal conditions deliver more water, gas, oil to cus-
tomers… Besides that such tubes do not roll asunder when stacking and it is easier
for repair robots to move inside of them… And this was April fool’s joke.
Nevertheless one tube-rolling mill wrote a formal letter to the newspaper editorial
office signed by chief engineer informing that they are ready to manufacture such
innovation tubes. That’s what can happen if one does not know the elements of fluid
dynamics set forth in this chapter. A tube is round because a circle as opposed to
other plane figures (square, rhomb, triangle etc.) has the least perimeter length at the
fixed surface.
Chapter 17
Cogeneration (CHP), Trigeneration
(CCHP) and Quadrogeneration (CCHPI),
or How Much of Mathematics Is
Contained in Thermal Engineering
Valery Ochkov
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
exceptions, were not connected to the existing municipal heating systems but
adapted to the individual “boiler” heat supply. Baku is a city in the south, but they
have quite strong and cold winds in the winter blowing from the Caspian Sea. Use
the household electric heaters for heating the apartments during this season is rather
expensive, and indeed electric grids would hardly sustain such a peak, and natural
gas is relatively cheap and available. However near the newly built houses there are
no gas boiler houses to supply heat to a house or an entire neighborhood, as is done,
for example, in our new suburbs in the vicinities of Moscow. In Baku the new
occupants have to supply their own heat. They have to install small gas boilers in
their new houses (Fig. 17.1).
These small boilers (in Baku they are called “combi”), intended for individual
cottages, are installed on the balconies of apartment buildings, where they can be
easily serviced, repaired or replaced. Figure 17.1 also represents a compressor-
condenser module of household air-conditioners (split units).
Of course, more rational and more “aesthetic” heating method is the construction
of a centralized gas boiler facility to be located near the house, in its basement or on
the roof. These boilers, by the way, are used to equip some of the new buildings in
Moscow and several other Russian cities. But in Baku this is not done for several
reasons, the main of which is that the builders are trying to minimize the costs for
the construction of a new house, providing for the future tenants the opportunity
(and obligation) to finish the construction of their houses by themselves. In Baku,
it’s not just “wallpaper, tiles and sanitary ware” (as in Moscow), but also heating.
By the way, to the “south of Baku”—in the Middle East, Central and South Asia,
and Africa, they often refuse from centralized power supply, and provide for
electric generators with internal combustion engines, which may be added to a
house as an annex or placed on a balcony of an apartment building.
Fig. 17.1 Small boilers on the balcony of an apartment building to heat individual apartments
222 V. Ochkov
This, “technical and economical nonsense” and to put it mildly, irrationality, can
be also seen in Russian cities where they build nice modern apartment buildings of
business class, which after tenants move in of are spotted with air conditioners.
Although it was possible to provide at once in these houses a centralized installation
of industrial air conditioners or, more rational systems for energy efficient
year-round hot water supply, heating in winter and cooling the apartments during
the hot season.
One can certainly argue about the economic feasibility of the centralized heating
and cooling system for modern apartment buildings. In Baku, for example, the cold
time with winds does not last long, and many apartments on “emergency” basis are
heated by electricity, rather than by hot water. By the way, that is done in Moscow
during the cold weather period when the “heating radiator is still not on.” But the
fact that all of these air-conditioners, small boilers (especially electric generators)
spoil the look of the house and degrade the environment is indisputable. Air con-
ditioners are generating additional noise and dripping water on passers-by. Now in
some Russian cities it is not allowed to mount air conditioners on the facades of
houses. The new designs provide for special niche on the balcony for room air
conditioners. In some European cities, open or “boiler” fuel combustion is not
allowed, not only for heating but also for cooking purposes. A person cannot even
have a barbecue on the coals in the garden of his single-family house. All this is
done only by means of safe and “green” electricity, without paying special attention
to the economic side of the issue. Use URL http://twt.mpei.ac.ru/MCS/Worksheets/
PTU/Vol-13.xmcd to find an open interactive design comparison of systems for
separate generation and cogeneration of electricity and heat. We can also discuss
quadrogeneration. Some IT-organizations have powerful servers with a large
self-contained power supply facilities equipped with gas turbines or reciprocating
engines of electric generators. Such organization can generate electricity, heat,
cooling and … its main product—information (combined cooling, heat, power and
information) (CCHPI).
A fact which is currently playing against the classical heat supply scheme
(CHP + peak hot water boilers + district heat networks) is as follows. At the time
when district heat networks were designed and constructed there was no such an
intense car traffic in the cities. Now, excavation and repair of main and local district
heating networks lead to huge traffic jams, where standing, or rather, “smoking”
machines “devour” the “equivalent fuel savings” that we receive due to district
heating.
Ideally, the solution should be the following (it has been already implemented in
many “advanced” countries caring about energy saving, ecology and aesthetics of
buildings): a combined cooling and heating plant connected to the mains intended
for the apartments is installed on the roof of a house, near a house or in its
basement. The electricity is generated far from the house—at the modern multiwatt
power plants—thermal ones (running on fossil or nuclear fuel), hydraulic, wind and
other power plants. Also this house-installed cooling and heating plant is addi-
tionally equipped with modern energy saving systems—heat pumps, heat/cold
accumulators (tanks with water, salt solution, or earth beds under the house), solar
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 223
collectors, etc. The apartments in such a house are supplied with electricity through
the wires, and, for example, an aqueous solution of salt—hot in the winter and cool
(with sub-zero temperatures) in the summer. Here one could still add a centralized
vacuum cleaner and a centralized removal of garbage sorted by tenants, but we
dream too much.
On the site of the book one can find interactive, open calculations created in
Mathcad and posted on the Internet, that shows the benefits of co-production of
different types of energy. You can see two examples in Figs. 17.2 and 17.3, the title
(see the Figure) of calculation and its main “result” operator.
Back to the aqueous solution of salt. Here, of course, one can consider other
coolants (freon coolants, information about their properties is available on the
described computation server), but salt is available, relatively inexpensive and,
most importantly, safe (non-toxic). Moscow, for example, is standing on a huge
underground salt lake, and the solution (brine) has long been used by Moscow
power industry to regenerate the filters used for softening the make-up water to
district heating systems.
So, to calculate the network of pipes through which the aqueous solution of
sodium chloride flows, it is necessary to know the physical properties of the
coolant, in particular its density and viscosity.
To obtain this information, the designer can run the online Google search engine,
for example, and make it search using key Density NaCl solution, adding to it the
word Mathcad. Without this word the search engine would reference to websites
where both traditional and quite unnecessary tables are published representing
dependency of density from temperature or constant. And what we need is a
function for calculating and such that it is available (“visible” as the programmers
say) in our Mathcad-calculation.
Fig. 17.2 Calculation of fuel economy on a CHP gas turbine compared to separate power
generation
224 V. Ochkov
Fig. 17.3 Determination of equivalent fuel economy for district heating with the turbine type T
(district heating) compared to separated production of the same amount of electricity and heat
Figure 17.4 shows that the search on the server by a key Density NaCl solution
Mathcad issued links the first of which leads to the calculation site shown in
Fig. 17.5.
On the site shown in Fig. 17.4 one can enter new values for temperature and
concentration of NaCl aqueous solution and obtain a new computed value of the
density of this solution. A visitor to the site can download a Mathcad file with a
function that returns the density of the salt solution depending on temperature and
concentration. But it is fashionable to do otherwise—to provide a link to the
Mathcad-file. To know its address—see Fig. 17.6.
If one right-clicks on the link shown in Fig. 17.5, a list will drop whose last item
stores information about the address of the desired file, see Fig. 17.6. This address
will need to be copied and moved to the Insert Reference window in Insert menu,
Mathcad 15, see Fig. 17.7.
Figure 17.7 displays command Insert|Reference in Mathcad 15 environment,
which enables to make a link to a cloud function called wspDT01MPaNaCl-
Water-Solution.xmcdz (D—density of NaCl water solution as function of tem-
perature T at pressure 0.1 MPa; suffix wsp is explained above).
226 V. Ochkov
By the way, if one will make a search in the Russian version of Google by the
keyword “online computations”, the first item in the list of links displayed will be a
link to the site, one of the pages of which is displayed in Fig. 17.4 here. This
demonstrates the high popularity of and demand for the MPEI calculation server.
As we have repeatedly noted, Mathcad’s environment has a very convenient tool—
a link from the current calculation to another one, to another file. After such a link all
the user’s variables and functions defined in the document to which the link is pro-
vided, get visible in the working document.
These links, as we have said, can be done to files stored not only on your work
station or in the local network, but on the Internet, in the “clouds”, and we are going
to use that when solving the problem of pipeline with heat and cool carrier—NaCl
aqueous solution.
Figure 17.9 shows the calculation in Mathcad 15 environment for a simple
piping system—determination of NaCl solution flow in its individual sections for a
specified geometry of pipes and specified pump head. After entering the input data
(geometric parameters of pipe sections: the length L, inner diameter d, roughness,
temperature and concentration of NaCl aqueous solution) give a link to the two
cloud functions returning the density and dynamic viscosity of the coolant and
stored in the “clouds” on the sites shown in Figs. 17.5 and 17.8. After such links
one can call the functions that returned the density (1140.175 kg/m3) and viscosity
(5.085 mPas) of NaCl solution with concentration 21 % and temperature –15 °C. It
should be emphasized that Mathcad can work with physical quantities and their
measuring units [8], which distinguishes it from other computer calculation pro-
grams—from Excel spreadsheets and programming languages (see Chap. 2).
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 227
Fig. 17.8 Website page displaying physical properties of aqueous solutions of salts
In the calculation shown in Fig. 17.9, there are two hidden areas. They are open
and shown in Figs. 17.10 and 17.11.
The first hidden area (Fig. 17.10) stores a function named λтp and the table itself
in the form of inserted picture for which this function has been created. This
function returns the fluid friction resistance coefficient as a function of Reynolds
number Re and the relative roughness Δ of the pipe inner surface [63]. The table
notes that in the range of Re = 2300–4000 (transient from laminar to turbulent
mode) a friction coefficient cannot be determined. However, a corresponding
function is given 0.03. This is done to perform the iterative calculation using this
function, and the “little trick” will be mentioned below. Incidentally, the described
site contains an interactive calculation (http://twt.mpei.ac.ru/MCS/Worksheets/
Hydro/La-De-Re.xmcd) of the friction coefficient in the transient zone (“The
Nikuradze spoon”).
The second hidden area of the calculation, open and shown in Fig. 17.11, stores
a function called Δh, which is used to calculate the differential (loss of) pressure of a
pipeline depending of its length L, inner diameter d and fluid volumetric rate q. In
this function, the program lines are a matrix with two columns and one line
(“horizontal” vector), whose first element is a comment, and the second element is a
formula. If the pipes had different roughnesses (failure of one of the assumptions),
the function Δh should be supplemented by another argument Δ. Function Δh uses
a function called λfriction, shown in Fig. 17.10. These are standard hydrodynamic
formulas, which can be found in the handbooks or on the Internet. We have already
used them in the previous chapter. The table shown in Fig. 17.10, was indeed found
on the Internet, copied and transferred as a picture into our calculation as a com-
ment. In this table, the fluid flow transient value λfriction is not defined. But in the
function λfriction it equals to 0.03. We repeat, that this is done in order to be able to
228 V. Ochkov
at solving such systems of equations will be unable to solve them. One will need to
use special “smart” techniques for solving such problems and highly specialized
software, created specifically for the calculation of pipelines. But these programs
will need information on the properties of liquids, available through a “cloud”
technique, as described above in this chapter and in other places in the book.
Blocks for solving algebraic equations in Mathcad’s environment can return not
only specific numerical values, but also be used to create user-defined functions,
which can be plotted. Now we use that feature, and construct a graph of a flow of
liquid in our hydraulic system (see Fig. 17.12) as a function of a pump head. For
this purpose, it is enough to make in the calculation, shown in this figure, small
changes, displayed in Fig. 17.14.
First, remove from the calculation the initial approximation operator for the
solution for the pump head. This operator is crossed out in Fig. 17.15. Second,
function Find should be assigned not to a vector, as shown in Fig. 17.15, but to a
user’s function named q. This function is not simple, but a vector-function. By
selecting its first element (in this calculation the system variable ORIGIN is equal
to one), one can plot the flow of fluid through the pump as function of its pressure.
Gaps in the line indicate that in these areas Mathcad could not solve a system of
nonlinear algebraic equations.
The Mathcad has the tools for solving not only the systems of nonlinear alge-
braic equations (see Figs. 17.7 and 17.12), but also for solving systems of linear
equations.
Our problem of the pipeline network has a direct equivalent in electrical engi-
neering: the pump is a source of electrical current (generator, battery, etc.) with a
certain electromotive force (EMF), expressed in volts, and the sections of pipelines,
valves, tees, angles, etc. appear to be electrical resistance, who’s value (ohms),
unlike the “hydraulic resistance,” depends little on the current. Therefore, solving a
typical problem of electrical engineering—finding the current values in some parts
of the circuit depending on the specified values of local resistances and EMF is
reduced to solving a system of linear algebraic equations. Figure 17.15 represents a
“remake” of our hydraulic problem, as shown in Fig. 17.12, into electrical problem.
Figure 17.15 shows how an electrical problem is solved with Solver block: enter
the initial approximations into the solution: one ampere is assigned to all the six
desired current values. But it is possible to assign other values. The selection of
these values is an important stage in the solution. Mathcad sometimes may not be
able to solve the system and asks to change the initial approximations. Engineer
solving on a computer his problem is to know about what decision he needs, and
provide a first approximation near this point. Next, in calculation in Fig. 17.15 they
record current balance equations (first Kirchhoff’s law), and the balance of “energy”
for individual closed branches of electrical network (second Kirchhoff’s law).
Function Find returned a value of six currents providing an approximate solution of
the original system. We have a problem with eight equations and six unknowns.
234 V. Ochkov
Fig. 17.14 Creating a characteristic curve of fluid flow through the pump versus its head
1
Mathematician’s comment. Not every problem in which the number of equations is greater than
the number of unknowns, is overdetermined. For example, if the equations are linearly dependent,
the described system can be defined (in fact this is shown in an example of calculation), and
uncertain. The ratio of the number of equations and the number of unknowns is not a criterion of
overdetermination. Probably, you should not generally pay as much attention to this delicate issue.
It might be better to say that anything can happen with the systems, but the package itself can (or
cannot!) figure it out.
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 235
2
Mathematician’s comment. If the problem is overdetermined, then why add another equation? Or
why they cannot be added if they are physically necessary? Something should be justified.
Probably, to refer to the theory of calculating circuits, or to something else. In any case, before
considering the numerical solution of the system one should substantiate the adequacy of a chosen
mathematical model (having the form of a system of linear algebraic equations) of a studied
physical model of electrical circuit.
236 V. Ochkov
Fig. 17.16 The solution of an overdetermined system of linear algebraic equations in Mathcad
Prime environment
3
Mathematician’s comment. Why is that? Is it always? Which system? Who or what is required to
specify them in this case? We need to know the specifics of the numerical method, hardwired into
the package.
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 237
Fig. 17.17 Analysis and solution of the “classical” system of linear algebraic equations in
Mathcad prime environment involving the basics of linear algebra
of the unknowns, and the unknowns and the right-side parts of the equations. Then
one can create a matrix of coefficients of the unknowns M and an absolute term
vector v. Function lsolve returned us a solution which is a repetition of the solution
in Fig. 17.15, but without any “initial approximations” that confused us.
From our system, we can remove the two extra equations,4 as we have noted
above. After that matrix M will become square, and it will be possible to calculate
the determinant and invert it (find the inverse matrix). The inverse matrix can be
multiplied by vector of free members and … we get a solution. That’s just where the
knowledge of linear algebra, which we acquired at the University, would be useful
for us! In Fig. 17.17 is solved a “classical” system of six linear algebraic equations
with six unknowns using a “classical” method—vector multiplication of inverted
(reverse) matrix of coefficients of the unknowns by free member vectors. The
function of rank that returns the rank of a square matrix M and an expanded matrix
4
Mathematician’s comment. Why are they superfluous, why these are superfluous and not others?
From the point of view of the theory of systems of linear algebraic equations, the extra equations
are found, and not in the only way, only after finding the rank of the augmented matrix. Here it is
probably more appropriate, without going into the mathematical details, to provide physical or
technical justification for allocating extra equations.
238 V. Ochkov
5
Mathematician’s comment-question. And you cannot get it? How is it determined?
6
Mathematician’s comment-question. What does it mean “similar”? If the solution is unique and
approximative (numerically) coincides with the desired, then for the technical calculation this is
not a nuisance, but the norm. If the mathematical theory proves the absence of a solution, which is
yet present (i.e. calculated), or the existence of an infinite number of solutions (unknowns are more
in number than the rank of the system), but only one is calculated, it means that the used model is
incorrect. And this is a very big trouble.
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 239
in solving such problems. But if you do not aim at reducing the number of equa-
tions, the solution of overdetermined algebraic system could be simply charged to a
computer.
Analysis of linear equation systems in our “non-linear” Study is not accidental.
The original hydraulic system is non-linear, but can be temporarily linearized, i.e.
take the rough assumption that the resistance of the pipeline does not depend on the
fluid flow rate. After that it will be possible to solve such a linear equation system
using the function lsolve (see Fig. 17.16) or multiplying the inverse matrix of
coefficients of the unknowns by a free members vector (see Fig. 17.17). For
completeness, here you can also add a Gauss method, that, by the way, in theory
easily and unambiguously copes with any number of equations and unknowns. The
resulting solution is then to be used as a first approximation for the function Find,
applied to a nonlinear original problem of a network of pipelines (see Fig. 17.12).
By the way, our problem on the electrical circuit (see Figs. 17.15, 17.16 and 17.17)
is also nonlinear. With the passage of electrical current the resistance gets heated.
And the value of the electric resistance, same as the value of thermal resistance see
Chap. 20), strongly depends on the temperature. If electric current is small (with a
slight heat resistance), these changes can be ignored, but in case of thermal elec-
trical installations one cannot do this. One can use a ohmmeter to measure the
resistance of an incandescent bulb, divide the square of voltage of the grid by the
measured resistance of the tungsten filament, and … do not get the value of the
bulb’s wattage, stamped on it. The electrical resistance of tungsten greatly increases
when it is heated. Otherwise, a tungsten filament of a light bulb would immediately
burn out when it is turned on. This is an interesting problem at the interface of heat
and electrical engineering: to calculate variation in time of a bulb filament tem-
perature when the bulb is switched on or off, as well as at power surges. This
problem belongs to the class of problems for differential equations, which we shall
consider in next chapter. An analogue of this non-standard heat exchanger problem
is shown in Fig. 18.10 in Chap. 18.
Let’s go back to the second title of the chapter. Immanuel Kant is credited with
the saying: “Every natural science contains as much truth as much mathematics it
contains.” In the Moscow Power Engineering Institute, where the author teaches a
course “Information Technology”, freshmen in parallel are taught a course of higher
mathematics, containing sections related to linear algebra (see above) and calculus
(mathematical analysis). In subsequent semesters they are also given a course of
numerical methods in mathematics. Without this triad it is impossible to proceed
with thermal engineering—thermodynamics and heat-mass exchange.
Foundations of linear algebra, more precisely, matrix calculus, with the original
and inverted matrices, matrix ranks, etc. were illustrated here by solving a problem,
as shown in Fig. 17.17. Function Find (see Figs. 17.12, 17.13 and 17.14) includes
various methods of a numerical solution of systems of linear and non-linear alge-
braic equations. (One could see the list, in Mathcad 15 environment right click on
the function Find.) Now we will look at the basics of calculus again using “heat
exchange” problems. However, the first of the problems discussed further is not a
calculus, since it does not contain a differential or integral calculus, any rows or
240 V. Ochkov
Fig. 17.19 Graphic analytic solution of underdetermined system of nonlinear algebraic equations
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 241
7
The problem of the dynamics of motion of a body in a viscous medium, we consider in the next
chapter.
8
Engineers have to solve not purely technical but technical and economic problems.
242 V. Ochkov
9
Present-day continuation of this “joke with the beard” is below. People on the balloon shout in
response that the batteries got low in their navigator. The man on the ground looks at his
smartphone with navigation and shouts back, “You are at the point of 0.9677 rad north latitude and
0.646 rad east longitude. Only mathematicians measure angles in radians, not degrees.
10
Caution! The computer can be pushed to the wrong decision.
17 Cogeneration (CHP), Trigeneration (CCHP) and Quadrogeneration … 243
Fig. 17.21 Analytical solution of the problem of the optimal size of the cylinder
the objective function with two unknowns r and h (see Fig. 17.20). However, you
can think a little with your computer and figure out that this problem has one
unknown quantity—the radius of the cylinder, for example. The minimum of a one
argument function can be found in the diagram (Fig. 17.22). But we can limit
ourselves to a solution with two unknown quantities11 (see Fig. 17.18) using the
KISS-principle: keep it simple, stupid! The main thing here is a skillful combination
of analytical and numerical methods for solving problems and … participation of
mathematicians in their critical analysis. Let me give a concrete example. Author of
the book sent this chapter to one good mathematician, who corrected the text and
11
In the transportation problem shown in Figs. 7.4, 7.5 and 7.6, we were guided by the
KISS-principle and introduced four unknowns, although it was possible to have only one
unknown.
244 V. Ochkov
Fig. 17.22 The combination of analytical and graphical methods in solving the problem of the
optimum size of the cylinder
a power plant may be as follows. Large cities in Central and Northern Russia, the
Urals and Siberia in the winter have significant problems with the disposal of snow
removed from the roads. This snow is not simple water, but is highly contaminated
with icing reagents, which are sprayed over the road after each snowfall. Our
municipal CHP in the winter period could store the snow on its territory and in
summer melt the snow separating icing reagents from it for their subsequent reuse
in the coming winter. The purified water can be used to make-up the power units
and cooling systems of turbine condensers. Efficiency of such a power plant, by the
way, would rise by reducing the temperature of cooling water flowing to the
condenser [66, 67]. Power plants have chemical shops where salt and other
impurities are removed from the natural water. They can be used to clean “last
year’s snow” melting water.
In many northern countries (Scandinavia, Canada), they never spray icing
reagents over the fallen snow on the country roads, the snow is simply raked from
the road or tamped straight on the road and sometimes dusted with granite chips,
which in the spring is collected, cleaned of the mud and used again. But the use of
granite chips on city streets is problematic, also due to the fact that they will clog
the sewers where the snow is dumped from snow melting plants, for example, or in
a natural way in spring or in a thaw.
For snow-melting in Moscow they started to apply not only sewage wastewater,
but also warm water from CHP heated in the turbine condensers with
straight-through cooling systems where water is taken directly from the river and
discharged back into the river. With this in mind, it should be recognized that the
scheme of dumping snow into CHP drains is similar to the old rejected scheme of
direct dumping of snow into the rivers. One “cosmetic” difference—the platforms
which are frightening the environmentalists and used to dump the snow from the
trucks are no longer built on embankments—now the snow is dumped into a river
remotely (by the way, this method have been used all the time at CHP when dealing
with snow on the territory of the power plant). But the CHP in Moscow may be
more actively and in a more civilized way involved in solving the problem of snow
in the city and addressing some of their own problems.
We ourselves need to clearly understand, that the snow containing various icing
regents cannot be discharged into the sewerage system. This snow, we repeat,
should be segregated from the regular snow when collected, then it should be
processed and separated from the icing reagents which may be reused, if the city
must spray roads with icing reagents. This should be done, not only considering the
sewage farm’s environment, but simply by a common cultural principle, “it is
necessary to clean up after yourself!”. Moscow CHP can be very useful, especially
those that are located closer to the center of Moscow.
Imagine such an ideal scheme for processing contaminated snow. Snow, more
precisely, the brown muddy mush, taken from the Moscow roads sprayed with icing
reagents, is brought into special bins at Moscow CHP sites. Maybe it is not even
worth melting the snow in winter, especially spending the fuel, and better to wait till
warm weather. In summer, the snow will be melting by itself and also due to the
CHP waste heat, increasing, by the way, the efficiency of the power plant due to the
246 V. Ochkov
Valery Ochkov
Abstract This chapter describes the common types of differential equations which
you will work when solving the problems of thermal engineering tasks. Described
the numerical methods for solving such problems. Continued analysis of finite
element methods applied to heat conduction problems.
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
and mass transfer.1 Water, by the way of spraying, may be enriched with oxygen.
Such spray systems are sometimes installed in dying reservoirs for water aeration,
cooling, and, consequently, slowing certain “harmful” biological processes therein.
The famous 140-m fountain on the Geneva Lake in Switzerland (see Fig. 18.7) is
not only a tourist attraction and advertising for the company that supplied the pump
for the fountain, but also a device for feeding the lake water with oxygen and
cooling it to some extent.
Spray cooling ponds are also widely used in industrial air-conditioning and
refrigeration systems (see Chap. 19).
Let’s fit the pipe end shown in Fig. 16.6 with several nozzles at different angles
to the horizon from which water will flow at the speed of v0, and see (calculate)
how water will flow out breaking into droplets. We’ll simplify the task—we’ll
simulate the flight of a single spherical drop of water under the force of inertia,
gravitational forces, and air resistance. Rather, we should say as follows: no force of
inertia exists, but there is the Newton’s law, according to which the sum of forces
acting on the body (on our droplet)—friction force and gravitational force—is the
product of body mass by its acceleration. In daily use, this product is sometimes
referred to as the force of inertia.
The problem will feature two pairs of required functions depending on time t:
droplet position in the space, or rather, in the xy plane (x axis is a horizontal line,
and y axis is a vertical: we solve the planar rather than spatial problem) and droplet
speed v using these two coordinates.
Speed is the first path derivative with respect to time, and acceleration is the
second path derivative with respect to time, or the first speed derivative with respect
1
There are also so-called dry cooling towers (Heller cooling towers) in which water and air are
segregated by a metal (usually aluminum) wall. Water cooling in such cooling towers is not as
inefficient, but there is no loss of cooling water due to evaporation and droplet entrainment. In
“wet” cooling towers and spray cooling ponds salt concentration in water is increased due to
evaporation, which might lead to deposition of scale in the tubes of the steam turbine condenser,
increase in steam condensation temperature and reduction in efficiency of the power unit. The
network calculation placed at http://twt.mpei.ac.ru/MCS/Worksheets/Water/acidation.xmcd,
allows us to calculate measures aimed at prevention of scale formation in the turbine condensers.
250 V. Ochkov
2
Temperature and pressure values used to determine the values of such thermal and physical
parameters, can be entered rather than values of water and air densities and air kinematic viscosity
value. We did so in other studies of the book, based on the “cloud” features and templates.
3
Wind affects spray cooling ponds operation. The above mentioned fountain on the Geneva Lake
is automatically deactivated at high wind speed, elsewise water from the fountain would irrigate
the lakeshore.
4
The wind blows from right to left. But the wind speed component could be taken into account in
the vertical direction. It takes place in the above mentioned cooling towers with the vertical air
draught. We can also switch to the spatial task, i.e. enter the wind speed component in the direction
perpendicular to x and y axes.
18 Differential Equations in Thermal Engineering 251
The technology of creating such animation is described in Chap. 12. The animation
itself can be viewed here: https://www.ptcusercommunity.com/videos/4530.
The system of equations shown in Fig. 18.2 should be, of course, much more
complicated. It should, for example, more accurately take into account the velocity
components along horizontal and vertical directions and absolute value of water
droplet velocity, as is done in the solution shown in Fig. 18.6. Such nuances are
discussed at forum https://www.ptcusercommunity.com/message/228797. But it is
possible to get confined to a simplified solution shown in Fig. 18.2, taking into
account the fact that our spherical water droplet in the actual flight will deform and
break up into smaller droplets. Our simplified solution gives a quite believable
animation. And this is sufficient for initial (quality) problem statement. The problem
can be further simplified: assume, for example, the air friction coefficient of the
droplet equal to 0.5—see the first formula and graph in Fig. 18.2.
The correctness of our model of the droplet flight can be verified qualitatively
rather than quantitatively: for example, compare two images—the flight curve of a
virtual droplet and type of real fountain of Geneva, mentioned above (Fig. 18.7).
Figures 18.3, 18.4, 18.5 and 18.6 showed statement and numerical solution of
the initial problem for the system or problem with the initial conditions of so-called
ordinary differential equations—see letter o in the title of the odesolve function. The
solution of such system is a function of one argument.
But in Mathcad 15 environment, tools are available that generate a function of
two arguments while solving differential equations. One such problem was just at the
author’s hand who wrote and corrected the text the tablet PC, the back wall of which
was warm in several places. This resulted in cooling down the tablet chips without
using noisy and energy consuming fans built into laptops and desktops.5 Figure 18.8
shows the final part of the Mathcad-document that calculates the temperature field of
three chips on the computer board with two dimensions of x and y.
The problem displayed in Fig. 18.8, uses the relax function built into Mathcad,
which returns the temperature values matrix (chip temperature in our problem)
depending on the input data—the chip coordinates, their specific heat evolution and
heat transfer conditions (arguments of the relax function). This problem is well
described in detail in the electronic manual in terms of heat and mass transfer by
Solodov [14]—http://twt.mpei.ac.ru/solodov. On this site, the reader will find other
useful materials related to the topic of this book, use of Mathcad and Internet for
thermotechnical calculations.
Based on the function of two variables rather than one variable, two partial
derivatives can be assumed—for the first argument and for the second argument.
The pdesolve function built into Mathcad 15 (p—partial) is designed for numerical
5
Elements of supercomputers are sometimes provided with water cooling. Water in such systems
can also be cooled down in spray cooling ponds, which were simulated in this chapter.
256 V. Ochkov
not quite true—the second derivative is based on the second measurement rather
than on the time), the finite-difference approximation in the grid nodes is shown
into which a flat plate is divided [60]. Using such approximation, one can make a
system of difference equations and solve it by the iterative method.
In Fig. 18.11, the Mathcad Prime 3 software is shown to solve this heat transfer
problem: given is a rectangular plate with the pre-set temperature T1, T2, T3, and
T4 in four angles of the plate. Along the plate edges, the temperature varies linearly
from one angle to another.6 Find temperature distribution in the plate—generate the
function T(x, y).
In Fig. 18.11, in addition to parameters of the plate itself, the following values
have been further defined: nx and ny—the number of grid divisions on the plate in
the horizontal (x) and vertical (y) directions, disparity of the solution (εT) and
maximum number of iterations (n). In the program shown in Fig. 18.11, operators
are intended for the following:
• setting the Temperature(T0, T1, x) function for linear temperature distribution
between two points;
• temperature calculation at the grid nodes at the plate edges: Ti,1, Ti,ny+1, T1,j и
Tnx+1,j; thus, the remaining (internal) matrix elements become zero;
• calculation of the horizontal (hx) and vertical (hy) grid spacing, as well as
auxiliary value of R7; moreover, the T matrix is copied into the TT matrix; the
auxiliary TT matrix will store temperature values at the grid nodes at the pre-
vious calculation; the first calculation does not contain the previous calculation,
6
These are the boundary conditions of the first kind. Under boundary conditions of the second
kind, heat flows at the body boundaries, and under boundary conditions of the third kind—heat
exchange conditions (heat exchange coefficient values).
7
This value should not be entered into calculation, and directly used in the iterative formula to
make it clearer to understand. But from the cycle, only that can be taken out which can be
calculated outside the loop. This speeds up calculation. And it might take us a lot of time, because
it involves nested loops with a lot of iterations.
18 Differential Equations in Thermal Engineering 259
thus, this matrix is introduced artificially here (these four operators are written in
one line in the form of a matrix with one line and four columns (undocumented
technique in the Mathcad environment) in the program for its compactness);
• in the loop with the s parameter, the IsεT feature is introduced, which is assigned
a single value; if this feature remains equal to unity (see below), this means that
solution has been found in the next approximation, since the temperature value
in each grid node is approximately equal to such value in the previous
approximation, rather, it differs in terms of modulus by no more than the pre-
determined εT value;
• further, in the loop with s parameter, approximations are made to solving the
problem: two nested loops with i and j parameters fill the T matrix—further
approximation to the solution, and then are verified if the value of the next
element of such matrix is close to the corresponding value of the TT matrix
260 V. Ochkov
Difeqs are time-tested methods and tools of research of surrounding real physical
processes. For example, body cooling or heating (see above), current passage
through electrical circuit or radioactive decay, body falling down (of our drop of
water) under the action of gravity, pendulum oscillation, etc.
Let’s talk about the recent past of students, postgraduate students, and engineers—
about school years. After all, it was the first time when familiarization with the
elements of differential equations theory begins. But this fact is persistently hidden
from students. Teachers of physics believe that it is more useful for a students to learn
the calculation formula than to understand where this formula is derived from.
Teachers of mathematics believe that it is enough to make a student learn how to define
a derivative using formal rules and tables than justify the need for study of this
mathematical concept.
Earlier, there were excuses—inaccessibility of computer tools intended for
solution of difeqs. Now, such excuses seem to be unsound. Almost all physical and
mathematical phenomena studied at school or university can be easily simulated
using differential equations solved by means of computer software with good
graphics and animation.
Once one of the authors of this book gave classes on informatics in a very
progressive Moscow Lyceum (high shcool). In terms of lessons, solution methods
were reviewed in the environment of mathematical software Mathcad for differ-
ential equations. The director of the Lyceum (incidentally, a physicist by education,
and part-time—Professor of the Physics Platform at the prestigious Moscow uni-
versity) said that the words “differential equations” namely “difeqs” are imposed a
strict secret ban at high school [68]. The students should at least understand simple
algebraic equations, but here they have to deal with the differential ones… But
when this director was shown what equations will be addressed at the lessons and
how they will be solved on the computer, the director changed his mind, and
expressed confidence that it will be interesting, understandable, and, what is most
important, useful for the students.
The Mathcad solves differential equations numerically, but not analytically. This
means that we obtain a set of numbers—table of arguments and values of the
desired function, rather than the formula-solution. The desired function can be
represented graphically and even animated. By changing the initial conditions, it is
possible to follow how the graphics will change. For example, solve the same
problem, but with a different initial condition. Moreover, such initial conditions can
be empirically selected which provide the desired behavior of the solution.
The main thing you need to understand that notorious difeqs are necessary and
useful simulation tool for a lot of physical and technical processes, rather than a
school or university “horror story” and “torture instrument”. That’s what we tried to
show in this chapter.
Chapter 19
Refrigeration Cycles
Volodymyr Voloshchuk
V. Voloshchuk (&)
National University of Water Management and Natural Resources, Rivne, Ukraine
e-mail: [email protected]
temperature is higher than that of the high temperature surroundings, and heat is
rejected naturally from the condensing refrigerant to the surroundings. The refrig-
erant liquid then passes to an expansion valve where the flow is throttled from high
to low pressure. The cycle is closed.
The efficiency with which a refrigerator or heat pump utilizes electricity to
produce either cooling or heating is known as the coefficient of performance. The
cooling coefficient of performance, COPR , is the cooling effect at the evaporator
divided by the work done at the compressor. The heating coefficient of perfor-
mance, COPHP , is the heating effect at the condenser divided by the work done at
the compressor [69–71]:
qL
COPR ¼ :
lC
qH
COPHP ¼ :
lC
www.wsp.ru and installing it on a computer,1 more than 300 functions not only on
the properties of working fluids, but also on some processes in thermodynamic
cycles become available for thermal engineering calculations.
Downloading functions from the Internet (or installing them on a computer from
a disk or other carrier if the workstation does not have access to the Internet) has
one essential drawback.
Computer programs, in particular, those for calculating the thermophysical
properties of individual substances and their mixtures are constantly upgraded and
improved. This is mainly due to the fact that new formulations (sets of formulas
with their description)2 emerge, that determine the procedures used to calculate the
concrete properties of certain substances.3 In addition, errors and inaccuracies in the
existing computer programs are corrected, their application domains are extended,
and their performance is improved (speed of operation, volume of computer
memory, etc.). Such programs are also continuously updated due to changes made
in the hardware and operating system software. Users of computer programs on the
properties of substances frequently fail to notice these changes and work with
outdated versions. Also, users face additional difficulties if they change their
computer and/or operating system: old computer programs cannot be installed and
started any longer.
In view of what was said above, and taking into account that at present almost all
computers on which engineering calculations are carried out (in particular, thermal
engineering ones) have constant highspeed connection to the Internet, we have
suggested a new technology for working with functions on the thermophysical
properties of working fluids, heat carriers/coolants, and materials used in thermal
power engineering that is based not on downloading computer programs but on
making reference to the functions stored on websites, or in “clouds”.
Figure 19.2 shows a part of the “cloud” calculation server of the Moscow Power
Engineering University’s (MPEI) National Research University (www.mpei.ru) and
Ltd Trieru (www.trie.ru) containing a list of links to a wide range of web-pages
with functions for calculating thermophysical properties of working substances
(including refrigerants). The links are located in the interactive handbook “Power &
1
Almost all software can be used for carrying out these calculations, such as the Excel spreadsheet
processor, Mathcad engineering calculator, Maple symbolic computation language, high level
programming languages C, BASIC, Pascal, and FORTRAN. In this paper, we limited ourselves to
considering only the Mathcad.
2
As regards water and steam (including solutions, e.g. seawater), which are used as the main
working fluids in power engineering, such formulations are developed and approved by the
International Association for the Properties of Water and Steam, IAPWS (see the website
www.iapws.org).
3
Here is a concrete example. In September 2011, the IAPWS approved a new formulation for
calculating the thermal conductivity of water and/or steam. Users who downloaded the computer
programs earlier will continue to work with their old versions for a long period of time, whereas
those who employ the new information technology described in this paper will immediately begin
to work with the software implementation of the new formulation on the thermal conductivity of
water and/or steam.
266 V. Voloshchuk
Fig. 19.3 The web-page on “cloud” functions of thermophysical properties of refrigerant R407c
with options of online computing, download and reference
Fig. 19.4 A copy of the internet-page for online calculation of specific enthalpy as a function of
pressure and temperature. Illustration of the determined point in h, p- and h, T-phase diagrams
with the link on-line calculation (see Fig. 19.3) to the calculation page shown on the
Fig. 19.4.
The user is able to change initial data (in this case pressure and temperature) and
to choose required dimensions. After clicking the button Recalculate he obtains the
results and graphical illustrations of the calculated point in h, p- and h, T-phase
diagrams.
The mathematical software Mathcad (the web-page, shown in Fig. 19.4, was
created with it—see Chap. 11) enables us to solve a wide range of scientific,
engineering and technical problems without using traditional programming.
Formulas are written in the Mathcad in the language of mathematics, which is
naturally readable and does not demand deep knowledge of programming language.
Mathcad gives us the possibility to see results clarified by plots and diagrams—a
feature that favorably distinguishes it from traditional programming languages and
spreadsheets. The computer application Mathcad also gives the possibility to use
units for checking the correctness of calculations and for presenting results in a
more convenient way, especially in engineering calculations. Mathcad has an
option to exchange computational applications with other computer programs (such
as MathSoft’s Axum and S-PLUS, Microsoft Excel, MATLAB). All that makes
Mathcad calculations more efficient.
The Mathcad has a convenient tool: references to another Mathcad document the
variables and functions of which become accessible or, as programmers say, visible
in the Mathcad document from which the corresponding reference is made. The
Mathcad user need not open another calculation document and insert it in his or her
268 V. Voloshchuk
engineers or scientists need to spend additional time for manual calculation and
input necessary properties into calculation or for search and integration into cal-
culation function on thermophysical properties of refrigerants. The resources of the
website www.vpu.ru/mas are proposed to help engineer to solve his working tasks.
If to analyze Fig. 19.3 a link R407cPSVT(T) can be observed, clicking on
which Mathcad-function titled as R407cPSVT(T) can be downloaded. The func-
tion R407cPSVT(T) determines pressure of vapor R407c on the saturation line as a
function of temperature. It is possible not to download and place this function into
user’s worksheet but just to do Internet-reference to the function.
If to move the pointer of the mouse to the link R407cPSVT(T), see Fig. 19.3,
and to click its right button a dialog window will appear (see Fig. 19.6) where the
corresponding URL address of the necessary “cloud” function can be seen and
copied.
To make this function activated in calculations of refrigerant cycles it is nec-
essary to make a reference to it in the MathCad worksheet. This procedure is shown
in Fig. 19.7: in Mathcad choose the menu items Insert then Reference and insert
the previously copied address shown in Fig. 19.6 in the appeared window. After
such procedure (see Fig. 19.5) it is possible to calculate the necessary pressures pev
and pcd of refrigerant R407c and to make the result visible: pev: = R407cPSVT
(t′′ev) = … (comment: pressure is shown in pascals but a user can fix other units, for
example megapascals (MPa) as demonstrated in Fig. 19.5). Pressures pev and pcd
are calculated at the specified temperatures on base of function relationships on
Fig. 19.8 Ts- and ph-diagrams for a simple vapor-compression refrigeration cycle
Fig. 19.9 Input data block and schematic of simple vapor-compression refrigeration cycle
prepared on technology of Mathcad calculation server
The web-page on this calculation can be found on the server in the interactive
reference handbook “Thermodynamic cycles”, Chap. “Heat pumps and refrigerators”.
As input data it is proposed: temperature of heat source fluid at evaporator outlet,
temperature of heat consumption fluid at condenser outlet, temperature difference at
cold end of the evaporator, temperature difference at warm end of the condenser,
isentropic efficiency of the compressor, refrigerant pressure losses in the condenser,
refrigerant pressure losses in the evaporator and a type of refrigerant. For determining
of thermodynamic properties of refrigerant references on “cloud” functions were used.
Ts-, ph- and hs property diagrams for the presented simple vapor-compression
cycle of a heat pump, plotted in Mathcad Calculation Server are shown in Fig. 19.10.
274 V. Voloshchuk
Fig. 19.10 Ts-, ph- and hs property diagrams for a simple vapor-compression refrigeration cycle,
plotted in Mathcad calculation server
Fig. 19.11 Results of on-line calculations on determining influence of temperatures of the source
and heat consumption on COP of a heat pump cycle: dashed lines—actual cycle; solid lines—
ideal cycle; working fluid—refrigerant R407c
19 Refrigeration Cycles 275
In the Fig. 19.11 results of on-line calculations made with the help of web-page
http://twt.mpei.ac.ru/MCS/Worksheets/PTU/Vv-27-eng.xmcd (see interactive ref-
erence handbook “Thermodynamic cycles”, Chap. “Heat pumps and refrigerators”,
“Thermodynamic calculations of Simple Vapor-Compression Refrigeration Cycle”)
are introduced. The results demonstrate influence of the temperature of low
potential heat source and outlet heating temperature on the coefficient of perfor-
mance of simple vapor-compression ideal and actual cycles of a heat pump. R407c
is chosen as a refrigerant.
Such dependences are widely spread in literature. But with the help of the
introduced recourses it is possible to investigate a quite wide variety of input
conditions (temperature of low potential heat source, outlet heating temperature,
compressor efficiency, pressure losses etc.) and, for example, to determine
numerically the most irreversibilities taking place in a heat pump cycles and so on.
In actual refrigeration cycles irreversible processes take place which decreases
their efficiency. Energy losses for friction due to moving of working fluid, heat
losses, pressure drops, temperature differences can be the reasons of these
irreversibilities.
In the calculations of the actual cycle contrary to ideal one 3 % pressure drops in
evaporator and condenser, minimum temperature differences in evaporator and
condenser Δt = 2 °C and 82 % isentropic efficiency of compression is taken into
account.
The Fig. 19.11 shows that for the chosen refrigerant, outlet heating temperature
t = 50 °C and increasing temperature of low potential heat source from −10 to 10 °C
coefficient of performance of ideal cycle of a heat pump increases from 3.7 to 5.6.
For actual cycle in the same temperature conditions COP of a heat pump increases
from 2.9 to 4.3 which is by 22…26 % lower in comparison with ideal cycle. Besides
COP increasing can be observed when temperature of heat consumption fluid is
decreased. For example, for temperature of low potential heat source t = 0 °C and
decreasing temperature of heat consumption fluid from +70 to +50 °C coefficient of
performance of actual cycle of a heat pump increases from 2.4 to 3.5. It should be
noted that compressor efficiencies are a strong function of pressure ratio. In the
introduced calculations efficiency of compression was set as constant for different
temperatures of the source and heat consumer and as a consequence for different
pressure ratios. In the proposed web-resource it is possible to change efficiency of
compression.
Figure 19.12 presents results of on-line calculations made with the help of the
same web-page http://twt.mpei.ac.ru/MCS/Worksheets/PTU/Vv-27-eng.xmcd for
the purpose of determining influence of working fluid type on COP at different
temperatures of heat source and outlet heating temperatures.
As it can be seen that for temperature changes of low potential heat source from
–10 to 10 °C and heating consumer from 50 to 70 °C replacing refrigerant R407c on
R134a results increasing COP of a heat pump by 0.3…0.6 units (Fig. 19.12).
There are a variety of ways that the refrigeration cycle can be tailored to suit an
application in a better way (not always necessarily resulting in a higher COP) than
the simple (basic) vapor compression cycle.
276 V. Voloshchuk
Fig. 19.12 Results of on-line calculations on determining influence of temperatures of the source
and heat consumption on COP of an actual heat pump cycle dashed lines—refrigerant R407c;
solid lines—refrigerant R134a
Consider the addition of a single heat exchanger to the basic vapor compression
cycle, exchanging heat between the fluid leaving the evaporator and the fluid
leaving the condenser, shown in the schematic of Fig. 19.13.
Figure 19.14 shows a heat pump cycle in Ts-, ph- и hs- property diagrams
plotted in Mathcad Calculation Server in accordance with input data introduced in
the Fig. 19.13. This web-page is located at http://twt.mpei.ac.ru/MCS/Worksheets/
PTU/Vv-28.xmcd (see interactive reference handbook “Thermodynamic cycles”,
Chap. “Heat pumps and refrigerators”, “Thermodynamic calculations of Simple
Vapor-Compression Refrigeration Cycle with Regeneration”). From ph- property
diagram it can be observed that due to addition of a regenerator (heat exchanger) the
specific enthalpy (and temperature) of the preexpansion state are reduced and
conversely the specific enthalpy (and temperature) of the pre-compression state are
increased. These differences point out two intended benefits of this cycle
modification.
First, since the specific enthalpy remains constant during expansion, a reduction
of the specific enthalpy prior to expansion results in a reduction of specific enthalpy
prior to evaporation. Therefore the unit will have more evaporative heat transfer to
provide more evaporator cooling capacity. Second, the state prior to compression is
further away from the saturated vapor line. For most compressors, it is imperative
that the state of the refrigerant prior to compression does not have any liquid in the
form of droplets or mist, since liquid entrained in a vapor undergoing compression
tends to damage the fast moving parts of a compressor, seriously degrading the
performance and working life span of the compressor. For this reason, it is usually
desirable for the refrigerant to enter the compressor as a superheated vapor, several
19 Refrigeration Cycles 277
Fig. 19.13 Input data block and schematic of vapor-compression refrigeration cycle with
regeneration prepared on technology of Mathcad calculation server
278 V. Voloshchuk
Fig. 19.14 Ts-, ph- and hs- property diagrams for a vapor-compression refrigeration cycle with
internal heat exchange, plotted in Mathcad calculation server
Fig. 19.15 Input data block and schematic of two-stage vapor compression refrigeration cycle
with flash chamber and separate vapor mixing intercooler prepared on technology of Mathcad
calculation server
hs-diagrams which were plotted with the help of interactive calculation in Mathcad
Calculation Server in accordance with initial data presented in Fig. 19.15.
One of the optimization task of thermodynamic analysis of two-stage vapor
compression refrigeration cycle is a task of determining intermediate pressure in
separate vapor mixing intercooler for which coefficient of performance has the
highest value.
For solving this task the same web-page http://twt.mpei.ac.ru/MCS/Worksheets/
PTU/Vv-30-eng.xmcd of the server www.vpu.ru/mas can be used. By changing
pressure in flash chamber pfc it is possible to calculate coefficient of performance
COP and to plot functional relationship COP = f(pfc).
Figure 19.17 illustrates such graphical relationship plotted in MathCad. From this
graph, it can be seen that there is an optimum value of pressure pfc_opt for which
COP is the highest. In the analyzed case according to data presented in Fig. 19.15,
19 Refrigeration Cycles 281
Fig. 19.16 Ts-, ph- and hs property diagrams for two-stage vapor compression refrigeration cycle
with flash chamber and separate vapor mixing intercooler, plotted in Mathcad calculation server
this pressure is equal pfc_opt = 1.5007 MPa. For determining pfc_opt instrument of
tracing is used and demonstrated in Fig. 19.17 which is possible in Mathcad.
Optimal value of intermediate pressure in separate vapor mixing intercooler
pfc_opt can also be determined with the help of other resources of Mathcad. For
example, using programming operators it is possible to create functional relation-
ship between coefficient of performance and pressure pfc. In this case a developer of
this relationship can use the algorithm which is introduced at http://twt.mpei.ac.ru/
MCS/Worksheets/PTU/Vv-30-eng.xmcd. Besides, for calculation of thermody-
namic properties of refrigerants in programming block of Mathcad-worksheet it is
proposed to do reference on “cloud” functions. This procedure is described above.
A fragment of developing functional relationship COP = f(pfc) in
Mathcad-worksheet with the help of programming block is shown in Fig. 19.18.
Having function COP = f(pfc) it is possible to determine its extreme using
Mathcad operators the choice of which depend on peculiarity of function. In the
analyzed case it is proposed to use operator Maximize(COP, pfc)—determining
local maximum. From Fig. 19.18 it can be seen that, on contrary of graphical
solution presented in Fig. 19.17, such numerical solution of the proposed task gives
more exact result: pfc_opt = 1.481 MPa.
In the other flash chamber/regenerative intercooling multistage subtype, there is
no separate mixing chamber.
282 V. Voloshchuk
Fig. 19.17 Graphical relationship COP = f(pfc) in two-stage vapor compression refrigeration
cycle with flash chamber and separate vapor mixing intercooler determined with the help of
web-page http://twt.mpei.ac.ru/MCS/Worksheets/PTU/Vv-30-eng.xmcd and plotted in Mathcad
Fig. 19.19 Input data block and schematic of two-stage vapor compression refrigeration cycle
with combination flash chamber/regenerative intercooler prepared on technology of Mathcad
calculation server
Fig. 19.20 Ts-, ph- and hs property diagrams for two-stage vapor compression refrigeration cycle
with combination flash chamber/regenerative intercooler, plotted in Mathcad calculation server
If to compare diagrams in Figs. 19.16 and 19.20 it can be seen that in cycle with
combination flash chamber/regenerative intercooler the state of vapor at compressor
inlet is saturated, but in cycle with flash chamber and separate vapor mixing
intercooler the state of vapor at compressor inlet is superheated. As a result, in the
first cycle (Fig. 19.16) performance conditions are better and effectiveness of
compression is higher.
Using resources of the server www.vpu.ru/mas located in the interactive refer-
ence handbook “Thermodynamic cycles”, Chap. “Heat pumps and refrigerators”
interactive calculations have been performed for determining coefficients of per-
formance of the four introduced above ideal cycles of heat pumps. The initial data
were the following: temperature of heat source at evaporator outlet t = 2 °C,
temperature of heat consumption fluid at condenser outlet t = 55 °C, temperature
difference at cold end of the evaporator δ = 2 °C, temperature difference at warm
end of the condenser δ = 2 °C, type of the refrigerant—R407c. The results of
calculations are presented in Fig. 19.21. It can be seen that for the same input
conditions addition of a heat exchanger to the basic vapor compression cycle,
exchanging heat between the fluid leaving the evaporator and the fluid leaving the
19 Refrigeration Cycles 285
condenser, results increasing COP from 3.841 to 4.055. Additional changes to the
basic heat pump cycle, e.g. two-stage compression, result increase of COP to 4.359
for cycle with combination flash chamber/regenerative intercooler and to 4.358 for
cycle with flash chamber and separate vapor mixing intercooler.
Chapter 20
Three-Layer Thermal Engineering Cake
or a Conclusion
Valery Ochkov
Abstract In the final chapter of the book described some methods for solving
problems of stationary thermal conductivity. Arguments about the relationship of
mathematics and thermal engineering, on the content of lectures on mathematics for
students are provided.
Figure 20.1—This is not merely a picture, not one of those drawings that just
illustrate the studies. The author has this figure printed in color, framed and hung on
the wall in his office among other “works of art” created with Mathcad. The author
brings his guests-mechanical (or thermal) engineers to the “picture” and asks them
to solve the heat exchange problem hidden in it.
Figure 20.1 is, as follows from its title, a Mathcad-document with the solution of
the problem of stationary thermal conductivity of a flat three-layer infinite wall.
Input data for calculation: the ambient temperature on both sides of the wall (Tb1
and Tb2), the heat transfer coefficients between the environment and the wall at both
sides (α1 and α2), the wall area F, the thickness of the layers of our heat engineering
cake δ1, δ2 and δ3, and also thermal conductivity values of the material of the layers
λ. It is necessary to determine the temperature at the edges of the three-layer wall
(T1 and T4) and at the boundaries of the layers (T2 and T3), and heat flow Q through
the “three-layer cake”. In this cake, sorry, calculation, there are three “flavors”.
The first “flavor”—this is how Mathcad can delicately solve the specific heat
flux balance equations Q in the form of one, so to speak, “five-layer” equation,
where the unknowns T1, T2, T3 and T4 are positioned in places corresponding to
the limits of the definite integrals. And we obtained the integrals because in our
V. Ochkov (&)
Moscow Power Engineering Institute, Moscow, Russia
e-mail: [email protected]
V. Ochkov
Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russia
calculation thermal conductivity of the wall’s material may be not only a constant
but a function of temperature. In the references the temperature dependence of
thermal conductivity of the construction, thermal insulation and other metals is
usually presented in the form of a polynomial of first or second degree. The
argument of this polynomial (a straight line or a parabola) may be any temperature
in Celsius or Kelvin (absolute temperature). In American references one can even
come across Fahrenheit or Rankine (Absolute) units. Temperature, but not as a
relative scale, but as a temperature unit, is also present in a thermal conductivity
unit. This all is very confusing for inexperienced mechanical and thermal engineers
and often leads to mistakes in their calculations.
Figure 20.2 shows an Internet site where the operations related to defining
thermal conductivity of insulating materials are automated.
A visitor to the web-site (see Fig. 20.2) selects from the drop-down list the
necessary material, specifies the temperature in Celsius, clicks Recalculate and
receives the answer: the thermal conductivity value of the selected material at a
given temperature, its density, the maximum temperature at which the material can
be used, and the formula used to calculate the thermal conductivity as function of
the temperature. In particular, thermal conductivity of the mineral wool grade 75
(see Fig. 20.2) is described by a linear function, wherein the temperature should be
Celsius. For other materials the formula may be in the form of a second order
polynomial (a parabola).
20 Three-Layer Thermal Engineering Cake or a Conclusion 289
1
“Stupid” Mathcad when displaying data on the screen simplifies and adds to the value of specific
energy or enthalpy the unit m2/s2 instead of the usual kJ/kg. True: simplicity is worse than robbery.
290 V. Ochkov
in the USA, and the concept of entropy relates with the absolute zero of temper-
ature, which is associated with the kelvin temperature unit in Europe or with
degrees Rankine in America. There is no way to break this duality tradition and no
reason indeed! Everything will be clear to a competent engineer, and the illiterate
may get entangled in any situation.
At the end of Fig. 20.3 one can see the operators translating values of temper-
atures to different scales and in different units. And whereas in the heads of our
mechanical engineers “sits” a constant 273.15, for their American counterparts, it is
459.67.
Surely it is more correct to give empirical dependence of thermal conductivity
and other thermophysical parameters on the temperature in the Kelvin, rather than
in Celsius, as shown in Fig. 20.2. Figure 20.4 shows a reference web-site with data
on the thermal conductivity of pure metals, where the visitor to the site can get
not only the desired thermal conductivity value, but also generate a polynomial of
n-th degree, describing the temperature dependence using kelvins, and not Celsius,
as in Fig. 20.2.
The author’s calculations server has reference pages for thermal conductivity of
various substances and materials, where one can change the temperature differential
range and obtain a simple (linear) dependence for the selected range—see, for
20 Three-Layer Thermal Engineering Cake or a Conclusion 291
instance, Fig. 20.5, where one can see the entire “dependency” or a selected portion
of it “under a magnifying glass.”
Back to Fig. 20.1. A visitor to this site has freedom to choose the units in the
formula: it can be a constant, linear function or a parabola with an argument either
in Celsius, or to an absolute scale.
Working with the temperatures is the second “flavor” of our three-layer thermal
engineering “pie”: on the site, as shown in Fig. 20.1, one can enter the temperature
dependence of thermal conductivity of the materials in the layers as constants or
polynomials.
As shown in Fig. 1.8, Chap. 1, if some materials cover the entire temperature
range it can lead to oscillations—“whipping” of the curve near the points of the
original conductivity table.
The third “flavor” of the three-layer cake, as shown in Fig. 20.1, is as follows.
This figure (paper sheet) is sandwiched between the glass and A4 cardboard frame.
So, if the temperature near the glass is equal to 20 °C, and near the cardboard
backing of the figure is minus 5 °C, the temperature in this “pie” would be dis-
tributed as shown in Fig. 20.1 with the figures and diagram: the thermal conduc-
tivity of glass is given by a polynomial of the second degree, which is taken from
the site shown in Fig. 20.2, and the thermal conductivity of paper and cardboard—
by constants found on the Internet. Heat flow through the picture’s area of 0.06 m2
is 60 W.
292 V. Ochkov
If the cube turns with its one side to the viewer, one could see the T, s-diagram of
a steam turbine cycle, by the second side—h, s-diagram and by the third side—T,
h-diagram. The enlarged view of the lower left corner of T, h-diagram can be seen
in Fig. 5.2 in Chap. 5.
This thermodynamic glass cube is a souvenir which the author gives as a gift to
his colleagues and partners, in particular the users of the program WaterSteamPro
(www.wsp.ru) and “Electronic Encyclopedia of Energy” (www.trie.ru). It is inter-
esting to give this cube to the chief engineers of thermal power plants and energy
companies.2 By the reaction of the people to this gift one can immediately see from
• Engineer Eifel.
20 Three-Layer Thermal Engineering Cake or a Conclusion 295
what “origin” this person has come up to the rank of top manager. Chief engineers
of thermal power plants are promoted from thermal engineers or electricians. So,
the chief engineer who is a thermal engineer immediately appreciates this gift, starts
to turn the cube over in his hands and talk about his student years, when he was
“fiddling” with the calculation of steam turbine cycle when writing a thesis project.
In the eyes of the chief engineer who is an electrician by education one can see at
once some confusion, the question—what is it? Here we have to explain that this is
a diagram of a steam-turbine power unit installed on your power plant, traced on the
thermodynamic surface of water and steam, and for you, electrical engineer, this is
like a hysteresis loop, for example. In order to smooth out some confusion, we have
to tell the old Soviet joke about the two Soviet violinists who were returning home
after an international competition.3 The first violinist took the second place and
grieved about it. The second violinist, which took the last place, was quite happy,
and he says to the first one to stop suffering, and that the second place is a very good
result, there were, in fact, three hundred participants. The first violinist objected to
that in the sense that the winner was given a Stradivarius violin as a gift. The second
“violinist” asks, “What is that?”—”Well, it’s like a Stalin gun for you”—this was
the answer of the first violinist. So, a glass cube with the steam turbine cycle—this
is the same cube, but with a hysteresis loop for an electrical engineer.
In addition to the flat (see Fig. 20.1) and volumetric (see Fig. 20.6) “heat engi-
neering works of art” the author has a more practical souvenir—a wall clock with a
dial shown in Fig. 20.8. On the clock, besides the numbers—typical
Mathcad-operators—there are three Internet addresses, which, as the author hopes,
will be “desktop” addresses for the readers of this book, will be bookmarked in your
favorites of the Internet browser. If the virtual clock is animated with moving arrows
attached to it, the clock can serve as a desktop picture on the computer screen. More
about this clock: https://www.ptcusercommunity.com/message/430540.
Our “three-layer thermal engineering cake with flavors” can be seen as a
“dessert” to this book, as a kind of conclusion. Conclusive remarks.
The author keeps repeating to his students that from the studies they need to try
getting… fun in the first place and only after that knowledge and skills as well.
Chapter is the work, very hard work, but without pleasure, without love for work,
even an unburdensome, prestigious and highly paid job can turn into hard labor.
Pleasure is unthinkable without jokes, without the light and ironic attitude toward
the most serious things. And what could be more “serious” than mathematics, or the
same heat engineering? Our secondary school with its cramming of answers for the
Unified State Exam completely discourages many students from craving knowledge
(Footnote 2 continued)
• This is like our father—an engineer?
• No, son, our dad is chief engineer.
Chief engineer mainly carries out not creative but administrative functions. Therefore he is
correctly renamed Technical Director.
3
Another joke—see the previous footnote.
296 V. Ochkov
and “getting pleasure from the process of learning.” This book can be seen as an
attempt to return to the students interest to this most important of human abilities—
to acquire knowledge. Many students aspire after a computer, but their future
engineering profession is not computer technology (IT), but all the same heat
engineering.
The authors hope that this book will help these students at least through the
computer and the Internet to restore the ability to receive pleasure from studying at
a technical college.
One of the authors of this and the Russian edition of the book Valery Ochkov
reads “Information Technology” course to the freshmen in MPEI (Moscow Power
Engineering Institute). In their second year they have a course in thermodynamics
delivered by Aleksey Alexandrov, another author of the Russian book. We hope
that this book will be a useful and, most importantly, a pleasant bridge between
these academic disciplines. And the basis of this “bridge” is mathematics, the
lectures in mathematics in MPEI are read by Elena Bogomolova, the third author of
the Russian book… “Construction and servicing of the bridge” is maintained by the
other authors of the Russian book.
More than half a century ago, S.M. Ulam, the “father” of the Monte Carlo
method, wrote [72] that each year about two hundred thousand theorems are
published in the mathematical journals. They are added to the already proven ones,
and the amount of mathematical knowledge (yet the same way as knowledge in the
other scientific disciplines) increases more and more. Although mathematics
teachers do not even try to explain in their lectures something mathematically new,
capable to modernize the engineering education, they barely have time even to
teach students the skills and methods which are two centuries old. The result is
20 Three-Layer Thermal Engineering Cake or a Conclusion 297
usually disappointing, the students poorly know the “old” mathematical techniques
and do not own the “new” one. Trying to answer the question “So what’s the point
of teaching mathematics to engineers?”, one inevitably comes to the conclusion that
the amount of mathematical knowledge received now at the Technical University
should be seriously reduced and its nature drastically changed.
Of course, the idea that every engineer has to be a good mathematician continues
to reign with us. Undoubtedly, as it was at the dawn of engineering sciences. Back
in the thirties of the last century, the outstanding scientist Norbert Wiener [73] read
in the University an annual course of combined lectures in mathematics and elec-
trical engineering. Now no one would think to combine these two disciplines into
one. Engineering and technological demands of society have changed dramatically.
The technical means used in solving problems has also changed. Foremost the
engineers, in fact, are consumers, rather than creators, of the finished mathematical
product, which is also “hidden” somewhere in the depths of computer engineering
developments. This fact must determine the current trends in mathematics education
of engineers.
In order to select the new priorities in teaching mathematics in a technical
college and be able to sacrifice a part for the sake of development and prosperity of
the whole, let’s have a look at the basic processes that have occurred in the last
half-century in the field of application of mathematics to engineering sciences. The
main thing we will see is that the object of the application of mathematical
knowledge has changed significantly.
The technologically weak computer of the twentieth century required that future
engineers should have good skills in algebraic transformations, differentiation and
integration, to facilitate the calculations carried out by means of a slide rule and old
math functions tables, and later with the help of calculators, slow IBM-360 and
low-power PC-286. In that time the development of all the world mathematics was
aimed at finding analytical solutions and creating suitable numerical methods for
their most accurate implementation.
Now with the advent of high-performance computing tools and environments the
analytical engineering research has moved into the background. And it can be
clearly seen from the examples in this book. For successful study of heat engi-
neering the current students should have good command of mathematical tools
relying on the potential of computing programs such as Mathcad features.
Many numerical methods of mathematics are implemented in the Mathcad, and
their list is constantly updated. The methods themselves get rid of the programming
errors that are inevitable in the beginning of work. For example, one of the first
versions of Mathcad plotted hyperbole as the letter N, which surprised many
students. Now diagrams created in Mathcad are flawless.
Referring students to software, it is important to consider that behind the ease of
use of finished programs there may be big problems hiding related to borders and
the peculiarities of their application. And if earlier each researcher, as a rule,
thoroughly studied the appropriate numerical method, wrote by himself a computer
program and knew its features, limitations and drawbacks, now, using the computer
transformations, which have become standard, the engineer often has no idea of the
298 V. Ochkov
wide range of their features. This is indicated by Mathcad users writing to the
forums, and there are a lot of tips in this book on how to avoid an error or correct a
mistake.
It is obvious that under present conditions of the developed applied technologies,
the vast majority of students do not need the mathematical knowledge, and prob-
ably that will not be useful to them in their future work. And then what will they
need? There are three global challenges put before a teacher of mathematics in a
technical college. The first—to provide the student with a minimum mathematical
tools technique required for the study of special disciplines. The second—to clarify
as strictly as possible those mathematical foundations on which the research and
computing processes are based. The third—to lay the fundamentals of a complex
method of formulating and solving problems, to teach the student to check and
analyze the result and provide a means of verification.
So, what is needed in this regard from mathematics by a mechanical engineer
and what heat engineering in a frame of Mathcad can give to teachers of mathe-
matics? We will be able to answer immediately to both questions by analyzing what
mathematical tools technique was used in the solution the heat engineering prob-
lems discussed in this book.
To perform thermal design calculations using Mathcad means the user shall have
at least basic skills of matrix and vector calculus. He should have a clear idea about
matrix, vector and scalar and their properties. Transposition, addition, multiplica-
tion and inversion of matrices, calculation of a rank of matrix—these are operations
a heating engineer should feel at home with. Of course, it is not required from him
to have skills to “manually” invert, for example, a fourth-order square matrix. But
he should have a good idea of the fact that if all the elements of a matrix are so
small that the computer perceives them as zeros, the matrix cannot be inverted,
though in theory the inverse matrix exists. In this regard, especially many questions
arise when solving systems of linear algebraic equations, which are one of the main
mathematical calculation tools for thermal problems. To select the most convenient
and the most accurate method of numerical solution of systems of linear algebraic
equations a mechanical engineer should know the strengths and weaknesses of both
Gauss method, and Cramer’s rule, and methods of solving matrix equations.
Knowing the general theory of linear equations systems, the Mathcad user can
either be confident of the correctness of the computational result, or envisage the
appearance of a computational “incident”, or even make discovery in engineering
when analyzing the resulting contradiction between the expected theoretical and
obtained computational result.
Another focal point of the study seems to be functions of one real variable. One
theoretically needs to know everything about that. In the thermal calculations in
Mathcad, symbolic operations on functions absolutely equally co-exist with
approximate methods of calculation applicable to functions. There is virtually no
chapter in this book, where actions relating to the functions of one variable were not
used in solving thermal problems. This is calculation of values, interpolation and
extrapolation, the search of roots and creating the inverse function, character or
numerical differentiation and integration, the search for extrema and study of the
20 Three-Layer Thermal Engineering Cake or a Conclusion 299
useful for all the students in the group. To be fair we note that the students working
with Mathcad have much less difficulty in solving math problems than the other
students.
Thus, it is evident that the structure, sequence and technique of sharing math-
ematical knowledge requires changes [74]. It is necessary to review the entire
teaching material on higher mathematics to eliminate or edit some math sections
and include others, more relevant to the future profession. This choice of material
should be conducted taking into account the practical importance of the studied
sections. As for the numerical methods, it is not necessary to teach methods or their
implementation (which is a task of computing software like Mathcad), but under-
standing of the nature and principles.
One of the teachers said, “If we do not know what to teach, let’s go back to the
meaning.”
And the last thing. In almost all the Studies of the book there is water. And there
is nothing surprising about that, water—is the source of life. No wonder the space
probes are searching at distant planets for water as the first sign of possible life on
them. In our book we discussed the technical principles for using water as the
working fluid and coolant. But, as the poet said, (http://www.stihi-rus.ru/1/Sluckiy/
40.htm), “For some reason physics are kept high, and lyricists are kept somewhat
down.” In [75] heat engineering, chemical, biological, economic, and other scien-
tific and technical aspects of water are supplemented with components of philo-
logical, art history, psychology and even religious nature…
References
1. Ochkov VF, Utenkov VF, Orlov KA (2000) Teplotekhnicheskiye raschety v srede Mathcad
(Thermal Engineering calculations with Mathcad). Teploenergetika J 2:73-78. http://twt.
mpei.ac.ru/ochkov/WSPHB/TE-2-2000.pdf. Accessed Sept 2015
2. Yankov GG (2011) Tekhnicheskaya kniga v sovremennoy Rossii: shans na vyjivanie
(Technical book in modern Russia: the chance of survival). University Book J 6:67–70.
http://www.unkniga.ru/vuz/101-tehknoga.html. Accessed May 2011
3. Invernizzi CM (2013) Closed power cycles: thermodynamic fundamentals and applications.
Springer, Berlin
4. Kirillin VA, Sychev V, Sheindlin AE (2008) Tekhnicheskaya termodinamika (Technical
Thermodynamics). Publishing House MPEI, Moscow
5. Sychev VV (2009) Slozhnyye termodinamicheskiye sistemy (Complex thermodynamic
systems). Publishing House MPEI, Moscow
6. Aleksandrov AA (2006) Termodinamicheskiye osnovy tsiklov teploenergeticheskikh
Ustanovok (The thermodynamic basis of cycles of thermal power plants). Publishing
House MPEI, Moscow
7. Andryushchenko A (1967) Osnovy tekhnicheskoy termodinamiki real’nykh protsessov
(Fundamentals of engineering thermodynamics of real processes). Higher School Publishing
House, Moscow
8. Ochkov VF (2002) Fizicheskiye i ekonomicheskiye velichiny v Mathcad i Maple (Physical
and economic value of Mathcad and Maple). Finance and statistics Publishing House,
Moscow
9. Korobov VI, Ochkov VF (2012) Chimichny rozrakhunki v sededovitschi Mathcad
(Chemical calculations in Mathcad). Dnipropetrowsk National University Publishing
House, Dnipropetrowsk
10. Ochkov VF (2009) Mathcad 14 dlia studentov i injenerov (Mathcad 14 for students and
engineers). BHV-Petersburg, Saint Petersburg
11. Korobov V, Ochkov VF (2011) Chemical kinetics with Mathcad and Maple. Springer, Berlin
12. Ochkov VF (2007) Mathcad 14 dlia studentov, injenerov I dizaynerov (Mathcad 14 for
students, engineers and designers). BHV-Petersburg, Saint Petersburg
13. Ochkov VF (2005) Mathcad 12 dlia studentov i injenerov (Mathcad 12 for students and
engineers). BHV-Petersburg, Saint Petersburg
14. Solodov AP, Ochkov VF (2004) Differential Models: An Introduction with Mathcad.
Springer
15. Vockoboynikov YE, Ochkov VF (2002) Programmirovaniye i resheniye zadach v pakete
Mathcad (Programming and problem solving in the package Mathcad). Publisher
House NSABU, Novosibirsk
16. Ochkov VF (2001) Sovety pol’zovatelyam Mathcad (Tips and tricks for Mathcad users).
Publishing House MPEI, Moscow. http://twt.mpei.ac.ru/ochkov/Units/Forword_book.htm.
Accessed Sept 2001
50. Ochkov VF (1997) Kak ya prodaval programmu (How I sold a program). ComputerPress J 7.
url http://twt.mpei.ac.ru/ochkov/paint.htm. Accessed Sept 2000
51. Ochkov VF, Gavrilenko SS (2010) Programmirovannyy vybor fil’trov (Programed choose of
filters). Water Treat Conditioning Supply J 7. url http://twt.mpei.ac.ru/ochkov/Opt-Filtr/
Filter-2.pdf. Accessed Nov 2010
52. Ochkov VF, Gurke S, Znamensky VE (2014) Rabota na forumakh — vazhnyy aspekt
inzhenernoy deyatel’nosti (Work on forums—an important aspect of engineering). Abstracts
of International Scientific Conference Informatization of engineering education Inforino. url
http://twt.mpei.ac.ru/ochkov/CloudUniOchkov.pdf. Accessed Sept 2014
53. Zorin VM (2008) Atomnyye elektrostantsii (Nuclear power plants. Main technological
process). Moscow Power Engineering Institute Publishing House, Moscow
54. Ochkov VF, Orlov KA, Ivanov EN, Makushin AA (2013) Calculation and visual displaying
of the water chemistry conditions in return cooling systems at thermal power stations. Therm
Eng J 60:465–471. http://twt.mpei.ac.ru/ochkov/WSPHB/ThermEng1307007OchkovKOR.
pdf. Accessed Sept 2013
55. Ochkov VF, Orlov KA, Ochkov AV, Znamenskiy VE (2013) Application of the cloud
technologies to perform heat engineering calculations. Therm Eng J 60(9):679–685. url
http://twt.mpei.ac.ru/ochkov/WSPHB/CloudFormulas.pdf. Accessed Sept 2013
56. Ochkov VF, Ustyuzhanin EE, Ko CK, Shishakov VV (2015) Thermophysical databases:
from tables to interactive internet resources and cloud templates. High Temp J 53(4):515–
520. url http://twt.mpei.ac.ru/ochkov/TBT-BD-Internet-eng.pdf. Accessed Nov 2015
57. Ochkov VF (2013) Zhivyye kinematicheskiye skhemy v Mathcad (Live kinematic diagrams
in Mathcad). Open Educ J 3:27–33. http://twt.mpei.ac.ru/ochkov/Mathcad-15/kinematic.html.
Accessed Sept 2013
58. Ochkov VF, Aleksandrov AA, Orlov KA (2007) Termodinamicheskiye tsikly: raschety v
Internete (Thermodynamic cycles: Online calculations). Vestnik MPEI J 1:43–50. http://twt.
mpei.ac.ru/ochkov/Therm_Cycle_Art/index.html. Accessed Sept 2015
59. Wang D, Zhi LH (2007) Symbolic-numeric computation. Springer, Berlin
60. Velichko VI, Solokhin VI (2014) Teploprovodnost’: zadachi i uprazhneniya (chislennoye
resheniye) (Heat conduction: problems and exercises (numerical solution). Publishing
House MPEI, Moscow
61. Ochkov VF, Orlov KA, Kondakova GY, Ochkov AV (2009) Trenazher dispetchera teplovoy
seti (Heat grid dispatcher training simulator). Energetic J 4:38–41. url http://twt.mpei.ac.ru/
ochkov/trenager/ThermoNet/page38.png (39, 40 and 41.png). Accessed Sept 2009
62. Ochkov VF, Mescheryakov SV, Orlov KA et al (2007) Trenazher dlya podgotovki personala
teplovykh setey (Training simulator for heat grid staff). New Russ Power Ind J 1:33–43.
http://twt.mpei.ac.ru/ochkov/trenager/ThermoNet/Article_NRE.html Accessed Feb 2007
63. Ochkov VF, Orlov KA, Chzho Ko Koa, Anokhin DN (2012) Oblachnyye funktsii dlya
inzhenernykh raschetov vodosnabzheniya (Clod functions for water supply engineering
calculations). Water Supply Sewerage J 9–10:68–74. url http://twt.mpei.ac.ru/GDHB/
CloudFunction.pdf. Accessed Sept 2012
64. Ochkov VF, Bogomolova EP (2015) Teaching mathematics with mathematical software.
J Humanistic Math 5(1):265–285. url http://twt.mpei.ac.ru/ochkov/teaching_of_math.pdf.
Accessed July 2015
65. Plis AI, Slivina NA (2003) Matematicheskiy praktikum dlya inzhenerov i ekonomistov
(Mathcad. Mathematics workshop for engineers and economists. Education guidance).
Finance and Statistics, Moscow
66. Ochkov VF (2005) Sneg v megapolise i energetika (Snow in a city and power industry).
Power Supply Water Treatment J 1. url http://twt.mpei.ac.ru/ochkov/Snow/index.html.
Accessed Sept 2005
67. OchkovVF (2010) Sneg v bol’shom gorode (Snow in a big city). Water Treat Conditioning
Supply J 3(27):72–74. url http://twt.mpei.ac.ru/ochkov/Snow/index-vvv.html. Accessed Sept
2010
References 307
68. Ochkov VF, Bogomolova EP (2015). Eto strashnoye slovo difury (This Terrible Word
“Difeqs”). Inf Sch J 1:55–58. url http://twt.mpei.ac.ru/ochkov/ODE.pdf. Accessed Sept 2015
69. Sokolov EY, Brodiansky VM (1981) Energeticheskiye osnovy transformatsii tepla i
protsessov okhlazhdeniya (Energy transformation basics of heat and cooling processes).
Energoatomizdat, Moscow
70. Cengel Y, Boles M (2006) Thermodynamics: an engineering approach, 5th edn.
McGraw-Hill Higher Education, Boston
71. Moroziuk TV (2006) Teoriya kholodil’nykh mashin i teplovykh nasosov (Theory of chillers
and heat pumps). Negoziant Publishing House, Odessa
72. Ulam S (2001) Priklyucheniya matematika (Mathematician’s adventures). RHD Publishing
House, Izhevsk
73. Wiener N (1964) I am a mathematician. MIT Press, Cambridge
74. Yartseva AD, Chernavsky AV (2012) O matematike: problemy prepodavaniya (About
mathematics, problems in teaching). Publishing House Znak, Moscow
75. Ochkov VF (2009) Voda v nauchnom, kul’turologicheskom i religioznom aspektakh (Water
in a scientific, cultural and religious aspects). Water Treat Condition Supply J 7:64–70. url
http://twt.mpei.ac.ru/ochkov/WSCR/WaterSCR.pdf. Accessed Sept 2009
76. Miyagawa K, Hill PG (2001) Rapid and accurate calculation of water and steam properties
using the tabular Taylor series expansion method. J. Eng. Gas Turbines Power 123(3):707–
712. doi:10.1115/1.1367340
77. Kunick M, Kretzschmar HJ, Gampe U (2008) Fast calculation of thermodynamic properties
of water and steam in process modelling using spline interpolation. In: Proceedings of the
15th international conference on the properties of water and steam, Berlin, 8-11 Sept 2008
78. Wagner W, Kretzschmar HJ (2008) International steam tables—properties of water and
steam based on the industrial formulation IAPWS-IF97. Springer, Berlin
79. Parry W, Bellows J, Gallagner J, Harwvey A (2000) ASME steam tables for industrial use.
ASME Press, New York
80. Tsanev SV, Burov VD, Remezov AN (2006) Gazoturbinnye i parogazovye ustanovki
teplovyh elektrostancii (Gas turbine and combined cycle gas turbine thermal power plants).
Publishing House MEI, Moscow
81. Andryushchenko AI (1998) O termodinamicheskoy effektinvosty slozhnyh tsiklov GTU v
parogazovyh ustanovkah (The thermodynamic efficiency of complex cycles of gas-turbine
units in combined-cycle installations). Teploenergetika J 3
82. Andryushchenko AI (1975) Osnovy termodinamiki tsiklov teploenergeticheskikh ustanovok
(The basis of thermodynamics of cycles of thermal power plants). Higher School, Moskow