ESTABLISHING REGULARITIES IN THE PROPAGATION OF PHASE TRANSFORMATION FRONT DURING TIMBER THERMAL MODIFICATION Текст научной статьи по специальности «Физика»

The creation of environmentally friendly protective materials for building structures made of wood could make it possible to influence the processes of stability and the physical-chemical properties at the thermal modification of hornbeam wood over a certain time. That necessitates studying the conditions for investigating phase transformations when the timber is exposed to high temperature, as well as establishing the mechanism of hornbeam wood thermal modification. Given this, a mathematical model of the phase transformation process during the transfer of heat flux to a sample was built. Based on the derived dependences, it was established that when hornbeam wood is exposed to temperature treatment, it undergoes endothermic phase transformations characterized by the heat absorption and change in the color of hornbeam wood. In particular, at a temperature of200 °C, the temperature in the wood decreases by 5 % due to the chemical changes in the structure of cell wall components (lignin, cellulose, and hemi-cellulose). It was found that the process of thermal modification is accompanied by the decomposition of hemicellulose and the amorphous part of cellulose, a decrease in moisture absorption, as well as a decrease in the volume of substances that are a medium for the development of fungi. In addition, lignin and the resulting pseudo lig-nin undergo a process of polymerization and redistribution throughout the cell volume. At the same time, they give the cell walls higher density, hardness, increase hydrophobicity (water repel-lency), thereby reducing the ability to absorb moisture and swell. It was established that the most effective parameter of phase transformations is the temperature and aging duration. The results of moisture absorption have been given; it has been found that over 6 hours of modified tim -ber exposure, its moisture absorption decreases by more than 10 times, which allows its application at facilities with high humidity

Keywords: thermally modified timber, modification efficiency, moisture absorption, diffusion, timber moisture resistance

UDC 614.842.5:349.211

pOI: 10.15587/1729-4061.2021.225310|

ESTABLISHING REGULARITIES IN THE PROPAGATION OF PHASE TRANSFORMATION FRONT DURING TIMBER THERMAL MODIFICATION

Yu. Tsapko

Doctor of Technical Sciences** V. D. Glukhovsky Scientific Research Institute for Binders and Materials Kyiv National University of Construction and Architecture Povitroflotsky ave., 31, Kyiv, Ukraine, 03037 E-mail: [email protected] O. Horbachova PhD*

E-mail: [email protected] А. Tsapko

Senior Research Fellow Department of Research on Quality and Storage Conditions of Petroleum Products and an Industrial Group of Goods Ukrainian State Research Institute "Resurs" Kasimira Malevicha str., 84, Kyiv, Ukraine, 03150 E-mail: [email protected] S. Mazurchuk PhD*

E-mail: [email protected] D. Zavialov Assistant* E-mail: [email protected] N. Buiskykh PhD*

E-mail: [email protected] *Department of Technology and Design of Wood Products** **National University of Life and Environmental Sciences of Ukraine Heroiv Oborony str., 15, Kyiv, Ukraine, 03041

Received date 21.12.2020 Accepted date 02.02.2021 Pu blished date 23.02.2021

Copyright © 2021, Yu. Tsapko, O. Horbachova, A. Tsapko, S. Mazurchuk, D. Zavialov, N. Buiskykh

This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0)

1. Introduction

Representing a building material, timber is widely used in construction and architecture due to its mechanical and operational properties; however, it is destroyed under the influence of atmospheric factors. Improving the operational level of facilities where timber structures are used is possi-

ble through thermal modification. That implies rendering the timber the capability to resist, for example, the effect of moisture, and prevent the destruction of timber thereby decelerating the process of destruction [1].

Given the above, the thermal modification of timber is associated with difficulties related to the use of technological modes such as time, temperature. This is primarily

because the structure and composition of wood differ so that the modification process is not completed while the use of low-heat modified timber leads to destruction. The knowledge of the physical-chemical properties of such substances, their quality indicators, the mechanism that affects materials, would make it possible to select them considering economic indicators, the duration and safety of application, environmental aspects, etc. [2, 3].

Therefore, the development of technological modes of timber modification, studying the thermal-physical transformations, the effect of the structure's features on this process, is an unresolved component of the production of durable building materials from wood. This predetermines the need to establish a mechanism of thermal modification.

2. Literature review and problem statement

In work [4], it was shown that the thermal modification of timber leads to changes in weight, wettability, and color, thereby increasing its resistance against destruction due to the chemical transformations of timber components. However, issues related to quality control, modeling, and studying the reasons for improving the properties of modified timber remain unresolved. This very approach was used in work [5] where it is shown that one of the best materials for cladding building structures is thermally modified timber because the change in temperature and humidity fields does not affect the intensive destruction. The use of properly selected and correct modification technology could slow down the processes of destruction and environmental consequences. However, the extent and modes of timber modification were not specified.

Extensive use of thermally modified timber [6] has led to the need to create a reliable quality control, which includes control over product deviations within certain limits. That allows a third-party control in the case of certification and regulation of consumer complaints and requirements. However, it is not specified what methods are needed to improve the targeted properties of modified timber during industrial production.

In particular, the authors of [7] devised a new technology of heat treatment in the presence of steam; it is, therefore, a standard hygrothermal treatment. When choosing wood species, one needs to deal with factors such as timber oxida-tive degradation. An option to overcome the corresponding difficulties may be to determine the reactions that occur through the presence of moisture.

A variant to overcome the relevant problems might be to study changes in the swelling and surface roughness of alder and elmwood after heat treatment at two different temperatures and durations [8]. The results showed that the parameters of swelling and surface roughness differed significantly for two temperatures and two durations of heat treatment. The reason for this may be those objective difficulties that are associated with the value of swelling and surface roughness, which decreased with an increase in the temperature and time of treatment.

A change in the coloration and reflectivity of timber surfaces, caused by artificial weathering obtained in a solar box chamber, which simulates external conditions and subsequent leaching of water, was assessed in [9]. As the weathering time increases, the untreated surfaces of timber samples get darker while the modified samples become

lighter. However, there are unresolved issues related to the possibility of having similar color or, in any case, reducing the starting chromatic difference observed at the beginning of tests for weather-related stability.

Many chemical reactions occur during heat treatment, leading to changes in the components of the primary cell wall of wood and the darkening of the material [10]. Other changes include the resistance of modified timber against fungal decay, making it suitable for use indoors and outdoors as cladding, flooring, floors, garden furniture, and window frames. However, the effectiveness of application was not determined for the most common types of timber, in particular, pine.

Laboratory tests have shown a positive effect of thermal modification on the durability, stability of size, as well as thermal conductivity of timber [11]. The monitoring results showed that the elements and windows made of thermally modified spruce have significantly lower moisture content in timber compared to windows made of unmodified timber. However, the effectiveness of application was not determined when compared with wax and other water repellents.

Article [12] reports a study into the use of artificial aging of timber, which plays an important role in assessing the results of work by reducing time compared to natural weathering. The approach is to protect the surface by using various types of commercially-available agents such as solvents contained in water, with large content of solids, powder coatings. However, the effectiveness of thermal modification of timber, its resistance to artificial aging is not shown.

One of the approaches to improve the durability of timber is a set of processes that provide the treated material with a better capability to cope with the damage caused by the external environment, by increasing the duration of treatment [13]. It was established that this process is also performed to enhance the physical, mechanical, or aesthetic properties of timber, and makes it possible to obtain products that are not harmful to users and the environment, similar to natural wood [14]. However, the mechanism of timber degradation caused by operation and its impact on destruction rates were not specified.

Thus, the scientific literature has revealed that the thermal modification of timber could reduce the destruction of a building structure and expand the scope of timber application. The above allows us to assert that it is expedient to conduct a study to determine the parameters that ensure the resistance against destruction, as well as define the mechanism affecting timber transformation exposed to thermal modification. Establishing the operational mechanism involved in the thermal modification of timber to render its resistance to moisture absorption has predetermined our research in this area.

3. The aim and objectives of the study

This work aims to establish the parameters of phase transformations during timber thermal modification. That could improve the technology of manufacturing such products and expand the scope of their application.

To achieve the set aim, the following tasks have been solved:

- to model the process of phase transformations' front propagation when timber is exposed to thermal influence due to high temperature;

- to establish features in the decrease of moisture absorption in order to confirm the termination of phase transformations during timber thermal modification.

tain time, the sample was weighed on scales; the amount of water absorbed was determined.

Based on the measured values, we registered the changes and determined the efficiency of timber thermal modification.

4. Materials and methods to study the process of timber thermal modification

4. 1. The examined materials and equipment used in the experiment

Hornbeam timber, which belongs to hardwood, is used for the manufacture of floors in rooms with high humidity (swimming pools, saunas); to reduce the rate of water absorption, it is advisable to conduct thermal modification.

Our research involved samples of raw hornbeam wood of 25x20x20 mm (Fig. 1).

The thermal modification of timber was conducted at a temperature of 200 °C over 1-6 hours with a 1-hour interval.

Fig. 1. A model sample of timber

To establish the degree of timber water absorption, the samples of thermally modified timber were used.

4. 2. Procedure for establishing the indicators of samples' properties

Since the phase transformations during timber thermal modification occur at the intercellular level, it is impossible to install a measuring device inside the timber to measure the temperature without damaging its structure; and the destruction of timber could lead to erroneous results. Our study, aimed at modeling the process of propagation of the phase transformations' front during timber thermal modification, was performed using the basic provisions from mathematical physics [15].

We determined the degree of moisture absorption by timber according to the working procedure whose essence was to experimentally establish the amount of the water absorbed by a sample at its exposure to 100 % moisture. To find the values of the amount of water absorbed by timber, specialized equipment was designed and manufactured (Fig. 2).

Fig. 2. Device for testing water absorption by timber

A tested sample was fixed in a special cuvette so that it was in a humid environment above the water. After a cer-

5. Modelling the process of timber thermal conductivity at its thermal modification

As a result of the thermal treatment of timber, under the action of heat flow, the direction of the decomposition of hemicellulose, the redistribution of lignin, etc. changes, towards the formation of volatile products and carbonized residue. This process is characterized by heat absorption during phase transformations and discoloration of hornbeam wood. Considering the above, the issue arises about studying the phase transformations of timber exposed to heat.

It should be noted that determining the thermal-physical characteristics of thermally modified timber is associated with some obstacles, namely the measurement of temperature in a timber layer that changes over time.

To establish the temperature of phase transformations within thermally modified timber, a method has been proposed to solve the problem of thermal conductivity for a plate with the thermal-physical properties that depend on temperature. Over the initial time, a time-constant heat flux qo is applied to the surface of a timber sample. That is, it is instantly heated to a temperature that is kept constant throughout the heating process; the temperature distribution passes to the center of the sample (Fig. 3).

Three regions were considered (Fig. 3):

- 1 - external environment, x<0;

- 2 - phase transition zone, 0<x<Z (Z is the hemicellu-lose conversion coordinate, m);

- 3 - timber (material of the sample of hard substance) (R-Z), m.

Front of phase transformations

T

Fig. 3. The scheme of the process of propagation of the phase

transformations' front during timber thermal modification: 1 — external environment; 2 — phase transition zone; 3 — timber

The differential equation of heat transfer in timber takes the following form:

T

d2T(x, t) 1 dT(x, t)

—^^---t-lJ- = 0, (x>Z(t); t>0),

dx a 9t

at the initial and boundary conditions

T ( x, £ )[o= 0, x ^ 0;

T M = z(t)=*(i)- f ^ °

T(x,£)| ^ x > 0,

(1)

(2)

(3)

(4)

r ( * ) = ^J

1

it-t

(x-Zt)

g ^ dT.

The unknown function y(x) can be found from a condition at the boundary x = Z ■ t:

dT.

However, the solution to equation (5) is cumbersome, so, by using a Laplace transform from [16], we can write down:

I— t l \ (x-Zt)2

t (x, p)=r e-p . e-i^M dT. dt=

Jt-T

(x-ZT)2

it-T

(7)

The original of this equation at any point x at moment t is:

T ( X ,0) = -M X-Z(/) ■

1 ; 2s[an I (t -t)3/2

(x-Z(t ))2 4a(t-T)

■ <(t )^t.

The solution to equation (8), taking into account the initial and boundary conditions in the region of L-transforms, takes the following form:

T( x,, } [43 - )

-fa

S

P

(13)

4a

Considering (11) to (13), let us write down equation (9) in the following form:

T (p) = T S z . e~Ta^.

4a

pVP--

where a is the coefficients of temperature timber conductivity, m2-s-1; Z(t) is the coordinate of timber phase transformations; 9(t) is the rate of timber phase transition changes; T is the temperature, °C; x is a coordinate, m.

A solution to the problem can be represented in the following form:

.1 Tl ■ Z ■ 2 p 4a

1

-X 4p

V"

4P -

(14)

4a

Using a ratio from [15], we can write

•Ja

(5)

(6)

4P\4P -

_Z_

4a

4a

xZ ( z Ja

■ erfc

x

4a - Z Jt l4~t 4a

(15)

x f-

-~T VP Va

4P [4p-

4a

>erfc

4a

2 4t

- - 4-r

a a J

■erfc

4a__^4t

24t 4a

(16)

Considering (11), (12), let us write down the solution to the problem in a time domain:

T ( x, t ) = To

3 -(X-Zi) [ X - Zt \ 1 r x —e a erfc \-— erfc—¡=

2 J , l40t J 2 J 240t

. (17)

(8)

At Z=0, when a timber phase transition change is finished, the solution will take the following form:

T ( x,t ) = To- erfc

24at

(18)

(9)

Thus, the argument of the action function 9 depends not on the parameter p but the difference:

p-4- 4P. (10)

Va

Since 9(t) is given by the following function:

<l>(i ) = V z (t ). (11)

Then

<(p) = \ (12) p

The derived equation allows calculating the temperature of phase transformations during timber thermal modification based on the experimental values of temperature, the geometric dimensions of a timber sample, and timber thermal conductivity.

Fig. 4 shows the calculation of temperature in a timber sample depending on the time of phase transformations, defined from equation (18) and by applying data from [15], namely the value of the thermal conductivity coefficient of timber at a temperature of about 200 °C, which is 1810-6 m2/s.

When timber is exposed to thermal treatment, it undergoes endothermic phase transformations, which are characterized by heat absorption and timber discoloration. In particular, at 200 °C, the temperature in the timber, due to the chemical changes in the structure of the components of the cell wall (lignin, cellulose, and hemicellulose), was 200 °C; according to Fig. 4, it was reduced, due to phase transformations, by 3 %.

o

196

194

3

13 1-

u a

s

£ 192

As can be seen from Fig. 6, the largest amount of absorbed moisture was recorded for the untreated hornbeam timber. Thermal modification gradually reduces the ability of the timber to absorb water.

190

0 2 4 6

Time of phase transformations, hours

Fig. 4. The dependence of a timber sample's temperature on the time of phase transformations (points correspond to the values calculated from equation (18), averaged by a trend's line)

Thus, in the process of thermal modification, hemicellulo-ses and the amorphous part of cellulose decompose, moisture absorption decreases, and the number of substances that are the medium for the development of fungi decreases. In addition, lignin and the formed pseudo lignin undergo a process of polymerization and redistribution throughout the volume of the cell, which gives the cell walls greater density, hardness. As a result, hydrophobicity (water repellency) increases, and the capacity to absorb moisture and swell decreases. Polymerized lignin fills the inner cavity of the cell thereby forming a closed porous structure with a low capability to bind water.

It is established ( Fig. 4) that the most effective parameter of phase transformations is temperature and exposure time.

Thus, we have derived the dependence that is necessary to predict the end time of phase transitions during timber thermal modification, which makes it possible to directly calculate the movement of timber phase transformations depending on the time of temperature impact.

6. Results of studying the features related to the decreased water absorption by the thermally modified timber

To establish the thermal-physical characteristics of thermally modified timber, we performed a study into the modification under the impact of a heating device at 200 °C over a certain time of modification (Fig. 5).

It should be noted that thermal modification should be performed at temperatures above 200 °C; however, it is necessary to take into consideration the ignition temperature of timber, which is from 215 °C to 225 °C, or carry out the process in a protective environment. Therefore, in a given case, taking into consideration safety requirements, the thermal modification was conducted at 200 °C.

As can be seen from Fig. 5, during the modification the hornbeam timber changed color from light to dark depending on the exposure time, respectively, and underwent phase transformations of the components and, accordingly, acquired other properties, in particular moisture absorption.

The results of moisture absorption by timber are given in Table 1.

Our results (Table 1) show that when the timber is thermally modified for 6 hours, the moisture absorption is reduced by more than 2.4 times, which allows it to be used at facilities with high humidity.

Fig. 6 shows the moisture absorption of hornbeam timber depending on the exposure time in a humid environment.

Fig. 5. Model samples of the hornbeam timber thermally modified at a temperature of 200 °C during: a — 1 hour; b — 2 hours; c — 3 hours; d — 4 hours; e — 5 hours; f — 6 hours

Table 1

Results of the experimentally determined timber water absorption

Sample No. Gain in the absorbed water, g

Time of exposure, day

0 1 2 6 8 12 21 26

0 6.80 8.1 8.35 8.62 8.73 8.75 8.76 8.76

1 7.55 8.03 8.31 8.59 8.57 8.63 8.70 8.72

2 7.34 7.73 7.99 8.3 8.38 8.33 8.38 8.48

3 7.38 7.67 7.97 8.24 8.30 8.34 8.44 8.44

4 7.11 7.32 7.69 7.9 8.00 8.07 8.11 8.12

5 7.26 7.46 7.74 7.89 7.93 7.99 7.99 8

6 7.19 7.39 7.6 7.8 7.79 7.86 7.86 7.86

7. Discussion of results of studying the process of propagation of the phase transformations' front at timber thermal modification

During the thermal modification of timber, as indicated by the research results, both theoretical equations (18) and Fig. 4 and the experimental study illustrated by Fig. 5, there is a natural process of the chemical transformations of hornbeam wood components. Namely, color changes under the impact of temperature and, accordingly, changes in structure, which, in turn, could lead to a decrease in moisture absorption. Phase transformations in the structure of wood under a temperature impact are characterized by an endothermic effect, which decreases over time. The thermally modified timber, at the same exposure temperature but over a longer time, is characterized by lower moisture absorption. Such a mechanism of thermal timber modification is probably a factor in regulating the extent of the formation of weather-resistant material. This agrees with the data known from [5, 8], the authors of which also link the effectiveness of protection against moisture during the thermal modification of hornbeam wood. In contrast to the results reported in [4, 7], our findings on the process of moisture absorption by thermally modified hornbeam timber and changes in the moisture insulating properties of hornbeam timber allow us to state the following:

a

e

- the main regulator of the process of moisture absorption is not only the formation of the coating layer but also the chemical transformations of the components of hornbeam wood, which, under the impact of temperature and humidity fields, provide for the resistance against moisture;

- a significant impact on the process of protection of thermally modified hornbeam timber from moisture fluctuations is exerted towards the formation of water-resistant capillary porous elements throughout the volume and at the surface of the natural material.

detection could make it possible to study the transformation of hornbeam timber, which moves in the direction of increased resistance to destruction, and to identify those variables that significantly affect the onset of the transformation of this process.

The results of the current study are limited to the use of the thermally modified hornbeam timber; it is necessary to refine the data when applying other types of wood.

9. Conclusions

e '3

O

= -0,0055x2 + 0,1824x + 7,5636

= -0,0032x2 + 0,1131x + 7,8628

= -0,0031x2 + 0,11x + 7,6016 -0,0031x2 + 0,1099x + 7,2737

= -0,003x2 + 0,1092x + 7,5781

= -0,0023x2 + 0,0811x + 7,4128 = -0,0022x2 + 0,0752x + 7,3289

6,5

10

15

20

Time of exposure, day

Fig. 6. The dependence of the amount of absorbed water on the exposure time of

hornbeam timber

Such conclusions can be considered reasonable from a practical point of view because they allow for a substantiated approach to determining the necessary technology for the thermal modification of hornbeam timber. From a theoretical point of view, they allow us to argue on determining the mechanism of phase transformation processes, which is a certain advantage of this study. The results of determining the moisture absorption during thermal modification (Table 1) indicate the ambiguous effect exerted on the nature of the change in humidity by the modified hornbeam timber. In particular, this implies the availability of data sufficient for the quality process of the inhibition of moisture diffusion, and the detection, on its basis, of the time when a decline in moisture resistance begins. This

1. We have modeled the process of propagation of the phase transformations' front, and established the dependence of temperature on the time of hornbeam timber exposure, as well as derived the estimation dependences that could help find a change in the temperature dynamics during phase transformations. At a temperature exposure, endothermic phase transformations take place in hornbeam timber, which are characterized by heat absorption and discoloration of hornbeam wood. In particular, at a temperature of 200 °C, the temperature in hornbeam timber, due to the chemical changes in the structure of the cell wall (lignin, cellulose, and hemi-cellulose), is reduced by up to 3 %.

2. Special features in the inhibition of the process of moisture advancement to the thermally modified hornbeam timber imply several aspects, namely, the formation of waterproof components, as well as capillary porous elements, which are characterized by the

formation, at the hornbeam timber surface, of a waterproof layer. Thermal modification of hornbeam timber reduces moisture absorption by more than 2.4 times within 6 hours, which allows it to be used at facilities with high humidity.

Acknowledgments

We acknowledge financial support to this work, performed within the funded budget No. 0118U002016, as well as for the development of scientific topics in the framework of the scientific cooperation program COST Action FP 1407 "Elucidating timber modification by using an integrated scientific and environmental approach" under the EU program HORIZON2020.

25

y

y

0

5

References

1. Tsapko, Y., Tsapko, A., Bondarenko, O. (2019). Effect of a flameretardant coating on the burning parameters of wood samples. Eastern-European Journal of Enterprise Technologies, 2 (10 (98)), 49-54. doi: http://doi.org/10.15587/1729-4061.2019.163591

2. Tsapko, Y., Lomaha, V., Bondarenko, O. P., Sukhanevych, M. (2020). Research of Mechanism of Fire Protection with Wood Lacquer. Materials Science Forum, 1006, 32-40. doi: http://doi.org/10.4028/www.scientific.net/msf.1006.32

3. Tsapko, Y., Lomaha, V., Tsapko, A., Mazurchuk, S., Horbachova, O., Zavialov, D. (2020). Determination of regularities of heat resistance under flame action on wood wall with fire-retardant varnish. Eastern-European Journal of Enterprise Technologies, 4 (10 (106)), 55-60. doi: http://doi.org/10.15587/1729-4061.2020.210009

4. Esteves, B. M., Pereira, H. M. (2008). Wood modification by heat treatment: A review. BioResources, 4 (1), 370-404. doi: http:// doi.org/10.15376/biores.4.1.370-404

5. Humar, M., Lesar, B., Krzisnik, D. (2020). Moisture Performance of Façade Elements Made of Thermally Modified Norway Spruce Wood. Forests, 11 (3), 348. doi: http://doi.org/10.3390/f11030348

6. Humar, M., Repic, R., Krzisnik, D., Lesar, B., Cerc Korosec, R., Brischke, C. et. al. (2020). Quality Control of Thermally Modified Timber Using Dynamic Vapor Sorption (DVS) Analysis. Forests, 11 (6), 666. doi: http://doi.org/10.3390/f11060666

7. Sandberg, D., Kutnar, A., Mantanis, G. (2017). Wood modification technologies - a review. iForest - Biogeosciences and Forestry, 10 (6), 895-908. doi: http://doi.org/10.3832/ifor2380-010

8. Aytin, A., Korkut, S. (2015). Effect of thermal treatment on the swelling and surface roughness of common alder and wych elm wood. Journal of Forestry Research, 27 (1), 225-229. doi: http://doi.org/10.1007/s11676-015-0136-7

9. Pelosi, C., Agresti, G., Lanteri, L., Picchio, R., Gennari, E., Lo Monaco, A. (2020). Artificial Weathering Effect on Surface of Heat-Treated Wood of Ayous (Triplochiton scleroxylon K. Shum). The 1st International Electronic Conference on Forests (IECF). Available at: https://www.researchgate.net/publication/345761222_Artificial_Weathering_Effect_on_Surface_of_Heat-Treated_ Wood_of_Ayous_Triplochiton_scleroxylon_K_Shum

10. Ugovsek, A., Subic, B., Rep, G., Humar, M., Lesar, B., Thaler, N., Brischke, C. et. al. (2016). Performance of Windows and façade elements made of thermally modified Norway spruce (Picea abies). in different climatic conditions. Proceedings of the WCTE 2016-World Conference on Timber Engineering, Vienna, 9.

11. Ugovsek, A., Subic, B., Starman, J., Rep, G., Humar, M., Lesar, B. et. al. (2018). Short-term performance of wooden windows and facade elements made of thermally modified and non-modified Norway spruce in different natural environments. Wood Material Science & Engineering, 14 (1), 42-47. doi: http://doi.org/10.1080/17480272.2018.1494627

12. Bonifazi, G., Serranti, S., Capobianco, G., Agresti, G., Calienno, L., Picchio, R. et. al. (2016). Hyperspectral imaging as a technique for investigating the effect of consolidating materials on wood. Journal of Electronic Imaging, 26 (1), 011003. doi: http://doi.org/ 10.1117/1.jei.26.1.011003

13. Jones, D., Sandberg, D., Goli, G., Todaro, L. (2019). Wood Modification in Europe: a state-of-the-art about processes, products and applications. Firenze University Press, 123. doi: http://doi.org/10.36253/978-88-6453-970-6

14. Janna, W. S. (2010). Engineering Heat Transfer. Boca Raton: CRC Press, 692.

15. Potter, M. C. (2018). Engineering analysis. New York: Springer, 444.

16. Temme, N. M. (1996). Special Functions: An Introduction to the Classical Functions of Mathematical Physics. Mathematics & Statistics. Applied Mathematics, 392. doi: http://doi.org/10.1002/9781118032572