^{1}

^{*}

^{2}

^{3}

^{4}

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.

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.

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 [

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

Earlier, in [

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 [

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 [

σ = ( 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 [

[ ( ε 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 [

∂ 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 [

φ ( 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) [

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. [

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 [

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” [

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 [

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 [

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 [

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 [

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 [

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 λ .

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.

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 [

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.

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

Let us represent the graphs of the same waves for the same parameters, but for the component of the cosine (

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 (

The exchange of information between cells in living organisms (electrolytes) occurs with the help of soliton-like pulses [

We consider the integrated waves for different values of the parameter m (

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.

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) (

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)

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.

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

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 [

Ψ ( 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.

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. 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.

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