Modifying wood by high intensive microwave pretreatment method is widely researched for the fabrication of wood-based nanocomposites, but the temperature uniformity and energy efficiency of microwave pretreatment have not reached the ideal state. In this study, the pretreated wood in rectangular cavity by high intensive microwave is theoretically studied by the finite element method based on the Maxwell electromagnetic field equations and the heat and mass transfer theory. The results show that the temperature uniformity and energy efficiency are related to the microwave feeding modes. Compared with the single-port and the two-port feeding mode, the four-port feeding mode is the best case on temperature uniformity and energy efficiency. The optimized parameters of cavity to pretreatment wood are achieved, which are that the height of cavities is between 0.08 m and 0.11 m in the four-port feeding mode when the thickness of wood is 0.06 m.
1. Introduction
Microwave pretreatment is a technique to improve process efficiency and product performance of wood in the forest product industry. High-intensity microwave pretreatment for wood modification originated from Australia [1, 2] at the end of last century, in which water in wood absorbed enough energy from the transient electromagnetic wave and turned into steam. High internal steam pressure destructed the weaker elements of wood and the microvoids in the radial-longitudinal planes of wood were created [3, 4], which puffed wood into loose structure, increased wood permeability, and improved the drying process. Some studies reported that several properties of wood were improved after microwave pretreatment. In Australia, Torgovnikov and Vinden’s studies [5–7] showed that, after high-intensity microwave pretreatment, the permeability of the wood increased up to 14-fold, the preservative penetration increased 10–14 times, the density of the wood reduced, and the elasticity modulus decreased 12–17%. As a result, a kind of wood with high permeability, low density, and good sound insulation was made, which was called puffing wood. And this treated wood could be used for new composites, such as the fabrication of wood-based nanocomposites and new functional materials [5–12]. To make the timber inorganic nanocomposites, direct infusion method is usually used, in which the composite compounds impregnate in wood cells cavity or walls through direct physical infusion method to make inorganic nanocomposites, such as various water-borne inorganic complex salt preservatives (boron preservative, boron-clay composite preservative, boron-silicon composite preservatives, basic zinc carbonate, copper oxide, sodium silicate, etc.) [13, 14], a variety of water-borne flame retardants (phosphorus-nitrogen boron-based flame retardants, amine-modified silicone-based flame retardants, phyllosilicates, etc.) [15–19]. This method of making the timber inorganic nanocomposites requires wood dried uniformly with large gap. Microwave pretreatment of wood can meet these requirements. As we know, the more even temperature profiles are, the better modification effect of wood will be. However, the most serious problem of the existing research is the nonuniformity of wood microwave modification. To solve this problem, we focused on related microwave pretreatment devices and the development of the microwave cavity with uniform distribution of microwave energy. The electromagnetic field distribution model and heat-mass transfer model are built for optimizing high-intensity microwave wood pretreatment process.
2. Material and Methods
In order to simulate the wood microwave heating process, some related models are created. There are two causes for the heat within wood during microwave heating. One part is the microwave absorbed by wood, which is described by absorption power of microwave, and the other is the heat transfer between each part of wood. Once the heat absorbed, heat transfer model, and its initial and boundary conditions are established, the temperature distribution within wood during microwave heating can be figured out. In the process of the simulation of high-intensity microwave heating, two parts are studied which are the distribution of electromagnetic wave and wood temperature in rectangle waveguide resonant cavity with different feeding modes. In order to simplify the calculation, some reasonable assumptions of the model are made as follows: the temperature and moisture are evenly distributed in the wood before heating. The air and wood in the cavity remain stationary during the heating process. Thirdly, wood is regarded as a homogeneous material. Finally, the moisture content and volume of wood remain unchanged while heating.
Finite element analysis is used to solve electromagnetic field and temperature distribution within wood. The microwave frequency is 0.915 GHz, and the length, width, and height of rectangular resonator cavity are 0.4953 m, 0.24765 m, and 0.08 m, respectively. The total input power of microwave excitation sources is set at 40 kW, which is averagely allocated to each excitation source. The waveguide of feeding port is a standard waveguide (Chinese national standard BJ9 model, length 0.24765 m, width 0.12382 m). Initial moisture content of wood is set at 60%, the length, width, and height of the wood are set at 0.74295 m, 0.2 m and 0.06 m, respectively, and the relative dielectric constant and dielectric loss factor are 7.40 and 1.42, respectively. The initial temperature and the density are 25°C and 570 kg/m3, respectively. The thermal conductivity λ (W/m °C) and specific heat C (J/kg °C) depended on the following equations:
(1)λ=0.23[1-0.72×(25-T)100],C=2650×(1+T100)0.2.
3. Theory
While building the electromagnetic field and temperature distribution models, thermodynamics and Maxwell’s electromagnetic theory are the foundation of the simulation.
3.1. The Distribution of the Electromagnetic Wave and Microwave Energy
According to Maxwell’s equations, the electromagnetic field satisfies the following equations during microwave heating progress in the rectangular cavity:
(2)∇×E→=-∂B→∂t,∇×H→=J→+∂D→∂t,∇·B→=0,∇·D→=ρc,
where E→ is the electric field intensity, H→ is the magnetic field intensity, B→ is the magnetic induction intensity, D→ is the electric displacement, J→ is the current density, and ρc is the free charge density, wherein J→=σE→,D→=εE→,B→=μH→,σ is the conductivity, ε is the dielectric constant, and μ is the magnetic permeability. According to (2), the distribution of the electromagnetic field inside the wood can be calculated [20].
In this study, the microwave excitation source is time-harmonic electromagnetic fields E→10, where E→=E→(xi→+yj→+zk→)e-iωt,H→=H→(xi→+yj→+zk→)e-iωt, B→ and D→ can be expressed as functions E→ and H→:
(3)∇×E→(xi→+yj→+zk→)=iωμH→(xi→+yj→+zk→),∇×H→(xi→+yj→+zk→)=(σ-iωε)E→(xi→+yj→+zk→)=-iωε*E→(xi→+yj→+zk→),
where ε* is the complex permittivity, which can be expressed as ε*=ε′+iε′′. Equation (4) can be obtained from (3) [20]:
(4)∇(E→·∇ε*ε*)+∇2E→+k2E→=0.
By solving (4), the electric field distribution in wood can be obtained. In this study, we assumed that the wood dielectric is constant value along the direction of the electric field and the first term of (4) is zero, so (4) can be simplified as follows:
(5)∇2E→+k2E→=0,
where k depends on the dielectric properties of wood and its expression is as follows:
(6)k=α+iβ,α=2πfcε′(1+tan2δ+1)2,β=2πfcε′(1+tan2δ-1)2,
where f is the frequency of the microwave radiation and tanδ is the wood loss factor; its value can be calculated from the following (7):
(7)tanδ=ε′′ε′.
In order to calculate the electric field distribution in wood, the solutions of (5) are obtained by the following boundary conditions:
(8)n×[E→1(xi→+yj→+zk→)-E→2(xi→+yj→+zk→)]=0,(9)n×[H→1(xi→+yj→+zk→)-H→2(xi→+yj→+zk→)]=0,
where the subscript 1 indicates air and the subscript 2 indicates wood. Equation (9) can also be expressed as a function of electric field intensity through the following expression:
(10)∂E→(x)∂xi→+∂E→(y)∂yj→+∂E→(z)∂zk→=iμ0ωH→(xi→+yj→+zk→).
The analytical solution of (5) is as follows:
(11)E=A1eikx2+y2+z2+B1e-ikx2+y2+z2,
where A1 and B1 depend on the boundary conditions (8) and (9). The electric field in wood can be used to calculate the electromagnetic field energy density by Poynting theorem:
(12)Q=12ωε0ε′′E→·E→*,
where ε0 is the permittivity of vacuum and E→* is the complex conjugate of the electric field.
3.2. Heat Transfer Model
In microwave pretreatment process, the microheat balance of wood is described by the following equation [20]:
(13)ρCP∂T∂t=∇·(kT∇T)+Q,
where ρ, CP, and kT are density, specific heat, and thermal conductivity of wood, respectively, and Q is electromagnetic energy density in the wood.
Boundary conditions of (13) are as follows [21]:
(14)t=0,T=Tini,0≤x≤L,t>0,-kT∂T∂x=h(T-Ta)+Lvapkm′(Cw,s-Cequi),000000000000000000000000000000000000000x=L,
where x=0 indicates the wood surface, L is the width of the wood; h is the heat transfer coefficient of the inner wood; Ta and Tini are air and wood initial temperature; Lvap is the heat of vaporization of the water; km′ is the mass transfer coefficient; Cw,s is the moisture content of the wood surface; Cequi is the equilibrium moisture content of the air.
3.3. Mass Transfer Model
To calculate the water concentration profile, the microscopic mass balance can be described by the following equation [22]:
(15)∂Cw∂t=∇(Dw∇Cw),
where Dw is the effective diffusivity coefficient of the wood.
The boundary conditions of (15) are as follows:
(16)t=0,Cw=Cw,ini,0≤l≤L,t>0,-Dw∂C∂x=km′(Cw,s-Cegui),l=0,l=L,
where Cw,ini is the wood initial moisture content.
4. Results and Discussion4.1. Effect of Feeding Mode on Temperature Distribution and Energy Efficiency
The temperature distribution and energy efficiency of each feeding mode are obtained by computational simulation. The spatial distribution of temperature in feeding modes is studied (Figure 1). After being pretreated by microwave with the same power and time, the wood temperature of the single-port, the two-port, and the four-port feeding mode gradually increases. And this is due to the microwave distribution in the wood. High-temperature region distribution of the four-port feeding mode is dispersed; the temperature uniformity is better than the single-port and the two-port, feeding mode. But the temperature uniformity of the single-port and the two-port feeding mode could not be derived directly from the figure. As it can be seen from the three kinds of feeding modes, the highest-temperature region is not on the surface but inside the wood. And this is caused by the penetration of microwave, which means the microwave passes through the surface of the wood and occurs in strong interaction with the inner wood. The temperature distribution of single-port feeding mode appeared to be a “U,” which is because the microwave feeding direction is oblique; the microwave energy distribution is uneven in the y-axis of wood. Figure 1 shows that temperature distribution in midsection is similar to the multisections. And most of the microwave energy is absorbed by the midsection of the wood. Thus, the midsection temperature distribution is studied specially.
The multisections temperature (°C) distribution within wood heated by microwave with different radiation methods. (a), (b), and (c) describe the single-port feeding mode, the two-port feeding mode, and the four-port feeding mode, respectively.
Figure 2 shows that temperature in the bottom part of wood is the highest in the single-port feeding mode. This is due to the deeper penetration of microwave in the wood when there is only one feeding port. In the two-port feeding mode, temperature distribution is symmetric along xy and xz plane; the highest temperature region is in the center. And this is due to the fact that the microwave interference is formed when microwaves in both directions with the same phase position come cross each other. In the four-port feeding mode, temperature distribution is disperse; high temperature region is distributed uniformly. And this is because there is microwave energy from all directions in the four-port feeding mode. Figure 2 shows the temperature distribution uniformity roughly; in order to analyze the temperature uniformity quantificationally, temperature coefficient of variation of the midsection temperature distribution within wood is studied. While energy efficiency is the ratio of energy absorbed in wood and whole input energy, which can reflect energy utilization and dissipation of resonator, the energy efficiency is very important on the performance evaluation of the cavity, which is taken into consideration too.
The midsection temperature (°C) distribution within wood heated by microwave with different radiation methods. (a), (b), and (c) describe the single-port feeding mode, the two-port feeding mode, and the four-port feeding mode, respectively.
Each point of the midsection temperature is evenly selected, which is used to figure out temperature coefficient of variation. Figure 3 shows the relation between temperature coefficient of variation, energy efficiency, and feeding mode. Figure 3 shows that energy efficiency of the single-port feeding mode is the highest, which is because the microwave feeding direction is oblique; microwave reflected several times in the microwave cavity; microwave can be fully absorbed by wood. But the microwave energy is relatively divergent, causing the uneven temperature distribution. Therefore, temperature coefficient of variation is larger. In the two-port feeding mode, temperature coefficient of variation is the largest, while energy efficiency dropped below 50%, which is not conducive to wood microwave pretreatment. In the four-port feeding mode, temperature coefficient of variation improved on the basis of the single-port feeding mode, reduced to nearly 20 percent. The temperature uniformity has been greatly improved, while energy efficiency is also high. Considering Figures 3 and 4 synthetically, high temperature region is distributed dispersedly and uniformly in the four-port feeding mode. And temperature coefficient of variation is small, which is in favor of rapid and uniform drying of wood. Thus, four-port feeding mode is the best choice in the three feeding modes. It also shows that the temperature uniformity do not improved as the feeding ports increasing, and the energy efficiency is not related to the feeding ports. The reason of temperature uniformity and energy efficiency is complicated, which is not only decided by the number of feeding ports.
The effect of microwave radiation methods on temperature coefficient of variation and energy efficiency.
The midsection temperature (°C) distribution within wood heated by microwave with the four-port feeding mode.
4.2. Effect of the Height of Cavities on Temperature Distribution and Energy Efficiency
Electromagnetic field distribution is influenced by the height of cavity (“h” for short), which has some influence on the heating effect. In order to determine the best height of cavity when thickness of the wood is 0.06 m, the respective temperature distribution of the wood is studied when the height of cavity is 0.07, 0.08, 0.09, 0.10, 0.11, 0.12, 0.13, and 0.14 m. Because the thickness of the actual wood always changes in the size, only the roughly height range of cavity needs to be determined.
Figure 4 shows the effect of the height of cavities on the midsection temperature distribution. High temperature region is gradually diverging from the center to the edge as the height of the cavity increases. When h≥0.11 m, the temperature of the edge area is gradually higher than in the central region, which is not conducive to wood microwave pretreatment. So the height of microwave cavity should be h<0.11 m. Figure 4 shows that when h≤0.08 m, high temperature region is too concentrated in the center, which might cause damage to the center of the wood. And the temperature distribution is uneven, so the height of cavity should be h>0.08 m.
Figure 5 shows the relation between temperature coefficient of variation, energy efficiency, and feeding mode as height of the cavity changes. Figure 5 shows that when h≤0.08 m, temperature coefficient of variation is very large; thus the temperature distribution is uneven. At the same time, energy efficiency is very low. As height of the cavity increases, energy efficiency which is determined by the shape of the cavity and the wavelength is also increasing. While considering the temperature distribution within the wood, h should not be too large. Using comprehensive analysis of temperature uniformity, drying efficiency, and energy efficiency, the ranges of the cavity height should be 0.08m<h<0.11m. It also shows that the energy efficiency improves with the cavity heightening, while the temperature uniformity remains almost unchanged. This is because when the volume of the cavity became larger, the distance of each feeding port became longer and electromagnetic leakage is less through feeding ports.
The effect of the height of cavities on temperature coefficient of variation and energy efficiency.
5. Conclusion
Effect of different feeding modes on wood microwave pretreatment is simulated and optimized parameters of cavity are achieved. Considering the wood temperature coefficient of variation and energy efficiency on different heating cavities, the following conclusions were made. Firstly, four-port feeding mode is the best among the three feeding modes on temperature uniformity and energy efficiency, and the height of cavities between 0.08 m and 0.11 m is the optimized parameter for wood whose thinness is 0.06 m. Temperature uniformity do not improve as the feeding ports increase, and energy efficiency is not related to the feeding ports. Lastly, energy efficiency improved with the heating cavity heightening, while temperature uniformity remains almost unchanged.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors acknowledge the support of the Special Fund for Forest Scientific Research in the Public Welfare (no. 201204708), National Natural Science Foundation of China (no. 31370564), New-Century Training Program Foundation for the Talents by the Ministry of Education of China (NCET-11-0979), and Graduate Student Innovation Fund of Central South University of Forestry and Technology (no. CX2014B20).
TurnerI. W.PuiggaliJ. R.JomaaW.A numerical investigation of combined microwave and convective drying of a hygroscopic porous material: a study based on pine wood199876219320910.1205/0263876985246222-s2.0-0031998685OloyedeA.GroombridgeP.Influence of microwave heating on the mechanical properties of wood20001001677310.1016/S0924-0136(99)00454-92-s2.0-0033884181LiX.ZhouY.YanY.CaiZ.FengF.A single cell model for pretreatment of wood by microwave explosion201064563363710.1515/HF.2010.0952-s2.0-77956192825YongDongZ.FengF.XianJunL.XiaoMeiJ.ZhiLinC.Effects of microwave treatment on residue growth stress and microstructure of Eucalyptus urophylla20093121461502-s2.0-75949084232TorgovnikovG.VindenP.Microwave wood modification technology and its applications201060217318210.13073/0015-7473-60.2.1732-s2.0-77956953533TorgovnikovG.VindenP.High-intensity microwave wood modification for increasing permeability200959484922-s2.0-67650074518PrzewlokaS. R.HannJ. A.VindenP.Assessment of commercial low viscosity resins as binders in the wood composite material Vintorg200765320921410.1007/s00107-006-0153-52-s2.0-34248995955SugiyantoK.VindenP.TorgovnikovG.PrzewlokaS.Microwave surface modification of Pinus radiata peeler cores: technical and cost analyses201060434635210.13073/0015-7473-60.4.3462-s2.0-78650082721TorgovnikovG.VindenP.New 300 kW plant for microwave wood modificationProceedings of the International Microwave Power Institutes 40th Annual Symposium2006Boston, Mass, USA260263VindenP.TorgovnikovG.HannJ.Microwave modification of Radiata pine railway sleepers for preservative treatment201169227127910.1007/s00107-010-0428-82-s2.0-79956214148GuangmingY.YuanL.YunchuH.Research on dispersing and surface treatment of inorganic-nanoparticles used for wood/inorganic-nano-composite20083211421452-s2.0-41349122754YuanG.LiuY.WuY.HuY.WuZ.Study on structure characterization and combining mechanism of Chinese fir/nano SiO2 composite200936959622-s2.0-71049162998ObandaD. N.ShupeT. F.BarnesH. M.Reducing leaching of boron-based wood preservatives—a review of research200899157312732210.1016/j.biortech.2007.12.0772-s2.0-44449096539KartalS. N.YoshimuraT.ImamuraY.Modification of wood with Si compounds to limit boron leaching from treated wood and to increase termite and decay resistance200963218719010.1016/j.ibiod.2008.08.0062-s2.0-58549106951BaysalE.AltinokM.ColakM.Kiyoka OzakiS.TokerH.Fire resistance of Douglas fir (Pseudotsuga menzieesi) treated with borates and natural extractives20079851101110510.1016/j.biortech.2006.04.0232-s2.0-33750302474PereyraA. M.GiudiceC. A.Flame-retardant impregnants for woods based on alkaline silicates200944449750310.1016/j.firesaf.2008.10.0042-s2.0-63449139411KartalS. N.HwangW.-J.YamamotoA.TanakaM.MatsumuraK.ImamuraY.Wood modification with a commercial silicon emulsion: effects on boron release and decay and termite resistance200760318919610.1016/j.ibiod.2007.03.0022-s2.0-34548862048CiudadA.LacastaA. M.HaurieL.FormosaJ.ChimenosJ. M.Improvement of passive fire protection in a gypsum panel by adding inorganic fillers: Experiment and theory20113117-183971397810.1016/j.applthermaleng.2011.07.0482-s2.0-84860407630KiliarisP.PapaspyridesC. D.Polymer/layered silicate (clay) nanocomposites: an overview of flame retardancy201035790295810.1016/j.progpolymsci.2010.03.0012-s2.0-77953608991AyappaK. G.Modelling transport processes during microwave heating: a review19971321692-s2.0-0002709805CampañoneL. A.PaolaC. A.MascheroniR. H.Modeling and simulation of microwave heating of foods under different process schedules20125273874910.1007/s11947-010-0378-52-s2.0-84856585625MarraF.De BonisM. V.RuoccoG.Combined microwaves and convection heating: a conjugate approach2010971313910.1016/j.jfoodeng.2009.09.0122-s2.0-71749115681