CC BY
35
DOI: 10.14529/ctcr160206
EVALUATION OF THE RADIATION PATTERNS
OF AIRCRAFT SLOT ANTENNAS USING ITERATIVE TECHNIQUES
I.M. Bigil'din, bigildin93@mail.ru,
A.B. Khashimov, xab@kipr.susu.ac.ru
South Ural State University, Chelyabinsk, Russian Federation
The rigorous electromagnetic formulation of the problem to define the pattern of slot antennas mounted on aircraft is used Fredholm's integral equation of II kind. For calculation of the slot antenna far-field pattern on perpendicular plane to aircraft axis is used an asymptotic correspondence of 3D and 2D problems in antenna engineering. We assume that slot antenna axis is parallel to aircraft axis. It is shown such correspondence is increasing the abilities of modeling the radiation pattern of slot antennas and arrays. High-precision Gauss quadrature formulas is used for matrix elements calculation. Turn from direct numerical tecniques of the integral equation solving to iterative procedures requires preliminary analysis of the using methods, especially for convergence of large-scale problems. On the base of multiple numerical experience for solving some widely-known test poblems we defined the most methods for solving such problems - conjugate gradients squared method and generalized minimum residual method. Numerical result shows that patterns of the aircraft slot antennas have a large number of directional nulls.
Keywords: slot antenna, aircraft, pattern, asymptotic correspondence, integral equations, iterative techniques.
The slot antennas often used as radiating system of radio equipment for an aircraft. An advantage of such antennas is feasible combining and mounting on an aircraft surface. The slot antennas design should be considered the influence of an aircraft surface. It means that the slot antenna is a part of an aircraft structure and surface, typically, is metal of high specific conductivity. Thus, primary electromagnetic field radiated the slot antenna to induce on an aircraft surface the current is essentially a surface current, it is a source of the secondary scattering field. Therefore, the exciting slot antenna and an aircraft surface form a single radiating system. Total radiation field is the result of interference of the primary and secondary fields. This leads to significant distortions of the slot antenna own pattern, as scattering field defines complex shape of the aircraft surface. Furthermore, for the aircraft slot antenna it is necessary to use cavity for providing optimal pattern. Thus, to own resonant characteristics of the slot antenna add loaded cavity-resonator Q, therefore band is significantly narrowed [1]. The slot antenna location on an aircraft surface depends on the desired pattern properties. For telemetry and command control systems is desirable to use patterns, close to uniform on a plane perpendicular to the axis of an aircraft. We consider dependence the slot antenna pattern on it location on an aircraft surface, assuming the linear polarization of it.
Applying the rigorous electromagnetic equations to form mathematical model of such problem we can consider the detail dependence of mutual coupling of all part an aircraft surface and the slot antennas or the slot array. Known approximate methods of physical optics [2] can provide useful evaluations of the pattern distortion for complex shape objects, but complete form of the scattered fields may be obtained only on basis of the solution of the corresponding rigorous electromagnetic problem.
We use an asymptotic correspondence of 3D vector and 2D scalar electromagnetic problems [3]. It is possible only for extensive aircraft surface and only for defining the patterns on a plane perpendicular to the axis of an aircraft. Thus, in this case the original slot antenna is substitute the filament of the equivalence magnetic current. The limit radius of the filament is a << X, where X is wavelength of the driving field. We assume that the slot antenna and magnetic current of the filament has equal driving amplitude and phase. Assuming a perfectly conducting cylindrical surface, the boundary conditions for the tangential components of the total electric field leads to integral equation of II kind [4]:
/ (p) + f J / (q)H(2)(К* ) ^^ ^ = 2H (p), (1)
2 L rpq
where jv is an electric current on the contour L, as cut of an aircraft surface S the plane ф = 0 (Fig. 1);
H(12)(krpq ) is the Hankel function first order and second kind for corresponding argument; rpq is
the radius-vector is linked points of view and source; Hl is tangent to contour L component of driving magnetic field, in this case source of it is the filament of the magnetic current. We assume known the given amplitude and phase distribution I0, у0 . Besides, for the aircraft cavity-resonator mounted slot antennas we assume that filament of the equivalence magnetic current is located above the perfectly conducting surface. Then, the current jv can define the far-field pattern of the total radiating system, but only for given cut plane.
Рис. 1. Typical versions to the aircraft slot antennas mounting
For numerical solution of the integral equation Eq. 1 is used a special case of the method of moments as the collocation method [5], leads Eq. 1 to the system of linear algebraic equations, often large-size. The noticeable feature of this system is the dominant main diagonal providing stability of numerical solution, whether direct or iterative. Calculation of the matrix elements is used the Gauss quadrature formulas, the choice of the number of nodes is based on the adaptive procedure to analyze the stability and accuracy of numerical solution, depending on the size and configuration of the contour. For most practically used solutions the choice of 8...12 nodes is sufficient to ensure the given relative accuracy of the numerical solution proportionally 10-3.
For solving the large-scale system of linear algebraic equations A • x = b it is necessary to use iterative methods. Matrix A must be square and should be large. The column vector b must have length n . For large-scale system integer n can greater then some ten thousands. We use several MATLAB iterative methods [6], but common feature of these methods is using an auxiliary handle function my fun to represent the numerical solution of the system as x = solver (A, b), where solver is using MATLAB
code. For example, function my _ fun (x,' notransp') must return column vector is equal to A • x . We consider some iterative methods to evaluate the efficiency of solving large-scale system.
1. qmr (Quasi-Minimal Residual method).
2. bicg (BiConjugate Gradient Conjugate Gradients).
3. bicgstab and bicgstabl (BiConjugate Gradients Stabilized method).
4. cgs (Conjugate Gradients Squared method).
5. gmres (Generalized Minimum Residual method).
6. lsqr (Least Squares method).
7. minres (Minimum Residual method).
8. pcg (Preconditioned Conjugate Gradients method).
9. symmlq (Symmetric LQ method).
10. tfqmr (Transpose-Free Quasi-Minimal Residual method).
The criterion for the efficiency of these methods is computing time to achieve the given accuracy. Comparison of the numerical results accuracy is obtained due to solve 2D test problem as finding the radiating field of the magnetic dipole on a perfectly conducting elliptic cylinder. In this case, the solution can be represent as an infinite series for Mathieu, Mathieu-Bessel, Mathieu-Hankel functions.
Бигильдин И. М., Хашимов А.Б.
Итерационные методы оценки поля излучения бортовых щелевых антенн
Using this analytical form of solving we can obtain the useful numerical results due to special adaptive procedure for choosing the size of summed series to provide given accuracy. For numerical solution of the Eq. 1 contour L is divided into n reasonable segments. For example, we assume that n = 7120. According to the results of numerical experiences is obtained the following conclusions.
1. Iterative methods pcg, minres, symmlq not provide a convergent solutions for reasonable computing time.
2. Most preferred iterative methods to solve Eq. 1 is qmr, cgs, gmres, Isqr. These methods provide a stable convergent solution, but gmres is the fastest method.
Fig. 2, 3 depict some patterns of the aircraft slot antennas and arrays on plane perpendicular to an aircraft axis (logarithmic scale in dB is used). We consider the several versions of the slot antennas location on an aircraft surface. Evidently, the patterns have a large number of directional nulls, but due to changes of an aircraft turns the interruptions of the radio channels will be short and quite acceptable.
Fig. 2. The missile pattern of phase quadrature driven slot array
Fig. 3. The aircraft slot antenna pattern disposed on the edge
References
1. Hansen R.C. Phased Array Antennas, J. Wiley & Sons, 2009. 560 p.
2. Ufimtsev P.Ya. Fundamentals of the Physical Theory of Diffraction, J. Wiley & Sons, 2007. 329 p.
3. Voitovich N.I., Khashimov A.B. On the Correspondence of Asymptotic Solutions to 2D and 3D Problems in Antenna Engineering. Journal of Communications Technology and Electronics, 2010, vol. 55, no. 12, pp. 1374-1379. Pleiades Publishing, Inc., 2010.
4. Khashimov A.B., Salikhov R.R., Al'metov R.S. [Using of High Precision Numerical Techniques for Antenna System Design]. Bulletin of the South Ural State University. Ser. Computational Mathematics and Software Engineering, 2014, vol. 3, no. 2, pp. 77-91. (in Russ.)
5. Computer Techiques for Electromagnetics. Pergamon Press, 1973. 488 p.
6. Barrett R., Berry M., Chan T.F. et al. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods. SIAM, Philadelphia, 1994. 368 p.
Received 11 February 2016
УДК 621.396.6(07), 537.8(07) DOI: 10.14529/^сг160206
ИТЕРАЦИОННЫЕ МЕТОДЫ ОЦЕНКИ ПОЛЯ ИЗЛУЧЕНИЯ БОРТОВЫХ ЩЕЛЕВЫХ АНТЕНН
И.М. Бигильдин, А.Б. Хашимов
Южно-Уральский государственный университет, г. Челябинск
Для строгой электродинамической формулировки задачи определения диаграммы направленности щелевых антенн, устанавливаемых на летательных аппаратах, используется интегральное уравнение Фредгольма II рода. Для расчета диаграммы направленности щелевой антенны в дальней зоне в плоскости, перпендикулярной оси летательного аппарата, используется асимптотическое соответствие трехмерных и двумерных задач в антенной технике. Предполагается, что ось щелевой антенны параллельна оси летательного аппарата. Показано, что такое соответствие расширяет возможности моделирования поля излучения щелевых антенн и антенных решеток. Для вычисления матричных элементов используются квадратурные формулы Гаусса высокой точности. Переход от прямых численных методов решения интегральных уравнений к итерационным методам требует предварительного анализа характеристик используемых методов, особенно вопросов сходимости в задачах большой размерности. На основе многочисленных вычислительных экспериментов решения хорошо известных тестовых задач определены наиболее приемлемые итерационные методы для решения сформулированных задач - квадратичный метод сопряженных градиентов и метод обобщенных минимальных невязок. Численные результаты решения показывают, что диаграммы направленности бортовых щелевых антенн имеют многочисленные глубокие узкие провалы.
Ключевые слова: щелевая антенна, летательный аппарат, диаграмма направленности, асимптотическое соответствие, интегральные уравнения, итерационные методы.
Литература
1. Hansen R.C. Phased Array Antennas / Hansen R.C. - J. Wiley & Sons, 2009. - 560 p.
2. Ufimtsev, P.Ya. Fundamentals of the Physical Theory of Diffraction. - J. Wiley & Sons, 2007. -329 p.
3. Voitovich, N.I. On the Correspondence of Asymptotic Solutions to 2D and 3D Problems in Antenna Engineering / N.I. Voitovich, A.B. Khashimov // Journal of Communications Technology and Electronics. - 2010. - Vol. 55, no. 12. - P. 1374-1379. Pleiades Publishing, Inc., 2010.
Бигильдин И. М., Хашимов А.Б.
Итерационные методы оценки поля излучения бортовых щелевых антенн
4. Хашимов, А.Б. Применение вычислительных схем повышенной точности в проектировании антенных систем / А.Б. Хашимов, Р.Р. Салихов, Р.С. Альметов // Вестник Южно-Уральского университета. Серия «Вычислительная математика и информатика». - 2014. - Т. 3, № 2. -
5. Computer Techiques for Electromagnetics /Ed. by R. Mittra. - Pergamon Press, 1973. - 488p.
6. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods / R. Barrett, M. Berry, T.F. Chan et al. - SIAM, Philadelphia, 1994. - 368p.
Бигильдин Ильдар Мансурович, магистрант кафедры конструирования и производства радиоаппаратуры, Южно-Уральский государственный университет, г. Челябинск; bigildin93@ mail.ru.
Хашимов Амур Бариевич, канд. физ.-мат. наук, доцент кафедры конструирования и производства радиоаппаратуры, Южно-Уральский государственный университет, г. Челябинск; xab@kipr.susu.ac.ru.
С.77-91.
Поступила в редакцию 11 февраля 2016 г.
ОБРАЗЕЦ ЦИТИРОВАНИЯ
FOR CITATION
Bigil'din, I.M. Evaluation of the Radiation Patterns of Aircraft Slot Antennas Using Iterative Techniques / I.M. Bigil'din, A.B. Khashimov // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2016. - Т. 16, № 2. - С. 63-67. DOI: 10.14529/ctcr160206
Bigil'din I.M., Khashimov A.B. Evaluation of the Radiation Patterns of Aircraft Slot Antennas Using Iterative Techniques. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2016, vol. 16, no. 2, pp. 63-67. DOI: 10.14529/ctcr160206
