Theoretical Electrical Engineering and Electrophysics
UDC 621.315 doi: 10.20998/2074-272X.2019.5.05
M.M. Rezinkina
CALCULATION OF ELECTROMAGNETIC FIELDS IN INHOMOGENEOUS MEDIA FOR SELECTION OF PROTECTIVE COATINGS
The aim of the work is to develop a method for high-frequency EMFs calculation in terms of the modified vector potential and to use it for choice of the parameters of the layered absorbing coatings that provide effective absorption of radio wave radiation. Methodology. A method for calculation of electromagnetic fields in inhomogeneous media is described. It is performed in terms of the modified magnetic potential at presented a plane electromagnetic wave as an incident and reflected one. To reduce the computational domain size, additional uniaxially perfectly matched layers ensuring fast and non-reflective attenuation of electromagnetic field (EMF) strengths are used at its boundaries. It is supposed that electrical parameters of such layers have not only resistive, but also magnetic losses. Results. Modification of the known methods for dividing the wave into the incident and reflected, as well as introducing uniaxially perfectly matched layers at the domain boundaries and formulation of the problem in terms of the modified vector potential, provide a possibility to solve the problem of the electromagnetic wave incidence onto heterogeneous media at reduced memory and counting time compared to the formulation tasks through the EMF strengths. Practical value. The described method was used to select parameters of radio-absorbing coatings, application of which on conductive bodies reduces the reflections from them of the incident electromagnetic waves. This provides a possibility to choose the electrical parameters of layered coatings with active losses, as well as magnetic losses, which provide a decrease of the reflected waves amplitude by an order of magnitude or more. Originality. For the first time, the problem of calculating the incidence of an electromagnetic wave on an inhomogeneous medium with active electrical and magnetic losses was formulated in terms of the modified vector potential. The parameters of a multilayer lossy coating, which ensures decrease of the reflected EMF wave strengths by an order of magnitude or more are determined. References 11, figures 3.
Key words: electromagnetic field, inhomogeneous media, modified vector magnetic potential, finite integration method, uniaxially perfectly matched layers.
Описаний метод розрахунку електромагттних потв в neodnopidnux середовищах, виконаний в термнах модифжованого магнтного потенцалу при подл плоскоТ електромагжтноТ хвилг на падаючу i вiдбuту. Для зменшення габаритгв розрахунковоТ областг на и границях введет дoдаткoвi одновкне iдeальno пoглunаючi шари, що забезпечують швидке i бeзвiдбuвne загасання напруженостей електромагштних пoлiв (ЕМП). Особлив^тю таких шаpiв е вuбip х електричних nаpамeтpiв такими, щоб забезпечити наявшсть в них не ттьки резистивних, але i магнiтнux втрат Математичне моделювання процеав при падшн ЕМП на npoвiднe середовище, перед яким помщен покриття з noдiбнuмu властивостями показало мoжлuвiсть ефективного загасання в них ЕМП радючастотного дiаnазoну. В результата проведених дослгджень обран параметри покриттгв, що забезпечують мтшальне вiдбuття падаючих eлeктpoмагнiтнuxхвиль. Бiбл. 11, рис. 3.
Ключовi слова: електромагнтне поле, нeoднopiднi середовища, модифшований векторний маг^тний noтeнцiал, метод кнцевого ттегрування, iдeальнo узгоджен noглuнаючi граничн шари.
Описан метод расчета электромагнитных полей в неоднородных средах, выполненный в терминах модифицированного магнитного потенциала при разделении плоской электромагнитной волны на падающую и отраженную. Для уменьшения габаритов расчетной области на ее границах введены дополнительные одноосно идеально поглощающие слои, обеспечивающие быстрое и безотражательное затухание напряженностей электромагнитных полей (ЭМП). Особенностью таких слоев является то, что их электрические параметры выбраны такими, чтобы обеспечить наличие в них не только резистивных, но и магнитных потерь. Математическое моделирование процессов при падении ЭМП на проводящую среду, перед которой помещены покрытия с подобными свойствами, показало возможность эффективного затухания в них ЭМП радиочастотного диапазона. В результате проведенных исследований выбраны параметры покрытий, обеспечивающих минимальное отражение падающих электромагнитных волн. Библ. 11, рис. 3.
Ключевые слова: электромагнитное поле, неоднородные среды, модифицированный векторный магнитный потенциал, метод конечного интегрирования, идеально согласованные поглощающие граничные слои.
Introduction. When solving a number of problems of theoretical electrical engineering and electrophysics, the problems of calculating electromagnetic fields (EMF) arise, including when a wave is incident on inhomogeneous media with flat interfaces. Such problems have to be solved, for example, when choosing the parameters of protective radar absorbing coatings of extended conductive objects from high-frequency radiation. Since in this case the application of analytical methods [1] is impossible, the development of numerical methods is necessary. Taking into account the flat geometry of the systems under consideration, the use of finite-difference methods seems most effective (see, for example, [2, 3]). To reduce the order of the system of equations being solved, the problem formulation in terms of the modified vector magnetic potential can be used [4]. In contrast to the traditional formulation of electromagnetic fields through EMF strengths, this approach allows to reduce the required computer resources due to the fact that in the calculation 3 components of the modified vector magnetic
potential are unknown, rather than 3 components of the electric field strength and 3 components of the magnetic field strength [4].
If there are several media interfaces with flat interfaces, it seems efficient to use the numerical finite integration method (FIM) [5]. The essence of this method is to impose a rectangular computational mesh on the considered region and integrate Maxwell equations over the volumes or faces of its cells. Thus, equations to be solved are obtained using conservation laws. Moreover, the formulation of the calculation problem using the finite integration method in terms of a modified vector magnetic potential with a special choice of the computational mesh ensures the automatic fulfillment of conditions at the media interfaces for the vectors of strengths and the flux density of the EMF [6, 7].
To limit the calculation region when determining the EMF in open areas when formulating the problem using EMF strengths, the so-called UPML (uniaxially perfectly
© M.M. Rezinkina
matched layers) are used [8]. To solve the problems associated with the propagation of high-frequency oscillations in inhomogeneous media with losses, it seems important to develop UPML in relation to the formulation of problems in terms of a modified vector magnetic potential.
The goal of the work is to develop a method for calculating high-frequency EMFs in lossy media using modified vector potential and UPML.
Definition of the problem of calculating the process of incidence of a plane electromagnetic wave on layered media. In order to reduce the number of equations in solving the problem of calculating EMF, we use the so-called modified vector magnetic potential A [4].
In the absence of external sources of electric charge, an additional condition when using the scalar electric (<) and the vector magnetic (a) potentials, instead of electric (E) and magnetic (H) field strengths, the gauging condition < = 0 can be chosen (see, for example, [1]). Then vectors of the magnetic flux density B and the electric field strength E are expressed through A in the form:
B = rot A = UU0 H ; (1)
e = -SA.
St
(2)
In order for the reflectionless EMF attenuation to occur upon falling onto the radar absorbing coating, it is necessary to ensure the presence of losses not only related to the active conductivity of the medium, but also magnetic losses.
The second Maxwell equation, taking into account magnetic losses, has the form [9, p. 69]:
„ , SB r0te = + 7^H
(3)
where 7 are the equivalent magnetic losses [Q/m].
Taking into account (1), rewrite the right-hand side of (3) in the form:
dB dA 7u
— + 7 UH = rot— rota .
dt dt uUo
Then write (3) as
rote = -rot
SA Y m
rotA,
dt UUo
whence, for piecewise homogeneous media, E can be expressed in terms of A as follows:
e = -—-y- a . St MMo Here, the first Maxwell equation
SD
(4)
rot H = vE+-,
dt
where D is the electric field induction; 7 is the electrical conductance takes the form:
S 2 A
rota SA
rot-= -Y^7-e0e—2
MMo St St2
YY, A -Y^SA . (5)
MMo MMo St To obtain a numerical solution for the vector potential, the considered region was divided into parallelepiped cells (Fig. 1) so that the nodes of the computational mesh (i, j, k) lay at the interfaces of the media. It was assumed that the electrical properties of the medium within each of the cells are homogeneous.
v>z Mz,jz,kz
Fig. 1. Cell of the calculation scheme
The equation to be solved was obtained by integrating (5) over sections perpendicular to the coordinate axes. Then the equations for the components of the vector A can be written as a result of integration of (5) over the sections Sx, Sy, Sz (see Fig. 1) and using the Stokes theorem. So for Ay - the y-th component of the vector A - we obtain:
rota
MMo
f SAy
dl = J (-y-^
S 2 Ay
St
-e0e-
St2
(6)
YYM Yueoe SAV M-AV -!-M0--rV)dS,
UUo " UUo dt where ly is the contour covering the area Sy.
As a result of integration of (5) over the corresponding areas: Sx for Ax and Sz for Az (see Fig. 1), similar expressions are written for the remaining components of the vector potential. Passing to the difference form of writing, we obtain a system of equations to be solved.
Features of the numerical calculation of EMF using a modified vector magnetic potential. When calculating electromagnetic processes when a plane electromagnetic wave having a finite length and harmonic filling is incident on objects, the problem arises of setting conditions for the strengths and potential of the EMP at the boundaries of the computational domain. To limit the computational domain when applying finite-difference methods of numerical computation, the introduction of so-called uniaxial perfectly matched layers (UPML) at the external boundaries of the computational domain is usually used, which ensure fast and reflectionless attenuation of the EMF. However, in this case, for the problem formulation under consideration, it is impossible to specify zero boundary conditions for the EMF potentials and strengths at the external UPML boundaries. The fact is that in addition to the reflected wave, which can really be assumed to be equal to zero at a sufficient distance from the object causing it, there is also an incident wave, which over time can be non-equal to zero and at a sufficient distance from the object. To solve this problem, a representation of the electromagnetic wave in the form of the sum of the incident and reflected waves is used. The propagation of both of these waves is described by the same Maxwell equations. In this case, the computational domain is also divided into two zones: the first (I), the inner, in which the studied object is located, and the second (II), the outer, located between the boundaries of the 1st region and the UPML inner borders. Here, in zone I, the propagation of the total wave is calculated, which contains both the incident and reflected waves. In zone II, only the propagation of the reflected wave is calculated. At the boundaries of zones I and II, the transition from calculating
J
S
V
V
only the reflected and total waves is taken into account. At the outer boundaries of UPML, zero conditions can be used for the potentials and strengths of the EMF, since they relate only to the reflected wave.
In [9], the solution of the problem of propagation and reflection of electromagnetic waves is presented using the described approach when using the Yee algorithm to calculate the components of the electric and magnetic field strengths. This approach was applied to solve the problem formulated in terms of a modified vector potential using the finite integration method. In this case, terms appearing in difference schemes containing the known values of the vector potential corresponding to the incident wave. The addition or subtraction of these terms when writing the equations to be solved for the nodes adjacent to the interface of subregions I and II allowed to implement the approach described above. Since the values of the incident wave at each moment of time and for each node of the calculation scheme are known, they are transferred to the right side of the equations. Such terms are added to the right sides of the difference equations written for nodes lying on the boundaries of subregions I and II, and are subtracted from the right sides of difference equations written for nodes located in II in planes in front of the boundaries of subregions I and II.
The introduction of UPML in calculating the propagation of an electromagnetic wave in terms of a modified vector magnetic potential requires a different approach than when solving through electric and magnetic field strengths. This is due to the fact that the parameter sM
s— = 1 +
r—
j®MMo
which is chosen equal to sS
ss= 1 + -
r
(7)
(8)
r
jmeeo
rM
(9)
(1o)
M
j®MMo
Then from (9, 10) we obtain:
rM = r
MMo eeo
Hence the relation between y and y^:
« 0.142-106r— /m]. (11)
SSo S
We obtain expressions for EMF in UPML layers by writing (6) through complexes (indicated by an underscore of the corresponding value) and taking into account (7), (8):
rota „ iv . . „ . n
—=dl = J (-j®râv - (j®) SosA
MoM
y
S„
rr—-Av - jv^-r—ây)dS =
Mo-
Mo M
= J [-(J®)2sosây (1 +
Y
Sv
jmsos
)-
- ja—y^Ay (1+-JL~)\dS = M00M j®e0e
= - j (ja)2e0eAySESMdS.
Sy
Dividing both parts of the resulting expression into Sp, we write:
j rot~ dl = - j (ja)2e0eAvssdS .
I MMm S
y y
In view of (9), (10), we transform the last expression to the form:
J2®°Ali = -J j®rAydS
L
rM
s„
yy
Then, in the last expression, we move from the complex form of writing to the functions of time:
■rot[5a / dt ]
dl = -
dA
rm
r-
dt
ds.
(12)
jaee0
when formulating the problem through A appears in the denominator of the equation to be solved. Therefore, its conversion during the transition from the frequency domain to the time domain is impossible. To use UpML in such a formulation, it is proposed to use a medium in which s-m and pm are much smaller than y and respectively, as absorbing layers. To ensure the constancy of the propagation speed of electromagnetic waves and the absence of reflections when the EMF falls on the UPML layers, it is necessary sM = se. At ^-m << y^ and s-m<<y:
A numerical solution was obtained by writing equations of the form (6), (12) for nodes of the computational mesh and replacing the differential operators in them by difference analogues. To solve the obtained system of equations, the sweep method was used (for more details see [6, 7]).
To evaluate the effectiveness of the proposed method of introducing UPML, numerical experiments were carried out to calculate electromagnetic processes in the case of a plane electromagnetic wave incident on a perfectly conducting sheet. It was assumed that the EMF with |£"|=1 V/m from time t = 0 falls in the direction of the Oz axis onto a conductive sheet whose thickness is d = 3-A (where A is the spatial step). It is believed that this sheet is located in the z = const plane and has dimensions in the x and y coordinates many orders of magnitude greater than the wavelength of the incident EMF. When calculating at the boundaries of the computational domain in the directions of the Ox, Oy, and Oz axes zero conditions were used for the vector potential of the reflected waves. The frequency of the incident EMF was assumed to be 10 GHz. In the calculations, the time step was set equal to At = 77200, A = c-77200 (c = 3-108 m/s). The UPML parameters were as follows: the number of layers N = 10, y varied according to the exponential law with an exponent of 3 [8] from 35 S/m at the inner UPML boundary, y^ varied depending on y in accordance with (11). The dimensions of the computational domain in the direction of the Oz axis are Zmax = 0.5-A (where X is the wavelength), and the number of steps in reflection zone II is 10.
Figure 2 shows the calculated time sweeps of the reflected (scattered) electric field strength |^scat| when a plane electromagnetic wave falls on a conductive sheet.
v
S
v
v
s
s
As can be seen from the Figure, due to the fact that UPML layers with the parameters described above are located at the boundaries of the computational domain, the relative error in calculating the strength |escat| does not exceed 3%, despite a significant truncation of the length of the computational domain.
Thus, a technique for calculating the propagation of an electromagnetic wave in media with magnetic losses in terms of a modified vector magnetic potential is described. This ensures the introduction of UPML at the boundaries of the computational domain, which was previously possible only when formulating the problem through EMF strengths. i
O.i 0.60.4 0.2 0-0.2 -0.4 -O.i -O.i
E, V/m A f\f \ f N ■ ! \ A \ ! \ / ^ !
j V i
/i,AA ' i ■■ 1/ f. 1 \ 1 \ I \ t A \ ■ : \ \ ; \ \ > / '' I'! / ' \ i \i f\ ! ' I K { y 1 . \ / \ / \ / \ / ^
i <1
1 i i: 1> r A i ' 1 \ ' / ■ i r I / \1 ■ h< i A i 1 !. 1 ; i i 1 ! \ 1 1 ' ' w \ \ \w \ i ; i p i ; \ ;
+ V / v JMf- \ ( V \ -4-
; \ T1—MT \ A i *I0A(-
♦10- (-10) *io"(-io ) 3*I0A(-10 )
2. The proposed calculation method is the theoretical basis for developing a method for selecting parameters of radar absorbing coatings with electric and magnetic losses, the application of which to conductive objects provides a significant reduction in the amplitude of the reflected radio waves.
Fig. 2. Calculated dependencies |Escat|(t) for nodes located at
different distances (z) from the conducting sheet: z = X (-),
-z = 0.5-X (---), -z = 0 (-----)
Mathematical modelling of EMF fall on layered coatings with dielectric and magnetic losses. In some cases, various dielectrics with low values of the reflection coefficient and high values of the absorption coefficient of radio waves are used as protective coatings from high-frequency radiation. Scopes of radar absorbing coatings are various. So, they are used to create screens for radio receivers, to protect biological objects from electromagnetic radiation, to equip special research chambers, etc. A large number of works have been devoted to this issue (see, for example, [10, 11]), but the problem of choosing coating parameters remains relevant, since the question of a significant decrease in the amplitude of waves reflected from conductive objects has not been completely resolved.
Figure 3 shows the calculated time sweeps of the electric field strength when a plane electromagnetic wave falls on a conductive sheet, before which a radar absorbing coating is placed with the following parameters: the 1st layer (closest to the conductive sheet) - e = f = 300, y = 315 S/m, the next layer - e = f = 100, y = 275 S/m, the last layer - e = f = 60, y = 17.5 S/m (yf changed depending on y in accordance with (11)). As can be seen from the Figure, due to the attenuation of the electromagnetic wave in the coating layers, the reflection from the conducting sheet in the steady state does not exceed 10% of the initial reflection level.
Conclusions.
1. A method for calculating the propagation of an electromagnetic wave in lossy media located in open areas at formulating the problem in terms of a modified vector magnetic potential and using UPML the implementation of which allows to reduce the memory and calculation time compared to formulation of the problem through EMF strengths is proposed and tested by comparison with an analytical solution.
How to cite this article:
*10*{-I0) 2*I0A(-10) 3*I0A(-10) 4* I CX-Fig. 3. Calculated dependencies |Escat|(t) in the presence in front of the conductive sheet of the radar absorbing coating for nodes located at different distances (z) from the conductive sheet: z = X (-), -z = 0.5-X (---), -z = 0 (-----)
REFERENCES
1. Stratton J.A. Electromagnetic theory. NJ, IEEE Press, 2007. 614 p.
2. Yee K.S., Chen Jei Shuan, Chang A.H. Conformal finite-different time-domain (FDTD) with overlapping grids. IEEE Transactions on Antennas and Propagation, 1992, vol.40, no.9, pp. 1068-1075. doi: 10.1109/8.166532.
3. Werner D.H., Mittra R. Frontiers in electromagnetics. New York, IEEE Press, 1999. 876 p.
4. Biro O., Preis K. On the use of the magnetic vector potential in the finite-element analysis of three-dimensional eddy currents. IEEE Transactions on Magnetics, 1989, vol.25, no.4, pp. 3145-3159. doi: 10.1109/20.34388.
5. Clemens M., Weiland T. Discrete electromagnetism with the finite integration technique. Progress in Electromagnetics Research, 2001, vol.32, pp. 65-87. doi: 10.2528/PIER00080103.
6. Rezinkina M.M., Rezinkin O.L. Modeling of the electromagnetic wavefront sharpening in a nonlinear dielectric. Technical Physics, 2011, vol.56, iss.3, pp. 406-412. doi: 10.1134/S1063784211030169.
ZRezinkina M.M. Modeling of the dendrite shape variation with applied electric field strength in poly(ethylene). Technical Physics Letters, 2000, vol.26, iss.3, pp. 196-198. doi: 10.1134/1.1262789.
8. Berenger J.-P. A perfectly matched layer for the absorption of electromagnetic waves. Journal of Computational Physics, 1994, vol.114, no.2, pp. 185-200. doi: 10.1006/jcph.1994.1159.
9. Taflove A., Hagness S.C. Computational Electrodynamics: the Finite-Difference Time-Domain Method. Boston - London, Artech House, 2000. 852 p.
10. Zhou Z., Chen K., Zhu B., Zhao J., Feng Y., Li Y. UltraWideband Microwave Absorption by Design and Optimization of Metasurface Salisbury Screen. IEEE Access, 2018, vol.6, pp. 26843-26853. doi: 10.1109/access.2018.2835815.
11. Bottauscio O., Chiampi M., Manzin A. Numerical analysis of magnetic shielding efficiency of multilayered screens. IEEE Transactions on Magnetics, 2004, vol.40, iss.2, pp. 726-729. doi: 10.1109/tmag.2004.825171.
Received 12.04.2019
M.M. Rezinkina, Doctor of Technical Science, Professor, 1 National Technical University «Kharkiv Polytechnic Institute», 2, Kyrpychova Str., Kharkiv, 61002, Ukraine, e-mail: [email protected]
Rezinkina M.M. Calculation of electromagnetic fields in inhomogeneous media for selection of protective coatings. Electrical engineering & electromechanics, 2019, no.5, pp. 31-34. doi: 10.20998/2074-272X.2019.5.05.