An experimental study on change of gas permeability depending on pore structures in three species (hinoki, Douglas fir, and hemlock) of softwood

The purpose of this study is to analyze the pore structure for the heartwood, intermediate wood, and sapwood sections in three species of softwood (hinoki, Douglas fir, and hemlock) and to investigate the correlation of gas permeability depending on pore structure. For this study, gas permeability and pore size were measured by capillary flow porometry, and classification of a novel method was performed to determine the type of pores (through pores, blind pores, and closed pores) based on International Union of Pure and Applied Chemistry (IUPAC). Gas permeability, through pore porosity, and pore size increased from heartwood to sapwood. The results of multiple regression analysis showed that through pore porosity, mean pore size, and bulk density were significant factors affecting gas permeability.


Introduction
The permeability of wood greatly influences drying rate because free water moves to the surface layer by capillary flow at the beginning of drying [1]. High-temperature and high-frequency drying methods generate water vapor that can move rapidly through the capillary passage [2].
In recent research, permeability of wood is recognized as an important parameter in various applications such as heat-treated wood [3][4][5], nanomaterial impregnation [6][7][8] and use of the transverse surface of wood as a porous sound-absorbing material [9][10][11] Applying an appropriate drying schedule can improve the permeability of wood [12,13]. Impregnation of aqueous metal nanoparticles into heat-treated solid woods can improve its permeability because some of the wood's extractives dissolve, producing perforations and pits [14]. Further, heat treatment of wood increases the permeability of wood due to loss of bound water, shrinkage at low temperatures, and physicochemical changes of cell walls at high temperatures [15].
As the diameter of a tree reaches a certain level during the growth process, the color between the inside and the outside of the cross-section appears different. The dark part of the inside is called heartwood and the light part of the outside is called sapwood [18].
Sapwood is a living tissue that acts as a pathway for water and nutrients. It contains a high level of moisture and comprised a soft texture. But heartwood is composed of dead cells with no physiological dysfunction and a low level of moisture content [19].
Over the years, many researches have reported that sapwood has greater gas permeability than heartwood [20][21][22][23] due to the difference in chemical and anatomical characteristics [24].
From a chemical point of view, cellulose and holocellulose (a combination of cellulose and hemicellulose) contained in heartwood are generally less abundant than those of softwood [25][26][27][28]. It is also known that heartwood has a large amount of extractives [28][29][30]. The content of the extractives has a significant effect on permeability, resulting in less permeability of heartwood than sapwood and difficulty in drying, impregnating, and using the wood pulp [25,26,28,30,31].
With respect to anatomical characteristics, the main cause of the permeability difference in softwood is extent of closure of the bordered pit pair in heartwood and of the extract deposits accumulated on the pit membrane [32]. Increased permeability in sapwood is attributed to the size and frequency of resin canals and parenchyma rays [33].
The permeability of wood is related to pore structure, which has been investigated using various methods.
Schneider and Wagner reported that the mean pore size of heartwood was smaller than that of softwood as measured by mercury intrusion porometry [34]. Plötze and Niemz reported pore distribution and pore size of 24 wood species using a mercury intrusion porosimeter and classified them as macro-pore, meso-pore, and micro-pore according to pore diameter [35]. Kang et al. [3] reported change in pore size of heat-treated malas (Homalium foetidum) by capillary flow porometry. Jang et al. reported that the pore size varies depending on method, between mercury intrusion porosimetry and capillary flow porometry [36]. Jang et al. [37] classified three pore shapes (through pore, blind pore, and closed pore) in heartwood, intermediate wood, and sapwood of cross-sectional yellow poplar and measured their contents. The results showed that pore size and through pore porosity increased from heartwood to sapwood [37]. In the same way, Jang and Kang analyzed the pore structure of heat-treated woods that are classified into three pore types based on IUPAC. They reported that, as the heat treatment temperature increased, the through pore porosity and pore size also increased [5].
Unlike hardwood, softwood is composed of simple structures referred to as tracheid, and its heartwood is more difficult to dry and impregnate with chemical than that of sapwood, which is commonly known to cause pit aspiration of the tracheid [19]. The tracheid of softwood is considered a type of through pore, while pit aspiration of the tracheid is considered closed pore in accordance with IUPAC's morphological classification.
Pore classification based on IUPAC between heartwood and sapwood was performed in a previous study [37]. However, in the previous study, pore structure was classified for only one species of hardwood (yellow poplar), and no comparison of species was performed. For such a classification to be useful in wood, it needs to be tested in various species. Also, the studies did not clearly demonstrate statistical explanations on how the gas permeability difference correlates with the type of pores in various species of wood. Therefore, this study selected three species of softwood [hinoki (Chamaecyparis obtusa (Siebold & Zucc.) Endl, Douglas fir (Pseudotsuga menziesii), and hemlock (Tsuga heterophylla)]. In addition, statistical analysis was performed to investigate the correlation between air permeability and pore structure of wood.
Specifically, this study analyzed pore size, gas permeability, and pore type in accordance with IUPAC's morphological classification of pores in three species of softwood. In addition, this study performed statistical analysis on correlations and associations among pore size, pore type, and gas permeability by Pearson correlation analysis and multiple regression. The logs were cut into cross-sections with thickness of 1 cm. In hinoki, the entire diameter was 20 cm, and the diameter of the heartwood was 14 cm. The entire diameter of the Douglas fir was 32 cm, and the diameter of the heartwood was 20 cm. In the hemlock, the entire diameter was 29 cm, and the diameter of the heartwood was 20 cm. From the cross-sectional discs, 10 cylindrical specimen (30 mm diameter × 10 mm thickness) of heartwood, intermediate wood, and sapwood were cut. Intermediate wood samples were distinguished by interfacing with approximately 50% of the area of the heartwood and sapwood. The cylindrical specimens were dried at 40 °C using a dryer for 1 week and controlled within a 10% moisture content (MC) to exclude the influence of MC between heartwood and sapwood to focus on pore structure.

Specimen preparation
To measure the true density of the wood substance, internal voids must be ignored. Wood samples from the sapwood and the heartwood were pulverized into fine wood powder [34,37,38].

Scanning electron microscope (SEM) imaging for cross-sections of specimens
To observe the pore morphology of a specimen, its cross-sections were cut into dimension of approximately 5 mm (radial) × 5 mm (tangential) × 4 mm (longitudinal), and the surfaces were smoothed using a microtome (model: HM400S, Microm GmbH, Germany). Further, specimens were coated with gold ions by an ion sputtercoater (model: SCM, Emcrafts, Korea) and observed at an acceleration voltage of 13 kV and 400× magnification using scanning electron microscopy (SEM; model: Genesis-1000, Emcrafts, Korea).

Measurement of bulk density and true density
Bulk density of the cylindrical specimens was measured by the dimension method. True density of cylindrical specimens and sawdust was measured by gas pycnometer (model: PYC-100A-1, Porous Material Inc., USA) [5,37,38]. Before measuring the specimen, the gas pycnometer was calibrated to a standard stainless-steel cylindrical sample from the manufacturer, the volume and mass of which were 15.4346 cm 3 and 41.1767 g, respectively. The volumes of the reference chamber and sample chamber were 35.5310 cm 3 and 33.5224 cm 3 , respectively. For convenience, the true density of sawdust intermediate wood was calculated as the average of heartwood and sapwood.

Measurement of gas permeability
Gas permeability was measured using a capillary flow porometer (model: CFP-1200AEL, Porous Materials, Inc., USA). To prevent air leakage from edge of a cylindrical specimen, the side surface was sealed with a silicon O-ring. Increasing pressure slowly, air flow through the pore was measured, and the Darcy permeability constant was calculated from Eq. (1) shown below: where C is the Darcy permeability constant; F is the flow; T is the sample thickness; V is the viscosity of air; D is the sample diameter; and P is the pressure.

Measurement of pore diameter
Pore diameter of specimens was measured using a capillary flow porometer (model: CFP-1200AEL, Porous Material Inc, USA). The ASTM F316-03 [39] method was used to selectively measure the constricted part of a through pore related with permeability and is widely used for analysis of pore structure of various porous materials [40]. It is also a possible way to selectively measure only through pores in line with permeability of wood [36,37]. As shown in Fig. 2, in the process of gas permeability estimation by flow rate under pressure, the graph of flow rate versus pressure is called a 'dry curve' . Next, the specimens were wetted in Galwick solution that has an extremely low surface tension (surface tension: 0.159 mN/m) and non-volatility to penetrate well into the cavity. Pressure was measured sequentially from the moment when the flow rate of the gas was detected with solution extrusion, and the pressure at the time of initial flow increase is called the 'bubble point' . In the wet state, the graph of flow rate versus pressure is called a 'wet curve' . Here, an imaginary curve having a slope of 1/2 of that of the 'dry curve' can be drawn and is called a 'half-dry curve' . The point where the 'wet curve' and 'halfdry curve' meet is called 'mean flow pore pressure' , and the pore size is determined by Eq. (2) as below: where D is the limiting diameter, τ is the surface tension, p is the pressure, and C is the constant of 2860 when p is in Pa, 2.15 when p is in cmHg, and 0.415 when p is in psi units.

Porosity of through pore, blind pore, and closed pore types
Following IUPAC [41], pore types of solid porous materials can be classified as through pores, blind pores, and closed pores (Fig. 3). Hence, this study measured the content of each type of pores in the three species of wood as in previous studies [5,37]. Specifically, content was measured for through pore porosity (φ through ), blind pore porosity (φ blind ), and closed pore porosity (φ closed ) of the cylindrical specimens using the following method. Total porosity φ total (%) was calculated from Eq. (3) measuring true volume (V true ) of the sawdust by gas pycnometry and bulk volume (V bulk ) of the cylindrical specimens by the dimensional method: Pore volume of the cylindrical specimen by gas pycnometer was obtained as the sum of through pore volume (V through ) and blind pore volume (V blind ). Thus, the sum of through pore porosity (φ through ) and blind pore porosity (φ blind ) was calculated as in Eq. (4): Finally, from the results of Eqs. (3) and (4), closed pore porosity (φ closed ) was obtained with Eq. (5): To distinguish through pore porosity (φ through ) from blind pore porosity (φ blind ), this study wetted cylindrical specimens in Galwick solution (surface tension: 0.159 mN/m) and intruded the liquid into the pores utilizing a vacuum pump. After that, the specimens were placed in a chamber of which both sides were sealed with O-rings, and the air pressure was increased in the longitudinal direction so that only Galwick in the through We obtained the mass of Galwick solution in blind pores by measuring the difference between specimen mass after extrusion (M after extrusion ) and dried specimen mass (M dried ). Also, the volume of blind pore (V blind ) was calculated as the difference between specimen mass (M after extrusion − M dried ) and specific gravity of Galwick solution (ρ Galwick ) using Eqs. (6, 7, and 8):

Statistical analyses
This study applied a statistical approach to corroborate reliability of the experimental results.
Physical properties (pore type, pore size, and density) of wood can affect gas permeability. Among the three types of pores for solid porous materials defined by IUPAC, only through pores can permit fluid inflow. Thus, through pore porosity may have a significant effect on gas permeability [5,37]. Moreover, as pore size increases, the flow of fluid generally increases; as density increases, the flow of fluid generally decreases. In short, pore size, bulk density, and through pore porosity are significant factors affecting gas permeability. Thus, this study statistically investigated the effects of through pore porosity, mean pore size, and bulk density on gas permeability.
Specifically, this study compared the average gas permeability, pore size, and porosity of three parts of heartwood, intermediate wood, and sapwood using analysis of variance (ANOVA).
In addition, this study analyzed statistical correlation by Pearson correlation analysis to understand univariate associations among the variables. Furthermore, this study utilized multivariate regression to analyze effects of physical factors on gas permeability to enable investigation of factor-specific effects on dependent variables after controlling the other factors. Multivariate regression has another advantage in that it provides the simultaneous effects of several variables on the dependent variable. The multivariate regression model was used to display the effects of through pore porosity, mean pore size, and bulk density on gas permeability as displayed in the equation below: where Y is the gas permeability; X 1 is the through pore porosity; X 2 is the mean pore size; X 3 is the bulk density; D 1 is the wood dummy variable (1: if the wood is hinoki, otherwise 0); D 2 is the wood dummy variable (1 (9) if the wood is Douglas fir, otherwise 0); α 0 is the constant (intercept term); and ε is the residuals (error term). This study uses Eq. (9) as an ordinary least square (OLS) regression model that formulates a linear function. Thus, α 0 in Eq. (9) represents a constant (intercept term) of the linear function model, while ε indicates residuals (error term) of the regression model. Basically, regression analysis such as Eq. (9) assumes that the expected value of the error term is zero [E(ε i ) = 0], resulting in an estimated value of 0 (zero) for ε.
If estimated coefficients of β 1 , β 2 and β 3 in Eq. (9) are significant, the three physical properties (X 1 , X 2 , and X 3 ) substantially affect gas permeability. Wood dummy variables (D 1 and D 2 ) were included in Eq. (9) to control the effect of wood species. Specifically, D 1 has a value of 1 if the wood is hinoki and 0 (zero) if the woods are other species. Likewise, D 2 has a value of 1 if the wood is Douglas fir and 0 (zero) if the woods are other species. Therefore, hemlock has a value of 0 (zero) for both D 1 and D 2 , and those two dummy variables (D 1 and D 2 ) control the effect of the difference among the three species of wood. Further, this study analyzes VIF (variance inflation factor) values of all independent variables in Eq. (9) to assess the multi-collinearity among the independent variables. Multi-collinearity is not serious if the estimated VIF value does not exceed 10.
All sample data were analyzed using both Pearson correlation analysis and regression analysis (n = 90). Also, measured values of all variables were winsorized at 2% of both sides to minimize the effect of extreme values.
All statistical analyses were performed using IBM SPSS statistics v25 software (IBM Corp., Armonk, NY, USA).

SEM observation
The cross-sectional surface of each specimen is presented in Fig. 4. Earlywood and latewood were observed in every specimen. However, it was difficult to distinguish characteristics of heartwood and sapwood in SEM images.

Results of bulk density and true density
The average bulk densities of hinoki were 0.47 g/cm 3 (standard deviation, SD = 0.01) of heartwood, 0.45 g/ cm 3  The average true densities of hinoki sawdust were 1.53 g/cm 3 of heartwood and 1.43 g/cm 3 of sapwood; those of the cylindrical specimens were 1.35 g/cm 3 (SD = 0.03) of heartwood, 1.35 g/cm 3 (SD = 0.02) of intermediate wood, and 1.39 g/cm 3 (SD = 0.03) of sapwood. The average true densities of Douglas fir sawdust were 1.42 g/cm 3 of heartwood and 1.53 g/cm 3 of sapwood, while those of the cylindrical specimen were 1.08 g/cm 3 (SD = 0.01) for heartwood, 1.26 g/ cm 3 (SD = 0.05) for intermediate wood, and 1.48 g/cm 3 (SD = 0.23) for sapwood. The average true densities of hemlock sawdust were 1.52 g/cm 3 for heartwood and 1.40 g/cm 3 for sapwood, while those of the cylindrical specimen were 1.26 g/cm 3 (SD = 0.10) for heartwood, 1.25 g/cm 3 (SD = 0.01) for intermediate wood, and 1.28 g/cm 3 (SD = 0.04) for sapwood.
It seems that the difference of true density between heartwood and sapwood in sawdust is caused by difference in extractive content. Further, true density tends to decrease when moving away from heartwood toward sapwood mainly due to the decrease in closed pore porosity.

Results of gas permeability
As shown in In all three softwoods, the permeability of sapwood was greater than that of heartwood, in agreement with prior studies [2,20,22]. The morphological pore characteristics of heartwood and sapwood were not significantly different, while the variance of gas permeability was high. SEM could observe only the cross-sectional surface of tracheid and could not observe the internal pathway. Hence, it is not feasible to predict the permeability of wood in the longitudinal direction by morphological observation of a cross-sectional surface of wood by SEM.

Results of through pore diameter
The pore size measured by the capillary flow porometer consists of the constricted part of the through pore. In a cross-section of softwood, the through pore is tracheid, and the constricted part is margo pit membrane.
As shown in Fig. 6, maximum pore sizes of hinoki were 0.58 μm (SD   [42][43][44]. Porosity analysis (through pore, blind pore, and closed pore) As shown in Fig. 7 Pit aspirations of tracheid in the heartwood are in the form of blind pore or closed pore, resulting in lower through pore porosity in heartwood than sapwood. Table 1 presents the results of Pearson correlation analysis among the variables. Gas permeability showed a highly positive relationship with through pore porosity. (0.709) at a 1% level of significance and a negative trend with bulk density (− 0.258) at a 5% level of significance.

Correlation and regression results
These results indicate that, as through pore porosity increases, gas permeability also increases. However, as bulk density increases, gas permeability decreases. No significant relation was found between mean pore size and gas permeability. Table 1 shows that the multi-collinearity problem in regression analysis was not significant because the correlations among the independent variables (porosity, pore size, and density) were statistically weak or not significant. Table 2 shows the regression results with Eq. (9). Panels A, B, and C present estimated results of Eq. (9) with each independent variables, and Panel D shows estimated results of regression with all independent variables. All estimated coefficients of porosity (β 1 ) and pore size (β 2 ) were positively significant at a 1% level in Panels A, B, and D. In addition, estimated coefficients of density (β 3 ) were negatively significant at a 1% level in Panels C and D. These results suggest that an increase of through pore porosity and mean pore size as well as a decrease of bulk density substantially increase gas permeability. All F-values of Panels A, B, C, and D were significant at a 1% level, indicating goodness-of-fit for the estimated regression results. Most adj. R 2 s were above 60%, especially that of 78.40% in Panel D. Hence, the explanatory power of the regression model was also acceptable. All VIFs (variance inflation factor) in Panel D were less than 10, which suggests multi-collinearity was not significant.
The results of Table 2 represent the three species (hinoki, Douglas fir, and hemlock) of softwood and should not be applied to other softwood species. They can be developed into a forecasting model after including and analyzing more species of softwood.

Conclusions
Summary of the results of the experiments on permeability and pore structure of heartwood and sapwood is as follows: 1. SEM was used to observe the cross-section of the wood and confirmed that the pore was slightly smaller in heartwood compared to sapwood. However, the difference was not significant, so we could not analyze detailed pore structure from heartwood to sapwood. 2. The gas permeability of hinoki, Douglas fir, and hemlock increased from heartwood to sapwood. 3. The pore structures of softwood was classified into through pore, blind pore, and closed pore following IUPAC definitions and observed using capillary flow porometer and gas pycnometer. Through pore porosity tended to increase from heartwood to sapwood. 4. Capillary flow porometry was used to evaluate pore size on the cross-sectional surface of the wood and showed that maximum pore size and mean flow pore size of sapwood were greater than those of heartwood. 5. Through pore porosity, mean pore size, and bulk density were the main factors affecting gas permeability in multivariate regression analysis.
Type of pore is related to permeability of the wood, and the quantity of pores affect drying, processing, and modifying of woods.