Movable shortcircuit technique to extract the relative permittivity of materials from a coaxial cell
M. G. Lountala^{1} , F. Moukanda Mbango^{2} , F. Ndagijimana^{3} , D. LilongaBoyenga^{4}
^{1, 2, 4}Electrical and Electronics Engineering Laboratory, Marien Ngouabi University, BP. 69, Brazzaville, Congo
^{3}Grenoble Electrical Engineering Laboratory, University of Grenoble Alpes, CNRS, Grenoble INP, G2Elab, F38000 Grenoble, France
^{2}Corresponding author
Journal of Measurements in Engineering, Vol. 7, Issue 4, 2019, p. 183194.
https://doi.org/10.21595/jme.2019.20925
Received 26 July 2019; received in revised form 23 November 2019; accepted 16 December 2019; published 31 December 2019
JVE Conferences
In recent years, industrial applications have been based on the use of intrinsic material properties that improve designs, processes, qualities and product controls. To get to those intrinsic parameters, various appropriate techniques are required. In this paper, a new technique has been developed and presented. It essentially puts the emphasis on the dielectric relative permittivity extraction from the principle of a movable shortcircuit through the coaxial transmissionline cell. This technique is aimed at drastically reducing the discontinuity impacts at the interface feed line (connector) and ideal line, solving the phase constant frequency limit, stopping the constraints bound to the higher mode propagations and improving the accuracy level when the frequency range has increased. The technique is based on the use of the sum of two different lengths of the cell by removing the first value of the phase constant in the frequency range of interest when it is negative. This new technique can be easily implemented; its focus is not on iterative principles, but on the use of the constant propagation of a QuasiTEM mode of the transmissionline. The biofood industry (semolina), environmental field (palm tree) and building trade (aquarium sand) were used to test the validity of the technique in 220 GHz.
Keywords: Discontinuity impacts, movable shortcircuit, propagation constant, relative permittivity, transmissionline.
1. Introduction
The material characterization has a wide variety of techniques [13], depending on the application domain [4]. These domains help to reveal why different materials show different properties and behaviors [57]. A large panel of material characterization techniques is described in the literature. Most of accurate techniques are developed in relation to the fixture’s geometry, the precise measurement apparatus and the kind and shape of the material to be tested [8]. The transmissionline technique is the wellknown and popular in the literature. That technique can be used in transmission [79] and/or in reflection with openend configuration [10, 11] or shortcircuit [12]. All structures terminated by an open circuit and/or shortcircuit use the iterative procedure [13]. Probe technique which is used for liquid through ColeCole theory [14], but also through combination of cavity mixed with probe [1] and [6] has a good accuracy according to the material thickness. For thick material, the freespace technique is quite often used [15, 16]. But the resonator technique, through the use of cavity is the best and is welldeveloped in [17]. The sample to be characterized defines the technique through its state and sharp [18]. All techniques have advantages and inconveniences. One of them is discontinuity that stops the scanned frequency range, promotes the spread of highorder modes [19] and increases measurement errors through the intrinsic material complex parameters [20]. In order to contribute to the techniques improvement, we have developed a new technique called Movable ShortCircuited Line (MSCL). This technique is broadband and belongs to the reflection transmissionline technique [12]. The use of two transmissionlines is a good way to reduce discontinuities instead of choosing an equivalent electric circuit as their representation [21]. Those discontinuities exist at the junction of connector and the ideal test cell interface. Discontinuity represents imperfections of the entire test cell. From the two transmissionline technique using the wave cascade matrix principal [7, 22], we have developed the mathematic formulation for the MSCL. It is based upon the sum of measurements when using two different lengths of the test cell, but in inserting the first value of the studied band. That has allowed sorting out the discontinuity impacts on the extraction of the relative permittivity. Both circular coaxial line fixtures are differentiated only by their lengths while diameters and conductors are identical. The material under test is trapped inside the fixture, and sparameter measurements are done in two configurations: in presence and absence of the sample under test (SUT). The use of the propagation constant parameter through measurement of the reflection coefficient is the foundation of the technique. This is a quick, simple and highly reliable technique when the fixture is completely full of the material under test (MUT). The technique is suitable for destructive (wafer, etc.) as well as nondestructive (liquid, powder, etc.) samples.
2. Theory and mathematic model
2.1. Propagation constant determination
The propagation constant $\gamma $ is the unique parameter we can extract, and it must be applied in this technique. The shortcircuit reflection coefficient ${S}_{11}^{sc}$ is given in Eq. (1) as:
Fig. 1 describes the test fixture, where $Z$ and $L$ are respectively the cell and the connector lengths.
Fig. 1. A coaxial test fixture design
The propagation constant is determined from the reflection coefficient through the following expression:
The propagation constant is a complex parameter as defined below.
where $\theta =\beta l$ is the electric length and $\alpha l$ the attenuation constant. We assume that only a QuasiTEM mode propagates in and the vacuum configuration is linked to the filled structure by:
where ${\beta}_{d}$ and ${\beta}_{vac}$ are respectively the phase constants of the sample under test (SUT) and vacuum filled coaxial. In that case, the material relative parameters (permittivity and/or permeability) are obtained from equation (4) as:
We have easily filled up the circular coaxial test fixture with the nondestructive material, the one we used to validate the technique implementation.
2.2. Deembedding technique
The work requires the good knowledge of the connector propagation constant ${\gamma}_{0}$ as it is shown in the Fig. 2. and ${l}_{0}$ the connector length.
Fig. 2. A simplified fixture roadmap
The discontinuity impedance ${Z}_{d}$ appears at the connectorfixture interface. From ${Z}_{in}^{sc}$ and ${Z}_{in}^{oc}$, respectively, the input impedance of the connector when the connector is in shortcircuit and in opencircuit configurations. The connector propagation constant is as follows:
If we call by ${S}_{11}^{g}$ the entire reflection coefficient and $\Gamma $ the input cell reflection coefficient, both are linked by Eq. (7) as:
Because the main and only goal of this technique is to determine the relative permittivity ${\epsilon}_{r}$ of the material as mentioned in introduction, the focus of our work is to get that parameters by using the phase constant of the structure filled of vacuum and/or material under test (MUT).Whereas it is necessary to neglect the material losses, we have come with a mathematic equation to solve the negative number of the phase constant and discontinuities as in imperfection between the contact interface and ideal connector line as follows:
where ${\left(\beta l\right)}_{mes}$ is the phase constant of the test cell, measured in any aforementioned configurations, and ${\left(\beta l\right)}_{mes}^{iv}$ its initial value once measured in the work frequency range. ${\left(\beta l\right)}_{id}$ is the ideal phase constant as $\beta l>0$.
3. Technique description
Due to the high order mode propagations in the test cell (at interface of connector – ideal cell that leads imperfections) and the phase constant behavior which is the high order mode propagation consequence, the frequency range of study is limited [79]. On the other hand, the use of the shortcircuit configuration involves the use of iterative method [12] and [21]; we have put forward this new technique. There are two methodologies we have developed for the relative permittivity ${\epsilon}_{r}$ or permeability ${\mu}_{r}$ determination. Both of them are concerned with a model based on the transmissionline propagation constant. The determination can be based on a variety of criteria.
3.1. Use of one transmissionline technique
Let us consider a transmissionline cell, terminated by a short circuit in the end of the line. From the input impedance as illustrated below:
Fig. 3. Ideal transmissionline ended by shortcircuit
Eq. (9) is provided to calculate the input impedance once the line propagation constant is known.
Taking into account the reflection coefficient after deembedding the entire structure, the input impedance is given by:
If we suppose that the transmission line is lossless, in that case, $\gamma l\approx j\theta $, Eq. (9) will be rewritten. At the same time, we respectively named by ${\left({Z}_{in}^{sc}\right)}_{v}$ and ${\left({Z}_{in}^{sc}\right)}_{MUT}$ the input impedances in absence and presence of the MUT, which are linked to Eq. (9) by:
In this context, the iterative technique may be used. This solution can be watered down by its complexity through time taken to get to the aim. That’s why we have started by using a technique based on the average of two different extracted relative permittivity from two different lengths of a transmission line. If the test cell is homogenous like the one we manufactured and use to validate the technique process, then Eq. (11) becomes
and for the second transmissionline, we used the following equation in Eq. (12b),
We have noticed that the real relative permittivity or relative permeability is obtained once the average operation is done, as follows Eq. (13a):
We have named ${\left({\beta}_{k}{l}_{k}\right)}_{mes}^{x}$ and ${\left({\left({\beta}_{k}{l}_{k}\right)}_{mes}^{x}\right)}_{iv}$ the phase constants when the test cell is filled up of MUT or not and measured ($x$ is MUT or vacuum and, $k=1$ or $2$, the number of the cell that we used). If the sample to characterize is not magnetic, the final Eq. (13a) will approximately be reduced to:
The combination of both transmission lines can be done in another way. This is the technique we have named “Movable ShortCircuit Line” (MSCL).
3.2. Movable shortcircuit transmissionline technique
In this section, we combined two transmission lines, having the same characteristic impedance, geometry and dimensions (inner and outer diameters), but different lengths. The shortcircuit position is considered movable because of the fixture’s length as shown in Fig. 4.
Fig. 4. Two transmissionlines with shortcircuit at the bottom
In this situation, discontinuities are the same and, we assume that higher propagation modes cannot spread or propagate in. The ideal line, which is the length’s difference, $\u2206l={l}_{2}{l}_{1}$, has a propagation constant given in Eq. (14).
From Eq. (8) and supposing that ${\left({\beta}_{1}{l}_{1}\right)}_{mes}^{iv}\cong {\left({\beta}_{2}{l}_{2}\right)}_{mes}^{iv}$, which are the initial values of the phase constant when using the test cell corresponding to length ${l}_{1}$ and ${l}_{2}$ in each configuration (with or without MUT), we can write:
Associating Eqs. (12) and (15), the relative permittivity or permeability can be estimated through Eq. (16), given below as:
4. Technique improvement
This technique presents the peculiarity to behave according to the resonances of the waveguide. Electromagnetic simulations allowed us during the deembedding steps, to redefine the reflection coefficient of Eq. (7) by adding “$j$” as given in the following equation:
So, all equations that we got previously are combined to the new formulation of the reflection coefficient. In that case, Eq. (5) will be written as:
where “$\mathcal{g}$“ is a constant to be determined. We note by $\u2206{T}_{v}$ and $\u2206{T}_{m}$ the linearised phase constants when the test cell is in filled up with vacuum and MUT before correction, respectively, and by $\u2206{T}_{v}^{j}$ and $\u2206{T}_{m}^{j}$ once having corrected. By the way, we defined two elements $B$ and $D$ such as $D=\frac{\u2206{T}_{m}^{j}}{\u2206{T}_{m}}$ and $B=\frac{\u2206{T}_{v}^{j}}{\u2206{T}_{v}}$. In that case, Eq. (18) is written as:
Eq. (19) says that $\mathcal{g}={\left(\frac{\sum B}{\sum D}\right)}^{2}$ changes according to the material under test and the study frequency range. In addition, the fixture dimensions limit the study’s bandwidth, especially in frequencies lower than $2$ GHz. In other words, as long as the gap between the inner diameter “$b$” of the outer conductor and the outer diameter “$a$” of the inner conductor is small, low frequencies than 2 GHz are reachable. That is the main reason in the choice of the test cell in terms of shape and dimensions. In point of fact, we validated the technique implementation with aquarium sand, semolina and palm trees with two species: vinifera and laurentiis. Further reasons to be sure that the results are right, we compared them with those obtained by using the two transmissionline technique (2TL) when the test cells are exactly the same. We have named by ${\phi}_{1}={\beta}_{1}{l}_{1}$ and ${\phi}_{2}={\beta}_{2}{l}_{2}$ the linearised phase constant when the technique is used without improvement and by ${\phi}_{1}^{j}={\beta}_{1}^{j}{l}_{1}$ and ${\phi}_{2}^{j}={\beta}_{2}^{j}{l}_{2}$ once improved. In that case:
5. Technique validation
Measurements have been done from 20 MHz to 20 GHz with a particular attention on the frequency range 220 GHz. We have validated the technique with experiment on nonmagnetic materials such as: palm tree raffia, semolina and aquarium sand. A vector network analyzer (VNA), connected to two circular coaxial test cells, which are manufactured in brass, where $b$ = 14.38 mm and $a$ = 4 mm as it is shown in Fig. 4, with 80 mm and 100 mm of lengths, is the RF measurement equipment we used.
Fig. 5. The test cell design
5.1. Electromagnetic simulations
We have evaluated the technique relative error from the phase constant when the structure is filled up with vacuum. Using the circular coaxial dimensions that we have designed and manufactured, we compared both in good situation (absence of connector) and the real situation (there is a connector placed at a tip of the test cell).
Fig. 6. Phase constants of ideal line after extraction from electromagnetic simulations: use of Eq. (20a)
We have extracted the phase constants in both cases and can notice that both sketch curves are nearly the same and negatives from 20 MHz up to 1 GHz. In that case, we have applied Eq. (8) to solve that.
Fig. 7. Relative error’s curve from ideal and deembedded phase constants
Once the incertitude of the error was evaluated, we observed, as the following figure shows that the relative error decreases when frequency increases. Our main goal is to have a technique with an error of less than 10 % in the frequency range 220 GHz. The limits generated by the propagation modes now are stopped up to 20 GHz in our case when the MUT is vacuum.
5.2. Experimental validation
To validate the movable shortcircuit line technique implementation, we applied the extraction procedure to the biotechnology food such as semolina in frequency up to 5 GHz as it shown below.
Fig. 8. Relative permittivity result comparisons using two different equations
The use of Eq. (19) gives good results from 1.7 GHz while both techniques (basic and improvement ones) start to approach each another from 3 GHz. We have compared the improvement technique results to those obtained with the use of two transmissionlines in reflection/transmission [9].
Fig. 9. Relative permittivity result comparisons using two different techniques
These results are closes to each another and prove that the movable shortcircuit line technique, in its improvement version, is suitable for material relative permittivity extraction. The appropriateness or suitability of the low frequency where the technique can be acceptable is around 1.6 GHz. Because of that, we have extracted other material relative permittivities such as aquarium sand and palm raffia in order to certify that assertion.
Fig. 10. Relative permittivity results using the MSCL technique from Eq. (19)
As we focused on frequency range 220 GHz, it is observed that the movable shortcircuit technique has allowed covering frequencies upper than 10.4 GHz. Theoretically, when only vacuum is inside of this fixture, high order modes are generated from 10.4 GHz. Eq. (20), given in [21] is satisfied.
where “$b$” and “$a$”, inner diameter of the outer conductor and outer diameter of inner conductor, respectively, are in mm. The next sketched curves are used to compare results between two different techniques: the two transmissionlines (2TL) and MSCL.
Fig. 11. Comparison of relative permittivity of a) aquarium sand from MSCL and 2TL techniques and b) Laurentiis from MSCL and 2TL techniques
a)
b)
The two transmissionline results are obtained with 46 % of relative error in the entire frequency range. The results are in conformity with those from electromagnetic simulation ones. The correction on the return loss coefficient allowed improving the extraction process and the technique implementation to cover frequencies 1.73.5 GHz.
6. Conclusions
In spite of the multitude of techniques which suffer from frequency limits because of the high order mode propagations that are generated by the test cell discontinuities, we have presented in this paper how to overcome that difficulty through the movable broadband shortcircuit technique. This new way to extract the relative permittivity and/or permeability parameter is well adapted for industrial applications. We avoided using iterative solutions and suggested how to calculate automatically the correction factor according to the MUT and the scanned frequency range. The electromagnetic simulations have proved that the error rate is inversely proportional to frequency. It means that the more the covered frequencies increase, the better is the accuracy, which covers 7% to 3% when using homogenous configuration. We validated the technique implementation with some biotechnology materials: two palm tree species (vinifera and laurentiis), aquarium sand and semolina in the frequency range 218 GHz with a good accuracy. We have proved that the technique result can be repeated as long as the calibration of the connector is well made, especially its short circuit in microwave domains. Nevertheless, despite the progress we have made, the technique is not suitable for frequencies lower than 1.7 GHz, due to the test cell (dimensions or shape) and/or the mathematical formulation. This means that developing a new fixture or a new technique in order to reach low frequencies will be suitable.
Acknowledgements
Authors would like to thank in a particular way F. Kuhlman for helpful suggestions and given samples to validate the new technique, also the anonymous reviewers for their attentive revision and thorough assessment.
References
 Zarral L., Djahli F., Ndagijimana F. Technique of coaxial frame in reflection for the characterization of single and multilayer materials with correction of air gap. International Journal of Antennas and Propagation, Vol. 15, Issue 5, 2014, p. 19. [Search CrossRef]
 Peyman A. Holden S. J., Watts S., Perrott R., Gabriel C. Dielectric properties of porcine cerebrospinal tissues at microwave frequencies: in vivo, in vitro and systematic variation with age. Physics in Medicine and Biology, Vol. 52, Issue 8, 2007, p. 22292245. [Publisher]
 Hand J. W. Modelling the interaction of electromagnetic fields (10 MHz – 10 GHz) with the human body. Physics in Medicine and Biology, Vol. 53, Issue 16, 2008, p. 243286. [Publisher]
 CabanesSempere M., et al. Characterization method of dielectric properties of free falling drops in a microwave processing cavity and its application in microwave internal gelation. Measurement Science and Technology, Vol. 24, Issue 9, 2013. [Publisher]
 La Gioia A., O’Halloran M., Porter E. Modelling the sensing radius of a coaxial probe for dielectric characterization of biological tissues. IEEE Access, Vol. 6, 2018, p. 4651646526. [Search CrossRef]
 Moukanda M., Ndagijimana F., Chilo J., Saguet P. Coaxial probe fixture used for complex permittivity measurement of thin layers. African Physical Review, Vol. 2, Issue 31, 2008, p. 6567. [Search CrossRef]
 ReynosoHernandez J. A. Unified method for determining the complex propagation constant of reflecting and nonreflecting transmission lines. IEEE Microwave and Wireless Components Letters, Vol. 13, Issue 8, 2003, p. 351353. [Publisher]
 Takach A. A., Moukanda M. F., Ndagijimana F., AlHusseini M., Jomaa J. Twoline technique for dielectric material characterization with application in 3Dprinting filament electrical parameters extraction. Progress in Electromagnetic Research M., Vol. 85, 2019, p. 195207. [Publisher]
 Moukanda Mbango F., Al Takach A., Ndagijimana F. Complex relative permittivity extraction technique of biotechnology materials in microwave domains. International Journal of Electronics, Communication and Instrumentation Engineering Research and Development, Vol. 9, Issue 1, 2019, p. 3342. [Publisher]
 Mosig J. R., Besson J. C. E., GexFarby M., Gardiol F. E. Reflection of an openended coaxial line and application to nondestructive measurement of materials. IEEE Transaction on Instrumentation and Measurement, Vol. 30, 1981, p. 4651. [Publisher]
 Seaman R., Burdete E., Dehaan R. Openended coaxial exposure device for applying RF/microwave fields to very small biological preparations. IEEE Transaction on Microwave Theory and Techniques, Vol. 37, 1989, p. 102111. [Publisher]
 James BakerJarvis Transmission/reflection and shortcircuit line permittivity measurements. National Institute of Standards and Technology, 1990, p. 6876. [Search CrossRef]
 Wolfson B. J., Wentworth S. M. Complex permittivity and permeability measurement using a rectangular waveguide. Microwave and Optical Technology Letters, Vol. 27, Issue 3, 2000, p. 180182. [Publisher]
 Marsland T., Evans E. Dielectric measurements with an openended coaxial probe. IEEE Proceeding, Vol. 134, 1987, p. 341349. [Search CrossRef]
 Ghodgaonkar D. K., Varadan V. V., Varadan V. K. Freespace measurement of complex permittivity and complex permeability of magnetic materials at microwave frequencies. IEEE Transaction on Instrumentation and Measurement, Vol. 39, Issue 2, 1990, p. 387394. [Publisher]
 Sagnard F., Bentabet F., Vignat C. In situ measurements of the complex permittivity of materials using reflection ellipsometry in the microwave band: Experiments (Part II). IEEE Transaction on Instrumentation and Measurement, Vol. 54, Issue 3, 2005, p. 12741282. [Publisher]
 Salahum E., Queffelec P., Le Floch M., Gelin Ph. A Broadband Permeameter for in Situ Measurements of Rectangular Samples. IEEE Transactions on Magnetics. Vol. 37, Issue 4, 2001, p. 27432745. [Publisher]
 Pozar D. M. Microwave Engineering, 3rd Edition. John Wiley and Sons, Inc., Hoboken, 2005. [Search CrossRef]
 Aleksandrov N. L., et al. Selective excitation of high order modes in circular waveguides. International Journal of Infrared and Millimeter Waves, Vol. 13, Issue 9, 1992, p. 13691385. [Publisher]
 Moukanda Mbango F, Ndagijimana F. Electric parameter extractions using a broadband technique from coaxial line discontinuities. International Journal of Scientific Research and Management, Vol. 7, Issue 5, 2019, p. 248253. [Search CrossRef]
 Moukanda Mbango F. Contribution à la Caractérisation Electrique de Matériaux Utilisés en Microélectronique Radiofréquence. Thesis of Université Joseph Fourier, Grenoble, 2008. [Search CrossRef]
 Janezic M. D., Jargon J. A. Complex permittivity determination from propagation constant measurements. IEEE Microwave and Guided Wave Letters. Vol. 9, Issue 2, 1999, p. 7678. [Publisher]
Cited By
IETE Journal of Research
F. Moukanda Mbango, M. G. Lountala, F. Ndagijimana, D. LilongaBoyenga

2021
