Abstract
The single exponential kinetics (SEK) and parallel exponential kinetics (PEK) models were fitted to kinetic sorption data of welded and unwelded Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.). Furthermore, diffusion coefficients of water vapour in wood were determined using two different Fickian diffusion solutions. The objective was to identify how well these models could represent the moisture contents of the specimens and to characterize differences between the sorption behaviour of welded and unwelded wood. This knowledge can be used to enhance the moisture resistance of welded wood, develop drying schedules, and improve the quality of timbers. The PEK and SEK models provided the most precise and the second most precise fits to the sorption kinetic data, respectively. The two Fickian models are equivalent when both the infinite series are truncated at
1 Introduction
1.1 Wood welding and water sorption
The studies related to structural applications of welded wood are few in part due to the vulnerability of the welded joint to damage from moisture (Vaziri 2011). Enhancing the water resistance of the welded joint requires a full understanding of the water sorption mechanism, and the present study contributes to increasing this knowledge.
Since the 1980s, the dynamic vapour sorption (DVS) technique has been used to characterize the sorption isotherms where the time-dependent changes in the mass for different steps of relative humidity (RH) are recorded. The sorption data are usually used for fitting against advanced computer-based kinetic models. Recently, considerable progress has been made regarding the understanding and modelling of moisture transport processes in wood (Glass et al. 2018; Thybring et al. 2019a; Zelinka et al. 2018). These models are based on analytical expressions and a few parameters that need to be estimated by inverse analysis.
The sorption kinetics data of wood are often interpreted by a Fickian diffusion law (Jalaludin et al. 2010; Papadopoulos and Hill 2003; Pfriem et al. 2010; Salin 2010). Fick (1855) was the first who made this kind of modelling possible by placing diffusion on a quantitative basis (Avramidis 2007). Thybring et al. (2019b) reviewed a handful of popular and recent theoretical models.
The parallel exponential kinetic (PEK) model has been comparatively recently introduced to the wood science literature. The PEK model, first proposed for cellulosic materials by Kohler et al. (2003, 2006, is an extension of a single kinetic exponential model (SEK) (Glass et al. 2017). It has been found that the PEK model provides precise fits to the sorption kinetic data of natural fibres (Belbekhouche et al. 2011; Hill et al. 2010a; Kohler et al. 2003; Xie et al. 2011), regenerated cellulose (Okubayashi et al. 2004; 2005a, b), and wood (Hill et al. 2010b, 2010c; Jalaludin et al. 2010; Popescu and Hill 2013; Popescu et al. 2014; Sharratt et al. 2010).
Abundant speculations have been derived from PEK model parameters that mainly are based upon the assumption that the model consists of two processes: (1) a fast process associated with the monolayer formation of strongly bound water to the polar groups within the cell wall, and (2) a slow process associated with the water sorption to the multilayer of weakly bound water in the larger water clusters (Kohler et al. 2003, 2006). However, none of the claims is experimentally proven, and it is not yet clear what these two processes represent (Hill et al. 2010d; Xie et al. 2011).
Thybring et al. (2019a) showed by multi-exponential decay analysis (MEDEA) that the sorption kinetic data to equilibrium frequently contains more than two characteristic exponential components. Glass and co-workers (Glass et al. 2017, 2018; Thybring et al. 2019a) claimed that fitting PEK to the sorption kinetic data severely mischaracterizes both the equilibrium moisture content (EMC) and the full kinetic behaviour. Therefore, the PEK model cannot be used to derive physically meaningful properties from water vapour sorption measurements.
On the other hand, the degree of accuracy that Glass and co-workers expect from DVS experiments seems idealistic, especially for a complex hygroscopic polymer such as wood. Furthermore, considering the above-mentioned material properties making such a conclusion must be supported by a much larger data set than only one sample.
To the best of our knowledge, this is the first study in which the time-dependent transition of MC after a stepwise change in RH from one level to another (i.e., sorption kinetics) of welded wood has been studied. The purpose of this study was to compare the sorption behaviours of welded and unwelded wood at the same experimental conditions rather than to characterize moisture content (MC) accurately.
2 Materials and methods
2.1 The experiment
In this work, the welded and unwelded Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) were studied, and the experimental procedure is described in detail in Vaziri et al. (2023).
Considering Glass and co-workers’ recommendation (Glass et al. 2017), the stop criterion (
The specimens were cut with identical dimensions, but due to the rigidity of the welded wood, some specimens got irregular edges, which could be one of the error sources in the DVS experiment.
It should be considered that errors in EMC measurement with DVS are associated with error sources such as stop criterion, weighing of samples, dimensional measurement of samples, temperature, RH control, etc. Some of these uncertainties in EMC measurements are discussed by Berger et al. (2020), Glass et al. (2018), Willems (2022), and Zelinka et al. (2018).
2.2 Diffusion problems – or Fick’s law
The theoretical background to Fickian diffusion of this paper is mainly relying on Crank (1975). Many recent publications inaccurately cite the works of Crank. Therefore, going back to the original source is preferable.
2.2.1 Some theoretical background
Diffusion of water in wood refers to the molecular movement of bound water through the gross wood system by transferring through the cell-wall substance and lumen void volume from a higher water concentration area to a lower one (Avramidis 2007).
The diffusion coefficient can be either constant or dependent upon the concentration, pressure, temperature, and composition of the medium (Crank 1975). In general, the diffusion coefficients of wood are a function of the above-mentioned factors and the method of their determination (Keey et al. 1999; Stamm 1964).
If no dimensional changes occur in wood, the diffusion coefficient
The diffusion problem is defined as
where,
Concentration is hard to measure, but the mass of the specimen can be measured in the DVS. The relative mass change of the diffusing water and the relative concentration change of it is however equivalent. Integrating the space-dependent
over specimen thickness
where
2.2.2 The two Fickian-based models of fractional change in moisture content,
E
(
t
)
In Crank (1975) two alternative solution forms of Equation (3) are presented. Both solutions to the partial differential equation are in the forms of infinite series. One solution is more accurate for short times, whereas the other one is more accurate for long times for finite numbers of series components. For
The first of the two solutions which is determined by the method of separation of variables (space and time) yields
where,
The second solution which can be achieved using the Laplace transform yields
is more accurate for short times (Equation 4.18, Crank 1975), i.e., right after the diffusion process has started from a steady-state condition after a stepwise change in RH.
It is not precisely stated in Crank (1975) for which time intervals the two different solutions are valid for a given
2.2.3 The simplified versions of
M
C
(
t
)
used for plotting and root-mean-square error (RMSE) estimations
For expressing
can be used.
For full accuracy, the entire series of Equations (4) and (5) should be used. It is common in the literature to cut the series already at
2.2.4 Determination of the diffusion coefficients
In this paper, the diffusion coefficients were assumed to be constant for each RH step and therefore were determined for each of the RH steps.
For the short times model – the initial rates of sorption and desorption (ISORAD) method determining
The method for determining
which is Equation (4) truncated at
knowing
For numerical comparison, the ISORAD method results in
and the sorption method which was proposed in Section 10.6.5, Equation (10.158), Crank (1975) results in
2.3 The exponential models
2.3.1 The parallel exponential kinetics (PEK) model
The PEK model is expressed as
where
Equilibrium moisture content (EMC) is theoretically defined as
but the exact values of
2.3.2 The single exponential kinetics (SEK) model
SEK is a simplified version of the PEK model with only one exponential term:
where
The idea of testing this model was to see if there is an artificial compromise exponential process that can mimic the two physical exponential processes behind PEK in a “good-enough” manner.
2.4 The fitting procedure
Curve fitting was performed in MATLAB® R2020a (Mathworks, Natick, Massachusetts, USA) using the in-built fit function for non-weighted and non-linear least-squares fitting.
The input data for fitting the PEK and SEK models consisted of experimental data of time
Due to the limitations of the curve-fitting function in Matlab, one cannot easily create a model that automatically makes sure that:
and, in turn, automatically determines a breaking point between the two parameters. This had to be implemented manually by introducing the variable limits
to the algorithm. The numbers 16 and 260 were empirically chosen by experimenting combined with some engineering knowledge of the processes behind the equation. In alignment with Equations (15) and (16), the boundary condition
was produced for the SEK fitting model.
For increasing RH values (adsorption), the exponential coefficients
Figure 1 shows the data from the first 3.0 min of the experiment for 10 % RH. Due to the curvature of the fitting curves at the first minutes of the sorption, some approaches have removed these data points and shifted the time axis for better curve fitting (Glass et al. 2018; Hill et al. 2010c). The approach used here left these data and the time axis unaltered. Unlike Thybring et al. (2019b), the fitting in this study was done for all parameters to minimize the total residual of the curve against measurement data. In Thybring et al. (2019b, p. 727), on the other hand,
The starting point of the fitting algorithm was set to 1 for the PEK and SEK time constants
The optional parameter “DiffMinChange” was set to
3 Results and discussion
In all results, the behaviours of the short-time and long-time Fickian models are almost identical all the time and for all specimens. The two Fickian models are equivalent when both the infinite series are truncated at n = 10. Some interesting examples are shown here for different RH steps, adsorbing and desorbing, comparing unwelded and welded wood. The examples are to some extent taken arbitrarily, since an all-embracing bouquet of representative figures would be unrealistic to produce for the paper.
The differences between pine and beech were so small that it could not be motivated showing similar examples for pine and welded pine as are already shown in Figures 2–5 for beech and welded beech.
3.1 Adsorption
3.1.1 RH 0 %–5 %
Figure 2 presents a typical example of the model fits to the kinetic sorption data in adsorption in an RH step of 0 %–5 % comparing beech and welded beech. In Figures 2A and B, the MC of beech and welded beech for RH 0 %–5 % are compared. The MC change is higher for beech than welded beech. The deviation in MC between the replicates is higher for welded beech than for beech. Such spread between the replicates can be seen also for welded pine but the spread is much smaller than for welded beech. It is not clear whether it is due to degradation of wood by welding or just that the variation in specimen’s geometry is higher for welded wood.
Figures 2C and D show the prediction errors (residuals) for the first 50 min of the sorption process. The residuals are lower in absolute value for welded beech. Their characteristics also vary somewhat between welded beech and beech. From all these figures (especially Figure 2C and D) it can be seen that the SEK and PEK models provide better and different fits to the sorption data than the two Fickian models. For the unwelded wood case, SEK and PEK are identical, which also is verified by Appendix Figure A1.
Figures 3A and B focus on the later parts of the above given sorption process. For beech, the Fickian models fit better than the SEK and PEK to the sorption data at the end of this sorption step. For welded beech, PEK and the Fickian models equally fit better than SEK to the sorption data in the later parts of the sorption step.
Figures 3C and D illustrate the first 25 min of the sorption process. The SEK and PEK fit better than the Fickian models to the sorption data for both welded and unwelded beech.
3.1.2 RH 90 %–95 %
In Figure 2, all replicates start with the same MC level but later at 90 %–95 % RH in Figure 4, the spread between the replicates has increased throughout the experiment. For each stepwise change in RH, the deviation between the initial MC values of the three replicates increase further.
The MC is generally lower, and the MC increase is also smaller for welded beech in Figures 4A and B. The residuals of the model fits in Figures 4C and D are also lower for welded beech. Moreover, the residuals have completely different characteristics for welded wood, implying that other phenomena might lie behind the dynamic vapour sorption process of the welded joints than of untreated wood.
3.2 Desorption – RH 80 %–75 %
Figures 5A and B show desorption curves of welded and unwelded beech replicates for 80%–75 % RH. The Fickian models overlap each other and deviate considerably from SEK and PEK which much better resemble the behaviour of the MC curves. SEK has however large residuals in the beginning. The Fickian and SEK challenges are common for desorption.
Figures 5C and D show the residuals of the first few minutes of the desorption process. The Fickian and PEK models reveal similar fits, but the SEK model does not fit precisely to the sorption curves, which is comparable to the Fickian long-term model truncated too early (
Figures 5E and F show the residuals of the model fits for the entire desorption process. Focus is on the non-Fickian curves, disregarding the early large residuals of SEK. The residuals are similar for welded and unwelded beech at the beginning and end of the fits but differ in the middle in pattern. The SEK (ignoring the early stage) and PEK models provide precise fits to the welded and unwelded sorption kinetic data. The residuals of PEK do not exceed ±0.10 percentage points of MC (in the beginning, c.f. Figures 5E and F) and most of the time not more than 0.02 % percentage points. The Fickian models generally had a better fit to the sorption data in adsorption than desorption (Figures 2, 4, and 5 – and a summary in Figure 7).
3.3 Aggregating the results
Appendix Figure A1 shows that the SEK and PEK outperform the Fickian models. For most of the cases, the two Fickian models perform equivalently, whereas it is only in the 0 %–5 % RH step that SEK and PEK perform equivalently. PEK with more fitting parameters should in the general case provide a better fit than SEK.
All models have more challenges in higher RH and desorption, and combining high RH and desorption makes that even harder. This can be seen in Appendix Figures A1–A4 for the four different wood types, and in Figure 7 where an average is taken over the wood types.
Moreover, in Figure 6, one can see that the RMSEs of welded wood are in general smaller than those of unwelded wood. In addition to that, the RMSEs of pine and welded pine are slightly higher than those of beech and welded beech. Welded wood seems to be more Fickian in its properties. One possible explanation is that the rigid structure of the welded bond line restricts its swelling. The diffusion behaviour of many swelling polymers such as wood cannot be explained by Fick’s law with constant boundary conditions (Section 11.1, Crank 1975; Van der Wel and Adan 1999).
In Figure 7 the average RMSEs of all specimen types are shown for different stages of RH-targets. In that figure (and in Figure 6), the vertical axis is linear and not logarithmic as in Appendix Figures A1–A4. Having linear vertical axes for the RMSE representations emphasizes even more where the errors are large or small.
The precision of the PEK model fits decreases at higher RH, especially in adsorption but it is also notable for desorption. This phenomenon can for example be seen in Figure 7. Crank (1975) and Willems (2022) might have an explanation for this phenomenon when the involved heat exchange of MC-changes limits the sorption rates. That is, however, out of the scope of this paper.
All models showed smaller RMSE in adsorption than in desorption even if PEK has the worst performance for high-RH-adsorption (Figure 7).
4 Conclusions
The PEK model fits more precisely than the other studied models to the MC curves, particularly in adsorption. Nevertheless, the slightly higher RMSEs for the PEK at higher RH reveal that probably other phenomena simultaneously happen with water adsorption that this model cannot fully cover.
Generally, it is hard for the models to represent the MC curves at a high RH level for desorption. The Fickian-based models perform significantly poorly at high RH levels for desorption.
Welded wood seems to be more Fickian in its properties. One possible explanation is that due to the rigidity of the welded bond line, the welded joints swell less. In general, welded wood seems more suitable for the models used. The residuals are smaller than for unwelded wood.
Moreover, the RMSEs for pine and welded pine are slightly larger than those for beech and welded beech.
It is recommended for future work to investigate how welding parameters can affect the vapor sorption properties of the welded wood.
-
Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
-
Research funding: None declared.
-
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
References
Avramidis, S. (2007). Bound water migration in wood. In: Perré, P. (Ed.), Fundamentals of wood drying. A.R.BO.LOR, Nancy, France, pp. 105–124.Search in Google Scholar
Belbekhouche, S., Bras, J., Siqueira, G., Chappey, C., Lebrun, L., Khelifi, B., Marais, S., and Dufresne, A. (2011). Water sorption behaviour and gas barrier properties of cellulose whiskers and microfibrils films. Carbohydr. Polym. 83: 1740–1748, https://doi.org/10.1016/j.carbpol.2010.10.036.Search in Google Scholar
Berger, J., Colinart, T., Loiola, B.R., and Orlande, H.R. (2020). Parameter estimation and model selection for water sorption in a wood fibre material. Wood Sci. Technol. 54: 1423–1446, https://doi.org/10.1007/s00226-020-01206-0.Search in Google Scholar
Crank, J. (1975). The mathematics of diffusion, 2nd ed. Oxford: Clarendon.Search in Google Scholar
Fick, A. (1855). On liquid diffusion. Phil. Mag. J. Sci. 10: 30–39, https://doi.org/10.1080/14786445508641925.Search in Google Scholar
Glass, S.V., Boardman, C.R., and Zelinka, S.L. (2017). Short hold times in dynamic vapour sorption measurements mischaracterize the equilibrium moisture content of wood. Wood Sci. Technol. 51: 243–260, https://doi.org/10.1007/s00226-016-0883-4.Search in Google Scholar
Glass, S.V., Boardman, C.R., Thybring, E.E., and Zelinka, S.L. (2018). Quantifying and reducing errors in equilibrium moisture content measurements with dynamic vapour sorption (DVS) experiments. Wood Sci. Technol. 52: 909–927, https://doi.org/10.1007/s00226-018-1007-0.Search in Google Scholar
Grönquist, P., Frey, M., Keplinger, T., and Burgert, I. (2019). Mesoporosity of delignified wood investigated by water vapour sorption. ACS Omega 4: 12425–12431, https://doi.org/10.1021/acsomega.9b00862.Search in Google Scholar PubMed PubMed Central
Hill, C.A.S., Norton, A., and Newman, G. (2010a). The water vapour sorption behaviour of flax fibers–analysis using the parallel exponential kinetics model and determination of the activation energies of sorption. J. Appl. Polym. Sci. 116: 2166–2173, https://doi.org/10.1002/app.31819.Search in Google Scholar
Hill, C.A.S., Norton, A., and Newman, G. (2010b). Analysis of the water vapour sorption behaviour of Sitka spruce [Picea sitchensis (Bongard) Carr.] based on the parallel exponential kinetics model. Holzforschung 64: 469–473, https://doi.org/10.1515/hf.2010.059.Search in Google Scholar
Hill, C.A.S., Norton, A., and Newman, G. (2010c). The water vapour sorption properties of Sitka spruce determined using a dynamic vapour sorption apparatus. Wood Sci. Technol. 44: 497–514, https://doi.org/10.1007/s00226-010-0305-y.Search in Google Scholar
Hill, C.A.S., Moore, J., Jallaludin, Z., Leveneu, M., and Mahrdt, E. (2010d). The influence of earlywood/latewood and ring position upon the water vapour sorption properties of Sitka spruce. Int. Wood Prod. J. 2: 12–19.10.1179/2042645311Y.0000000001Search in Google Scholar
Jalaludin, Z., Hill, C.A.S., Samsi, H.W., Husain, H., and Xie, Y. (2010). Analysis of water vapour sorption of oleo-thermal modified wood of Acacia mangium and Endospermum malaccense by a parallel exponential kinetics model and according to the Hailwood–Horrobin model. Holzforschung 64: 763–770, https://doi.org/10.1515/hf.2010.100.Search in Google Scholar
Keey, R.B., Langrish, T.A.G., and Walker, J.C.F. (1999). Kiln-drying of lumber. Springer-Verlag, Berlin, Heidelberg, New York, p. 326.10.1007/978-3-642-59653-7Search in Google Scholar
Kohler, R., Renate, D., Bernhard, A., and Rainer, A. (2003). A numeric model for the kinetics of water vapour sorption on cellulosic reinforcement fibers. Compos. Interfaces 10: 255–257, https://doi.org/10.1163/156855403765826900.Search in Google Scholar
Kohler, R., Alex, R., Brielmann, R., and Ausperger, B. (2006). A new kinetic model for water sorption isotherms of cellulosic materials. Macromol. Symp. 244: 89–96, https://doi.org/10.1002/masy.200651208.Search in Google Scholar
Okubayashi, S., Griesse, U.J., and Bechtold, T. (2004). A kinetic study of moisture sorption and desorption of lyocellfibers. Carbohydr. Polym. 58: 293–299, https://doi.org/10.1016/j.carbpol.2004.07.004.Search in Google Scholar
Okubayashi, S., Griesser, U.J., and Bechtold, T. (2005a). Water accessibilities of man-made cellulosic fibers effects of fiber characteristics. Cellulose 12: 403–410, https://doi.org/10.1007/s10570-005-2179-y.Search in Google Scholar
Okubayashi, S., Griesser, U.J., and Bechtold, T. (2005b). Moisture sorption/desorption behaviour of various manmade cellulosic fibers. J. Appl. Polym. Sci. 97: 1621–1625, https://doi.org/10.1002/app.21871.Search in Google Scholar
Papadopoulos, A.N. and Hill, C.A.S. (2003). The sorption of water vapour by anhydride modified softwood. Wood Sci. Technol. 37: 221–231, https://doi.org/10.1007/s00226-003-0192-6.Search in Google Scholar
Pfriem, A., Zauer, M., and Wagenfuhr, A. (2010). Alteration of the unsteady sorption behaviour of maple (Acer pseudoplatanus L.) and spruce (Picea abies (L.) Karst.) due to thermal modification. Holzforschung 64: 235–241, https://doi.org/10.1515/hf.2010.029.Search in Google Scholar
Popescu, C.M. and Hill, C.A.S. (2013). The water vapour adsorption– desorption behaviour of naturally aged Tilia cordata Mill. wood. Polym. Degrad. Stabil. 98: 1804–1813, https://doi.org/10.1016/j.polymdegradstab.2013.05.021.Search in Google Scholar
Popescu, C.M., Hill, C.A.S., Curling, S., Ormondroyd, G., and Xie, Y. (2014). The water vapour sorption behaviour of acetylated birch wood: how acetylation affects the sorption isotherm and accessible hydroxyl content. J. Mater. Sci. 49: 2362–2371, https://doi.org/10.1007/s10853-013-7937-x.Search in Google Scholar
Salin, J.G. (2010). Problems and solutions in wood drying modeling: history and future. Wood Mater. Sci. Eng. 5: 123–134, https://doi.org/10.1080/17480272.2010.498056.Search in Google Scholar
Sharratt, V., Hill, C.A.S., Zaihan, J., and Kint, D.P.R. (2010). Photodegradation and weathering effects on timber surface moisture profiles as studied using dynamic vapour sorption. Polym. Degrad. Stab. 95: 2659–2662, https://doi.org/10.1016/j.polymdegradstab.2010.07.011.Search in Google Scholar
Sparr, G. and Sparr, A. (2000). Kontinuerliga system. Studentlitteratur, Sweden.Search in Google Scholar
Stamm, A.J. (1964). Wood and cellulose science. The Ronald Press Company, New York, p. 549.Search in Google Scholar
Thybring, E.E., Boardman, C.R., Glass, S.V., and Zelinka, S.L. (2019a). The parallel exponential kinetics model is unfit to characterize moisture sorption kinetics in cellulosic materials. Cellulose 26: 723–735, https://doi.org/10.1007/s10570-018-2134-3.Search in Google Scholar
Thybring, E.E., Glass, S.V., and Zelinka, S.L. (2019b). Kinetics of water vapour sorption in wood cell walls: state of the art and research needs. Forests 10: 704, https://doi.org/10.3390/f10080704.Search in Google Scholar
Van der Wel, G.K. and Adan, O.C.G. (1999). Moisture in organic coatings: a review. Prog. Org. Coat. 37: 1–14, https://doi.org/10.1016/s0300-9440(99)00058-2.Search in Google Scholar
Vaziri, M. (2011). Water resistance of Scots pine joints produced by linear friction welding. Ph.D. thesis. Skellefteå, Luleå University of Technology.Search in Google Scholar
Vaziri, M., Dreimol, C., Abrahamsson, L., Niemz, P., and Sandberg, D. (2023). Water-vapour sorption of welded bond-line of European beech and Scots pine. Holzforschung 77, It is published now.10.1515/hf-2022-0012Search in Google Scholar
Willems, W. (2022). Improved boundary conditions for models of coupled heat and water vapour transfer between an adsorbent and the ambient. Int. J. Therm. Sci. 172: 107296, https://doi.org/10.1016/j.ijthermalsci.2021.107296.Search in Google Scholar
Xie, Y., Hill, C.A.S., Jalaludin, Z., and Sun, D. (2011). The water vapour sorption behaviour of three celluloses: analysis using parallel exponential kinetics and interpretation using the Kelvin–Voigt viscoelastic model. Cellulose 18: 517–530, https://doi.org/10.1007/s10570-011-9512-4.Search in Google Scholar
Zelinka, S.L., Glass, S.V., and Thybring, E.E. (2018). Myth versus reality: do parabolic sorption isotherm models reflect actual wood-water thermodynamics? Wood Sci. Technol. 52: 1701–1706, https://doi.org/10.1007/s00226-018-1035-9.Search in Google Scholar
© 2023 the author(s), published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.