Original Article

Water Sorption/Desorption Kinetics and Convective Drying of Eucalyptus globulus Wood

Mahyoub AMER2,, Bousselham KABOUCHI2, Salah El ALAMI2, Brahim AZIZE2, Mohamed RAHOUTI3, Abderrahim FAMIRI4, Abdelwahed FIDAH4
Author Information & Copyright
2Laboratory of Condensed Matter and Interdisciplinary Sciences, Faculty of Sciences, Mohammed V University in Rabat, 4 Avenue Ibn Battouta B.P. 1014 RP, Morocco
3Center of Plant and Microbial Biotechnologies, Biodiversity and Environment, Faculty of Sciences, Mohammed V University in Rabat, 4 Avenue Ibn Battouta B.P. 1014 RP, Morocco
4Physics Mechanics of Wood Laboratory, Research Centre of Forestry in Rabat, Charia Omar Ibn Al Khattab. B.P, 763 Agdal-Rabat 10050 Morocco
Corresponding author: Mahyoub AMER (e-mail: mahyoub51973@gmail.com, ORCID: 0000-0003-0547-8970)

© The Korean Society of Wood Science & Technology.

Received: Apr 05, 2019; Accepted: Aug 10, 2019

Published Online: Sep 25, 2019

Abstract

Radial and tangential water diffusion in Eucalyptus globulus wood was investigated using three mature trees from a forest in Khemis Sahel (North Morocco). Absorption and desorption kinetics experiments were conducted at ambient temperature (25 °C) and 30 °C, respectively, and a relative humidity of 60%. The diffusion coefficients in the two directions were determined under imposed hygrothermal conditions; they were greater in the radial direction for the absorption as well as desorption processes. Convective drying under load, preceded by reconditioning and followed up by balancing, revealed the drying conditions that corresponded to the appropriate drying schedules for E. globulus wood. This was verified by measuring the cracks and bowsbefore and after drying of boards.

Keywords: Eucalyptus globulus; absorption kinetic; desorption kinetic; diffusion coefficient; convective drying

1. INTRODUCTION

Eucalyptus trees have a major socio-economic role in Morocco, due to their plasticity and rapid growth capacity (Famiri et al., 2001). Indeed, Eucalyptus wood is used in the industry as pulp, logs, firewood and charcoal. However, the uses of eucalyptus as a timber remain limited because the appearance of cracks and/or splits due to the growth stresses releasing during cutting, sawing and drying. In addition, the eucalyptus exhibit high shrinkage and collapses (Maziri et al., 2010; Wessels et al., 2016; Amer et al., 2017).

Beneficial use of eucalyptus such Eucalyptus globulus wood needs searching methods to reduce these defects. This requires a study of water transfer in this kind of wood and control the drying rate. Diffusivity of water is an important phenomenon in wood drying studies (Peralta and Bangi, 2003; Kouchade, 2004; Ghazil, 2010).

Wood is a hygroscopic material. In storage, it is exposed to both absorption and desorption processes of moisture from the surrounding air (Baronas et al., 2001; Okoh, 2014). Moisture content significantly influences wood physical and mechanical properties, which are of great importance for wood uses in different applications areas (Shi, 2007; Noorolahi et al., 2008; Sonderegger et al., 2011; Franke et al., 2016; Lee and Lee, 2018; Pang and Jeong, 2019).

Water transfer in wood is a very complex process. It diffuses in three different forms: free water in wood above fiber saturation point (FSP), bound water and water vapor in the wood below FSP point (Vasic et al., 2012; Engelund et al., 2013; Okoh, 2014; Bennani et al., 2017). Water diffusion in wood depends on numerous factors, such as wood moisture content, diffusion direction, time, temperature, relative humidity and type of wood (Jannot et al., 2006; Tamme et al., 2011; Sandoval-Torres et al., 2014; Kim et al., 2017).

Drying eucalyptus wood can be done naturally or industrially, but conventional kiln drying is the most common practice among the artificial drying methods used in the wood industry (Travan et al., 2010; Chang et al., 2017; Kang and Kang, 2018). The conventional kilns typically operate in temperatures ranging from 40 to 80 °C (Batista et al., 2013).

In this work, we studied the water diffusion in E. globulus wood in radial and tangential directions at a fixed temperature and relative humidity values by analyzing absorption and desorption kinetics of water and estimating the diffusion coefficients. We also investigated convective drying tests in order to reduce defects and get a good quality of drying for the considered wood.

2. MATERIALS and METHODS

2.1. Materials

The vegetal material used in this study consists of three trees of E. globulus Labill. from the Khmiss Sahel forest in Larache region (North of Morocco). The selected trees are more than 30-years-old. The corresponding circumferences are ranged between 100 and 118 cm. These trees have a good straightness and are free of defects. In Table 1, we reported dendrometric characteristics measurements of sampled trees.

Table 1. Dendrometric measurements of the studied trees
Tree ht (m) C1.3 (m) Cbase (m) MC (%)
1 15 1.18 1.90 91
2 14 1.03 1.82 89
3 13 1.00 1.80 89

where ht: total tree height, C1.3: circumference at 1.3 m, Cbase: circumference at the base, and MC: moisture content.

Download Excel Table
2.2. Preparation of test specimens

The wood samples used in this study were taken from the heartwood along the same longitudinal bars and the first log of each tree. The selected samples have a parallelepipedic form with dimensions (20 × 20 × 15) mm3 (Fig. 1). We have prepared 40 specimens for studying absorption kinetics in the radial direction (20 samples) and tangential direction (20 samples). The same number (40 specimens) is used in the desorption study. The samples were covered with aluminum tape on the four sides. The water diffusion in radial or tangential directions is studied along the two other open sides.

wood-47-5-557_F1
Fig. 1. Test specimen preparation for sorption and desorption kinetic studies.
Download Original Figure
2.3. Experimental

Water diffusion in E. globulus wood was carried out by using a water basin for absorption study and a climatically chamber for desorption measurement (Fig. 2). The climatic conditions (temperature T and relative humidity Hr) surrounding the samples were controlled during desorption experiments.

wood-47-5-557_F2
Fig. 2. Experimental material used in desorption experiments.
Download Original Figure
2.4. Measurement method

To follow water absorption kinetics in the wood, the samples which have initially moisture content about 10% were submerged in the water basin at ambient temperature (25°C). The changes in the samples weight are recorded at regular intervals time (each 2 min) until saturation in 24 hours.

In the study of desorption kinetic, the moisture content in the wood samples was almost 60% initially. The samples were placed in climatically chamber at temperature 30 °C and relative humidity 60%, their weights were measured at regular time intervals during 24 hours until equilibrium state.

The relative humidity and samples weights were measured respectively by a digital hygrometer of accuracy (0.01 %) and an electronic balance with precision (0.001 g) controlled with a computer.

The wood samples were placed in an oven at 103 °C during 48 hours to obtain the dry mass and to determine the moisture content.

2.5. Mathematical formulations

The wood moisture content denoted (MC) is equal to the quantity of water contained in the wood, expressed as a percentage of dry weight:

MC % = W h W o W o × 100
················ (1)

where Wh: wet wood mass (g), and Wo: dry wood mass (g).

Driving force of transfer is moisture gradient. The following approach treats the internal transfer mechanism of moisture diffusion in wood. According to Fick's second law, the diffusion equation is given by (Houngan et al., 2015):

C t = d dx D dC dx
················ (2)

where C = C(x, t): wood moisture concentration (g cm-3) at the x thickness (cm) in time t (s), and D: the diffusion coefficient (m2s-1). By defining a dimensionless water content C+:

C + = C C eq C in C eq MC MC eq MC in MC eq
················ (3)

where Ceq and MCeq: average samples moisture concentration (g cm-3) and moisture content (%) at equilibrium according to the hygro-thermal imposed conditions, Cin and MCin: average initial moisture concentration (g cm-3) and moisture content (%) in the samples, C and MC: average moisture concentration (g cm-3) and moisture content (%) in wood samples at the time t (s).

Assuming that D is constant, equation (2) becomes:

C + t = D d 2 C + dx 2
················ (4)

Analytic solution for equation 4 is given by the following relation (Kouchade, 2004):

C + t = 4 e D π t
················ (5)

Referred to gain or loss of the dimensionless mass by:

m t = Δ m t Δ m max = W t W in W eq W in
················ (6)

Where Δm and Δmmax represent the variation of the sample mass (g) at time t (s) and at equilibrium; W(t): weight at any time (g); Win: initial weight (g), and Weq: weight at equilibrium (g).

The experimentally measured magnitude m*(t) and the calculated analytically C+(t) being dimensionless. Thus, C+(t) can be written C+(t) ≡ m*(t). The slope linear part in curve (Fig. 3) represent the variation of m*(t) depending √t . We can determine D by using the expression (5).

wood-47-5-557_F3
Fig. 3. Variation of m* (t) depending √t for a test sample
Download Original Figure

Value of diffusion coefficient D is calculated by the direct method of half-rise time t0,5 which is corresponding to m*(t) = 0.5. In this time, the diffusion coefficient D is expressed by the following equation (Kouchade, 2004):

D = 0.049 e 2 t 0.5
················ (7)

where e: sample thickness (cm).

End splits of wood boards are determined by a parameter called splitting indicator Ie, which is given by the following relation (Hakam et al., 2005):

Ie max = f max L
················ (8)
Ie aver = 1 L 1 n f n n
················ (9)
Δ Ie k = Ie after Ie before
················ (10)

where L: boards length (cm); n: split number; fn: depth of nth split (cm); fmax: depth of maximal split (cm); Iebefore: splitting indicator before drying; Ieafter: splitting indicator after drying, and k: average or maximum.

2.6. Drying tests

Drying tests are realized in a traditional wood dryer type Cathild installed at the Forest Research Center (FRC) in Rabat. The heating in the dryer was performed by electrical resistances. The drying conditions were controlled automatically by connection to a computer with specific software, which permits visualization and control in real-time drying parameters evolution, such as relative humidity, temperature, air velocity and equilibrium moisture. The full automatic control system is used for the experiments. The dryer is filled with boards stacked as horizontal rows and separated by stickers with a section of 2 × 2 cm2, six stickers per row (Hakam et al., 2005). Drying tests were carried out on boards with length 2.50 m, width from 25 to 30 cm (according to the diameter of the log) and thickness of 3 cm.

Drying conditions are based on drying techniques of hardwood and softwood species which were subject to several studies conducted by TCWFI (Technical Center for Wood and Furnishing Industries) and other authors for the eucalyptus wood species (Kantay et al., 2002, Hakam et al., 2003), which have obtained wood products with a good quality and lower deformation.

Many experiments were done on drying E. globulus wood to achieve optimum drying conditions, minimize damage and distortions (collapse, splits, cracks, etc.) and improve the wood quality. The best results of these experiments are shown in Tables 2 and 3. In the first test of drying conditions, the temperature varied from 46 to 60 °C, while in the second test temperature varied between 40 and 70 °C.

Table 2. Drying schedule test 1
Sequence V (%) Vt 50% R & 50% L T gradient (°C/h) t (h) MC (%) T (°C) RH (%) EMC (%)

S1 Warming 1 99 2 6 46 90 17.0
S2 Warming 2 99 4 6 48 92 17.5
99 > 50 48 87 16.5
S3 Drying above FSP 99 50 - 40 48 83 15.5
90 40 - 35 48 79 14.2
85 35 - 30 48 77 13.3
75 30 - 27 48 69 11.5
70 27 - 24 50 63 10.2
65 24 - 21 52 60 9.5
S4 Drying below FSP 60 9 21 - 18 54 54 8.1
55 18 - 15 56 48 7.5
50 15 - 12 58 45 6.8
40 12 - 9 60 41 6.2
S5 Equalizing 40 4 12 12 50 36 6.0
S6 Conditioning 40 4 40
Download Excel Table
Table 3. Drying schedule test 2
Sequence V (%) Vt 50% R & 50% L T gradient (°C/h) t (h) MC (%) T (°C) RH (%) EMC (%)

S1 Warming 1 49 2 8 40 88 17.8
S2 Warming 2 70 4 4 44 84 16.4
90 > 50 48 82 15.2
S3 Drying above FSP 90 9 50 - 40 48 79 14.1
90 40 - 35 50 77 13.2
85 35 - 30 50 73 12.4
80 30 - 27 55 70 11.0
70 27 - 24 58 63 9.6
60 24 - 21 60 55 8.1
S4 Drying below FSP 55 9 21 - 18 62 52 7.4
50 18 - 15 66 46 6.4
49 15 - 12 68 39 5.5
49 12 - 9 70 35 4.8
S5 Equalizing 49 4 8 12 55 40 6.0
S6 Conditioning 49 5 40

where V: ventilation velocity (%) of the maximum value, Vt : ventilation time (hour); 50% right & 50% left, T: temperature (°C), t: time of sequence (hour), RH: relative humidity (%), EMC: equilibrium moisture content (%), FSP: fiber saturation point, MC: wood moisture content (%). T gradient: upwards or downwards slope of the temperature (°C/hour).

Download Excel Table

3. RESULTS and DISCUSSION

3.1. Absorption and desorption phenomena
3.1.1. Absorption kinetics in the radial and tangential directions

General profiles of water sorption kinetics in radial and tangential directions for wood samples are almost similar (Fig. 4). The obtained results of absorption kinetics in the radial and tangential directions show that the moisture of wood samples in the radial direction increased from 10 to 32.5 % during 24 hours. While in the tangential direction it increased from 10 to 30.8 % in the same period. By comparing the curves of absorption kinetics in the two directions, we can deduce that the E. globulus wood is slightly more diffusive in the radial direction. In this direction, the vascular rays are the main source of the water diffusion; especially their structure, arrangement and positioning, which promote the phenomenon of diffusivity (Kouchade, 2004). In the tangential direction, the early-wood and late-wood layers are arranged in parallel, which promotes more the diffusion in the radial direction.

wood-47-5-557_F4
Fig. 4. Moisture absorption as a function in time for the radial and tangential directions.
Download Original Figure
3.1.2. Desorption kinetics in the radial and tangential directions

Fig. 5 shows that the moisture content in the two directions decreases continuously with time. In the curves of desorption, it can be seen that the moisture content decreases rapidly in the initial time. Within the first six hours, in the radial direction, the moisture decreased from 60 to 22%, whereas it diminished from 60 to 27% in the tangential direction. At the end of the desorption process (24 hours), the moisture content reached almost 10% and 14% in the radial and tangential directions respectively.

wood-47-5-557_F5
Fig. 5. Moisture desorption variation with time in radial and tangential directions.
Download Original Figure

In the desorption curves, we can distinguish three phases. During the first phase, the samples moisture content loss is very fast; this can be explained by the elimination of the free water through the lumens and pits of the wood structure. The second phase is characterized by a reduction of loss velocity of the moisture content; the wood samples are in the hygroscopic region (below the FSP), where the bound water and water vapor are transferred. The third phase is characterized by a very slowdown of desorption which practically becomes zero. This is mainly due to the nearness of the moisture content of samples with that of the surrounding air.

3.1.3. Diffusion coefficient

According to the radial and tangential directions in absorption and desorption processes, the water diffusion coefficient of wood samples was calculated from equation (6). The mean values of the water diffusion coefficient in two directions in both processes and the average of their variation coefficients were reported in Table 4. The results show that the average values of the diffusion coefficient in desorption are higher than those of absorption in the two directions. We also noted that these coefficients are larger in the radial direction than in the tangential direction for both diffusion processes. These results are similar to those reported by other authors for other wood species (Siau, 1995; Khouya, 2008).

Table 4. Radial and tangential diffusion coefficient (D)
Radial direction
D × 10-9
(m2s-1)
CV
(%)
Tangential direction
D × 10-9
(m2s-1)
CV
(%)
Absorption 1.02 2.00 0.96 2.10
Desorption 2.30 2.20 1.19 1.70
Download Excel Table
3.2. Eucalyptus globulus wood drying

Fig. 6 shows the drying kinetics of E. globulus wood concerning the two drying schedules described above (tests 1 and 2). It presents the evolution of the average moisture of wood as a function of time according to the conditions imposed for each test. The results obtained show that the program of drying test 1 reduce the average moisture content in wood boards from 56 to 12 % during 222 h. However, the second drying test was more accelerated to reduce the wood moisture content from 56 to 12% in the drying cycle of 174 h. It should be noted that a good drying process must take into account the quality of the final product and the energy consumed. Therefore, drying E. globulus wood according to the conditions in test 1 provides the best results concerning dried wood quality.

wood-47-5-557_F6
Fig. 6. Drying kinetics for E. globulus wood according to drying tables
Download Original Figure

The measurement of the boards splits before and after drying showed that the drying tests 1 and 2 reduce the deformation and degree of cracks in E. globulus wood. The value of splitting indicator (ΔIe) in test 1 (0.25%) is less than that obtained in test 2 (10.43%).

The analysis of the values of different bows classes, before and after drying, is reported in Table 5. This table indicates that the drying tests affect bows values. The more important variations are observed for test 1.

Table 5. Bows distribution before and after lumber drying
Classes Bows (mm) Test 1 Test 2
Before (%) After (%) Difference (%) Before (%) After (%) Difference (%)
0-3 34.32 42.80 8.48 40.20 42.60 2.40
3-6 25.20 29.50 4.30 31.63 33.05 1.42
6-9 16.00 14.45 -2.00 12.50 13.25 0.70
9-12 9.18 4.25 -4.93 6.77 8.40 1.63
12-15 5.80 5.50 -0.30 2.50 2.35 -0.15
> 15 9.50 3.50 -6.00 4.00 5.14 1.14
Download Excel Table

4. CONCLUSION

We have investigated the absorption and desorption kinetics of water in the E. globulus wood in the radial and tangential directions. We also have determined the diffusion coefficients for this wood. The results show that E. globulus wood has a slightly higher diffusivity in the radial direction than in the tangential direction. The average values of the diffusion coefficient in the radial direction are larger than that of the tangential direction in both of absorption and desorption process. Drying tests done on E. globulus wood shows that the appropriate drying conditions for this wood correspond to the drying schedule 1 (Table 2). This later also reduce the deformations during the drying process.

References

1.

Amer M., Kabouchi B., Rahouti M., Famiri A., Fidah A. Determination of growth stresses indicator, moisture profiles and basic density of clonal eucalyptus wood. Journal of the Indian Academy Wood Science. 2017; 14(1):91-98.

2.

Baronas R., Ivanauskas F., Juodeikien L., Kajalavičius A. Modelling of moisture movement in wood during outdoor storage. Nonlinear Analysis: Modelling and Control. 2001; 6(2):3-14.

3.

Batista D.C., Klitzke R.J., Da Rocha M.P., Bolzon de Muñiz G.I., Batista T.R. Volume loss as a tool to assess kiln drying of Eucalyptus Wood. Floresta e Ambiente. 2013; 20(2):250-256.

4.

Bennani L., Elkouali M., Talbi M., Ainane T. Modelling the absorption process of water in wood in the transient regime. International Journal of Chemical Science. 2017; 15(2):137-147.

5.

Chang Y.S., Shin H.K., Kim S., Han Y., Kim M.J., Eom C.D., Lee Y.G., Shim K.B. Evaluation of drying properties and yields of domestic Quercus species for enhancing utilization. Journal of the Korean Wood Science and Technology. 2017; 45(5):622-628.

6.

Engelund E.T., Thygesen L.G., Svensson S., Hill C.A.S. A critical discussion of the physics of wood-water interactions. Wood Science and Technology. 2013; 47:141-161.

7.

Famiri A., Kabouchi B., Hakam A., Gril J. Sawing and growth stresses in green wood of Eucalyptus grandis & E. gomphocephala. Forest Science, Bulgaria. 2001; 1/2:45-50.

8.

Franke B., Franke S., Schiere M., Müller A. Moisture diffusion in wood-Experimental and numerical investigations. 2016; In: WCTE World Conference on Timber Engineering; 22-25 August (2016), Vienna-Austria p. 1-8.

9.

Ghazil A. Etude de la migration des fluides dans le bois. PhD. Thesis. 2010 France: Université Henri Poincaré, Nancy 1. .

10.

Hakam A., Dikrallah A., Kabouchi B., Famiri A., Walia-Allah M., El Abid A. Eucalyptus wood drying. Journal of Physics IV France. 2005; 123:327-330.

11.

Hakam A., Dikrallah A., Walia-Allah M., Famiri A., Kabouchi B., El Abid A., Fechtal M. Contrôle et quantification de la dégradation mécanique du bois pendant le séchage. Cas d’eucalyptus. Physical Chemistry News. 2003; 13:110-114.

12.

Houngan C.A., Awanto C., Houndedako S., Anjorin M., Vianou A. Mass diffusivity determination of Teak wood (Tectona grandis) used as building material. Int. Conference on Computational Heat and Mass Transfer, Procedia Engineering. 2015; 127:201-207.

13.

Jannot Y., Kanmogne A., Talla A., Monkam L. Experimental determination and modelling of water desorption isotherms of tropical woods: afzelia, ebony, iroko, moabi and obeche. Holz als Roh. und Werkstoff. 2006; 64:121-124.

14.

Kang C.W., Kang H.Y. Development of a Kiln dry schedule for Lindera erythrocarpa grown in Hongsung, Chungnam Province, Korea. Journal of the Korean Wood Science Technology. 2018; 46(1):10-16.

15.

Kantay R., Unsal O., Korkut S. Drying problems of fast growing tree species: Evaluation of Maritime pine (Pinus pinaster Ait) and Eucalyptus (Eucalyptus camaldulensis Dehn) wood. 2002; In: Conference, Meeting Management of Fast Growing Plantations (IUFRO 2002), Izmit-Turkey p. 208-212.

16.

Khouya A. Contribution aux études expérimentale et numérique d’un processus de séchage du bois. PhD. Thesis. 2008 Tétouan, Morocco: Université Abdelmalek Essaadi. .

17.

Kim H., Han Y., Park Y., Yang S.Y., Chung H., Eom C.D., Lee H.M., Yeo H. Finite difference evaluation of moisture profile in Boxedheart large-cross-section Square timber of Pinus densiflora during high temperature drying. Journal of the Korean Wood Science and Technology. 2017; 45(6):762-771.

18.

Kouchadé A.C. Détermination en routine de la diffusivité massique dans le bois par méthode inverse à partir de la mesure électrique en régime transitoire. PhD. Thesis. 2004 Nancy, France: Ecole Nationale du génie Rural des Eaux et des Forets. .

19.

Lee J.M., Lee W.H. Dimensional stabilization through heat treatment of thermally compressed wood of Korean Pine. Journal of the Korean Wood Science and Technology. 2018; 46(5):471-485.

20.

Maziri A., El Ghorba M., Chergui M., Famiri A., Ziani M., Kabouchi B. Etude des contraintes de croissance chez l’Eucalyptus camaldulensis et leur relation avec les fentes d’abattage. Physical Chemistry News. 2010; 53:15-21.

21.

Noorolahi S., Khazaei J., Jafari S. Modeling Cyclic water absorption and desorption characteristics of three varieties of wood. 2008; In: World conference on agricultural information and IT; 24-27 August (2008), Tokyo, Japan p. 13-22.

22.

Okoh E.T. Water absorption properties of some tropical timber species. Journal of Energy and Natural Resources. 2014; 3(2):20-24.

23.

Pang S.J., Jeong G.Y. Effects of density, temperature, size, grain angle of wood materials on nondestructive moisture meters. Journal of the Korean Wood Science and Technology. 2019; 47(1):40-50.

24.

Peralta P.N., Bangi A.P. A nonlinear regression technique for calculating the average diffusion coefficient of wood during drying. Wood Fiber Sci. 2003; 35(3):401-408.

25.

Sandoval-Torres S., Hernández-Bautista E., Rodríguez-Ramírez E., Carrillo Parrab E. Numerical simulation of warm-air drying of Mexican softwood (Pinus pseudostrobus): An empirical and mechanistic approach. Chemical and Biochemical Engineering Quartely. 2014; 28:125-133.

26.

Shi S.Q. Diffusion model based on Fick’s second law for the moisture absorption process in wood fiber-based composites: is it suitable or not. Wood Science and Technology. 2007; 41(8):645-658.

27.

Siau J. Wood influence of water on physical properties. PhD. Thesis. 1995Virginia Polytechnic Institute and State University. .

28.

Sonderegger W., Vecellio M., Zwicker P., Niemz P. Combined bound water and water vapour diffusion of Norway spruce and European beech in and between the principal anatomical directions. Holzforschung. 2011; 65:819-828.

29.

Tamme V., Muiste P., Mitt R., Tamme H. Determination of effective diffusion coefficient and mechanical stress of pine wood during convective drying. Baltic Forestry. 2011; 17:110-117.

30.

Travan L., Allegretti O., Negri M. Eucalyptus drying process: qualitative comparison of different clones cultivated in Italy. 2010; In: The Final Conference of COST Action E53; 4-7 May (2010), Edinburgh, UK p. 57-70.

31.

Vasic M., Radojevic Z., Grbavcic Ž. Calculation of the effective diffusion coefficient during the drying of clay samples. Journal of the Serbian Chemical Society. 2012; 77(4):523-533.

32.

Wessels C.B., Crafford P.L., Toit B. Du., Grahn T., Johansson M., Lundqvist S.O., Sall H., Seifert T. Variation in physical and mechanical properties from three drought tolerant Eucalyptus species grown on the dry West coast of Southern Africa. European Journal of Wood and Wood Products. 2016; 74:563-575.