Next Article in Journal
Goal-Oriented Tuning of Particle Filters for the Fault Diagnostics of Process Systems
Next Article in Special Issue
Intake Valve Profile Optimization for a Piston-Type Expander Based on Load
Previous Article in Journal
Design and Experimental Analysis of an Adhesive Joint for a Hybrid Automotive Wheel
Previous Article in Special Issue
Modeling of an Organic Rankine Cycle Integrated into a Double-Effect Absorption System for the Simultaneous Production of Power and Cooling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computing the Thermal Efficiency of Autoclaves during Steaming of Frozen Prisms for Veneer Production at Changing Operational Conditions

1
Faculty of Forest Industry, University of Forestry, Kliment Ohridski Blvd., 10, 1797 Sofia, Bulgaria
2
Institute for Building Materials, ETH Zürich, CH 8093 Zürich, Switzerland
3
Department of Wood Science and Engineering, Luleå University of Technology, 931 77 Skellefteå, Sweden
*
Author to whom correspondence should be addressed.
Processes 2023, 11(3), 822; https://doi.org/10.3390/pr11030822
Submission received: 17 February 2023 / Revised: 6 March 2023 / Accepted: 7 March 2023 / Published: 9 March 2023
(This article belongs to the Special Issue Advances in Thermal Process Engineering and Simulation)

Abstract

:
A methodology for the computation of the thermal energy efficiency of modes for the heat treatment of frozen wooden prisms in an autoclave with saturated water vapor at changing operational conditions has been proposed. The methodology includes computer simulations with two own-coupled unsteady models: one to calculate the 2D temperature distribution in the cross-section of prismatic wood materials during their heat treatment, and the second to determine the heat balance of industrial autoclaves for such wood treatment. Simulations were carried out in order to determine the duration, energy consumption, and thermal efficiency of different modes, caused by changed operational conditions, for the autoclave steaming of frozen beech prisms with industrial parameters in the absence and presence of dispatcher intervention. The influence of nine combinations between the time of dispatcher intervention and the degree of reduction of the constant maximum temperature from the 130 °C of the basic mode on the thermal efficiency of the autoclave was investigated. The results show that all studied dispatching interventions cause an increase in both the duration and the thermal efficiency of the modes. This efficiency in the modes at changing operational conditions has values between 68.7% and 74.6%, while the efficiency in the basic steaming mode is equal to 68.0%.

1. Introduction

It is well known that for the decorative coating of details intended for the production of high-quality furniture as well as for the furnishing of office premises, the most frequently utilized practice is to use a veneer, obtained from variable wooden species, with different thickness, texture and color than the wood being heat treated [1,2,3,4]. In the manufacture of veneer, wood materials are subjected to heat treatment with saturated steam or hot water in various facilities in order to plasticize them before cutting the veneer [5,6,7,8,9,10,11,12,13,14,15,16]. This is determined by the circumstance that the heated wood has an increased deformation capability and is susceptible to cutting.
To reduce the duration and energy consumption of heat treatment, instead of using equipment operating at atmospheric pressure, autoclaves are used [11,12,13,14,16,17,18,19,20,21,22,23,24,25,26,27].
Heat treatment modes are usually presented in the literature sources, which are applied in practice to ensure the minimum possible duration of the wood heating process and maximum productivity of the equipment [5,6,7,8,9,10,11,12,13,14,28,29,30,31,32,33,34,35,36]. These sources completely lack data on modes in which dispatcher intervention increases their duration in order to ensure the necessary wood plasticity is achieved later than the usual time due to changed operational conditions. Only in [37,38] are results given for a computer-simulated study of the duration of modes for the steaming of non-frozen beech prisms in an autoclave and of the energy needed for heating the wood itself, when applied to the operations of various dispatcher interventions. It was established that during the steaming of prismatic materials with a cross-section of 0.4 × 0.4 m, u = 0.6 kg·kg−1 and tw0 = 0 °C in a basic mode with tm-max = 130 °C = const and a duration of 9.15 h, maximum productivity of the autoclave is ensured at an energy consumption of 65.35 kWh·m−3, which is required only for warming up of the wood. The dispatcher’s lowering of tm from 130 °C to 100 °C in the 3rd, 5th and 7th h of the basic mode causes an increase in its duration to 14.15, 13.15 and 12.15 h, respectively, while the indicated energy consumption is reduced in all studied cases to about 58.0 kWh·m−3. Information only on the duration of analogous modes with dispatcher intervention, intended for the heat treatment of frozen beech prisms with saturated water vapor, is published in [39]. The partial results are reflected below in the last column of Table 1.
In [14], it was found that during the steaming of beech materials with cross-sectional dimensions of 100 × 100 mm and u = 0.6 kg·kg−1 in an autoclave for 7 h on separate modes with tm = 100, 120 and 140 °C, the consumption of heat energy changes from 97.8 to 155.4 kWh·m−3 at initial wood temperature tw0 = 0 °C; and from 135.1 to 196.3 kWh·m−3 at tw0 = −20 °C. References [9,10] note that the efficiency of the heat treatment of wood materials with saturated water vapor in pits during veneer production does not typically exceed 18%.
Studying the impact of time- and size-varying dispatcher interventions on the duration and energy parameters of wood heat treatment modes has significant scientific and practical interests. In the event of organizational or technical problems in the production lines for the cutting and drying of veneer, it is required to extend the modes for the heat treatment of the wood until the problem is solved, while at the same time ensuring the necessary plasticity of the material. In such cases, it is necessary for the dispatcher or automated control system to intervene in a timely manner and to change the temperature-time parameters of the current mode in an appropriate way.
The proposal and application of a methodology for studying the influence of different dispatcher interventions in wood heat treatment modes will allow, in the future, for the development of software for advanced model-based systems for the automatic control of various types of such treatments [21,25,40,41,42,43]. Solving a task with such a great complexity and multifactorial nature [16,39] can only be done using adequate multiparameter temperature-time-energy models.
This paper presents a methodology for the application of two personal coupled non-stationary models for computing the thermal energy efficiency of different modes of autoclave steaming on frozen wood prisms intended for the production of veneer for the cases of time- and magnitude-variable dispatch interventions.

2. Materials and Methods

2.1. Materials for Research

Two coupled models, created and verified in [20], and the methodology suggested in [44] were used for simulated research of the thermal efficiency of modes for the steaming of frozen wooden prisms in an autoclave in cases of dispatcher intervention in the temperature-time modes’ parameters.
The study was carried out with ice-containing beech (Fagus sylvatica L.) prisms above the hygroscopic range, which are commonly used in veneer production.
During numerical simulations with the mathematical models presented below, the following parameters of the prisms, which influence τmode and η, were set: d × b × l = 0.4 × 0.4 × 2.5 m, tw0 = −20 °C, ρb = 560 kg·m−3, u = 0.6 kg·kg−1, and u fsp 293.15 = 0.31   kg · kg 1 . Beech prisms of such dimensions and with a moisture content above the hygroscopic range are relatively often subjected to steaming in veneer production in practice. Such prisms, with an initial temperature of −20 °C, contain significant amounts of both free and bound water in a frozen state, the thawing of which will favor the increase of the differences in the durations between the individual steaming modes studied.
The simulations were performed under the following industrial parameters of a well-insulated wood steaming autoclave used in practice: D = 2.4 m, L = 9.0 m, γ = 50%, and qsource = 500 kW. Detailed structural and thermo-physical characteristics of such autoclaves, as well as their application in the woodworking industry, is described in [12,14,20,22,23].

2.2. Modelling of the 2D Unsteady Temperature Change in Prisms

The following 2D model of the temperature change in prismatic wood materials during their heat treatment with saturated water vapor and subsequent cooling (conditioning) in an air medium before veneer cutting was created and verified in [20]:
c w-eff 1 , 2 , 3 · ρ w T τ = div ( λ w-cr grad   T )
at   T ( x , y , 0 ) = T w 0
and at the following prism surface temperatures:
during the steaming process:
T ( x , 0 , τ ) = T ( 0 , y , τ ) = T m ( τ )
during the subsequent conditioning process:
T ( x , 0 , τ ) y = α w-cond ( x , 0 , τ ) λ w-cr ( x , 0 , τ ) [ T ( x , 0 , τ ) T air-cond ( τ ) ]
T ( 0 , y , τ ) x = α w-cond ( 0 , y , τ ) λ w-cr ( 0 , y , τ ) [ T ( 0 , y , τ ) T air-cond ( τ ) ]
The defrosting process of ice-containing wood materials during steaming for the purpose of plasticization in the production of veneer takes place in three stages [14,20,32]. Figure 1 shows these three stages, as well as the symbols of the thermo-physical characteristics of the wood in each of them, for which it is necessary to have a mathematical description when solving the model (1)–(5).
The effective specific heat capacities of the wooden prisms during their defrosting, cw-eff1,2,3, which participate in Equation (1), are described mathematically as follows [6,20,26,32,45]:
First   stage :   c w-eff 1 = c w-fr + c ice-fw
Sec ond   stage :   c w-eff 2 = c w-nfr + c ice-fw
Third   stage :   c w-eff 3 = c w-nfr
where
c w-fr = K c-fr 526 + 2.95 T + 0.0022 T 2 + 2261 u + 1276 · u fsp 272.15 1 + u
K c-fr = 1.06 + 0.04 u + 0.00075 ( T 272.15 ) u fsp 272.15
c ice-bw = ( 69.344 T + 119.183 T · ln T 273.15 ) · ( u fsp 272.15 0.12 ) · exp [ 0.0567 ( T 272.15 ) ] 1 + u
c w-nfr = 1 1 + u · ( 2862 u + 2.95 T + 5.49 u · T + 0.0036 T 2 + 555 )
c ice-fw = 3.34 × 10 5 u u fsp 272.15 1 + u
u fsp 272.15 = u fsp 293.15 + 0.021
The mathematical descriptions of the wood thermal conductivity, λw-cr, wood density, ρw, and heat transfer coefficient, αw-cond, given in [20,26,32,46,47] were used in solving models (1)–(5).
It can be noted that the experimentally established data for λw = f(t,u) and cw = f(t,u) in the dissertations [6,48] were used as a basis in the mathematical descriptions of λw-cr, cw-fr, and cw-nfr, which were involved in Equations (1) and (4)–(12). These data are often used in both the European [8,10,11,12,13,14,28,47,49] and the American [29,50,51,52,53,54,55,56] literature sources, which consider approaches for calculating processes of the defrosting and/or heating of different wood materials.

2.3. Modelling of Thermal Efficiency of Modes for Steaming of Wooden Prisms in Autoclaves

To calculate this efficiency, it is necessary to have a mathematical description of the energy consumption of the entire autoclave and separately of the part of it that is used to heat the wood placed in the equipment.
In the simulations, the total energy consumed by the autoclave, Q a n , was calculated using the following unsteady model of its thermal balance, which was proposed and experimentally verified in [20]:
Q a n = Q w n + Q mb n + Q il n + Q e n + Q fv n + Q cw n
Dependence of all the components of the thermal balance on the multitude of influencing factors have been considered in [14,16,20].
The thermal energy consumption used for heating of the frozen prismatic wood materials at any time, n∙Δτ, Q w n , was calculated by the following equation [20,32,57]:
Q w n = ρ w 3.6 × 10 6 S w · { S w c w-eff 1 , 2 , 3 n at T i , k n + c w-eff 1 n at T w 0 2 · ( T i , j n T w 0 ) d S w }
where
S w = d · b 4
The thermal energy efficiency, η, of the separate steaming modes of prisms in an autoclave is equal to
η = 100 Q w-max n Q a-max n
where Q w-max n and Q a-max n are calculated by the models as the maximum values of Q w and Q a for each of the autoclave steaming modes given in Table 1.
Table 1. Change in the temperature tm1 and duration τsteam of the studied autoclave steaming modes.
Table 1. Change in the temperature tm1 and duration τsteam of the studied autoclave steaming modes.
Steaming
Modes
Ist Stage of tm1,
°C
IInd Stage of tm1,
°C
Ist Stage of tm1
τI, h
IInd Stage of tm1
τII, h
τsteam =
τI + τII, h
Mode 013013.913.9
Mode 1
Mode 2
Mode 3
130
130
130
120
120
120
3.0
7.0
11.0
13.2
8.7
4.2
16.2
15.7
15.2
Mode 4
Mode 5
Mode 6
130
130
130
110
110
110
3.0
7.0
11.0
15.7
10.7
5.7
18.7
17.7
16.7
Mode 7
Mode 8
Mode 9
130
130
130
100
100
100
3.0
7.0
11.0
17.7
12.1
6.5
20.7
19.1
17.5

2.4. Change in the Steaming Medium Temperature Tm of Modes in Cases of Absence and Presence of Dispatcher Intervention

In the numerical simulations with the models (1)–(5) and (15), the change of Tm shown in Figure 2, whose mathematical description is presented in [14,20,58], is assumed.
The following values of the parameters of the modes, whose symbols are given in Figure 2, were set: Tm0 = 253.15 K (i.e., tm0 = −20 °C); Tm1 = 403.15 K (i.e., tm1 = 130 °C); ΔTm1-a = 10 K = 10 °C, ΔTm1-b = 20 K = 20 °C, and ΔTm1-c = 30 K = 30 °C; Tm2 = 383.15 K (i.e., tm2 = 110 °C); Tm3 = 353.15 K (i.e., tm3 = 80 °C); Tm-cond = 293.15 K (i.e., tm-cond = 20 °C); Δτa = 3 h, Δτb = 7 h, and Δτc = 11 h. These values of the Tm parameters fully correspond to those applied in the regimes for the industrial steaming of wood materials in the production of veneer [12,14,22].
With these values, one basic mode without dispatcher intervention (Mode 0) was calculated, as well as 9 modes with dispatcher intervention, which in Table 1 are described as Mode 1 to Mode 9.
Table 1 also provides the temperatures and durations of the two stages (before and after the dispatcher’s intervention) of Mode 1 to Mode 9, and also the total duration of the introduction of steam into the autoclave, τsteam, for all of the investigated modes.
For joint numerical solving of the experimentally verified coupled models presented above and aimed at computation of the duration, thermal energy consumption, and energy efficiency of the modes given in Table 1, a personal software program was created in the Visual FORTRAN computing environment. Using this program, the modes shown in Table 1 were developed.
An explicit finite-difference scheme was used to transform the individual model equations into a FORTRAN-friendly programming form, which excludes any simplifications of the models.
The development of the modes consisted in selecting the values of the parameters shown in Figure 2 and Table 1, and that at the end of the conditioning of the steamed prisms, to ensure the required good plasticity of the wood before cutting the veneer. The degree of plasticity of the steamed prisms is assessed when the temperature distribution in their central cross-section falls completely within the optimal limits recommended for beech wood in veneer production topt-min = 62 °C and topt-max = 90 °C [14,59].
Simultaneously with the solution of the models, computations of tw-avg, Qw, and Qa were carried out. After determining the maximum values of Qw and Qa, the energy efficiencies η of each of the studied modes were calculated.

3. Results

3.1. Computing the 2D Unsteady Temperature Change in Prisms during Studied Modes

In Figure 3, as an example, the calculated change in ts, tw-avg, and t of 2 representative points t1 (with coordinates d/8, b/8) and t2 (with coordinates d/2, b/2) of the prisms during Mode 0, Mode 5, and Mode 8 is presented. Analogous figures for modes 1, 3, 4, 6, 7, and 9 can be seen in [39].

3.2. Computing the Qa, Qw, and η for the Cases of Absence and Presence of Dispatcher Intervention in Steaming Modes

In Figure 4, Figure 5 and Figure 6, the computed values of the energies Qa and Qw during the all studied modes is presented.
In Table 2, the change in τ2 = τsteam, and τ4 = τmode (see Figure 2) of all modes, and also of the tw-avg, Qw-max, Qa-max, and η is given.
Figure 7 shows the change in Qw--max and Qa-max for all steaming modes.
In Figure 8 the change in η of all studied modes is presented.

4. Discussion

It can be seen in Figure 3 that the temperature at the representative points of all the prisms changes along extremely complex curves, both during the steaming modes and during the subsequent conditioning of the heated prisms in an air environment. The temperature of the processing medium in the autoclave at the beginning of all modes rises gradually and reaches the maximum value of tm1 = 130 °C after 2.65 h.
The violation of the smoothness of the line tm = f(τ) around the 1st h at the beginning of the modes is caused by the phenomenon of the uneven melting of the ice formed by the free water in the wood, which is adequately reflected in the mathematical models.
Table 2. Change in τsteam, τmode, tw-avg at τ2, Qw-max, Qa-max, and η of the studied steaming modes, depending on Δτ.
Table 2. Change in τsteam, τmode, tw-avg at τ2, Qw-max, Qa-max, and η of the studied steaming modes, depending on Δτ.
Steaming
Modes
Δτ,
h
τ2 = τsteam,
h
τ4 = τmode,
h
tw-avg
at τ2, °C
Qw-max,
kWh·m−3
Qa-max,
kWh·m−3
η,
%
Mode 0013.917.491.797.95144.0868.0
Mode 1
Mode 2
Mode 3
3
7
11
16.2
15.7
15.2
18.7
18.2
17.7
90.2
90.0
90.5
96.44
96.34
96.69
137.19
137.06
137.49
70.3
70.3
70.3
Mode 4
Mode 5
Mode 6
3
7
11
18.7
17.7
16.7
20.2
19.2
18.2
87.3
87.1
87.9
94.01
93.86
94.50
129.65
129.44
131.21
72.5
72.5
72.0
Mode 7
Mode 8
Mode 9
3
7
11
20.7
19.1
17.5
22.2
20.6
19.0
81.8
81.7
83.0
89.60
89.53
90.55
120.11
119.98
131.75
74.6
74.6
68.7
When this phenomenon occurs, a huge amount of heat is utilized to melt the aforementioned ice, which disturbs the thermal balance of the autoclave in the considered case of limited power of the steam generator, qsource, feeding it, and slows down the rise of tm.
The duration of the basic Mode 0, in which there is no dispatcher intervention, is equal to 17.4 h, and the duration of the supply of steam to the autoclave in this mode is equal to 13.9 h. The earlier such an intervention is carried out and the greater the reduction in the maximum value of tm1, the longer these durations are compared to those in the basic mode.
Figure 4, Figure 5 and Figure 6 show that the increase in the steaming time causes a gradual smooth increase of Qw and Qa in the basic mode for the prisms’ steaming. Due to the influence of the last 5 terms in the right-hand side of Equation (15), the rate of increase of Qa is greater than that of Qw. This means that the differences between the corresponding graphs Qa and Qw are equal to the sum of the energies that provide the energy consumptions described by the last 5 terms in Equation (15).
The increase of Qw is analogous to that of tw-avg of the prisms in all investigated steaming modes, with and without dispatcher intervention.
The smoothness of increasing Qw and Qa during the steaming process is disturbed at the moments of application of dispatcher intervention in the steaming modes of the prisms. The earlier such intervention occurs or the greater the temperature change Δtm1 in the mode, the more significant the difference between the values of Qw and Qa of the basic mode when compared to the corresponding mode with dispatcher intervention.
The thermal efficiency η in the modes at changing operational conditions has values between 68.7% and 74.6%, while in the basic steaming mode, efficiency is equal to 68.0%. If the dispatcher’s intervention was provided in the 3rd or 7th h of the modes, with a more significant reduction in tm1, there was a greater increase in η compared to η of the basic mode.
When applying a dispatcher interference closer to the end of the modes, for example at the 11th h of the modes, this dependence is broken in the mode, with the smallest investigated decrease in tm1 being from 130 to 120 °C.
When the efficiency η is equal to 0.67, i.e., it is quite a bit greater than that of the basic steaming mode. The reason for this is the much higher Qa-max value at the 11th h compared to the Qa-max values at the 3rd and 7th h of the modes with a decrease in tm1 from 130 to 120 °C.

5. Conclusions

This paper considers a methodology for computing the energy consumption and thermal efficiency of autoclaves during the treatment with saturated water vapor of frozen prisms in veneer production at changing operational conditions. When such conditions occur, the parameters of the autoclave’s steaming modes have to be changed by the dispatcher or control system in such a way as to ensure optimal plasticity of the wood immediately before cutting the veneer.
The change in the maximum values of Qw and Qa, respectively Qw-max and Qa-max, and with their change, also that of of η, calculated during computer simulations with two own coupled models for each of the investigated modes for autoclave steaming of beech prisms with industrial parameters, led to the following conclusions:
  • At the moment τ2 = τsteam = 13.9 h, when the introduction of water vapor into the autoclave ends, the greatest values Qw-max = 97.95 kWh∙m−3 and Qa-max = 144.08 kWh∙m−3 are established in the basic mode, which takes place at tm1 = 130 °C = const. These values of Qw-max and Qa-max determine the presence of the lowest value of η = 68.0% of the basic mode compared to the thermal efficiency of all modes with dispatcher intervention.
  • When, upon application of dispatcher intervention, the temperature of the processing medium in the autoclave is reduced from tm1 = 130 °C to tm1 = 120 °C, the energies Qw and Qa reach their maximum values at moments τ2 = τsteam, which depend on the occurrence times of this intervention. Then, they are equal to about Qw-max ≈ 96.5 kWh∙m−3 and Qa-max ≈ 137.3 kWh∙m−3, respectively. As a result, the energy efficiency turns out to be the same, equal to 70.3% for all three such modes investigated.
  • When, after dispatcher intervention, the temperature tm1 is reduced from 130 to 110 °C, the energies Qw and Qa reach maximum values at the moment τ2 = τsteam only at Δτa = 3 h and Δτb = 7 h. Then, they are equal to about Qw-max ≈ 93.9 kWh∙m−3 and Qa-max ≈ 129.5 kWh∙m−3, respectively, resulting in η ≈ 72.5%. In the case when Δτc = 11 h, the maximum values of Qw-max = 94.6 kWh∙m−3 and Qa-max = 131.2 kWh∙m−3 are reached at the moment of application of the dispatcher intervention, and this causes a reduction of η to η = 72.0%. In this case τsteam = 16.7 h.
  • When, after dispatcher intervention, tm1 is reduced from 130 to 100 °C, Qw and Qa reach maximum values at the moment τ2 = τsteam also only at Δτa = 3 h and Δτb = 7 h. Then, they are equal to approximately Qw-max ≈ 89.6 kWh∙m−3 and Qa-max ≈ 120.1 kWh∙m−3, respectively, resulting in η ≈ 74.6%. In the case when Δτc = 11 h, the maximum values of Qw-max = 90.6 kWh∙m−3 and Qa-max = 131.8 kWh∙m−3 are reached at the time of application of the dispatcher intervention and this causes a reduction of η to η = 68.7%.
The presented methodology can be applied in the creation of system software for model-based energy-efficient automatic control of technologies for the water vapor treatment of frozen and non-frozen wood materials with a desired duration of modes, set by a dispatcher.

Author Contributions

Conceptualization, N.D., P.N., and D.A.; methodology, N.D., P.N., and D.A.; software, N.D.; validation, N.D., P.V., and N.T.; formal analysis, P.N. and N.D.; investigation, N.D., P.N., P.V., and D.A.; resources, N.D., P.N., and D.A.; data curation, N.D. and N.T.; writing—original draft preparation, N.D. and D.A.; writing—review and editing, N.D., P.N., and D.A.; visualization, N.T., P.V., and N.D.; supervision, N.D.; project administration, N.D.; funding acquisition, N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Symbols
bwidth of the wooden prisms, m
cspecific heat capacity, J·kg−1·K−1
dthickness of the prisms, m
Ddiameter of the steaming autoclave, m
llength of the prisms, m
Llength of the autoclave, m
qthermal power, kW
Qthermal energy, kWh·m−3
Saria, m2
Ttemperature, K
ttemperature, °C: t = T − 273.15
umoisture content, kg·kg−1 = %/100
xcoordinate along d
ycoordinate along b
αconvective heat transfer coefficient, W·m−2·K−1
γloading of the autoclave, m3·m−3 = %/100
ηenergy efficiency, %
λthermal conductivity, W·m−1·K−1
ρdensity, kg·m−3
τtime, s
Δτstep along τ, s
Subscripts:
aautoclave
adanatomical direction (for wood)
avgaverage
bbasic (for density or for steaming mode)
bwbound water
crcross sectional to the fibers
cwcondensed water (for autoclave)
eemission (for autoclave)
eff1effective (for c of wood with frozen bound water)
eff2effective (for c of wood with frozen free water)
eff3effective (for c of non-frozen wood)
frfrozen
fvfree volume (for autoclave)
fwfree water
iceice
ilinsulating layer
hheat
imesh point along x
jmesh point along y
mmedium
mbmetal body (for autoclave and trolleys in it for placing of wood materials)
nfrnon-frozen
ssurface
wwood
0initial
Superscripts:
ntime level: n = 0, 1, 2, 3, …, τend/Δτ
272.15 at 272.15 K, i.e., at −1 °C
293.15 at 293.15 K, i.e., at 20 °C

References

  1. Zhukov, E.V.; Onegin, V.I. Technology of Protective and Decorative Coatings of Wood and Wood Materials; Ecologia: Moscow, Russia, 1993; 304p. [Google Scholar]
  2. Rüdiger, A. Grundlagen des Möbel- und Innenausbaus; DRW-Verlag GmbH & Co. KG: Leinfelden-Echterdingen, Germany, 1995; 306p. [Google Scholar]
  3. Jaić, M.; Živanović, R. Surface Processing of Wood—Theoretical Base and Technological Processes; University in Belgrade: Belgrade, Serbia, 2000; 400p. [Google Scholar]
  4. Kavalov, A.; Angelski, D. Technology of Furniture; University of Forestry: Sofia, Bulgaria, 2014; 390p. [Google Scholar]
  5. Vorreiter, L. Holztechnologisches Handbuch, Band I; Verlag Georg Fromme & Co.: Wien, Austria, 1949; 547p. [Google Scholar]
  6. Chudinov, B.S. Theoretical Research of Thermo-Physical Properties and Thermal Treatment of Wood. Ph.D. Thesis, SibLTI, Krasnoyarsk, Russia, 1966. [Google Scholar]
  7. Kollmann, F.F.; Côté, W.A., Jr. Solid wood. In Principles of Wood Science and Technology; Springer: New York, NY, USA; Berlin/Heidelberg, Germany, 1984; 592p. [Google Scholar]
  8. Shubin, G.S. Drying and Thermal Treatment of Wood; Lesnaya Promyshlennost: Moscow, Russia, 1990; 337p. [Google Scholar]
  9. Sohor, M.; Kadlec, P. Hydrothermal Treatment of Wood for Production of Veneer. Drevo 1990, 2. [Google Scholar]
  10. Lawniczak, M. Hydrothermal and Plasticizing Treatment of Wood. Part I. Boiling and Steaming of Wood; Agricultural Academy: Poznan, Poland, 1995; 149p. [Google Scholar]
  11. Trebula, P.; Klement, I. Drying and Hydrothermal Treatment of Wood; Technical University in Zvolen: Zvolen, Slovakia, 2002; 449p. [Google Scholar]
  12. Videlov, C.H. Drying and Thermal Treatment of Wood; University of Forestry: Sofia, Bulgaria, 2003; 335p. [Google Scholar]
  13. Pervan, S. Technology for Treatment of Wood with Water Steam; University in Zagreb: Zagreb, Croatia, 2009. [Google Scholar]
  14. Deliiski, N.; Dzurenda, L. Modelling of the Thermal Processes in the Technologies for Wood Thermal Treatment; Technical University in Zvolen: Zvolen, Slovakia, 2010; 224p. [Google Scholar]
  15. Niemz, P.; Sonderegger, W. Holzphysik: Physik des Holzes und der Holzwerkstoffe; Carl Hanser Verlag GmbH & Company KG: Munich, Germany, 2017; 580p. [Google Scholar]
  16. Deliiski, N.; Dzurenda, L.; Angelski, D.; Tumbarkova, N. Influence of selected factors on the duration and energy efficiency of autoclave steaming regimes of non-frozen prisms for veneer production. Energies 2021, 14, 7433. [Google Scholar] [CrossRef]
  17. Burtin, P.; Jay-Allemand, C.; Charpentier, J.P.; Janin, G. Wood Colour and Phenolic Composition under Various Steaming Conditions. Holzforschung 2000, 54, 33–38. [Google Scholar] [CrossRef]
  18. Riehl, T.; Welling, J.; Frühwald, A. Druckdämpfen von Schnittholz; Arbeitsbericht 2002/01; Institut für Holzphysik, Bundesforschungsanstalt für Forst- und Holzwirtschaft: Hamburg, Germany, 2002. [Google Scholar]
  19. Bekhta, P.; Niemz, P. Effect of High Temperature on the Change in Color, Dimensional Stability and Mechanical Properties of Spruce Wood. Holzforschung 2003, 57, 539–546. [Google Scholar] [CrossRef]
  20. Deliiski, N. Modelling and Technologies for Steaming Wood Materials in Autoclaves. Ph.D. Thesis, University of Forestry, Sofia, Bulgaria, 2003; 358p. [Google Scholar]
  21. Deliiski, N. Modelling and Automatic Control of Heat Energy Consumption Required for Thermal Treatment of Logs. Drvna Ind. 2004, 55, 181–199. [Google Scholar]
  22. Deliiski, N.; Sokolovski, S. Autoclaves for Intensive Resource Saving Steaming of Wood Materials. In Proceedings of the 2nd International Scientific Conference “Woodworking Techniques”, Zalesina, Croatia, 11–15 September 2007; pp. 19–26. [Google Scholar]
  23. Sokolovski, S.; Deliiski, N.; Dzurenda, L. Constructive Dimensioning of Autoclaves for Treatment of Wood Materials under Pressure. In Proceedings of the 2nd International Scientific Conference “Woodworking techniques”, Zalesina, Croatia, 11−15 September 2007; pp. 117–126. [Google Scholar]
  24. Dagbro, O.; Torniainen, P.; Karlsson, O.; Morén, T. Colour Responses from Wood, Thermally Modified in Superheated Steam and Pressurized Steam Atmospheres. Wood Mater. Sci. Eng. 2010, 5, 211–219. [Google Scholar] [CrossRef]
  25. Deliiski, N. Model Based Automatic Control of the Wood Steaming Process in Autoclaves. In Proceedings of the 4th International Science Conference “Woodworking Techniques”, Prague, Czech Republic, 7−10 September 2011; pp. 67–72. [Google Scholar]
  26. Deliiski, N. Transient Heat Conduction in Capillary Porous Bodies. In Convection and Conduction Heat Transfer; Ahsan, A., Ed.; InTech Publishing House: Rijeka, Croatia, 2011; pp. 149–176. [Google Scholar]
  27. Deliiski, N.; Dzurenda, L.; Trichkov, N.; Tumbarkova, N. Computing the 2D Temperature Field in Non-Frozen Logs at Changing Atmospheric Temperature and during their Subsequent Autoclave Steaming. Acta Fac. Xylologiae Zvolen 2020, 62, 47–59. [Google Scholar]
  28. Chudinov, B.S. Theory of the Thermal Treatment of Wood; Nauka: Moscow, Russia, 1968; 255p. [Google Scholar]
  29. Steinhagen, H.P. Heat Transfer Computation for a Long, Frozen Log Heated in Agitated Water or Steam—A Practical Recipe. Holz Roh Werkstoff 1991, 49, 287–290. [Google Scholar] [CrossRef]
  30. Dzurenda, L.; Deliiski, N. Mathematical Model for Calculation Standard Values for Heat Energy Consumption during the Plasticization Process of Wood Logs and Prisms by Hot Water in Pits. Acta Fac. Xylologie Zvolen 2011, 53, 25–36. [Google Scholar]
  31. Shirazinia, M.; Moya, R.; Muñoz, F. Properties of Laminated Curves Manufactured with Steamed Veneers from Fast-Growth Tropical Wood in Costa Rica. Madera Bosques 2011, 17, 85–101. [Google Scholar] [CrossRef] [Green Version]
  32. Deliiski, N. Modelling of the Energy Needed for Heating of Capillary Porous Bodies in Frozen and Non-Frozen States; Lambert Academic Publishing, Scholars’ Press: Saarbrücken, Germany, 2013; 106p. [Google Scholar]
  33. Šprdlík, V.; Brabec, M.; Mihailović, S.; Rademacher, P. Plasticity Increase of Beech Veneer by Steaming and Gaseous Ammonia Treatment. Maderas Cienc. Tecnol. 2016, 1, 91–98. [Google Scholar] [CrossRef] [Green Version]
  34. Deliiski, N.; Dzurenda, L.; Tumbarkova, N.; Angelski, D. Mathematical Description of the Latent Heat of Bound Water in Wood during Freezing and Defrosting. Acta Fac. Xyilologiae Zvolen 2020, 62, 41–53. [Google Scholar]
  35. Konopka, A.; Chuchala, D.; Orlowski, K.A.; Vilkovská, T.; Klement, I. The Effect of Beech Wood (Fagus sylvatica L.) Steaming Process on the Colour Change versus Depth of Tested Wood Layer. Wood Mater. Sci. Eng. 2021, 17, 420–428. [Google Scholar] [CrossRef]
  36. Moya, R.; Tenorio, C.; Torres, J.D.C. Steaming and Heating Dipteryx panamensis Logs from Fast-Grown Plantation: Reduction of Growth Strain and Effects on Quality. For. Prod. J. 2021, 71, 3–10. [Google Scholar] [CrossRef]
  37. Deliiski, N.; Vitchev, P.; Angelski, D.; Tumbarkova, N. Change in the duration of the autoclave steaming regimes of non-frozen prisms at dispatching interferences in the production of veneer. Innov. Woodwork. Ind. Eng. Des. 2022, 11, 25–32. [Google Scholar]
  38. Deliiski, N.; Dzurenda, L.; Vitchev, P.; Tumbarkova, N.; Angelski, D. Change of the energy for warming up of wooden prisms during autoclave steaming at dispatching interferences in production of veneer. Chip- Chipless Woodwork. Processes 2022, 12, 123–128. [Google Scholar]
  39. Deliiski, N.; Niemz, P.; Vitchev, P.; Angelski, D.; Tumbarkova, N. Computing the duration of regimes for autoclave steaming of frozen wooden prisms under variable operating conditions in veneer production. Wood Mater. Sci. Eng. 2022, 17, 421–428. [Google Scholar] [CrossRef]
  40. Hadjiski, M.; Deliiski, N.; Grancharova, A. Spatiotemporal Parameter Estimation of Thermal Treatment Process via Initial Condition Reconstruction using Neural Networks. In Intuitionistic Fuzziness and Other Intelligent Theories and Their Applications; Hadjiski, M., Atanasov, K.T., Eds.; Springer: Cham, Switzerland, 2019; pp. 51–80. [Google Scholar]
  41. Hadjiski, M.; Deliiski, N. Advanced Control of the Wood Thermal Treatment Processing. Cybern. Inf. Technol. 2016, 16, 179–197. [Google Scholar] [CrossRef] [Green Version]
  42. Hadjiski, M.; Deliiski, N.; Tumbarkova, N. Intelligent Hybrid Control of Thermal Treatment Processes of Wood. In Proceedings of the IEEE 10th International Conference on Intelligent Systems, IS 2020, Varna, Bulgaria, 28–30 August 2020; pp. 482–489. [Google Scholar]
  43. Sgurev, V.; Hadjiski, M.; Deliiski, N. Multicriteria Optimal Control of Industrial Thermal Processes with Distributed Parameters under Variable Operational Conditions. In Complex Systems: Spanning Control and Computational Cybernetics: Applications. Studies in Systems, Decision and Control; Shi, P., Stefanovski, J., Kacprzyk, J., Eds.; Springer: Cham, Switzerland, 2022; Volume 415, pp. 333–356. [Google Scholar]
  44. Deliiski, N.; Dzurenda, L.; Angelski, D.; Tumbarkova, N. An Approach to Computing Regimes for Autoclave Steaming of Prisms for Veneer Production with a Limited Power of the Heat Generator. Acta Fac. Xylologiae Zvolen 2018, 60, 101–112. [Google Scholar]
  45. Deliiski, N. Mathematische Beschreibung der spezifischen Wärmekapazität des aufgetauten und gefrorenen Holzes. In Proceedings of the VIII International Symposium on Fundamental Research of Wood, Warsaw, Poland, 8–12 October 1990; pp. 229–233. [Google Scholar]
  46. Deliiski, N. Mathematical Description of the Thermal Conductivity Coefficient of Non-frozen and Frozen Wood. In Proceedings of the 2nd International Symposium on Wood Structure and Properties ’94, Zvolen, Slovakia, 5−9 September 1994; pp. 127–134. [Google Scholar]
  47. Hrčka, R.; Babiak, M. Wood Thermal Properties. In Wood in Civil Engineering; Consu, G., Ed.; InTechOpen: London, UK, 2017; pp. 25–43. [Google Scholar]
  48. Kanter, K.R. Investigation of the Thermal Properties of Wood. Ph.D. Thesis, MLTI, Moscow, Russia, 1955. [Google Scholar]
  49. Požgaj, A.; Chovanec, D.; Kurjatko, S.; Babiak, M. Structure and Properties of Wood, 2nd ed.; Priroda a.s.: Bratislava, Slovakia, 1997; 486p. [Google Scholar]
  50. Khattabi, A.; Steinhagen, H.P. Numerical Solution to Two-dimensional Heating of Logs. Holz Roh Werkstoff 1992, 50, 308–312. [Google Scholar] [CrossRef]
  51. Khattabi, A.; Steinhagen, H.P. Analysis of Transient Non-linear Heat Conduction in Wood Using Finite-difference Solutions. Holz Roh Werkstoff 1993, 51, 272–278. [Google Scholar] [CrossRef]
  52. Khattabi, A.; Steinhagen, H.P. Update of “Numerical Solution to Two-dimensional Heating of Logs”. Holz Roh Werkstoff 1995, 53, 93–94. [Google Scholar] [CrossRef]
  53. Steinhagen, H.P. Thermal Conductive Properties of Wood, Green or Dry, from –40 to +100 °C: A Literature Review; General Technical Report FPL-09; Forest Products Laboratory, Forest Service: New York, NY, USA, 1977. [Google Scholar]
  54. Steinhagen, H.P. Computerized Finite-difference Method to Calculate Transient Heat Conduction with Thawing. Wood Fiber Sci. 1986, 18, 460–467. [Google Scholar]
  55. Steinhagen, H.P.; Lee, H.W. Enthalpy Method to Compute Radial Heating and Thawing of Logs. Wood Fiber Sci. 1988, 20, 415–421. [Google Scholar]
  56. Steinhagen, H.P.; Lee, H.W.; Loehnertz, S.P. LOGHEAT: A Computer Program of Determining Log Heating Times for Frozen and Non-Frozen Logs. For. Prod. J. 1987, 37, 60–64. [Google Scholar]
  57. Deliiski, N.; Dzurenda, L.; Angelski, D.; Tumbarkova, N. Computing the Energy for Warming up of Prisms for Veneer Production during Autoclave Steaming with a Limited Power of the Heat Generator. Acta Fac. Xilologiae Zvolen 2019, 61, 63–74. [Google Scholar]
  58. Hadjiski, M.; Deliiski, N.; Angelski, D. Computing the Processing Medium Temperature and Heat Fluxes in the Beginning of Regimes for Autoclave Steaming of Frozen Wood Materials. In Proceedings of the International Conference Automatics and Informatics (ICAI), Varna, Bulgaria, 30 September–2 October 2021; pp. 393–397. [Google Scholar]
  59. Mörath, M. Das Dämpfen und Kochen in der Furnier– und Sperrholzindustrie. Holztechnik 1949, 29, 129–134. [Google Scholar]
Figure 1. Stages of the wood defrosting process and thermo-physical characteristics of wood used in them.
Figure 1. Stages of the wood defrosting process and thermo-physical characteristics of wood used in them.
Processes 11 00822 g001
Figure 2. Change of Tm during steaming modes in absence and presence of dispatcher intervention.
Figure 2. Change of Tm during steaming modes in absence and presence of dispatcher intervention.
Processes 11 00822 g002
Figure 3. Changes in ts, tw-avg, t1, and t2 of the prisms during steaming Mode 0, Mode 5, and Mode 8, depending on τ.
Figure 3. Changes in ts, tw-avg, t1, and t2 of the prisms during steaming Mode 0, Mode 5, and Mode 8, depending on τ.
Processes 11 00822 g003
Figure 4. Changes in tm, Qa, and Qw during steaming modes 0, 1, 2, and 3, depending on τ.
Figure 4. Changes in tm, Qa, and Qw during steaming modes 0, 1, 2, and 3, depending on τ.
Processes 11 00822 g004
Figure 5. Changes in tm, Qa, and Qw during steaming modes 0, 4, 5, and 6, depending on τ.
Figure 5. Changes in tm, Qa, and Qw during steaming modes 0, 4, 5, and 6, depending on τ.
Processes 11 00822 g005
Figure 6. Changes in tm, Qa, and Qw during steaming modes 0, 7, 8, and 9, depending on τ.
Figure 6. Changes in tm, Qa, and Qw during steaming modes 0, 7, 8, and 9, depending on τ.
Processes 11 00822 g006
Figure 7. Change in Qw-max (a) and Qa-max (b) of the studied steaming modes, depending on tm1 and τ.
Figure 7. Change in Qw-max (a) and Qa-max (b) of the studied steaming modes, depending on tm1 and τ.
Processes 11 00822 g007
Figure 8. Change in η of the studied modes, depending on tm1 and τ.
Figure 8. Change in η of the studied modes, depending on tm1 and τ.
Processes 11 00822 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deliiski, N.; Niemz, P.; Angelski, D.; Vitchev, P.; Tumbarkova, N. Computing the Thermal Efficiency of Autoclaves during Steaming of Frozen Prisms for Veneer Production at Changing Operational Conditions. Processes 2023, 11, 822. https://doi.org/10.3390/pr11030822

AMA Style

Deliiski N, Niemz P, Angelski D, Vitchev P, Tumbarkova N. Computing the Thermal Efficiency of Autoclaves during Steaming of Frozen Prisms for Veneer Production at Changing Operational Conditions. Processes. 2023; 11(3):822. https://doi.org/10.3390/pr11030822

Chicago/Turabian Style

Deliiski, Nencho, Peter Niemz, Dimitar Angelski, Pavlin Vitchev, and Natalia Tumbarkova. 2023. "Computing the Thermal Efficiency of Autoclaves during Steaming of Frozen Prisms for Veneer Production at Changing Operational Conditions" Processes 11, no. 3: 822. https://doi.org/10.3390/pr11030822

APA Style

Deliiski, N., Niemz, P., Angelski, D., Vitchev, P., & Tumbarkova, N. (2023). Computing the Thermal Efficiency of Autoclaves during Steaming of Frozen Prisms for Veneer Production at Changing Operational Conditions. Processes, 11(3), 822. https://doi.org/10.3390/pr11030822

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop