A new method for exact computation of the differential phase shift in the circular waveguide, loaded with an azimuthally magnetized ferrite cylinder and a dielectric toroid
A novel numerical technique for calculating the differential phase shift, brought forth by the circular waveguide, containing a ferrite cylinder of azimuthal magnetization and a dielectric toroid that supports normal TE<sub>0n</sub> modes, is worked out. It is based on the method of successive approximations and yields the roots a<sub>1</sub> and a<sub>2</sub> |a<sub>1</sub> |a<sub>2</sub> of a biquadratic equation, connecting the off diagonal ferrite Polder permeability tensor element a and the eigenvalue spectrum of the configuration. The spectrum depends on the relative permittivities of inner and outer medium e<sub>r</sub> and e<sub>d</sub> resp. the factor r, determining the degree of filling of geometry with anisotropic material, the element a and the normalized regarding frequency and the constant e<sub>r</sub> guide radius r<sub>0</sub> It is calculated by the roots of the structure's characteristic equation, derived formerly through complex Kummer confluent hypergeometric and real Bessel and Neumann functions. For chosen values of the quantities e<sub>r</sub> e<sub>d</sub> r, |a| and r<sub>0</sub> the imaginary part of the complex first parameter of confluent functions is varied, until the modulus of counted root a<sub>1</sub> (or a<sub>2</sub> coincides with the selected numerical equivalent of the element a within the range of the prescribed accuracy. If |a| |a<sub>1</sub> |a| |a<sub>2</sub> for the normalized as pointed out above_phase constant b of the electromagnetic field it holds b b<sub>(1)/sub> (b= b<sub>(1)/sub> where b<sub>1</sub> |a<sub>2</sub> (b<sub>(2)/sub> |a<sub>l</sub> (b<sub>(1)/sub> b(2) For each set of the parameters mentioned the procedure is repeated twice for positive (counterclockwise) and negative (clockwise) direction of the magnetization of its anisotropic load towards that of the wave propagation. The difference between the two computed constants gives the phase shift sought. Numerical outcomes for three geometries are obtained in case of TE<sub>01</sub> mode, assuming e<sub>r</sub> =e<sub>d</sub> The approach is rather complicated, too cumbersome and very time consuming, but it provides extraordinarily high accuracy of the results (more than ten decimal places) comparable with the one, due to the recently developed for the same purpose method, invented also by the authors. The combined application of both schemes makes possible to disclose completely the entire picture of the phase shifting behaviour of transmission line considered.
2016 International Conference on Electromagnetics in Advanced Applications (ICEAA)