Research Article  Open Access
Raúl V. HaroBáez, Juan Córcoles, Jorge A. RuizCruz, José R. MontejoGarai, Jesús M. Rebollar, "HigherOrder Mode Electromagnetic Analysis of a Material Sample between Two Flanged Coaxial Probes for Broadband Modelling of Dielectric Measurement Setups", Advances in Mathematical Physics, vol. 2019, Article ID 6404812, 17 pages, 2019. https://doi.org/10.1155/2019/6404812
HigherOrder Mode Electromagnetic Analysis of a Material Sample between Two Flanged Coaxial Probes for Broadband Modelling of Dielectric Measurement Setups
Abstract
The characterization of the dielectric properties of a material requires a measurement technique and its associated analysis method. In this work, the configuration involving two coaxial probes with a material for dielectric measurement between them is analyzed with a modematching approach. To that effect, two models with different complexity and particularities are proposed. It will be shown how convergence is sped up for accurate results by using a proper choice of higherorder modes along with a combination of perfect electric wall and perfect magnetic wall boundary conditions. It will also be shown how the frequency response is affected by the flange mounting size, which can be, rigorously and efficiently, taken into account with the same type of approach. This numerical study is validated through a wide range of simulations with reference values from another method, showing how the proposed approaches can be used for the broadband characterization of this wellknown, but with a recent renewed interest from the research community, dielectric measurement setup.
1. Introduction
The behavior of a material under the exposure to electromagnetic fields is determined by its dielectric properties. The accurate knowledge of these properties is essential in science and engineering, since they are used in a variety of industries [1, 2]. The properties of a material can be tracked during an industrial process, and, thus, they can be used to monitor a product, incorporating that information into an improved quality control. They are also important in diverse fields such as remote sensing, heating, molecular analysis, and creation of liquid standards [3, 4].
More specifically, in microwave engineering [5, 6], the dielectric properties of materials are key data for evaluating the performance of circuit substrates, printed antennas, radomes, dielectric resonators, etc. [7]. In modern computer aided design, these properties are incorporated in the socalled fullwave methods (techniques to solve the Maxwell equations starting from the full electromagnetic model) in order to provide accurate predictions for each application. In the telecommunications industry, with a massive use of the microwave spectrum, the precise characterization of the materials involved in the microwave components is essential to optimize and make designs with stringent electrical specifications, which are becoming progressively more demanding to address new systems with higher data capacity.
Regardless of the specific area where the dielectric properties are required, the process to characterize them is based on two techniques: (i) a measurement technique extracting a parameter related directly or (more commonly) indirectly with the dielectric properties and (ii) an electromagnetic analysis technique (circuitbased or fullwave) modelling the measurement, i.e., providing estimations of the results of that measurement for a given material.
With respect to the measurement techniques, there is a plurality of options, such as those based on coaxial probes, transmission lines, free space methods, use of resonators (as in microwave filter characterization), etc. [1]. This variety of measurement techniques provides different alternatives for dealing with constrains in the frequency range, the requirements of preparation for the sample, type of expected (low/high) losses in the material, etc. In this work, the coaxial probe will be used, since it provides a good tradeoff between many of these features, can be used for biological tissues, and can be considered a broadband method, which is one of the goals of this paper. In addition, its use is very simple and wellknown, especially in microwave engineering, where it has been used for a long time [8–11].
With respect to the analysis technique used in the dielectric characterization process, when the electromagnetic analysis is based on a fullwave approach, the results will be more precise, at the expense of a higher computational burden than simpler circuit analyses. This paper will use an accurate fullwave electromagnetic method based on modematching methods [12, 13], which provide more accurate estimations than circuital methods. In comparison with other fullwave methods [14], it also shows a great efficiency, which will be further improved by exploiting the specific features of the geometry associated with the twoport coaxial probe measurement problem.
In fact, this kind of problem has received very recent attention for different types of materials and with different types of probes for the input/output [15–17]. In this paper, two different models based on modematching along the longitudinal direction are proposed for characterizing the twoport problem associated with the coaxial probe for dielectric measurement of homogeneous samples. The first model will benefit from different types of boundary conditions that can be used to model the simplest structure, made up of two coaxial lines with the material to be measured inserted in a virtual circular waveguide that contains the sample. This circular waveguide will be modelled with perfect electric or magnetic wall boundary conditions and with a combination of both cases, which will be shown to provide faster convergent results. The second model will also take the flanges of the setup into account. The convergence of these models with respect to different parameters will be studied, resulting in a very efficient tool that can be used for the broadband analysis of the twoport scattering problem of a material sample placed between two coaxial probes for dielectric measurement.
2. The Coaxial Probe for Dielectric Measurement
2.1. OnePort Open Coaxial Probe
The coaxial probe open at one end, whose basic schematic version is shown in Figure 1, is one of the preferred techniques for measuring materials. It is a very simple setup, where the sample material can be inserted in the test unit with almost no preparation, filling the space just after the coaxial probe. This can be easily done for a large variety of malleable materials such as liquids or semisolids whose shape can be adjusted. However, for solids, contact surfaces must be flat enough and tighten to the test unit in order to avoid air gaps between the sample and the probe that will degrade the measurement.
(a)
(b)
From the electromagnetic point of view, the problem in Figure 1 is a coaxial probe radiating in a semiinfinite space (the orange block in the figure) whose dielectric properties are those of the material under test (MUT), which may have losses: complex dielectric permittivity ϵ=ϵ′jϵ′′ and complex magnetic permeability μ=μ′jμ′′. The coaxial transmission line has two conductors of internal and external radius a and b, respectively, with a supporting dielectric (blue color in Figure 1). These parameters are typically adjusted for providing 50 Ω characteristic impedance, as in classic microwave equipment. In the analysis technique, a semiinfinite space is assumed for the region related to the MUT, and, thus, the thickness (d) and the transversal size of the sample (R_{sample}) must be electrically large enough in order to confine the incoming electromagnetic field excited by the probe within the MUT limits. Within the same modelling, the metallic flange is assumed to be large enough for avoiding backward radiation.
The measurement technique for this case extracts the reflected electromagnetic wave at the oneport problem. The signal excited by the generator connected to the coaxial propagates through the transmission line and reaches the material, where it is partially reflected back. This reflection coefficient depends on the electromagnetic properties of the dielectric material. Thus, from this measurement, the values of the dielectric properties of the material in contact with the probe can be estimated. However, there is not a direct closed formula relating the ϵ, μ of the MUT (this is the goal of the dielectric characterization process) to the measured reflection coefficient, and this is where the different analysis techniques, from equivalent circuits to fullwave approaches, play a key role. It is important to note that this is also a classic problem in microwave engineering when designing antenna feed systems, and, thus, there has been continuous research in this subject from different perspectives during the last decades [8–11, 18–21]. These models have evolved between the complexity of the lumped and distributed elements in the initial equivalent circuits, to the way of approximating the electromagnetic field functions in the more recent fullwave approaches.
2.2. Proposed Models for the Characterization of the TwoPort Problem
The twoport measurement is a method directly related to the coaxial probe [22–24], where a thin sample of the MUT is inserted between two coaxial lines as in Figure 2, allowing a quick measurement of the return and insertion loss as a function of frequency. The analysis methods are very closely related to the coaxial probe discussed in the previous subsection. Nevertheless, different formulations have been the subject of research in very recent years [15–17].
(a)
(b)
In this work, the two models shown in Figure 3 will be proposed for MUTs with homogeneous dielectric properties. The simplest model A is shown in Figures 3(a) and 3(b). It will allow a straight analysis through modematching, since it involves a problem with two waveguide discontinuities between coaxial and circular waveguides [25] surrounded by perfect electric wall (PEW) at its boundaries. A variation will be introduced to account for the particularities of the circular waveguide that contains the sample. The computations will be also done with a waveguide having perfect magnetic wall (PMW) at its boundaries and with a combination of both PEW/PMW. Their convergence properties will be tested, allowing us to decide the faster method for the twoport problem characterization.
An alternative model B is shown in Figures 3(c) and 3(d), which takes the finite flanges associated with the input/output coaxial lines into account. This model can also be simulated with modematching, as long as the problem of waveguide bifurcations is included in the formulation. This type of approach has been used in the past for probes and, more recently, to improve the modematching accuracy in radiation problems [11, 26, 27]. In this work, the model which incorporates the finite flange information will be proposed for the twoport coaxial holder. It will exploit, as model A, the advantage of the circular revolution, since only modes with no circular variation will be needed to solve this problem.
3. Electromagnetic Analysis of the Proposed Models through ModeMatching
The proposed two models are analyzed with the modematching method, comprised of two main stages. First, the modes of all the waveguides involved in the problem are calculated. Then, the modal series to expand the fields in each waveguide are matched at each side of the step discontinuity, providing the Generalized Scattering Matrix (GSM) of each step. Finally, the GSMs of each step are cascaded [14, 28]. The main steps are briefly described now.
3.1. Characterization of the Waveguides
The electromagnetic field inside homogeneous waveguides can be described in terms of transversal electromagnetic (TEM), transversal electric (TE), and transversal magnetic (TM) modes [29]. The electromagnetic field of a TE or TM mode is derived from a scalar function which is the solution to the Helmholtz equation with Neumann or Dirichlet boundary conditions:where is the transversal to z (two dimensional) nabla differential operator, S is the crosssection with normal , and C is its contour. In our case for the models in Figure 3, only coaxial and circular waveguides are needed, as shown in Figure 4. In (1), function is the longitudinal zcomponent, up to a constant factor (not depending on the spatial point), of the magnetic or electric modal field, of the TE or TM mode, respectively. The cutoff frequency of the mode is directly related to the cutoff wavenumber of the mode [7, 29].
For the coaxial waveguide, the TEM mode is also necessary to build the complete mode spectrum of the waveguide, and it can be easily computed analytically by solving the Laplace equation, where represents the potential function that provides the electric field: Once is solved for TE, TM, and TEM modes, the transversal components are obtained using basic gradient operations [7, 29]. When the crosssection of the waveguide belongs to a canonical case (such as rectangular or circular), the function is readily obtained by applying the technique of the separation of variables. In our case, it is wellknown that the longitudinal field of the TE and TM modes for the waveguides in Figure 4 is written in terms of the Bessel function of the first () and second kind ():where is related to the power carried by the mode and is thus used as a normalization factor. For a hollow circular waveguide, in (3). The cutoff wavenumbers of the modes for a PEW circular waveguide are , , and are the zeros of and , respectively. For a PMW boundary at the circular waveguide, the dual solution with , is obtained [29]. This case for PMW is also found in [27], where the virtual waveguide with PMW is used for improving the field approximation in horn problems. For a coaxial waveguide in (3), is a combination of Bessel functions whose amplitudes and the cutoff wavenumber are obtained by imposing the boundary condition at the contour [7, 29].
It is important to note that the cutoff wavenumbers for the coaxial waveguide must be computed for each specific aspect ratio b/a. Nevertheless, in our problem, since the structure will be excited with a TEM mode, which does not have angular variation and has PMW symmetry at both xz and yz planes, all the higherorder modes generated in the structure are modes in both the circular and coaxial waveguides.
3.2. ModeMatching Formulation
As every waveguide in the two models of Figure 3 is surrounded by PEW (or PMW), the transverse electromagnetic field in these structures can be expanded in a set of orthogonal TE, TM, and TEM modes formulated in the previous section (although TE modes will not be needed as discussed in previous section, they will be kept in this formulation for the sake of the generality). With them, the electromagnetic field in each point of the structure is represented as a modal series.
In the reference plane (see Figure 5), the fields can be written as a sum of orthonormal modes with the incident and reflected amplitudes at the discontinuity plane. In each waveguide (index g=w, for the larger, g=s for the smaller), the following can be written:where , are the transverse modal vectors corresponding to the mth TE, TM, or TEM mode of the guide (g), which can be propagating or evanescent. These fields take the lateral boundary conditions into account while the electric and magnetic modal functions are related by the mode impedance as usual, which are also typically normalized to unity as follows. The electric and magnetic field boundary conditions (EFBC and MFBC, respectively) for this problem relate the transversal fields in (4) at both sides of the discontinuity.These equations are imposed by means of a classic Galerkin method where the EFBC is tested by magnetic modal fields from the larger waveguide, and the MFBC is tested with electric modal fields in the smaller waveguide. This fieldmatching procedure also requires truncating the series (4) to terms and collecting the modal amplitudes in column vectors. The final equations expressed in matrix form becomewhere matrix is filled with the interaction between the electric and magnetic modal functions from both waveguides in the step:The equation system presented in (7) provides, after its resolution, the Generalized Scattering Matrix (GSM) of the problem in the form:with:This formulation has been used in the past many times. Examples of the innercross product computation for different waveguide steps can be found in [25, 30, 31]. After this computation, the cascading of the GSM of each step in the problem is carried out systematically [14, 28], as it will be explained in the next section for the two proposed models.
(a)
(b)
4. Results for the Proposed Models
The modematching technique described in the previous section has been applied for the two proposed models in Figure 3, further detailed in this section. The first model uses two coaxial waveguides attached to both sides of a virtual circular waveguide filled by the sample material, which is denoted as model A. The second model takes into account the finite flanges used in the measurement, and it is called model B.
The test cases will be evaluated under the two models A and B, for the broadband frequency range from 0.1 GHz to 20 GHz, varying their constituent parameters. The coaxial waveguide used for the measurement has an inner and an outer radius of a=0.46mm and b=1.75mm, respectively, filled by Teflon () for providing 50 Ω. The sample material for the tests has , () and , (). This material is prepared with a thickness of d=0.2mm and d=1mm. These main data are shown in Table 1.

The reference values are provided by the method from the Dielectric Assessment Kit (DAK) product line manufactured by the Swiss company SPEAG [32], used in highprecision systems for measurements of dielectric properties. These instruments allow dielectric parameter measurements over a very broad frequency range.
4.1. Model A with the Step between the Coaxial and the Large Circular Waveguide Containing the Sample
Model A is detailed in Figure 6. This model will be also divided into three cases (A1, A2, A3), depending on the boundary conditions used for the outer boundary of the virtual circular waveguide containing the sample. In all these cases, from the analysis point of view, the model only requires the characterization of two GSMs, called and according to Figure 6(b). However, and are only different in the port numbering, but both represent a step between the same types of waveguides with the same dimensions; the computational cost reduces to computing a single GSM, since once is obtained, is also available. Then, they have just to be cascaded through wellknown operations to get the GSM representing the full structure.
4.1.1. Case Model A1: Large Virtual Circular Waveguide with PEW
Since the problem has no angular variation, under the TEM excitation by the coaxial line, the series (4) in the modematching will only require the TEM mode at the coaxial lines, and modes at the coaxial and the circular waveguides. The field in the circular waveguide will be always represented in both models A and B by modes. However, in order to see the importance of the higherorder modes at the input/output, we will first compute in Figure 7 the results by only using the TEM mode to approximate the field within the coaxial lines, as initial pioneer works with probes did.
In the simulation in Figure 7, the magnitude of the S_{11} and S_{21} parameters is shown for different values of the radius R_{sample} of the virtual circular waveguide between the input and output coaxial lines. The analysis starts with R_{sample} equal to 20 times the outer radius of the input/output coaxial, which shows a ripple that starts to decrease when this circular waveguide radius is enlarged. In this case, when R_{sample} increases, the results in both reflection and transmission are flatter and converge to the larger case of R_{sample}=200b. However, the results do not reproduce the results of the DAK reference values, showing a relevant difference in the curves.
In order to overcome this issue, Figure 8 shows the same type of analysis as Figure 7, but using also modes in the coaxial transmission lines in addition to the TEM mode. Although the ripple also appears as in the previous case and decreases when R_{sample} increases as well, now the results converge to those of the reference case obtained by DAK. Thus, the response improves the results obtained from the previous case, having a very good agreement. The phase is also shown in this graph, showing also a very good agreement. The phase reference is taken at the interface between the coaxial and the circular waveguide in all the cases.
The same type of analysis is now done with the sample having a thickness of 1mm. From now on, the analysis will be always performed with TEM and modes in the coaxial lines, since its need has been clearly justified in Figures 7 and 8. The results for the new thickness d=1mm are shown in Figure 9, where the same type of performance as with the thinner sample can be seen: ripple in the curves related to the size of the large virtual circular waveguide containing the sample. This ripple reduces when the size of R_{sample} increases till achieving convergent results matching those of the DAK reference.
4.1.2. Case Model A2: Large Virtual Circular Waveguide with PMW
In model A, we are expecting that the electromagnetic field at the lateral limit of the circular waveguide (i.e., for large ρ) is small enough to represent a suitable model for the actual measurement, where it is assumed that the field is confined within the MUT limits. In order to confirm that, the model should not provide very different results when the lateral boundary conditions of the circular waveguide are changed to PMW and the waveguide steps in Figure 6(b) are recomputed with these new modes. This is confirmed with the results shown in Figure 10 for a sample thickness of 1mm, where the magnitude of S_{11} and S_{21} has a very similar behavior to that of the case of model A1 with PEW in Figure 9. There is also a ripple in the modematching results, more important for smaller R_{sample}, which disappears when R_{sample} is large enough.
In all the simulations done so far, we have presented the results with a number of modes large enough to achieve convergent results with respect to the number of terms in the series (4). In this type of problems, as a general rule to overcome the relative convergence problem [14], the number of modes in the modal series is selected following a criterion related to the size of the involved waveguides. One possible solution is keeping a relation between the modes of the different waveguides as the square root of their areas. This means that for a step between circular waveguides, the ratio of the modes between two waveguides in a step discontinuity would be the ratio of their radii.
The convergence has been verified with several analyses, varying the number of modes while keeping a fixed geometry. One of those analyses is shown in Figure 11. It shows the results of this last case in Figure 10, although the behavior would be very similar for other geometries and for both models A or B. In Figure 11, the black solid line represents the reference and the other lines represent the results for different number of modes in the coaxial lines N_{coax} and in the large circular waveguide N_{circ}.
4.1.3. Case Model A3: Large Virtual Circular Waveguide with Combination of PEW and PMW
Figure 12 shows the results of models A1 and A2 with R_{sample}=10b plotted within the same graph. It can be seen that the ripple of the responses for the two cases (PEW/PMW at the outer lateral boundary of the circular waveguide, respectively) are interleaved. This would also happen in the other results of Figures 9 and 10, although it is more difficult to see when the ripple starts to decrease for large R_{sample}. In addition, since both boundary conditions in the circular waveguide are approximations of the actual measurement problem, the combination of both results could be also another suitable analysis model. Moreover, this type of strategy has also been used in [26].
Following this rationale, Figure 12 also shows the combination of the PEW analysis in model A1 with those using the PMW in model A2, which will be the model A3 proposed in this subsection: In order to verify model A3, Figure 13 shows the results for different values of . It can be seen how the combination of the two boundary conditions converges to the reference results with smaller , in both magnitude and phase. This translates into a higher computational efficiency, since a smaller number of modes will be needed to get the results using this smaller circular waveguide (in comparison with models A1, A2) containing the sample.
(a)
(b)
4.2. Model B with the Steps for Taking the Finite Flanges into Account
The second model to compute the twoport problem is shown in Figure 14. This geometry takes the finite flanges into account, as well as the outer size of the coaxial lines involved in the measurement. This may be relevant for some cases where the assumption that all the fields are confined within the MUT limits starts to be not fully compliant with the experimental data. The structure has been plotted symmetrically with respect to the middle plane for the sake of the simplicity, but data from the input coaxial could be made different from the output coaxial.
From the modematching point of view, the new geometry involves two new virtual coaxial transmission lines at each side of the problem, associated with the incorporation of the finite flange with radius R_{flange} and length to the model. They have an inner radius of R_{o} and R_{flange}, respectively, both with outer radius of R_{sample}. They are filled by air.
This geometry involves a new step discontinuity between the two virtual airfilled coaxial lines associated with the outer part of the structure, as is shown in Figure 14(b) with the block []. Model B also involves a bifurcation between, at one side, the Teflon coaxial (radii a and b) and one airfilled coaxial (radii R_{flange} and R_{sample}) and, at the other side, the larger circular waveguide of radius R_{sample}. In this problem, two additional virtual ports appear, carrying the signal not confined within the MUT limits and not taken into account by model A. Once the GSM of the steps and bifurcations are computed, they are also easily cascaded.
For the following tests, in addition to the data in Table 1, the flange length will be L_{flange}=1mm, with a cover for the input/output coaxial of 0.5mm thickness, i.e., R_{o}=b+0.5mm. As was done for model A, the first simulation in Figure 15 is done using a fixed geometry with finite R_{flange} =5.4b, varying the radius of the R_{sample} associated with the circular waveguide between 10b and 100b. In this case, although the results are very similar for the different R_{sample}, the finite size of the flange is leading to results different from those of the DAK reference, regardless of the size of the virtual circular waveguide required in the modematching simulations. The DAK values corresponded to the typical geometry with a large enough flange and electromagnetic field fully confined within the MUT limits, which is coherent with results of model A.
This effect of the size of the flange is further analyzed in Figure 16, where the simulations are run with different number of modes and different values for R_{flange}. The ripple is now associated with R_{flange} since these simulations have a fixed value of R_{sample}. It is also important to note that additional terms in series (4) are needed for having convergent results in terms of number of modes, since now the structure is more complex involving more waveguide steps.
In this case, since the new airfilled coaxial lines have a large size, many modes will be needed to approximate the field within them, with a high aspect ratio between the outer and inner conductor. Thus, their modes will have to be computed carefully to avoid instabilities related to the behavior of the function in (3) when it is evaluated close to zero.
Figure 17 shows the results of model B for different values of R_{flange} in Figure 14(a), with a fixed value of . Results show that for a large flange, the DAK reference is achieved. Figure 17 also shows the backward and forward radiation for different values of the flange, characterized by the Sparameter between the input (or output) Teflon coaxial port and the virtual airfilled coaxial ports 3 and 4 in Figure 14. As expected, for a smaller flange, large peaks of signal at specific frequencies are not confined within the MUT limits. Figure 18 shows the same type of analysis of Figure 17 for a thickness of 1mm, where it can be seen that these peaks of forward and backward radiation are less significant.
5. Conclusions
The mathematical formulation of the modematching method applied to the twoport scattering problem arising in a material placed between two coaxial probes has been presented. This wellknown configuration has received renewed attention in very recent years. In this line, this paper has developed two different approaches to characterize this dielectric measurement setup for broadband applications. The analysis of the problem with the MUT considering both PMW and PEW has been studied, while their combination has been shown to improve the efficiency of the method. In all cases, the use of higherorder modes in the coaxial waveguide at the MUT interface is essential to obtain fast and accurate results. The effect of the finite flanges has been also efficiently and rigorously assessed, proving that a flatter frequency response is achieved when the field is confined within the MUT, i.e., when large flanges are considered. Two setups have been analyzed and their results have been compared with the reference data provided by the DAK characterization method, showing excellent agreement between them.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported in part by the Spanish Government (Agencia Estatal de Investigación, Fondo Europeo de Desarrollo Regional) under Grant TEC201676070C31/2R (AEI/FEDER, UE) and in part by the program of Comunidad de Madrid under Grant S2013/ICE3000. The authors want to acknowledge SPEAG company, Zurich, Switzerland, especially Pedro CrespoValero and Ferenc Muranyi, for their assistance with the DAK results.
References
 Keysight Technologies, “Basics of measuring the dielectric properties of materials,” Application Note 59892589EN, USA, March 7, 2017. View at: Google Scholar
 V. Komarov, S. Wang, and J. Tang, “Permittivity and measurement,” in Encyclopedia of RF and Microwave Engineering, pp. 1–20, John Wiley & Sons, 2005. View at: Google Scholar
 U. Kaatze and Y. Feldman, “Broadband dielectric spectrometry of liquids and biosystems,” Measurement Science and Technology, vol. 17, no. 2, pp. R17–R35, 2006. View at: Publisher Site  Google Scholar
 A. P. Gregory and R. N. Clarke, “A review of RF and microwave techniques for dielectric measurements on polar liquids,” IEEE Transactions on Dielectrics and Electrical Insulation, vol. 13, no. 4, pp. 727–743, 2006. View at: Publisher Site  Google Scholar
 R. W. Bruce, “New Frontiers in the Use of Microwave Energy: Power and Metrology,” in Proceedings of the Materials Research Society Symposium, vol. 124, pp. 3–15, 1990. View at: Publisher Site  Google Scholar
 L. F. Chen, C. K. Ong, C. P. Neo, V. V. Varadan, and V. K. Varadan, Microwave Electronics: Measurement and Materials Characterization, John Wiley & Sons, Ltd, 2004.
 D. M. Pozar, Microwave Engineering, John Wiley, Hoboken, NJ, USA, 2011.
 R. W. P. King and L. D. Scott, “The Cylindrical Antenna as a Probe for Studying the Electrical Properties of Media,” IEEE Transactions on Antennas and Propagation, vol. 19, no. 3, pp. 406–416, 1971. View at: Publisher Site  Google Scholar
 E. C. Burdette, F. L. Cain, and J. Seals, “In vivo probe measurement technique for determining dielectric properties at VHF through microwave frequencies,” IEEE Transactions on Microwave Theory and Techniques, vol. 28, no. 4, pp. 414–427, 1980. View at: Publisher Site  Google Scholar
 J. R. Mosig, J.C. E. Besson, M. GexFabry, and F. E. Gardiol, “Reflection of an openended coaxial line and application to nondestructive measurement of materials,” IEEE Transactions on Instrumentation and Measurement, vol. IM30, no. 1, pp. 46–51, 1981. View at: Publisher Site  Google Scholar
 J. D. BañosPolglase and J. D. Rebollar, “Coaxial probes with finite ground plane,” IEEE Electronics Letters, vol. 24, no. 5, pp. 291292, 1988. View at: Publisher Site  Google Scholar
 P. J. B. Clarricoats and K. R. Slinn, “Numerical method for the solution of waveguidediscontinuity problems,” IEEE Electronics Letters, vol. 2, no. 6, pp. 226–228, 1966. View at: Publisher Site  Google Scholar
 A. Wexler, “Solution of waveguide discontinuities by modal analysis,” IEEE Transactions on Microwave Theory and Techniques, vol. 15, no. 9, pp. 508–517, 1967. View at: Publisher Site  Google Scholar
 T. Itoh, Numerical Techniques for Microwave and MillimeterWave Passive Structures, John Wiley, New York, NY, USA, 1989.
 M. W. Hyde and M. J. Havrilla, “A broadband, nondestructive microwave sensor for characterizing magnetic sheet materials,” IEEE Sensors Journal, vol. 16, no. 12, pp. 4740–4748, 2016. View at: Publisher Site  Google Scholar
 M. W. Hyde, M. J. Havrilla, and A. E. Bogle, “Nondestructive determination of the permittivity tensor of a uniaxial material using a twoport clamped coaxial probe,” IEEE Transactions on Microwave Theory and Techniques, vol. 64, no. 1, pp. 239–246, 2016. View at: Publisher Site  Google Scholar
 M. H. Hosseini, H. Heidar, and M. H. Shams, “Wideband nondestructive measurement of complex permittivity and permeability using coupled coaxial probes,” IEEE Transactions on Instrumentation and Measurement, vol. 66, no. 1, pp. 148–157, 2017. View at: Publisher Site  Google Scholar
 D. K. Misra and M. Chabbra, “A quasistatic analysis of openended coaxial lines,” IEEE Transactions on Microwave Theory and Techniques, vol. 35, no. 10, pp. 925–928, 1987. View at: Publisher Site  Google Scholar
 J. BakerJarvis, M. Janezic, P. Domich, and R. Geyer, “Analysis of an openended coaxial probe with liftoff for nondestructive testing,” IEEE Transactions on Instrumentation and Measurement, vol. 43, no. 5, pp. 711–718, 1994. View at: Publisher Site  Google Scholar
 D. V. Blackham and R. D. Pollard, “An improved technique for permittivity measurements using a coaxial probe,” IEEE Transactions on Instrumentation and Measurement, vol. 46, no. 5, pp. 1093–1099, 1997. View at: Publisher Site  Google Scholar
 B. GarcíaBaños, J. M. CataláCivera, A. J. Canós, and F. PeñarandaFoix, “Design rules for the optimization of the sensitivity of openended coaxial microwave sensors for monitoring changes in dielectric materials,” Measurement Science and Technology, vol. 16, no. 5, pp. 1186–1192, 2005. View at: Publisher Site  Google Scholar
 R. A. Weck, “ThinFilm Shielding for Microcircuit Applications and a Useful Laboratory Tool for PlaneWave Shielding Evaluations,” IEEE Transactions on Electromagnetic Compatibility, vol. 10, no. 1, pp. 105–112, 1968. View at: Publisher Site  Google Scholar
 A. R. Ondrejka and J. W. Adams, “Shielding effectiveness (SE) measurement techniques,” in Proceedings of the IEEE National Symposium on Electromagnetic Compatibility, pp. 1–6, 1984. View at: Google Scholar
 J. BakerJarvis and M. D. Janezic, “Analysis of a twoport flanged coaxial holder for shielding effectiveness and dielectric measurements of thin films and thin materials,” IEEE Transactions on Electromagnetic Compatibility, vol. 38, no. 1, pp. 67–70, 1996. View at: Publisher Site  Google Scholar
 A. P. Orfanidis, G. A. Kyriacou, and J. N. Sahalos, “A modematching technique for the study of circular and coaxial waveguide discontinuities based on closedform coupling integrals,” IEEE Transactions on Microwave Theory and Techniques, vol. 48, no. 5, pp. 880–883, 2000. View at: Publisher Site  Google Scholar
 Z. Shen and R. H. MacPhie, “Input admittance of a multilayer insulated monopole Antenna,” IEEE Transactions on Antennas and Propagation, vol. 46, no. 11, pp. 1679–1686, 1998. View at: Publisher Site  Google Scholar
 L. PoloLopez, J. A. RuizCruz, J. R. MontejoGarai, and J. M. Rebollar, “Analysis of waveguide discontinuities with lateral and transverse perfect magnetic wall boundary conditions,” in Proceedings of the 2016 URSI International Symposium on Electromagnetic Theory, EMTS 2016, pp. 811–814, Finland, August 2016. View at: Google Scholar
 J. M. Rebollar and J. A. Encinar, “Field theory analysis of multiportmultidiscontinuity structures: application to shortcircuited Eplane septum,” IEE Proceedings Part H Microwaves, Antennas and Propagation, vol. 135, no. 1, pp. 1–7, 1988. View at: Publisher Site  Google Scholar
 R. E. Collin, Field Theory of Guided Waves, IEEE Press, New York, NY, USA, 1991. View at: MathSciNet
 Z. Shen and R. H. MacPhie, “Scattering by a Thick OffCentered Circular Iris in Circular Waveguide,” IEEE Transactions on Microwave Theory and Techniques, vol. 43, no. 11, pp. 2639–2642, 1995. View at: Publisher Site  Google Scholar
 A. MoránLópez, J. Córcoles, J. A. RuizCruz, J. R. MontejoGarai, and J. M. Rebollar, “Electromagnetic Scattering at the Waveguide Step between Equilateral Triangular Waveguides,” Advances in Mathematical Physics, vol. 2016, Article ID 2974675, pp. 1–16, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 DAK, “DAK  Dielectric Assessment Kit Product Line,” https://speag.swiss/products/dak/dielectricmeasurements. View at: Google Scholar
Copyright
Copyright © 2019 Raúl V. HaroBáez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.