Modeling of Plane Arrays Using a Variational Approach

Affiliation(s)

^{1}
Department of Numerical Methods of Mathematical Physics, Pidstryhach Institute for Applied Problems of Mechanics
and Mathematics, NASU, Lviv, Ukraine.

^{2}
Department of Informatics, Lviv State University of Internal Affairs, Lviv, Ukraine.

Abstract

The variational statement of synthesis problem is generalized in order to account
the additional requirements to the synthesized radiation pattern (RP)
and field distribution in the specified points of near zone. For this aim, the
minimizing functional is supplemented by term providing the possibility to
minimize the values of field in these points; creating the deep zeros in the RP
for the certain angular coordinates is realized too. The approach foresees reduction
of an explicit formula for field values in a near zone. The results of
computational modeling testify the possibility to create zeros in the given RP
and to minimize the values of field in a near zone of plane arrays in a great
extent.

Keywords

Synthesis Problem, Variational Statement, Plane Array, Near Zone, Euler’s Equation, Numerical Modeling

Synthesis Problem, Variational Statement, Plane Array, Near Zone, Euler’s Equation, Numerical Modeling

1. Introduction

A lot of modern antenna systems are functioning by condition of a different series of various requirements to their radiation characteristics. One of such requirements is forming the deep minimums or zeros in the given points of near or far zones of antenna. This requirement corresponds to issues of the electromagnetic compatibility (EMC). Of course, all antenna systems are working at the condition that is characterized by the existence of various types of radio electronic systems and devices in its neighboring environment. Many engineering problems appear in this connection. The solution to these problems should provide the effective functioning neighbor antenna radiating systems that are placed closely. Therefore, the need for solving two interrelated problems of EMC appears. They are distribution of frequencies while working in the assigned frequency range, and a constructive improvement of the radio electronic system’s individual elements.

The EMC problems were stated and discussed since 80^{th} of last century. In this regard, it should be noted pioneer monograph [1] , and many journal papers (see, for example [2] [3] [4] [5] and references therein). The last papers appurtenant to this problem show that the development of the methods and means for providing the certain requirements to the field characteristics of antenna systems in the near and far zones is actual now. This leads to necessity of development of the modeling tools and respective software for solving the practical engineering problems related to EMC of antenna systems.

Some general recommendations how to improve EMC issues for radiating systems are set out in [6]. Different technical solutions and measurements for improving the EMC requirements in the neighboring environment are discussed in papers [7] [8] [9]. Application of numerical techniques in order to provide the series of restrictions on the radiation characteristics of antennas on the whole [10] [11] , to RP only [12] [13] or to characteristics radiation in a near zone [14] gives possible to elaborate the general vision to EMC problem solving. However, a comprehensive solution to the EMC problem for radiating systems and antennas is far from complete. The generalized variational approach for solving the array’s synthesis problems, proposed in this paper, is one of the means that allows to solve this problem at a more fundamental level. The mathematical statement of synthesis problem including the EMC specifications for plane arrays and results of modeling are discussed in the presented paper.

The paper is organized as follows.

The generalization of usual synthesis problem for plane array allowing to prescribe the restrictions on the radiation characteristics in a near zone is proposed and explanations to terms of formulated functional are given in Section 2. The formula for RP of array (array factor) is discussed with introducing the generalized angular coordinates. In Section 3, the procedure of obtaining the nonlinear integral equation for optimal distribution of currents in array elements, corresponding to minimizing functional, is described. Reduction of explicit formula for components of electromagnetic (EM) field in a near zone of array is made in Section 4. Section 5 and Section 6 are devoted to methods of creating zeros (deep gaps) in RP of array and values of EM field in a near zone. The computational results of modeling for rectangular and hexagonal arrays are presented and discussed in Sections 7; Section 8 contains concluding remarks and proposal for the future investigations.

2. Variational Statement of Synthesis Problem

The functional aimed to take into account a series of requirements to both the RP and field distribution in a near zone is formulated as sum of several terms [15] [16] [17] , each of them ensures the best approximation to the array’s radiation characteristics [18]

$\begin{array}{c}\sigma ={\displaystyle \underset{0}{\overset{\text{\pi}}{\int}}{\displaystyle \underset{0}{\overset{\text{2\pi}}{\int}}p\left(\theta ,\phi \right){\left[F\left(\theta ,\phi \right)-\left|f\left(\theta ,\phi \right)\right|\right]}^{2}\text{d}\phi \text{d}\theta}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\displaystyle \underset{0}{\overset{\text{\pi}}{\int}}{\displaystyle \underset{0}{\overset{\text{2\pi}}{\int}}q\left({\theta}^{\prime},{\phi}^{\prime}\right){\left[U\left({\theta}^{\prime},{\phi}^{\prime}\right)-c\left|u\left({\theta}^{\prime},{\phi}^{\prime}\right)\right|\right]}^{2}\text{d}{\phi}^{\prime}\text{d}{\theta}^{\prime}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\displaystyle \underset{n=-N,m=-M}{\overset{N,M}{\sum}}{t}_{nm}{\left|{I}_{nm}\right|}^{2}},\end{array}$ (1)

function $F\left(\theta ,\phi \right)$ is prescribed amplitude RP (by the condition of problem, it is real and positive), $\left|f\left(\theta ,\phi \right)\right|$ characterizes the module of RP created (synthesized) by array, $p\left(\theta ,\phi \right)\ge 0$ is a weight function that aimed to formulate specific demands to the module of $f\left(\theta ,\phi \right)$ , in particular, to create zeros.

The second term in (1) provides an opportunity to minimize the values of field’s amplitude in the fixed points of array’s near zone. Such an assignment is fulfilled by means of increasing the $q\left({\theta}^{\prime},{\phi}^{\prime}\right)$ value in the above points in comparison with the rest of them. Value $U\left({\theta}^{\prime},{\phi}^{\prime}\right)$ is amplitude of given field in a near zone, by this $u\left({\theta}^{\prime},{\phi}^{\prime}\right)$ is field created by array. The value of additional optimizing parameter c is determined as $\text{d}\sigma /\text{d}c=0$ .

The last term in (1) is applied to restrict the amplitude of currents in the array’s elements. This provides to solve the problem to make away with the effect of super-directivity [19] [20] [21] for array. The additional optimizing values ${t}_{nm}$ are applied to restrict the amplitude of currents in the array’s radiators.

We use for RP (array factor) of rectangular array the next formula [18]

$f\left({s}_{1},{s}_{2}\right)=AI\equiv {\displaystyle \underset{m=-{M}_{1}}{\overset{{M}_{1}}{\sum}}{\displaystyle \underset{n=-{N}_{1}}{\overset{{N}_{1}}{\sum}}{f}_{nm}\left({s}_{1},{s}_{2}\right){I}_{nm}{\text{e}}^{i\left({c}_{1}n{s}_{1}+{c}_{2}m{s}_{2}\right)}}}$ , (2)

the RP of hexagonal array is determined similarly [22]. ${I}_{nm}$ are the currents in array’s radiators, ${f}_{nm}\left({s}_{1},{s}_{2}\right)$ are the RPs of separate radiators.

We will use below the generalized angular coordinates ${s}_{1},{s}_{2}$ instead of usual $\theta ,\phi $ . This provides ability to present the dependence of RP of array on the one generalized parameter that includes the geometrical size of array and length of wave or frequency simultaneously. Such coordinates are ${s}_{1}=\mathrm{sin}\theta \mathrm{cos}\phi /\mathrm{sin}{\alpha}_{1}$ , ${s}_{2}=\mathrm{sin}\theta \mathrm{sin}\phi /\mathrm{sin}{\alpha}_{2}$ where ${c}_{1}=k{d}_{1}\mathrm{sin}{\alpha}_{1}$ , ${c}_{2}=k{d}_{2}\mathrm{sin}{\alpha}_{2}$ , $k=2\text{\pi}/\lambda $ is wavenumber, $\lambda $ is length of wave, values ${d}_{1}$ and ${d}_{2}$ correspond to distance between elements of array in Ox and Oy directions. It is assumed that the given RP F differs from zero within the angles $\left[\pm {\alpha}_{1},\pm {\alpha}_{2}\right]$ .

3. Numerical Minimization of Functional

In order to find the minimum of (1), we apply the approach of the variation calculus that yields in receiving the formula for the gradient of functional. Using such formula, we pass to non-linear integral equation [18] with respect to unknown currents in the array’s elements.

The notations can be simplified in a great extent if to use the operator form as $f\left({s}_{1},{s}_{2}\right)=AI$ for synthesized RP. For this goal, the notion of adjoint operators ${A}^{*}$ and ${B}^{*}$ is applied too.

The nm-th component for gradient of (1) on the currents ${I}_{nm}$ can be presented as

$\begin{array}{l}{z}_{nm}=t{I}_{nm}+{A}^{*}\left[p\left(A{I}_{nm}-F{\text{e}}^{i\mathrm{arg}A{I}_{nm}}\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{B}^{*}\left[q\left({c}^{2}B{I}_{nm}-cU{\text{e}}^{i\mathrm{arg}B{I}_{nm}}\right)\right],\left(m=\stackrel{\xaf}{-{M}_{1},{M}_{1}},n=\stackrel{\xaf}{-{N}_{1},{N}_{1}}\right).\end{array}$ (3)

Here, operator B presents explicit relation for the field’s components in a near zone. It is determined similarly to A, the difference is that the additional multiplier, depending on the distance r to array is presented. Equating gradient to zero, we receive the non-linear vector-matrix equation for the optimal currents in array’s elements

$tI+{A}^{*}\left[p\left(AI-F{\text{e}}^{i\mathrm{arg}AI}\right)\right]+{B}^{*}\left[q\left({c}^{2}BI-cU{\text{e}}^{i\mathrm{arg}BI}\right)\right]=0.$ (4)

In practice, we do not deal with solving (4) because in the iterative procedure of its solving there is necessary to calculate the actions of A, B, ${A}^{*}$ and ${B}^{*}$ operators. If one to act on (4) by operator A, we have the nonlinear integral equation of Hammerstein’s type with respect to RP f [18] , because the $A{A}^{*}(\xb7)$ , $A{B}^{*}(\xb7)$ are included into respective exponent.

The explicit relation for operator ${A}^{*}$ is [18]

${\left[{A}^{*}f\right]}_{nm}=\frac{{c}_{1}{c}_{2}}{4{\text{\pi}}^{2}}{\displaystyle \underset{-1}{\overset{1}{\int}}{\displaystyle \underset{-1}{\overset{1}{\int}}f\left({s}_{1},{s}_{2}\right){f}_{{}_{nm}}^{*}\left({s}_{1},{s}_{2}\right)\times {\text{e}}^{-i\left({c}_{1}n{s}_{1}+{c}_{2}m{s}_{2}\right)}\text{d}{s}_{1}\text{d}{s}_{2}}}.$ (5)

In the last formula procedure of integration in limits $\left[0,2\text{\pi}\right]$ and $\left[0,\text{\pi}\right]$ is replaced by twice integration over $\left[-1,1\right]$ . This is result of use of new coordinates ${s}_{1}$ and ${s}_{2}$ instead of initial $\theta $ and $\phi $ . Determination of operator ${B}^{*}$ is realized similarly.

If one to use the subsequent application of A, ${A}^{*}$ , as well as B, ${B}^{*}$ operators and operator form (2) of RP, the numerical calculations are much simple. This allows to receive the nonlinear integral equation relatively RP $f\left({s}_{1},{s}_{2}\right)$ and explicit formulas for the kernel of this equation [18].

4. Determination of Field Components in Near Zone

In order to carry out with explicit relation for the field components in a near zone using the currents ${I}_{nm}$ in the radiators of array, formula for the electric vector potential [23]

${\stackrel{\to}{A}}_{Fr}^{e}\left(r,\theta ,\phi \right)=\frac{{\text{e}}^{-ikr}}{4\text{\pi}R}{\displaystyle \underset{V}{\int}I\left({x}^{\prime},{y}^{\prime},{z}^{\prime}\right){\text{e}}^{ik\left[{r}^{\prime}\mathrm{cos}\alpha +\frac{{{r}^{\prime}}^{2}}{2r}\left(1-{\mathrm{cos}}^{2}\alpha \right)\right]}\text{d}V}$ (6)

is applied. Here

$R\approx r\left[1-\frac{{r}^{\prime}}{r}\mathrm{cos}\alpha +\frac{{{r}^{\prime}}^{2}}{2{r}^{2}}\left(1-{\mathrm{cos}}^{2}\alpha \right)\right]$ (7)

that is we will restrict ourselves to the value of order $O\left(1/{r}^{2}\right)$ . In the case of plane array

${r}^{\prime}\mathrm{cos}\alpha ={x}^{\prime}\mathrm{sin}\theta \mathrm{cos}\phi +{y}^{\prime}\mathrm{sin}\theta \mathrm{sin}\phi $ (8)

so far as the array is placed in the ${x}^{\prime}O{y}^{\prime}$ plane.

For plane array, ${\stackrel{\to}{A}}_{Fr}^{e}\left(r,\theta ,\phi \right)$ can be written as

$\begin{array}{l}{\stackrel{\to}{A}}_{Fr}^{e}\left(r,\theta ,\phi \right)=\frac{{\text{e}}^{-ikr}}{4\text{\pi}R}{\displaystyle \underset{m=-{M}_{1}}{\overset{{M}_{1}}{\sum}}{\displaystyle \underset{n=-{N}_{1}}{\overset{{N}_{1}}{\sum}}{f}_{nm}\left({s}_{1},{s}_{2}\right){\stackrel{\to}{I}}_{nm}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \mathrm{exp}\left[i\left(nk{d}_{1}\mathrm{sin}\theta \mathrm{cos}\phi +mk{d}_{2}\mathrm{sin}\theta \mathrm{sin}\phi \right)\right]{v}_{nm}\left(r,{r}^{\prime},\theta ,\phi \right),\end{array}$ (9)

where

${v}_{nm}\left(r,{r}^{\prime},\theta ,\phi \right)=\mathrm{exp}\left[i\left({x}_{nm}^{2}+{y}_{nm}^{2}-{\left({x}_{nm}\mathrm{sin}\theta \mathrm{cos}\phi +{y}_{nm}\mathrm{sin}\theta \mathrm{sin}\phi \right)}^{2}\right)/2r\right].$ (10)

Here ${\stackrel{\to}{I}}_{nm}$ represent the vector current distributions of array’s elements; by this radius ${r}^{\prime}$ is replaced by $\sqrt{{x}_{nm}^{2}+{y}_{nm}^{2}}$ for all array’s radiators. Without loss of generality, one can assume the separation for ${\stackrel{\to}{I}}_{nm}$ currents onto two axes $O{x}^{\prime}$ and $O{y}^{\prime}$ , and to consider values of ${\stackrel{\to}{A}}_{Fr}^{e}\left(r,\theta ,\phi \right)$ for respective components independently.

Following the above assumptions, we will accept the same relations for the E and H field components in a near zone. As result, we have

${H}_{r}=0$ ${H}_{\theta}=ik{A}_{\phi}^{e}$ ${H}_{\phi}=-ik{A}_{\theta}^{e}$ (11)

${E}_{r}=0$ ${E}_{\theta}=-i\omega \mu {A}_{\theta}^{e}$ ${E}_{\phi}=-i\omega \mu {A}_{\phi}^{e}$ (12)

In such a way, we have obtained relations (9), (11) or (9), (12) for determination E and H field components. These subsequent formulas can be considered as the linear operator B. This is because of linear dependence of E and H on the distribution of currents ${\stackrel{\to}{I}}_{nm}$ . In such consideration, operator ${B}^{*}$ is determined by introducing the respective Hilbertian space, as it was done to obtain ${A}^{*}$ [24].

5. Creation of Zeros in the RP

We rewrite (2) as

$f\left({s}_{1},{s}_{2}\right)=AI\equiv {\displaystyle \underset{m=-{M}_{1}}{\overset{{M}_{1}}{\sum}}{\displaystyle \underset{n=-{N}_{1}}{\overset{{N}_{1}}{\sum}}{f}_{nm}\left({s}_{1},{s}_{2}\right){I}_{nm}{\text{e}}^{i\left({c}_{1}{d}_{n}{s}_{1}+{c}_{2}{d}_{m}{s}_{2}\right)}}}$ , (13)

where

$\begin{array}{l}{d}_{n}=n-\left({N}_{1}+1-{\delta}_{1}/2\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{N}_{1}=\mathrm{entier}\left(N/2\right),\\ {d}_{m}=m-\left({M}_{1}+1-{\delta}_{2}/2\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{M}_{1}=\mathrm{entier}\left(M/2\right).\end{array}$ (14)

Function $f\left({s}_{1},{s}_{2}\right)$ is periodical, namely

$f\left({s}_{1}+\frac{2\text{\pi}\nu}{{c}_{1}},{s}_{2}+\frac{2\text{\pi}\mu}{{c}_{2}}\right)={\text{e}}^{\text{i\pi}\left({\delta}_{1}\nu +{\delta}_{2}\mu \right)}f\left({s}_{1},{s}_{2}\right)$ , (15)

that is important for consideration below.

Formulas (13)-(15) are used for presentation of RP using the Kotelnikov’s formula [25] and forming zeros in the given points $\left({s}_{1}^{*},{s}_{2}^{*}\right)$ . For this goal, we present (13) as

$f\left({s}_{1},{s}_{2}\right)={\displaystyle \underset{n=1}{\overset{N}{\sum}}{\phi}_{n}\left({s}_{2}\right){\text{e}}^{i{c}_{1}{d}_{n}{s}_{1}}}$ , (16)

where

${\phi}_{n}\left({s}_{2}\right)={\displaystyle \underset{m=1}{\overset{M}{\sum}}{I}_{nm}{\text{e}}^{i{c}_{2}{d}_{m}{s}_{2}}}$ . (17)

We use the greed

$\left\{{\Delta}_{1}=\frac{2\text{\pi}}{{c}_{1}{L}_{1}},{\Delta}_{2}=\frac{2\text{\pi}}{{c}_{2}{L}_{2}}\right\},{L}_{1}\ge N,{L}_{2}\ge M$ in the region $\Omega =\{\left[-\text{\pi}/{c}_{1},\text{\pi}/{c}_{1}\right]\times \left[-\text{\pi}/{c}_{2},\text{\pi}/{c}_{2}\right]$ of ${s}_{1}$ and ${s}_{2}$ variables.

Kotelnikov’s formula for (17) is

${\stackrel{\xaf}{f}}_{n}\left({s}_{2}\right)=\frac{1}{{L}_{2}}{\displaystyle \underset{{d}_{\mu}=1}{\overset{{L}_{2}}{\sum}}{\phi}_{n}\left({s}_{2}^{\left(\mu \right)}\right)\frac{\mathrm{sin}\left[{L}_{2}{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]}{\mathrm{sin}\left[{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]}}$ , (18)

where ${s}_{2}^{\left(\mu \right)}={d}_{\mu}{\Delta}_{2}$ .

Substituting (18) into (13), we receive

$f\left({s}_{1},{s}_{2}\right)=\frac{1}{{L}_{2}}{\displaystyle \underset{{d}_{\mu}=1}{\overset{{L}_{2}}{\sum}}\frac{\mathrm{sin}\left[{L}_{2}{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]}{\mathrm{sin}\left[{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]}}\times {\displaystyle \underset{n=1}{\overset{N}{\sum}}{\phi}_{n}\left({s}_{2}^{\left(\mu \right)}\right){\text{e}}^{i{c}_{1}{d}_{n}{s}_{1}}}.$ (19)

Finally, the Kotelnikov’s formula for plane array is

$\begin{array}{c}\stackrel{\xaf}{f}\left({s}_{1},{s}_{2}\right)=\frac{1}{{L}_{1}{L}_{2}}{\displaystyle \underset{{d}_{\nu}=1}{\overset{{L}_{1}}{\sum}}{\displaystyle \underset{{d}_{\mu}=1}{\overset{{L}_{2}}{\sum}}{f}_{0}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)}}\times \frac{\mathrm{sin}\left[{L}_{1}{c}_{1}\left({s}_{1}-{s}_{1}^{\left(\nu \right)}\right)/2\right]}{\mathrm{sin}\left[{c}_{1}\left({s}_{1}-{s}_{1}^{\left(\nu \right)}\right)/2\right]}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot \frac{\mathrm{sin}\left[{L}_{2}{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]}{\mathrm{sin}\left[{c}_{2}\left({s}_{2}-{s}_{2}^{\left(\mu \right)}\right)/2\right]},\end{array}$ (20)

${s}_{1}^{\left(\nu \right)}={d}_{\nu}{\Delta}_{1}$ , ${f}_{0}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)$ are the values of initial PR in the interpolation points. Of course,

$\stackrel{\xaf}{f}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)={f}_{0}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)$ . (21)

Let us minimization point $\left({s}_{1}^{*},{s}_{2}^{*}\right)$ coincides with one of interpolation point, as example $\left({s}_{1}^{\left({\nu}_{0}\right)},{s}_{2}^{\left({\mu}_{0}\right)}\right)$ . We determine the function $\stackrel{\xaf}{f}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)$ in the interpolation points as

$\stackrel{\xaf}{f}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right)=\{\begin{array}{l}0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\nu ={\nu}_{0},\mu ={\mu}_{0},\\ {f}_{0}\left({s}_{1}^{\left(\nu \right)},{s}_{2}^{\left(\mu \right)}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\nu \ne {\nu}_{0},\mu \ne {\mu}_{0}.\end{array}$ (22)

The problem to synthesize zero is formulated as

$\sigma =\Vert \stackrel{\xaf}{f}-{f}^{*}\Vert ={\displaystyle \underset{\Omega}{\iint}{\left|\stackrel{\xaf}{f}\left({s}_{1},{s}_{2}\right)-{f}^{*}\left({s}_{1},{s}_{2}\right)\right|}^{2}\text{d}{s}_{1}\text{d}{s}_{2}}$ , (23)

where ${f}^{*}\left({s}_{1},{s}_{2}\right)$ is RP with formed zero.

The condition of stationarity of (23) leads to formula for currents ${I}_{nm}^{*}$

${I}_{nm}^{*}=\frac{{c}_{1}{c}_{2}}{4{\text{\pi}}^{2}}{\displaystyle \underset{\Omega}{\iint}f\left(s,s\right){\text{e}}^{-i\left({c}_{1}{d}_{n}{s}_{1}+{c}_{2}{d}_{m}{s}_{2}\right)}\text{d}{s}_{1}\text{d}{s}_{2},}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(n=\stackrel{\xaf}{1,N},\text{\hspace{0.17em}}m=\stackrel{\xaf}{1,M}\right)$ . (24)

In fact, values (24) are the coefficients of discrete Fourier transform of function $f\left({s}_{1},{s}_{2}\right)$ . More simple formula for optimal currents ${I}_{nm}^{*}$ is [18]

${I}_{nm}^{*}={I}_{nm}^{0}-\frac{f\left({s}_{1}^{\left({\nu}_{0}\right)},{s}_{2}^{\left({\mu}_{0}\right)}\right)}{{L}_{1}{L}_{2}}{\text{e}}^{-i\left[{c}_{1}{d}_{\nu}{s}_{1}^{\left({\nu}_{0}\right)}+{c}_{2}{d}_{\mu}{s}_{2}^{\left({\mu}_{0}\right)}\right]}$ . (25)

In the case, when point $\left({s}_{1}^{*},{s}_{2}^{*}\right)$ does not coincide with any interpolation point and ${L}_{1}>N$ , ${L}_{2}>M$ , optimal ${I}_{nm}^{*}$ can be determined by the iterative process

${I}_{nm}^{*\left(j+1\right)}={I}_{nm}^{*\left(J\right)}-\frac{{f}^{\left(j\right)}\left({s}_{1}^{\left({\nu}_{0}\right)},{s}_{2}^{\left({\mu}_{0}\right)}\right)}{{L}_{1}{L}_{2}}{\text{e}}^{-i\left[{c}_{1}{d}_{\nu}{s}_{1}^{\left({\nu}_{0}\right)}+{c}_{2}{d}_{\mu}{s}_{2}^{\left({\mu}_{0}\right)}\right]}$ , (26)

where j is iteration number.

Iterative process (26) converges because [26]

$\left|{f}^{\left(j+1\right)}\left({s}_{1}^{\left({\nu}_{0}\right)},{s}_{2}^{\left({\mu}_{0}\right)}\right)\right|<\left|{f}^{\left(j\right)}\left({s}_{1}^{\left({\nu}_{0}\right)},{s}_{2}^{\left({\mu}_{0}\right)}\right)\right|\left(1-\frac{NM}{{L}_{1}{L}_{2}}\right)$ . (27)

6. Creation of Zeros in a Near Zone

Let us ${E}_{\phi}\left({\theta}^{\prime},{\phi}^{\prime}\right)=u\left({\theta}^{\prime},{\phi}^{\prime}\right)$ . The weight function $q\left({\theta}^{\prime},{\phi}^{\prime}\right)$ is used to create zero in the certain point $\left({\theta}^{*},{\phi}^{*}\right)$ of near zone. As a rule, neighbor $\left({\theta}^{*}\pm \Delta \theta ,{\phi}^{*}\pm \Delta \phi \right)$ of this point is used, and integration within $\left[0,\text{\pi}\right]\times \left[0,2\text{\pi}\right]$ is replaced by integration within $\left[{\theta}^{*}\pm \Delta \theta ,{\phi}^{*}\pm \Delta \phi \right]$ .

Minimizing (1), we minimize ${E}_{\phi}\left({\theta}^{\prime},{\phi}^{\prime}\right)$ too, and we have optimal currents ${I}_{nm}^{*}$ [27]. The components of electric vector potential (9) are determined by known ${I}_{nm}^{*}$ . After this, we determine the rest of components of EM field. The radiation power density (RPD) is determined as [22]

${W}_{av}\left(x,y,z\right)=\frac{1}{2}\mathrm{Re}\left[E\times {H}^{*}\right]\text{\hspace{0.17em}}\left(\text{W}/{\text{m}}^{\text{2}}\right)$ . (28)

7. Computational Modeling

7.1. Creating Zeros in RP

Firstly, we show the computational results for the hexagonal array that has 127 elements. By this we assume known RPs ${f}_{nm}\left({s}_{1},{s}_{2}\right)$ of array radiators. The computations are carried out using free MathWorks software.

Results are shown for frequency equal to 11.99 GHz, distance between radiators in the linear subarrays is equal to 0.0078 m, and distance between the linear subarrays is equal to 0.0087 m. The prescribed amplitude RP $F\equiv 1$ in the given range $F\equiv 1$ (Figure 1). In Figure 2, Figure 3, the results that show forming zeros in the vicinity of the first sidelobe are shown for E-plane of array. Using the weight function $p\left(\theta ,\phi \right)$ in functional (1), we form zeros in points $\phi =\text{\pi}/2,\text{3\pi}/2$ and $\theta =\pm {26.4}^{\circ}$ .

Figure 1. Front view of the prescribed RP.

Figure 2. The initial synthesized amplitude RP.

Figure 3. The synthesized RP with zero in the area of first sidelobes, $p=10.0$ .

The level of the first sidelobe was reduced to −52.6 dB (Figure 3). The obtained results testify that such a way to create zeros has some limitation. Use of Kotelnikov’s procedure (25) or (26) is more complicate, but it can provide much small zeros.

7.2. Creating Zeros in a Near Zone

Secondly, the results for antenna with rectangular placement of radiator are considered. Array that consists of 121 radiators is examined, and distance d between them is equal to 0.01 m, that is square of array is equal to 0.01 m^{2}; the work frequency is equal to 5.976 GHz, by this the condition
$d<\lambda /2$ is fulfilled that allows to form the main lobe of RP without several maximums (see [22] ). In Figure 4, values of the RPD are shown without forming zeros at the distance 0.5 m of array.

For this case, RPD changes from 0.16 W/m^{2} to 0.207 W/m^{2}. The requirement to form zero was stated for the distance
$r=0.5\text{\hspace{0.17em}}\text{m}$ of array at the next points:
$\theta =\pm \text{\pi}/4$ ,
$\phi =\text{\pi}$ . To obtain minimal values of RPD in these points, the field
$U\left({\theta}^{\prime},{\phi}^{\prime}\right)$ in (1) was equated to zero; herewith values of the weight function
$q\left({\theta}^{\prime},{\phi}^{\prime}\right)$ were decreased on one order greater in the comparison of rest values of
${\theta}^{\prime}$ and
${\phi}^{\prime}$ angles.

Such modification of initial data allows to decrease the RPD from quantity of 0.1831 W/m^{2} to quantity of 0.0373 W/m^{2} that is 21.7% of its initial value. In the rest range of
${\theta}^{\prime}$ and
${\phi}^{\prime}$ values, the RPD increased not more than on 10.2% (see Figure 5). The obtained results testify that the field values and its structure outside the considered
${\theta}^{\prime}$ and
${\phi}^{\prime}$ for such r is close to constant value.

In Figure 6, the results related to minimization of the RPD in one point are shown. Similarly to the previous example, the RPR values diminish significantly.

If distance r to array’s plane increases the values of the RPD decrease quickly. The respective results are shown in Figure 7 for the distance that is equal to $r=10.0\text{\hspace{0.17em}}\text{m}$ to array.

The presented results testify that the RPD diminishes in a great extent for such parameters of array; by this its values are more two order lower than for values corresponding to $r=0.5\text{\hspace{0.17em}}\text{m}$ . For example, the maximum of RPD for $r=6.0\text{\hspace{0.17em}}\text{m}$ reaches $1.37\times {10}^{-2}\text{W}/{\text{m}}^{2}$ . It confirms that area of small values of the RPD is placed at distance equal about to 6.0m from the plane of array. This is important data for the creating zeros for the procedure of design of such

Figure 4. The initial distribution of the RPD at $r=0.5\text{\hspace{0.17em}}\text{m}$ .

Figure 5. The RPD with forming zero in two points at $r=0.5\text{\hspace{0.17em}}\text{m}$ .

Figure 6. The RPD with forming zero in one point at $r=0.5\text{\hspace{0.17em}}\text{m}$ .

Figure 7. The RPD at the distance 10.0 m.

type of arrays. The form and values of the main and sidelobes of the array’s RP are varied few in comparison with the case when the synthesis problem is considered without the creating of zeros. The main lobe is wider at level of 0.5 power of the RP on 4.7%, and the level of sidelobes wherein increased on 3.7 dB.

One more requirement to increase the array’s performance consists of the ability to take into account the additional limitations on the amplitude and phase distribution of the currents ${\stackrel{\to}{I}}_{nm}$ in the array elements. This is can be carry out by the prescribing the weight multipliers ${t}_{nm}$ instead of fixed t [28]. For goal to decrease ${\stackrel{\to}{I}}_{nm}$ , the values ${t}_{nm}$ should be increased. The numerical calculations demonstrate that such a way is suitable to avoid the superdirectivity of array’s elements to a large degree. The weight multipliers ${t}_{nm}$ are prescribed so that the amplitude of currents ${\stackrel{\to}{I}}_{nm}$ is close to constant distribution.

8. Conclusions

The variational approach for solving the antenna synthesis problem by the given amplitude RP is modified for the goal to take into account the additional restrictions on the array radiation in the near and far zones. The optimizing functional is formularized in such a way that the deep zeros can be created in the points of near and far zones.

The results of the numerical calculations demonstrate that the level of formed zeros less than −80 dB in the area of far sidelobes of the RP can be reached. By this, the array’s RP in the main lobe varies insignificantly, and difference of the obtained and initial RP’s values does not exceed 3.7 dB. The minimization of the field at the distance about 10.0 m from array does not lead to the considerable variation of the array’s RP.

The future research in the topic under investigation is designated for consideration of vector character of array’s RP and non-transversal character of EM field in a near zone of array, taking into account the mutual coupling of array’s radiators, as well as at development of powerful numerical methods and respective software for modeling the array’s radiation properties in the real physical environment.

Cite this paper

Andriychuk, M. and Kuleshnyk, Y. (2018) Modeling of Plane Arrays Using a Variational Approach.*Open Journal of Antennas and Propagation*, **6**, 93-105. doi: 10.4236/ojapr.2018.64009.

Andriychuk, M. and Kuleshnyk, Y. (2018) Modeling of Plane Arrays Using a Variational Approach.

References

[1] Keiser, K.L. (2000) Electromagnetic Compatibility Handbook. CRC Press, Boca Raton, London, New York, Washington DC.

[2] Hansen, D. (2003) Review of EMC Main Aspects in Fast PLC Including Some History. 2003 IEEE International Symposium on Electromagnetic Compatibility (EMC ‘03), 1, 184-192.

https://doi.org/10.1109/ICSMC2.2003.1428226

[3] Lezynski, P. (2018) Random Modulation in Inverters with Respect to Electromagnetic Compatibility and Power Quality. IEEE Journal of Emerging and Selected Topics in Power Electronics, 6, 782-790.

https://doi.org/10.1109/JESTPE.2017.2787599

[4] Lambrecht, N., Pues, H., De Zutter, D. and Ginste, D.V. (2017) Modeling of Contact Bounce in a Transient Electromagnetic Compatibility Test for the Analysis and Optimization of Nonlinear Devices. IEEE Transactions on Electromagnetic Compatibility, 59, 541-544.

https://doi.org/10.1109/TEMC.2016.2618345

[5] Trainotti, V. (2017) Electromagnetic Compatibility (EMC) Antenna Gain and Factor. IEEE Transactions on Electromagnetic Compatibility, 59, 1006-1015.

https://doi.org/10.1109/TEMC.2016.2642833

[6] Montrose, M.L. (2008) EMC and Propagating Fields. IEEE Antennas and Propagation Magazine, 50, 224-226.

https://doi.org/10.1109/MAP.2008.4653720

[7] Liang, C.-H., Dang, X.-J., Wang, N. and Yuan, H.-B. (2011) Generalized Isolation between Antennas for EMC Problems in Complex EM Environments. IEEE Transactions on Electromagnetic Compatibility, 53, 645-652.

https://doi.org/10.1109/TEMC.2011.2140375

[8] Foged, L.J., Scialacqua, L., Mioc, F., Saccardi, F., Iversen, P.O., Shmidov, L., Braun, R., Quijano, J.L.A. and Vecchi, G. (2013) Echo Suppression by Spatial-Filtering Techniques in Advanced Planar and Spherical Near-Field Antenna Measurements [AMTA Corner]. IEEE Antennas and Propagation Magazine, 55, 235-242.

https://doi.org/10.1109/MAP.2013.6735525

[9] Perez, J.R. and Basterrechea, J. (2007) Comparison of Different Heuristic Optimization Methods for Near-Field Antenna Measurements. IEEE Transactions on Antennas and Propagation, 55, 549-555.

https://doi.org/10.1109/TAP.2007.891508

[10] Bogdanov, F.G., Karkashadze, D.D., Jobava, R.G., Gheonjian, A.L., Yavolovskaya, E.A., Bondarenko, N.G. and Ullrich, C. (2010) Validation of Hybrid MoM Scheme with Included Equivalent Glass Antenna Model for Handling Automotive EMC Problems. IEEE Transactions on Electromagnetic Compatibility, 52, 164-172.

https://doi.org/10.1109/TEMC.2009.2036003

[11] Parrikar, R.P. and Gupta, K.C. (1998) Multiport Network Model for CAD of Electromagnetically Coupled Microstrip Patch Antennas. IEEE Transactions on Antennas and Propagation, 46, 475-483.

https://doi.org/10.1109/8.664110

[12] Foged, L.J., Scialacqua, L., Saccardi, F., Quijano, J.L.A. and Vecchi, G. (2014) Application of the Dual-Equation Equivalent-Current Reconstruction to Electrically Large Structures by Fast Multipole Method Enhancement [AMTA Corner]. IEEE Antennas and Propagation Magazine, 56, 264-273.

https://doi.org/10.1109/MAP.2014.6971966

[13] Ferrieres, X., Parmantier, J.-P., Bertuol, S. and Ruddle, A.R. (2004) Application of a Hybrid Finite Difference/Finite Volume Method to Solve an Automotive EMC Problem. IEEE Transactions on Electromagnetic Compatibility, 46, 624-634.

https://doi.org/10.1109/TEMC.2004.837837

[14] Quijano, J.L.A. and Vecchi, G. (2010) Near-and Very Near-Field Accuracy in 3-D Source Reconstruction. IEEE Antennas and Wireless Propagation Letters, 9, 634-637.

https://doi.org/10.1109/LAWP.2010.2055032

[15] Voitovich, N.N. (1985) Synthesis of Plane Closed Antenna with the Field Restrictions in the Near Zone. Radiotechnika i Elektronika, 30, 458-462. (In Russian)

[16] Schnattinger, G., Mauermayer, R.A.M. and Eibert, T.F. (2014) Monostatic Radar Cross Section Near-Field Far-Field Transformations by Multilevel Plane-Wave Decomposition. IEEE Transactions on Antennas and Propagation, 62, 4259-4268.

https://doi.org/10.1109/TAP.2014.2323429

[17] Galindo-Israel, V. and Rahmat-Samii, Y. (1981) A New Look at Fresnel Field Computation Using the Jacobi-Bessel Series. IEEE Transactions on Antennas and Propagation, 29, 885-898.

https://doi.org/10.1109/TAP.1981.1142680

[18] Andriychuk, M.I., Voitovich, N.N., Savenko, P.A. and Tkachuk, V.P. (1993) Antenna Synthesis According to the Amplitude Radiation Pattern. Naukova Dumka, Kiev. (In Russian)

[19] Salt, H. (1977) Practical Realization of Superdirective Arrays. Radio and Electronic Engineer, 47, 143-156.

https://doi.org/10.1049/ree.1977.0021

[20] Yang, T.C. (2018) Performance Analysis of Superdirectivity of Circular Arrays and Implications for Sonar Systems. IEEE Journal of Oceanic Engineering, 43, 1-11.

https://doi.org/10.1109/JOE.2018.2801144

[21] Diao, J. and Warnick, K.F. (2018) Practical Superdirectivity with Resonant Screened Apertures Motivated by a Poynting Streamlines Analysis. IEEE Transactions on Antennas and Propagation, 66, 432-437.

https://doi.org/10.1109/TAP.2017.2772929

[22] Balanis, C.A. (1997) Antenna Theory: Analysis and Design. John Wiley & Sons, New York, Chichester, Brisbane, Toronto, Singapore.

[23] Markov, G.T., Petrov, B.M. and Grudinskaya, G.P. (1979) Electrodinamics and Propagation of Radiowaves. Sov. Radio, Moscow. (In Russian)

[24] Andriychuk, M.I. (2000) Investigation of Solution of the Nonlinear Synthesis Problem for the Waveguide Array. International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory, Tbilisi, 47-51.

[25] Kotelnikov, V.A. (2006) On the Capacity of “Ether” and Wire in Electro-Communication (Appendix). UFN, 176, 762-770.

https://doi.org/10.3367/UFNr.0176.200607h.0762

[26] Andriychuk, M.I. and Savenko, P.O. (2000) Synthesis of a Waveguide Array with Due Regard for the Mutual Coupling of Radiarors. International Conference on Mathematical Methods in Electromagnetic Theory, Kharkiv, 12-15 September 2000, Vol. 2, 604-606.

[27] Savenko, P.O. and Anokhin, V.J. (1997) Synthesis of Amplitude-Phase Distribution and Shape of Plane Antenna Aperture for a Given Power Pattern. IEEE Transactions on Antennas and Propagation, 45, 744-747.

https://doi.org/10.1109/8.564102

[28] Andriychuk, M.I. and Voitovich, N.N. (2013) Antenna Synthesis According to Power Radiation Parttern with Condition of Norm Equality. International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory, Lviv, September 23-26 2013, 137-140.