ESTABLISH THE DIELECTRIC PROPERTIES OF A NUMBER OF WOOD SPECIES FOR A RANGE OF FREQUENCIES, TEMPERATURES AND PRESSURES.

Transcription

1 ESTABLISH THE DIELECTRIC PROPERTIES OF A NUMBER OF WOOD SPECIES FOR A RANGE OF FREQUENCIES, TEMPERATURES AND PRESSURES Georgiana Daian Doctor of Philosophy 2005 SWINBURNE UNIVERSITY OF TECHNOLOGY

2 Abstract The present research offers an assessment of a new way of approaching the dielectric properties of wood. The concept takes into consideration a probable influential factor which has never been considered and refers to the elevated pressure which develops within wood during microwave processing. In addition, in order to understand and analyse the wood dielectric properties dependency on the internal temperature and pressure, a custom-made device and the appropriate measuring procedure, as well as a numerical model which calculates complex dielectric permittivity of an anisotropic wood structure is presented. Initially, the fundamentals of microwave wood interaction including polarisation mechanism, factors affecting the dielectric response of wood, microwave related characteristics of wood and its microwave heating process are examined. The importance of measuring the wood dielectric properties is also posed. Subsequently, an insight into the main aspects of the techniques used in other industries to measure the materials dielectric properties is provided. The report ends with a comparative discussion of the measurement techniques applicability. A method for measuring the dielectric permittivity of wood is then presented together with a description of the software required to extract the dielectric parameters from the measured quantities. The experimental results of the dielectric parameters of the main Australian wood species are also presented. Besides, a numerical model built to describe the 3-D wood structure is put forward. This model was introduced into an efficient solver that calculates the effective dielectric constant of any three dimensional structure of dielectric materials. The analysis of the calculated results for the permittivity and the comparison to the measured values shows a very practical qualitative and quantitative agreement between all the measured and calculated values. Finally, the importance of the research, the key results and the recommendations for further work are emphasized. Keywords: dielectric properties of wood, microwave processing, 3-D wood model, wood internal temperatures and pressures, measurements techniques. i

3 Acknowledgement I would like to express my gratitude to my academic supervisor, Dr. Alex Taube of the Swinburne University of Technology, for his kind assistance and invaluable guidance throughout this study. I am grateful to Professor Gregory Torgovnikov of Melbourne University for serving on my Ph.D. committee and for his support and valuable feedback. I would like to extend my gratitude to the Cooperative Research Center for Wood Innovations, especially to Professor Peter Vinden for providing motivation behind this research program and the funding to make it possible. Additionally, I thank to the Industrial Research Institute Swinburne for all the facilities that were made available to me in order to undertake this research. I am deeply grateful to Dr Amikam Birnboim of University of Haifa for his effective collaboration and valuable guidance in performing the computer-aided simulation. Sincere appreciation is also extended to Dr. Yury Shramkov (Industrial Research Institute Swinburne) for providing considerable technical support. Special thanks go to Professor Elias Siores for having confidence in my research abilities and offering encouragement on numerous occasions. I would like to thank my parents for their continuous moral support and for their simply being there despite the distance which separated us, and I would like to thank my husband for his unfailing support. ii

4 DECLARATION This thesis contains no material which has been accepted for the award of any other degree or diploma, except where due reference is made in the text of the thesis. To the best of my knowledge contains no material previously published or written by another person except where due reference is made in the text of the thesis. Signed: Georgiana Daian Dated: August 2005 iii

5 Table of Contents Abstract...i Acknowledgement...ii Table of Contents...iv List of Figures...vii List of Tables...xi Nomenclature...xiii Abbreviations...xv Introduction...1 PART I LITERATURE REVIEW CHAPTER 1: MICROWAVE WOOD INTERACTION FUNDAMENTALS Microwave History, Applicability and Definition Maxwell s Equations and Electromagnetic Waves Polarisation of Wood Time and frequency dependent polarisation Wood dielectric response as the result of dipole polarisation Factors Affecting the Dielectric Properties of Wood Structural Components of Wood Wood Moisture Content Wood Anisotropy Wood Density Wood Internal Temperature Microwave Frequency Wood Characteristics Related to Microwave Heating Heterogenous and anisotropic Poor thermal, moderate electrical conductor and non permeable under microwave energy Microwave Heating Process of Wood The Importance of Measuring the Dielectric Properties of Wood...27 CHAPTER 2: TECHNIQUES FOR MEASURING THE DIELECTRIC PROPERTIES OF MATERIALS AT MICROWAVE FREQUENCIES Classifications of the Dielectric Properties Measurement Techniques Transmission Line Techniques Transmission lines fundamentals and electromagnetic field distribution Operating frequency domain; measurements errors and results accuracy; advantages and disadvantages in using the different transmission line techniques Understanding the Fundamental Principles of the Transmission Line Measurements Measurement system set-up...43 iv

6 Reflection and transmission parameters Transmission lines terminations. Fundamentals of electromagnetic waves for each particular case Algorithms for calculating the dielectric parameters Transmission-Reflection method algorithms Reflection methods algorithms Transmission line calibrations Free Space Technique Operating frequency domain; measurements errors and results accuracy; advantages and disadvantages in using the technique Measurement system set-up Free-space calibration Algorithms for calculating the dielectric parameters Open-End Coaxial Probe Technique Operating frequency domain; measurements errors and results accuracy; advantages and disadvantages in using the technique Coaxial probe calibration Measurement system set-up Models for calculating the dielectric parameters Resonant Cavity Techniques Resonant structures fundamentals and electromagnetic field distribution Operating frequency domain; measurements errors and results accuracy; advantages and disadvantages in using resonant cavities techniques Measurement system set-up Cavity perturbation method: the principle for calculating the dielectric parameters Comparison of the Permittivity Measurement Techniques and the Review Outcome...77 PART II RESEARCH METHODOLOGY, DEVELOPED TECHNIQUES, RESULTS AND ANALYSIS CHAPTER 3: RESEARCH METHODOLOGY Defining the Research Question and Objectives Approach to the Research Problem...83 CHAPTER 4: MEASURING DIELECTRIC PROPERTIES OF WOOD AT MICROWAVE FREQUENCIES, ELEVATED TEMPERATURES AND PRESSURES Measurement Technique and Experimental System Set-Up Measurement Device Calibration Experimental Method Numerical Procedure Verification of Experimental Method and Numerical Procedure...95 v

7 CHAPTER 5: MODELING DIELECTRIC PROPERTIES OF WOOD AT MICROWAVE FREQUENCIES Structural Design of the Wood Modelling Software Method for Calculating the Permittivity and the 3-D Wood Model Software s Input Instructions To create wood model geometry To attribute different densities and moistures To calculate the structure s permittivity Dielectric parameters of air Dielectric parameters of free water Dielectric parameters of wood substance/moist wood substance Dielectric parameters of rays CHAPTER 6: RESULTS AND DISCUSSION Complex Permittivity of Australian Wood Species Experimental results and analysis Modeling results and analysis Discussion CHAPTER 7: CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK References APPENDIX A: Additional Data Related to the Fundamental Principles of the Transmission Lines and Resonant Cavities APPENDIX B: The Input File Interface of the Software which Extracts the Dielectric Parameters from the Measured Quantities Via the Solution of a Transcedental Equationin in the Complex Plane APPENDIX C: The Input File Interface of the Software which Predicts the Dielectric Properties of Wood APPENDIX D: Additional Data Related to Dielectric Parameters of Free Water, Bound Water and Cell Wall Substance in Wood List of Publications vi

8 List of Figures Figure 1: The potential of microwave applications in wood industry...6 Figure 2: Propagation of a plane wave...8 Figure 3: Electromagnetic spectrum and frequencies used in microwave processing...9 Figure 4: Frequency dependence of the various polarisation mechanisms...12 Figure 5: Permittivity of the CWS at 10 9 GHz and high temperatures...15 Figure 6: Dielectric parameters of water as a function of temperature, at f = 2.45 GHz...17 Figure 7: Typical structure of hardwood...18 Figure 8: Dielectric parameters of wood as a function of density and moisture content, at f = 2.45 GHz and T = 20 C...20 Figure 9: Dielectric parameters of wood (ρ 0 = 0.7 g/cm 3 ) as a function of moisture content and temperature, at f = 2.45 GHz...21 Figure 10: Microwave frequency response of the wood dielectric parameters at 20 C.22 Figure 11: Principal directions on wood stem...23 Figure 12: Physical model of microwave heating of wood...25 Figure 13: Microwave drying periods for wood...26 Figure 14: The scheme for scaling-up a microwave process...28 Figure 15: Exploratory research description...28 Figure 16: Exponential decay of the electric field in a lossy dielectric...30 Figure 17: Measurement techniques versus frequency and dielectric loss of MUT...34 Figure 18: Equivalent circuit of a transmission line...35 Figure 19: Waveguide types...36 Figure 20: Coaxial transmission line...37 Figure 21: Field distributions in coaxial line...38 Figure 22: Parallel plate transmission line...38 Figure 23: Microstrip line - configuration and field distribution...39 Figure 24: Strip line configuration...40 Figure 25: Transmission line measurement system...43 vii

9 Figure 26: Definition of the two-port network...44 Figure 27: S-parameters...44 Figure 28: Measuring S-parameters...45 Figure 29: Familiar forms of the reflection/transmission parameters...45 Figure 30: Reflection parameters...46 Figure 31: Transmission parameters...46 Figure 32: Transmission line terminated with Z Figure 33: Transmission line terminated with Short...48 Figure 34: Transmission line terminated with Open...48 Figure 35: Standing waves...49 Figure 36: A dielectric sample in a transmission line (waveguide/coaxial line) and the incident and reflected electric field distributions...54 Figure 37: Short-circuited waveguide with two-layer medium...56 Figure 38: Error model for one-port measurements...59 Figure 39: Forward two-port error model...60 Figure 40: Schematic figure of the free-space measurement system...63 Figure 41: Coaxial probe technique; reflection from the MUT...68 Figure 42: Coaxial probe model...69 Figure 43: A parallel resonant circuit formed from conventional circuit elements...70 Figure 44: A rectangular cavity resonator and dielectric field distributions for the TE 101 and TE 102 resonance modes...71 Figure 45: Cylindrical cavity resonator...72 Figure 46: Microwave resonant cavities...74 Figure 47: Transmitted power response of unloaded resonator at resonant frequency...76 Figure 48: Microwave pressure windows availability and their applications...86 Figure 49: Designed high-pressure cavity...86 Figure 50: Measurement system set-up...87 Figure 51: Schematic representation of the experiment...89 Figure 52: Three dimensional plot of minus the imaginary part of tanh(z)/z...94 Figure 53: Structural design of the wood modeling software...99 viii

10 Figure 54: Algorithm principle: a three dimensional box divided into a mesh of small cells Figure 55: Typical structure of hardwood Figure 56: A cross section through the 3-D wood model Figure 57: A radial section through the 3-D wood model Figure 58: A tangential section through the 3-D wood model Figure 59: Specific gravity of bound water Figure 60: 3-D ray s cell model Figure 61: Measured dielectric parameters of Blue Gum as a function of moisture content and structural direction, at C Figure 62: Dielectric constant of different Australian wood species at 2.45GHz, elevated temperatures and pressures Figure 63: Loss factor of different Australian wood species at 2.45GHz, elevated temperatures and pressures Figure 64: Modeling vs. measurement. Dielectric constant of Mountain Ash, ρ 0 =0.75g/cm Figure 65: Modeling vs. measurement. Loss factor of Mountain Ash, ρ 0 =0.75g/cm Figure 66: Modeling vs. measurement. Dielectric constant of Mountain Ash, ρ 0 =0.60g/cm Figure 67: Modeling vs. measurement. Loss factor of Mountain Ash, ρ 0 =0.60g/cm Figure 68: Modeling vs. measurement. Dielectric parameters of Messmate Figure 69: Modeling vs. measurement. Dielectric constant of Blue Gum Figure 70: Modeling vs. measurement. Loss factor of: Blue Gum Figure 71: Estimated dielectric parameters of bound water in wood at 2.45 GHz and room temperature Figure 72: Calculated dielectric parameters of pure water (or free water in wood) at 2.45GHz, as function of temperature Figure 73: Dielectric parameters of cell wall substance for fiber at 2.45GHz, as function of temperature Figure 74: Estimated dielectric parameters of bound water in wood at 2.45 GHz, as a function of temperature ix

11 Figure A-1: Field distributions and key expressions of calculation for modes in rectangular waveguide (NMAB, CETS, and NRC, 1994) Figure A-2: Field distributions and key expressions of calculation for modes in cylindrical waveguide (NMAB, CETS, and NRC, 1994) Figure A-3: Instantaneous values of incident and reflected waves on an open-ended line (Integrated Publishing 1998) Figure A-4: Electromagnetic-field configuration of TE mnk mode in rectangular cavity resonator (Pozar 1998) Figure A-5: Electromagnetic field configuration of TE mnk mode in rectangular cavity resonator (Pozar 1998) Figure A-6: Field configuration of TE 011 and TE 111 mode in cylindrical cavity resonator (Pozar 1998) Figure A-7: Field configuration of TM 011 and TM 111 mode in cylindrical cavity resonator (Pozar 1998) x

12 List of Tables Table 1: Measurement techniques according to the operation in the broad/narrow band frequency...33 Table 2: Measurement techniques classified upon their destructive features...34 Table 3: Comparison of the permittivity measurement techniques...78 Table 4: Complex solutions of the transcendental equation together with the values of the dielectric constant derived from each solution...95 Table 5: Dielectric permittivity of Soda Lime Glass measured at various positions of the short...97 Table 6: Fractional volumes formulas Table 7: Calculation of the specific gravities of wood Table 8: Relaxation parameters of bound water Table 9: Oven-dry densities of wood pieces from a set of six Table 10: Measured dielectric permittivity of Blue Gum at 2.45GHz and C Table 11: Measured dielectric properties of different Australian wood species at 2.45GHz, elevated temperatures and pressures Table 12: Anatomical dimensions of the macroelements in Eucalyptus globulus Table 13: Proportional volume of wood components for the six Blue Gum samples..123 Table 14: Dielectric parameters of moist wood substance and ray material for the six Blue Gum samples Table 15: Calculated and measured dielectric permittivity for Blue Gum wood samples (data obtained and calculated at temperature 20 o -25 o C) Table 16: Calculated and measured dielectric permittivity for wood samples at elevated temperatures and pressures Table D-1: Free water relaxation parameters function of temperature Table D-2: Free water dielectric parameters function of frequency and temperature..176 Table D-3: Dielectric parameters of the cell wall substance for fibers as function of frequency and temperature Table D-4: Dielectric parameters of the cell wall substance for rays as function of frequency and temperature xi

13 Table D-5: Estimated dielectric parameters of bound water in wood at 2.45GHz Table D-6: Estimated dielectric parameters of bound water in wood at 0.915GHz xii

14 Nomenclature a a 1, a2 B b 1, b 2 C Radius of a dipole molecule Incident waves Magnetic field [T] Reflected waves Capacitance c 8 Speed of light in free space ( = m/s) c speed of electromagnetic waves propagating through a dielectric [m/s] D Electric flux density [C/m 2 ] D p Penetration depth [m] d E E f G 0 Sample length [m] Electric field strength [V/m] Voltage [V] Oscillation frequency of the fields at a given point in space [Hz] Resonant frequency [Hz]. Chapter 2, Eq 54 to 59 Conductance accounting for radiation losses H I Magnetic field intensity [A/m] Current [A] J Current density [A/m 2 ] 23 k Boltzmznn s constant ( = J/K) k Volumetric shrinkage coefficient at fiber saturation point FSP Q S 11, S22 S 12, S 21 SG, SG, SG 0 x Quality factor Reflection scattering parameters (S-parameters) Transmission scattering parameters (S-parameters) FSP SG bw, x, SG bw, FSP Specific gravity of wood at 0% moisture content, x% moisture content and fiber saturation point, respectively. Specific gravity of bound water at x% moisture content and fiber saturation point, respectively xiii

15 T t Temperature [ C] Transmission coefficient. Chapter 2 Time [s] tan δ Loss tangent tan δ L, δ R tan δ T tan, Loss tangent values in longitudinal, radial and tangential direction V Resonant cavity volume [m 3 ] V ws, V bw, V fw, a ir V Fractional volume of wood substance, bound water, free water and air v Sample volume [m 3 ] Y Z 0 Admittance Characteristic impedance[ω] Z in Z L Line impedance [Ω] Load impedance [Ω] α Relaxation time distribution constant. Chapter 1, Eq 13 Attenuation constant [Np/m]. Chapter 1, Eq 16, 17 and 18 δ Phase angle [rad] ε Permittivity of free space ( = C 2 /N m 2 ) 0 ε ε eff Complex permittivity Effective permittivity ε Dielectric constant of the dielectric medium ε, ε, ε Dielectric constant in longitudinal, radial and tangential direction L ε ε s R T Dielectric constant at very high frequency (optic) Dielectric constant at very low frequency (static) ε Dielectric loss factor Φ Γ γ 0 Phase of the transmission coefficient [rad] Reflection coefficient Propagation constant of the empty waveguide γ Propagation constant of the filled waveguide η Viscosity [N s/m 2] ϕ Phase of the reflection coefficient [ rad] xiv

16 λ Wavelength in free space ( = m) 0 λc Cutoff wavelength [m] λg Wavelength in waveguide [m] µ Effective permeability eff 7 µ Permeability of free space ( = 4π 10 Wb/A m) 0 ρ 3 Volume charge density [C/m ] Magnitude of the reflection coefficient. Chapter 2 ρ Oven-dry density of wood [g/cm 3 ] 0 σ e Electronic conduction σ i τ ω Ionic conduction Relaxation time [s] Magnitude of the transmission coefficient. Chapter 2 Angular frequency [rad/s] Abbreviations CWS Cell wall substance DUT Device under test FSP Fiber saturation point (%) mc Moisture content (%) MUT Material under test NRW Nicolson-Ross-Weir procedure RL Power loss TE Transverse electric waves TEM Transverse electromagnetic waves TM Transverse magnetic waves TR Transmission-reflection VNA Vector Network Analyser VSWR Voltage Standing Wave Ratio xv