Technique for the electric and magnetic parameter measurement of powdered materials

Similar documents
CHARACTERIZATION OF MAGNETICALLY LOADED MICROWAVE ABSORBERS

Study of Specific Absorption Rate (SAR) in the human head by metamaterial attachment

Design of Composite Electromagnetic Wave Absorber Made of Fine Aluminum Particles Dispersed in Polystyrene Resin by Controlling Permeability

ELECTROMAGNETIC ENVIRONMENT GENERATED IN A TEM CELL FOR BIOLOGICAL DOSIMETRY APPLICATIONS

Electromagnetic Waves

Measurements of Electromagnetic Properties of Ferrites above Curie Temperature and impact on kicker performance

Transmission-Reflection Method to Estimate Permittivity of Polymer

Measuring Complex Permittivity, Permeability, or Sheet Impedance of Materials at Microwave Frequencies with a Free-Space Focused Beam

Spectral Domain Analysis of Open Planar Transmission Lines

USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION

<CMC Beads as a New EM Absorber >

An Approximate Approach to Determining the Permittivity and Permeability near λ/2 Resonances in Transmission/Reflection Measurements

PHY3128 / PHYM203 (Electronics / Instrumentation) Transmission Lines

UNIT I ELECTROSTATIC FIELDS

MEASURE THE COMPLEX PERMEABILITY OF FER- ROMAGNETIC THIN FILMS: COMPARISON SHORTED MICROSTRIP METHOD WITH MICROSTRIP TRANS- MISSION METHOD

Measuring Complex Permittivity, Permeability, or Sheet Impedance of Materials at Microwave Frequencies with a Free-Space Spot Probe System

A Modified Formula for Microwave Measurement of Dielectric Loss Using a Closed Cylindrical Cavity Dielectric Resonator

A New Method for VHF/UHF Characterization of Anisotropic Dielectric Materials

A Polynomial Approximation for the Prediction of Reflected Energy from Pyramidal RF Absorbers

A Microwave Method for Complex Permittivity Extraction of Thin Materials

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

ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER

A Broadband Flexible Metamaterial Absorber Based on Double Resonance

SCITECH Volume 4, Issue 1 RESEARCH ORGANISATION November 09, 2017

Research Article. Received 1 February 2016; Accepted 10 April 2016

Characterising Properties and Loss in High Powered Metamaterials

Lecture 36 Date:

Progress In Electromagnetics Research M, Vol. 20, 73 80, 2011

NASA Contractor Report. Application of FEM to Estimate Complex Permittivity of Dielectric Material at Microwave Frequency Using Waveguide Measurements

Exam in TFY4240 Electromagnetic Theory Wednesday Dec 9, :00 13:00

CONTROL OF MICROWAVE HEATING IN RECTANGULAR WAVEGUIDE

The Co-Conical Field Generation System A 40 GHz Antenna Test Cell

From Active Metamaterials to Transformation Electromagnetics: AMULET from the academic's perspective

A CALIBRATION PROCEDURE FOR TWO-PORT VNA WITH THREE MEASUREMENT CHANNELS BASED ON T-MATRIX

APPLICATION OF NON-BIASED FERRITES

INTRODUCTION TO TRANSMISSION LINES DR. FARID FARAHMAND FALL 2012

CHAPTER 9 ELECTROMAGNETIC WAVES

IPC-TM-650 TEST METHODS MANUAL

Electromagnetic Theory for Microwaves and Optoelectronics

Non-reciprocal passive protection of HF transmitters from reflection failure by the help of semiconductor isolators

Computation of Electromagnetic Energy Absorption in the Human Body Tissues by High Frequency Structure Simulator

A SIMPLE APPROACH FOR EVALUATING THE RECI- PROCITY OF MATERIALS WITHOUT USING ANY CAL- IBRATION STANDARD. U. C. Hasar

FDTD implementations for electromagnetic shields

Electromagnetic-Thermal Analysis Study Based on HFSS-ANSYS Link

Advances in Radio Science

1 The formation and analysis of optical waveguides

Lecture 12. Microwave Networks and Scattering Parameters

Contents. Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation. Measurements

ECE357H1S ELECTROMAGNETIC FIELDS TERM TEST March 2016, 18:00 19:00. Examiner: Prof. Sean V. Hum

Power Absorption of Near Field of Elementary Radiators in Proximity of a Composite Layer

A MATLAB GUI FOR SIMULATING THE PROPAGATION OF THE ELECTROMAGNETIC FIELD IN A 2-D INFINITE SPACE

SCATTERING CROSS SECTION OF A META-SPHERE

Current EMC Research at NIST

Chap. 1 Fundamental Concepts

Measurement of Complex Permittivity of Large Concrete Samples with an Open-Ended Coaxial Line

Engineering Electromagnetics

Kimmo Silvonen, Transmission lines, ver

New Aspects of Old Equations: Metamaterials and Beyond (Part 2) 신종화 KAIST 물리학과

ECE 107: Electromagnetism

5 RCS Management of Edge Diffracted Waves

Microstrip Coupler with Complementary Split-Ring Resonator (CSRR)

ELECTRODYNAMICS OF CONTINUOUS MEDIA

ELECTRICAL CALCULATIONS SOFTWARE. R.Ubic 1996

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

COLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE

Plane Waves GATE Problems (Part I)

Four Ports. 1 ENGN4545/ENGN6545: Radiofrequency Engineering L#16

Contents. Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation. Measurements

DESIGN OF MULTILAYER MICROWAVE BROADBAND ABSORBERS USING CENTRAL FORCE OPTIMIZATION

Grade 7 Revision sheet Key Answer. Kingdom of Bahrain Ministry of Education Ahlia School -ABCD

DETERMINATION OF THE COMPLEX PERMITTIVITY VALUES OF PLANAR DIELECTRIC SUBSTRATES BY MEANS OF A MULTIFREQUENCY PSO-BASED TECH- NIQUE

Eigenvalue Analysis of Waveguides and Planar Transmission Lines Loaded with Full Tensor Anisotropic Materials

Coupling Impedance of Ferrite Devices Description of Simulation Approach

4.4 Microstrip dipole

Electromagnetic materials characterization of ferromagnetic composites at microwave frequencies

On-Wafer Characterization of Electromagnetic Properties of Thin-Film RF Materials

This section reviews the basic theory of accuracy enhancement for one-port networks.

ELECTROMAGNETISM. Second Edition. I. S. Grant W. R. Phillips. John Wiley & Sons. Department of Physics University of Manchester

Finite Element Method (FEM)

A Methodology of Metamaterial Effective Permittivity and Permeability Value Measurement

Microwave Phase Shift Using Ferrite Filled Waveguide Below Cutoff

MATERIAL STUDY IN GHZ FREQUENCY DOMAIN OF A DIVALENT LIQUID CRYSTAL

FINAL EXAM IN FYS-3007

High Directivity Horn Antenna of Metamaterial in Terahertz Xiangjin Quan, Shiquan Zhang, Hui Li

Microstrip Antennas. Prof. Girish Kumar Electrical Engineering Department, IIT Bombay. (022)

ELECTROMAGNETIC INTERFERENCE (EMI) ANALYSIS FOR OBLIQUE INCIDENCE OF EM WAVES IN DOUBLE SHIELDS

Some Remarks on Shielding. Herbert Kapitza (FLA) (using slides from a talk by Mike Thuot) DESY,

Omar M. Ramahi University of Waterloo Waterloo, Ontario, Canada

Shielding Effectiveness of Modern Energy Saving Windows

Graduate Diploma in Engineering Circuits and waves

Measurement of Dielectric Substrate Parameters Using Split-Post Dielectric Resonator Assisted by 3D Simulators

General Appendix A Transmission Line Resonance due to Reflections (1-D Cavity Resonances)

TM-Radiation From an Obliquely Flanged Parallel-Plate Waveguide

E. Verney DIOM Laboratory Saint-Etienne University 23 rue Paul Michelon, St. Etienne, Cedex 42023, France

Cambridge International Examinations Cambridge International Advanced Level

Electromagnetic Modelling Process to Improve Cabling of Power Electronic Structures

Left-handed and right-handed metamaterials composed of split ring resonators and strip wires


ECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM. 28 April Examiner: Prof. Sean V. Hum. Duration: hours

Analysis and Design of the CRLH SICL Unit Cell using Effective Parameters

Transcription:

Computational Methods and Experimental Measurements XIV 41 Technique for the electric and magnetic parameter measurement of powdered materials R. Kubacki,. Nowosielski & R. Przesmycki Faculty of Electronics, Military University of Technology, Poland Abstract A measurement technique for the electric and magnetic properties of powdered ferrites has been developed. This technique allows one to determine relative permittivity and permeability of powdered materials. Measurements were taken in a coaxial transmission line, which guarantees the broadband frequency measurements. The calculations were based on evaluation of the scattering parameters of measured powdered material and two plastic walls. The values of relative permittivity (ε, ε ) and permeability (µ, µ ) have been done for ferrite powder in the frequency range from 00 MHz to 100 MHz. Keywords: absorbing materials, permittivity and permeability measurements. 1 Introduction With the increasing quantity of devices emitting electromagnetic radiation, there is a need to develop and introduce into practice materials and techniques to protect against unwanted radiation. This is important for EMC purposes as well as for the protection of people against the harmful radiation. There are two ways in which materials can shield against radiation: by reflecting or by absorbing the incident electromagnetic energy. The ideal absorber should have low reflectivity and a high value of absorption of the incident electromagnetic energy. The electromagnetic field incident to the boundary surface of any material is reflected and the level of reflected energy is a function of the internal parameters of the material. In general, the shielding effectiveness of a material or configuration of materials is a measure of its ability to attenuate electromagnetic energy. This ability depends on both its reflection and absorption properties. Energy not being reflected or absorbed by the material is transmitted from one side to the other. doi:10.495/cmem0901

4 Computational Methods and Experimental Measurements XIV For a given material, the amount of energy transmitted has a complex dependence upon the angle of incidence and polarization of the incoming electromagnetic wave. There are a few ways to obtain the minimal value of reflected power; however, the most often applied technique is to cover the protected walls with graphite materials having a pyramidal layer structure. Unfortunately, such structures operate correctly at rather higher (microwave) frequencies. At frequencies lower than 10 MHz, such a technique with a pyramidal layer structure does not guarantee the decrease of the level of reflected energy because, generally speaking, the wavelength is much bigger than the highest of the pyramids. Another way to obtain low reflectivity of material is to use materials having not only electric properties, but also magnetic properties. Materials with such properties are ferrite materials. Nowadays, walls in anechoic chambers are covered by broadband absorbers composed of two layers. The first one is pyramidal structure polyurethane foam loaded with graphite and the second layer is a sintered ferrite. Commercially available monolithic ferrite materials guarantee good absorption properties and a low level of reflected field; however, in many practical uses such a solid form of absorber cannot be introduced. From this point of view there is a need to develop absorbers that can fulfill the following requirements: - operate at desirable frequencies or over a wide frequency band; - be a solid or flexible form of material; - have the desirable level of reflectivity and absorption. New materials can be used for specific studies, e.g. stealthness, but in this case rather composite materials should be developed. It is possible to simulate the properties of an absorber having the above properties; however the constitutive material parameters, e.g. relative permittivity ε r and permeability µ r, of such materials or any component of these materials should be determined. In many cases components are in powder form. In fact, most of absorbing material is usually composed of powder components joined together with an epoxy bonder. The important thing is to assess the electric and magnetic parameters of the components of the mixture, being in powder form before their mixing. Network analyzer measurements There are many methods that can be used to determine the permittivity ε r and the permeability µ r of materials. Resonant cavity techniques offer high accuracy of measurements but they can only be used mainly for one frequency. On the other hand, waveguide methods make possible measurements over a wide frequency band but they have some disadvantages when measured materials are lossy. Nevertheless, measurements in waveguides are based on measuring of the scattering parameters (S ik ). In this paper measurements of scattering parameters were taken according to the approach given in Baker-Jarvis [1]. This approach refers to toroidal material

Computational Methods and Experimental Measurements XIV 43 s material under test 1 m Figure 1: Coaxial transmission line with a sample of material under test. under test inserted into a transmission coaxial line fig. 1. In this method the specimen is treated as homogeneous. Baker-Jarvis [1] described the technique of determining values of the propagation constant (γ m ) and characteristic impedance (Z m ), taking into account the measured scattering parameters. However, in order to determine the values of γ m and Z m of the material under test (m), the scattering parameters of this sample should be extracted from the scattering parameters of whole coaxial section (s). In this case airline sections having lengths 1 and must be compensated. This is an important task to minimize measurement errors due to uncontrolled displacement of the sample. To minimize these errors a method was developed taking into account the nonlinear optimization of the calculation with inserted airlines sections [1]. Measurements of the scattering parameters of the sample inserted into a coaxial line were taken using a vector network analyzer system. However, special additional accessories are necessary for the calibration process. Effective auxiliary tools for calibration were developed and described in [, 3]. These tools are in fact three airlines with different lengths depending on the frequency range measurements. This special calibration process guarantees the highest accuracy of measurements. Taking into account calibration evaluation, the computational program can compensate for errors caused by any potential displacement of the sample. The above-mentioned measurement technique is also called the reflection/transmission method. This method is well adapted to broadband characterization of isotropic and homogeneous materials; however, there are some disadvantages to this technique. On one hand, at frequencies where ε µ n λ m r r = (n integer) S 11 vanishes and in such case noise affects measurements. Consequently, the network analyzer phase uncertainty is large. On the other hand, when material is lossy (or thickness of a sample is large) the value of S 1 can also vanish. So, the reflection/transmission method should be applied to the characterization of several materials, where permittivity and permeability are approximately known.

44 Computational Methods and Experimental Measurements XIV 3 Measurement method of the relative permittivity and permeability of powdered materials There are few waveguide methods that allow one to measure the complex relative permittivity ε r of non-solid materials. In the case of liquid materials, measurement methods are easy to realize, especially when they are lossy. For example, a typical waveguide setup used for measuring liquid biological tissues consists of a coaxial open-ended line applied to a liquid surface (Kubacki et al. [4]). In this case the measurement of the complex reflection coefficient is sufficient to determine the values of two parameters ε r (ε, ε ) of the liquid. However, when the material also has magnetic properties µ r (µ, µ ), such measurements based on determining only the reflection coefficient are not sufficient. In this case the transmission coefficient should also be measured. This could be done using waveguides or coaxial lines where values of two complex scattering parameters can be measured. The method described in this paper allows one to measure complex relative permittivity ε r and permeability µ r of non solid (mainly powder) materials over a wide frequency range from a few khz to a few GHz. The technique of measurement has been based on a method developed for coaxial transmission lines where the material under test must have toroidal form. In the proposed method, the material under test is inserted into a coaxial line between two plastic walls as is shown in fig.. This allows one to measure powdered materials. The values of electric (ε, ε ) and magnetic parameters (µ, µ ) of measured powder are extracted in a mathematical evaluation from the measured data of scattering parameters of the whole sample (S 11, S, S 1, S 1 ). Mathematical evaluation that allows one to extract the electromagnetic data of the measured powder from measured data of the whole sample has been derived on the basis of scattering parameters describing two-port networks as is presented in fig.3. In this figure, the scattering parameters of plastic walls were defined as E ik while the scattering parameters of measured powdered material are defined as S ik. measured powder plastic walls w w Figure : The sketch of the coaxial line with a sample of material under test between two plastic walls.

Computational Methods and Experimental Measurements XIV 45 E1 A S1 E1 E11 Eik E S11 Sik S E11 Eik E E1 S1 E1 A Figure 3: Scattering parameters of plastic walls (E ik ) and powder under test (S ik ). Taking into account the reciprocal and source-free networks, the scattering matrix satisfies the following condition: S ik =S ki. In our case for symmetrical twoport networks: S 11 =S, S 1 =S 1 and also E 11 =E, E 1 =E 1. For each network it is possible to formulate the following scattering equations: b1 S11 S1 a1 = (1) b S1 S a Scattering matrix (1) refers to a single network and for the sample composed of three layers such matrix becomes much more complicated. To receive mathematical formulas that allow one to extract scattering parameters of networks representing the powder substrate, the so called graph method of calculation has been introduced. The graph method has been derived from the theory of power flow in the network branches. As an example of how to use graph method in practice the reflection coefficient at the gate A-A (fig. 3) has been presented in fig. 4. A S1 ΓA S11 Sik S E11 Eik S1 A Figure 4: The graph method used for the example analysis of the reflection coefficient at the gate A-A. The final formulas of the scattering parameters Sw ik of the whole sample (material under test and two plastic walls) are as follows: Sw 11 11 1 ( 1 E11S11) S11 + E11S ( 1 E S ) E S = E + E () 11 11 11 1 1

46 Computational Methods and Experimental Measurements XIV Sw E1S1 1 = (3) ( 1 E11S11) E11S1 where Sw ik represents the scattering parameters of whole sample. Formulas () and (3) allow one to determine the scattering parameters of the specimen (S 11, S 1 ) from the scattering data of whole sample (Sw 11, Sw 1 ). 4 Experiments and results Measurements of complex relative permittivity and permeability were carried out using a vector network analyzer (VNA). This system consisted of a 7 mm coaxial airline, equipped with measurement cables and PC7 connectors. The center conductor of the coaxial airline is 3,04 mm in diameter to receive the 50 Ω characteristic impedance of the holder. The system measures the magnitude and phase of four scattering parameters (S 11, S, S 1, S 1 ). First of all the empty coaxial airlines were employed to calibrate the system. The additional tools consisted of three auxiliary airlines that were used to realize the calibration. These airlines are sample holders as well as calibration tools. They have been presented in fig. 5. The whole measurement setup is presented in fig. 6. Taking into account that various components of the VNA introduce magnitude and phase uncertainties, the calibration process is important in order to obtain the highest accuracy of measurements [1]. Calibration removes the systematic measurement uncertainties of the system. The calibration coefficients were determined by solving a set of simultaneous equations generated from the linear fractional transformation. After calibration, when the system was operated with the errors correction, the measurements were updated by the calibration coefficients. Data of calibration process and measurements of the sample were acquired using the MultiCal program, than they were processed with a special program developed for extracting values of permittivity and permeability. Figure 5: Coaxial airlines as sample holders and/or calibration tools.

Computational Methods and Experimental Measurements XIV 47 Figure 6: S-parameter test of the vector network analyzer and a coaxial line with sample. E 0 [ (1-ρ )e -γ ] E 0 ρ E 0 [ -ρ(1-ρ )e -γ ] E 0 [ ρ (1-ρ )e -3γ ] E 0 [ -ρ 3 (1-ρ )e -4γ ] E 0 [ ρ 4 (1-ρ )e -5γ ] E 0 Figure 7: The infinite reflections and transmissions according to the theory of optical rays. The final formulas of the scattering parameters of a powder (S 11, S 1 ) were extracted from the scattering parameters of the whole sample () and (3). On the other hand, these scattering parameters could be determined by analyzing the electric field at the specimen interface. In order to determine the material properties from scattering data there is a need to introduce the relationship of the electric field incident at the boundary of the specimen and propagate through the materials. The reflection and transmission coefficients can be interpreted as the infinite summation of optical rays when the incident electromagnetic field is treated as a plane wave. Such a situation is typical in a coaxial line where the TEM mode is propagated. The infinite reflections and transmissions are shown in fig. 7.

48 Computational Methods and Experimental Measurements XIV The reflection and transmission of rays have been presented according to the theory of optics (fig. 7), taking into account two boundary conditions: air/sample as well as sample/air. The reflection coefficient (S 11 ) consists of infinite rays in which the first one is reflected at the first interface air/sample (ρ) and all the others are transmitted, arising from multiple reflections inside the sample. On the other hand, the transmission coefficient (S 1 ) is the infinite summation of the rays transmitted to the second interface sample/air. This infinite sum of components is in fact a geometric series yielding to the following final sums: γ 1 e S11 = ρ (4) γ 1 ρ e γ ( 1 ρ ) e S1 = (5) γ 1 ρ e where, in the coaxial transmission line with the fundamental mode TEM: µ r ε r ρ - reflection coefficient: ρ = µ + ε γ - propagation constant: γ = j ω ε rµ r Formulas (4) and (5) are functions of unknown parameters of ε r and µ r of the specimen. Values of these parameters (ε, ε, µ µ ) can be obtained by solving the above complex equations (4) and (5). In this case the iterative procedure yields to a very stable solution for a specimen of arbitrary length. For the purpose of the validation of the solution of the sample of solid material, the measurement in the configuration with two plastic walls was used. Then, the values of ε, ε, µ µ were determined using the proposed method. The obtained data were confirmed by measuring this solid sample (without plastic walls) using the reflection/transmission method described in [1, 3]. The measured values of ε, ε, µ and µ by the proposed method were within good accuracy with the reflection/transmission method. As an example of powdered material measurements, the values of complex permittivity and permeability of ferrite powder have been measured using the proposed method. Ferrite powder was received in a ball mill from the solid state material (spinel class ferrite, named G-175) to receive granulate (0,-) µm. For this powdered material, relative values of permittivity (ε, ε ) versus frequency were presented in fig. 8, while relative values of permeability (µ, µ ) were presented in fig. 9. The measured values of permittivity (ε, ε ) and permeability (µ, µ ) of the ferrite powder versus frequency were presented in figs. 8 and 9. However, the obtained data for the powdered ferrite were significantly lower than that of solid state ferrite. This is because the magnetic domains of ferrite were destroyed during the milling process. Further investigation should be undertaken to characterize the relationship between the volume of granulates and the decreasing permeability. r r

Computational Methods and Experimental Measurements XIV 49 ε ε 6 5 4 3 1 0 00 300 400 500 600 700 800 900 1000 1100 100 f [MHz] Figure 8: Values of the relative permittivity of ferromagnetic powder in the function of frequency. µ µ 3 1 0 00 300 400 500 600 700 800 900 1000 1100 100 f [MHz] Figure 9: Values of the relative permeability of ferromagnetic powder in the function of frequency. In the described measurements all grains of powder were relatively small, compared to the wavelength, to receive homogeneous material. Measurement problems can appear in the measurement transmission line when material under test is heterogeneous. This can happen when materials contain inclusions or when grains of powder are not small compared to the wavelength. In these conditions the higher order modes can be excited in the sample by diffraction on the heterogenities. So, the parameters ε eff and µ eff of such a sample cannot be deduced in classical calculation because it supposes that only the dominant mode propagates. Moreover, in such structures with local inclusions, the local resonances are excited. This causes a lack of energy to be measured when higher order modes appear. This is due to metallic losses in the coaxial line and inclusions.

50 Computational Methods and Experimental Measurements XIV 5 Conclusions The measurement technique of electric and magnetic properties of powdered ferrites has been developed. This technique seems to be useful for lossy materials. Based on this technique, values of relative permittivity (ε, ε ) and permeability (µ, µ ) of ferromagnetic powders can be measured in the frequency range from 00 MHz to 100 MHz. The aim of this work was to investigate the measurement technique of powdered materials, because it could be an important tool for simulating the final properties of absorbing materials composed of conductor and ferrite materials. However, it turned out that in the case of ferrite powder such an attempt falls down because the values of relative permeability (µ, µ ) of powder was significantly lower than that of solid state ferrite. Further investigation should be directed to find the difference in permeability between the powdered and solid state ferrites. Acknowledgement The authors would like to thank Ryszard Frender from the Telecommunications Research Institute, Warsaw, Poland for providing ferrite powders for the measurements. References [1] Baker-Jarvis J., Transmission/reflection and short-circuit line permittivity measurements, NIST Technical Note, 005. [] Marks R.B., A multiline method of network analyzer calibration, IEEE Trans. Microwave Theory Tech., 4, pp. 105-115, 1991. [3] Wiatr W., Frender R., Żebrowski M., Characterization of microwave absorbing materials using a wideband T/R measurement technique, Conference Proceedings on International Conference on Microwaves, Radar and Wireless Communications, Wrocław, Poland,, pp. 475-478, 008. [4] Kubacki R., Sobiech J., Wardak K., The comparison of dielectric properties of young and mature animal tissues in microwaves, Przegląd Elektrotechniczny, 5, pp. 43-45, 006.