Next Article in Journal
A Cost Analysis of Internet of Things Sensor Data Storage on Blockchain via Smart Contracts
Next Article in Special Issue
Low-Profile and Closely Spaced Four-Element MIMO Antenna for Wireless Body Area Networks
Previous Article in Journal
Marching On-In-Time Unstructured PEEC Method for Electrically Large Structures with Conductive, Dielectric, and Magnetic Media
Previous Article in Special Issue
A Compact Semi-Circular and Arc-Shaped Slot Antenna for Heterogeneous RF Front-Ends
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complex Bianisotropy Effect on the Propagation Constant of a Shielded Multilayered Coplanar Waveguide Using Improved Full Generalized Exponential Matrix Technique

1
Department of Electrical Engineering, University 20 Aout 1955-Skikda, Skikda 21000, Algeria
2
Department of Electronics, University of Ferhat Abbas, Sétif -1-, Sétif 19000, Algeria
3
School of Electrical Engineering and Computer Science, University of Bradford, Bradford BD7 1DP, UK
4
Instituto de Telecomunicações, Campus Universitário de Santiago, 3810-193 Aveiro, Portugal
5
Faculty of Computing, Engineering and Science, University of South Wales, Pontypridd CF37 1DL, UK
6
Department of Computer and Informatics, Faculty of Electrical Engineering and Informatics, Technical University of Kosice, Letná 1/9, 040 01 Košice, Slovakia
7
Department of Electronics, University Mentouri, Constantine 1, Constantine 25000, Algeria
*
Author to whom correspondence should be addressed.
Electronics 2020, 9(2), 243; https://doi.org/10.3390/electronics9020243
Submission received: 23 December 2019 / Revised: 14 January 2020 / Accepted: 28 January 2020 / Published: 2 February 2020
(This article belongs to the Special Issue Recent Technical Developments in Energy-Efficient 5G Mobile Cells)

Abstract

:
A theoretical study of the electromagnetic propagation in a complex medium suspended multilayer coplanar waveguide (CPW) is presented. The study is based on the generalized exponential matrix technique (GEMT) combined with Galerkin’s spectral method of moments applied to a CPW printed on a bianisotropic medium. The analytical formulation is based on a Full-GEMT, a method that avoids usual procedures of heavy and tedious mathematical expressions in the development of calculations and uses matrix-based mathematical expressions instead. These particularities are exploited to develop a mathematical model for the characterization of wave propagation in a three-layer shielded suspended CPW structure. This study is based on the development of mathematical formulations in full compact matrix-based expressions resulting in Green’s functions in a matrix form. The implemented method incorporates a new accelerating procedure developed in the GEMT which provides an initial value used to speed up searching for the exact solution in the principal computation code. This helped us to obtain accurate solutions with tolerable computing time. Good agreements have been achieved with the literature in terms of accuracy and rapid convergence. The results for different cases of bianisotropy have been investigated, and particularly, the effect on the dispersion characteristics is presented and compared with the isotropic case.

1. Introduction

To establish strong foundations for the development of modern and mass-market applications in the field of telecommunications, microwave designers have to develop further efficient devices, which aim to meet the specific needs of modern telecommunication systems, particularly the 5G technology. In addition to providing far better levels of reliability and performances by offering new services, the new technology should be fully consistent with traditional services including 2G, 3G, 4G, Wi-Fi, and other relevant wireless systems, which requires suitable and high performance microwave devices in terms of miniaturization and ultra-wideband characteristics that are proven to be the most challenging issues of all time.
An important class of existing microwave devices is that exploiting the particular properties of bianisotropic media [1] for the development of special and innovative devices that may respond to the needs of modern technologies [2]. In this class of promising materials, we may mention, for example, non-reciprocal, gyrotropic, ferrites, chirals, metamaterials, metasurfaces, etc. [1,3]. Over the last three decades, the electromagnetism of bianisotropic media have gained a great deal of interest from scientists and researchers within the frame of artificial media with new and exciting properties [4]. However, practical exploitation of these phenomena did not develop on a large scale; many physics and engineering problems needed to be solved. Recently, as the science of materials has tremendously advanced, the concept of bianisotropic media has substantially reemerged as a field of importance in microwaves and optics technology [5,6,7].
The particular properties of bianisotropic media arise from a coupling between the electric and magnetic fields that can be explicitly described by general constitutive relations. Due to their diversity, they have found many potential applications from microwaves to optical frequencies such as polarization transformers, directional couplers, antenna and transmission line substrates, antenna radomes, radar systems, chiral waveguides and others [8,9,10,11,12,13].
The electromagnetic properties of bianisotropic media should be analyzed to perceive their exotic characteristics. Several studies have been conducted to characterize the electromagnetic behavior of bianisotropic media [14,15,16,17,18], ferrites [19], metamaterials [20], chiral [21,22], nonreciprocal [23] for simple and complex dielectric based microwave planar structures using numerical and analytical methods [8,14,17,20,24,25,26,27,28,29,30,31,32,33,34,35]. In [8] and [35], the method of lines is used to analyze planar transmission lines with conductor losses and to analyze integrated optical waveguide structures, respectively. In [17] and [20], the transmission line matrix (TLM) method is used for modeling dispersive chiral media and the analysis of dispersion in metamaterials. In [25], a fast computation of planar microstrip lines using the generalized equivalent circuit method of moments is presented. The finite difference technique and the iterated moment method are utilized for the analysis of dielectric [33] and optical waveguides [34]. In [27], a complex image method based on genetic algorithm (GA) is proposed to calculate the Green’s functions of a coplanar waveguide structure. Anisotropic based multilayers, microstrip and waveguide structures are treated in [28,29,31,32]. Recently, in [14] and [24], Karma et al. studied microstrip transmission lines with anisotropic and uniaxial anisotropic substrates using the discrete mode matching method.
To extract the effective constitutive parameters of bianisotropic materials, various techniques such as stepwise method, S-parameters method, resonator method, coaxial probe method, free-space characterization method, rectangular waveguide measurements and recursive algorithms have been employed [36,37,38,39].
By knowing the intrinsic physical properties of complex media, designers can predict the response of microwave components for the development of inventive devices. However, the complexity of mathematical modeling of bianisotropic media is a real challenge in the characterization of microwave components. This has recently become a central area of research in microwaves and optics.
This paper introduces a mathematical modeling of complex media characterized by full 3 × 3 bianisotropic tensors of permittivity, permeability and magneto-electric parameters. The objective of this work is to sufficiently develop predictive mathematical models to judiciously characterize the propagation of electromagnetic waves in a suspended shielded bianisotropic three-layer CPW structure using a Full-GEM technique [28,40,41]. Three primary considerations have been exploited to accomplish new and promising results. These include: the consideration of the most general reciprocal and non-reciprocal chiral and achiral (Tellegen) complex media using the Full-GEMT, estimation of the effective permittivity constant to be used as an initial value to search for the exact solution, and the use of Muller’s method for the extraction of the complex solution of the associated propagation constant in both chiral and achiral media.
The efficient spectral Galerkin-based method of moments (SGMoM) is extensively used to analyze microwave planar structures [18,22,25,27,29,30,41,42,43,44,45,46,47,48]. In our recent work [42], we presented an analytical modeling of a shielded microstrip line based on an anisotropic medium with full 3 × 3 permittivity and permeability tensors. To reduce calculations, some conditions on the permittivity and permeability constants were considered ensuring the decoupling of the TE and TM modes. On the other hand, the computational time drastically increases with the required accuracy because of the slow convergence of integrals and series summations making this approach time intensive [43]. For accelerated convergence and efficient computation, several techniques have been used [44,45,46,47]. The time intensive part in the SGMoM is the evaluation of the matrix elements and the determinant calculation, as the matrix size is large in most cases. In [48], the SGMoM calculation of the propagation constant (β) for a shielded microstrip line is accelerated using asymptotic expansions for the Bessel’s and the Green’s functions with the aid of super convergent series in the approximation of the summation of the leading terms. In [49], for the multilayered shielded microstrip analysis, series summation calculation are accelerated using the Levin’s transformation in the spectral method. In [28], the Green’s function-based volume integral equation computation is accelerated using the fast Fourier transform technique. In [44,45], the impedance matrix elements based on Sommerfeld-type infinite double integral Green functions evaluation is accelerated by converting the infinite double integral of the impedance-matrix elements into a finite one-dimensional integral by using the asymptotic Green’s functions and triangular basis functions with edge condition. In [47], integral-equation formulation in the SGMoM is accelerated by extracting suitable half-space parts of the kernels which leads to an exponentially decaying integrand functions. The integrals of the extracted parts are expressed as combinations of proper integrals and fast converging improper integrals. In [50], an efficient quasi-static analysis is presented, which can be used for speeding up full-wave SGMoM computations as well. Accelerated versions of full-wave spectral domain approach are also reported in [51,52,53,54].
In this paper, we propose a novel approach for the numerical acceleration of the SGMoM for the analysis of bianisotropic medium-based microstrip structures by accelerating and fixing problems of the convergence of the series summation in the elements of the Galerkin’s matrix based on Green’s functions.
The herein considered complex medium is bianisotropic with non-zero magneto-electric tensors (ξ ≠ 0 and η ≠ 0 at the same time). Our previous study, presented in [41], did not treat the case of bianisotropic media with both magneto-electric tensors; only one tensor was considered non zero. The relative simplicity of this case of medium does not require a more complex resolution technique or longer calculation time. Note that this technique failed to provide accurate solutions for general complex bianisotropic media due to the round-off errors and the highly oscillatory fields behavior. A new procedure is implemented to improve the technique and expand it to support the general case of bianisotropy for a CPW structure. Due to the complexity of the considered bianisotropic medium, the resolution method has required a more efficient technique to overcome the drawbacks in terms of non-convergence or considerable calculation time for a tolerable accuracy. An improvement is achieved by introducing an intermediate calculation procedure based on the GEMT to retrieve an approximated initial value of the relative effective permittivity of the three layer-structure as a function of ε r ,   μ r ,   ξ and η : the bianisotropic layer constitutive parameters. This value is used for searching for the exact solutions of the normalized complex constant of propagation ((β/κ0)2 and (α/κ0)2).

2. Exponential Matrix Technique Formulation

The general CPW geometry and the appropriate coordinate system with the z-axis as the direction of propagation are shown in Figure 1. The considered structure is based on a complex bianisotropic medium (region 1) characterized by full 3 × 3-magneto-electric tensors expressing the cross coupling between electric and magnetic fields.
Bianisotropic materials, in their general form, are characterized by the following constitutive relations [30,41,45].
D = ε 0 [ ε ] E + ε 0 μ 0 [ ξ ] H B = μ 0 [ μ ] H + ε 0 μ 0 [ η ] E
The tensors of the relative permittivity [ ε ] , relative permeability [ μ ] and magneto-electric elements [ ξ ] and [ η ] are represented in the Cartesian coordinate system as follows:
ψ = [ ψ x x ψ x y ψ x z ψ y x ψ y y ψ y z ψ z x ψ z y ψ z z ] ,
where ψ stands for [ ε ] , [ μ ] , [ ξ ] , or [ η ] .
Starting from Maxwell’s equations and using the GEMT in the spectral domain, we come to four coupled first-order differential equations for the transverse electromagnetic field components as functions of their derivatives [40,41] given in the Fourier domain:
[ f ˜ ( i ) ( α , β , z ) ] z = [ P ( i ) ] 4 × 4 [ f ˜ ( i ) ( α , β , z ) ] ,
α and β are the Fourier variables corresponding to the space domain wavenumbers κ x and κ y ,
with
[ f ˜ ( i ) ( α , β , z ) ] = [ E ˜ x ( i ) ( α , β , z ) E ˜ y ( i ) ( α , β , z ) H ˜ x ( i ) ( α , β , z ) H ˜ y ( i ) ( α , β , z ) ] ,
and
[ P ] 4 × 4 = j κ 0   { [ M ] + [ N ] [ Q ] [ R ] } ,
where
[ M ] = [ η y x η y y Ζ 0 μ y x Ζ 0 μ y y η x x η x y Ζ 0 μ x x Ζ 0 μ x y Y 0 ε y x Y 0 ε y y ξ y x ξ y y Y 0 ε x x Y 0 ε x y ξ x x ξ x y ] = [ [ η T ] Ζ 0 [ μ T ] Y 0 [ ε T ] [ ξ T ] ] ,
[ N ] = [ ( η y z + κ x n ) Ζ 0 μ y z ( η x z κ y n ) Ζ 0 μ x z Y 0 ε y z ( ξ y z κ x n ) Y 0 ε x z ( ξ x z + κ y n ) ] ,
[ Q ] = 1 ( ε z z μ z z ξ z z η z z ) [ Ζ 0 μ z z ξ z z η z z Y 0 ε z z ] ,
[ R ] = [ Y 0 ε x z Y 0 ε y z ( ξ x z κ y n ) ( ξ y z + κ x n ) ( η x z + κ y n ) ( η y z κ x n ) Ζ 0 μ x z Ζ 0 μ y z ] ,
κ x n = κ x κ 0 ,
  κ y n = κ y κ 0 ,
Ζ 0 = 1 Y 0 = μ 0 ε 0 ,
with
[ η T ] = [ η y x η y y η x x η x y ] ,   [ ε T ] = [ ε y x ε y y ε x x ε x y ] ,   [ ξ T ] = [ ξ y x ξ y y ξ x x ξ x y ] ,   [ μ T ] = [ μ y x μ y y μ x x μ x y ] .
where κ 0 is the free space wavenumber and ω is the angular frequency.
This study is essentially based on the development of mathematical formulations in compact matrix-based forms; this is deemed as a promising approach, since it avoids excessive and complex calculation developments. This can dramatically reduce the complexity of wave propagation modeling in complex media.
The matrix [ P ] (Equation (5)) is the first foundation for this technique associated with the studied bianisotropic-medium based CPW structure. Its elements are given as functions of the constitutive tensors elements. In previous works [30,45], for particular cases of media calculations were explicitly developed, which is not obvious with heavy mathematical calculations case that characterize the herein studied complex bianisotropic structure.
Equation (3) admits a general solution of the form:
[ E x ( z ) E y ( z ) H x ( z ) H y ( z ) ] = T ( κ x , κ y , z )   [ E x ( 0 ) E y ( 0 ) H x ( 0 ) H y ( 0 ) ] ,
with
T ( κ x , κ y , z ) = exp ( [ P ] z ) ,
and
f ˜ ( i ) ( α , β , z ( i ) ) = T ( κ x , κ y , z ) f ˜ ( i ) ( α , β , z ( i 1 ) ) .
The 4 × 4 transfer matrix T ( κ x , κ y ; z ) is calculated in the formulation of the GEMT by means of the Cayley Hamilton theorem for the determination of the complex function roots [40]. It is expressed in the following polynomial form:
T ( z ) = a 0 [ I ] + a 1 [ P ] + a 2 [ P ] 2 + a 3 [ P ] 3 ,
where aj are scalar expansion coefficients, determined by solving the Vandermode linear algebraic system [40], and [ I ] is a 4 × 4 identity matrix. The transfer matrix is easily obtained by multiplying the different transfer matrices related to the different layers of the structure, this constitutes the main advantage of this new technique. By imposing the appropriate boundary conditions between the heterogeneous medium layers, the appropriate Green’s tensor which models the CPW structure is derived. Details can be found in [41]. This technique exhibits a compact matrix form with the advantage of being easily inserted in the calculation code.

2.1. Implementation of the Acceleration Procedure

To overcome the drawbacks of the resolution method used in [41], when applied for a general complex bianisotropic medium, the resolution method has to be improved in terms of convergence and accuracy for a tolerable computing time. A new procedure is introduced in the calculation technique. This latter is based on the GEMT technique (detailed calculations can be found in [40]) used to retrieve an approximated value for the effective relative permittivity of the whole inhomogeneous structure to be used as an initial value to search for the exact solution of the propagation constant. This value is evaluated for each frequency point by extracting the eigenvalues of matrix [ P ] by resolving Equation (12) for κx = 0 and κy = 0. By applying this procedure, analytical expressions of the effective relative permittivity maybe obtained as functions of the constitutive parameters of the bianisotropic layer ε r ,   μ r ,   ξ and η . The application of this procedure allows the acceleration of the Matlab® [55] calculation code and provides a better solution accuracy.
The expansion coefficients a i   ( i = 0 ,   1 ,   2 ,   3 ) in Equation (10) are determined by solving the Vandermode linear algebraic system:
[ 1 λ 0 λ 0 2 λ 0 3 1 λ 1 λ 1 2 λ 1 3 1 λ 2 λ 2 2 λ 2 3 1 λ 3 λ 3 2 λ 3 3 ] [ a 0 ( z ) a 1 ( z ) a 2 ( z ) a 3 ( z ) ] = [ exp ( λ 0 z ) exp ( λ 1 z ) exp ( λ 2 z ) exp ( λ 3 z ) ] ,
λi ( i = 0 ,   1 ,   2 ,   3 ) : are eigenvalues of [ P ] which correspond to propagating waves [56] defined by:
det ( λ [ I ] [ P ] ) = λ 4 + α 1 λ 3 + α 2 λ 2 + α 3 λ + α 4 = 0 .
The coefficients α i   ( i = 0 ,   1 ,   2 ,   3 ) are given in terms of the matrix [ P ] , explicit expressions may be found in [40].

2.2. Derivation of the Initial Value Expression of the Effective Relative Permittivity

a. Isotropic case
The initial effective permittivity value expression is calculated in terms of the medium constitutive parameters ε r ,   μ r ,   ξ and η using the total transfer matrix (Equation (10)) for κx = 0 and κy = 0, which permits the extraction of an approximated value. As an example of calculations, we present the derived analytical expressions of the initial effective permittivity of the isotropic and some bianisotropic cases. For the isotropic case, the derived matrix [ P ] is given by:
[ P ] = [ 0 0 0 μ r Z 0 κ 0 0 0 μ r Z 0 κ 0 0 0 ε r κ 0 / Z 0 0 0 ε r κ 0 / Z 0 0 0 0 ] .
In this case, the normalized eigenvalues of matrix [ P ] , with respect to κ 0 , are:
λ n 0 = ε r μ r ,   λ n 1 = ε r μ r ,   λ n 2 = ε r μ r   and   λ n 3 = ε r μ r ;
functions of the relative permittivity and permeability of the isotropic medium. For the fundamental propagating mode, the numerical maximal value is taken as the initial value.
ε r e f f 0 = λ n 0 2 = ε r μ r
b. Uniaxial anisotropic case
In this example, we consider the following special case of biaxial anisotropy
[ ε ] = [ ε x 0 0 0 ε y 0 0 0 ε z ] ,   [ μ ] = [ μ x 0 0 0 μ y 0 0 0 μ z ] ,   [ ξ ] = 0 ,   [ η ] = 0 ,
the derived matrix [ P ] is:
[ P ] = [ 0 0 0 μ y Z 0 κ 0 0 0 μ x Z 0 κ 0 0 0 ε y κ 0 / Z 0 0 0 ε x κ 0 / Z 0 0 0 0 ] ,
where the four normalized eigenvalues are found to be
λ n 0 = μ x μ z ( ε y μ z 1 ) ,   λ n 1 = μ x μ z ( ε y μ z 1 ) ,   λ n 2 = ε x ε z ( ε z μ y 1 ) ,   λ n 3 = ε x ε z ( ε z μ y 1 )
and
ε r e f f 0 = max ( λ n i 2 ) = max ( μ x μ z ( ε y μ z 1 ) ,   ε x ε z ( ε z μ y 1 ) )
c. Diagonal bianisotropy case
As an example, we take the following case:
[ ε ] = ε r [ I ] ,   [ μ ] = μ r [ I ] ,   [ ξ ] = [ ξ x x 0 0 0 0 0 0 0 0 ] ,   [ η ] = [ η x x 0 0 0 0 0 0 0 0 ] ,
where [ I ] is a 3 × 3 identity matrix. The derived [ P ] is:
[ P ] = [ 0 0 0 μ r Z 0 κ 0 η x x κ 0 0 μ r Z 0 κ 0 0 0 ε r κ 0 / Z 0 0 0 ε r κ 0 / Z 0 0 ξ x x κ 0 0 ] ,
and the four normalized eigenvalues of [ P ] are:
λ n 0 = ε r μ r + ε r μ r ξ x x η x x ,   λ n 1 = ε r μ r + ε r μ r ξ x x η x x λ n 2 = ε r μ r ε r μ r ξ x x η x x ,   λ n 3 = ε r μ r ε r μ r ξ x x η x x .
The initial value depends on the constitutive parameters ε r , μ r , ξ x x and η x x .
ε r e f f 0 = max ( λ n i 2 ) = ε r μ r + ε r μ r ξ x x η x x
d. Gyrotropic bianisotropy case
Considering the following medium case:
[ ε ] = ε r [ I ] ,   [ μ ] = μ r [ I ] ,   [ ξ ] = [ 0 0 ξ x z 0 0 0 ξ z x 0 0 ] ,   [ η ] = [ 0 0 η x z 0 0 0 η z x 0 0 ] .
The derived corresponding [ P ] is:
[ P ] = [ 0 0 0 μ r Z 0 κ 0 0 0 κ 0 Z 0 ε r ( ε r μ r ξ z x η x z ) 0 0 ε r κ 0 / Z 0 0 0 κ 0 μ r Z 0 ( ε r μ r ξ x z η z x ) 0 0 0 ] ,
which gives as solutions:
λ n 0 = ε r μ r ξ z x η x z ,   λ n 1 = ε r μ r ξ z x η x z λ n 2 = ε r μ r ξ x z η z x ,   λ n 3 = ε r μ r ξ x z η z x ,
and
ε r e f f 0 = max ( λ n i 2 ) = max ( ε r μ r ξ x z η z x ,   ε r μ r ξ z x η x z ) .
It is medium-case dependent.
Notice that if ξ x z ,   ξ z x ,   η x z and η z x are taken so that ξ x z η z x ξ z x η x z , two different solutions are obtained corresponding to bifurcating modes [21]. This may constitute an independent issue, which is outside the scope of this work. This shows the efficiency of the procedure, without which, the calculation code would diverge to adjacent solutions or give spurious ones [57], since the electromagnetic fields of bianisotropic media are highly oscillatory [56]. The result presented by Equation (15d) shows that the use of this approach has not only made it feasible to get the optimal initial value, in some cases, but also to more accurately infer the appropriate approximation, mainly in the presence of the bifurcating modes phenomenon, in which two neighboring modes are excited.
In order to show the benefits of using the new procedure, two examples of the studied cases of bianisotropy are considered. Case1: ξ z x , 1 = ξ x z , 1 = 0.75 ε r , η z x , 1 = ξ x z , 1 and η z x , 1 = η x z , 1 , and Case2: η x z , 8 = ξ x z , 8 = j ε r , ξ z x , 8 = ξ x z , 8 and η z x , 8 = η x z , 8 . By the introduction of the new procedure, the computing time for Case1, for a frequency point, is reduced by about 33% from 1.937 s to 1.302 s. In Case2, with the aid of the procedure, we get a solution in 1.442 s while without the procedure, the technique failed to find a solution and gave a spurious value instead, after a long execution time. This is due to the oscillating behavior of the series summations of the manipulated complex Galerkin’s matrix.

3. Method of Solution

By applying the boundary conditions, the expressions of the electric and magnetic tangential components are evaluated at the interface air-dielectric in terms of the tangential current densities on the strips j ˜ x and j ˜ y . A matrix of the Green’s tensor elements G i j ( α n , β ) for the CPW structure is achieved. It is arranged in the following system of equations
[ j ˜ x j ˜ y ] = 1 Δ G [ G 22 ( α n , β ) G 12 ( α n , β ) G 21 ( α n , β ) G 11 ( α n , β ) ] [ E ˜ x E ˜ y ] ,
with
Δ G = G 11 G 22 G 12 G 21 .
and α n : the discrete Fourier transform variable with n the Fourier number of terms n = 1, 2, 3, …, N.
For the resolution of the problem, the SGMoM is applied, the spectral electric field components are expanded in terms of trigonometric basis function sets [58,59]. A homogeneous system of linear equations arranged in a compact matrix form is derived [58]:
[ M ( β ) ]   [ a 1   a 2   a M   b 1   b 2 b N ] T = [ 0 ] ,
where
[ M ( β ) ] =   [ M q , p 1 , 1 M q , q 1 , 2 M p , p 2 , 1 M p , q 1 , 1 ] ,
and
M q , p 1 , 1 ( β ) = n 1 Δ G G 22 ( α n , β )   J ˜ x , p J ˜ y , q * ;
M q , q 1 , 2 ( β ) = n 1 Δ G G 12 ( α n , β )   J ˜ y , q J ˜ y , q * ;
M p , p 2 , 1 ( β ) = n 1 Δ G G 12 ( α n , β )   J ˜ x , p J ˜ x , p * ;
M p , q 2 , 2 ( β ) = n 1 Δ G G 11 ( α n , β )   J ˜ y , q J ˜ x , p * .
The system admits nontrivial solutions when det [M(β)] = 0 [41,45,58,60], from which the frequency dependent propagation constant can be determined. For lossy media, a complex constant solution is expected.
Using the new technique, original results for the dispersion characteristics of complex bianisotropic chiral and achiral media are obtained through the ratio (β/κ0)2, discussed and compared with the isotropic case (ξ = η = 0) using the technique in [41].
Due to the great number of possible medium cases, we have restricted our analysis to highlighting the main results of achiral media that have been less addressed in the literature. Accurate solutions of the determinant roots in Equation (18) are obtained within a tolerance of 10−12.

4. Results and Discussions

In order to validate our calculations and test the efficiency of the proposed method, three magnetic anisotropic cases have initially been considered. Numerical results have been computed and compared with available literature [61,62] (Figure 2) and good agreements are observed. On the other hand, a rapid convergence has been achieved with a reduced Fourier number (N = 500) and basis functions (K = 8) against (N = 3000) used by Khodja et al. [61] for the same number of basis functions.
In this study, we consider a suspended three-layer CPW structure implanted on a complex bianisotropic dielectric material (Figure 1) with the following geometrical dimensions a = 10 mm, d1 = 4.5 mm, d2 = 1 mm, d3 = 4.5 mm, w = 1 mm, s = 1 mm, εr = 2.53, μr = 1. Different sub-figures are differentiated by the included legends where only the sign of the constitutive element changes respectively to the previous case in the same figure.
In order to examine the effect of the magneto-electric parameters on the dispersion characteristics, we first start with the examination of the axial bianisotropy effect. The magneto-electric elements, whether they are real, imaginary, positive or negative directly affect the phase constant as well as the attenuation factor.
The obtained results, treat two principal cases of diagonal bianisotropic medium: achiral and chiral. In each case, the magneto-electric pair ( ξ i j ,   η i j ) is considered non-zero, so that the new original results of the achiral medium case can be validated and compared with the chiral case. In addition, in this parametric study, we examined the effects of the gyrotropic elements of the magneto-electric tensors on the complex propagation coefficient. Results are grouped in figures according to the constitutive parameters effects.

4.1. Effect of Diagonal Bianisotropy

For diagonal bianisotropy three cases of magneto-electric elements are considered:
  • ξ i i , 1 = a ε r , η i i = ξ i i ;
  • ξ i i , 2 = a ε r , η i i = ξ i i ;
  • ξ i i , 3 = j a ε r , η i i = ξ i i .
where (i = x, y, z) and a = (−1, −0.75, −0.5, −0.25, 0.25, 0.5, 0.75, 1). Two main cases are distinguished: achiral ( ξ i i = a ε r ) and chiral ( ξ i i = j a ε r ).
In Figure 3, the effect of element ξ y y , for achiral and chiral media cases is presented. An identical effect is observed on the ratio (β/κ0)2 (Figure 3a), with reciprocal effect (curves superposition for ξ = ± a ε r ), for both chiral and achiral cases. These media cases show low losses for achirality (case (i)) with η y y = ξ y y and almost zero losses for chirality with η y y = ξ y y (case (ii)) (Figure 3b). It can be concluded that the effect of a reciprocal chirality is almost the same as for an achiral medium with equal magneto-electric elements. In these cases (ii = yy) propagating modes are excited in the guiding structure, however, for achiral with η i i = ξ i i and chiral with η i i = ξ i i , no solutions are obtained for ii = xx and ii = zz, hence, the medium does not support any propagating modes.
Figure 4 illustrates the effect of diagonal bianisotropic media for the case achiral with η i i = ξ i i . Unlike the previous case, the fundamental propagating mode is excited for all diagonal magneto-electric elements (ii = xx, yy, zz). However, each of the elements has its own effect. According to Figure 4, for (ii = xx), a non-reciprocity for the achiral case (Figure 4a) is observed on (β/κ0)2. Higher losses are observed for a 0.5 (Figure 4b). The effect of element ξ y y , 2 is presented in Figure 4c and d. These cases are non-reciprocal and exhibit relatively lower losses. The effect on (β/κ0)2 is almost identical for both cases ξ y y , 2 and ξ z z , 2 (Figure 4c,e).

4.2. Effect of Gyrotropic Bianisotropy

For the gyrotropic elements, five achiral cases are considered:
  • ξ i j , 1 = a ε r , η i j = ξ i j , ξ j i = ξ i j ,   η j i = η i j ;
  • ξ i j , 2 = a ε r ,   η i j = ξ i j ,   ξ j i = ξ i j ,   η j i = η i j ;
  • ξ i j , 3 = a ε r ,   η i j = ξ i j ,   ξ j i = ξ i j ,   η j i = η i j ;
  • ξ i j , 4 = a ε r ,   η i j = ξ i j ,   ξ j i = ξ i j ,   η j i = η i j ;
  • ξ i j , 5 = a ε r ,   η i j = ξ i j ,   ξ j i = ξ i j ,   η j i = η i j ;
In Figure 5a–d, five cases were grouped, each differs from the other by a single change in sign of the magneto-electric element. It can be seen that the combination of the constitutive parameters shows non-reciprocity and a distinct effect on (β/κ0)2 and (α/κ0)2 parameters. In Figure 5a,b, for | ξ x y , 1 | and | ξ x y , 2 | 0.25 ε r , we observe a weak effect on the ratio (β/κ0)2 compared to the isotropic case with lower losses for ξ x y , 2 (case (ii)). The sign change between η x y , 1 = η y x , 1 (case (i)) and η x y , 2 = η y x , 2 (case (ii)) keeps the same effect on (β/κ0)2 (Figure 5a). Only the ratio (α/κ0)2 is affected with the appearance of non-reciprocity (Figure 5b). As shown in Figure 5c,d, for the combinations ξ x y , 3 , ξ x y , 4 and ξ x y , 5 , a sign change has only a weak effect on (β/κ0)2 and (α/κ0)2.
In Figure 6 are presented results of the element ξ x z combinations. The remarkable effect is that the ratio (β/κ0)2 is almost constant with respect to frequency and decreases with ξ x z , 1 to reach the unity for | ξ x z , 1 | = 0.75 ε r and zero for | ξ x z , 1 | = ε r (Figure 6a), all with negligible losses (Figure 6b,d). The elements, ξ x z , 2 and ξ x z , 3 have the same effect on both (β/κ0)2 and (α/κ0)2 (Figure 6c,d). In this case, only 3 combinations: ξ x z , 1 , ξ x z , 2 and ξ x z , 3 support propagating modes, with the appearance of the non-reciprocal effect. For the two other cases ξ x z , 4 and ξ x z , 5 , no solutions are obtained. The (β/κ0)2 and (α/κ0)2 variations of both cases (ii) and (iii) are completely different from case (i). The exchange of sign between ξ x z , 3 = ξ z x , 3 and η x z , 2 = η z x , 2 preserved the same effect ((ii) and (iii) curves are superposed), as shown in Figure 6c,d.
Figure 7a,b presents the effect of sign exchange between the magneto-electric element of achiral medium cases. The positive sign ξ y z , 2 = η y z , 2 and ξ y z , 3 = η y z , 3 reveals a non-reciprocity on (β/κ0)2, while it shows a significant effect on (α/κ0)2 and reduced losses compared to isotropic case.

5. Conclusions

In this work, an analytical modeling of a three-layer CPW structure implanted on a complex medium is presented. This study is based on the Full-GEMT developed in a matrix form for the characterization of the bianisotropic-substrate CPW structure. This resulted in compact matrix form expressions of the Green’s tensor. The implemented resolution method includes a new accelerating procedure developed in the GEMT that contributed to accomplishing accurate solutions with improved computing time. The computing time, for one frequency point, has been reduced by 33% for some calculation cases and more for others. A satisfactory calculation convergence is achieved with a significantly reduced Fourier number compared to literature. The presented technique can dramatically reduce the complexity of wave propagation characterization, in highly complex media, in terms of mathematical modelling and computational solution method.
According to our calculations, cases of complex media have been achieved, where the ratio (β/κ0)2 is close to unity such as for cases ξ x z , 1 , ξ x z , 2 . and ξ x z , 3 . This characteristic is vital for the realization of media with permittivity close to unity for better use in the design of radiating antenna structures.
For cases ξ x y , 1 and ξ x y , 2 , where there is only one element that changes sign between one case and another, the results of (β/κ0)2 are similar (reversed cases with similar effects), however, (α/κ0)2 presents a different variation either in form or in magnitude.
It is noted that for the case of achiral medium when ξ = η , ξ i j = ξ j i and η i j = η j i ( ξ x y , 2 ), the medium is reciprocal and the effect is well reversed (non-reciprocal) when ξ i j = ξ j i ( ξ x y , 3 ). For the cases ξ x y , 3 and ξ x y , 5 , one had to have inverted cases, whereas, it is not the case. This can be explained by the properties imposed by the geometry of the studied structure.
An original result that should be taken into consideration is the losses, which show changes with each sign change between the magneto-electric elements. On the other hand, it is found that losses in achiral media are of the same magnitude as those in non-reciprocal chiral media, which must be taken into account. Furthermore, it is worth noting that the transverse elements ξ x z are the most influential on the phase constant in bianisotropic a CPW structure, and this may be attributed to the geometry of the studied structure.
Investigation of some achiral media cases has shown new results such as the notion of achiral media with a relative permittivity approaching unity with reduced losses. This new finding could serve as a valuable concept from which designers of unusual synthetic materials may benefit to enhance the media intrinsic properties for future innovative applications use.
Finally, the technique discussed in this paper may be extended to deal with propagation of bifurcated modes and enclosed multilayer microstrip structures.

Author Contributions

Design and concept, D.S. and C.Z.; methodology, C.Z. and I.E.; investigation, I.O., A.U.; resources, I.E. and J.R.; writing—original draft preparation, C.Z. and D.S.; writing—review and editing, I.E., A.U., H.A. and R.A.-A.; supervision, R.A.-A.; project administration, J.R., formal analysis H.A. and F.B. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement H2020-MSCA-ITN-2016 SECRET-722424. This work is also funded by the FCT/MEC through national funds and when applicable co-financed by the ERDF, under the PT2020 Partnership Agreement under the UID/EEA/50008/2019 project.

Acknowledgments

This work is supported by the European Union’s Horizon 2020 Research and Innovation program under grant agreement H2020-MSCA-ITN-2016-SECRET-722424. This work is also funded by the FCT/MEC through national funds and when applicable co-financed by the ERDF, under the PT2020 Partnership Agreement under the UID/EEA/50008/2019 project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Capolino, F. Theory and Phenomena of Metamaterials; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2009. [Google Scholar]
  2. Asadchy, V.S.; Díaz-Rubio, A.; Tretyakov, S.A. Bianisotropic metasurfaces: Physics and applications. Nanophotonics 2018, 7, 1069–1094. [Google Scholar] [CrossRef]
  3. Ra’di, Y.; Grbic, A. Magnet-free nonreciprocal bianisotropic metasurfaces. Phys. Rev. B 2016, 94, 195432. [Google Scholar] [CrossRef] [Green Version]
  4. Sihvola, A.; Semchenko, I.; Khakhomov, S. View on the history of electromagnetics of metamaterials: Evolution of the congress series of complex media. Photonics Nanostructures-Fundam. Appl. 2014, 12, 279–283. [Google Scholar] [CrossRef] [Green Version]
  5. Novitsky, A.; Shalin, A.S.; Lavrinenko, A.V. Spherically symmetric inhomogeneous bianisotropic media: Wave propagation and light scattering. Phys. Rev. A 2017, 95, 053818. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, N.; Wang, G.P. Effective medium theory with closed-form expressions for bi-anisotropic optical metamaterials. Opt. Express 2019, 27, 23739–23750. [Google Scholar] [CrossRef] [PubMed]
  7. Sihvola, A.; Lindell, I.V. Bianisotropic materials and PEMC. In Theory and Phenomena of Metamaterials; CRC Press: Boca Raton, FL, USA, 2017; pp. 26-1–26-7. [Google Scholar]
  8. Kesari, V.; Keshari, J.P. Hybrid-mode analysis of circular waveguide with chiral dielectric lining for dispersion characteristics for potential application in broadbanding a gyro-traveling-wave tube. J. Electromagn. Waves Appl. 2019, 33, 204–214. [Google Scholar] [CrossRef]
  9. Crowgey, B.R.; Tuncer, O.; Tang, J.; Rothwell, E.J.; Shanker, B.; Kempel, L.C.; Havrilla, M.J. Characterization of Biaxial Anisotropic Material Using a Reduced Aperture Waveguide. IEEE Trans. Instrum. Meas. 2013, 62, 2739–2750. [Google Scholar] [CrossRef]
  10. Wu, B.; Wang, W.; Pacheco, J.; Chen, X.; Lu, J.; Grzegorczyk, T.; Kong, J.A.; Kao, P.; Theophelakes, P.A.; Hogan, M.J. Anisotropic metamaterials as antenna substrate to enhance directivity. Microw. Opt. Technol. Lett. 2006, 48, 680–683. [Google Scholar] [CrossRef]
  11. Bodnar, D.G.; Bassett, H.L. Analysis of an anisotropic dielectric radome. IEEE Trans. Antennas Propag. 1975, 23, 841–846. [Google Scholar] [CrossRef]
  12. Zebiri, C.; Sayad, D.; Elfergani, I.; Iqbal, A.; Mshwat, W.F.; Kosha, J.; Rodriguez, J.; Abd-Alhameed, R. A compact semi-circular and arc-shaped slot antenna for heterogeneous RF front-ends. Electronics 2019, 8, 1123. [Google Scholar] [CrossRef] [Green Version]
  13. Meshram, M.R.; Agrawal, N.K.; Sinha, B.; Misra, P.S. Characterization of M-type barium hexagonal ferrite-based wide band microwave absorber. J. Magn. Magn. Mater. 2004, 271, 207–214. [Google Scholar] [CrossRef]
  14. Kamra, V.; Dreher, A. Efficient analysis of multiple microstrip transmission lines with anisotropic substrates. IEEE Microw. Wirel. Compon. Lett. 2018, 28, 636–638. [Google Scholar] [CrossRef]
  15. Buzov, A.L.; Buzova, M.A.; Klyuev, D.S.; Mishin, D.V.; Neshcheret, A.M. Calculating the Input Impedance of a Microstrip Antenna with a Substrate of a Chiral Metamaterial. J. Commun. Technol. Electron. 2018, 63, 1259–1264. [Google Scholar] [CrossRef]
  16. Klyuev, D.S.; Minkin, M.A.; Mishin, D.V.; Neshcheret, A.M.; Tabakov, D.P. Characteristics of Radiation from a Microstrip Antenna on a Substrate Made of a Chiral Metamaterial. Radiophys. Quantum Electron. 2018, 61, 445–455. [Google Scholar] [CrossRef]
  17. Zhou, Z.; Keller, S.M. The Application of Least-Squares Finite-Element Method to Simulate Wave Propagation in Bianisotropic Media. IEEE Trans. Antennas Propag. 2019, 67, 2574–2582. [Google Scholar] [CrossRef]
  18. Zebiri, C.; Lashab, M.; Benabdelaziz, F. Effect of anisotropic magneto-chirality on the characteristics of a microstrip resonator. IET Microw. Antennas Propag. 2010, 4, 446–452. [Google Scholar] [CrossRef]
  19. Balbastre, J.V.; Nuño, L. Modelling the propagation of electromagnetic waves across complex metamaterials in closed structures. J. Comput. Appl. Math. 2019, 352, 40–49. [Google Scholar] [CrossRef]
  20. Xiong, Y.; Russer, J.A.; Che, W.; Shen, G.; Han, Y.; Russer, P. Dispersion analysis of a fishnet metamaterial based on the rotated transmission-line matrix method. IET Microw. Antennas Propag. 2015, 9, 1345–1353. [Google Scholar] [CrossRef]
  21. Aib, S.; Benabdelaziz, F.; Zebiri, C.; Sayad, D. Propagation in diagonal anisotropic chirowaveguides. Adv. Optoelectron. 2017, 2017, 9524046. [Google Scholar] [CrossRef] [Green Version]
  22. Zebiri, C.; Daoudi, S.; Benabdelaziz, F.; Lashab, M.; Sayad, D.; Ali, N.T.; Abd-Alhameed, R.A. Gyro-chirality effect of bianisotropic substrate on the operational of rectangular microstrip patch antenna. Int. J. Appl. Electromagn. Mech. 2016, 51, 249–260. [Google Scholar] [CrossRef]
  23. Alù, A.; Krishnaswamy, H. Artificial nonreciprocal photonic materials at GHz-to-THz frequencies. MRS Bull. 2018, 43, 436–442. [Google Scholar] [CrossRef]
  24. Kamra, V.; Dreher, A. Discrete mode matching method for the analysis of microstrip lines on uniaxial anisotropic substrates. In Proceedings of the International Applied Computational Electromagnetics Society Symposium ACES, Firenze, Italy, 26–30 March 2017; pp. 1–2. [Google Scholar]
  25. Oueslati, N.; Aguili, T. An Improved MoM-GEC Method for Fast and Accurate Computation of Transmission Planar Structures in Waveguides: Application to Planar Microstrip Lines. Prog. Electromagn. Res. 2016, 48, 9–24. [Google Scholar] [CrossRef] [Green Version]
  26. Ardakani, H.H.; Mehrdadian, A.; Forooraghi, K. Analysis of graphene-based microstrip structures. IEEE Access 2017, 5, 20887–20897. [Google Scholar] [CrossRef]
  27. Han, D.; Lee, C.; Kahng, S. Formulation of the Green’s Functions for Coplanar Waveguide Microwave Devices as Genetic Algorithm-Based Complex Images. J. Electr. Eng. Technol. 2017, 12, 1600–1604. [Google Scholar]
  28. Hu, Y.; Fang, Y.; Wang, D.; Zhan, Q.; Zhang, R.; Liu, Q.H. The scattering of electromagnetic fields from anisotropic objects embedded in anisotropic multilayers. IEEE Trans. Antennas Propag. 2019, 67, 7561–7568. [Google Scholar] [CrossRef]
  29. Daoudi, S.; Benabdelaziz, F.; Zebiri, C.; Sayad, D.; Abdussalam, F.M.; Abd-Alhameed, R. Dispersion characteristics of a gyro-chiro-ferrite shielded multilayered microstrip line using the generalized exponential matrix technique. In Proceedings of the ITA Internet Technologies and Applications, Wrexham, UK, 12–15 September 2017; pp. 293–298. [Google Scholar]
  30. Zebiri, C.; Lashab, M.; Benabdelaziz, F. Rectangular microstrip antenna with uniaxial bianisotropic chiral substrate-superstrate. IET Microw. Antennas Propag. 2011, 5, 17–29. [Google Scholar] [CrossRef]
  31. Kamra, V.; Dreher, A. Analysis of Anisotropic Inhomogeneous Dielectric Waveguides With Discrete Mode Matching Method. In Proceedings of the IMS International Microwave Symposium, Boston, MA, USA, 2–8 June 2019; pp. 24–27. [Google Scholar]
  32. Kamra, V.; Dreher, A. Multilayered Transmission Lines on Quasi-planar Substrates With Anisotropic Medium. Adv. Radio Sci. 2019, 17, 77–82. [Google Scholar] [CrossRef]
  33. Horikis, T.P. Dielectric waveguides of arbitrary cross sectional shape. Appl. Math. Model. 2013, 37, 5080–5091. [Google Scholar] [CrossRef]
  34. She, S.X. Iterated-moment method for the analysis of optical waveguides of arbitrary cross section. JOSA A 1989, 6, 1031–1037. [Google Scholar] [CrossRef]
  35. Yang, W.D.; Pregla, R. The method of lines for analysis of integrated optical waveguide structures with arbitrary curved interfaces. J. Lightwave Technol. 1996, 14, 879–884. [Google Scholar] [CrossRef]
  36. Hasar, U.C.; Barroso, J.J.; Sabah, C.; Kaya, Y.; Ertugrul, M. Stepwise technique for accurate and unique retrieval of electromagnetic properties of bianisotropic metamaterials. J. Opt. Soc. Am. B 2013, 30, 1058–1068. [Google Scholar] [CrossRef]
  37. Hasar, U.C.; Barroso, J.J.; Bute, M.; Muratoglu, A.; Ertugrul, M. Boundary effects on the determination of electromagnetic properties of bianisotropic metamaterials from scattering parameters. IEEE Trans. Antennas Propag. 2016, 64, 3459–3469. [Google Scholar] [CrossRef]
  38. Hasar, U.C.; Muratoglu, A.; Bute, M.; Barroso, J.J.; Ertugrul, M. Effective Constitutive Parameters Retrieval Method for Bianisotropic Metamaterials Using Waveguide Measurements. IEEE Trans. Microw. Theory Tech. 2017, 65, 1488–1497. [Google Scholar] [CrossRef]
  39. Hasar, U.C.; Buldu, G.; Kaya, Y.; Ozturk, G. Determination of Effective Constitutive Parameters of Inhomogeneous Metamaterials With Bi-anisotropy. IEEE Trans. Microw. Theory Tech. 2018, 66, 3734–3744. [Google Scholar] [CrossRef]
  40. Tsalamengas, J.L. Interaction of electromagnetic waves with general bianisotropic slabs. IEEE Trans. Microw. Theory Tech. 1992, 40, 1870–1878. [Google Scholar] [CrossRef]
  41. Daoudi, S.; Benabdelaziz, F.; Zebiri, C.; Sayad, D. Generalized Exponential Matrix Technique Application for the Evaluation of the Dispersion Characteristics of a Chiro-Ferrite shielded Multilayered Microstrip Line. Prog. Electromagn. Res. M 2017, 61, 1–14. [Google Scholar] [CrossRef] [Green Version]
  42. Sayad, D.; Zebiri, C.; Daoudi, S.; Benabdelaziz, F. Analysis of the Effect of a Gyrotropic Anisotropy on the Phase Constant and Characteristic Impedance of a Shielded Microstrip Line. Adv. Electromagn. 2019, 8, 15–22. [Google Scholar]
  43. Bianconi, G.; Mittra, R. Efficient Numerical Techniques for Analyzing Microstrip Circuits and Antennas Etched on Layered Media via the Characteristic Basis Function Method. In Computational Electromagnetics; Springer: New York, NY, USA, 2014; pp. 111–148. [Google Scholar]
  44. Park, S.-O.; Balanis, C.A. Analytical technique to evaluate the asymptotic part of the impedance matrix of Sommerfeld-type integrals. IEEE Trans. Antennas Propag. 1997, 45, 798–805. [Google Scholar] [CrossRef] [Green Version]
  45. Sayad, D.; Benabdelaziz, F.; Zebiri, C.; Daoudi, S.; Abd-Alhameed, R.A. Spectral domain analysis of gyrotropic anisotropy chiral effect on the input impedance of a printed dipole antenna. Prog. Electromagn. Res. M 2016, 51, 1–8. [Google Scholar] [CrossRef] [Green Version]
  46. Jain, S.; Song, J. Accelerated spectral domain approach for shielded microstrip lines by approximating summation with super convergent series. In Proceedings of the Digests of the 14th Biennial IEEE Conference on Electromagnetic Field Computation, Chicago, IL, USA, 9–12 May 2010; p. 1. [Google Scholar]
  47. Lucido, M. A new high-efficient spectral-domain analysis of single and multiple coupled microstrip lines in planarly layered media. IEEE Trans. Microw. Theory Tech. 2012, 60, 2025–2034. [Google Scholar] [CrossRef]
  48. Jain, S.; Song, J.; Kamgaing, T.; Mekonnen, Y.S. Acceleration of spectral domain approach for generalized multilayered shielded microstrip interconnects using two fast convergent series. IEEE Trans. Compon. Packag. Manuf. Technol. 2013, 3, 401–410. [Google Scholar] [CrossRef]
  49. Xu, H.; Jain, S.; Song, J.; Kamgaing, T.; Mekonnen, Y.S. Acceleration of spectral domain immitance approach for generalized multilayered shielded microstrips using the Levin’s transformation. IEEE Antennas Wirel. Propag. Lett. 2014, 14, 92–95. [Google Scholar] [CrossRef]
  50. Medina, F.; Horno, M. Quasianalytical static solution of the boxed microstrip line embedded in a layered medium. IEEE Trans. Microw. Theory Tech. 1992, 40, 1748–1756. [Google Scholar] [CrossRef]
  51. Tsalamengas, J.L.; Fikioris, G. Rapidly converging spectral-domain analysis of rectangularly shielded layered microstrip lines. IEEE Trans. Microw. Theory Tech. 2003, 51, 1729–1734. [Google Scholar] [CrossRef]
  52. Railton, C.J.; McGeehan, J.P. A rigorous and computationally efficient analysis of microstrip for use as an electro-optic modulator. IEEE Trans. Microw. Theory Tech. 1989, 37, 1099–1104. [Google Scholar] [CrossRef] [Green Version]
  53. Tao, J.-W.; Angenieux, G.; Flechet, B. Full-wave description of propagation and losses in quasi-planar transmission lines by quasi-analytical solution. IEEE Trans. Microw. Theory Tech. 1994, 42, 1246–1253. [Google Scholar] [CrossRef]
  54. Cano, G.; Medina, F.; Horno, M. Efficient spectral domain analysis of generalized microstrip lines in stratified media including thin, anisotropic and lossy substrates. IEEE Trans. Microw. Theory Tech. 1992, 40, 217–227. [Google Scholar] [CrossRef]
  55. MATLAB 9.4 (R2018a), The MathWorks, Inc.: Natick, MA, USA, 2018.
  56. Yang, H.Y. A numerical method of evaluating electromagnetic fields in a generalized anisotropic medium. IEEE Trans. Microw. Theory Tech. 1995, 43, 1626–1628. [Google Scholar] [CrossRef]
  57. Mariotte, F.; Pelet, P.; Engheta, N. A review of recent study of guided waves in chiral media. Prog. Electromagn. Res. 1994, 9, 311–350. [Google Scholar]
  58. Krowne, C.M. Electromagnetic propagation and field behavior in highly anisotropic media. Adv. Imaging Electron Phys. 1995, 92, 79–214. [Google Scholar]
  59. Nguyen, C. Analysis Methods for RF, Microwave, and Millimeter-Wave Planar Transmission Line Structures; John Wiley & Sons: New York, NY, USA, 2003. [Google Scholar]
  60. Mirshekar-Syhakal, D. Spectral Domain Method for Microwave Integrated Circuits; Wiley: New York, NY, USA, 1990. [Google Scholar]
  61. Maze-Merceur, G.; Tedjini, S.; Bonnefoy, J.L. Analysis of a CPW on electric and magnetic biaxial substrate. IEEE Trans. Microw. Theory Tech. 1993, 41, 457–461. [Google Scholar] [CrossRef]
  62. Khodja, A.; Yagoub, M.C.E.; Touhami, R.; Baudrand, H. Practical Recurrence Formulation for Composite Substrates: Application to Coplanar Structures with Bi-Anisotropic Dielectrics. In Proceedings of the 18th Mediterranean Microwave Symposium (MMS), Istanbul, Turkey, 31 October–2 November 2018; pp. 341–344. [Google Scholar]
Figure 1. Geometry of the shielded suspended 3-layer CPW structure.
Figure 1. Geometry of the shielded suspended 3-layer CPW structure.
Electronics 09 00243 g001
Figure 2. (β/κ0)2 for the dominant mode of a shielded CPW with magnetic anisotropy, (2a = 3.556 mm, 2w = 2s = 0.7112 mm, d1 = 2.8448 mm, d2 = 0.7112 mm, d3 = 3.556 mm and εr = 3).
Figure 2. (β/κ0)2 for the dominant mode of a shielded CPW with magnetic anisotropy, (2a = 3.556 mm, 2w = 2s = 0.7112 mm, d1 = 2.8448 mm, d2 = 0.7112 mm, d3 = 3.556 mm and εr = 3).
Electronics 09 00243 g002
Figure 3. Effect of diagonal bianisotropy ((i) ξyy real with ηyy = ξyy and (ii) ξyy imaginary with ηyy = −ξyy) on (a): (β/κ0)2 and (b): (α/κ0)2.
Figure 3. Effect of diagonal bianisotropy ((i) ξyy real with ηyy = ξyy and (ii) ξyy imaginary with ηyy = −ξyy) on (a): (β/κ0)2 and (b): (α/κ0)2.
Electronics 09 00243 g003
Figure 4. Effect of diagonal bianisotropy elements, ξxx on: (a) (β/κ0)2 and (b) (α/κ0)2; ξyy on: (c) (β/κ0)2 and (d) (α/κ0)2; ξzz on: (e) (β/κ0)2 and (f) (α/κ0)2 with ηxx = −ξxx.
Figure 4. Effect of diagonal bianisotropy elements, ξxx on: (a) (β/κ0)2 and (b) (α/κ0)2; ξyy on: (c) (β/κ0)2 and (d) (α/κ0)2; ξzz on: (e) (β/κ0)2 and (f) (α/κ0)2 with ηxx = −ξxx.
Electronics 09 00243 g004aElectronics 09 00243 g004b
Figure 5. Effect of gyrotropic achiral bianisotropy for different ξxy, ξyx, ηxy and ηyx combinations on: (a) (β/κ0)2, (b) (α/κ0)2, (c) (β/κ0)2 and (d) (α/κ0)2.
Figure 5. Effect of gyrotropic achiral bianisotropy for different ξxy, ξyx, ηxy and ηyx combinations on: (a) (β/κ0)2, (b) (α/κ0)2, (c) (β/κ0)2 and (d) (α/κ0)2.
Electronics 09 00243 g005
Figure 6. Effect of gyrotropic achiral bianisotropy for different ξxz, ξzx, ηxz and ηzx combinations on: (a) (β/κ0)2, (b) (α/κ0)2, (c) (β/κ0)2 and (d) (α/κ0)2.
Figure 6. Effect of gyrotropic achiral bianisotropy for different ξxz, ξzx, ηxz and ηzx combinations on: (a) (β/κ0)2, (b) (α/κ0)2, (c) (β/κ0)2 and (d) (α/κ0)2.
Electronics 09 00243 g006
Figure 7. Effect of gyrotropic achiral bianisotropy for different ξyz, ξzy, ηyz and ηzy combinations on: (a) (β/κ0)2 and (b) (α/κ0)2.
Figure 7. Effect of gyrotropic achiral bianisotropy for different ξyz, ξzy, ηyz and ηzy combinations on: (a) (β/κ0)2 and (b) (α/κ0)2.
Electronics 09 00243 g007

Share and Cite

MDPI and ACS Style

Sayad, D.; Zebiri, C.; Elfergani, I.; Rodriguez, J.; Abobaker, H.; Ullah, A.; Abd-Alhameed, R.; Otung, I.; Benabdelaziz, F. Complex Bianisotropy Effect on the Propagation Constant of a Shielded Multilayered Coplanar Waveguide Using Improved Full Generalized Exponential Matrix Technique. Electronics 2020, 9, 243. https://doi.org/10.3390/electronics9020243

AMA Style

Sayad D, Zebiri C, Elfergani I, Rodriguez J, Abobaker H, Ullah A, Abd-Alhameed R, Otung I, Benabdelaziz F. Complex Bianisotropy Effect on the Propagation Constant of a Shielded Multilayered Coplanar Waveguide Using Improved Full Generalized Exponential Matrix Technique. Electronics. 2020; 9(2):243. https://doi.org/10.3390/electronics9020243

Chicago/Turabian Style

Sayad, Djamel, Chemseddine Zebiri, Issa Elfergani, Jonathan Rodriguez, Hasan Abobaker, Atta Ullah, Raed Abd-Alhameed, Ifiok Otung, and Fatiha Benabdelaziz. 2020. "Complex Bianisotropy Effect on the Propagation Constant of a Shielded Multilayered Coplanar Waveguide Using Improved Full Generalized Exponential Matrix Technique" Electronics 9, no. 2: 243. https://doi.org/10.3390/electronics9020243

APA Style

Sayad, D., Zebiri, C., Elfergani, I., Rodriguez, J., Abobaker, H., Ullah, A., Abd-Alhameed, R., Otung, I., & Benabdelaziz, F. (2020). Complex Bianisotropy Effect on the Propagation Constant of a Shielded Multilayered Coplanar Waveguide Using Improved Full Generalized Exponential Matrix Technique. Electronics, 9(2), 243. https://doi.org/10.3390/electronics9020243

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop