Charge method in mathematical modeling of printed radiators

Published: 08/18/2010
Original: Computer-aided design of microwave devices and systems // Interuniversity digest of scientific papers. (Moscow, MIREA), 1982, p.p.138…148
© 1982, V. S. Filippov, A. A. Sapozhnikov
© 2010, EDS–Soft,,   E-mail:

The method of researching the electrodynamic characteristics of the PAA radiating curtain described in this work applies to antennas with rectangular printed radiators and also other radiators that can be presented as an array of rectangular items. However, the charge method is applied to researching a printed vibrator radiator (fig.1) being part of an antenna array in this particular example.

As it is known [1], some freedom in choosing the integral representation of electromagnetic fields in diffraction problems makes it possible to get various integral equations. The essence of the proposed method comes down to determining the intermediate characteristic — the charge distribution over the metal plate and then by means of integration — the electric current. The vectors of the electromagnetic field are determined via the vector potential of electric currents using the known equations:


where .

The radiator is regarded as part of an infinite periodic antenna array, which makes it possible to proceed to the analysis of the field within one period. Due to this, it is reasonable to use a representation in the form of expansion in plane waves for the vector potential [2]:



and is the vector of the electric current volume density, (x,y,z), (x',y',z') are the coordinates of the point of observation and the point of integration respectively.


Presuming that the metal radiators are ideal conductors and they are also infinitely thin, the vector of the electric current volume density can be written down in the following way:


and it is reasonable to replace the effect of the screen with a mirror view of the radiator:


where is the vector of the electric charge surface density, is the Dirac delta function.

By integrating expression (2) using the axial coordinate z and taking into account (3) and (4), it is possible to get the values of the vector potential harmonic in two characteristic areas of the structure:


The wave of the potential is presented in the form of the superposition of waves related to Е- and Н-waves of the spatial waveguide:


The specified waves are determined by the following expressions:


where .

The above relations are obtained for the case of homogeneous space above the array (z>0). Now taking into account the influence of obstacles in the form of the border between dielectric layers, it is possible to transform expression (6) in the following way for the area z0<z<z1, where z0 and z1 are distances from the screen to the metal plate and to the upper dielectric layer respectively:


In this relation the quantities:


are the functions of the factors of the reflection of Е- and H-waves from the dielectric interface z=z1. The reflector factors can be derived from a system of equations corresponding to the boundary conditions on the dielectric interface:



and , are the conductivities of the (m,nth harmonic of the E- and H-wave respectively.

The component of the vector quantities in (8) can be obtained by substituting (5) into (7):




and the symbol means the X or Y coordinate.

Integrating the last expression by parts and taking into account the boundary conditions for the normal components of the current on the radiator edges and the continuity equation:


we will get


where are the components of the electric charge surface density related to the corresponding components of the current density.

Thus, the field vectors (1) expressed via electrodynamic potentials are completely determined by the so far unknown distribution of the electric charge over the radiator surface. The required charge can be found with the help of the integral equation method. To do it, we use a boundary condition for the complete electric field on the surface of an ideally conducting radiator:



is an extraneous electric field and E0 is the known value of the electric field intensity in the gap between vibrator arms. By integrating the first equation from (1) using transverse coordinates, we can get:


where С is the integration constant.

By using the expressions found above in (16) for electric potentials, it is possible to get a system of two integral equations of the first type:


The presence of two components of the charge actually requires the solution of the system of equations (17). In the particular case when the radiators are narrow vibrators with one nonzero component of the current whose transverse distribution is known, the system of equations (17) is reduced to one-dimensional integral equation:





The kernel (19) of the integral equation (18) is presented as a sum of two summands (20). In case the arguments coincide, the summand has an integrable singularity and, if we isolate it [3], we can reduce the integral equation of the first type (18) to the equation of the second type:



the numerical solution of which is a well-defined problem [4]. In the last expression, the summand of the kernel:


is a smooth function of position. Where R is the radius of a circle with its center in the singular point x=x', and

The orthogonal function series expansion coefficients for the isolated singularity of the type can be found in the following form:


where , , are zero-order and first-order Bessel functions, , are zero-order and first-order Struve functions respectively.

The function in the expression (22) is a results of isolating the singularity:


Thus, we get the integral equation (22) and the algorithms of its numerical solution are well known. In this problem, it is convenient to use the Krylov-Bogolyubov interpolation method according to which the integral equation is reduced to a system of linear algebraic equations in case of the piecewise constant approximation of the unknown function. In order to do it, the vibrator is divided into N finite elements that go along only axis X because of the above mentioned small vibrator width. If we consider the unknown function is constant within each finite element and coordinate the solution in their middle points, we can come to the system of equations:



Solving this system of equations, we can find the electric charge distribution over the metallic vibrator that is used to find all integral characteristics of a radiator in the array: the partial radiation pattern, the input impedance and others. In case the finite elements are not very small, the matrix of the system will be well-conditioned because the isolation of the singularity of the integral equation kernel leads to the absolute domination of diagonal elements over the rest of the elements in the matrix of the system.

It should be noted that the characteristics of a radiator are determined in this method via the charge distribution over the vibrator surface, instead of current as it is done in many related problems. The main advantage of the method comes down to the faster convergence of the necessary radiator characteristics in case of the specified calculation accuracy. That is due to the fact that the polynomial approximant selected for describing the charge distribution is equivalent to the polynomial for the current whose power is greater by one. It is known that the integral characteristics of radiators have better convergence as compared to the current distribution in case of the numerical solution of boundary value problems. In this case, the current distribution itself is the integral characteristic of the charge distribution, which determines the advantages of the method.

The developed methods of analyzing antenna arrays with printed radiators serve as the basis for creating the program for calculating the characteristics of a printed-circuit vibrator-type PAA. In fig.2 you can see the results of the numerical solution of the test problem for determining the electric charge on the vibrator of an array radiator with the following parameters: dx= 0,6, dy= 0,25, a= 0,5, b= 0,03, c= 0,01, z0= 0,15, ε1= 4, ε2=ε3= 1, φ= 900, θ= 50.




1. E.N. Vasiliev Algorithmization of Diffraction Problems Based on Integral Equations. - In Applied Electrodynamics Moscow: Higher School, 1977, #1.
2. G.T. Markov Antennas. – Moscow/Leningrad: Gosenergoizdat, 1960.
3. A.S. Ilyinsky, A.G. Sveshnikov Numerical Methods in Diffraction Theory. Moscow: MGU, 1975.
4. A.N. Tikhonov, V.Ya. Arsenin Solution Methods for Ill-Defined Problems. Moscow: Science, 1979.