Open access peer-reviewed chapter - ONLINE FIRST

Waveguide Port Approach in EM Simulation of Microwave Antennas

Written By

Faik Bogdanov, Irina Chochia, Lily Svanidze and Roman Jobava

Submitted: January 29th, 2022 Reviewed: February 1st, 2022 Published: March 24th, 2022

DOI: 10.5772/intechopen.102996

IntechOpen
Microwave Technologies Edited by Ahmed Kishk

From the Edited Volume

Microwave Technologies [Working Title]

Dr. Ahmed Kishk and Dr. Kim Ho Yeap

Chapter metrics overview

21 Chapter Downloads

View Full Metrics

Abstract

This chapter generalizes a recently proposed MoM-based approach to waveguide port excitation (WPE) problems on arbitrary conducting and composite geometries. This approach combines the canonical aperture coupling approach with the EFIE-PMCHWT formulation for composite structures. Each WPE problem in this approach is divided into equivalent sub-problems for internal and external regions, which are solved using the MoM. Internal WPE problems are solved using waveguide modal expansion in the port plane, while external problems are solved using the equivalence principle to reduce these problems to the systems of algebraic equations for unknown electric and magnetic currents. The developed approach is validated on radiation and coupling problems for coaxial ports by comparing simulated results with those obtained by other approaches and measurements. An excellent agreement between the simulated and measured results is demonstrated. Finally, this approach is applied to practical EMC problems for microwave antennas fed by coaxial ports.

Keywords

  • coupling problem
  • coaxial port
  • equivalence principle
  • method of moments (MoM)
  • waveguide modal expansion
  • waveguide port

1. Introduction

The excitation problem is of increasing importance at microwave frequencies [1]. Microwave antennas and other microwave devices are often fed by waveguides (rectangular, circular, and coaxial) or transmission lines (such as microstrips). In general, these devices are composite structures consisting of both conductive and dielectric elements. Therefore, the appropriate modeling of the waveguide excitation of such structures is of great interest. Such modeling in numerical methods is usually done by truncating the feed waveguide to create a waveguide port and formulate suitable boundary conditions (BC) imposed on the port. Such BCs should be able to launch an incident wave into the waveguide and absorb the reflected (in active mode) or received (in passive mode) wave without spurious reflections [2].

To date, most approaches to solving the waveguide port excitation (WPE) problem are based on volume discretization methods, such as the finite-element method (FEM) [1, 2], finite difference time domain (FDTD) [3, 4, 5], discontinues Galerkin time-domain (DGTD) [6], contour integral method (CIM) [7], etc. Most of these works use various modal absorbing boundary conditions (MABC) [4, 5], developed for time-domain methods as termination conditions imposed on the port.

At present, many electromagnetic (EM) problems are solved using surface integral equations (SIE) together with the method of moments (MoM) [8]. Within the framework of SIE, the WPE problem was first formulated as an aperture coupling problem for a conducting geometry, and the MoM solution for magnetic currents was obtained in the presence of a short-circuited conductive sheet [9]. This approach was then modified using the pseudo-image method for magnetic currents in addition to electric currents [10]. Further, MoM was applied to the waveguide port problems [11] and antenna radiation problems with aperture port excitation [12]. However, until recently, a MoM-based solution to the WPE problem for arbitrary geometries has been poorly represented in the literature. In our recent works [13, 14, 15], such a solution was obtained for radiation and coupling problems for various types of geometries.

This chapter generalizes the recently proposed MoM-based approach to WPE problems [13, 14, 15] on arbitrary conducting and composite geometries. The obtained approach combines the canonical aperture coupling approach with the EFIE-PMCHWT formulation for composite structures [16, 17, 18, 19, 20, 21, 22]. Each WPE problem in this approach is divided into equivalent sub-problems for internal and external regions, which are solved using the MoM. The internal WPE problems are solved using waveguide modal expansion in the port plane, while the external problems are solved using the equivalence principle to reduce these problems to the systems of algebraic equations for unknown electric and magnetic currents. The obtained solution also considers the problem of material junctions between adjacent surfaces, considered in [19, 20, 21, 22].

The developed approach is validated on radiation and coupling problems for coaxial ports by comparing the simulated results with those from other approaches and measurements. In addition, this approach is applied to practical EMC problems for microwave antennas fed by coaxial ports. The MoM calculations were performed using the TriD numerical code incorporated in the EMCoS Studio software package [23].

Advertisement

2. Waveguide port approach for conducting geometry

2.1 Dividing the original problem into equivalent problems

Figure 1a illustrates a canonical waveguide port problem for conducting geometry. This geometry consists of a semi-infinite waveguide 1with perfect electric conducting (PEC) walls and a microwave structure 2, which is yet supposed to be conductive. We intend to create port Pin waveguide 1to divide the geometry into two regions (Aand B) to truncate the mesh in the region Aand impose appropriate termination conditions in the port plane.

Figure 1.

(a) Waveguide port problem for conducting geometry; (b) Equivalence for the internal regionA; (c) Equivalence for the external regionB.

For this purpose, we follow the classical approach for the aperture problem [9] to divide the original problem into two equivalent problems, as shown in Figure 1b and c. We introduce a perfectly conducting surface Sainto the port plane Pto separate regions Aand Band consider two equivalent sub-problems: internal (for region A) and external (for region B). In addition, we introduce equivalent magnetic currents and Mon both sides of Sato restore the tangential electric fields on the boundary surface Sa. Let us consider these equivalent problems separately.

2.2 Formulation of the internal equivalent problem

Consider an internal equivalent problem for the region A. The total EM field in the region Ais composed of the incident field Einc,Hincand the reflected field Eref,Hrefgenerated by magnetic currents in the presence of a conductor. According to the equivalence principle [24], these currents are related to the total electric field EASaon the port surface Saby the relation:

M=n×EASa=n0×EASaE1

where nis the internal normal in the region A, and n0=nis the propagation direction of the incident wave.

Equation (1) relates the total electric field at the port surface Sato magnetic currents depending on the geometric and material properties of the external region B. The internal equivalent problem is to find the modal expansion of the total EM field at the port surface Sathrough these currents.

2.3 MoM solution of the internal equivalent problem

The total EM field in the region Aon the port surface Sacan be generally written as the sum of the incident (+) and reflected (−) TEM (if exists), TE, and TM modes [2, 13]:

EASa=a0++a0e0TEM+s=1NTEas++asesTE+s=1NTMbs++bsesTME2
HASa=1Za0+a0n0×e0TEM+s=1NTE1ZsTEas+asn0×esTE+s=1NTM1ZsTMbs+bsn0×esTME3

where as±and bs±are the mode amplitudes, e0TEM, esTEand esTMare the transverse modal functions of TEM, TE and TM waves with wave impedances Z=μ/ε, ZsTE=ωμ/γsand ZsTM=γs/ωε, respectively, εand μare the permittivity and permeability of the medium, Γsis the propagation constant, and NTEand NTMare the numbers of accounted TE and TM modes, respectively.

Next, we use the MoM to relate the mode amplitudes of the reflected fields in (2), (3) to the magnetic currents M. The port surface Sais discretized into planar patches, and the unknown magnetic currents are approximated as

M=n=1NaMnfnE4

where fnare linear independent basis functions (BFs), Mnare unknown expansion current coefficients, and Nais the number of these BFs on the surface Sa. Substituting now (4) and (2) into (1), multiplying both sides by n0×e0TEM, n0×esTEand n0×esTM, respectively, and integrating over the port surface Sa, we relate the amplitudes of the reflected waves with those of the incident waves and magnetic current coefficients:

as=as+n=1NaMnTnsRs,bs=bs+n=1NaMnTnsRss=012E5

where

Tn0=Safnn0×e0TEMdS,Tns=Safnn0×esTEdS,Tns=Safnn0×esTMdSE6
R0=San0×e0TEMn0×e0TEMdS,Rs=San0×esTEn0×esTEdS,Rs=San0×esTMn0×esTMdS.E7

Substitution of (5) into (2), (3) determines the total electric and magnetic fields on the port surface in region A through the still unknown magnetic currents.

2.4 Formulation of the external equivalent problem for conducting geometry

Consider now an external equivalent problem for conducting geometry. The scattered EM field in the external region Bin Figure 1c of the conducting geometry is produced by electric currents Jflowing over surfaces Saand Scand equivalent magnetic currents Mat the surface Sa, which can be written as:

EBscJM=LEJJ+LEMME8
HBscJM=LHJJ+LHMME9

where LEJ, LEM, LHJand LHMare the linear integro-differential operators of electric and magnetic fields applied to the electric and magnetic currents, respectively. Applying the boundary conditions for the tangential electric and magnetic fields on the surfaces Saand Sc, we obtain the following system of integral equations for the unknown electric and magnetic currents Jand M

EBscJMtanSa+Sc=0just outsideSaE10
HBscJMtanSa+Sc=HASatanjust insideSaE11

2.5 MoM solution of the external equivalent problem for conducting geometry

To obtain the MoM solution to the BC (10) and (11), we consider, along with Eq. (4), the following expansion for an unknown electric current J:

J=n=1Na+NcInfn,E12

where fnare the BFs taken the same as for the expansion of magnetic currents in (4), Inare the unknown expansion current coefficients on the surfaces Saand Sc, and Naand Ncare the numbers of these BFs on these surfaces. Substitution of expansions (4) and (12) in (8) and (9) gives the following expressions for the EM field in region B:

EBscJM=n=1Na+NcInLEJfn+n=1NaMnLEMfn,E13
HBscJM=n=1Na+NcInLHJfn+n=1NaMnLHMfn.E14

Substituting now (3), (5), (13) and (14) into (10) and (11), introducing the boundary operators L̂JJ=LEJSa+Sc, L̂JM=LEMSa+Scoutside, L̂MJ=LHJSainsideand L̂MM=LHMSaand testing the resulting equations with appropriate weighting functions w1r, w2r,…, wmrleads to the following system of linear algebraic equations

ZmnJJZmnJMZmnMJZmnMMInMn=0VmWE15

with elements defined as:

ZmnJJ=wmL̂JJfn,ZmnJM=wmL˜JM+12n×fn,E16
ZmnMJ=wmL˜MJ+12n×fn,ZmnMM=wmL̂MMfn+QmnW,E17
QmnW=T̂m0WTn0R0+s=1NTET̂msWsTMTnsRs+s=1NTMT̂msWsTMTnsRs,E18
VmW=2T̂m0Wa0++s=1NTET̂msWsTMas++s=1NTMT̂msWsTMbs+,E19

where L˜JMand L˜MJare the regular parts of the boundary operators L̂JMand L̂MJ, the notation wf=swfdSis used for the scalar product, and

T̂m0=Sawmn0×e0TEMdS,T̂ms=Sawmn0×esTEdS,T̂ms=Sawmn0×esTMdSE20

In the case of Galerkin’s procedure wm=fm, and coefficients (20) and (6) become the same. The MoM system (15) determines the solution to the waveguide port problem in the conducting geometry.

2.6 Validation of the developed approach for conducting geometry

The developed approach has been validated to simulate the scattering characteristics of a flanged coaxial line as proposed in [25, 26, 27]. Such structures are frequently used in biomedical engineering for non-destructive testing of various materials [25, 26, 27].

When modeling a coaxial line, it is convenient to choose the port plane at the output of the line to provide fast damping of evanescent waves. In this case, it can be assumed NTE=NTM=0in (2) and (3) to take taken into account only fundamental, TEM mode with the modal function e0TEM=eρ/ρlnD/d, where ρis the radial distance, eρis the unit radial vector, and Dand dare the outer and inner diameters of the coaxial waveguide.

Figure 2 shows a flanged coaxial line consisting of a coaxial waveguide section with an outer radius D/2 = 4.725 mm, an inner radius of d/2 = 1.4364 mm, and a length L= 10 mm, ended with a circular disc with a diameter 2R= 200 mm. The bottom plane of the waveguide is accepted as a waveguide port, and the structure is excited in this port by TEM mode. To validate the developed approach for conducting geometry, we analyze the case when both the waveguide and outer space have the same permittivity εr=2.05.

Figure 2.

Geometry of open-ended coaxial line flanged with a circular disc:D/2 = 4.725 mm,d/2 = 1.4364 mm,L= 10 mm, 2R= 200 mm,εr=2.05inside and outside the line.

Figure 3a and b show the magnitude and phase of the reflection coefficient at the end of the coaxial line used as the reference plane. We compare the simulation results obtained using the developed approach, the mode-matching technique [25], and the matrix pencil method [26] with measurement data [27]. Note that infinite flanges are assumed in [25, 26]. Phase data conforms to the time convention expt.

Figure 3.

(a) Magnitude and (b) phase of the reflection coefficient versus the frequency of excitation at the end of the flanged coaxial line, calculated for various approaches.

Comparison of various results shows excellent agreement between them. However, the phase characteristics obtained by our approach agree somewhat more accurately with the measurement data. Thus, the obtained results validate the developed approach to modeling a coaxial waveguide port for conducting geometries.

Advertisement

3. Waveguide port approach for composite geometry

3.1 Equivalent problems for composite geometry

Figure 4 shows the geometry of the problem, consisting of a composite structure composed of k-1homogeneous regions Di, i =1,2,…, k−1, located in the free space region D0and exposed to waveguide excitation from the waveguide region B, which will be considered as k-th region of the problem. The region Dkis a finite section of the waveguide, confined by the PEC walls, the port surface Sa, and the dielectric surface Skd, through which the structure is fed. The port surface Saseparates the region Dk(B) from the semi-infinite waveguide region Awith incident waveguide excitation. In addition, each region Diis excited, in general, by the incident field Eiinc, Hiinc.

Figure 4.

Waveguide port problem for composite geometry with waveguide excitation.

To formulate the waveguide port excitation problem through the port surface Sa, we first consider the aperture coupling problem between the waveguide regions Aand B[9]. Thus, we cover the port surface Sawith a PEC sheet and introduce equivalent magnetic currents Mand Mon both sides of Sato divide the excitation problem into two different equivalence problems: the internal problem for region A, and the external problem for region B(Dk), as done in Section 2.1. Then, the internal equivalent problem is identical to that formulated in Section 2.2 and solved in Section 2.3. The external equivalent problem requires consideration of equivalent problems for each boundary surface in regions Di, i =1,2,…, k, including the port surface Sa.

3.2 Formulation of the external equivalent problem for composite geometry

An external equivalent problem for composite geometry is reduced to a set of equivalent problems for each conducting and dielectric boundary Sicand Sidof free space region D0(i= 0), composite structure regions Di(i= 1,…,k-1), and finite waveguide region Dk(i= k). In turn, each surface Sidcomprises a set of boundary surfaces sij=DiDj(ij), being the interfaces between the regions Diand Dj.

Per the equivalence principle [24], the total EM field inside the i-th region Dican be expressed as the sum of the incident field Eiinc, Hiincand that induced by the total surface currents distributed over its boundary surface Siand radiating into a homogeneous medium with constitutive parameters εiand μiof the region Di. The total electric currents Jion the boundary surface Siconsist of conducting currents Jic, flowing on the inner sides of conducting boundaries Sic, and equivalent electric currents Jid, flowing on the inner sides of dielectric boundaries Sid. Magnetic currents in the region Diare equivalent currents Mid, flowing on dielectric boundaries Sid. In addition, in the waveguide region Dkthere are equivalent magnetic currents Mon the port surface Sa.

Unknown electric and magnetic currents can be found using the boundary conditions at the conducting boundaries of the composite structure:

Eiinc+EiscJicJidMidtanSic=0,i=0,1,,k1E21

dielectric boundaries of regions Di(i,j=0,1,,k,ij):

Eiinc+EiscJicJidMidMδiktansij=[Ejinc+Ejsc(JjcJjdMjdMδjk)]tansij,E22
Hiinc+HiscJicJidMidMδiktansij=[Hjinc+Hjsc(JjcJjdMjdMδjk)]tansij,E23

and on the port surface Saand the conducting boundary Skcof the k-th region:

Ekinc+EkscJkcJkdMkdMtanSa+Skc=0just outsideSaE24
Hkinc+HkscJkcJkdMkdMtanSa=HASatanjust insideSa,E25

where δikis the Kronecker delta, which shows that magnetic currents Mradiate only in a waveguide region Dk. The magnetic field on the right-hand side of (25) is expressed by Eq. (3). The scattered EM fields in (21)–(25) can be expressed in terms of electric and magnetic currents Jiand Miin the dielectric region Dias

EiscJiMi=LiEJJiLiEMMiE26
HiscJiMi=LiHJJiLiHMMiE27

where LiEJ, LiEM, LiHJand LiHMare linear integro-differential operators of EM fields applied to currents radiated in the i-th region. It can also be shown [19, 20, 21, 22] that the equivalent currents on opposite sides of the dielectric boundaries are related as:

Jid=Jjd,Mid=MjdonsijE28

Equation (21)(25) together with relations (26)(28) and expansions (3) represent the general (EFIE-PMCHWT) form of integral equations for a composite structure with an arbitrary excitation, including the waveguide port.

3.3 MoM solution of the external equivalent problem for composite geometry

To solve the coupled system of integral Eqs. (21)(28), we use the MoM to discretize the geometry of all boundary surfaces of the regions Di(i=1,…,k) into the planar patches and to consider the following expansions for the unknown currents:

Jkc=n=1Na+NkcIncfn,M=n=1NaMn,Jici=0k1=n=1NCInCfn,E29
Jidi=0k=n=1NdIndfn,Midi=0k=n=1NdMndfn,E30

where fnare the suitable BFs, Inc, Mn,InC, Indand Mndare the unknown expansion current coefficients, and Na, Nkc, NCand Ndare the numbers of BFs on the surfaces Sa, Skc, Sici=0k1, if any, and Sidi=0k1, respectively. Expansions (29) and (30) take into account relations (28) for unknown equivalent currents on opposite sides of the dielectric boundaries. They also consider the ratios for adjacent currents at material junctions, which are the boundaries between several media [22].

Substituting (29) and (30) into (21)(25) taking into account (3), (5), (26)(28) and testing the resulting equations with weighting functions w1r, w2r,…, wmr, defined in the range of the respective boundary operators, we obtain the following MoM system of linear algebraic equations:

ZmnJcJcZmnJcM0ZmnJcJdZmnJcMdZmnMJcZmnMM+QmnW0ZmnMJdZmnMMd00ZmnJCJCZmnJCJdZmnJCMdZmnJdJcZmnJdMZmnJdJCZmnJdJdZmnJdMdZmnMdJcZmnMdMZmnMdJCZmnMdJdZmnMdMdIncMnInCIndMnd=VmcVmM+VmWVmCVmdVmHdE31

where the matrix elements are defined as Zmnαβ=wmL̂αβfn, L̂αβis the respective boundary integral operator, superscripts α,β=JcMJCJdMd; Vmc=wmEkinc, VmM=wmHkinc, VmCi=wmEiinc, Vmdij=wmEiincEjinc, VmHdij=wmHiincHjincare the voltage elements due to the incident wave in i-th and j-th media, and the elements QmnWand VmWare the same as those expressed by (18) and (19) and determine the additional inclusions in the matrix and voltage elements due to the waveguide ports.

The MoM system (31) generalizes the solution (15) of the canonical waveguide port problem to the case of composite geometry. In the structure of the MoM matrix of this solution, blocks of waveguide excitation, complex structure, and couplings between these objects through dielectric interfaces are clearly seen.

3.4 Validation of the developed approach for composite geometry

The developed approach has been validated to simulate the scattering characteristics of a single monopole antenna, fed by waveguide excitation from a flanged coaxial line with dielectric filling. Figure 5a shows a schematic view of such antenna with a height ha= 10 mm placed above a square metallic plate of 20 mm × 20 mm, which serves as a reflector. The coaxial line has an outer diameter D= 6.98 mm, an inner diameter d= 2 mm, and a length hb= 15 mm. The line bottom end is accepted as a waveguide port, and the input impedance of the antenna at this port is simulated for various dielectric fillings of the line.

Figure 5.

Single monopole antenna fed by a flanged coaxial line: (a) with dielectric filling; (b) without dielectric filling.

Figure 6 shows a comparison of the input impedances, calculated by the developed approach for the model of Figure 5a with εr=1.0001, by the WPE approach for the conducting model of Figure 5b, and by discontinuous Galerkin time-domain (DGTD) method [28]. An excellent agreement between the obtained results is seen, which confirms the equivalence and correctness of both WPE approaches (for conducting and composite geometries) for very low dielectric fillings of coaxial lines.

Figure 6.

Comparison of the input impedances of a monopole antenna in the port plane, calculated by the MoM forεr=1.0001andεr=1with DGTD method.

Figure 7 shows a comparison of the input impedances, calculated for the model of Figure 5a with εr=2.25using the developed approach and DGTD method. An excellent agreement between both results is seen, which validates our approach to treat arbitrary dielectric and geometric parameters of composite structures with waveguide port excitation.

Figure 7.

Comparison of the input impedances of a monopole antenna in the port plane for the dielectric filling of a coaxial lineεr=2.25calculated by MoM and DGTD method.

Comparison of Figures 6 and 7 shows that the use of dielectric filling of the coaxial line shifts the resonances of the input impedance to lower frequencies. In addition, this leads to a change in the line’s characteristic impedance from 75 Ω in Figure 6 to 50 Ω in Figure 7. Thus, the developed WPE approach for composite geometries covers a wider area of geometries and provides more control over the characteristics of the analyzed structures.

Advertisement

4. Waveguide port approach in coupling problems

4.1 Problem formulation

Consider the coupling problem between several composite structures fed by waveguide excitations. Although each structure can be formed from an arbitrary number of dielectric regions, for simplicity, we will consider only one-region structures with composite (dielectric and conducting) boundaries. Figure 8 shows the geometry of the problem consisting of Nwaveguides Wiradiating into dielectric regions Di, i= 1,2,…N, surrounded by closed surfaces SDiwith partially conducting boundaries SDicand inward unit normal nDi. Waveguides Wiare filled, in general, by dielectrics with permittivities εiand permeabilities μi, and the regions Diare filled by dielectrics with parameters εDiand μDi. An outer space region D0is a free space with material parameters ε0, μ0.

Figure 8.

Geometry of the problem.

The waveguide ports Piin cross-sections Siadivide the waveguides Wiinto semi-infinite regions Aiand finite regions Bito truncate the mesh in regions Aiwith incident waveguide excitation and act as excitation sources of composite regions Dithrough the dielectric boundaries SDiBidbetween the regions Diand Bi. Each region Bi, Diand D0is also excited, for generality, by the impressed EM field Eαinc, Hαinc, α=Bi,Di,D0.

To formulate the waveguide port excitation problems through the port surfaces Sia, we consider the aperture coupling problems between the regions Aiand Bito divide an original problem into two sets of equivalence problems: internal problems for regions Aiand external problems for regions Bi, Diand D0. For this purpose, we cover the port surfaces Siawith PEC sheets and introduce equivalent magnetic currents Miand Mion both sides of Siato restore tangential electric fields on the port surfaces Sia.

4.2 Solution of the internal equivalent problem

The internal equivalent problems for the considered geometry are similar to those formulated in Section 2.2 and implemented in Section 2.3. According to the equivalence principle [24], the magnetic currents in the regions Aiare related to the total electric field EAiSiaon the port surface Siaby the relation:

Mi=ni×EAiSia=n0i×EAiSiaE32

where niis an inward normal in the regionAi, and n0i=niis the propagation direction of the incident wave. Thus, the solution of the internal problem is expressed by formulas analogous to those obtained in Section 2.3 with adding the index i, when necessary.

4.3 Formulation of the external equivalent problem

When considering the external equivalent problem, let SBicbe the conducting boundary of the region Bi, including the inner sides of the waveguide walls and the conductive part of the boundary surface between the regions Biand Di; SDicis the conductive part of the boundary surface SDi, and SD0cis the conducting boundary of the region D0, including the outer sides of the waveguide walls and all conducting boundaries between the regions D0and Di. Further, SDiBidis the dielectric boundary between the regions Diand Bi, and SDiD0dis the dielectric boundary between the regions Diand D0. Per the equivalence principle [24], the dielectric boundaries between different regions can be replaced by oppositely directed equivalent electric and magnetic currents flowing on both sides of the dielectric interfaces.

The EM field in the waveguide region Biis created by electric currents JBicflowing along the port surface Siaand conducting surface SBic, equivalent electric and magnetic currents JDiBidand MDiBidflowing along the dielectric interfaces SDiBid, and equivalent magnetic currents Miflowing along the port surface Sia. The EM field in the region Diis created by electric currents JDicflowing along the conducting surfaces SDic, equivalent currents JDiBidand MDiBidflowing on dielectric boundaries SDiBidbetween the regions Diand Bi, and equivalent currents JDiDd0and MDiDd0flowing on dielectric boundaries SDiD0dbetween the regions Diand D0. The field in the free space region D0is created by electric currents SD0cflowing along the total conducting boundary of the region D0, and equivalent currents JDiDd0and MDiDd0at dielectric boundaries between the regions Diand D0.

The unknown currents JBic,Mi,JDic,JDiBid,MDiBid,JD0c,JDiD0d,MDiD0dcan be found from the boundary conditions on the port surface and the conducting boundaries of the waveguide region Bi:

EBiinc+EBiscJBicMiJDiBidMDiBidtanSia+Sic=0just outsideSiaE33
HBiinc+HBiscJBicMiJDiBidMDiBidtanSia=HAiSiatanjust insideSia,E34

and the boundary conditions on the conducting and dielectric boundaries of the regions Diand D0:

EDiinc+EDiscJDicJDiBidMDiBidJDiD0dMDiD0dtan=0onSDicE35
EBiinc+EBiscJBicMiJDiBidMDiBidtan=EDiinc+EDiscJDicJDiBidMDiBidJDiD0dMDiD0dtanonSDiBidE36
HBiinc+HBiscJDicMiJDiBidMDiBidtan=HDiinc+HDiscJDicJDiBidMDiBidJDiD0dMDiD0dtanonSDiBidE37
ED0inc+i=1NED0scJD0cJDiD0dMDiD0dtan=0onSD0cE38
EDiinc+EDiscJDicJDiBidMDiBidJDiD0dMDiD0dtan=[ED0inc+ED0scJD0c+i=1NED0sc(JDiD0dMDiD0d)]tanonSDiD0dE39
HDiinc+HDiscJDicJDiBidMDiBidJDiD0dMDiD0dtan=[HD0inc+HD0scJD0c+i=1NHD0sc(JDiD0dMDiD0d)]tanonSDiD0dE40

The scattered EM fields in (33)(40) are related to the equivalent electric and magnetic currents by Eqs. (26) and (27). After substituting (26) and (27) into (35)(40), Eqs. (35)(40) represent a coupled system of integral equations in terms of unknown currents for solving the coupling problem between several composite structures.

4.4 MoM solution of the external equivalent problem

To solve the boundary problem (35)(40), we use the following MoM expansions for the unknown currents:

JBici=1N=n=1Na+NBcIncBfn,Mii=1N=n=1NaMnfn,JDici=1N=n=1NDcIncDfn,JD0c=n=1ND0cIncD0fn,E41
JDiBidJDiD0di=1N=n=1NdIndfn,MDiBidMDiD0di=1N=n=1NdMndfn,Mii=1N=n=1NaMnfnE42

where fnare the suitable BFs, IncB,IncD,IncD0,Ind,Mndand Mnare the unknown expansion current coefficients, and Na,NBc,NDc,ND0candNdare the numbers of these BFs on the surfaces Siai=1N, SBici=1N, SDici=1N,SD0cand SDiBidSDiD0di=1N, respectively. Substituting now (41) and (42) into (35)(40) with an accounting of (3), (5), and (26) and (27) for each i-th region and testing the obtained equations with weighting functions w1r, w2r,…, wmr, defined in the range of the respective boundary operators, we obtain the following MoM system of linear algebraic equations:

ZJBcJBcZJBcM00ZJBcJdZJBcMdZMJBcZMM+QW00ZMJdZMMd00ZJDcJDc0ZJDcJdZJDcMd000ZJD0cJD0cZJD0cJdZJD0cMdZJdJBcZJdJDcZJdJDcZJdJD0cZZJdJdZJdMdZMdJBcZMdMZMdJDcZMdJD0cZMdJdZMdMdIcBMIcDIcD0IdMd=VcBVM+VWVcDVcD0VdVHdE43

where the elements of the block matrices are defined as: Zmnαβ=wmL̂iαβfn, L̂αβis the respective boundary integral operator, superscripts α,β=JBcMJDcJD0cJdMd, voltage elements are defined in the same way as in Eq. (31), and the elements of the block matrices QWand VWare expressed by (18) and (19) for each i-th feeding waveguide and determine the additional inclusions in the matrix and voltage elements due to the waveguide ports. The MoM system (43) defines a solution to the coupling problem between several composite geometries. In the structure of the MoM matrix of this solution, blocks of waveguide excitations, complex geometries, and couplings between them are clearly seen.

4.5 Validation of the developed approach for coupling problems

The developed approach has been validated on a two-element antenna array fed by coaxial waveguide ports by comparing the simulation results obtained using the developed MoM approach and the DGTD method [28]. Figure 9 shows a schematic view of two identical monopole antennas flanged over the PEC plate and fed by coaxial waveguides with generally different diameters and dielectric fillings. The monopoles located at a distance La= 40 mm from each other have the same height ha= 10 mm above the PEC plate with a width W= 40 mm and a length L= 80 mm, which serves as a reflector. Coaxial waveguides have the same inner diameter d1 = d2 =2 mm, but generally different outer diameters D1 and D2 and relative permittivities ε1 and ε2. The depth of each coaxial waveguide under the flange is hb= 15 mm, and its end is taken as the reference plane of the waveguide port.

Figure 9.

Schematic view of an array of two identical monopole antennas fed by coaxial waveguides and flanged above the PEC plate.

Figure 10 shows the real and imaginary parts of the transmission coefficient S21=a02/a01+between waveguide ports 1 and 2 with the same radii and dielectric fillings: D1/2 = D2/2 = 6.65 mm, and εr1 = εr2 = 5.17, which leads to the same characteristic impedances: Zc1=Zc2= 50 Ω. The developed MoM approach and the DGTD method are compared. The first antenna in these simulations is considered active, and the second is passive. Comparison of these results shows very good agreement between them over a wide frequency range from 1 GHz up to 10 GHz. This validates the developed approach in modeling coupling problems for coaxial waveguide ports with the same characteristic impedance.

Figure 10.

Transmission coefficient between the of the antenna array waveguide ports with the same parameters of the feeding coaxial waveguides.

Figures 11 and 12 show a comparison of the transmission coefficient S21=a02/a01+Zc1/Zc2between waveguide ports 1 and 2, calculated by the MoM and DGTD method for different parameters of coaxial waveguides. Figure 11 is made for the same fillings of waveguides: εr1 = εr2 = 2.25, but with different outer radii: D1/2 =3.49 mm and D2/2 = 6.52 mm, while Figure 12 is performed for different fillings: εr1 = 4 and εr2 = 1.78, but with the same outer radii D1/2 = D2/2 = 5.3 mm. Both cases result in characteristic impedances of waveguides Zc1= 50 Ω and Zc2= 75 Ω. Comparison of the MoM and DGTD results again shows very good agreement between both simulated results, which validates the developed approach to modeling coupling problems for coaxial waveguide ports with different characteristic impedances.

Figure 11.

Transmission coefficient between the waveguide ports of the antenna array for the same permittivitiesεr1 =εr2 = 2.25, but different outer radii:D1/2=3.49 mm,D2/2 = 6.52 mm.

Figure 12.

Transmission coefficient between the waveguide ports of the antenna array for the same outer radiiD1/2= D2/2 = 5.3 mm, but different permittivities:εr1 = 4 andεr2 = 1.78.

Advertisement

5. Application of waveguide port approach

The obtained approach has been applied to practical EMC problems for microwave antennas fed by coaxial waveguides. Such waveguides are the most commonly used to excite microwave antennas and electronic devices. This excitation usually uses microwave coaxial connectors, such as BNC and SMA.

5.1 Modeling of two branches feeding large printed UWB antenna

First, based on the measurement data [28], a printed ultra-wideband (UWB) antenna is considered. Figure 13 shows a schematic view of a large printed UWB antenna with a two-branch-feed, the bottom of which is connected to the core of a 50 Ω SMA connector with waveguide excitation, the covering of which is connected to a metal plate serving as a reflector. The bottom end of the connector is accepted as a waveguide port, and the input impedance of the UWB antenna at the waveguide port is measured and simulated.

Figure 13.

Schematic view of a large printed UWB antenna with two branch feed connected to a coaxial waveguide port.

The model of a printed UWB antenna is a square metal patch with a length La= 40 mm and a width Wa= 40 mm, printed on a dielectric substrate with a length Lb= 43 mm, a width Wb= 47.5 mm, a thickness t= 1.5 mm and material parameters εrd= 4.4 and tanδd= 0.02. The antenna is connected to a two-branch-feeding strip with a total width t1 = 15 mm, a distance between the branches t2 = 11 mm, and a height of the branches h1 = 3.5 mm. The UWB antenna is placed at a height h2 = 3 mm above a metallic plate of a length L= 275 mm and a width W= 207 mm and is connected to the SMA connector. The model of the SMA connector is represented by a coaxial waveguide with an outer radius D/2 = 2.125 mm, inner radius d/2 = 0.635 mm, and a length Lcon = 6.8 mm, filled with a polyethylene dielectric with relative permittivity εr= 2.24 and loss tangent tan δ = 0.005.

Figure 14 shows a comparison of simulated input impedances of a printed UWB antenna at the waveguide port with measurement results [28]. Comparison of the simulated results with measurement data shows a good agreement between them in a wide frequency range from 1 to 10 GHz. This validates the developed approach to modeling the composite antenna geometries fed by coaxial waveguides with dielectric filling.

Figure 14.

Comparison of the simulated and measured input impedances of a printed UWB antenna.

5.2 Coupling problem between GPS and SDARS antennas

In conclusion, based on the measurement data [28], the coupling between the GPS and SDARS patch antennas was analyzed in the frequency range from 1 GHz to 3 GHz. Figure 15 shows the measurement setup (a) and its schematic view (b) for studying the coupling between active GPS and passive SDARS antennas, separated by a distance of d= 4 cm. Both antennas are fed by 50 Ohm coaxial lines with standard SMA connectors with parameters described in Section 5.1.

Figure 15.

Measurement setup (a) and its schematic view (b) for the coupled GPS and SDARS patch antennas.

The parameters of the setup are the following. The SDARS antenna is a square metallic patch of 32 mm × 32 mm size with two opposite cut corners, printed on a dielectric substrate with dimensions 34 mm × 34 mm × 3.25 mm and εr= 4.1. The GPS antenna is constructed by a square metallic patch of 21 mm × 21 mm size with truncated corners, printed on a 25 mm × 25 mm × 4 mm dielectric substrate with εr1 = 20.34. Both patch antennas are mounted on a 190 mm × 145 mm metal plate.

Figure 16 shows a comparison of the transmission coefficient between active GPS and passive SDARS patch antennas, obtained by the developed MoM approach and measurements. A pretty good agreement between the simulated results and measured data in the frequency range of 1–3 GHz is observed. This comparison validates the developed waveguide port approach with measurements to model coupling problems between different composite geometry antennas with coaxial waveguide ports.

Figure 16.

Transmission coefficient between GPS and SDARS patch antennas.

Advertisement

6. Conclusion

The MoM-based waveguide port approach was developed to model waveguide port excitation problems on arbitrary conducting and composite geometries. The developed approach was validated for modeling radiation and coupling problems for coaxial ports by comparing the simulated results with those obtained by other approaches and measurements. The approach has been applied to practical EMC problems for microwave antennas fed by coaxial connectors. A good agreement between the simulated and measured results has been demonstrated. The efficiency of the developed approach for solving various complex problems with waveguide excitation has been verified.

References

  1. 1. Lou Z, Jin JM. An accurate waveguide port boundary condition for the time-domain finite-element method. IEEE Transactions on Microwave Theory and Techniques. 2005;53(9):3014-3023. DOI: 10.1109/APS.2005.1551498
  2. 2. Jin JM, Riley DJ. Finite Element Analysis of Antennas and Arrays. Hoboken: Wiley; 2008. 467p
  3. 3. Alimenti F, Mezzanotte P, Roselli L, Sorrentino R. Modal absorption in the FDTD method: A critical review. International Journal of Numerical Model. 1997;10:245-264
  4. 4. Alimenti F, Mezzanotte P, Roselli L, Sorrentino R. A revised formulation of modal absorbing and matched modal source boundary conditions for the efficient FDTD analysis of waveguide structures. IEEE Transactions on Microwave Theory and Techniques. 2000;48(1):50-59. DOI: 10.1109/22.817471
  5. 5. Flisgen T, Heller J, Rienen U. Time-domain absorbing boundary terminations for waveguide ports based on state-space models. IEEE Transactions on Magnetics. 2014;50(2):145-148. DOI: 10.1109/TMAG.2013.2283065
  6. 6. Li P, Jiang LJ, Bağci H. Transient analysis of dispersive power-ground plate pairs with arbitrarily shaped antipads by the DGTD method with wave port excitation. IEEE Transactions on Electromagnetic Compatibility. 2017;59(1):172-183
  7. 7. Štumpf M, Šeděnka V, Kadlec P. On modeling of excitation ports in the time-domain contour-integral method. In: Proceedings of the International Conference on Electromagnetics in Advanced Applications (ICEAA’17); 11–15 September 2017; Verona, Italy. IEEE; 2017. pp. 1292-1294. DOI: 10.1109/ICEAA.2017.8065509
  8. 8. Harrington RF. Field Computation by Moment Methods. New York: Macmillan; 1968
  9. 9. Harrington RF, Mautz JR. A generalized network formulation for aperture problems. IEEE Transactions on Antennas and Propagation. 1976;24(6):870-873. DOI: 10.1109/TAP.1976.1141420
  10. 10. Wang T, Harrington R, Mautz JR. Electromagnetic scattering from and transmission through arbitrary apertures in conducting bodies. IEEE Transactions on Antennas and Propagation. 1990;38(11):1805-1814. DOI: 10.1109/8.102743
  11. 11. Bunger R. Moment-method analysis of arbitrary 3-D metallic N-port waveguide structures. IEEE Transactions on Microwave Theory and Techniques. 2000;48(4):531-537
  12. 12. Wen D, Chen H, Zhang J. Radiation Characteristic Analysis of Ridged Horn by Aperture Port Excitation. Xi’an, China: IEEE; 2012. pp. 759-761
  13. 13. Bogdanov F, Chochia I, Svanidze L, Jobava R. Incorporation of MoM-based waveguide port model into the mixed conducting and dielectric geometry. In: Proceedings of the International Symposium on Electromagnetic Compatibility (EMC Europe’17); 4–8 September 2017; Angers, France. IEEE; 2017. pp. 1-6. DOI: 10.1109/EMCEurope.2017.8094718
  14. 14. Bogdanov F, Svanidze L, Gheonjian A, Eremyan D, Kutchadze Z, Jobava R. Application of MoM-based waveguide port approach to the analysis of EM coupling problems related to microwave antennas. In: Proceedings of the European Conference on Antennas and Propagation (EUCAP’18); 9–13 April 2018; London, England. IEEE; 2018. pp. 1-5. DOI: 10.1049/cp.2018.1194
  15. 15. Bogdanov F, Chochia I, Svanidze L, Jobava R. Validation of MoM-based solution of waveguide port problem for composite structures applied to microwave antenna and PCB geometries. In: Proceedings of the European Microwave Conference (EuMC’19); 1–3 October 2019; Paris, France. IEEE; 2019. pp. 642-645. DOI: 10.23919/EuMC.2019.8910683
  16. 16. Arvas E, Rahhal-Arabi A, Sadigh A, Rao SM. Scattering from multiple conducting and dielectric bodies of arbitrary shape. IEEE Antennas and Propagation Magazine. 1991;33(2):29-36. DOI: 10.1109/74.88184
  17. 17. Goggans PM, Kishk AA, Glisson AW. Electromagnetic scattering from objects composed of multiple homogeneous regions using a region-by-region solution. IEEE Transactions on Antennas and Propagation. 1994;42(6):865-871. DOI: 10.1109/8.301713
  18. 18. Shin J, Glisson AW, Kishk AA. Analysis of combined conducting and dielectric structures of arbitrary shapes using an E-PMCHW integral equation formulation. In: Proceedings of IEEE-AP-S Symposium (AP-S’00). Vol. 4. 2000. pp. 2282-2285
  19. 19. Kolundzija BM. Electromagnetic modeling of composite metallic and dielectric structures. IEEE Transactions on Microwave Theory and Techniques. 1999;47(7):1021-1032
  20. 20. Yiёa-Oijala P, Taskinen M, Sarvas J. Surface integral equation method for general composite metallic and dielectric structures with junctions. Progress In Electromagnetics Research. 2005;52:81-108. DOI: 10.2528/PIER04071301
  21. 21. Carr M, Topsakal E, Volakis JL. A procedure for modeling material junctions in 3-D surface integral equation approaches. IEEE Transactions on Antennas and Propagation. 2004;52(5):1374-1379. DOI: 10.1109/TAP.2004.827247
  22. 22. Bogdanov F, Svanidze L, Jobava R. MoM solution to scattering problem on multi-region composite structures with various type material junctions. In: Proceedings of International Seminar/Workshop (DIPED’18); 24–27 September 2018; Tbilisi. pp. 13-18. DOI: 10.1109/DIPED.2018.8543281
  23. 23. EMCoS. EMCoS Studio 2021 [Internet]. Available from:http://www.emcos.com
  24. 24. Harrington RF. Time-Harmonic Electromagnetic Fields. New York: McGraw-Hill; 1961
  25. 25. Cho YH. Analytic and numerically efficient scattering equations for an infinitely flanged coaxial line. PIERS Letters. 2012;28:149-158. DOI: 10.2528/PIERL11110901
  26. 26. Tan W, Shen Z. Efficient analysis of open-ended coaxial line using Sommerfeld identity and matrix pencil method. IEEE Microwave and Wireless Components Letters. 2008;18(1):7-9. DOI: 10.1109/LMWC.2007.911971
  27. 27. Mosig JR, Besson JE, Gex-Fabry M, Gardiol FE. Reflection of an open-ended coaxial line and application to non-destructive measurement of materials. IEEE Transactions on Instrumentation and Measurement. 1981;30(1):46-51. DOI: 10.1109/TIM.1981.6312437
  28. 28. Eremyan D, Kutchadze Z, Badzagua I, Gheonjian A, Jobava R. Validation of discontinuous Galerkin method. EMCoS Internal Report. Tbilisi, Georgia; 2013

Written By

Faik Bogdanov, Irina Chochia, Lily Svanidze and Roman Jobava

Submitted: January 29th, 2022 Reviewed: February 1st, 2022 Published: March 24th, 2022