Applications and Solving Techniques of Propagated Wave in Waveguides Filled with Inhomogeneous Dielectric Materials

This chapter presents techniques to solve problems of propagation along the straight rectangular and circular waveguides with inhomogeneous dielectric materials in the cross section. These techniques are very important to improve the methods that are based on Laplace and Fourier transforms and their inverse transforms also for the discontinuous rectangular and circular profiles in the cross section (and not only for the continuous profiles). The main objective of this chapter is to develop the techniques that enable us to solve problems with inhomogeneous dielectric materials in the cross section of the straight rectangular and circular waveguides. The second objective is to understand the influence of the inhomogeneous dielectric materials on the output fields. The method in this chapter is based on the Laplace and Fourier transforms and their inverse transforms. The proposed techniques together with the methods that are based on Laplace and Fourier transforms and their inverse transforms are important to improve the methods also for the discontinuous rectangular and circular profiles in the cross section. The applications are useful for straight waveguides in the microwave and the millimeter-wave regimes, for the straight hollow waveguide and for infrared field, also in the cases of inhomogeneous dielectric materials in the cross section.


Introduction
The methods of straight waveguides have been proposed in the literature. Review of numerical and approximate methods for the modal analysis of general optical dielectric waveguides with emphasis on recent developments has been published [1]. Examples of important methods have been proposed such as finite-difference method, the integral-equation method, and methods based on series expansion. Full-vectorial matched interface and boundary method for modal analysis of dielectric waveguides has been proposed [2]. The method distinguishes itself with other existing interface methods by avoiding the use of the Taylor series expansion and by introducing the concept of the iterative use of low-order jump conditions. A review of the hollow waveguide and the applications has been presented [3,4]. A review of hollow waveguides, infrared transmitting, and fibers has been presented [5]. Hollow waveguides with metallic and dielectric layers have been proposed to reduce the transmission losses. A hollow waveguide can be made from any flexible or rigid tube, such as glass, plastic or glass, and the inner hollow surface is covered by a metallic layer and a dielectric overlayer. The structure of the layer enables to transmit both the TE and TM polarization with low attenuation [6,7].
Selective suppression of electromagnetic modes in rectangular waveguides has been presented [8] by using distributed wall losses. Analytical model for the corrugated rectangular waveguide has been extended to compute the dispersion and interaction impedance [9].
A Fourier operator method has been used to derive for the first time an exact closed-form eigenvalue equation for the scalar mode propagation constants of a buried rectangular dielectric waveguide [10]. Wave propagation in an inhomogeneous transversely magnetized rectangular waveguide has been studied with the aid of a modified Sturm-Liouville differential equation [11]. A fundamental and accurate technique to compute the propagation constant of waves in a lossy rectangular waveguide has been proposed [12]. This method is based on matching the electric and magnetic fields at the boundary and allowing the wavenumbers to take complex values.
A method that relates to the propagation constant for the bound modes in the dielectric rectangular waveguides has been proposed [13]. An analysis of rectangular folded-waveguide slow-wave structure has been developed using conformal mapping by using Schwarz Christoffel transformation [14]. A simple closed form expression to compute the time-domain reflection coefficient for a transient TE 10 mode wave incident on a dielectric step discontinuity in a rectangular waveguide has been presented [15]. In this paper, an exponential series approximation was provided for efficient computation of the reflected and transmitted field waveforms.
The electromagnetic fields in rectangular conducting waveguides filled with uniaxial anisotropic media have been characterized [16]. In this paper, the electric type dyadic Green's function due to an electric source was derived by using eigenfunctions expansion and the Ohm-Rayleigh method. An improved generalized admittance matrix technique based on mode matching method has been proposed [17]. The generalized scattering matrix of waveguide structure and its discontinuity problems is obtained with relationship equations and reflection coefficients.
A full-vectorial boundary integral equation method for computing guided modes of optical waveguides has been proposed [18]. Method for the propagation constants of fiber waveguides of arbitrary cross section shape has been proposed [19]. In this paper, the proposed techniques used to solve problems of scattering by irregularly shaped dielectric bodies, and in the static limit, for solving the problem of an irregular dielectric or permeable body in an external field.
The rectangular dielectric waveguide technique for the determination of complex permittivity of a wide class of dielectric materials of various thicknesses and cross sections has been described [20]. In this paper, the technique has been presented to determine the dielectric constant of materials. The fields and propagation constants of dielectric waveguides have been determined within the scalar regime [21] by using two-dimensional Fourier series expansions. Propagation of modes in some rectangular waveguides using the finite-difference time-domain method has been proposed [22]. Analysis of rectangular waveguide using finite element method has been proposed [23].
Wave propagation and dielectric permittivity reconstruction in the case of a rectangular waveguide have been studied [24]. According to this paper, we study the electromagnetic wave propagation in a rectangular waveguide filled with an inhomogeneous dielectric material in the longitudinal direction. Light propagation in a cylindrical waveguide with a complex, metallic and dielectric function has been proposed [25]. Advancement of algebraic function approximation in eigenvalue problems of lossless metallic waveguides to infinite dimensions has been investigated [26]. The method of algebraic function approximation in eigenvalue problems of lossless metallic waveguides such as a closed uniform cylindrical waveguides has been proposed [27]. Analysis of longitudinally inhomogeneous waveguides using Taylor's series expansion has been proposed [28]. Analysis of longitudinally inhomogeneous waveguides using the Fourier series expansion has been proposed [29]. The method of external excitation for analysis of arbitrarily-shaped hollow conducting waveguides has been proposed [30].
A circular metallic hollow waveguide with inner dielectric multilayers has been designed by Miyagi and Kawakami [31] with the emphasis on low-loss transmission of the HE 11 mode for the infrared. According to this paper, the transmission losses of the dielectric-coated metal waveguides are drastically reduced when a multiple dielectric layer is formed instead of a single dielectric layer. The simplest and most efficient multilayer structure is a three-dielectriclayer stack deposited on a metal layer.
A transfer matrix function for the analysis of electromagnetic wave propagation along the straight dielectric waveguide with arbitrary profiles has been proposed [32]. According to this paper, the method is based on the Laplace and Fourier transforms, and the inverse Laplace and Fourier transforms. A rigorous approach for the propagation of electromagnetic fields along a straight hollow waveguide with a circular cross section has been proposed [33]. The cross section is made of a metallic layer, and only one dielectric layer upon it. The separation of variables is obtained by using the orthogonal-relations. The longitudinal components of the fields are developed into the Fourier-Bessel series. The transverse components of the fields are expressed as functions of the longitudinal components in the Laplace plane and are obtained by using the inverse Laplace transform by the residue method.
In order to solve more complex problems of coatings in the cross section of the dielectric waveguides, such as rectangular and circular profiles, then it is important to develop in each modal an improved technique to calculate the dielectric profile, the elements of the matrix and its derivatives of the dielectric profile.
The main objective of this chapter is to develop the techniques that enable us to solve problems with inhomogeneous dielectric materials in the cross section of the straight rectangular and circular waveguides. The second objective is to understand the influence of the inhomogeneous dielectric materials on the output fields. Thus, we need to develop the technique and a particular application to calculate the profiles in the cross section. Namely, we need to calculate the dielectric profile, the elements of the matrix and its derivatives of the dielectric profile in the cases of the straight rectangular and circular waveguides. The proposed techniques are important to improve the methods that are based on Laplace and Fourier transforms and their inverse Laplace and Fourier transforms also for the discontinuous rectangular and circular profiles in the cross section (and not only for the continuous profiles).

Formulation of the problem
In this chapter, we present techniques for solving discontinuous problems of dielectric materials in the cross section of the straight waveguide for applications in the microwave and millimeter-wave regimes and in the cases of infrared regime. The proposed techniques are very effective in relation to the conventional methods because they allow the development of expressions in the cross section only according to the specific discontinuous problem. In this way, the mode model method becomes an improved method to solve discontinuous problems in the cross section.
Three examples of inhomogeneous dielectric materials in the cross section of the straight waveguides are shown in Figure 1   Emerging Waveguide Technology circular profile in the cross section of the straight rectangular waveguide. Figure 1(c) shows an example for three dielectric layers and a metallic layer in the cross section of the straight hollow waveguide.
In order to solve these inhomogeneous dielectric materials, we need to calculate the dielectric profile, the elements of the matrix and its derivatives of the dielectric profile. In the next section, we explain the techniques to solve problems with inhomogeneous dielectric materials in the cross section of the straight rectangular waveguide (Figure 1(a) and 1(b)) and also in the cross section of the straight circular waveguide (Figure 1(c)).

The technique to solve inhomogeneous dielectric profiles
The particular application is based on the ω ε function [34]. The ω ε function is used in order to solve discontinuous problems of the rectangular profile, and circular profile in the cross section of the straight waveguide. The ω ε function is defined as where C ε is a constant, and Ð ω ε r ðÞ dr ¼ 1. In the limit ε ! 0, the ω ε function is shown in Figure 2.
The technique that based on ω ε function is very effective to solve complex problems, in relation to the conventional methods, especially when we have a large numbers of dielectric layers and a metal layer, as shown in Figure 1(c). We will demonstrate how to use with the proposed technique for all the cases that are shown in the examples in Figure 1  The dielectric profile gx ; y ðÞ is given according to ε x; y ðÞ ¼ ε 0 1 þ gx ; y ðÞ ðÞ .AccordingtoFigure 3(a) and 3(b) and for gx ; y ðÞ ¼ g 0 ,weobtain Ày 12 exp Àjk x x þ k y y ÀÁ ÀÁ dy ) : (1) If y 11 and y 12 are not functions of x, then the dielectric profile is given by cos k x x ðÞ dx ð y 12 y 11 cos k y y ÀÁ dy: The derivative of the dielectric profile in the case of y 11 and y 12 are functions of x, is given by where g x x; y ðÞ ¼ 1=e x;y ðÞ ðÞ de x; y ðÞ =dx ðÞ , e x; y ðÞ ¼ e 0 1 þ gx ;y ðÞ ðÞ , k x ¼ nπx ðÞ =a, and k y ¼ mπy ðÞ =b. Similarly, we can calculate the value of g y n; m ðÞ , where g y x; y ðÞ ¼ 1=e x; y ðÞ ðÞ de x; y ðÞ =dy ðÞ .
For the cross section as shown in Figures 1(a) and 3(b), the center of the rectangle is located at (0.5 a, 0.5 b), y 12 = b/2 + c/2 and y 11 = b/2Àc/2. Thus, for this case, y 12 À y 11 = c and y 12 þ y 11 =b. The equation of the circle is given by x À a=2 ðÞ 2 þ y À b=2 ðÞ 2 ¼ r 2 . The center of the circle is located at (0.5 a, 0.5 b), as shown in Figure 1(b). We obtain two possibilities without this The dielectric profile for the circle is given where the center is located at (0.5 a, 0.5 b) (Figure 1 where else gx ; y ðÞ = 0. The radius of the circle is given by r ¼ The derivatives of the dielectric profile for the circle are given where the center is located at (0.5 a, 0.5 b), as shown in Figure 1 else g x = 0, and g y =0.
The elements of the matrices for the circular profile are given where the center is located at (0.5 a, 0.5 b) (Figure 1 Applications and Solving Techniques of Propagated Wave in Waveguides Filled with Inhomogeneous Dielectric… http://dx.doi.org/10.5772/intechopen.76793 g y n; m ðÞ ¼ À 2g 0 ab where r ¼ The proposed techniques in this subsection and the subsection 3.1 relate to the method for the propagation along the straight rectangular metallic waveguide [32]. The techniques and the particular applications to solve the rectangular and circular profiles in the cross section of the straight rectangular waveguide are important in order to improve the mode model [32]. The method is based on the Laplace and Fourier transform and their inverse transforms. Laplace transform is necessary to obtain the comfortable and simple input-output connections of the fields. The output transverse field profiles are computed by the inverse Laplace and Fourier transforms.
The matrix G is given by the form.
Similarly, the G x and G y matrices are obtained by the derivatives of the dielectric profile. These matrices relate to the method that is based on the Laplace and Fourier transforms and their inverse [32].
Several examples will demonstrate in the next section in order to explore the effects of the rectangular and circular materials in the cross section (Figure 1(a) and 1(b)) along the straight waveguide on the output field. All the graphical results will be demonstrated as a response to a half-sine (TE 10 ) input-wave profile and the rectangular and circular materials in the cross section of the straight rectangular waveguide.
3.3. The technique based on ω ε function for the circular profile in the cross section of the straight circular waveguide The cross section of hollow waveguide (Figure 1(c)) is made of a tube of various types of three dielectric layers and a metallic layer. The internal and external diameters are denoted as 2b, 2b 1 ,2 b 2 , 2a, and 2(a+δ m ), respectively, where δ m is the thickness of the metallic layer. In addition, we denote the thickness of the dielectric layers as d 1 , d 2 , and d 3 , respectively, where The refractive index in the particular case with the three dielectric layers and the metallic layer in the cross section of the straight hollow waveguide (Figure 1(c)) is calculated as follows: nr ðÞ¼ where the parameters ε 1 , ε 2 , ε 3 , and ε 4 are very small [e.g., ε 1 = ε 2 = ε 3 = a À b ½ =50, ε 4 = δ m /50]. The refractive indices of the air, dielectric and metallic layers are denoted as n 0 , n 1 , n 2 , n 3 , and n m , respectively. In this study, we suppose that n 3 > n 2 > n 1 .
The transverse derivative of the dielectric profile is calculated as follows: The proposed technique in this subsection relates to the theoretical model that based on Maxwell's equations, Fourier-Bessel series, Laplace transform, and the inverse Laplace transform [33]. In this theoretical model, the longitudinal components of the fields are developed into the Fourier-Bessel series. The transverse components of the fields are expressed as functions of the longitudinal components in the Laplace plane and are obtained by using the inverse Laplace transform by the residue method. The separation of the variables is obtained by using the orthogonal relations.
The output transverse components of the fields of the straight hollow waveguide are finally expressed in a form of transfer matrix function, and the derivation has been already explained in detail in Ref. [33]. The contribution of the proposed technique to calculate the refractive indices (nr ðÞ ) and the transverse derivative of the dielectric profile (g r ) for three dielectric layers is important to improve the method that is based on the mode model [33]. This improved method is important to reduce the transmission losses of the dielectric coated metal waveguides.

Numerical results
Several examples for all geometry of the rectangular and circular waveguides and the dielectric profile are demonstrated in this section for three cases, as shown in Figure 1(a)-(c).

Numerical results for the rectangular dielectric material in the cross section of the straight rectangular waveguide
The analytical method for the dielectric slab [35] is shown in Figure 4(a). The slab profile in the cross section is based on transcendental equation, as follows: where ν The criterion for the convergence of the solution is given by , (18) for N ≥ 1, where 2N þ 1 ðÞ 2 is the number of the modes. By increasing the order N, then E y N ðÞ approaches E y .    (18). Between N = 7 and N = 9 the value of the criterion is equal to À2. According to Figure 4(b) and 4(c), the comparison has shown good agreement.
All the next graphical results are demonstrated as a response to a half-sine (TE 10 ) input-wave profile. Figure 5(a)-(f) relates to discontinuous problem according to Figure 1(a). Figure 5(a)-(d) shows the results of the output field as a response to a half-sine (TE 10 ) input-wave profile. In this case, c=d=3.3mm and the center of the rectangle is located at the point (0.5 a, 0.5 b), as shown in Figure 1 (a) for e r =3,5,7,and10,respectively.Theotherparametersarea=b=2 cm ,z=15cm, k 0 =1671=m, λ =3.75cm,andβ =581=m. The output fields are strongly affected by the input wave profile (TE 10 mode), the rectangular profile, and the location of the center of the rectangle (0.5 a, 0.5 b).  By increasing only the value of e r of the rectangular dielectric profile, the TE 10 wave profile decreased, the Gaussian shape of the output field increased, and the relative amplitude decreased. In addition, by increasing only the value of e r , the width of the Gaussian shape decreased.
The output profiles and the amplitudes for N = 1, 3, 5, and 7 are shown in Figure 5(f), for e r = 10. By increasing only the parameter of the order N, the output field approaches to the final output field.
According to the results, we see that the output fields are strongly affected by the input wave profile (TE 10 mode), the rectangular profile, and the location of the center of the rectangle (0.5 a, 0.5 b).

4.2.
Numerical results for the circular dielectric material in the cross section of the straight rectangular waveguide Figure 6(a)-(f) relates to discontinuous problem according to Figure 1(b). Figure 6(a)-(d) demonstrates the results of the output field as a response to TE 10 input-wave profile, for e r = 3, 5, 7, and 10, respectively. The radius of the circle in this case is equal to 2.5 mm and the center of the circle is located at the point (0.5 a, 0.5 b), as shown in Figure 1(b).
The other parameters are a = b = 2 cm ,z=15cm, k 0 =1671=m, λ =3.75cm,andβ =581=m.  Emerging Waveguide Technology 308 The output profiles are shown in Figure 6 By increasing only the parameter e r of the circular dielectric profile from 3 to 10, the Gaussian shape of the output transverse profile of the field increased, the TE 10 wave profile decreased, and the relative amplitude of the output field decreased.
The output fields of Figure 6(a)-(f) are strongly affected by the input wave profile (TE 10 mode), the circular profile, and the location of the center of the circle (0.5 a, 0.5 b).
It is interesting to see a similar behavior of the output results in the cases of rectangular profiles ( Figure 5(a)-(f)) that relate to Figure 1(a) and in the cases of circular profiles (Figure 6(a)-(f)) that relate to Figure 1(b), for every value of e r , respectively. According to these output results, we see the similar behavior for every value of e r , but the amplitudes of the output fields are different.

Numerical results for the three dielectric layers and a metallic layer in the cross section of the straight hollow waveguide
The comparison between our theoretical result of the output power density with the laboratory result in the case of the straight hollow waveguide with one dielectric layer and a metallic layer is demonstrated in Figure 7(a) and 7(b). In this example, the diameter (2a) of the waveguide is 2 mm, the thickness of the dielectric layer [d AgI ðÞ ]is0:75 μm, and the minimum spot-size (w 0 )is 0.3 mm. The length of the straight waveguide is 1 m. The refractive indices of the air, the dielectric layer (AgI) and the metallic layer (Ag) are n 0 ðÞ ¼ 1, n AgI ðÞ ¼ 2:2, and n Ag ðÞ ¼ 13:5 À j75:3, respectively. The value of the refractive index of the material at a wavelength of λ=10.6 μm is taken from the table compiled by Miyagi et al. [36].
The results of the output power density jS av j ðÞ (e.g., Figure 7(a)) show the behavior of the solutions for the TEM 00 mode in excitation. The comparison between our theoretical result (Figure 7(a)) and the published experimental data [37], as shown also in Figure 7(b) shows good agreement of a Gaussian shape as expected, except for the secondary small propagation mode. The experimental result is taken into account the roughness of the internal wall of the waveguide, but our theoretical model is not taken the roughness.
Miyagi and Kawakami showed that transmission losses of the dielectric-coated metal waveguides are drastically reduced when a multiple dielectric layer is formed instead of a single dielectric layer [31]. The simplest and most efficient multilayer structure is three dielectric layers that deposited on a metal layer.
Thus, in this subsection, we present the output results that relate to the proposed technique and the particular application for the cross section of the straight hollow waveguide with three dielectric layers (and not only with one dielectric layer) and a metallic layer. In this case, we can improve the method [33], in the case of the three dielectric layers and a metallic layer in the cross section of the straight hollow waveguide.

Emerging Waveguide Technology
By changing only the spot size from w 0 = 0.1 to 0.3 mm, with the same other parameters, the output power density is changed, as shown in Figure 8(a) and 8(b). The results of the output power density jS av j ðÞ show the behavior of the solutions for the TEM 00 mode in excitation. By changing only the spot size from w 0 = 0.3 to 0.1 mm, the width of the output Gaussian becomes more narrow. The output field results are strongly affected by the spot size and the structure of the three layers and the metallic layer in the cross section of the straight hollow waveguide.
This mode model can be a useful tool to predict the relevant parameters in the case of the hollow waveguide with three dielectric layers and a metallic layer (Ag) for practical applications (output fields, output power density and output power transmission), before carrying out experiments in the laboratory.

Conclusions
The main objective of this chapter was to develop the techniques that enable us to solve problems with inhomogeneous dielectric materials in the cross section of the straight rectangular and circular waveguides. The proposed techniques are very effective in relation to the conventional methods because they allow the development of expressions in the cross section only according to the specific discontinuous problem. In this way, the mode model methods become an improved methods to solve discontinuous problems in the cross section (and not only for continuous problems).
Three examples of inhomogeneous dielectric profiles in the cross section of the straight waveguides were shown in Figure 1(a)-(c). Figure 1(a) and 1(b) show the rectangular and circular profiles in the cross section of the straight rectangular waveguide, respectively. Figure 1(c) shows three dielectric layers and a metallic layer in the cross section of the straight hollow waveguide. The second objective is to understand the influence of the inhomogeneous dielectric materials on the output fields.
The proposed techniques are important to improve the methods that are based on Laplace and Fourier transforms and their inverse transforms also for the discontinuous rectangular and The result of the comparison between the theoretical models with the analytical solution [35] is shown in Figure 4(b) and the convergence is shown in Figure 4(c). The comparison has shown good agreement.
The results for the rectangular profile in the cross section of the straight rectangular waveguide (Figure 1(a)) were shown in Figure 5(a)-(f). Figure 5(a)-(d) shows the results of the output field as a response to a half-sine (TE 10 ) input-wave profile. The output fields are strongly affected by the input wave profile (TE 10 mode), the rectangular profile, and the location of the center of the rectangle (0.5 a, 0.5 b).
By increasing only the value of e r of the rectangular dielectric material, the Gaussian shape of the output field increased, the TE 10 wave profile decreased, the relative amplitude decreased, and the width of the Gaussian shape decreased. The output field approaches to the final output field, by increasing only the parameter of the order N, as shown in Figure 5(f).
The results for the circular profile in the cross section of the straight rectangular waveguide (Figure 1(b)) were shown in Figure 6(a)-(f). Figure 6(a)-(d) shows the results of the output field as a response to a half-sine (TE 10 ) input-wave profile.
By increasing only the value of the parameter e r of the circular dielectric material (Figure 1(b)) in the rectangular cross section from 3 to 10, the Gaussian shape of the output transverse profile of the field increased, the TE 10 wave profile decreased, and the relative amplitude of the output field decreased.
The output fields of Figure 5 It is interesting to see a similar behavior of the output results in the cases of rectangular profiles ( Figure 5(a)-(f))thatrelatetoFigure 1(a) andinthecasesofcircularprofiles( Figure 6(a)-(f)) that relate to Figure 1(b), for every value of e r , respectively. According to these output results, we see the similar behavior for every value of e r , but the amplitudes of the output fields are different.
The comparison between our theoretical result (Figure 7(a)) and the published experimental data [37], as shown also in Figure 7(b) shows good agreement of a Gaussian shape as expected, except for the secondary small propagation mode. The experimental result is taken into account the roughness of the internal wall of the waveguide, but our theoretical model is not taken the roughness.
The output power density for the straight hollow waveguide with three dielectric layers (Figure 1(c)) is shown in Figure 8(a)-(c). The output field results are strongly affected by the spot size and the structure of the three layers and the metallic layer in the cross section of the straight hollow waveguide.
These models are useful to predict the structure of the output fields for rectangular and circular profiles in straight waveguides in the cases of microwave and millimeter-wave regimes and in