Journal of Electromagnetic Analysis and Applications
Vol.10 No.02(2018), Article ID:82351,21 pages
10.4236/jemaa.2018.102002

Non-Monochromatic Electromagnetic Radiation of Inhomogeneous Media

Nicolay N. Grinchik1, Vladimir G. Andrejev2, Galina M. Zayats3, Yury N. Grinchik4

1A. V. Luikov Heat and Mass Transfer Institute of the National Academy of Sciences of Belarus, Minsk, Belarus

2Ryazan State Radio Engineering University, Ryazan, Russia

3Institute of Mathematics of the National Academy of Sciences of Belarus, Minsk, Belarus

4Minsk, Belarus

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 1, 2017; Accepted: February 6, 2018; Published: February 9, 2018

ABSTRACT

A coordinated physicomathematical model for the propagation of a soliton-like electromagnetic pulse in a heterogeneous medium is developed in the presence of strong discontinuities in the electromagnetic field. The model is based on the reduction of Maxwell’s equations to the well-studied wave equation. When the electromagnetic pulse was specified, its amplitude modulation was taken into account, as was the nonstationary broadening of the spectral line. Conditions for matching the momentum for the first initial boundary-value problem are obtained. The time dispersion of the electrical induction is taken into account in terms of the function of signal conditioning which takes account of the broadening of its spectral line and integration over the continuous spectrum. With this approach, it is not necessary to neglect spatial derivatives, and also to use spatial nonlocal relations to take account of the effect of surface charge, surface current, and spatial dispersion of electrical induction at the interfaces of adjacent media.

Keywords:

Strong Discontinuities in the Electromagnetic Field, Nonstationary Broadening of the Spectral Line, Surface Charge, Surface Current, Electromagnetic Wave, Heterogeneous Medium, Soliton-Like Electromagnetic Pulse, Induced Surface Charge

1. Introduction

In many fields of technology, for example, radar and biophysics, it is of interest to study signal propagation in a heterogeneous medium. Real pulse signals in most cases are not strictly monochromatic and are characterized by the width of spectral lines. In this case, in a heterogeneous medium, at the interfaces of adjacent media, as a rule, the function of the first kind is discontinuous, i.e., a strong discontinuity in the electromagnetic field.

In radio-physics, one often uses, for the problem of propagation of electromagnetic waves in the presence of “perfectly conducting” bodies, a surface impedance of the form as the boundary condition in determining outside the conductor (M.A. Leontovich, 1948); the fulfillment of the condition of continuity of the tangential components of the electric and magnetic field strength, i.e., the absence of a surface current, is assumed. Not only do the induced surface charge and current character the properties of the surface, but they are also functions of the process; therefore, surface impedances hold for conditions under which they have been determined. Under other experimental conditions, one cannot use surface impedances, since surface impedances depend on the structure of the determined field and are not boundary conditions in the ordinary sense. When the notion “surface impedance” is used there are a number of constraints and assumptions: the depth of penetration of an electromagnetic wave into a medium and the wavelength in it must be small compared to the wavelength in the ambient space, compared to the distances from the field’s sources, and compared to the radii of curvature of the body’s surface; the changes in the relative permittivity along the body’s length and in the magnetic permeability at distance of the order of the wavelength (or on a length equal to the penetration depth) are small. Furthermore, if we have no ware absorption or it is minor, it is difficult to use the existing impedance conditions, since in the body, there can be waves traveling not only from the surface but from the body as well. By virtue of the aforesaid, the scattering of electromagnetic waves by a group of closely spaced bodies assumes the use of exact boundary conditions with introduction of a surface charge and a surface current.

2. Statement of the Problem

Let us consider the interface between two media with different electrophysical properties. On both sides, the vectors of magnetic field strength and magnetic induction, as well as electric field and electric displacement vectors are finite and continuous, but on the surface S they can undergo a discontinuity of the first kind. Furthermore, induced surface charges s and surface currents i (vectors lying in the tangent plane to the surface S) arise at the interface under the action of an external electric field.

When studying the electric field interacting with the material medium, we use Maxwell’s equations (1857)

j total = × H , D = ρ , (1)

B t = × E , B = 0, (2)

where j total = λ E + D t , B = μ μ 0 H , and D = ε ε 0 E .

The system of equations is supplemented, on the interface S, with the boundary conditions

D n 1 D n 2 = σ , (3)

E τ 1 E τ 2 = 0 , (4)

B n 1 B n 2 = 0 , (5)

H τ 1 H τ 2 = [ i τ n ] . (6)

The subscripts n and t denote the normal and tangential components of the vectors to the surface S, and the subscripts 1 and 2, the adjacent media with different electrophysical properties. By t we mean any direction tangential to the discontinuity surface. At the same time, there is no closing relation for the induced surface charge s, which leads to the necessity of introducing an impedance matrix, which is determined experimentally or in some cases theoretically from quantum representations.

Not only does the induced surface charge s characterize the properties of the surface, but it is also is a function of the process, i.e., s ( E ( E / t , H ( H / t ) ) ) . Therefore the surface impedances are valid for conditions under which they were determined. Under other experimental conditions, it is impossible to use surface impedances. The problem of determining the surface charge and the surface current at the interfaces is further complicated by the study and modeling of nonstationary electrochemical processes and of pulsed electrolysis, since the surface charge is influenced by heat and mass fluxes. Let us show that s can be calculated from the phenomenological macroscopic equations of the Maxwell electromagnetic field and the law of conservation of electric charge taking account of the specific nature of the interface of adjacent media.

The Maxwell Equations (1), (2) represent a system of first-order differential equations containing unknown components of the field E ( x , y , z ; t ) and H ( x , y , z ; t ) , where the spatial coordinates x , y , z and the time parameter t vary within the ranges ( x , y , z ) D R 3 and 0 t . However, to determine the unknown vector functions in a rapidly varying field, it is expedient to exclude one field from Equations (1), (2) and to go to higher-order equations containing only the unknown field E or H . This is due to the fact that the resulting equations in many cases are classical and well-studied.

We assume that the parameters of the medium e, m and g depend only on the coordinates x , y , z and have sufficient smoothness, i.e., we consider piecewise homogeneous media.

A Generalized Wave Equation for E and Conditions at the Interfaces in the Presence of Strong Discontinuities of the Electromagnetic Field. Formulation of the Physicomathematical Model of Propagation of the Electromagnetic Field in a Layered Medium

We multiply the left and right sides of the equation for the total current (1) by μ μ 0 and differentiate it with respect to time. We apply the curl operator to the left and right sides of Equation (2) and, setting μ 1 -const, obtain

j total t = 1 μ μ 0 2 E 1 μ μ 0 grad ( div E ) . (7)

In Cartesian coordinates (7), will have the form

j total x t = 1 μ μ 0 ( 2 E x x 2 + 2 E x y 2 + 2 E x z 2 ) 1 μ μ 0 x ( E x x + E y y + E z z ) , (8)

j total y t = 1 μ μ 0 ( 2 E y x 2 + 2 E y y 2 + 2 E y z 2 ) 1 μ μ 0 y ( E x x + E y y + E z z ) , (9)

j total z t = 1 μ μ 0 ( 2 E z x 2 + 2 E z y 2 + 2 E z z 2 ) 1 μ μ 0 z ( E x x + E y y + E z z ) . (10)

On the interface, there is also the relation

div i τ + I q x 1 I q x 2 = σ t , (11)

which reflects the law of conservation of charge at the interface of adjacent media. Equations (3)-(6) can be written in the Cartesian coordinate system:

D x 1 D x 2 = σ , (12)

E y 1 E y 2 = 0 , (13)

E z 1 E z 2 = 0 , (14)

B x 1 B x 2 = 0 , (15)

H y 1 H y 2 = i z , (16)

H z 1 H z 2 = i y . (17)

where i τ = i y j + i z k is the surface current density, while the x coordinate is directed along the normal to the interface. By the densities i y , i z of surface currents, we mean the quantity of electricity flowing per unit time through a unit length of a segment located on the surface along which the current flows and perpendicular to the direction of the current.

The order of the system of differential Equations (8)-(10) is 18. Therefore, at the interface S, generally speaking, nine boundary conditions must be specified. In addition, three more conditions (12), (16), and (17) containing unknown (before the solution) quantities must be satisfied at this boundary. Therefore, the total number of conjugation conditions at the interface S must be equal to 12 for correct solution of the problem.

Differentiating the expression of Equation (12) with respect to time and taking into account relation (11), at the interface, we obtain the condition for the normal components of the total current:

div i τ + j total x 1 = j total x 2 , (18)

which makes it possible to exclude from consideration the surface charge density.

Due to the equality of tangential projections of the electric field in z and y, according to Equations (13), (14), the expressions for the surface current densities i z and i y have the form

i z = λ ¯ E z | x = ξ , i y = λ ¯ E y | x = ξ , (19)

where

λ ¯ = 1 2 ( λ 1 + λ 2 ) | x = ξ (20)

is the mean value of electrical conductivity at the interfaces z of adjacent media in accordance with the Dirichlet theorem for a piecewise smooth and piecewise- differentiable function. Consequently, it follows from Formulas (19) and (20) that

[ div i τ ] | x = ξ = 0. (21)

Relation (21) and, as a consequence, the equality of the normal components of the total current were obtained by another method in 1948 in the work of G.A. Grinberg and V.A. Fock. It was shown in [1] - [6] that Equation (21) leads to the equality of the derivatives of the electric field strength along the normal to the surface:

[ E x x ] | x = ξ = 0.

Earlier, in [1] - [6] , we obtained and formulated all fifteen conditions at the interface of adjacent media necessary for solving the complete system of Equations (8)-(10):

1) The equality of normal components of the total current;

2) The equality of tangential projections of the vortex of the electric field;

3) The law of conservation of electric charge;

4) The equality of tangential components of the electric field and their derivatives in the tangential direction;

5) The equality of the derivatives of normal components of the total current in a direction tangential to the interface of adjacent media, with account of the effect of surface currents without introducing explicitly the surface charge s.

These conditions are valid in each cross section of the sample under study. For numerical modeling of the interaction of an electromagnetic field and a heterogeneous medium, it is most expedient to use, in our opinion, the through-counting schemes [1] - [6] , since the continuity condition of the total current must be satisfied at the interfaces of adjacent media. Indeed, in solving electro-dynamics problems, according to Equations (1)-(6), it is necessary to set six conditions on each face. Therefore, in 2D we must already have 24 conditions. Furthermore, on “sharp” faces, i.e., at angular points in Cartesian coordinates, it is necessary to additionally specify the Sommerfeld conditions in the asymptotic approximation of Malyuzhinets. Thus, to solve the problem even for one object, we must have 28 conditions. It is clear that for a group of bodies, it is technically quite difficult to realize the task.

Furthermore, Maxwell’s boundary conditions do not contain any closing relations for the induced surface charge and surface current. The use of known Leontovich-Shchukin conditions at interfaces, for example, for nanoobjects, is incorrect, since these conditions assume that the electrophysical properties and characteristics of the object vary little over distances of the order of a wavelength. The surface charge due to the passage of a direct current through the interface of adjacent media can be determined according to [1] [4] [5] [6] [7] :

σ = ( U / R S ) ε 0 [ ( ε 1 / λ 1 ) ( ε 2 / λ 2 ) ] . (22)

As follows from (22), the induced electric charge s is determined by the magnitude of the current, and also by a factor that takes into account the properties of the medium. In textbooks [8] - [14] , it is assumed that the surface charge is zero and is associated with electrostatics. As can be seen from Formula (22a), this assertion is valid only if

[ ( ε 1 / λ 1 ) ( ε 2 / λ 2 ) ] = 0.

The reformulation of Maxwell’s equations and their reduction to a wave equation make it possible to exclude from consideration surface charge and current [1] - [6] , but in this case the number of boundary conditions becomes significantly larger, since at each boundary it is required to specify nine first derivatives now with respect to coordinates, as well as E x , E y , and E z , and their derivatives with respect to time, i.e. 15 conditions, and the problem also becomes technically unrealizable. The way out is to use the continuity condition of the total current and through-counting schemes [1] - [6] . The calculation method consists in that the investigated region containing various inhomogeneous (angular, layered) structures is conventionally placed in the “casing” and we specify the boundary conditions only on it; usually, they are functions of the form E x = E x ( t ) sin ω t and E y = E y ( t ) cos ω . The medium inside the “casing” is considered together with inhomogeneous inclusions as a single medium, when its properties depend on the coordinates, i.e., ε ( r ) , λ ( r ) . This approach was used by us earlier to simulate the interaction of an electromagnetic wave and macro objects using the Comsol Multiphysics package [4] [6] . The characteristic dimension of the “casing” must be greater than the characteristic dimension of the object so that the reflected wave does not affect the boundary conditions. The proposed method also proved to be workable as applied to nanoobjects measuring 10 - 50 nm, wavelength ~400 - 500 nm (light). Note that the Comsol Multiphysics software package specifically focused on solution of electrodynamics problems. Unfortunately, using this package, as well as other software packages for solving electrodynamics problems, it is possible to obtain unreliable results, since the solution does not require matching conditions at angular points for the initial boundary-value problem. For clarity, we consider the one-dimensional wave equation of oscillations of a string:

2 u t 2 a 2 2 u x 2 = f ( x , t ) , ( x , t ) D ¯ , (23)

u | t = 0 = φ ( x ) , u t | t = 0 = ψ ( x ) , 0 x l , (24)

u | x = 0 = μ 1 ( t ) , u | x = l = μ 2 ( t ) , t 0 . (25)

With given functions f x x C ( D ¯ ) , φ C 2 ( 0 x l ) , ψ C 1 ( 0 x l ) , and μ 1 C 2 ( 0 x ) , you need to find a function u C 2 ( D ¯ ) , that satisfies Equation(23) in the region D ¯ , initial conditions (24), and boundary conditions of the first kind (25). Problem (23)-(25) describes the process of oscillations of a homogeneous string of length l stretched along a segment 0 x l . The first initial condition of (24) defines the graph u = φ ( x ) of the string at the initial time t = 0 , and the quantity ψ ( x ) from the second initial condition of (24) is the initial velocity of the string at the point with a coordinate x.

Note that when formulating the problem (23)-(25), certain restrictions must be imposed on the given functions φ , ψ , and μ 1 . In particular, matching conditions must be satisfied at the angular points of the region [15] [16] :

φ ( 0 ) = μ 1 ( 0 ) , φ ( l ) = μ 2 ( 0 ) , μ 1 ( 0 ) = ψ ( 0 ) , μ 2 ( 0 ) = ψ ( l ) . (26)

These conditions are necessary conditions for the continuous differentiability of the solution u ( x , t ) in a closed region D ¯ . Since the solution is u C 2 ( D ¯ ) , then in addition to conditions (26), second-order conditions must be satisfied:

μ 1 ( 0 ) a 2 φ ( 0 ) = f ( 0 , 0 ) , μ 2 ( 0 ) a 2 φ ( l ) = f ( l , 0 ) . (27)

Indeed, we differentiate conditions (25) twice with respect to t, and the first condition of (24), twice with respect to x, then

2 u t 2 | x = 0 , t = 0 = μ 1 ( 0 ) , 2 u t 2 | x = l , t = 0 = μ 2 ( 0 ) ,

2 u x 2 | x = 0 , t = 0 = φ ( 0 ) , 2 u x 2 | x = l , t = 0 = φ ( l ) . (28)

Substituting the values of the derivatives at the corresponding points to Equation (23), we obtain the required conditions (27) [15] [16] .

In the future, we will be interested in waves generated by sources that perform modulated oscillations:

A cos w t ± m A 2 cos w 1 t ± m A 2 cos w 2 t , (29)

where w 1 = w + Ω and w 2 = w + Ω . Expression (29) is easily transformed into

S = A [ 1 ± m cos ( Ω t K x ) ] cos ( w t K x ) (30)

or

S = A [ 1 ± m cos ( Ω t K x ) ] sin ( w t K x ) , (31)

where K = 0.5 ( k 1 k 2 ) , K = 0.5 ( k 1 + k 2 ) , and Ω w , Ω is the frequency and k 1 and k 2 are the wave numbers.

In fact, the wave generated by our modulated source does not differ from the superposition of the three waves that would be created by three independent sources. Let us imagine that waves of the form (30) and (31), which were formed with the help of special “tuning forks”, i.e. technical means, begin to propagate in a continuous medium. In this case, because of the interaction of a “direct” wave of the form (30) and (31) with the backward wave, the spectral line broadens. The broadening of the spectral line is explained by many effects: spin-spin interaction, signal absorption by the medium, exchange interaction, collision broadening, Doppler effect, etc. [17] [18] [19] . As is known, the general solution of the wave Equation (23) is the d’Alembert solution: the sum of the forward and backward waves, with the backward wave having a broadening of ± 2 Δ . As a result, we obtain for an electromagnetic wave

E x ( t ) = A [ 1 + m cos Ω t ] [ 2 sin w t sin ( w 2 Δ ) t sin ( w + 2 Δ ) t ] = 4 A [ 1 ± m cos Ω t ] sin 2 Δ t sin w t , (32)

where the frequency w corresponds to a wave moving to the right, and the frequency w ± 2 Δ corresponds to a wave moving to the left. A similar expression can also be obtained for E y :

E y ( t ) = A [ 1 ± m cos Ω t ] [ 2 cos w t cos ( w 2 Δ ) t cos ( w + 2 Δ ) t ] = 4 A [ 1 ± m cos Ω t ] sin 2 Δ t cos w t . (33)

In deriving (32) and (33), we used the well-known trigonometry formulas [20] :

sin α sin β sin γ = 1 4 [ sin ( α + β γ ) + sin ( β + γ α ) + sin ( γ + α β ) sin ( α + β + γ ) ] ,

sin α sin β cos γ = 1 4 [ cos ( α + β γ ) + cos ( β + γ α ) + cos ( γ + α β ) cos ( α + β + γ ) ]

with the condition α = β .

The broadening of the spectral lines is due to the interaction of the emitting atom with the surrounding particles: other atoms and molecules, ions and electrons. Therefore, functions of the form (32) and (33) continuously fill the frequency band w 2 Δ w w + 2 Δ . In addition, it takes some time to establish a signal (transient process), finally, for the components E x and E y we have

E x ( 0 , t ) = A [ w 2 Δ w + 2 Δ ( 1 2 e λ sin 2 φ t + e λ sin 2 φ t ) sin 2 φ t d φ ] × ( 1 m cos σ t ) sin w t 4 Δ , (34)

E y ( 0, t ) = A [ w 2 Δ w + 2 Δ ( 1 2 e λ sin 2 φ t + e λ sin 2 φ t ) sin 2 φ t d φ ] × ( 1 m c o s σ t ) c o s w t 4 Δ . (35)

It is easy to verify that E x ( 0, t ) and E z ( 0, t ) satisfy the necessary matching conditions, since E x ( 0 , t ) = 0 , E x ( 0 , t ) = 0 , E z ( 0 , t ) = 0 , and E z ( 0 , t ) = 0 , when the pulse propagates into a medium with zero initial conditions. A group of waves is a kind of oscillatory circuit with distributed parameters, in which “forced” oscillations are not established immediately, but only after a certain time after “including an external emf” [8] .

At the initial time moment, when t 0 , the expression in parentheses in the integrand expressions (34), (35) exactly corresponds to the “establishment function” in the oscillatory circuit under the action of an external sinusoidal electromotive force [8] .

Below we confine ourselves to linear electrodynamics. In this case the constitutive equations are linear relations. For constant fields, the constitutive equations corresponding to Equations (1), (2) can be written in the form

D i = ε i j E j , B i = μ i j H j .

Here ε i j and μ i j are the tensors of permeability and magnetic permeability. They are determined by the specific properties of the medium and thereby characterize its electromagnetic properties. Such constitutive equations hold only for rather slowly varying fields. For fast-changing fields that change rapidly compared to the characteristic relaxation time in the medium or to the periods of characteristic natural oscillations of the medium, the situation becomes somewhat more complicated. In this case, the state of the medium turns out to depend not only on the field at a given instant t, but also on its values at the preceding instants of time. This can be understood if we take into account, for example, the fact that the relaxation process that has begun in the medium under the action of a field that appeared at the time t and vanished after a certain time, which is much shorter than the relaxation time, will still flow after the field disappears. Therefore, for high-frequency fields, it is necessary to use constitutive equations of the form [21]

D i ( t ) = t d t ε ^ i j ( t t ) E j ( t ) , B i ( t ) = t d t μ ^ i j ( t t ) H j ( t ) . (36)

Relations (36) take into account the influence of prehistory on the electromagnetic properties of the medium. In this case, we usually speak of the time or frequency dispersion.

It is quite natural to raise the question of why, for example, the electric induction in relations (36) is considered to depend only on the electric field. In principle, one could speak of the dependence of D on the magnetic induction B . However, in the latter case, using the field equation

rot E = 1 c B t ,

it would be possible to eliminate the quantity B . In this case, although the electric induction would again be expressed only by the electric field strength, the spatial derivatives of E would arise in this constitutive equation. Therefore, we can say that the constitutive Equations (36) are valid only when the spatial derivatives can be neglected. In other words, such constitutive equations are valid for fields that vary slowly in space [21] .

When functions of the form (34) and (35) are used the relaxation nature of the process in a rapidly varying field is taken into account by us as a transition function, which is an integral not only in frequency but in time and space after numerical integration of the wave equation. This remark applies not only to temporal, but also to spatial dispersion. According to [21] , for fields that change sharply in space, it is obviously necessary to take into account the influence of the field at remote points on the electromagnetic properties of the medium at a given point in space. Indeed, for example, because of the transport processes, the state at a certain point in the medium will be assigned not only by the value of the field at this point, but also by the field in entire regions of the medium, from which the field effect is transferred as a result of the transfer of matter. Therefore, instead of the constitutive Equations (36), it is necessary to use spatially nonlocal relations that take into account not only temporal, but also spatial dispersion. For homogeneous isotropic and nonisotropic media, such relations can be written in the following form [21] :

D ( r , t ) = t d t d r ε ^ ( t t , | r r | ) E ( r , t ) , B ( r , t ) = t d t d r μ ^ ( t t , | r r | ) H ( r , t ) . (37)

Thus, the electromagnetic properties of such a medium are determined by two functions dependent on r and t. When spatial dispersion is taken into account, in accordance with (37), the vector functions D and B can be determined by solution of the integral-differential Equations (1), (2), and (37). In this case, only harmonic functions are considered in the paper without taking into account the modulation and broadening of the spectral line. According to the popular opinion, any uniformly continuous function can be reduced to a trigonometric series by means of an expansion in a Fourier series, i.e., to harmonic functions. However, the functions (34) and (35) containing the inverse hyperbolic cosine are not actually expanded in a Fourier series (message of the Wolfram Mathematics package). Note that for the Fourier expansion we are to calculate integrals of the form:

π π f ( x ) c o s ( n x ) d x , π π f ( x ) s i n ( n x ) d x ,

where f ( x ) is defined by Formulas (34) and (35); numerically they are calculated problematically, in view of the properties of the integrand, there is no analytical solution.

The most significant applications of the theory of electromagnetic field have been obtained in radio engineering, when the field functions are harmonicly dependent on time, i.e.,

E ( r , t ) = E 1 ( r ) c o s w t + E 2 ( r ) s i n w t ,

H ( r , t ) = H 1 ( r ) c o s w t + H 2 ( r ) s i n w t .

Passing to the complex domain and introducing complex quantities E ( r ) = E 1 ( r ) + i E 2 ( r ) and H ( r ) = H 1 ( r ) + i H 2 ( r ) and the complex dielectric constant of the medium, we exclude from Maxwell’s equations the time parameter t and determine the real electromagnetic field by means of the formulas

E = R e ( E e i w t ) , H = R e ( H e i w t ) ,

where R e ( ) is the real part of the complex quantity in brackets.

Note that for an arbitrary function E ( t ) , for example, (34) and (35), it is difficult to represent it in the form E ( t ) = a ( t ) cos φ ( t ) , since the amplitude a ( t ) and phase φ ( t ) cannot be uniquely determined and it is not clear how to decompose E ( t ) into the factors a and cos φ . Even greater difficulties arise in the case of transition to the complex representation W ( t ) = U ( t ) + i V ( t ) when the real oscillation E ( t ) is supplemented with the imaginary part V ( t ) . Problems arising in this case have been considered in detail in [22] . In this paper, it is emphasized that without defining uniquely amplitude, phase, and frequency, some methods using complex representation and claiming increased accuracy become meaningless.

We do not know of papers using (37) and in which solutions of the problems even for a homogeneous medium are given. However, it is not clear how to take into account temporal and spatial dispersion at the interfaces of adjacent media in the presence of strong discontinuities in the electromagnetic field. With this phenomenological approach in constitutive equations, it is necessary to take into account, for example, the dependence not only on frequency, but also the broadening of the spectral line Δ , as well as the parameter that determines the establishment of the vibration λ .

3. Numerical Simulation

The method of the “non-sharp” boundary is substantiated using the Dirichlet theorem for piecewise smooth functions and the condition of continuity of the total current. Numerical simulation of the propagation of a soliton-like pulse signal in a medium with strong electromagnetic field discontinuities is performed with account taken of the broadening of the spectral line.

3.1. Features of Calculating the Propagation of Electromagnetic Waves in Angular Structures. Analysis of a Soliton-Like Pulse Signal

Electromagnetic phenomena arising from the incidence of plane electromagnetic waves on the interface of different media play an important role in engineering, since all real devices are bounded by surfaces and are inhomogeneous in space. With the help of the proposed physicomathematical model, it is possible to investigate the transmission of an electromagnetic wave by a layered angular structure.

In connection with the fact that in every cross section of the layered medium, the condition of continuity of the total current is valid, we use the through-counting schemes without identifying explicitly the interface of adjacent media. The procedure of calculation of E x at the interface is as follows. Since E x 1 E x 2 , E x ( x ) experiences a discontinuity of the first-kind function (strong electromagnetic field discontinuity). Let us determine the value of the electric field strength at the point of discontinuity for x = ξ on condition that E x ( x ) is piecewise smooth and piecewise differentiable, i.e., it has finite one-sided derivatives E x ( x ) and E x ( x ) . At the points of discontinuity x i , we have

E x + ( x i ) = lim Δ x i + 0 E ( x i + Δ x i ) E ( x i + 0 ) Δ x i ,

E x ( x i ) = lim Δ x i 0 E ( x i + Δ x i ) E ( x i 0 ) Δ x i .

In this case, according to the Dirichlet theorem for a uniformly continuous function E ( x ) , at the discontinuity point ξ , we have:

E x = ξ = 1 2 [ E ( ξ 0 ) + E ( ξ + 0 ) ] . (38)

The Dirichlet condition also has a physical interpretation. When two solid conductors, dielectrics, and electrolytes come in contact in different combinations: metal-electrolyte, dielectric-electrolyte, metal-vacuum, etc., a double electric layer (DEL) is always formed at the interface of adjacent media, whose structure is usually unknown, but it significantly affects the electrokinetic phenomena, the rate of electrochemical processes, etc. It is important to note that in actual practice the electrophysical characteristics of E(x) in the DEL are continuously changing, so Equation (38) is valid for the case where the thickness of the DEL, i.e., the thickness of the interphase boundary, is much smaller than the characteristic dimension of the homogeneous medium. In the case of a composite, for example, a metal with impregnations of dielectric balls, at a sufficiently high concentration of both components and smallness in their characteristic dimensions, the interphase boundaries overlap and condition (38) may be violated.

If the thickness of the DEL is much smaller than the characteristic dimensions L of the objects under consideration, then Equation (38) also follows from the condition of linear variation in E ( x ) in the region of the DEL. In actual practice, the thickness of the DEL depends on the kind of contacting substances and can be tens of angstroms. According to modern concepts, the outer shell of the DEL consists of two parts: the first is formed by ions closely drawn to the surface of the metal (“dense” or “Helmholtz” layer of thickness h), and the second, by ions located at distances from the surface exceeding the radius of the ion, with the number of these ions decreasing with distance from the interface (“diffuse layer”). The distribution of the potential in the dense and diffuse parts of the DEL is in fact exponential [23] , i.e., the linearity condition E ( x ) is violated. At the same time, the sum of the charges of the dense and diffuse parts of the outer shell of the DEL is equal to the charge of the inner lining of the DEL of the metal surface. However, if the thickness of the DEL is much smaller than the characteristic dimension of the object, then the expansion of E ( x ) in a power series is valid, and we can confine ourselves to a linear approximation. According to the more general Dirichlet theorem (1829), the physical interpretation and knowledge of the E ( x ) dependences in the region of the DEL are not required in the justification of Equation (38). Nevertheless, the above-mentioned known physical features of the DEL confirm the validity of fulfillment of condition (38).

The concept of surface charge has only an auxiliary meaning and was introduced for convenience in considering the field of charges at distances sufficiently large in comparison with the dimensions of the charges themselves. By studying the field at the adjacent media interfaces, i.e., near or inside the charges, we must return to the concept of the volume distribution of charges when div E 0 , and the vector E is finite and continuous. In real physical processes, there is a double electrical layer at the interfaces of adjacent media and there is no disruption of the normal component of the electric field strength.

Thus, with account taken of the foregoing and of the validity of formula (38), in each cross section of the layered medium, for numerical solution it is expedient to use cross-counting schemes and at the same time to discretize the medium in such a way that the boundaries of the layers have common nodes.

For numerical simulation, we used the Wolfram Mathematics software package. We split the medium into finite elements in such a way that the nodes of the finite element mesh, lying on the interface, simultaneously belonged to media with different electrophysical properties. In this case, the condition of equality of the total currents and equality of the charge flows must be satisfied at the interface.

3.2. Analysis of Pulsed Soliton-Like Signals

We have considered waves of the following form, which will henceforth be called integrated

E x ( t ) = ( 1 m c o s θ t ) s i n ω t 4 Δ ω 2 Δ ω + 2 Δ ( 1 2 e λ s i n 2 φ t + e λ s i n 2 φ t ) s i n 2 φ t d φ , (39)

E y ( t ) = ( 1 m c o s θ t ) c o s ω t 4 Δ ω 2 Δ ω + 2 Δ ( 1 2 e λ s i n 2 φ t + e λ s i n 2 φ t ) s i n 2 φ t d φ , (40)

as well as waves that will be called ordinary:

E x ( t ) = ( 1 2 e λ sin 2 ω t + e λ sin 2 ω t ) sin 2 ω t ( 1 m cos θ t ) sin ω t , E y ( t ) = ( 1 2 e λ sin 2 ω t + e λ sin 2 ω t ) sin 2 ω t ( 1 m cos θ t ) cos ω t , (41)

where ω = 10 9 , θ = 10 6 , Δ = 10 4 , and λ = 0.1 . Let us consider the difference between the integrated and non-integrated waves, and also consider the integrated waves for different parameters m.

The graph of ordinary waves (component of the sine wave) with t [ 2 × 10 6 ,5 × 10 5 ] is shown in Figure 1(a). As can be seen from the graph, these waves satisfy the principle of the transient process and are of an impulse nature. It is also worth noting that near zero, they are very close to 0. Next, we consider the integrated waves on the same segment (component of the sine) (Figure 1(b)) and also consider in detail the shape of the wave crests (component of the sine), reducing the segment to t [ 2 × 10 6 ,5 × 10 5 ] .

Let us represent the graphs of the same waves for the same parameters, but for the component of the cosine (Figure 2). Obviously, ordinary waves have a sinusoidal shape upon close examination. As could be seen from the graphs above, the integrated waves change their shape. Waves become “impulsive”, which is a significant characteristic of waves, for example, for radar problems. In this case necessary conditions for coordination are fulfilled. The penetrating power of such a pulse in complex media, for example, the ionosphere and an

Figure 1. Graph of ordinary waves (component of the sine) (a) and graphs of integrated waves (components of the sine), m = 1, different time scales (b and c).

ust be much higher than for a rectangular and saw-tooth signal or other form of signal when necessary conditions for coordination are not fulfilled and as a consequence, we have the expenditure of energy on “noise” due to discontinuities in the solution of differential equations at the interface of adjacent media.

Let us consider the energy of ordinary waves and, by reducing the scale, the energy of the integrated waves is more detail (Figure 3). It can be seen from the figure that the qualitative characteristics of waves after integration are modified. The rate of increase in the electromagnetic pulse of the integrated waves also becomes several orders of magnitude higher. This means that the penetrating power of the integrated waves will be much higher than that of ordinary waves.

The exchange of information between cells in living organisms (electrolytes) occurs with the help of soliton-like pulses [24] , therefore, it can be expected that these pulses may prove to be promising for signal transmission in the ionosphere and plasma.

We consider the integrated waves for different values of the parameter m (Figure 4). The figure shows only the components of the sine for m = 0.7, 1.0, and 1.3. As noted earlier, at m = 1.0, waves after a certain period are damped (go to zero). As the parameter m decreases, the wave in the damping regions is smoothed out and no longer has a total attenuation, but only decreases to a certain limit. When the parameter m increases at the damping points, a new wave crest appears which breaks the damping region into two new ones, smaller

Figure 2. Graph of ordinary waves (component of the cosine) (a) and graphs of integrated waves (components of the cosine), m = 1, different time scales (b and c).

Figure 3. Graph of the energy of ordinary waves (a) and graphs of integrated waves for different time scales (b and c).

in size. It is also worth noting that the drawn conclusions are valid for the components of the cosine and for the energy and all components for ordinary waves.

3.3. Numerical Simulation of an Electromagnetic Field in a Heterogeneous Medium

The obtained waves can be used to solve the problem of their propagation in a medium with strong field discontinuities. Imagine that we have a square domain with dimensions [ 0 ; 1 ] × [ 0 ; 1 ] . At the center of the domain, there is a copper bar (parallelepiped) (Figure 5).

It is necessary to determine how the intensity of the electromagnetic field varies in the given region. Earlier, we obtained a mathematical model that describes this process using differential equations, namely:

λ E t + ε ε 0 2 E t 2 = 1 μ μ 0 ( 2 E grad ( div E ) ) , (42)

E | t = 0 = 0 , (43)

E t | t = 0 = 0 , (44)

Figure 4. Graphs of integrated waves (components of the sine) for different values of the parameter m: (a) m = 0.7; (b) 1.0; and (c) 1.3.

E | Γ = φ ( t ) . (45)

Equations (43) and (44) are the initial conditions for the problem, (45) is the boundary condition (39, 40), and Γ is the boundary of the domain. The function E depends on three variables and is of the form E ( t , x , y ) . It is also a vector E = ( E x , E y ) . The boundary Γ of our domain is nothing but the boundary of a square. To solve this problem, we used the Wolfram Mathematics package.

Figure 5. Problem solution domains: 1) copper bar; 2) domain with air.

The results obtained at m = 1 for the multiplier ( 1 cos θ t ) and various instants of time are shown in Figure 6.

In conclusion, we note that the attempt to compare the results of numerical simulation of this problem in the Wolfram Mathematics and COMSOL Multiphysics environment was unsuccessful, because the task using the COMSOL Multiphysics package was not considered when specifying functions of the form (39), (40) on the boundaries. The calculation results are also difficult to compare, for example, with calculations based on difference schemes for the equations of oscillations of a string [25] , since the problem is considered on a time interval 0 < t T , i.e., the matching problem for t = 0 is not considered. Attempts to calculate using the existing wave packets for a group of waves also failed, since the wave packet [26]

Ψ ( x , t ) = C ( x , t ) e i ( ω 0 t k 0 x ) = 2 C ( k 0 ) sin ( [ ( d ω d k ) 0 t x ] Δ k ) ( d ω d k ) 0 t x e i ( ω 0 t k 0 x )

does not satisfy the matching conditions at the initial time. Therefore, the Wolfram Mathematics package informs about the error in formulating the initial boundary-value problem. We note, however, that the numerical calculations in both packages qualitatively (increased concentration of the electromagnetic field on sharp edges) correspond to the experimental data provided that COMSOL Multiphysics uses harmonic functions at the boundary of the region without taking into account the time of signal establishment, but without transient processes.

4. Conclusions and Suggestions

A coordinated physicomathematical model for the propagation of a soliton-like electromagnetic pulse in a heterogeneous medium is developed in the presence

Figure 6. Distributions of the intensity in the region m = 1 for the multiplier (1-cosθt) and different instants of time: top, t = 4 × 109; center, t = 6 × 109; bottom, t = 8 × 109.

of strong discontinuities in the electromagnetic field. The model is based on the reduction of Maxwell’s equations to the well-studied wave equation. When the electromagnetic pulse was specified, its amplitude modulation was taken into account, as was the nonstationary broadening of the spectral line. The results obtained: the condition of continuity of the total current in a heterogeneous medium, a consistent soliton-like pulse of the electromagnetic field, allow us to justify the use of through-counting schemes for numerical analysis of the interaction of an electromagnetic field and a piecewise homogeneous body.

The method of the “non-sharp” boundary is substantiated using the Dirichlet theorem for piecewise smooth functions and the condition of continuity of the total current. Numerical simulation of the propagation of a soliton-like pulse signal in a medium with strong electromagnetic field discontinuities is performed with account taken of the broadening of the spectral line.

The work was supported by the Belarusian Republican Foundation for Basic Research, project No. F16K-075.

Cite this paper

Grinchik, N.N., Andrejev, V.G., Zayats, G.M. and Grinchik, Yu.N. (2018) Non-Monochromatic Electromagnetic Radiation of Inhomogeneous Media. Journal of Electromagnetic Analysis and Applications, 10, 13-33. https://doi.org/10.4236/jemaa.2018.102002

References

  1. 1. Grinchik, N.N. (2008) Simulation of Electrophysical and Thermophysical Processes in Porous Media. Nauka, Minsk, Belarus. (In Russian)

  2. 2. Grinchik, N.N. (2000) Interaction of Thermal and Electric Phenomena in Polariged Media, Matemat. Modelirovanie, 12, 67-76.

  3. 3. Grinchik, N.N. and Tsurko, V.A. (2002) Problem of Modeling of the Interaction of Nonstationary Electric, Thermal, and Diffusion Feilds in Layered Media. Journal of Engineering Physics and Thermophysics, 75, 693-699.

  4. 4. Grinchik, N.N., et al. (2011) Fundament Problems of the Electrodynamics of Heterogeneous Media with Boundary Conditions Corresponding to the Total-current Continuity. Electromagnetic Wave Propagation in Complex Matter, INTECH, Rijeka.

  5. 5. Grinchik, N.N. and Grinchik, Yu.N. (2012) Fundamenttal Problems of the Electrodynamics of Heterogeneous Media. Physics Research International, 2012, Article ID 185647.

  6. 6. Grinchik, N.N. (2014) Electrodynamics of Inhomogeneous (Laminated, Angular) Structures. Journal of Electromagnetic Analysis and Applications, 6, 57-105. https://doi.org/10.4236/jemaa.2014.65009

  7. 7. Stratton, J.A. (1948) The Theory of Electromagnetism. OGIZ, Moscow-Leningrad. (In Russian)

  8. 8. Kalashnikov, S.G. (1985) Electricity. Nauka, Moscow. (In Russian)

  9. 9. Tamm, I.E. (2003) Foundations of the Theory of Electricity. Fizmatgiz, Moscow. (In Russian)

  10. 10. Purcell, E. (1971) Electricity and Magnetism. Nauka, Moscow. (In Russian)

  11. 11. Irodov, I.E. (1983) Fundamental Lows of Electromagnetism. Vysshaya Shkola, Moscow. (In Russian)

  12. 12. Zil’berman, G.E. (1970) Electricity and Magnetism. Nauka, Moscow. (In Russian)

  13. 13. Telesnin, R.V. and Iakovlev V.F. (1969) A Course in Physics. Electricity. Prosveshchenie, Moscow. (In Russian)

  14. 14. Sivukhin, D.V. (1996) Electricity, in Two Parts. Nauka, Moscow. (In Russian)

  15. 15. Erofeenko, V.T. and Kozlovskaya, I.C. (2001) Partial Equations with Applications in Economics: A Course of Lectures. BGU, Minsk, Belarus. (In Russian)

  16. 16. Sveshnikov, A.G., Bogolubov, A.N. and Kravtsov, V.V. (2004) Lectures on Mathematical Physics: A Textbook for University Students. Nauka, Moscow. (In Russian)

  17. 17. Sobel’man, I.I. (1954) On the Theory of Width of Atomic Spectral Lines. Physics-Uspekhi, 54, 551-586. https://doi.org/10.3367/UFNr.0054.195412c.0551

  18. 18. Vol’kenstein, M.V. (1951) Molecular Optics. Gostekhizdat, Moscow-Leningrad. (In Russian)

  19. 19. DeWitt, C. and Blandin, A., Eds. (1966) Quantum Optics and Electronics: Lectures at Summer School of Theoretical Physics of the University of Grenoble, Les Houches, France. Mir, Moscow. (In Russian)

  20. 20. Bronshtein, I.N. (2010) Reference Book of Mathematics for Engineers and University Students. Publishing House Lan, St. Petersburg. (In Russian)

  21. 21. Silin, V.P. and Rukhadze A.A. (2013) Electromagnetic Properties of Plasma and Plasmalike Media. 3rd Edition, Publishing House Librokom, Moscow. (In Russian)

  22. 22. Vakman, D.E. and Vainshtein L.A. (1977) Amplitude, Phase, Frequency—Fundamental Concepts of Oscillation Theory. Physics-Uspekhi, 20, 1002-1016. https://doi.org/10.3367/UFNr.0123.197712f.0657

  23. 23. Evans, D.F. and Matesich, M.A. (1977) Measurement and Interpretation of Electrical Conductivity. In: Yeager, E. and Salkind, A.J., Eds., Measurement Methods in Electrochemistry, Vol. 2, Mir, Moscow. (In Russian)

  24. 24. Antonov, V.F., Chernysh, A.M., Pasechnik, V.I., et al. (2000) Biophysics. Gumanitarnyi Izdatel’skii Tsentr “Vlados”, Moscow. (In Russian)

  25. 25. Samarskii, A.A. (1989) Theory of Difference Schemes. Nauka, Moscow. (In Russian)

  26. 26. Blokhintsev, D.I. (2015) Quantum Mechanics: Lectures on Selected Issues. 2nd Edition, Lenand, Moscow. (In Russian)