# Modeling of Coaxial Slot Waveguides Using Analytical and Numerical Approaches: Revisited

Authors: Kok Yeow You; Jamaliah Salleh; Mohd Fareq Abd Malek; Zulkifly Abbas; Cheng Ee Meng; Lee Kim Yee
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

Our reviews of analytical methods and numerical methods for coaxial slot waveguides are presented. The theories, background, and physical principles related to frequency-domain electromagnetic equations for coaxial waveguides are reassessed. Comparisons of the accuracies of various types of admittance and impedance equations and numerical simulations are made, and the fringing field at the aperture sensor, which is represented by the lumped capacitance circuit, is evaluated. The accuracy and limitations of the analytical equations are explained in detail. The reasons for the replacement of analytical methods by numerical methods are outlined.


---

## Body

## 1. Introduction

Since the 1940s, coaxial slot waveguides have been used as antennas and as cables for electric power. Recently, this type of waveguide has been used in dielectric measurements and the treatment of cancer. Hence, a theoretical formula was needed to design and model the waveguide. Many theoretical formulas have been developed and derived in the past 70 years. Indeed, in the past, the analytical formulas were an important tool for waveguide design because few numerical methods were available and advance computer equipment and software were not yet available to reduce the costs associated with experimental design. Until the 1960s, many reliable numerical techniques were developed and applied for solving electromagnetic problems, such as the finite difference time-domain (FDTD) method [1], the moment method (MoM) [2], and the finite element method (FEM) [3]. At that time, 2D numerical techniques were used to analyze coaxial waveguides because of their simplicity and the symmetrical shapes they produced.In reality, no completely analytical solutions exist due to the difficulty of deriving the formulas and the many unknown environmental factors that are not included in the formulation. In the author’s view, analytical formulas referred to equations that were based on theoretical concepts and that had no unknown variables. Although the complicated analytical formulas were solved with the aid of numerical methods, such as series expansion and the trapezium rule, the solutions were still obtained by an analytical approach. When a theoretical formula has a matrix of unknown variables and is solved by using a modern numerical routine, such as FEM and MoM, we will define it as a numerical approach. In fact, modern numerical methods can produce results that are close to practical measurements because they divide the problem into many small segments and solve them one by one.However, the development of powerful numerical techniques does not mean that analytical formulas are not still useful in predicting waveguide performance, and, in addition, many of numerical models are based on analytical formulations. For instance, many online broad-band dielectric measurements involve the prediction of the dielectric properties of samples based on measurements of electric signals, but, unfortunately, those numerical techniques are not suitable for solving this type of inverse problem. Therefore, until now, analytical formulas have continued to be important tools for the analysis of coaxial sensors. Furthermore, the application of numerical methods to electromagnetic fields is very challenging because the method strongly affects the operational wavelength and the size of the discrete mesh grid.Most of the origin publications for those analytical formulations do not describe in detail what are the weaknesses and limitations of their equations. However, among are some analytical formulas which are quite reliable and accurate for low-frequency modeling. In this paper, the detailed results of our investigation of the restrictions and limitations of frequency-domain analytical models for coaxial waveguides are presented, and the reasons for using numerical methods to replace analytical methods are discussed. In addition, the strengths and weaknesses of numerical methods are analyzed and discussed.


## 2. Analytical Methods

### 2.1. Open-Ended Coaxial Line

The earliest work [4, 5] concerned the homogeneous case, in which an air-filled coaxial line with an infinite conducting flange is radiated into the infinite half-space. Simultaneously, the variation expressions for normalized admittance aperture, Ỹ, are given by(1)Ỹ=jk2k1ln⁡(b/a)∫0∞dζζ(ζ2-k22)1/2|Jo(ζa)-Jo(ζb)|2=G(0)Yo+jB(0)Yo.
The real part and the imaginary part in (1) are called normalized conductance, G(0)/Yo, and susceptance, B(0)/Yo, respectively [4, 5]:(2)G(0)Yo=k2k1ln⁡(b/a)∫0π/2dθsin⁡θ[Jo(k2asin⁡θ)-Jo(k2bsin⁡θ)]2,(3)B(0)Yo=k2πk1ln⁡(b/a)∫0πdϕ{sin⁡(ϕ2)2Si(k2a|u|)-Si(2k2a|sin⁡(ϕ2)|)-Si(2k2b|sin⁡(ϕ2)|)},
where u2=a2+b2-2abcosϕ and Si are the sine integral or Si(x)=∫0xsint/tdt. Presently, a simplified version of (1) introduced in [6] is the most commonly used equation to calculate the admittance of the open-ended coaxial sensor(4)Ỹ=jk22πk1ln⁡(b/a)∫ab∫ab∫0πexp⁡(-jk2r)rcos⁡ϕ′dϕ′dρ′dρ.The discontinuity at the open end of the coaxial line is commonly interpreted using admittance expressions (1) and (4) due to the naturally capacitive properties that dominate at the open-end surface (z=0). Finally, analytical equations (1) and (4) can be solved using several classical numerical approaches, such as the series expansion method and the Gaussian quadrature routines. However, the analytical equations only are suitable for ideal conditions in which the component electric field, Ez, has been neglected in (1) and (4) [4], as shown in Figure 1(a). In actual situations, in addition to components Eρfield, the component Ez also is present on the surface z=0, as shown in Figure 1(b).Propagated wave from coaxial line (z<0) into free space (z>0) through the ground plane (z=0).
(a)(b)Typically, componentEz is neglected to facilitate solving the problem, since components Ez involve unknown variables in its formulation [7]. So, to get accurate results from the calculations, modern numerical methods have played an important role in determining the values of unknown variables.In the early 1980s, the admittance model was simplified by using an equivalent circuit, in which the circuit consists of two parallel capacitance terms and one conductance term,Go [8]:(5)Ỹin=1Y0[jω(C1+εrC2)+Goω5/2].
The two capacitances, C1 and C2 in (5), are the capacitive properties inside the coaxial line near the aperture probe and the capacitive properties in an external sample under test, respectively. However, the two capacitance terms are not sufficient to represent the equivalent circuit for the aperture probe, which interacts with lossy materials at high-frequency operation. In this study, we propose that inductance, L, and resistance, R, should be taken into account in the circuit model, as shown in Figure 2. In this work, the composition for the C, R, and L elements refer to the impedance before it is converted to function admittance. Thus, the normalized input impedance, Z̃in, and admittance, Ỹin, are given as(6a)Z̃in=Yo[1jω(C1+εrC2)+jω(L1+εrL2)+R],(6b)Ỹin=1Z̃in.The capacitances C1 and C2 in (6a) have the same meaning as in (5). The L1 and L2 are the inductive properties inside the coaxial line near the aperture probe and the inductive properties in an external sample under test, respectively. The R is the resistive properties at the aperture probe. The values of components C1, L1, L2, and R are obtained by optimization between the normalized input impedance, Z̃in, and the finite element simulation results, while the C2 values in (5) and (6a) are calculated from C2=2.38εo(b-a) [9].Comparison of (a) normalized conductance,G(0)/Yo, and (b) normalized susceptance, B(0)/Yo, for water at 25°C by considering size of probe with a=0.65 mm and b=2.05 mm.
(a)(b)The comparison of calculated normalized conductance,G(0)/Yo, and susceptance, B(0)/Yo, using several models and the simulation results (using COMSOL Multiphysics [10], a finite element analysis simulation software) and measurements for water at room temperature are shown in Figures 2(a) and 2(b), respectively. The calculations for admittance in (4) were conducted by using a 5×5× 6-order Gaussian triple integral method, while (2) and (3) were solved by using the series expansion method with 25 series terms. The dispersive properties of the water in simulation and all admittance calculations were obtained from the Cole-Cole model [11] with εs=78.6, ε∞=4.22, τ=8.8ps, and α=0.013. The results measured over the frequency range from 0.3 GHz to 18 GHz were obtained by using a Teflon-filled coaxial sensor with a=0.65 mm and b=2.05 mm.From Figures2(a) and 2(b), the closed form equation (5) with C1 = 0.1 pF, C2 = 0.0293 pF, and Go=9×10-29 is applicable only up to 6 GHz, while (6a) shows that the calculated values are in good agreement with the corresponding simulation results up to 40 GHz. The values for the components C1, L1, L2, and R in (6a) were given as C1=0.1 pF, C2=0.0293 pF, L1=7 pH, L2=0.9 pH, and R=2.8 Ω. Equations (2), (3), and (4) were solved analytically, but with different approaches. The series expansion solution of (2)-(3) agrees with those corresponding admittance results for frequencies only up to 18 GHz and becomes unstable for higher frequencies. This unstable condition may be due to an inaccuracy of the decimal numeral in the series terms. However, the series expansion solutions can be used efficiently to minimize the run time of the program compared to the Gaussian integral method.


### 2.2. Monopole Driving from Coaxial Line

For the extended-conductor case, the normalized input impedance,Zin, can be derived using the induced electromagnetic field (EMF) method [12] as(7)Z̃in=j(0.5)k1k2ln⁡(b/a)sin2(k2h)∫0hsin⁡[k2(h-z)]×[e-jk2R1R1+e-jk2R2R2-2cos⁡(kh)e-jk2r′r′]dz.
While, in [12], listed out the input impedance, Z̃infor monopole driven from coaxial line can be written as(8)Z̃in=j(0.5)k1k2ln⁡(b/a)sin2(k2h)×∫0hsin⁡[k2(h-z)]×[e-jk2R1R1-cos⁡(kh)e-jk2r′r′-zsin⁡(k2h)e-jk2r′(jr′2+1k2r′3)]dz,
where εcis the relative permittivity of lossless material filled in coaxial line (at region z≤0). Equations (7) and (8) can be solved easily by Simpson’s rules. For extended conductor cases, the energy is released through the extended conductor from the end of the coaxial line, thus the inductive properties are dominant at driving point (z=0). Hence, the analysis and the prediction of the performance of an extended conductor waveguide are often referred as input impedance, Zin=Rin+jXin, rather than input admittance, Yin=Gin+jBin.Figure3 shows the calculated normalized resistance, Rin/Zo, and the normalized reactance, Xin/Zo, for air at room temperature using (7) and (8). It was found that (7) and (8) give similar resistance values, Rin/Zo, in the frequency range of 300 kHz to 20 GHz. Meanwhile, both equations give les consistent values of reactance, and this may be due to the fact that the field distribution near the source point was different between the two models. The impedance integral equations are important in analytical calculations and for numerical solutions. Numerical solutions, such as MoM, cannot be used if the problems have no established analytical integral equation. Moreover, the accuracy of the MoM method is dependent on the rigor of the analytical integral equation. For instance, (1) and (8) have been used as the governing equations in the MoM method [7, 13, 14]. However, in the MoM solution, the current, Iz, is an unknown parameter that must be solved by a numerical routine. Instead, in analytical calculations, the current, Iz, is usually replaced by the sinusoidal current model [12, 13, 15] in the integral equation before it is solved.Comparison of calculated normalized resistance,Rin/Zo, and normalized reactance, Xin/Zo, for air at 25°C  using (7) and (8) (a=0.65 mm, b=2.05 mm, and h=14.32 mm).
(a)(b)The distribution current,Iz, substituted in (7) and (8) was assumed to have a convenient sinusoidal form, such as [12, 13, 15]:(9)Iz=I(0)sin⁡k2(h-z)sin⁡(k2h),
where I(0) is the amplitude of the driving point current at z=0, and it was canceled out in the derivations in (7) and (8). Until now, no completely actual distribution current, Iz, formulation has existed for an arbitrarily sized extended conductor. The current calculated using (9) is clearly the real value, but, in an actual case, the current is a complex value due to the azimuthal current, Iϕ, contributed by the size of the radius of conductor a, as shown in Figure 4(b). This means that (9) is valid only for thin extended conductors. In addition, it cannot be guaranteed that the current distribution will always be sinusoidal in all cases of the extended conductor.(a) Ideal analytical line form current; (b) Actual situations of current distribution on a cylindrical extended conductor.
(a)(b)Figure5 shows the variational of distribution current, Iz along the length, z, of monopole. The unknown complex current, Iz, is determined by using point matching MoM with Pocklington’s integral equation and frill-generator source [14, 16, 17]. For cylindrical monopole structure, the magnetic field, H⃗, around the surface monopole which is generated by a current, Iz, is in an azimuthal, ϕ, direction, and the relationship is given as Hϕ=Iz/2πaϕ̂. This means that the imaginary part of current, Iz, is referring to the current flow in z-direction (perpendicular to the ϕ direction). While, the real part of Iz is represented in the azimuthal current. From Figure 5, it is clear that, at low frequencies (f<3 GHz), the value of current, Iz, is approaching the real number. When the monopole operated at higher frequencies, the distribution current, Iz, along the length of monopole become complex number.Simulated distribution current,Iz, on the 1.436 cm of monopole driving from coaxial line with a=0.65 mm and b=2.05 mm.
(a)(b)(c)(d)(e)(f)In addition to the uncertainty of the distribution current, (7) and (8) also have not considered the fringing fields contributed by the electric field, Eρ, near the driving point (z=0) and the end of the terminate conductor (z=h). For a short conductor, (h<10a), driven from a coaxial line, the input impedance is affected significantly by the fringing field at the end of the coaxial line, especially for high-frequency operation [18]. This means that the analytical impedance formulas of (7) and (8) are accurate only for the low-frequency condition. For the reasons stated, the full wave (Hϕ, Eρ, and Ez) analysis of the waveguide must be conducted by using a numerical method. In practice, the effects of fringing fields at the driving point are always empirically corrected by capacitance, Cf, element circuits.In addition to integration models, the input impedance of the monopole also can be represented by using the lumped-element model [19] and expressed as(10)Z=j(ωL1-1ωC1)︸Lowfrequencies+1jωC2+(1/(jωL2+R1))︸Highfrequencies+1jωC3+[1/(jωL3+(50jωL4/(jωL4+50)))]︸.DrivingfringingfieldsFrom Figures7(a) and 7(b), the input impedance show that only two peaks of curve line occurred from 300 kHz to 20 GHz, the three terms of (10) are sufficient to model the impedance properties cover the frequency range. The first term in (10) describes the characteristics of the impedance at low frequencies, while the second term contributes to the modeling of impedance for high frequencies (>10 GHz). The third term is used to model the equivalent circuit, which is near the driving point for the monopole. This means that if there are three peaks of impedance curve line (See Figure 14) over 300 kHz to 20 GHz, up to four terms of lumped-element expression are required. The equivalent circuit and all optimized values of resistance, R, capacitance, C, and inductance, L, elements in (10) for air are shown in Figure 6(b), and the values are accurate up to 20 GHz. The RLC values in (10) are obtained by optimizing (10) to the measurement results. At low frequencies, the capacitance term, C1, in (10) has played an important role in the modeling, since the input impedance at driving point (z=0) shows the nature of capacitive properties (Xin/Zo in negative (−) sign) at a frequency of less than 4 GHz, as shown in Figure 7(b).Figure 6
Dimensions of monopole and simulated electric fields contour; (b) equivalent circuit for input impedance of monopole atz=0.Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for air at 25°C by considering size of monopole driving from coaxial line with a=0.65 mm, b=2.05 mm, and h=14.32 mm.
(a)(b)For instance, a monopole withh=1.436 cm driven from the coaxial line, as shown in Figure 6(a), was tested in this work. The comparison of calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, using several models, moment method [14, 16, 17], simulation results (using COMSOL Multiphysics [10] and CST Microwave Studio [20]), and measurements for air at room temperature are shown in Figures 7(a) and 7(b), respectively. Different from the coaxial probe, for a monopole, the fringing field capacitance, Cf, elements first refer to the admittance before being converted to the impedance function, as follows:(11a)Ỹ=1Equation(7)+jωCf,(11b)Z̃corrected=1Ỹ.In this work, the value ofCf was 5.5 pF. Figure 7 shows that a 1.436 cm monopole has been matched with a standard, 50 Ω cable at a frequency of 4.6 GHz and 14.4 GHz.


### 2.3. Coupling Monopole Driving from Coaxial Line

When two different lengths of monopoles are placed closed to each other, as shown in Figure8(a), the combination of two electromagnetic fields occurs. Now, the input impedance of the monopole is required to consider the self-radiation and mutual-radiation effects.(a) Arbitrary lengths of two parallel, mutually coupled monopoles; (b) actual configuration of coupling monopole with one excitation port.
(a)(b)The mutual normalized impedance,Z̃12, is well expressed as [21]:(12)Z̃12=j(0.5)k1k2ln⁡(b/a)sin⁡(k2h1)sin⁡(k2h2)×∫0h2sin⁡[k2(h2-z)]×[e-jk2R1′R1′+e-jk2R2′R2′-2cos⁡(k2h1)e-jk2R′R′]dz,
where R1′=D2+(z-h1)2, R2′=D2+(z+h1)2, and R′=D2+z2. Finally, the input impedance, Z̃in_port-1, of monopole 1 can be calculated as [15, 21, 22]:
(13)Z̃in_port-1=Z̃11-Z̃122Z̃22,
where the term Z̃11=(7) is assumed to be equal to Z̃22, and Z̃21 is equal to Z̃12. Similarly, the input impedance, Zin_port-1, of the coupling monopole also can be represented by using the lumped-element model and expressed as(14)Zin_port-1=j(ωL1-1ωC1)︸Lowfrequencies+1jωC2+(1/(jωL2+R1))︸Highfrequencies+1jωC3+[1/(jωL3+(50jωL4/(jωL4+50)))]︸.Drivingfringingfields
The equivalent circuit for (14) and its resistance, R, capacitance, C, and inductance, L, elements for air, which are accurate up to 20 GHz, are shown in Figure 9. The RLC values in equation (14) are obtained by optimizing (14) to the measurement results. From the comparison of calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, using several models, the CST simulation results and measurements for air at room temperature are shown in Figures 10(a) and 10(b), respectively.Figure 9
Equivalent circuit for input impedance of coupling monopole at port-1 (excitation port).Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, of the coupling monopole at excitation port-1 (a=0.65 mm, b=2.05 mm, h1=h2=14.32 mm, and D=10 mm).
(a)(b)The equation for the coupling input impedance,Zin_port-2, at port-2 can be expressed in terms of input impedance, Zin_port-1, at excitation port-1 as(15a)T=(1-Equation(14)-ZoEquation(14)+Zo)e-(α+jβ)d′,(15b)Zin_port-2=2Zo(1-TT),where T is the transmission coefficient at parasitic port-2. The exponential term, e-(α+jβ)d′, is the transmission factor of the transmitted waves from port-1 to port-2. The symbols α and β are the attenuation constant and phase constant, respectively. The equation d′=h1+h2+D+Δ gives the length of the transmission wave from port-1 to port-2. The normalized input impedances, Zin_port-2=Rin+jXin, calculated by using (15b), are plotted and compared with measurement and CST simulation results, as shown in Figure 11. The normalized resistance, Rin/Zo, and reactance, Xin/Zo, in Figures 11(a) and 11(b) are calculated with α=3.18×10-5f neper/m,β=2πf/c rad/m, d′=0.044m and Zo=50Ω.Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, of the coupling monopole at parasitic port-2 (a=0.65 mm, b=2.05 mm, h1=h2=14.32 mm, and D=10 mm).
(a)(b)


### 2.4. Coated Conductor Driving from Coaxial Line

A coated monopole is a bare monopole enclosed by a thin cylindrical low-loss dielectric material with an outer radius,b, as shown in Figure 12.(a) Two-dimensional and (b) three-dimensional configurations of a coated monopole.
(a)(b)In general, the transmission line formulas are adapted easily to the analysis of a coated antenna. The input impedance,Z̃in, of the coated monopole at the drive point (z=0) [23–25] is expressed as(16)Z̃in=-jZcZocot(kLh′),
where Zo and Zcare the characteristic impedance in the coaxial line and the characteristic impedance for the coated monopole transmission line, respectively, which can be given as:
(17)Zc=γkL2πkc[ln⁡(ba)+(kck2)2Ho(2)(k2b)k2bH1(2)(k2b)].
The complex propagation constant, kL, is expressed as:(18)kL=kck2[Ho(2)(k2b)+k2bln⁡(b/a)H1(2)(k2b)kc2Ho(2)(k2b)+k23bln⁡(b/a)H1(2)(k2b)]1/2.
The transmissionline model requires that the propagation constant, k2, of the external medium be much larger than the propagation constant, kc, of the insulation [23–25]. Thus, the coated monopole is almost always used in an ambient medium with higher dielectric properties than its insulation, such as soil, seawater, or biological tissue. Moreover, coating material, such as Teflon, functions as a hygienic protector when the monopole is immersed in a biological sample. However, for the case in Figure 12, the fringing effects at the top end of the monopole are required to consider when its input impedance is calculated using (16). In this study, the effective length, h’, of the monopole was used to correct the fringing effects, and it is given as(19)h′=h+κf,
where h and f are the actual length of the coated monopole and the operational frequency, respectively. The symbol κ is a coefficient value which depends on the dimensions of the monopole. The fringing effects are assumed to be inversely proportional to the square root of frequency. Similarly, the equivalent circuit can be used to represent the input impedance properties of the coated monopole, as shown in Figure 12, and the corresponding formulations are expressed as (20):(20)Zin=j(ωLT-1ωCT)+∑n=1m[1jωCn+(1/(jωLn+Rn))]+1jωCB+[1/(jωLB+(50jωLB′/(jωLB′+50)))].
Figures 13 and 14 show the calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, respectively, using several models and the COMSOL simulation results and measurements for air and water at room temperature (25±1°C) over the frequency range from 300 kHz to 20 GHz. The results were measured using a Teflon-coated monopole driven from coaxial line with a=0.65 mm, c=2.05 mm, and h=13.92 mm. The values for component resistor, R, reactor, L, and capacitor, C, in (20) for air and water are listed in Table 1, which provides the accurately calculated impedance up to 20 GHz. The RLC values in (20) are obtained by optimizing (20) to the measurement results. As expected, compared to the water sample, the calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, for air using (16) do not agree well with the measurements and simulation results. The propagation constant, γ, at the boundary between the coated medium and the external medium is γ=±koεr-εc, thus the tendency properties for the air case (k2<kc) are different those for the water case (k2>kc). When k2<kc, the waves decay in the normal direction to the interface between the coated dielectric and external medium due to the fact that the corresponding propagation constant,γ, is imaginary.Table 1
RCL component values in (20) for air and water samples.
RCL componentsAirWaterCB0.33 pF0.70 pFCT0.40 pF2.50 pFC10.54 pF0.65 pFC2—0.45 pFLB1 nH1 nHLB′1.1 nH2.4 nHLT0.5 nH0.35 nHL10.2 nH0.32 nHL2—0.16 nHR12.9 Ω5.56 ΩR2—4.0 ΩComparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for air at 25°C by considering size of probe with a=0.65 mm, c=2.05 mm, and h=13.92 mm.
(a)(b)Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for water at 25°C by considering size of probe with a=0.65mm, c=2.05 mm, and h=13.92 mm.
(a)(b)


## 2.1. Open-Ended Coaxial Line

The earliest work [4, 5] concerned the homogeneous case, in which an air-filled coaxial line with an infinite conducting flange is radiated into the infinite half-space. Simultaneously, the variation expressions for normalized admittance aperture, Ỹ, are given by(1)Ỹ=jk2k1ln⁡(b/a)∫0∞dζζ(ζ2-k22)1/2|Jo(ζa)-Jo(ζb)|2=G(0)Yo+jB(0)Yo.
The real part and the imaginary part in (1) are called normalized conductance, G(0)/Yo, and susceptance, B(0)/Yo, respectively [4, 5]:(2)G(0)Yo=k2k1ln⁡(b/a)∫0π/2dθsin⁡θ[Jo(k2asin⁡θ)-Jo(k2bsin⁡θ)]2,(3)B(0)Yo=k2πk1ln⁡(b/a)∫0πdϕ{sin⁡(ϕ2)2Si(k2a|u|)-Si(2k2a|sin⁡(ϕ2)|)-Si(2k2b|sin⁡(ϕ2)|)},
where u2=a2+b2-2abcosϕ and Si are the sine integral or Si(x)=∫0xsint/tdt. Presently, a simplified version of (1) introduced in [6] is the most commonly used equation to calculate the admittance of the open-ended coaxial sensor(4)Ỹ=jk22πk1ln⁡(b/a)∫ab∫ab∫0πexp⁡(-jk2r)rcos⁡ϕ′dϕ′dρ′dρ.The discontinuity at the open end of the coaxial line is commonly interpreted using admittance expressions (1) and (4) due to the naturally capacitive properties that dominate at the open-end surface (z=0). Finally, analytical equations (1) and (4) can be solved using several classical numerical approaches, such as the series expansion method and the Gaussian quadrature routines. However, the analytical equations only are suitable for ideal conditions in which the component electric field, Ez, has been neglected in (1) and (4) [4], as shown in Figure 1(a). In actual situations, in addition to components Eρfield, the component Ez also is present on the surface z=0, as shown in Figure 1(b).Propagated wave from coaxial line (z<0) into free space (z>0) through the ground plane (z=0).
(a)(b)Typically, componentEz is neglected to facilitate solving the problem, since components Ez involve unknown variables in its formulation [7]. So, to get accurate results from the calculations, modern numerical methods have played an important role in determining the values of unknown variables.In the early 1980s, the admittance model was simplified by using an equivalent circuit, in which the circuit consists of two parallel capacitance terms and one conductance term,Go [8]:(5)Ỹin=1Y0[jω(C1+εrC2)+Goω5/2].
The two capacitances, C1 and C2 in (5), are the capacitive properties inside the coaxial line near the aperture probe and the capacitive properties in an external sample under test, respectively. However, the two capacitance terms are not sufficient to represent the equivalent circuit for the aperture probe, which interacts with lossy materials at high-frequency operation. In this study, we propose that inductance, L, and resistance, R, should be taken into account in the circuit model, as shown in Figure 2. In this work, the composition for the C, R, and L elements refer to the impedance before it is converted to function admittance. Thus, the normalized input impedance, Z̃in, and admittance, Ỹin, are given as(6a)Z̃in=Yo[1jω(C1+εrC2)+jω(L1+εrL2)+R],(6b)Ỹin=1Z̃in.The capacitances C1 and C2 in (6a) have the same meaning as in (5). The L1 and L2 are the inductive properties inside the coaxial line near the aperture probe and the inductive properties in an external sample under test, respectively. The R is the resistive properties at the aperture probe. The values of components C1, L1, L2, and R are obtained by optimization between the normalized input impedance, Z̃in, and the finite element simulation results, while the C2 values in (5) and (6a) are calculated from C2=2.38εo(b-a) [9].Comparison of (a) normalized conductance,G(0)/Yo, and (b) normalized susceptance, B(0)/Yo, for water at 25°C by considering size of probe with a=0.65 mm and b=2.05 mm.
(a)(b)The comparison of calculated normalized conductance,G(0)/Yo, and susceptance, B(0)/Yo, using several models and the simulation results (using COMSOL Multiphysics [10], a finite element analysis simulation software) and measurements for water at room temperature are shown in Figures 2(a) and 2(b), respectively. The calculations for admittance in (4) were conducted by using a 5×5× 6-order Gaussian triple integral method, while (2) and (3) were solved by using the series expansion method with 25 series terms. The dispersive properties of the water in simulation and all admittance calculations were obtained from the Cole-Cole model [11] with εs=78.6, ε∞=4.22, τ=8.8ps, and α=0.013. The results measured over the frequency range from 0.3 GHz to 18 GHz were obtained by using a Teflon-filled coaxial sensor with a=0.65 mm and b=2.05 mm.From Figures2(a) and 2(b), the closed form equation (5) with C1 = 0.1 pF, C2 = 0.0293 pF, and Go=9×10-29 is applicable only up to 6 GHz, while (6a) shows that the calculated values are in good agreement with the corresponding simulation results up to 40 GHz. The values for the components C1, L1, L2, and R in (6a) were given as C1=0.1 pF, C2=0.0293 pF, L1=7 pH, L2=0.9 pH, and R=2.8 Ω. Equations (2), (3), and (4) were solved analytically, but with different approaches. The series expansion solution of (2)-(3) agrees with those corresponding admittance results for frequencies only up to 18 GHz and becomes unstable for higher frequencies. This unstable condition may be due to an inaccuracy of the decimal numeral in the series terms. However, the series expansion solutions can be used efficiently to minimize the run time of the program compared to the Gaussian integral method.


## 2.2. Monopole Driving from Coaxial Line

For the extended-conductor case, the normalized input impedance,Zin, can be derived using the induced electromagnetic field (EMF) method [12] as(7)Z̃in=j(0.5)k1k2ln⁡(b/a)sin2(k2h)∫0hsin⁡[k2(h-z)]×[e-jk2R1R1+e-jk2R2R2-2cos⁡(kh)e-jk2r′r′]dz.
While, in [12], listed out the input impedance, Z̃infor monopole driven from coaxial line can be written as(8)Z̃in=j(0.5)k1k2ln⁡(b/a)sin2(k2h)×∫0hsin⁡[k2(h-z)]×[e-jk2R1R1-cos⁡(kh)e-jk2r′r′-zsin⁡(k2h)e-jk2r′(jr′2+1k2r′3)]dz,
where εcis the relative permittivity of lossless material filled in coaxial line (at region z≤0). Equations (7) and (8) can be solved easily by Simpson’s rules. For extended conductor cases, the energy is released through the extended conductor from the end of the coaxial line, thus the inductive properties are dominant at driving point (z=0). Hence, the analysis and the prediction of the performance of an extended conductor waveguide are often referred as input impedance, Zin=Rin+jXin, rather than input admittance, Yin=Gin+jBin.Figure3 shows the calculated normalized resistance, Rin/Zo, and the normalized reactance, Xin/Zo, for air at room temperature using (7) and (8). It was found that (7) and (8) give similar resistance values, Rin/Zo, in the frequency range of 300 kHz to 20 GHz. Meanwhile, both equations give les consistent values of reactance, and this may be due to the fact that the field distribution near the source point was different between the two models. The impedance integral equations are important in analytical calculations and for numerical solutions. Numerical solutions, such as MoM, cannot be used if the problems have no established analytical integral equation. Moreover, the accuracy of the MoM method is dependent on the rigor of the analytical integral equation. For instance, (1) and (8) have been used as the governing equations in the MoM method [7, 13, 14]. However, in the MoM solution, the current, Iz, is an unknown parameter that must be solved by a numerical routine. Instead, in analytical calculations, the current, Iz, is usually replaced by the sinusoidal current model [12, 13, 15] in the integral equation before it is solved.Comparison of calculated normalized resistance,Rin/Zo, and normalized reactance, Xin/Zo, for air at 25°C  using (7) and (8) (a=0.65 mm, b=2.05 mm, and h=14.32 mm).
(a)(b)The distribution current,Iz, substituted in (7) and (8) was assumed to have a convenient sinusoidal form, such as [12, 13, 15]:(9)Iz=I(0)sin⁡k2(h-z)sin⁡(k2h),
where I(0) is the amplitude of the driving point current at z=0, and it was canceled out in the derivations in (7) and (8). Until now, no completely actual distribution current, Iz, formulation has existed for an arbitrarily sized extended conductor. The current calculated using (9) is clearly the real value, but, in an actual case, the current is a complex value due to the azimuthal current, Iϕ, contributed by the size of the radius of conductor a, as shown in Figure 4(b). This means that (9) is valid only for thin extended conductors. In addition, it cannot be guaranteed that the current distribution will always be sinusoidal in all cases of the extended conductor.(a) Ideal analytical line form current; (b) Actual situations of current distribution on a cylindrical extended conductor.
(a)(b)Figure5 shows the variational of distribution current, Iz along the length, z, of monopole. The unknown complex current, Iz, is determined by using point matching MoM with Pocklington’s integral equation and frill-generator source [14, 16, 17]. For cylindrical monopole structure, the magnetic field, H⃗, around the surface monopole which is generated by a current, Iz, is in an azimuthal, ϕ, direction, and the relationship is given as Hϕ=Iz/2πaϕ̂. This means that the imaginary part of current, Iz, is referring to the current flow in z-direction (perpendicular to the ϕ direction). While, the real part of Iz is represented in the azimuthal current. From Figure 5, it is clear that, at low frequencies (f<3 GHz), the value of current, Iz, is approaching the real number. When the monopole operated at higher frequencies, the distribution current, Iz, along the length of monopole become complex number.Simulated distribution current,Iz, on the 1.436 cm of monopole driving from coaxial line with a=0.65 mm and b=2.05 mm.
(a)(b)(c)(d)(e)(f)In addition to the uncertainty of the distribution current, (7) and (8) also have not considered the fringing fields contributed by the electric field, Eρ, near the driving point (z=0) and the end of the terminate conductor (z=h). For a short conductor, (h<10a), driven from a coaxial line, the input impedance is affected significantly by the fringing field at the end of the coaxial line, especially for high-frequency operation [18]. This means that the analytical impedance formulas of (7) and (8) are accurate only for the low-frequency condition. For the reasons stated, the full wave (Hϕ, Eρ, and Ez) analysis of the waveguide must be conducted by using a numerical method. In practice, the effects of fringing fields at the driving point are always empirically corrected by capacitance, Cf, element circuits.In addition to integration models, the input impedance of the monopole also can be represented by using the lumped-element model [19] and expressed as(10)Z=j(ωL1-1ωC1)︸Lowfrequencies+1jωC2+(1/(jωL2+R1))︸Highfrequencies+1jωC3+[1/(jωL3+(50jωL4/(jωL4+50)))]︸.DrivingfringingfieldsFrom Figures7(a) and 7(b), the input impedance show that only two peaks of curve line occurred from 300 kHz to 20 GHz, the three terms of (10) are sufficient to model the impedance properties cover the frequency range. The first term in (10) describes the characteristics of the impedance at low frequencies, while the second term contributes to the modeling of impedance for high frequencies (>10 GHz). The third term is used to model the equivalent circuit, which is near the driving point for the monopole. This means that if there are three peaks of impedance curve line (See Figure 14) over 300 kHz to 20 GHz, up to four terms of lumped-element expression are required. The equivalent circuit and all optimized values of resistance, R, capacitance, C, and inductance, L, elements in (10) for air are shown in Figure 6(b), and the values are accurate up to 20 GHz. The RLC values in (10) are obtained by optimizing (10) to the measurement results. At low frequencies, the capacitance term, C1, in (10) has played an important role in the modeling, since the input impedance at driving point (z=0) shows the nature of capacitive properties (Xin/Zo in negative (−) sign) at a frequency of less than 4 GHz, as shown in Figure 7(b).Figure 6
Dimensions of monopole and simulated electric fields contour; (b) equivalent circuit for input impedance of monopole atz=0.Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for air at 25°C by considering size of monopole driving from coaxial line with a=0.65 mm, b=2.05 mm, and h=14.32 mm.
(a)(b)For instance, a monopole withh=1.436 cm driven from the coaxial line, as shown in Figure 6(a), was tested in this work. The comparison of calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, using several models, moment method [14, 16, 17], simulation results (using COMSOL Multiphysics [10] and CST Microwave Studio [20]), and measurements for air at room temperature are shown in Figures 7(a) and 7(b), respectively. Different from the coaxial probe, for a monopole, the fringing field capacitance, Cf, elements first refer to the admittance before being converted to the impedance function, as follows:(11a)Ỹ=1Equation(7)+jωCf,(11b)Z̃corrected=1Ỹ.In this work, the value ofCf was 5.5 pF. Figure 7 shows that a 1.436 cm monopole has been matched with a standard, 50 Ω cable at a frequency of 4.6 GHz and 14.4 GHz.


## 2.3. Coupling Monopole Driving from Coaxial Line

When two different lengths of monopoles are placed closed to each other, as shown in Figure8(a), the combination of two electromagnetic fields occurs. Now, the input impedance of the monopole is required to consider the self-radiation and mutual-radiation effects.(a) Arbitrary lengths of two parallel, mutually coupled monopoles; (b) actual configuration of coupling monopole with one excitation port.
(a)(b)The mutual normalized impedance,Z̃12, is well expressed as [21]:(12)Z̃12=j(0.5)k1k2ln⁡(b/a)sin⁡(k2h1)sin⁡(k2h2)×∫0h2sin⁡[k2(h2-z)]×[e-jk2R1′R1′+e-jk2R2′R2′-2cos⁡(k2h1)e-jk2R′R′]dz,
where R1′=D2+(z-h1)2, R2′=D2+(z+h1)2, and R′=D2+z2. Finally, the input impedance, Z̃in_port-1, of monopole 1 can be calculated as [15, 21, 22]:
(13)Z̃in_port-1=Z̃11-Z̃122Z̃22,
where the term Z̃11=(7) is assumed to be equal to Z̃22, and Z̃21 is equal to Z̃12. Similarly, the input impedance, Zin_port-1, of the coupling monopole also can be represented by using the lumped-element model and expressed as(14)Zin_port-1=j(ωL1-1ωC1)︸Lowfrequencies+1jωC2+(1/(jωL2+R1))︸Highfrequencies+1jωC3+[1/(jωL3+(50jωL4/(jωL4+50)))]︸.Drivingfringingfields
The equivalent circuit for (14) and its resistance, R, capacitance, C, and inductance, L, elements for air, which are accurate up to 20 GHz, are shown in Figure 9. The RLC values in equation (14) are obtained by optimizing (14) to the measurement results. From the comparison of calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, using several models, the CST simulation results and measurements for air at room temperature are shown in Figures 10(a) and 10(b), respectively.Figure 9
Equivalent circuit for input impedance of coupling monopole at port-1 (excitation port).Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, of the coupling monopole at excitation port-1 (a=0.65 mm, b=2.05 mm, h1=h2=14.32 mm, and D=10 mm).
(a)(b)The equation for the coupling input impedance,Zin_port-2, at port-2 can be expressed in terms of input impedance, Zin_port-1, at excitation port-1 as(15a)T=(1-Equation(14)-ZoEquation(14)+Zo)e-(α+jβ)d′,(15b)Zin_port-2=2Zo(1-TT),where T is the transmission coefficient at parasitic port-2. The exponential term, e-(α+jβ)d′, is the transmission factor of the transmitted waves from port-1 to port-2. The symbols α and β are the attenuation constant and phase constant, respectively. The equation d′=h1+h2+D+Δ gives the length of the transmission wave from port-1 to port-2. The normalized input impedances, Zin_port-2=Rin+jXin, calculated by using (15b), are plotted and compared with measurement and CST simulation results, as shown in Figure 11. The normalized resistance, Rin/Zo, and reactance, Xin/Zo, in Figures 11(a) and 11(b) are calculated with α=3.18×10-5f neper/m,β=2πf/c rad/m, d′=0.044m and Zo=50Ω.Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, of the coupling monopole at parasitic port-2 (a=0.65 mm, b=2.05 mm, h1=h2=14.32 mm, and D=10 mm).
(a)(b)


## 2.4. Coated Conductor Driving from Coaxial Line

A coated monopole is a bare monopole enclosed by a thin cylindrical low-loss dielectric material with an outer radius,b, as shown in Figure 12.(a) Two-dimensional and (b) three-dimensional configurations of a coated monopole.
(a)(b)In general, the transmission line formulas are adapted easily to the analysis of a coated antenna. The input impedance,Z̃in, of the coated monopole at the drive point (z=0) [23–25] is expressed as(16)Z̃in=-jZcZocot(kLh′),
where Zo and Zcare the characteristic impedance in the coaxial line and the characteristic impedance for the coated monopole transmission line, respectively, which can be given as:
(17)Zc=γkL2πkc[ln⁡(ba)+(kck2)2Ho(2)(k2b)k2bH1(2)(k2b)].
The complex propagation constant, kL, is expressed as:(18)kL=kck2[Ho(2)(k2b)+k2bln⁡(b/a)H1(2)(k2b)kc2Ho(2)(k2b)+k23bln⁡(b/a)H1(2)(k2b)]1/2.
The transmissionline model requires that the propagation constant, k2, of the external medium be much larger than the propagation constant, kc, of the insulation [23–25]. Thus, the coated monopole is almost always used in an ambient medium with higher dielectric properties than its insulation, such as soil, seawater, or biological tissue. Moreover, coating material, such as Teflon, functions as a hygienic protector when the monopole is immersed in a biological sample. However, for the case in Figure 12, the fringing effects at the top end of the monopole are required to consider when its input impedance is calculated using (16). In this study, the effective length, h’, of the monopole was used to correct the fringing effects, and it is given as(19)h′=h+κf,
where h and f are the actual length of the coated monopole and the operational frequency, respectively. The symbol κ is a coefficient value which depends on the dimensions of the monopole. The fringing effects are assumed to be inversely proportional to the square root of frequency. Similarly, the equivalent circuit can be used to represent the input impedance properties of the coated monopole, as shown in Figure 12, and the corresponding formulations are expressed as (20):(20)Zin=j(ωLT-1ωCT)+∑n=1m[1jωCn+(1/(jωLn+Rn))]+1jωCB+[1/(jωLB+(50jωLB′/(jωLB′+50)))].
Figures 13 and 14 show the calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, respectively, using several models and the COMSOL simulation results and measurements for air and water at room temperature (25±1°C) over the frequency range from 300 kHz to 20 GHz. The results were measured using a Teflon-coated monopole driven from coaxial line with a=0.65 mm, c=2.05 mm, and h=13.92 mm. The values for component resistor, R, reactor, L, and capacitor, C, in (20) for air and water are listed in Table 1, which provides the accurately calculated impedance up to 20 GHz. The RLC values in (20) are obtained by optimizing (20) to the measurement results. As expected, compared to the water sample, the calculated normalized resistance, Rin/Zo, and reactance, Xin/Zo, for air using (16) do not agree well with the measurements and simulation results. The propagation constant, γ, at the boundary between the coated medium and the external medium is γ=±koεr-εc, thus the tendency properties for the air case (k2<kc) are different those for the water case (k2>kc). When k2<kc, the waves decay in the normal direction to the interface between the coated dielectric and external medium due to the fact that the corresponding propagation constant,γ, is imaginary.Table 1
RCL component values in (20) for air and water samples.
RCL componentsAirWaterCB0.33 pF0.70 pFCT0.40 pF2.50 pFC10.54 pF0.65 pFC2—0.45 pFLB1 nH1 nHLB′1.1 nH2.4 nHLT0.5 nH0.35 nHL10.2 nH0.32 nHL2—0.16 nHR12.9 Ω5.56 ΩR2—4.0 ΩComparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for air at 25°C by considering size of probe with a=0.65 mm, c=2.05 mm, and h=13.92 mm.
(a)(b)Comparison of (a) normalized resistance,Rin/Zo, and (b) normalized reactance, Xin/Zo, for water at 25°C by considering size of probe with a=0.65mm, c=2.05 mm, and h=13.92 mm.
(a)(b)


## 3. Conclusions

In this work, the coaxial waveguide was used as an example of the problem of linking electromagnetic theory with practical modeling, since the coaxial slot waveguides have been used as antennas over the past of 70 years. Recently, many scientific applications have involved this kind of waveguide. For instance, an open-ended coaxial probe was applied as a dielectric probe to measure the dielectric properties of the material being tested. In addition, the dielectric-coated monopole was used for hyperthermia treatment, and the coupler monopole was designed to be an array antenna. Hence, many models have been developed for the coaxial waveguide, and those models have been modified for use in modeling other devices. In particular, the sinusoidal current model is also used for planar waveguides. In this study, the accuracy of frequency-domain analytical models was tested by acquiring measurements and numerical simulation results. We found that the semiempirical equivalent circuit modeling worked successfully, covered a wide frequency range, and was very useful in the design of circuits that matched the waveguides. Although the analytical models are less accurate compare to numerical method, it provides significant rapid and economize computation. Implicitly, the analytical models still retain the academy valuable, especially for who has preliminary study of the antenna modeling.



---

*Source: 701532-2012-04-26.xml*