=Paper=
{{Paper
|id=Vol-2131/paper9
|storemode=property
|title=Usage of the MATHCAD Framework for Energy Saving Simulation
|pdfUrl=https://ceur-ws.org/Vol-2131/paper9.pdf
|volume=Vol-2131
|authors=Alexander G. Gorokhovsky,Elena E. Shishkina,Natalia R. Vukovic
}}
==Usage of the MATHCAD Framework for Energy Saving Simulation==
https://ceur-ws.org/Vol-2131/paper9.pdf
Usage of the MATHCAD Framework for Energy Saving Simulation
Alexander G. Gorokhovsky Elena E. Shishkina
Ural State Forestry Engineering University Ural State Forestry Engineering University
620100, Russia, Ekaterinburg, 620100, Russia, Ekaterinburg,
goralegr@yandex.ru elenashishkina@yandex.ru
Natalia R. Vukovic
Ural State Forestry Engineering University
620100, Russia, Ekaterinburg,
shpak@usfeu.ru
Abstract
This study explores the use of energy savings technologies in timber pro-
duction for wood housing construction. Methods applied were based on
the creation of system of optimized regimes of drying of sawn timber was
carried out in two stages: computing experiment with the tools of
MathCAD and production experiment. The findings discovered the fol-
lowing tendency: however the minimum cost of energy is achieved at the
minimum time of drying that is explained by an essential difference in the
cost of thermal and electric energy. The regimes of steeples structure re-
ceived by computer modeling allow receiving sawn timber that fully
meets the requirements of the consumer without application of moisture
treatment and the conditioning processing. The developed technique of
formation of the steeples regimes of drying of sawn timber allows deter-
mining structure and sizes of parameters of the regime depending on the
required quality of drying and energy costs on its carrying out. Calcula-
tions of technical and economic efficiency have shown that the total an-
nual economic effect in more than 17 dollars per 1 m3 and it can be ex-
tended by including on power component.
Introduction
Glued wall beam is one of main constructional materials in wooden housing construction. In this case the initial ma-
terial in production of glued wall beam is sawn timber. The conducted researchers [1] showed that electric power con-
sumption for chamber drying of sawn timber 1.5-2 times exceed those for their sawing. In addition, in the structure of
complete cost of chamber drying of sawn timber the share of energy component is quite considerable about 60%, and in
the total cost of drying of sawn timber up to 30%. Therefore, creation of energy-saving technology of drying of sawn
timber is an economically efficient task.
This problem is highly actual for modern scientist. So the compression strengths are compared with steel and timber
wall stud strengths and shown to be suitable for residential building applications. The combined plain channel and stiff-
ened channel experimental data covers a broad range of section slenderness values, and design models are developed to
predict their compression strength are investigated in the articles of Bambach, M. R. (2018) [2]. Also The European
Committee for Standardization (CEN) thus developed horizontal standards to enable the sustainability assessment of
construction works over their entire life cycle, which were analyzed in the works of Achenbach, H., Wenker, J. L., &
Rüter, S. (2018) and in it for the categories GWP and AP, around 30% of the impacts originate from the prefabrication
of the building elements, their transport and the processes at the construction site [3].
A tangible reduction of energy consumption is possible in the following directions [1]:
-costs on heating of outside air on the basis of air exchange between the chamber and surrounding atmosphere;
-electric power consumption for the drive of fan.
It should be noted that development of drying technology goes now mainly by improvement of the modes of drying
on the basis of modern methods of computer modeling and optimization [2]. Surely that optimization of the modes of
drying should not only increase their energy efficiency, but also provide quality of drying of sawn timber allowing mak-
ing further both a glued wall beam and other types of products for housing construction.
The issues concerning effect of drying regimes on the quality of drying of sawn timber were in detail investigated in
the 1950s-1980s [5,7, etc.].
P.S. Sergovsky, one of the founders of Russian science on wood drying [6,7,8]notes that not all quality indicators but
only two of them depend on drying
regime: integrity of the material caused by the size of internal tension in the wood and level of preservation of
strength in the wood caused by the level and duration of temperature impacts on it.
Creation of drying regimes should be carried out in such a way that during the whole process the maximum values of
internal tensions in wood do not exceed the maximum permissible value. The regime is characterized by a safety coeffi-
cient:
пр. р.
В (1)
макс
where ζпр.р. – calculated strength of wood;
ζмаx – the maximum size of internal tension.
If В < 1, then the regime does not provide integrity of material, if В > 1, the greatest possible intensity of process is
not reached. However, the Guidance Technical Materials on physical and mechanical properties of wood for strength of
wood give a variation factor of 10%. At the same time efficiency of drying can be completely guaranteed at B = 1.3
(with probability of p = 99.73%). Respectively, at B = 1.2 it is guaranteed with probability p = 95%, and at B = 1.1, p =
90%. The set of parameters of regime uniquely determines the size of indicators of efficiency and quality of dried up
wood [9].At the same time efficiency of drying can be unambiguously estimated by total power consumption per 1m3 of
sawn timber. Consequently, the task on increasing of data values of these indicators can be considered as optimization.
Methodology
Creation of system of optimized regimes of drying of sawn timber was carried out in two stages:
1. Computing experiment
2. Production experiment
Removal of moisture out of wood during the drying process is rather complex physical and chemical process accom-
panied by heat mass exchange (HME).
For mathematical description of the process of low-temperature convective drying of unlimited sawn timber, A.V.
Lykov [10] proposes the following system of differential equations in private derivatives (DEPD HME).
t u ,
a 2 t (2)
c
u
a м 2 u a м 2 t , (3)
For unlimited sawn timber the initial and boundary conditions of the III kind have the form:
t x0 ,0 f x , (4)
ux,0 x , (5)
t R,
x
t c t R, 1 б m uR, u p 0 (6)
uR, t R,
am
x
a m
x
m u R, u p 0 (7)
Symmetry condition:
t 0, u 0, (8)
0
x x
Where t – temperature, oC;
u – humidity;
η – time;
a – coefficient of thermal diffusivity, m2/c;
am - coefficient of moisture conductivity, m2/c;
ε – coefficient of phase transformation;
ρ – density of wood, kg/m3;
c – thermal capacity of wood, kJ;
δ – thermogradient coefficient;
x – coordinate in the direction of thickness of a plate, m;
R – a half of thickness of a plate, m;
λ – coefficient of thermal conductivity, W/(м·град);
α – coefficient of heat exchange, W/(м2·град);
αm – coefficient of moisture exchange, м/с;
ρб – basic density of wood, кг/м3;
tс – temperature of the medium, 0С;
uр – equilibrium humidity of wood.
For the solution of the system (2) – (10) software was developed in MathCAD computing environment [11] on the
basis of the implicit method [12]. In addition, software was used to calculate internal stresses [13] on the basis of multi-
rod model of the board [14]. Using the above – mentioned software, a computer experiment was implemented. During
the experiment, constant factors were the following:
1. Type of sawn timber – conditional (pine, section 40x150mm);
2. Type of regime – step less
- temperature of processing medium
u н u
t c t н t к t н ; (9)
u н 0,1
- equilibrium humidity
u р u рк u рн u рк е е
b0 b1u
, (10)
where tн, tк – respectively, initial and final temperature of the agent of drying, 0С;
uн, u – respectively, initial and current humidity of wood;
uрн, uрк – respectively, initial and final values of equilibrium humidity ;
b0, b1 – coefficients.
Expression (12) represents function of desirability [15] which is characterized by two transition values u, denoted re-
spectively uп1 and uп2. Moreover, uп2 = 0,35 and uн = 0,6 remained constant in all experiments. Variable factors during
the experiment are:
uрн (х1), uрк (х2), uп1 (х3), tн (х4), tк (х5).
The factors varied at three levels, their values in encoded and natural terms are presented in Table 1.
Table 1: Variable factors during a computing experiment
Values of factors at levels
№ factors low main upper
coded natural. coded natural coded natural
1 uрн (х1) - 0,1 0 0,14 + 0,18
2 uрк (х2) - 0,02 0 0,03 + 0,04
3 uп1 (х3) - 0,1 0 0,15 + 0,2
4 tн (х4) - 60 0 70 + 80
5 tк (х5) - 80 0 90 + 100
Output parameters:
η1 (у1) – duration of drying of sawn timber to humidity of W = 12 %;
η2 (у2) – duration of drying of sawn timber to humidity of W = 7 %;
Sт (у3) – difference of humidity on board thickness;
Sw (у4) – an average square deviation of humidity ;
Bmin (у5) – the minimum value of safety criterion of the regime during each drying;
Qт (y6) – costs of thermal energy for drying of 1m3 of sawn timber;
Сэ (у7) – total costs of energy (electric and thermal) spent on drying of 1m3 of sawn timber.
In the course of experiment Hartly plan [15] consisting of 27 issues (experiments) was implemented. As a result, the
dependences of each output parameter on each input one were obtained in the form of polynoms of the second range.
Then optimization on each of output parameters was carried out (η1, η2, Sт, Sw, Bmin).
Formulation of optimization problems was the following:
y1 min (11)
1 X 1
y 2 min (12)
1 X 1
y3 min (13)
1 X 1
y 4 min (14)
1 X 1
y5 min (15)
1 X 1
x1
x2
where
X x3
x4
x5
Production tests were carried out for the purpose of check of practical suitability of received optimized regimes of
drying.
As the main experiment plan of the kind B3 was performed.
Calculation of consumption of thermal and electric energy was carried out by a standard technique [4]. The price of
energy unit was taken as average for the Ural region of the Russian Federation.
Results and Analysis
Results of optimization performed in MathCAD computing system are given in Table 2.
However results of optimization of the regime on the required category of quality are of the greatest interest (table 3).
Table 2: Results of optimization of the regime of drying on private criteria
operating factor Values of operating factors for criteria of optimality
№
η1 (hour) η2 (hour) ST Sw Bmin
1 uрн 0,1 0,1 0,18 0,18 0,18
2 uрк 0,02 0,02 0,04 0,04 0,02
3 uп1 0,2 0,174 0,1 0,1 0,1
4 tн, 0С 80 80 60 60,5 80
5 tк, 0С 92 100 100 99,5 100
Value of criterion of optimality
90,8 120 0,011 0,088 2,039
Table 3: Results of optimization of drying regimes on categories of quality
№ operating factor Values of operating factors
/ criterion of quality criteria of optimality
I II III
1 uрн 0,168 0,11 0,1
2 uрк 0,036 0,033 0,038
3 uп1 0,1 0,1 0,2
4 tн, 0С 66,25 74,7 80
5 tк, 0С 100 100 100
6 η2 (hour) 247 179 143
7 Bmin 1,703 1,299 1,263
8 ST 0,019 0,03 0,035
9 Sw 0,01 0,015 0,02
Formulation of the problem was the following:
Ist category of quality:
2 min
1 X 1
(16)
Bmin 1,3
S w 0,01
S T 0,02
II nd category of quality:
2 min
1 X 1
(17)
Bmin 1,2
S w 0,015
S T 0,025
III rd category of quality:
2 min
1 X 1
Bmin 1,2
S w 0,02
S T 0,035
At the third stage, optimization of regimes in energy efficiency parameters was carried out. Two problems of optimi-
zation of these indicators are formulated on the basis of the data obtained at the previous stages of optimization for dry-
ing on the second (II) category of quality:
(19)
}
(20)
}
Results of optimization are shown in Table 4
Table 4 : Results of optimization of drying regimes in terms of energy efficiency parameters
Values of operating factors in terms of optimization on
operating factor
№ parameters
/ criterion of quality
Qт Сэ
1 uрн 0,11 0,1
2 uрк 0,038 0,0385
3 uп1 0,11 0,1
4 tн, 0С 62,4 71,6
5 tк, 0С 81 98,6
6 η1, hour 133,3 122,6
7 Qт, GJ/m3 1,867 1,882
8 Qэ, GJ/m3 0,581 0,533
9 Qт + Qэ, GJ/m3 2,448 2,415
10 Сэ, rubles 860 823
Note – Qэ – amount of electricity
The values of parameters of the regime received as a result of optimization by solving a compromise task by the
method of conditional centre of masses rather closely coincide with the regime parameters for the second (II) category
of quality received by analytical optimization (table 5).
Table 5: Results of analytical and experimental optimization of drying regimes for conditional sawn timber
Values
parameters of the regime analytical optimization
№ /values of output parame- On the second on Qт
experimental optimization
ter (II) category on Сэ
of quality
1 uрн 0,11 0,11 0,1 0,114
2 uрк 0,033 0,038 0,0385 0,033
3 tн, 0С 74,7 62,4 71,6 73,4
4 tк, 0С 100 81 98,6 95
5 Τdrying (Wк = 12 %), hour 125 133,3 122,6 1,38
Sw, % (category of
6 1,5 (II) - - 0,835 (I)
quality)
inside tensions, (category
7 1,376 (II) - - I – II
of quality)
Consumption of energy
8 - 2,448 2,415 2,769
on drying, GJ/m3
Also the drying time is quite close: in the experiment, it is 9.4% higher that can be attributed primarily to idealization
of conditions of drying in a computing experiment. It should be noted that experimentally received valued of an average
quadratic deviation of wood humidity is significant, it is nearly 1.8 times less than those, obtained analytically. At the
same time, according to the experiment this indicator of quality of drying completely corresponds to the first (I) catego-
ry of quality. Obtained discrepancies in our opinion should be attributed to the error of method of analytical determina-
tion of Sw.
The fact is that this technique is calculated on final quantity of stages of the regime of drying and for the steeples re-
gimes we considered that, hypothetically, the quantity of steps was 10. Obviously, that for increase in accuracy, it is
necessary to increase quantity of these hypothetical steps to 20 -25. Besides, it should be noted that results of theory and
the experiment rather closely coincide during drying by standard 3-staged regimes. Thus, drying time in the experiment
differs from theoretical only by 5.5%. As in the previous case, there are discrepancies between theory and the experi-
ment concerning Sw. It is true, to be fair, it should be noted that in this case these discrepancies are much less and their
amount is about 16%. It also can be explained by the fact that in theory and experiment the number of steps (stages) of
regime was still the same -3. It should also be noted that the application of standard regimes, both in theory and in prac-
tice, allowed obtaining quality of drying corresponding only to the third (III) category of quality. As for energy con-
sumption, in experimental optimization they are more than in analytical one by 10-12% that confirms a rather reliability
of carried out computing experiment.
Acknowledgement
The work is carried out based on the task on fulfilment of government contractual work in the field of scientific ac-
tivities as a part of base portion of the state task of the Ministry of Education and Science of the Russian Federation to
Ural State Forest Engineering University (the # 26.8660.2017/8.9 "The Research Methodology of Forms of Economic
and Technological Reality in the Aspect of Sustainable Forest Management").
Conclusions
Optimization of drying regimes in parameters of energy efficiency has shown a rather close coincidence of parame-
ters of the regime by optimization on the consumption of thermal energy and the total cost of energy. However the min-
imum cost of energy is achieved at the minimum time of drying that is explained by an essential difference in the cost
of thermal and electric energy. The regimes of steeples structure received by computer modeling allow receiving sawn
timber that fully meets the requirements of the consumer without application of moisture treatment and the conditioning
processing. The developed technique of formation of the steeples regimes of drying of sawn timber allows determining
structure and sizes of parameters of the regime depending on the required quality of drying and energy costs on its car-
rying out. Calculations of technical and economic efficiency have shown that the total annual economic effect at the
volume of drying of 3000 m3 can reach 1 million rubles, including on a power component – more than 750 thousand
rubles, which is equivalent to more than 17 dollars per 1 m3 and it can be extended by including on power component.
Bibliographic list
[1]. Shishkina E.E. Drying of sawn timber in low-power chambers with natural air circulation: dis. ... cand. Tech.
Sciences / Shishkina Elena Evgenievna. St. Petersburg: SPbGLTA them. CM. Kirova, 2006. - 160 p.
[2]. Bambach, M.R. Geometric optimization and compression design of natural fibre composite structural channel
sections. Composite Structures, 2018 - 185, P. 549-560.
[3]. Achenbach, H., Wenker, J.L., & Rüter, S. Life cycle assessment of product-and construction stage of prefabri-
cated timber houses: a sector representative approach for Germany according to EN 15804, EN 15978 and EN 16485.
European Journal of Wood and Wood Products, 2018. - 76(2), P. 711-729.
[4]. Shishkina E.E. Energy-saving technology of convective drying of lumber based on controlled moisture transfer
in wood: dis. ... Dr. techn. Sciences: 05.21.05 / Shishkina Elena Evgenievna. Arkhangelsk: SAFU them. M.V. Lomono-
sov Moscow State University, 2005. - 336 p.
[5]. Sergovsky P.S., Fakhretdinov A.A. New modes of drying of aspen lumber / / Woodworking industry, 1991. - №
1. P. 4 - 7.
[6]. Sergovsky P.S. Hydrothermal processing of wood. Moscow: Forest Industry, 1975. - 360 p.
[7]. Krasoukhina L.P. On rational regimes for drying birch lumber in batch chambers // Derevoobrabatyvayushchaya
promyshlennost.1988. - № 6. P. 5 - 7.
[8]. Sergovsky P.S. On rational regimes for drying lumber in air chambers of batch action, Derevoobrabatyvayush-
chaya promyshlennost, 1969. - No. 2, pp. 1 - 4; № 3. P. 1 - 4.
[9]. Gorokhovsky A.G. Technology of drying lumber based on modeling and optimization of heat and mass transfer
in wood: diss. ... Dr. techn. / Gorokhovsky Alexander Grigorievich. St. Petersburg. 2008. - 263 p.
[10]. Lykov A.V. On systems of differential equations of heat and mass transfer in capillary-porous bodies // Engi-
neering-Physics Journal. 1974. T. XXVI. № 1. P. 18 - 25.
[11]. Kiryanov D.V. Mathcad - 12. St. Petersburg: BHV - Petersburg. 2005.
[12]. Gavrilova R.I. Investigation of the drying process with variable coefficients of heat and mass transfer // Engi-
neering and Physics Journal. 1964. - T. VII. No. 8. P. 37 - 42.
[13]. Ugolev BN, Lapshin Yu.G., Krotov E.V. Stress control during drying of wood. Moscow: Forest Industry, 1980.
- 206 p.
[14]. Pizhurin AA, Rosenblit M.S. Basics of modeling and optimization of woodworking processes. M.: Forest In-
dustry. 1988.
[15]. Pen R.Z. Statistical modeling and optimization of the pulp and paper industry. Krasnoyarsk: Publishing house
of Krasnoyarsk University. 1982.