Next Article in Journal
Variably Sized and Multi-Colored Silica-Nanoparticles Characterized by Fluorescence Correlation Methods for Cellular Dynamics
Previous Article in Journal
Skew Rolling of Bimetallic Rods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling and Experimental Study of Moisture Transport in Wood, Taking into Account Diffusion and the Accompanying Adsorption of Water Vapour by Cell Walls

Civil Engineering and Architecture Faculty, Opole University of Technology, 45-061 Opole, Poland
*
Author to whom correspondence should be addressed.
Materials 2021, 14(1), 17; https://doi.org/10.3390/ma14010017
Submission received: 17 November 2020 / Revised: 16 December 2020 / Accepted: 21 December 2020 / Published: 22 December 2020
(This article belongs to the Section Materials Physics)

Abstract

:
The paper presents a model of moisture transport in wood taking into account diffusion and the accompanying adsorption of water vapour through the skeleton. A two-parameter form of the source term was proposed, depending on the distance of the current moisture content (MC) from the equilibrium state. The tests on cubic samples with a side of 2 cm were carried out which allowed to determine the coefficients of the proposed model on the basis of the reverse method. The tests were performed for pine, larch, oak and ash in all directions of orthotropy. Tests on thin samples were also performed to verify the source term.

1. Introduction

Wood is widely used in construction as a structural and finishing material. The properties of wood such as strength, dimensional stability and durability are strongly dependent on its humidity. This means that the spatial distribution of moisture and its changes versus time play a key role in predicting its behaviour. The most important from the point of view of material engineering is the transport of moisture in the range of the moisture content of the wood corresponding to dry conditions and the moisture content occurring at full saturation of the fibres, i.e., in the moisture content from 0% to 30% [1,2]. In this range of moisture in wood there is water in the form of water vapour and water bound by surface forces, and changes in its content lead to shrinkage or swelling of the wood [3].
One of the first models for predicting moisture distributions in the hygroscopic range under isothermal conditions are based on the balance of total moisture diffusion. These models, although they do not fully explain the physical problem and the results obtained by means of the diffusion equation cannot always be matched to experimental data, which was analysed by Shi [4], are attractive due to the small number of parameters needed to describe the phenomenon. On the other hand, if the diffusion coefficient is assumed to depend on the moisture content, a very good fit of the model to experimental data is obtained [5,6].
Recently, more advanced models have also been developed, which analyse the combined flow of vapour and bound water [7,8,9,10,11,12]. These models assume that vapour and bound water are transported by diffusion. However, the form of the diffusion coefficient is not entirely clear. Siau [13] assumed that the resistance to vapour flow is influenced by both cell lumens and cell walls. Therefore, the coefficient can only be considered as an apparent diffusion coefficient. In fact, it is assumed that the vapour passes through the cell walls. On the one hand, it is adsorbed on the cell wall and then desorbed on the other. In this case, the diffusion coefficient should depend on the moisture content, e.g., in [14]. However, this is only one theory of the mechanism of vapour transmission. According to Dinwoodie [15], the main vapour transport route leads through pits in cell walls. A diffusion coefficient independent of moisture content is then to be expected. This case was assumed and confirmed by the following analyses.
The issue of the transport of bound water itself is also not fully explained. According to Kärger and Valiullin [16] diffusion in mesopores is the result of the interplay of several mechanisms. Thus, at partial pore-space saturation one may distinguish between surface diffusion and the process of mass transfer through the pore gas phase. On the other hand Hozjan and Svensson [11] claim that the diffusion of bound water in the porous structure is very slow compared to the diffusion of vapour. Therefore, for moisture transport in larger wood samples, the transport of bound water is minor importance and may be omitted [11,17]. However, it should take into account that this assumption is valid only for relative air humidity below 65% [18].
Another important issue is the possibility of surface resistance on boundary surfaces of wooden samples. This implies the adoption of more complex boundary conditions, e.g., boundary conditions of the third type instead of the first. This issue was considered experimentally by Rosen [19] analysing the effect of air velocity on water vapour adsorption for both longitudinal and transverse directions. He noted that for air velocities above 3 m/s, this effect is negligible and can be ignored.
The last issue, and probably the most important one, is the form of the source term. Its form has a very large impact on the results obtained from the model. A wrong one may give the apparent impression that, for example, the vapour diffusion coefficient is variable or that the transport of bound water is of great or minor importance.
The literature review shows that the topics related to moisture transport in the hygroscopic range are not fully explained and should be further analysed. This article proposes a model of moisture transport in wood that takes into account the diffusion and accompanying adsorption of water vapour through the skeleton. A new form of the source term has been proposed, which is an original contribution to modelling moisture flows in wood compared to those found in the literature. Its form was verified on thin pine samples for the relative humidity range from 30% to 50% and from 50% to 70%. The constancy of water vapour diffusion coefficients in the model for each of the orthotropic directions, boundary conditions of the first type for water vapour and boundary conditions with a source of moisture sorption for bound water were adopted. On the basis of measurements of changes in the mass of samples by minimising the target function, material coefficients were determined for the proposed model. The tests were performed on samples of the following dimensions 2 cm × 2 cm × 2 cm of pine, larch, oak and ash in the relative humidity range from 25% to 45%. The proposed model allowed to describe the process of moisture diffusion occurring in wood with a good fit for experimental research taking into account only three material parameters, i.e., the diffusion coefficient and two parameters describing the source term. Such a small number of model parameters allows to optimise the reverse methods of determining coefficients in wood and wood-based materials.

2. Model of Moisture Transport in Wood

2.1. Transport Equations

The moisture flow in wood in the hygroscopic range is mainly carried out by water vapour diffusion and sorption at the phase boundary between the air in the pores of the material and its skeleton. There is also a slight diffusion of bound water. However, it is much slower compared to the diffusion of water vapour in a porous structure [11] and thus negligible as mentioned in the introduction. For a one-dimensional case, without the influence of diffusion of bound water, the description of the phenomenon comes down to solving the following system of transport equations:
ε ρ v t = D v 2 ρ v x 2 m ˙ ,
ρ w d r y C b w t = m ˙ ,
in which ε is the porosity of the material [m3/m3], ρ v is the water vapour density [kg/m3], D v is the water vapour diffusion coefficient in the material in one of the orthotropic directions [m2/s], m ˙ is the source term related to the moisture sorption process [kg/m3s], C b w is the moisture content (mass of water divided by mass of dry material) [kg/kg] and ρ w d r y is the density of dry wood [kg/m3].
For Equations (1) and (2) uniform initial conditions are adopted:
ρ v ( x , t = 0 ) = ρ v 0 ,
C b w ( x , t = 0 ) = C b w 0  
and boundary conditions of the first type for Equation (1):
ρ v ( x = 0 , t ) = ρ v a i r ,
ρ v ( x = d , t ) = ρ v a i r ,
where ρ v 0 is the water vapour density corresponding to the initial relative air humidity in which the wood was conditioned [kg/m3], C b w 0 is the moisture content corresponding to the initial relative air humidity in pores [kg/kg], ρ v a i r is the water vapour density corresponding to the relative air humidity of the surrounding air [kg/m3], d is the sample thickness [m].
Such boundary conditions are correct for ambient air velocities greater than 3 m/s [18], which in the studies carried out in this paper was met (3.5–4 m/s).
For Equation (2), no boundary conditions are defined because the model does not assume bound water transport. Its changes are related only to the adsorption of water vapour, which takes into account the source term.
The change in average moisture content can be determined directly by measuring the mass of water that enters the material. It is also possible to determine the spatial distribution of the local moisture content in the sample. Isotopic or NMR methods are used [10,20]. The average moisture content is linked to the local moisture content by a relationship:
C ^ b w ( t ) = 1 d 0 d C b w ( x , t )   d x .

2.2. Moisture Equilibrium

Moisture adsorption is the process responsible for the bounding of water vapour molecules on the inner surface of the material. The amount of absorbed moisture is greater the higher the relative air humidity and less the higher the temperature [21]. By studying the changes in the absorbed moisture with respect to relative air humidity, a so-called sorption isotherm is obtained, which describes how much the material is capable of absorbing water at a constant relative humidity and a constant temperature to move to a state of equilibrium. This equilibrium is determined by comparing the chemical potential of water vapour and bound water. For vapour, this potential is well defined on the basis of thermodynamic considerations and has the form (e.g., in [21]):
μ v ( T , ρ v ) = R v   T   l n ( ρ v   R v   T p v s a t ) ,
in which R v is the gaseous constant of water [J/(kg K], T is temperature [K] and p v s a t is the saturation pressure of water vapour [Pa]. Several forms of this potential have already been proposed for water in porous materials [22,23,24], but it does not seem that any of them is particularly preferred. This article proposes a new form:
μ b w ( C b w ) = 1 a ( C b w m a x C b w 1 ) 1 b ,
in which C b w m a x is the maximum moisture content corresponding to fibre saturation point [kg/kg] while a and b are material constants expressed in [kg/J] and [-] respectively. By comparing Equations (8) and (9) and assuming that, for a state of equilibrium, the moisture content C b w becomes an equilibrium moisture content C b w e q an expression is obtained for the moisture content at which the equilibrium between moisture balance between the material and the environment is achieved:
C b w e q ( T , ρ v ) = C b w m a x ( a   R v   T   l n ( ρ v   R v   T p v s a t ) ) b + 1

2.3. Source Term

As sorption stops at equilibrium, its rate should be a function of the difference between the equilibrium and the current moisture content:
m ˙ = k   ( C b w e q C b w )
where k is the coefficient describing the rate of adsorption [kg/(m3s].
However, the study [8] showed that the adsorption rate cannot be constant. This article adopts the original form of a source term similar to the function describing the chemical potential of associated water, with the difference that the maximum moisture content was replaced by an equilibrium moisture content:
m ˙ = k 0 ( C b w e q C b w 1 ) n = k 0 ( C b w e q C b w ) n 1 C b w   n ( C b w e q C b w )
where k 0 is the coefficient of absorption rate [kg/(m3s] and n is the material coefficient [-]. This form of function may be interpreted as depending on the difference between the equilibrium and the current moisture content with a variable adsorption coefficient or as depending on the distance of the current moisture content to the equilibrium state.

2.4. Water Vapour Diffusion Coefficient

The coefficient D v is the effective coefficient of water vapour diffusion, which may be related to the coefficient of water vapour diffusion in air according to the relationship:
D v = D v a ψ
in which ψ is the dimensionless coefficient taking into account the resistance of the wood due to its porous structure and D v a is the coefficient of water vapour diffusion in air. The coefficient D v a adopted was as in [25] given as [26]:
D v a ( T ) = 2.16 · 10 5     ( T 273.15 ) 1.8

3. Verification of the Adopted Forms of Sorption Isotherm and the Source Term

The verification of the functions describing the sorption isotherm and source term was carried out on wood chips obtained from pine wood with density 459.6 kg/m3 (Figure 1). They were about 0.1 mm thick. This allowed, in the case of the determination of the sorption curve, to significantly shorten the testing time in relation to the testing of samples of several centimetres and, in the case of the source term, to minimise the influence of diffusion on the adsorption process. This type of testing is also used for other fibrous materials [27].
Initially, the samples were dried at 60 ° C . After drying, they were moved to a climatic chamber with constant climatic conditions T = 25 ° C and R H = 30 % . After the mass of the samples was stabilized (the mass changes of the samples were less than 0.002 g per 24 h), the relative humidity in the chamber was changed to R H = 50 % , and then the samples were weighed several times for 5 h. After the mass of the samples was stabilised again (about a week), the relative humidity in the chamber was changed again to R H = 70 % and the procedure was as in the previous stage. Once the mass was stabilised and the samples were weighed, they were placed in a exsiccator with relative air humidity R H 98 % and temperature T = 25 ° C . This made it possible to determine the close to maximum moisture content that can be achieved by the wood still in the hygroscopic range.
On the basis of the measurements of the change of masses, the moisture content at sequent moments were calculated according to the formula:
C b w ( t ) = m w ( t ) m w d r y m w d r y
where m w and m w d r y are the wet and dry masses of the sample. They were used to determine the parameters of the equation describing the sorption isotherm and the parameters of the model of moisture adsorption by wood.
The coefficients of Equation (10) were determined by the Levenberg-Marquardt algorithm implemented in the Matlab package. The sorption isotherms for the determined parameters are shown in Figure 2.
A very good fit of the model curve to the data obtained from the measurements confirmed that Equation (9) describes the potential of bound water very well. In addition, it can be seen that in the range R H = 20 ÷ 60 % the dependence of the moisture content on relative air humidity is almost linear, which has been used in the calculation (shown in Section 6) to determine the equilibrium moisture content.
At this stage of research, a verification of the adopted source term was also carried out. For such thin samples, the water vapour density is almost immediately equalised in the entire pore volume. With a constant relative humidity of the surrounding air, we obtain, therefore, a homogeneous in the whole volume and constant in time distribution of the equilibrium moisture content C b w e q = C b w e q ( ρ v a i r ) = c o n s t . This means that Equation (2), taking into account (12), is reduced to the form of:
C b w t = h 0 ( C b w e q = c o n s t . C b w 1 ) n ,   h 0 = k 0 ρ w d r y
There is not analytical solution for the above equation. The solution was obtained using the finite difference method (FDM) with an explicit scheme. The coefficients of the equation were determined so that the values obtained according to Equation (16) were as close as possible to those obtained during the experiment. This was achieved by finding the minimum of the function F by means of a domain search. This function was in the form of:
F ( x ) = k = 1 N ( C b w , k m , k C b w , k c , k C b w , k m , k ) 2
where N is the number of measurements taken, C b w , k m and C b w , k c are the moisture contents measured and calculated according to the model, and x is the vector of sought parameters h 0 and n .
The graph of the change of the moisture content versus time, obtained on the basis of the determined parameters, is shown in Figure 3.
In Figure 3 we can see a very good match between the values measured and obtained from the model in the humidity range from 30 % to 70 % . The material tested differed in terms of structure from the samples used in the main testing, but the results obtained confirm that the adopted form of the source term describes well the process of moisture adsorption by wood.

4. Description of Main Experiment

The main tests were performed on cubes with the following dimensions 2 cm × 2 cm × 2 cm, sourced from 4 types of wood (pine, larch, oak and ash), 36 pieces for each type. The samples were cut so that the growth rings are orthogonally aligned to their sides (Figure 4).
Initially, the samples were dried at 60   ° C to determine the mass of dry material. They were then transferred to a climatic chamber with constant climatic conditions T = 25   ° C and R H = 25 % . Once the masses were stabilized (the mass changes of samples were less than 0.002 g per 7 days), the samples were insulated on four sides in such a way as to force a one-dimensional moisture flow. In this way 12 samples were obtained each of the three anatomical directions: radial, tangential and longitudinal (Figure 5) and weighed again to calculate the mass of insulation used.
The wooden cubes prepared in this way were placed in the climatic chamber where, after the masses were stabilised again, the air humidity was increased to R H = 45 % . The change in the mass of the samples on the first day was measured every several hours. The time between sequent successive measurements was gradually extended to 4 days in the final stage of measurements. The measurements of mass change took 39 days (total approx. 3 months). The measured masses were used to determine the average moisture content in the samples according to the formula:
C ^ b w ( t ) = m w ( t ) m w d r y m w d r y
where m w and m w d r y are the current mass and mass of dry sample respectively. Then, for each direction, 9 cubes were selected for which the results were the least divergent (only these were taken into account in the calculations). For these nine cubes the average concentration of the samples was calculated for each time moment according to the formula:
C ¯ b w = i = 1 P C ^ b w i P
where P = 9 is the number of samples in the series.
For selected cubes the average apparent density of dry wood was calculated ρ w d r y . The density of the wood substance for all types of wood is almost identical and depending on the research method [28] oscillate around 1500 kg/m3 (assumed as in [29] citing as in [30]). This allows to calculate the approximate porosity of the tested wood from the formula:
ε = 1 ρ w d r y 1500 [ m 3 / m 3 ]
The average densities and approximate porosities of the tested wood depending on its species are shown in Table 1.
The data obtained in this way made it possible to determine the coefficients of the model proposed in the article.

5. Numerical Solution

By discretising Equations (1) and (2), using the implicit and centred finite difference scheme, we obtained:
ε ( ρ v ) j + 1 i ( ρ v ) j i Δ t = D v ( ρ v ) j + 1 i 1 2 ( ρ v ) j + 1 i + ( ρ v ) j + 1 i + 1 Δ x 2 m ˙ j + 1 i
ρ w d r y ( C b w ) j + 1 i ( C b w ) j i Δ t = m ˙ j + 1 i
where Δ t is the time increase, Δ x is the spatial increase, ( j ) is the index referring to the current moment in time for which the saturation values are known, ( j + 1 ) is the index referring to the future moment in time for which the saturation values are searching and ( i 1 ,   i ) and ( i + 1 ) are the indices referring to spatial nodes.
In order to obtain a solution to Equations (21) and (22), it is more advantageous to present in a matrix form:
A α β ( y j + 1   ) y j + 1   y j   Δ t + B α β ( y j + 1   )   y j + 1   = C α ( y j + 1   )  
where α , β = ρ , C , but A α β , B α β , C α are matrixes of coefficients depending on state variables (in the matrix B α β boundary conditions are taken into account), y j   = [ ( ρ v ) j     ( C b w ) j   ] is the vector of state variables in all spatial nodes for the ( j ) -th known point in time at which the initial conditions are taken into account, y j + 1   is the saturation vector in all spatial nodes for the ( j + 1 ) -th search moment in time.
The value of the sought solution y j + 1   is the argument A α β , B α β , C α . Therefore it is necessary to approximate the solution. This approximation was performed using the Newton-Raphson iterative method. The values of the variable in the ( k + 1 ) -th iteration in the ( j + 1 ) -th time step has a form [31]:
y j + 1 k + 1 = y j + 1 k ( Ψ y ( y j + 1 k ) ) 1 Ψ ( y j + 1 k )  
where:
Ψ ( y j + 1 k ) = A α β ( y j + 1   ) y j + 1   y j   Δ t + B α β ( y j + 1   )   y j + 1   C α ( y j + 1   )  
The condition for the end of the calculation is in the form Ψ ( S j + 1 k ) < ε where ε is a small value that determines the permissible error or assuming in advance the number of required iterations.
Based on the above algorithm, in the Matlab environment an original program solving the system of Equations (1)–(6) was written.

6. Optimization Procedure

The proposed model of moisture transport in wood has three parameters: two source term parameters k 0 , n and the water vapour diffusion coefficient D v . However, the tests lasted 39 days, which is too little time for the samples to reach equilibrium. The final moisture content C b w e q (equilibrium moisture content) that the samples would achieve after a very long time was also required. So eventually, the number of unknowns rose to four parameters.
The model coefficients were determined in such a way that the differences between the moisture content determined during the experiment and calculated using the model were as small as possible. The target function was minimised by the domain method. In the case of a single series, the target function F was in the form of:
F ( x ) = i = 1 N ( C ¯ b w m , i C ¯ b w c , i C ¯ b w m , i ) 2  
where N is the number of measuring moments, C ¯ b w m and C ¯ b w c are the moisture contents measured and calculated according to the model in time t i and x is a vector for the sought parameters: k 0 , n , ψ = D v a D v , Δ C b w = C b w e q C b w l a s t ( C b w l a s t —last measured value).
Matching errors were calculated to evaluate the fitting of model results to the measured ones. Local error according to the formula:
e l ( t i ) = | C ¯ b w m , i C ¯ b w c , i | C ¯ b w m , i  
and the global matching error according to the formula used in [32]:
e g = i = 1 N ( C ¯ b w m , i C ¯ b w c , i ) 2 i = 1 N ( C ¯ b w m , i ) 2  
The diagrams in Figure 6 show a very good fit of the model to the results obtained from measurements for all three directions of orthotropy. However, it can be seen that the coefficients of the source term shown in Table 2 differ slightly (although they should be the same because the source term is independent of direction). The reasons include measurement errors and heterogeneity of wooden samples. In order to determine the parameters k and n which would describe with a good approximation the moisture flow in all directions, the next step was to look for the best fit of the results for the three directions simultaneously.
Unfortunately in this case there are 8 coefficients to be determined. For this reason, this stage was divided into two parts. In the first one, for pre-determined diffusion coefficients and final moisture content, the parameters of the source term were determined by the domain method. In the second part, for the parameters set k 0 and n , the remaining coefficients for the best fit were sought. These calculations were repeated several times until the search values stopped changing. As shown in Table 3, for the results obtained in this way, the matching errors hardly changed compared to step one. This is due to the fact that the first stage already gives a very good approximation of the final results. The algorithm for determining coefficients of model is also shown in Figure 7.
In this way, calculations were made for all four sample types and the final results are presented in the next section.

7. Results and Discussion

As a result of calculations carried out in accordance with the scheme presented in Section 4, the parameter values of the proposed model of moisture transport in wood were obtained. They are presented in Table 4, Table 5, Table 6 and Table 7. For greater transparency and easier comparison of results, the diffusion coefficients are replaced by the vapour resistance factor for the successive orthotropic directions calculated according to the transformed formula (13):
ψ R = D v a D v , R ,   ψ T = D v a D v , T ,   ψ L = D v a D v , L
The results in graphical form are presented in Figure 8, Figure 9, Figure 10 and Figure 11. A model curve of the changes of moisture content C ¯ b w versus time and its measured values together with the standard deviation was plotted. The obtained results were also compared with the results for the multi-Fickian model (model m-F) presented in [10] for parameters: C 1 = 3.8 · 10 3 [1/s], C 3 = 15 [-], C 4 = 5.94 · 10 7 [1/s], C 21 = 3.58 [-], C 22 = 2.21 [-], C 23 = 1.59 · 10 3 [-], C 24 = 14.98 [-]. Standard deviation was calculated according to the formula:
s = i = 1 P ( C ^ b w i C ¯ b w ) 2 P 1  
where P is the number of samples in the series. Maximum and minimum values of standard deviation are presented in Table 8.
All the diagrams show a very good match between the curves obtained from the model and the measured values. This is confirmed by very small values of matching errors: local ( max ( e l ) = 0.77 % ÷ 1.71 % ) and global ( e g = 0.37 % ÷ 0.97 % ). In addition, the model’s matching to the measurements is similar for the entire course of the test, which indicates that the model describes the process well for moisture levels far from and close to equilibrium.
The results of several researchers presented by Time [33] show that for spruce at R H = 25 % the values of diffusion resistance in the longitudinal and transverse directions should be in the ranges ψ L 1.7 ; 5.0 and ψ T 11 ; 150 . Similarly Krabbenhoft and Damkilde [8] citing for Siau [13] show that ψ L 1.7 ; 2.8 and ψ T 56 ; 156 , although for higher humidity. Diffusion resistances obtained from the model presented in Table 5, Table 6, Table 7 and Table 8 do not differ much from the values determined by these researchers.
The results of the calculations also confirm the correct selection of the form of the source term. The coefficients k 0 and n calculated for the three directions separately (Table 2) are very similar. A small difference is due to the heterogeneity of individual samples. In turn, the values calculated for the three directions simultaneously also describe the process very well for each direction. An incorrect selection of the form of the mass source would take effect of different values of parameters k and n for each direction or in the impossibility of obtaining, with constant diffusion coefficients, a satisfactory match between the calculation and the measurement.
For comparison, Figure 8 also shows the curves determined on the basis of the multi-Fickan model presented in [10]. This model also describes the process quite well (especially for the radial and tangential directions), but the source term in this model has as many as seven parameters, which makes it very difficult to use it for the reverse method. The form of the source term proposed here has only two parameters, but the fit to the experimental data is even better. All this confirms the suitability of the proposed model for modelling moisture transport in wood.

8. Conclusions

The model of moisture transport in wood proposed in the paper describes very well the diffusion process in the range of moisture present in the studies. This is confirmed by very small matching errors not exceeding for local 2 % and for global 1 % . The source term adopted in the model, although describing a very non-linear adsorption process and in a large humidity range, has only two parameters. The small number of unknowns makes it easy to determine the water vapour diffusion coefficients in wood from the reverse method. The main difference between this model and most of the literature is that the diffusion coefficients are constant and independent of the humidity of the wood and the surrounding air. However, the main tests were performed for a fairly narrow humidity range, and although tests performed on thin samples give a high probability of model accuracy for higher humidity ranges, it is necessary to extend the tests to a larger humidity range to allow more far-reaching conclusions to be drawn.

Author Contributions

Methodology: A.K. and K.P.; Measurements: A.K. and K.P.; Writing—original draft preparation: A.K. and K.P.; Writing—review and editing: A.K. 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.

Data Availability Statement

he data are included in this article (Figure 3, Figure 8, Figure 9, Figure 10 and Figure 11)

Conflicts of Interest

The authors declare no conflicts interest.

References

  1. Kang, W.; Chung, W.Y.; Eom, C.-D.; Yeo, H. Some considerations in heterogeneous nonisothermal transport models for wood: A numerical study. J. Wood Sci. 2008, 54, 267–277. [Google Scholar] [CrossRef]
  2. Wang, J.; Mukhopadhyaya, P.; Morris, P. No title. J. Build. Phys. 2014, 37, 327–347. [Google Scholar] [CrossRef]
  3. Derome, D.; Rafsanjani, A.; Hering, S.; Dressler, M.; Patera, A.; Lanvermann, C.; Sedighi-Gilani, M.; Wittel, F.K.; Niemz, P.; Carmeliet, J. The role of water in the behavior of wood. J. Build. Phys. 2013, 36, 398–421. [Google Scholar] [CrossRef] [Green Version]
  4. 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 Sci. Technol. 2007, 41, 645–658. [Google Scholar] [CrossRef]
  5. Olek, W.; Perré, P.; Weres, J. Inverse analysis of the transient bound water diffusion in wood. Holzforschung 2005, 59, 38–45. [Google Scholar] [CrossRef]
  6. Perkowski, Z.; Świrska-Perkowska, J.; Gajda, M. Comparison of moisture diffusion coefficients for pine, oak and linden wood. J. Build. Phys. 2016, 41, 135–161. [Google Scholar] [CrossRef]
  7. Salin, J.-G. Mass Transfer from Wooden Surfaces and Internal Moisture Non-Equilibrium. Dry. Technol. 1996, 14, 2213–2224. [Google Scholar] [CrossRef]
  8. Krabbenhoft, K.; Damkilde, L. A model for non-Fickian moisture transfer in wood. Mater. Struct. 2004, 37, 615–622. [Google Scholar] [CrossRef]
  9. Frandsen, H.; Damkilde, L.; Svensson, S. A revised multi-Fickian moisture transport model to describe non-Fickian effects in wood. Holzforschung 2007, 61, 563–572. [Google Scholar] [CrossRef]
  10. Dvinskikh, S.V.; Henriksson, M.; Mendicino, A.L.; Fortino, S.; Toratti, T. NMR imaging study and multi-Fickian numerical simulation of moisture transfer in Norway spruce samples. Eng. Struct. 2011, 33, 3079–3086. [Google Scholar] [CrossRef]
  11. Hozjan, T.; Svensson, S. Theoretical analysis of moisture transport in wood as an open porous hygroscopic material. Holzforschung 2011, 65, 97–102. [Google Scholar] [CrossRef]
  12. Eitelberger, J.; Svensson, S. The Sorption Behavior of Wood Studied by Means of an Improved Cup Method. Transp. Porous Media 2011, 92, 321–335. [Google Scholar] [CrossRef]
  13. Siau, J.F. Transport Processes in Wood; Springer: Berlin/Heidelberg, Germany, 1984. [Google Scholar]
  14. 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. [Google Scholar] [CrossRef] [Green Version]
  15. Dinwoodie, J.M. Timber: Its Nature and Behaviour; Taylor & Francis: London, UK; New York, NY, USA, 2000. [Google Scholar]
  16. Kärger, J.; Valiullin, R. Mass transfer in mesoporous materials: The benefit of microscopic diffusion measurement. Chem. Soc. Rev. 2013, 42, 4172–4197. [Google Scholar] [CrossRef]
  17. Thybring, E.E.; Thygesen, L.G.; Svensson, S.; Hill, C.A.S. A critical discussion of the physics of wood–water interactions. Wood Sci. Technol. 2013, 47, 141–161. [Google Scholar] [CrossRef] [Green Version]
  18. Konopka, D.; Kaliske, M. Transient multi-Fickian hygro-mechanical analysis of wood. Comput. Struct. 2018, 197, 12–27. [Google Scholar] [CrossRef]
  19. Rosen, H.N. The influence of external resistance on moisture adsorption rates in wood. Wood Fiber 1978, 10, 218–228. [Google Scholar]
  20. Gezici-Koç, Ö.; Erich, S.J.F.; Huinink, H.P.; Van Der Ven, L.G.J.; Adan, O.C.G. Bound and free water distribution in wood during water uptake and drying as measured by 1D magnetic resonance imaging. Cellulose 2017, 24, 535–553. [Google Scholar] [CrossRef] [Green Version]
  21. Siau, J.F. Chemical potential as a driving force for nonisothermal moisture movement in wood. Wood Sci. Technol. 1983, 17, 101–105. [Google Scholar] [CrossRef]
  22. Bénet, J.-C.; Ramirez-Martinez, A.; Ouedraogo, F.; Cherblanc, F. Measurement of the chemical potential of a liquid in porous media. J. Porous Media 2012, 15, 1019–1029. [Google Scholar] [CrossRef]
  23. Abuku, M.; Ogura, D.; Hokoi, S. Measurement and modeling of water sorption isotherm for salt-laden autoclaved aerated concrete. Jpn. Arch. Rev. 2018, 1, 167–174. [Google Scholar] [CrossRef] [Green Version]
  24. Pavlík, Z.; Žumár, J.; Medved, I.; Ĉerný, R. Water Vapor Adsorption in Porous Building Materials: Experimental Measurement and Theoretical Analysis. Transp. Porous Med. 2012, 91, 939–954. [Google Scholar] [CrossRef]
  25. Zhang, Z.; Kim, C.S.; Kleinstreuer, C. Water Vapor Transport and Its Effects on the Deposition of Hygroscopic Droplets in a Human Upper Airway Model. Aerosol Sci. Technol. 2006, 40, 1–16. [Google Scholar] [CrossRef] [Green Version]
  26. Vargaftik, N.B. Tables on Thermophysical Properties of Liquids and Gases; Hemisphere: Washington, DC, USA, 1975. [Google Scholar]
  27. Kosiński, P.; Brzyski, P.; Duliasz, B. Moisture and wetting properties of thermal insulation materials based on hemp fiber, cellulose and mineral wool in a loose state. J. Nat. Fibers 2018, 17, 199–213. [Google Scholar] [CrossRef]
  28. Stamm, A.J. Density of Wood Substance, Adsorption by Wood, and Permeability of Wood. J. Phys. Chem. 1929, 33, 398–414. [Google Scholar] [CrossRef]
  29. Krzysik, F. Wood Science; PWN: Warsaw, Poland, 1975. (In Polish) [Google Scholar]
  30. Trendelenburg, R.; Mayer-Wegelin, H. Das Holz Als Rohstoff; Hanser: Munchen, Germany, 1955. [Google Scholar]
  31. Owen, D.R.J.; Hinton, E. Finite elements in plasticity: Theory and practice. Appl. Ocean Res. 1981, 3, 149. [Google Scholar] [CrossRef]
  32. Olek, W.; Rémond, R.; Weres, J.; Perré, P. Non-Fickian moisture diffusion in thermally modified beech wood analyzed by the inverse method. Int. J. Therm. Sci. 2016, 109, 291–298. [Google Scholar] [CrossRef]
  33. Time, B. Hygroscopic Moisture Transport in Wood. Ph.D. Thesis, Department of Building and Construction Engineering, Norwegian University of Science and Technology, Trondheim, Norway, 1998. [Google Scholar]
Figure 1. Thin wood samples used in the tests.
Figure 1. Thin wood samples used in the tests.
Materials 14 00017 g001
Figure 2. Sorption isotherm of the tested wood.
Figure 2. Sorption isotherm of the tested wood.
Materials 14 00017 g002
Figure 3. Change in moisture content in thin samples versus time.
Figure 3. Change in moisture content in thin samples versus time.
Materials 14 00017 g003
Figure 4. The cubic wood sample.
Figure 4. The cubic wood sample.
Materials 14 00017 g004
Figure 5. Insulated wood samples.
Figure 5. Insulated wood samples.
Materials 14 00017 g005
Figure 6. Diagrams of changes of moisture content for the pine samples.
Figure 6. Diagrams of changes of moisture content for the pine samples.
Materials 14 00017 g006
Figure 7. Algorithm for determining model coefficients.
Figure 7. Algorithm for determining model coefficients.
Materials 14 00017 g007
Figure 8. Changes in the moisture content of the pine samples versus time.
Figure 8. Changes in the moisture content of the pine samples versus time.
Materials 14 00017 g008
Figure 9. Changes in the moisture content of the larch samples versus time.
Figure 9. Changes in the moisture content of the larch samples versus time.
Materials 14 00017 g009
Figure 10. Changes in the moisture content of the oak samples versus time.
Figure 10. Changes in the moisture content of the oak samples versus time.
Materials 14 00017 g010
Figure 11. Changes in the moisture content of the ash samples versus time.
Figure 11. Changes in the moisture content of the ash samples versus time.
Materials 14 00017 g011
Table 1. Average apparent densities and porosities of the dry wood used in the tests.
Table 1. Average apparent densities and porosities of the dry wood used in the tests.
Type of WoodAverage Apparent Density of Dry Wood ρ w d r y   [ kg / m 3 ] Porosity   ε
[m3/m3]
Pine459.60.69
Larch558.40.63
Oak739.80.51
Ash680.80.55
Table 2. Coefficients obtained for the pine samples.
Table 2. Coefficients obtained for the pine samples.
Anatomical Direction Diffusion   Resistance   Factor   ψ [ ] Distance from Equilibrium Moisture Content
Δ C b w   [ kg / kg ]
Source Term Parameters
k 0 · 10 3   [ kg / ( m 3 s ] n [–]
Radial R 680.00508.44.5
Tangential T 820.00486.44.3
Longitudinal L 2.10.00437.44.4
Table 3. Matching errors of the final results and of the first calculation stage for the pine samples.
Table 3. Matching errors of the final results and of the first calculation stage for the pine samples.
Anatomical DirectionMatching Error [%]
Stage 1Final Results
Maximum Local ErrorGlobal ErrorMaximum Local ErrorGlobal Error
Radial R 0.860.420.890.43
Tangential T 1.00.370.990.38
Longitudinal L 1.360.621.460.63
Table 4. Model coefficients and matching errors for pine.
Table 4. Model coefficients and matching errors for pine.
Anatomical DirectionSource Term Parameters Vapour   Resistance   Factor   ψ [ - ] Error [%]
k 0 · 10 3   [ kg / ( m 3 s ] n [ - ] Max .   Local   e l , m a x Global   e g
Radial R 6.94.4 ψ R 650.860.42
Tangential T ψ T 831.00.37
Longitudinal L ψ L 2.01.460.63
Table 5. Model coefficients and matching errors for larch.
Table 5. Model coefficients and matching errors for larch.
Anatomical DirectionSource Term Parameters Vapour   Resistance   Factor   ψ [ - ] Error [%]
k 0 · 10 3   [ kg / ( m 3 s ] n [ - ] Max .   Local   e l , m a x Global   e g
Radial R 8.25.4 ψ R 1571.380.58
Tangential T ψ T 1141.130.59
Longitudinal L ψ L 2.40.900.49
Table 6. Model coefficients and matching errors for oak.
Table 6. Model coefficients and matching errors for oak.
Anatomical DirectionSource Term Parameters Vapour   Resistance   Factor   ψ [ - ] Error [%]
k 0 · 10 3   [ kg / ( m 3 s ] n [ - ] Max .   Local   e l , m a x Global   e g
Radial R 7.04.5 ψ R 1351.580.62
Tangential T ψ T 1650.770.42
Longitudinal L ψ L 6.81.100.44
Table 7. Model coefficients and fitting errors for ash.
Table 7. Model coefficients and fitting errors for ash.
Anatomical DirectionSource Term Parameters Vapour   Resistance   Factor   ψ [ - ] Error [%]
k 0 · 10 3   [ kg / ( m 3 s ] n [ - ] Max .   Local   e l , m a x Global   e g
Radial R 13.44.8 ψ R 1601.710.97
Tangential T ψ T 1781.710.95
Longitudinal L ψ L 5.41.220.58
Table 8. Maximum and minimum values of standard deviation.
Table 8. Maximum and minimum values of standard deviation.
Type of WoodAnatomical DirectionStandard Deviation [10−2 kg/kg]
MinimumMaximum
PineRadial R 0.0470.076
Tangential T 0.0970.146
Longitudinal L 0.0580.107
LarchRadial R 0.0610.096
Tangential T 0.0340.110
Longitudinal L 0.0300.059
OakRadial R 0.0380.090
Tangential T 0.0520.123
Longitudinal L 0.0520.088
AshRadial R 0.0690.106
Tangential T 0.0480.093
Longitudinal L 0.0660.097
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kucharczyk, A.; Pawlik, K. Modelling and Experimental Study of Moisture Transport in Wood, Taking into Account Diffusion and the Accompanying Adsorption of Water Vapour by Cell Walls. Materials 2021, 14, 17. https://doi.org/10.3390/ma14010017

AMA Style

Kucharczyk A, Pawlik K. Modelling and Experimental Study of Moisture Transport in Wood, Taking into Account Diffusion and the Accompanying Adsorption of Water Vapour by Cell Walls. Materials. 2021; 14(1):17. https://doi.org/10.3390/ma14010017

Chicago/Turabian Style

Kucharczyk, Andrzej, and Kamil Pawlik. 2021. "Modelling and Experimental Study of Moisture Transport in Wood, Taking into Account Diffusion and the Accompanying Adsorption of Water Vapour by Cell Walls" Materials 14, no. 1: 17. https://doi.org/10.3390/ma14010017

APA Style

Kucharczyk, A., & Pawlik, K. (2021). Modelling and Experimental Study of Moisture Transport in Wood, Taking into Account Diffusion and the Accompanying Adsorption of Water Vapour by Cell Walls. Materials, 14(1), 17. https://doi.org/10.3390/ma14010017

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