Preprint
Communication

Antenna Array on a Multilayer Dielectric Structure as a Nanostructure

Altmetrics

Downloads

107

Views

37

Comments

0

This version is not peer-reviewed

Submitted:

27 February 2024

Posted:

28 February 2024

You are already at the latest version

Alerts
Abstract
The work contains a detailed analysis of the electromagnetic field originating from a point source placed above the layered structure.The surface impedance of the multilayer dielectric structure was determined. Based on the boundary conditions, a method for determining the surface impedance of a layered structure with a perfectly conductive screen and or a magnetic screen is presented. The resonance conditions and the conditions for the excitation of surface waves are derived. The integral expression defining the generalized impedance was analyzed with a view to adapting it to numerical implementation. This allowed the development of a computer program based on which calculations were made of the effect of excitation of waves on multilayer dielectric substrates.
Keywords: 
Subject: Engineering  -   Electrical and Electronic Engineering

1. Introduction

Printed antennas have been one of the most innovative fields of antenna technology for over a dozen years. In the early 1950s, radiation from strip lines and microstrip lines was observed. However, from the point of view of the use of transmission lines, it was an undesirable phenomenon, so apart from a few articles suggesting its use in the construction of antennas, it did not arouse much interest. Only Munson’s article [1] showed that microstrip antennas can have practical applications. In 1979, the first international conference was held in La Cruses, devoted to a comprehensive approach to antennas on dielectric substrates, i.e. materials, construction plans, configuration systems and theoretical foundations.
Microstrip antennas have many interesting features, such as:
• accurate representation on the surface,
• low manufacturing cost,
• high repeatability of workmanship,
• insignificant volume,
• masking the operating frequency,
• simplicity of production provided that relatively advanced technologies are used,
• flat shape and low mass allow the use of dielectric antennas on fast flying objects without fear of deterioration of their aerodynamic properties.
Microstrip antennas are characterized by many interesting features. However, the single- or multi-layer dielectric substrate used favors the excitation of surface waves, outgoing waves, and guided waves, which propagate along the dielectric plane and disrupt the normal operation of the antenna.
Other disadvantages of microstrip antennas are:
• narrow operating band,
• limited power load.
These antennas allow for miniaturization of the antenna system, thus increasing its density. This causes mutual couplings that change the field distributions on aperture antennas and current distributions in linear antennas. This state of affairs, in turn, causes a change in the spatial characteristics of the antenna radiation and their input impedance.
The basic design of a microstrip antenna consists of a metallic strip printed on a thin, grounded (shield) dielectric substrate [2,3]. A single element can be powered by: a coaxial line run perpendicularly through the substrate (coaxial probe fed antenna) [4,5] or a microstrip line run on a dielectric substrate in the plane of the antenna [6,7]. This is particularly important when implementing planar antenna arrays (periodic antennas) due to the possibility of achieving high repeatability of execution and reducing production costs.
The first studies and publications concerned single radiating elements, which have large transverse dimensions of the antenna radiation characteristics and, therefore, small antenna gain. These defects can be removed by using periodic antenna arrays, taking into account their mutual coupling (through their mutual impedance).
There are several methods for determining mutual couplings in the periodic antenna wall. Among them, the mutual impedance method is a commonly used method. Within this approach, it is assumed that the antenna supply currents clearly determine the current distribution on the antenna dipoles and hence its radiation field (or scattered field) Es. We use the method of moments to determine these currents. Using the method developed by Richmond [8] for obtaining the elements of the impedance matrix based on the concept of reaction between sections of the base and testing functions. The power supply currents and voltages at the input of one antenna element (cell) are related to their counterparts at the input of the other antenna element through the mutual impedance matrix of the antenna system composed of many interacting elements. Hence, the supply currents of all system elements are related to the supply voltages of all system elements through the mutual impedance matrix.
Therefore, to determine the mutual impedance, it is necessary to know the electromagnetic field generated by the Hertz dipole, which leads to the need to know the Green’s function for the multilayer dielectric structure. We determine it by solving the non-homogeneous Helmholtz equation.
The spectral response of a periodic antenna structure placed in a dielectric homogeneous medium depends on the antenna geometry, medium parameters, incidence angle, polarization and excitation field geometry [17]. However, in many applications the above number of parameters determining the operation of the antenna is not sufficient to model its desired directional characteristics. Increasing the number of parameters of the antenna structure can be achieved by introducing a multilayer dielectric medium with a specific number of metalized periodic surfaces placed on flat boundaries between dielectric layers.

2. Surface impedance of the layered structure. reflectivity. transverse resonance conditions

2.1. Surface impedance of a layered structure

We assume the geometry of the system and the excitation independent of the y coordinate (Figure 1), thus reducing the issue to a two-dimensional problem
Surface impedance binds the appropriate components of the electromagnetic field tangential to the surface z=0 according to the formula [9]
n _ × n _ × E _ ( k x ) = z 0 η ( k x ) n _ × H _ ( k x )
where:
z 0 = μ 0 ε 0 specific vacuum impedance z>0,
n _ vector normal to the impedance surface (Figure 1),
η ( k x ) = 0 η e 1 / η h 0 relative (relative to the medium) surface impedance z=0.
For a TM type field (E-polarization) with components H _ = 0 , H y , 0 E _ = E x , 0 , E z
from the definition of surface impedance (1) we obtain
E ~ x = z 0 η e H ~ y ,
from Maxwell’s equations
E ˜ x = i z o k H ˜ y z ,
by comparing expression (2) and (3) we obtain the differential impedance condition for the TM problem.
H ˜ y z + i k η e H ˜ y = 0 .
We assume the following form of the field at the boundary of the centers
H ˜ y = e i ( k x x k z z ) + R e e i ( k x x + k z z ) .
From the condition of continuity (boundary) of the tangential components of the field in the z = 0 plane, we obtain the expression for the reflection coefficient
R e ( k x ) = k z k η e ( k x ) k z + k η e ( k x ) .
For a TE type field (H-polarization) with components . By definition we get
H ~ x = Y o η h E ~ y ,
and from Maxwell’s equations
H ˜ x = i Y o k E ˜ y z
from (7) and (8) we obtain the boundary admittance condition for the TE problem
E ~ z z + i k η h E ~ y = 0 .
Assuming the field form analogous (as for the TM problem), we obtain the reflection coefficient from the condition of continuity in the z = 0 plane
R h ( k x ) = k z k η h ( k x ) k z + k η h ( k x ) .

2.2 Reflectance and surface impedance of a dielectric layer lying on a perfectly conductive plane

Figure 2. Dielectric layer with an electric shield.
Figure 2. Dielectric layer with an electric shield.
Preprints 99982 g002
We consider the TM field (E-polarization) - an even field. To determine the reflection coefficient appearing in equations (6) and (10), it should be noted that it does not depend on the angle ψ in the (x, y) plane. Thus, we can analyze the issue in the selected plane from here y = 0 from where ρ = x ; k ρ = k x k y = 0 ..
For the considered TM field with components E _ = [ E x , 0 , E z ] H _ = [ 0 , H y , 0 ]
we assume the following form of the field
H ~ y ( k x ) = I d l 8 π 2 e i k x ± e ± i k z z + R e ( k x ) e i k z z d l a z > z o A cos χ ( z + z o + b ) cos χ b z b < z < z o .
From Maxwell’s equations we get:
E ˜ x = i ω ε ε o H ˜ y z
ε = 1 f o r z > z o ε ( r e l a t i v e ) z o b < z < z o .
Meaning that
E ~ x ( k x ) = i ω ε o I d l 8 π 2 e i k x x i k z e ± i k z z + R e ( k x ) e i k z z d l a z > z o A χ ε sin χ ( z + z o + b ) cos χ b z o b < z < z o
where:
χ ( k x ) = χ ( k ρ ) = k 2 n 2 k ρ 2 .
Equations (11) and (12) take into account the condition of zeroing of the electric field component on the screen z = - z o - b . We will call the obtained solution an even TM field, i.e., symmetrical with respect to the plane of the screen (perfectly conductive).
From the conditions of continuity of components in the plane we obtain:
i k z ( e i k z z + R p e ( k x ) e i k z z o ) = A χ ε tg ( χ b ) .
Hence the reflection coefficient for the even field TM
R p e ( k ρ ) = e 2 i k z z k z + i χ ε tg ( χ b ) k z i χ ε tg ( χ b ) = e 2 i k z z k z k η p e ( k ρ ) k z + k η p e ( k ρ )
and surface impedance
η p e ( k ρ ) = i χ k ε tg ( χ b ) .

2.3 Odd TE reflection coefficient, surface impedance

We consider the TE (H-polarization) odd field. We are analyzing the problem
in the plane from here x = 0   ρ = y , k ρ = k y . For the considered field TE with components E _ = [ E x , 0,0 ] H _ = [ 0 , H y , H z ]
we adopt its following representation.
E ~ x ( k y ) = I d l 8 π 2 ω μ 0 1 k z e i k y y e ± i k z z + R n ( h ) ( k y ) e i k z z d l a z > z o A sin χ ( z + z o + b ) sin ( χ b ) z o b < z < z o
from Maxwell’s equations for the medium
H ~ x = i ω μ μ o E ~ x z
μ = 1 d l a z > z o μ z o b < z < z o
H ~ x ( k y ) = i I d l 8 π 2 1 k z e i k y y ( ± e ± i k z z + R n ( h ) ( k y ) e i k z z ) i k z A χ μ cos χ ( z + z o + b ) sin ( χ b ) .
Considering the condition of zeroing of the electric field component Ex on the screen z = z 0   b and the continuity of the electromagnetic field components Ex, Hy in the plane z = z 0 , we obtain
( e i k z z + R n h ( k y ) ) e i k z z o = A
i k z ( e + i k z z + R n h ( k y ) e i k z z o ) = A χ μ ctg ( χ b )
Hence, the reflection coefficient for the odd field and TE polarization is described by the relationship.
R n h ( k ρ ) = e 2 i k z z k z k η n h ( k ρ ) k z + k η n h ( k )
and surface admittance
η n h ( k ρ ) = i k χ μ ctg ( χ b ) .

2.4. Reflection coefficient for TE field (H-polarization) even field, surface impedance

We analyze the issue in the plane, hence   x = 0 , k x = 0 , stąd ρ = y ,   k ρ = k y
Figure 3. Dielectric layer with a magnetic screen.
Figure 3. Dielectric layer with a magnetic screen.
Preprints 99982 g003
For the considered TE field with components: E _ = [ E x , 0,0 ] i H _ = [ 0 , H y , H z ] We adopt the following representation.
E ~ x ( k y ) = + I d l 8 π 2 ω μ 0 e i k y y k z e ± i k z z + R p h ( k y ) e i k z z d l a z > z o A cos χ ( z + z o + b ) cos ( χ b ) z o b < z < z o
from Maxwell’s equation for the medium
H ˜ y = i ω μ μ 0 E ˜ x z
μ = 1 z > z 0 μ r z 0 b < z < z 0 .
In numerical calculations, we assume μ=1 for a perfect dielectric
H ~ y ( k y ) = i I d l 8 π 2 e i k y y k z ± e ± i k z z + R p h ( k y ) e i k z z d l a z > z o A χ μ sin χ ( z + z o + b ) cos ( χ b ) z o b < z < z o .
Considering the condition of zeroing of the magnetic field component on the screen z = z 0 b and the continuity of the electromagnetic field components Ex, Hy in the plane z = z 0 , we obtain
e i k z z + R p h ( k y ) = A
i k z ( e i k z z + R p h ( k y ) e i k z z = A χ μ tg ( χ b )
Hence, the reflection coefficient for the even TE field is described by the expression
R p h = e 2 i k z z k z + i χ μ tg ( χ b ) k z i χ μ tg ( χ b ) = e 2 i k z z k z k η p h ( k ρ ) k z + k η p h ( k ρ )
and surface admittance
η p h ( k ρ ) = i k χ μ tg ( χ b ) .

2.5. Reflection coefficient for the TM field (E - polarization), odd field, surface impedance.

Carrying out analogous considerations as in the previous points, we obtain:
R n e = e 2 i k z z k z k η n e ( k ρ ) k z + k η n e ( k ρ )
η n e ( k ρ ) = i k χ ε ctg ( χ b ) .

2.6. Transverse resonance conditions for the dielectric layer on the screen - surface wave poles for the TM and |TE fields

Zeroing of the denominator of the reflection coefficient R(kρ) for a specific polarization and field symmetry
k z + k η e ( h ) ( k ρ ) = 0
gives us the condition for transverse resonance of the dielectric layer, now treated as a waveguide structure. The solution of this condition determines first order pole type singularities in the integral expressions in the field representation. Singularities of this type located on the physical Rieman lobe Im(kz)>0 correspond to surface waves carried by an open waveguide structure. For H polarization we obtain the following resonance conditions for the even part of the field
k z i χ μ tg ( χ b ) =
for the odd part of the field
k z + i χ μ ctg ( χ b ) = 0 .
For E- polarization for the even part of the field
k z i χ ε tg ( χ b ) = 0
for the odd part of the field
k z + i χ ε ctg ( χ b ) = 0 .
The occurrence of this phenomenon should be taken into account when analyzing linear antennas on a dielectric substrate.

2.7. Reflection coefficient, surface impedance for the dielectric layer at a distance b from the perfectly conducting plane (screen)

We analyze the system presented in Figure 4. We consider the TM field (E polarization), an even field with a component E _ = E x 0 , E z H _ = 0 , H y , 0 .
We analyze the problem in the y=0 plane, hence ρ=x , kρ=ky. We assume the following field representation H ˜ y
H ~ y ( k x ) = I d l 8 π 2 e i k x x ( 1 + R e ( k x ) ) e i k z z z > z 0 e i k z z + R e ( k x ) ) e i k z z z 0 < z < 0 A ( k x ) sin ( χ b ) { sin χ ( z + z 0 + b ) f ( k x ) sin χ ( z + z 0 ) } z o b < z < z 0 A ( k x ) f ( k x ) cos k z ( z + z 0 + b + b 1 ) cos ( k z b ) z 0 b b 1 < z < z b .
Using Maxwell’s equations for medium III
E ˜ x ( k x ) = i ω ε ε o H ˜ y z
we designate:
E ~ x ( k x ) = i I d l e i k x x 8 π 2 ω ε 0 i k z ( 1 + R e ( k x ) ) e i k z z z > z 0 i k z [ e i k z z + R e ( k x ) ] e i k z z z 0 < z < 0 A ( k x ) [ χ cos χ ( z + z 0 + b ) f ( k x ) χ cos χ ( z + z 0 ) ] sin ( χ b ) z o b < z < z 0 A ( k x ) f ( k x ) k z sin k z ( z + z 0 + b + b 1 ) cos ( k z b 1 ) z 0 b b 1 < z < z b
where:
f ( x ) = χ ε sin ( χ b ) χ ε ctg ( χ b ) k z tg ( k z b 1 ) 1 .
Taking into account the continuity of the field H ˜ y i E ˜ x in the i plane z = z 0 b and z = z 0 , we obtain an expression for the reflection coefficient
R e ( k ρ ) = e 2 i k z z 0 k z k η e ( k ρ ) k z + k η e ( k ρ )
and surface impedance
η e ( k ρ ) = η n e η p e + η o p η n e + η o p
where : are η p e   a n d   η n e   described by dependencies (18) and (36), respectively
η o p ( k ρ ) = i k z k tg ( k z b 1 ) .
The condition of transverse resonance is expressed by the following relationship:
k z + k η e = η o p ( k ρ ) ( k z + k η n e ) + η m e ( k z + k η n e ) η n e + η o p

2.8. Reflectance and surface impedance of the TE field (H-polarization) even field

We analyze the system shown in Figure 5
We consider the TE field - an even field with components E _ = E x , 0,0   i   H _ = [ 0 , H y , H z ]
We assume the following field representation ρ=y , kρ=ky. We assume the following field representation Ex(ky)
E ~ x k y = I d l 8 π 2 ω μ 0 e k z i k y y 1 + R h k y e i k z z z > z 0 e i k z z + R h ( k y ) ) e i k z z z 0 < z < 0 A k y sin χ z + z 0 + b g k y sin χ z + z 0 sin ( χ b ) z o b < z < z 0 A k y g k y cos k z z + z 0 + b + b 1 cos ( k z b 1 ) z 0 b b 1 < z < z b .
They use Maxwell’s equations (28) to calculate the area H ˜ y ( k y )
H ~ y ( k y x ) = i I d l 8 π 2 e i k x x k z i k z ( 1 + R h ( k y x ) ) e i k z z z > z 0 i k z [ e i k z z + R h e ( k y x ) ] e i k z z z 0 < z < 0 A ( k y ) χ cos χ ( z + z 0 + b ) g ( k y x ) χ cos χ ( z + z 0 ) sin ( χ b ) z o b < z < z 0 A ( k y ) g ( k y ) k z sin k z ( z + z 0 + b + b 1 ) cos ( k z b 1 ) z 0 b b 1 < z < z b .
where:
g ( k y ) = χ μ sin ( χ b ) χ μ ctg ( χ b ) k z tg ( k z b 1 ) 1
Taking into account field continuity H ˜ y , E ˜ x in the plane z = z 0 b i z = z 0 we obtain an expression for the reflectance coefficient
R h ( k ρ ) = e + 2 i k z z 0 k z k η h ( k ρ ) k z + k η h ( k ρ )
and surface admittance of the surface structure
η h ( k ρ ) = η n h η p h + η o p η n h + η o p
The resonance condition is expressed by the following relationship:
k z + k η h = η o p ( k z + k η n h ) + η n h ( k z + k η p h ) η n h + η o p = 0 .

2.8. Reflectance coefficient and surface impedance for TE fields (H - polarization) odd field and TM fields (E - polarization) odd field

Carrying out the considerations as in the previous section, we obtain an odd field for the TE field (H-polarization).
Reflectance coefficient
R h ( k ρ ) = e 2 i k z z k z k η h ( k ρ ) k z + k η h ( k ρ )
surface admittance
η h = η n h η p h + η o n η n h + η o n
where
η o n = i k z k ctg ( k z b 1 ) ,
η n h = i χ k ctg ( χ b ) ,
η p h = i χ k tg ( χ b ) .
Resonance equation
k z + k η h = η o n ( k z + k η n h ) + η n h ( k z + k η p h ) η n p + η o p
and for E polarization
R e ( k ρ ) = e 2 i k z z 0 k z k η e ( k ρ ) k z + k η e ( k ρ )
surface impedance
η e = η n e η p e + η o n η n e + η o n
where: η o n , η p e , η n e are described by dependencies (58), (18) and (26), respectively. And the resonance equation has the form:
k z + k η e = η o n ( k z + k η p e ) + η n e ( k z + k η p e ) η n e + η o n = 0 .

3. Distribution coefficients Generalized impedance zmn

There are two complementary approaches to the analysis of an antenna structure placed in a dielectric homogeneous medium. In the first, the composite antenna system is analyzed by constructing supermodes of the entire structure, in the second, the system is considered as a cascade assembly of flat discrete elements, i.e. boundaries between two dielectrics, periodic metalized planes and dielectric layers [10,11]. The second approach leads to the definition of the scattering matrix, transmission or impedance of the entire structure by cascading appropriate matrices assigned to individual discrete elements of the antenna structure. It is particularly useful in modeling dielectric multilayer antenna walls, where the stored data regarding one planar antenna element can be repeatedly used in the analysis of various antenna systems with modified parameters of other discrete elements of the structure [18,19]. We consider the composite periodic antenna structure shown in Figure 6. The system consists of many layered elements limited at the top and bottom by half-spaces (upper and lower) of homogeneous dielectric media. The third axis (y) of the orthogonal coordinate system (x,y,z) is perpendicular to the drawing plane. The basic components of the structure are: a metalized, infinitely periodic in the (x, y) plane, an antenna surface, a boundary between two dielectrics and a homogeneous dielectric layer.
After determining the Green’s function for the adopted form of the base current, further considerations come down to presenting the generalized impedance in the form of a double integral, where the integrating function will be presented in an explicit form [12]. We obtain the following expression describing the generalized impedance: z m n = 2 π 4 F G k x k x ; k y k y F J o k x k x
F J o k x k x F J o k x e i k x k x x m k x x n i k y k y y m i k y y n + i k x x + i k y y d k y d k x d k y d k x d y d x
The integral described by expression (66) converges slowly, moreover, it contains poles related to the occurrence of a surface wave and requires great caution when calculating in the vicinity of the source [20]. There are two ways to calculate the integral. When calculating in the spatial domain, the spectral variables kx and ky are converted into variables in the polar coordinate systemϕ and K. Using the properties of the spectral Green’s functions, integration over the variable Preprints 99982 i001 is performed analytically. The integral (66) therefore reduces to the form:
x y x 0 y 0 0 J m x G x , y , x 0 , y 0 , K J n x 0 , y 0 d K d y 0 d x 0 d y d x
By integrating the function G over the variable K, we obtain the field at point (x, y) produced by the source placed at point (xo, yo), hence the method is similar to the standard method of moments used to solve linear antennas [13] and can be used in calculations previously learned and developed techniques (e.g. use of Toeplitz matrix symmetry). The main inconvenience of this method, which makes the calculation of antennas with very thin substrates very difficult or even impossible, is the singularity that occurs when the point (x,y) approaches the source in the integration process. Despite this, it is used by some authors [14].
The second approach to calculating the integral requires moving to the spectral domain. The integration over the spatial variables in (66) can then be performed analytically. Occurrence. exponential factors means the Fourier transform of the functions Jm and Jn. Then the integral (66) reduces to
F m k x , k y G _ k x , k y F n k x , k y d k y d k x
where Fm and Fn are the Fourier transforms of the basis and test functions, respectively.
This integral is calculated numerically, but it is the most time-consuming process due to the slow convergence of the integral [15,16]. Slow convergence is caused by the fact that the singularities of the integral in the spatial domain have been removed by converting the domain to a spectral one and have become "scattered" in this domain. Calculating the integral (66) in the spectral domain allows it to be easily extended to the case of an infinite array of antennas [21]. Further analysis is carried out in the spectral domain.We wzorach tych are the reflectance coefficients described by the relations:
  • for single-layer system (18, 26, 33, 35)
  • for a multi-layer system, dependencies (18, 26, 33, 36, 48, 58)
In the above expressions (based on the solution of equations (38÷41), the quantity kz may assume purely imaginary values. This is due to the fact that the dielectric substrate is non-lossy tg δ 1 0 4 But the size ε 0 , μ 0 , ε , μ , x m x n . y m y n , k e f , k n k , L N , b , i ω are real, non-negative parameters that change discreetly within certain ranges.
We will now examine the behavior of the integrand functions in the above formulas. Studies of the integrand functions appearing in the formulas indicate the existence of singular points of the type - branches for:
k 2 k x 2 k y 2 = 0
or
k 2 n 2 k x 2 k y 2 = 0
-
apparent pole for:
k x 2 k e f 2 = 0 (this is a removable singularity)
-
isolated poles (which may occur)
for   1 + ε k 2 k x 2 k y 2 k 2 n 2 k x 2 k y 2 ε k 2 k x 2 k y 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n k x 2 k y 2 = 0
or
for   1 μ k 2 k x 2 k y 2 k 2 n 2 k x 2 k y 2 μ k 2 k x 2 k y 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n k x 2 k y 2 = 0
hence, after transformation, we get that for
k 2 k x 2 k y 2 > 0   or   k 2 n 2 k x 2 k y 2 < 0
poles do not exist.
It is possible for a pole to occur only for
( k n ) 2 k x 2 k y 2 > 0       i       k 2 k x 2 k y 2 < 0 k < k x 2 + k y 2 < k n
In order to determine the value of isolated poles for R p e and R n h , the relationship should be examined:
R p e : 1 + ε k 2 k x 2 k y 2 k 2 n 2 k x 2 k y 2 ε k 2 k x 2 k y 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n k x 2 k y 2
R n h : 1 μ k 2 k x 2 k y 2 k 2 n 2 k x 2 k y 2 μ k 2 k x 2 k y 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n k x 2 k y 2 = 0
From the above it follows that if
a )               k 2 k x 2 k y 2 > 0 b )               k 2 n 2 k x 2 k y 2 < 0   c )   k 2 k x 2 k y 2 < 0         i k 2 n 2 k x 2 k y 2 > 0
If conditions a) and b) are met, the poles do not occur, but if condition c) is met, then the necessary condition for the presence of poles is met.
Let us therefore examine the sufficient condition from the condition mentioned above (74) point. c it follows that:
k < k x 2 + k y 2 < k n
For the situation under consideration, the equations will be written in the form:
R p e : 1 + i ε k x 2 + k y 2 k 2 k 2 n 2 k x 2 k y 2 i ε k x 2 + k y 2 k 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n 2 k x 2 k y 2 = 0
R n h : 1 i μ k x 2 + k y 2 k 2 k 2 n 2 k x 2 k y 2 i μ k x 2 + k y 2 k 2 + k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n 2 k x 2 k y 2 = 0
Hence, after transformations we have:
1 + e 2 i a r c t g ε k x 2 + k y 2 k 2 k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n 2 k x 2 k y 2 = 0
or
1 e 2 i a r c t g μ k x 2 + k y 2 k 2 k 2 n 2 k x 2 k y 2 e 2 i b 1 k 2 n 2 k x 2 k y 2 = 0
For a multilayer system, the equation determining the singularities (poles) of the function takes the form:-
dla   R e tg χ b 1 = 2 χ k z ε χ 2 + ε k z 2 1 + χ 2 ε 2 k z 2 χ 2 + ε 2 k z 2 e 2 b 2 k z k ρ 2 = x e   is   a   sin gular   point :   k 2 < x e < k 2 n 2
- for R h
tg χ b 1 = 2 χ k z μ χ 2 + μ 2 k z 2 1 χ 2 μ r 2 k z 2 χ 2 + μ r 2 k z 2 e 2 b 2 k z k ρ 2 = x h   is   a   sin gular   point :   k 2 < x h < k 2 n 2 k z = k ρ 2 k 2
The remaining patterns remain as for a single-layer system.
The appropriate branches of the elements are specified in individual ranges
in the following way:10
dla   k ρ 2 > k 2 n 2 > k 2 k z = ± i k ρ 2 k 2 χ = ± i k ρ 2 k 2 n 2 R e , h = e 2 z k p 2 k 2 = e 2 b k p 2 k 2 n 2 + η e , h 1 + η e , h e 2 b k p 2 k 2 n 2
η e = ε k p 2 k 2 k p 2 k 2 n 2 ε k p 2 k 2 + k p 2 k 2 n 2 η h = μ k p 2 k 2 k p 2 k 2 n 2 μ k p 2 k 2 + k p 2 k 2 n 2
20
2 0   for   k 2 n 2 > k ρ 2 > k 2       k ρ ( k n , k ) ( k , k n )       k z = ± i k ρ 2 k 2       κ = ± i k 2 n 2 k ρ 2 R e , h = e 2 z k p 2 k 2 = e 2 i b k p 2 k 2 n 2 + η e , h 1 + η e , h e 2 b k p 2 k 2 n 2
η e = i ε k p 2 k 2 k 2 n 2 k p 2 i ε k p 2 k 2 + k 2 n 2 k p 2 η h = i μ k p 2 k 2 k 2 n 2 k p 2 i μ k p 2 k 2 + k 2 n 2 k p 2
3 0   dla   k 2 n 2 > k 2 > k ρ 2       k ρ ( k , + k ) k z = ± i k ρ 2 k 2         κ = ± i k 2 n 2 k ρ 2 R e , h = e 2 z k p 2 k 2 = e 2 i b k p 2 k 2 n 2 + η e , h 1 + η e , h e 2 b k p 2 k 2 n 2
η e = ε k 2 k p 2 k 2 n 2 k p 2 ε k 2 k p 2 + k 2 n 2 k p 2 η h = μ k 2 k p 2 k 2 n 2 k p 2 μ k 2 k p 2 + k 2 n 2 k p 2
Numerical analysis of submitted above expression was made. Results of calculation for typical dielectrics are shown in Figure 7.

Conclusions

Microstrip antennas combine field and peripheral issues and require the use of analytical methods with a high degree of complexity. Therefore, at present, there are no standard methods that can be used in engineering practice. The work is a step towards filling these gaps. It covers all issues related to the analysis of coupling of arbitrarily placed linear antennas on a multilayer dielectric and generalizes the obtained results to antenna arrays.
The problem was reduced to an exact solution of the Hertzian dipole radiation for a layered system. The obtained solution was used to synthesize a solution for the case of finite dimensions of the antenna, using the standard formalism of Green’s functions, representing the field coming from a point source. These functions satisfy the inhomogeneous Helmholtz equation, the Sommerfeld radiation condition and the appropriate boundary conditions. Using Fourier transforms and their properties, an analytical expression for the generalized impedance was determined in the form of a double integral. Due to the fact that there are singular points in the integrand function, the behavior of this function was examined. The conditions for the generation of surface waves were determined. Surface waves in microstrip antennas can be eliminated by selection of base type. It can be done by minimization of b / ε r proportion. Even three-layer dielectric system causes minimization of this proportion but it requires high-guide work. Condition of surface waves minimization is at variance with the antenna bandwidth requirements.

References

  1. Munson R.: „Conformal microstrip antennas and microstrip phased arrays”, IEEE Transactions on Antennas and Propagation Vol. 22 No.11 November 1974. [CrossRef]
  2. Minatti, G.; Faenzi, M.; Martini, E.; Caminita, F.; de Vita, P.; Gonz’alez-Ovejero, D.; Sabbadini, M.; Maci, S. “Modulated metasurface antennas for space: Synthesis, analysis and realizations”. IEEE Transactions on Antennas and Propagation 2014, 63, 1288–1300.
  3. Sharawi S. M. “Printed Multi-Band MIMO Antenna Systems and Their Performance Metrics” IEEE Antennas and Propagation Magazine 2013, 218 - 232.
  4. Oltman G.H. "Electromagnetically coupled microstrip dipoles", IEEE Transactions on Antennas and Propagation AP-29, No 1, 1981.M. N. Sadiku, “Elements of electromagnetics”. Oxford University Press, 2014.
  5. Malaisamy K., Haridoss, Pandian A., Rima S. “Design of Crossed Dipole Yagi–Uda MIMO Antenna for Radar Applications Internationa” Journal of Antennas and Propagation 2022 Volume12.
  6. Haobo Y, Hongguang Z , Hongwei L, Xiaojie D, and Xi C “Computing Mutual Impedance of Antennas by Spherical Harmonic Transform” IEEE Transactions on antennas and propagation, VOL. 67, NO. 7, JULY 2019 4377.
  7. Aourik S. , Erkik A. , Oukaira A. , Dhaou S. Jamal, Z. Lakhssassi A. „An Advanced Array Configuration Antenna Based on Mutual Coupling Reduction” Electronics 2023, 12(7), 1707;
  8. Richmond J.H., Otto D.V.: „Digital computer solutions of the rigorous equations of electromagnetics”, IEEE Transactions on Antennas and Propagation Vol.17, 1969.
  9. Liu S., Sun S., and Chew W. C., “A potential-based integral equation method for low-frequency electromagnetic problems,” IEEE Transactions on Antennas and Propagation, vol. 66, no. 3, pp. 1413–1426, 2018.
  10. Gao, X.; Han, X.; Cao, W.-P.; Li, H.O.; Ma, H.F.; Cui, T.J. “Ultrawideband and high-efficiency linear polarization converter based on double v-shaped metasurface”. IEEE Transactions on Antennas and Propagation. 2015, 63, 3522–3530. [Google Scholar] [CrossRef] . [CrossRef]
  11. Cao Y. F. ; Zhang X. Y. ; Xue Q. “Compact Shared-Aperture Dual-Band Dual-Polarized Array Using Filtering Slot Antenna and Dual-Function Metasurface” IEEE Transactions on Antennas and Propagation, Publication Year: 2022, Page(s):1120 – 1131 04.
  12. Wnuk M. “Multilayer dielectric periodic antenna structure in a cascade view” Appl. Sci. 2022, Volume 12, Issue 9, pp.1-14.
  13. Liu, Q.S.; Sun, S.; Chew, W.C. “A potential-based integral equation method for low-frequency electromagnetic problems”. IEEE Transactions on Antennas and Propagation. 2018, 66, 1413–1426. [Google Scholar] [CrossRef. [CrossRef]
  14. Yaghjian, A.D. “Reflections on Maxwell’s treatise.” Prog. Electromagnetic. Res. 2014, 149, 217–249.
  15. Khaleghi S. S. M., Moradi G. , Jafargholi R. S. S. A. “Microstrip Line Impedance Matching Using ENZ Metamaterials”, Design, and Application Transactions on antennas and propagation, Vol. 67, No. 4, 2019 2243.
  16. Radavaram S., Pour M. , “Wideband Radiation Reconfigurable Microstrip Patch Antenna Loaded With Two Inverted” IEEE Transactions on Antennas and Propagation, Vol. 67, No. 3, 2019 1501-1505.
  17. Balanis, C.A. “Antenna Theory: Analysis and Design”; John Wiley & Sons: Hoboken, NJ, USA, 2016.
  18. Kong, J.A. “Electromagnetic Wave Theory”; EMW Publishing: Cambridge, MA, USA, 2008.
  19. Sadiku, M.N. “Elements of Electromagnetics”; Oxford University Press: Oxford, UK, 2014.
  20. Bender, C.M.; Orszag, S.A. “Advanced Methods for Scientists and Engineers I: Asymptotic Methods and Perturbation Theory”; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013.
  21. Oppenheim, A.V.; Schafer, R.W. “Electromagnetic Field Theory Discrete-Time Signal Processing”; Pearson Education: London, UK, 2014.
Figure 1. Surface impedance.
Figure 1. Surface impedance.
Preprints 99982 g001
Figure 4. Multilayer dielectric structure with an electric shield.
Figure 4. Multilayer dielectric structure with an electric shield.
Preprints 99982 g004
Figure 5. Multilayer dielectric structure with a magnetic screen.
Figure 5. Multilayer dielectric structure with a magnetic screen.
Preprints 99982 g005
Figure 6. Plane periodic antenna structure infinite in the (x,y) plane, multi-layer structure.
Figure 6. Plane periodic antenna structure infinite in the (x,y) plane, multi-layer structure.
Preprints 99982 g006
Figure 7. Distribution of singular points with thickness of dielectric layer.
Figure 7. Distribution of singular points with thickness of dielectric layer.
Preprints 99982 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated