# An Efficient Analysis Method for Cylindrical Conformal Microstrip Antenna Fed by Microstripline

Authors: Chengyou Yin; Mengzhong Hu
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

Firstly, a domain-division solution is provided in this paper, which can solve Hankel function’s singular problem of Sommerfeld integral in calculating the spatial domain Green function for cylindrically stratified media and realize the computation for all elements ofZ-matrix in the method of moment. Then, the arbitrary shape cylindrical conformal microstrip antenna (CCMA) fed by microstripline is theoretically analysed using the RWG basis function. Numerical results are presented in the form of mutual coupling between two current modes and electromagnetic characteristics of a CCMA fed by microstripline. The results show that, comparing with the planar counterpart microstrip antenna, the reflection coefficient of the CCMA increases, the current on the microstripline varies sharply, and the backward radiation field appears. However, the frontward radiation field of the CCMA is similar to the planar case.


---

## Body

## 1. Introduction

Microstrip antennas find many applications in satellite and mobile communication because of their many advantages, such as low profile, light weight, easy fabrication, and conformability to mounting hosts. With the development of technology, people find it is not enough to do researches only on planar microstrip antennas, so they begin to study and analyse conformal microstrip antennas, such as cylindrical, spherical, conical, or paraboloidal conformal microstrip antennas. This paper just focuses on cylindrical conformal microstrip antenna (CCMA). It is found that the development of analysis method for CCMA is not as mature as the planar case. At present, researchers usually use the transmission line model, cavity model [1–3], or their improved models to analyse CCMA in engineering projects. However, these models are only accurate for regular shape (such as rectangle, triangle, or circle) CCMA with thin substrate and cannot be applied to analyse arbitrary shape CCMA. As full-wave solutions are more accurate and applicable to many structures, some researchers begin to study the MoM-based solution for CCMA, such as [4–13] in recent years. Reading these references, one can find that the spatial domain Green function of cylindrically stratified media is the focused problem. A quasistatic images extraction technique has been proposed to solve the spatial domain Green function of electric field and mixed potential, respectively, in [4, 5], but when the source and observation points are located at the same radial distance from the axis of conductive cylinder, it is not valid along axial line (defined as ρ=ρ′ and φ=φ′). An eigenfunction solution of the Green function for arbitrary source and observation locations has been provided in [6–8]; however it needs a lot of computer resources, time consuming, and does not adapt to analyse big size or complex shape CCMA. The spatial and spectral MoM for CCMA is studied in [8–11], but all of them used the rooftop basis function which is not the general basis function compared with RWG for complex shape CCMA. In addition, recent relative MoM-based electromagnetic simulation software (such as IE3D, FEKO, or ADS) does not realize the calculation of Green’s function of cylindrically stratified media too. Therefore, based on the former research works, this paper will provide a method which can analyse arbitrary shape CCMA because of the RWG basis function used in MoM. The paper is documented in the following structure. In Section 2, theory and formulation of the domain-division solution for the Green’s function for CCMA will be introduced in detail firstly. Then, based on RWG-MoM, theoretical analysis of arbitrary shape CCMA fed by microstripline will be done. In Section 3, a simulation example for investigating the effectiveness of the domain-division solution is presented and the electromagnetic characteristics of a CCMA fed by microstripline are gotten in another simulation example, in which results are compared with those of the planar case. A time convention of ejωt has been adopted in this paper.


## 2. Theory and Formulation

As we know, analysis of CCMA using MoM needs calculating the spatial Green function when source and observation points are located at the same radial distance from the axis of cylinder (ρ=ρ′), but, in this case, integral kernel function converges very slowly for electrically large conductive cylinders, because the Green’s function involves Bessel and Hankel functions along with their derivatives. So the expressions of Green function gotten from [12, 13] cannot be directly used in the analysis of CCMA and some special work should be done. Combining the methods of [4–11], a domain-division solution can be provided to solve this problem which divides the computation of Green’s function into two domain’s problems. The first one is the unaxial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|>0.05), and the other one is the axial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|<0.05).


### 2.1. Solution of the Unaxial Line Problem

Based on the formulation of spectral domain Green’s function in [12, 13], the expressions of the spectral domain Green function components for the nth Harmonic wave can be rewritten as
(1)G~uvEnkze3=(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)fuv(n,kz).
For different uv, the values of e1, e2, and e3 in the above formula are show in Table 1.Table 1
The values ofe1,e2, and e3 for different uv.
u
v
e
1
e
2
e
3
z
z
1
0
0
φ
z
0
1
1
z
φ
0
1
1
φ
φ
0
2
0Hence, the expressions of the spectral domain Green function for all tangential electric field components in cylindrical coordinates can be written as(2)G~uvEkze3=-14ω∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×fuv(n,kz)ejn(φ-φ′).Using the expressions of the spectral domain Green function in [12, 13], fuv (u and v can be z or φ) in the above formula is easy to get.In (2), when n→∞, it is found that fuv(n,kz) will converge to a function which is constant with respect to n,
(3)limn→∞fuv(n,kz)≈Cuv(kz).
With the aid of the series expression of H0(2)(kρi|ρ-ρ′|), given by
(4)∑n=-∞∞Hn(2)(kρiρ)Jn(kρiρ′)ejn(φ-φ′)=H0(2)(kρi|ρ-ρ′|),
then (2) can be expressed as：(5)G~uvEkze3=-14ω∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×(fuv(n,kz)-Cuv(kz))ejn(φ-φ′)+-14ωCuv(kz)(kρi2)e1F1uv[H0(2)(kρi|ρ-ρ′|)].
Expressions of F1uv[H0(2)(kρi|ρ-ρ′|)]in (5) can be gotten using (4) and its derivatives, which can be found in [5]. Computation of F1uv[H0(2)(kρi|ρ-ρ′|)]involves Bessel and Hankel functions’ derivatives. With the aid of formula (1.2.34) in [14], their derivatives can be calculated easily.From (5), one can find that the spectral domain Green function now includes two parts after extracting the asymptotic part with respect to n. The first part is still a summation of infinite series with respect to n, but it is now more rapidly convergent and the limits of the infinite summation can be truncated at relatively small values Nt. The second part is the asymptotic part of the spectral domain Green function which now is not an infinite summation but a function with respect to H0(2)(kρi|ρ-ρ′|).To study the convergence of formulae (2) and (5), the spectral domain Green function is computed by (2) and (5), respectively. The results are shown in Figures 1, 2, and 3.Figure 1
Comparison ofG~zzE imaginary part versus Nt using (2) and (5).Figure 2
Comparison ofG~φzE imaginary part versus Nt using (2) and (5).Figure 3
Comparison ofG~φφE imaginary part versus Nt using (2) and (5).Parameters used in the above computation are set as follow:Δφ=0.05, a0=3λ0, a1=3.06λ0, εr=3.25, kz=0 for computation of G~zzE and kz=ks(1+jT1) for computation of G~φzE, G~φφE. As seen in Figures 1–3, it can be found that results of the spectral domain Green function using (2) converge difficultly even for Nt=300. Especially for G~φφE shown in Figure 3, the numerical results seem to be divergent. But the results of the spectral domain Green function using (5) converge very quickly and one can make Nt=200 to get accurate results for the spectral domain Green function.Whenkz→∞, it is found that Cuv(kz) is convergent to Cuv∞, so, extracting the asymptotic part of the spatial domain Green function components with respect to kz, the expressions for calculating the spatial domain Green function of electric field can be written as
(6)GuvE=(j∂∂z)e3{-18ωπ∫-∞∞(g~uv1E+g~uv2E)=(j∂∂z)e3=(j∂∂z)e3×cos[kz(z-z′)]dkz-18ωπ∫-∞∞(g~uv1E+g~uv2E)}+-18ωπCuv∞F2uv,g~uv1E=∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×(fuv(n,kz)-Cuv(kz))ejn(φ-φ′),g~uv2E=(Cuv(kz)-Cuv∞)F1uv[H0(2)(kρi|ρ-ρ′|)].
Using the following equations:
(7)I1=∫-∞∞H0(2)(kρi|ρ-ρ′|)e-jkz(z-z′)=2je-jki|r-r′||r-r′|,I2=∫-∞∞kρi∂[H0(2)(kρi|ρ-ρ′|)]∂(kρi|ρ-ρ′|)e-jkz(z-z′)=2j(-jks|r-r′|-1)|ρ-ρ′|e-jks|r-r′||r-r′|3,I3=∫-∞∞kρi2∂2[H0(2)(kρi|ρ-ρ′|)]∂(kρi|ρ-ρ′|)2e-jkz(z-z′)=2j[-ks2|ρ-ρ′|2|r-r′|2-(z-z′)2-j(z-z′)2×ks|r-r′|+2|ρ-ρ′|2+2j|ρ-ρ′|2ks|r-r′|]×e-jks|r-r′||r-r′|5.
Expressions for calculating F2uv in (6) can be deducted. Analysing GuvE in (6), one can find that GuvE also includes two parts. The first part is an infinite integral with respect to kz which is now fast decaying with respect to even for very small |φ-φ′| values and can be solved by the method of deformed integration path proposed in [13]. The second part is the asymptotic part of the spatial domain Green function which has closed form expressions.


### 2.2. Solution for Axial Line Problem (|ρ-ρ′|=0 and |φ-φ′|<0.05)

Although a twice-extracting technique is proposed to solve the unaxial line problem in the above section, as to axial line problem,H0(2)(kρi|ρ-ρ′|) in (6) shows more singular and is also difficult to calculate. Fortunately, [7] does much research on this problem and obtains an approximated representation of spatial Green’s function for this case by performing Watson’s transformation.WhenS→0 (S is shown in Figure 6), the method of [7] is not valid because of high-order singular problem in this condition. Even using the asymptotic extraction approach in [15] convergence is still difficult in this case. However, the Green function of planar stratified media can be used instead, because cylindrical stratified media can be treated as planar stratified media approximately when S<0.2λ0(λ0is the free space wave length). Therefore, to sum up, computation of the spatial Green function for CCMA can be divided into two parts which are shown in Figure 4. The first part is about the unaxial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|>0.05), and the other is about the axial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|<0.05). The axial line problem involves two subdomains. One is S>0.2λ0, and the other is S<0.2λ0. The method proposed in [7] is used to attack the problem in S>0.2λ0domain, and the Green function of planar stratified media is used to solve the problem in S<0.2λ0 domain. Computation for the Green function of planar stratified media can be found in [14–17].Figure 4
Calculation of the spatial Green function for MMCA.


### 2.3. Theoretical Analysis of Arbitrary Shape CCMA Fed by Microstripline

Considering CCMA fed by microstripline,Eρs does not need to be considered, so one can get
(8)[EφsEzs]=〈[G-φφEG-φzEG-zφEG-zzE][JφJz]〉.
If the source current Ji=[JφiJzi]T is excited at the end of microstripline as shown in Figure 5, then excited electric field is written as
(9)[EφiEzi]=〈[G-φφEG-φzEG-zφEG-zzE][JφiJzi]〉.Figure 5
CCMA fed by microstripline.Figure 6
Real and imaginary parts of the mutual impedance (Z12) between two identical z-directed current sources versus separation S when α=700° (unaxial line problem).Using the Galerkin testing technique, on the patch of CCMA, we arrive at(10)〈fw,Es|t〉=-〈fw,Ei|t〉,
where fw is the test basis function in MoM. Upon substituting (8) and (9) into (10), one obtains
(11)∫S∫S′[fφw·(G-φφE·Jφ+G-φzE·Jz)]dS′dS=-∫S∫Si[fφw·(G-φφE·Jφi+G-φzE·Jzi)]dSidS,∫S∫S′[fzw·(G-zφE·Jφ+G-zzE·Jz)]dS′dS=-∫S∫Si[fzw·(G-zφE·Jφi+G-zzE·Jzi)]dSidS.The current distribution on the patch of the antenna can be approximated in terms offn as
(12)Js=∑nInfn,
where fn is the RWG basis function in this paper, which is different from the traditional method in [8–10], and one can find the define about RWG in [18]. From (11) and (12), if fw is equal to fm, current coefficient In can be gotten from the following identity:
(13)In∫S∫S′[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][fn|φfn|z]dS′dS=-∫S∫Si[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][JφiJzi]dSidS.
Considering m,n=1,2,…N, formula (13) can be written in matrix form as
(14)[Z11Z12⋯Z1NZ21Z22⋯Z2N⋮⋮⋮⋮ZN1ZN2⋯ZNN][I1I2⋮IN]=[V1V2⋮VN],
where
(15)Zmn=∫S∫S′[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][fn|φfn|z]dS′dS,Vm=-∫S∫Si[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][JφiJzi]dSidS.In order to calculateVm in (15), Ji can be treated as a half-RWG basis function added at the end of the microstripline which is shown in Figure 5.After getting all the current coefficients in (14), one can obtain the current distributions on the patch and microstripline. Using GPOF method [19], current distributions on the microstripline can be expressed as
(16)I(z)≈∑i=1KIie(αi+jβi)z,
where Ii, αi, and βi are parameters approximated by GPOF method. It is found that αi≪βi, so (16) can be approximated as follows:
(17)I(z)=I1e-jβ1z+I2ejβ1z⋯IKejβKz.In (17), I1 and I2 represent the two opposite travelling waves I+ and I-. The others (I3,I4,…IK) are the high-order travelling waves which can be ignored generally. So the reflection coefficient can be found as
(18)S11=I2I1,dB(S11)=20log10|I2I1|,
and the impedance of the antenna is given by
(19)Zin=1+S111-S11Zc,
where Zc is the characteristic impedance of microstripline. In addition, to calculate the far field of the CCMA fed by microstripline, the method in [19] can be used. However, it should be noticed that this paper chooses the RWG basis function to approximate the current, so formulas for calculating current in [19] should be changed properly.


## 2.1. Solution of the Unaxial Line Problem

Based on the formulation of spectral domain Green’s function in [12, 13], the expressions of the spectral domain Green function components for the nth Harmonic wave can be rewritten as
(1)G~uvEnkze3=(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)fuv(n,kz).
For different uv, the values of e1, e2, and e3 in the above formula are show in Table 1.Table 1
The values ofe1,e2, and e3 for different uv.
u
v
e
1
e
2
e
3
z
z
1
0
0
φ
z
0
1
1
z
φ
0
1
1
φ
φ
0
2
0Hence, the expressions of the spectral domain Green function for all tangential electric field components in cylindrical coordinates can be written as(2)G~uvEkze3=-14ω∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×fuv(n,kz)ejn(φ-φ′).Using the expressions of the spectral domain Green function in [12, 13], fuv (u and v can be z or φ) in the above formula is easy to get.In (2), when n→∞, it is found that fuv(n,kz) will converge to a function which is constant with respect to n,
(3)limn→∞fuv(n,kz)≈Cuv(kz).
With the aid of the series expression of H0(2)(kρi|ρ-ρ′|), given by
(4)∑n=-∞∞Hn(2)(kρiρ)Jn(kρiρ′)ejn(φ-φ′)=H0(2)(kρi|ρ-ρ′|),
then (2) can be expressed as：(5)G~uvEkze3=-14ω∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×(fuv(n,kz)-Cuv(kz))ejn(φ-φ′)+-14ωCuv(kz)(kρi2)e1F1uv[H0(2)(kρi|ρ-ρ′|)].
Expressions of F1uv[H0(2)(kρi|ρ-ρ′|)]in (5) can be gotten using (4) and its derivatives, which can be found in [5]. Computation of F1uv[H0(2)(kρi|ρ-ρ′|)]involves Bessel and Hankel functions’ derivatives. With the aid of formula (1.2.34) in [14], their derivatives can be calculated easily.From (5), one can find that the spectral domain Green function now includes two parts after extracting the asymptotic part with respect to n. The first part is still a summation of infinite series with respect to n, but it is now more rapidly convergent and the limits of the infinite summation can be truncated at relatively small values Nt. The second part is the asymptotic part of the spectral domain Green function which now is not an infinite summation but a function with respect to H0(2)(kρi|ρ-ρ′|).To study the convergence of formulae (2) and (5), the spectral domain Green function is computed by (2) and (5), respectively. The results are shown in Figures 1, 2, and 3.Figure 1
Comparison ofG~zzE imaginary part versus Nt using (2) and (5).Figure 2
Comparison ofG~φzE imaginary part versus Nt using (2) and (5).Figure 3
Comparison ofG~φφE imaginary part versus Nt using (2) and (5).Parameters used in the above computation are set as follow:Δφ=0.05, a0=3λ0, a1=3.06λ0, εr=3.25, kz=0 for computation of G~zzE and kz=ks(1+jT1) for computation of G~φzE, G~φφE. As seen in Figures 1–3, it can be found that results of the spectral domain Green function using (2) converge difficultly even for Nt=300. Especially for G~φφE shown in Figure 3, the numerical results seem to be divergent. But the results of the spectral domain Green function using (5) converge very quickly and one can make Nt=200 to get accurate results for the spectral domain Green function.Whenkz→∞, it is found that Cuv(kz) is convergent to Cuv∞, so, extracting the asymptotic part of the spatial domain Green function components with respect to kz, the expressions for calculating the spatial domain Green function of electric field can be written as
(6)GuvE=(j∂∂z)e3{-18ωπ∫-∞∞(g~uv1E+g~uv2E)=(j∂∂z)e3=(j∂∂z)e3×cos[kz(z-z′)]dkz-18ωπ∫-∞∞(g~uv1E+g~uv2E)}+-18ωπCuv∞F2uv,g~uv1E=∑n=-∞∞(kρi2)e1ne2Hn(2)(kρiρ)Jn(kρiρ′)×(fuv(n,kz)-Cuv(kz))ejn(φ-φ′),g~uv2E=(Cuv(kz)-Cuv∞)F1uv[H0(2)(kρi|ρ-ρ′|)].
Using the following equations:
(7)I1=∫-∞∞H0(2)(kρi|ρ-ρ′|)e-jkz(z-z′)=2je-jki|r-r′||r-r′|,I2=∫-∞∞kρi∂[H0(2)(kρi|ρ-ρ′|)]∂(kρi|ρ-ρ′|)e-jkz(z-z′)=2j(-jks|r-r′|-1)|ρ-ρ′|e-jks|r-r′||r-r′|3,I3=∫-∞∞kρi2∂2[H0(2)(kρi|ρ-ρ′|)]∂(kρi|ρ-ρ′|)2e-jkz(z-z′)=2j[-ks2|ρ-ρ′|2|r-r′|2-(z-z′)2-j(z-z′)2×ks|r-r′|+2|ρ-ρ′|2+2j|ρ-ρ′|2ks|r-r′|]×e-jks|r-r′||r-r′|5.
Expressions for calculating F2uv in (6) can be deducted. Analysing GuvE in (6), one can find that GuvE also includes two parts. The first part is an infinite integral with respect to kz which is now fast decaying with respect to even for very small |φ-φ′| values and can be solved by the method of deformed integration path proposed in [13]. The second part is the asymptotic part of the spatial domain Green function which has closed form expressions.


## 2.2. Solution for Axial Line Problem (|ρ-ρ′|=0 and |φ-φ′|<0.05)

Although a twice-extracting technique is proposed to solve the unaxial line problem in the above section, as to axial line problem,H0(2)(kρi|ρ-ρ′|) in (6) shows more singular and is also difficult to calculate. Fortunately, [7] does much research on this problem and obtains an approximated representation of spatial Green’s function for this case by performing Watson’s transformation.WhenS→0 (S is shown in Figure 6), the method of [7] is not valid because of high-order singular problem in this condition. Even using the asymptotic extraction approach in [15] convergence is still difficult in this case. However, the Green function of planar stratified media can be used instead, because cylindrical stratified media can be treated as planar stratified media approximately when S<0.2λ0(λ0is the free space wave length). Therefore, to sum up, computation of the spatial Green function for CCMA can be divided into two parts which are shown in Figure 4. The first part is about the unaxial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|>0.05), and the other is about the axial line problem (defined as |ρ-ρ′|=0 and |φ-φ′|<0.05). The axial line problem involves two subdomains. One is S>0.2λ0, and the other is S<0.2λ0. The method proposed in [7] is used to attack the problem in S>0.2λ0domain, and the Green function of planar stratified media is used to solve the problem in S<0.2λ0 domain. Computation for the Green function of planar stratified media can be found in [14–17].Figure 4
Calculation of the spatial Green function for MMCA.


## 2.3. Theoretical Analysis of Arbitrary Shape CCMA Fed by Microstripline

Considering CCMA fed by microstripline,Eρs does not need to be considered, so one can get
(8)[EφsEzs]=〈[G-φφEG-φzEG-zφEG-zzE][JφJz]〉.
If the source current Ji=[JφiJzi]T is excited at the end of microstripline as shown in Figure 5, then excited electric field is written as
(9)[EφiEzi]=〈[G-φφEG-φzEG-zφEG-zzE][JφiJzi]〉.Figure 5
CCMA fed by microstripline.Figure 6
Real and imaginary parts of the mutual impedance (Z12) between two identical z-directed current sources versus separation S when α=700° (unaxial line problem).Using the Galerkin testing technique, on the patch of CCMA, we arrive at(10)〈fw,Es|t〉=-〈fw,Ei|t〉,
where fw is the test basis function in MoM. Upon substituting (8) and (9) into (10), one obtains
(11)∫S∫S′[fφw·(G-φφE·Jφ+G-φzE·Jz)]dS′dS=-∫S∫Si[fφw·(G-φφE·Jφi+G-φzE·Jzi)]dSidS,∫S∫S′[fzw·(G-zφE·Jφ+G-zzE·Jz)]dS′dS=-∫S∫Si[fzw·(G-zφE·Jφi+G-zzE·Jzi)]dSidS.The current distribution on the patch of the antenna can be approximated in terms offn as
(12)Js=∑nInfn,
where fn is the RWG basis function in this paper, which is different from the traditional method in [8–10], and one can find the define about RWG in [18]. From (11) and (12), if fw is equal to fm, current coefficient In can be gotten from the following identity:
(13)In∫S∫S′[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][fn|φfn|z]dS′dS=-∫S∫Si[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][JφiJzi]dSidS.
Considering m,n=1,2,…N, formula (13) can be written in matrix form as
(14)[Z11Z12⋯Z1NZ21Z22⋯Z2N⋮⋮⋮⋮ZN1ZN2⋯ZNN][I1I2⋮IN]=[V1V2⋮VN],
where
(15)Zmn=∫S∫S′[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][fn|φfn|z]dS′dS,Vm=-∫S∫Si[fm|φfm|z]·[G-φφEG-φzEG-zφEG-zzE][JφiJzi]dSidS.In order to calculateVm in (15), Ji can be treated as a half-RWG basis function added at the end of the microstripline which is shown in Figure 5.After getting all the current coefficients in (14), one can obtain the current distributions on the patch and microstripline. Using GPOF method [19], current distributions on the microstripline can be expressed as
(16)I(z)≈∑i=1KIie(αi+jβi)z,
where Ii, αi, and βi are parameters approximated by GPOF method. It is found that αi≪βi, so (16) can be approximated as follows:
(17)I(z)=I1e-jβ1z+I2ejβ1z⋯IKejβKz.In (17), I1 and I2 represent the two opposite travelling waves I+ and I-. The others (I3,I4,…IK) are the high-order travelling waves which can be ignored generally. So the reflection coefficient can be found as
(18)S11=I2I1,dB(S11)=20log10|I2I1|,
and the impedance of the antenna is given by
(19)Zin=1+S111-S11Zc,
where Zc is the characteristic impedance of microstripline. In addition, to calculate the far field of the CCMA fed by microstripline, the method in [19] can be used. However, it should be noticed that this paper chooses the RWG basis function to approximate the current, so formulas for calculating current in [19] should be changed properly.


## 3. Numerical Results

### 3.1. Mutual Impedance between the Current Modes

In order to assess the validity and accuracy of the method proposed in this paper, numerical results for the mutual impedance between two tangential electric current modes are obtained using the domain-division solution and compared with the results of [6, 7]. The numerical results are shown in Figures 6 and 7, and the parameters of the cylinder and substrate in these figures are set as a0=3λ0, a1=3.06λ0, εr=3.25. Figure 6 shows that results obtained by the twice-extracting technique in this paper are in good agreement with the results in [6], so the technique is effective for the unaxial line problem. As seen in Figure 7, the method of [7] is used successfully to solve axial line problem in this paper.Figure 7
Real and imaginary parts of the mutual impedance (Z12) between two identical z-directed current sources versus separation S when α=90° (axial line problem).


### 3.2. Analysis of a CCMA Fed by Microstripline

Using the method in this paper, a CCMA fed by microstripline illustrated in Figure8 is analysed. Sizes of the antenna are set as φ- and z-direction lengths of patch are 0.0615 m and 0.0604, respectively, microstripline length is 0.1278 m, the feed point is 0.02696 m away from the bottom edge’s centre of the patch in φ-direction, thickness of the substrate is 0.00234 m, the microstripline wide is 0.1 mm (Zc=184), εr=2.32, and a0=0.187451 m. Results of the normalised current distribution on microstripline, reflection coefficient, and normalised radiation patterns are presented in Figures 9–12.Figure 8
CCMA fed by microstripline.Figure 9
Comparison of the current distribution on microstripline.Figure 10
Comparison of the reflection coefficient.Figure 11
Normalised radiation field inxoy-plane (dB).Figure 12
Normalised radiation field inxoz-plane (dB).Figures9 and 10 show that resonance frequency of the antenna in planar case is 1.575 GH, but in cylindrical case, the reflection coefficient does not satisfy the request dB(S11) < −10 dB at 1.575 GHz. The reason for it may be that characteristic impedance of the microstripline in cylindrical case does not match the impedance of the antenna. So the sizes of microstripline should be designed again to match the impedance of the CCMA. If the microstripline wide is designed to be 2.1 mm (Zc=71Ω), then the reflection coefficient of the CCMA will decrease, which is shown in Figure 10.Figures11 and 12 show that backward radiation field appears in cylindrical case, but the frontward radiation field is similar to the planar case.


## 3.1. Mutual Impedance between the Current Modes

In order to assess the validity and accuracy of the method proposed in this paper, numerical results for the mutual impedance between two tangential electric current modes are obtained using the domain-division solution and compared with the results of [6, 7]. The numerical results are shown in Figures 6 and 7, and the parameters of the cylinder and substrate in these figures are set as a0=3λ0, a1=3.06λ0, εr=3.25. Figure 6 shows that results obtained by the twice-extracting technique in this paper are in good agreement with the results in [6], so the technique is effective for the unaxial line problem. As seen in Figure 7, the method of [7] is used successfully to solve axial line problem in this paper.Figure 7
Real and imaginary parts of the mutual impedance (Z12) between two identical z-directed current sources versus separation S when α=90° (axial line problem).


## 3.2. Analysis of a CCMA Fed by Microstripline

Using the method in this paper, a CCMA fed by microstripline illustrated in Figure8 is analysed. Sizes of the antenna are set as φ- and z-direction lengths of patch are 0.0615 m and 0.0604, respectively, microstripline length is 0.1278 m, the feed point is 0.02696 m away from the bottom edge’s centre of the patch in φ-direction, thickness of the substrate is 0.00234 m, the microstripline wide is 0.1 mm (Zc=184), εr=2.32, and a0=0.187451 m. Results of the normalised current distribution on microstripline, reflection coefficient, and normalised radiation patterns are presented in Figures 9–12.Figure 8
CCMA fed by microstripline.Figure 9
Comparison of the current distribution on microstripline.Figure 10
Comparison of the reflection coefficient.Figure 11
Normalised radiation field inxoy-plane (dB).Figure 12
Normalised radiation field inxoz-plane (dB).Figures9 and 10 show that resonance frequency of the antenna in planar case is 1.575 GH, but in cylindrical case, the reflection coefficient does not satisfy the request dB(S11) < −10 dB at 1.575 GHz. The reason for it may be that characteristic impedance of the microstripline in cylindrical case does not match the impedance of the antenna. So the sizes of microstripline should be designed again to match the impedance of the CCMA. If the microstripline wide is designed to be 2.1 mm (Zc=71Ω), then the reflection coefficient of the CCMA will decrease, which is shown in Figure 10.Figures11 and 12 show that backward radiation field appears in cylindrical case, but the frontward radiation field is similar to the planar case.


## 4. Conclusion

The development of the researches on CCMA is unmatured, and the Green functions of cylindrically stratified media are very difficult to calculate accurately and efficiently. Based on the former researches [1–14], this paper proposes a domain-division solution to solve Hankel function’s singular problem of Sommerfeld integral in calculating the spatial domain Green function of cylindrically stratified media and realize the computation for all elements of Z-matrix in MoM. Then, theoretical analysis of arbitrary shape CCMA fed by microstripline is done by MoM with the RWG basis function which is more general than the rooftop basis function used in the traditional method. At last, numerical results in the form of mutual impedance between two current modes reveal that the proposed method is valid and efficient. In addition, comparison of the results of a CCMA and a planar microstrip antenna fed by microstripline depicts that the current distribution on microstripline, reflection coefficient, and normalised radiation patterns are affected by antenna patch’s curvature. It also shows that, comparing with the planar counterpart microstrip antenna, the reflection coefficient of the CCMA increases, the current on the microstripline varies sharply, and the backward radiation field appears. However, the frontward radiation field is similar to the planar case. In order to improve the performance of the CCMA, the sizes of microstripline should be designed properly to match the impedance of CCMA.



---

*Source: 629748-2012-11-26.xml*