# Fundamental Problems of the Electrodynamics of Heterogeneous Media

Authors: N. N. Grinchik; Yu. N. Grinchik
Journal: Physics Research International
 (2012)

---

## Abstract

The consistent physic-mathematical model of propagation of an electromagnetic wave in a heterogeneous medium is constructed using the generalized wave equation and the Dirichlet theorem. Twelve conditions at the interfaces of adjacent media are obtained and justified without using a surface charge and surface current in explicit form. The conditions are fulfilled automatically in each section of counting schemes for calculations. A consistent physicomathematical model of interaction of nonstationaly electric and thermal fields in a layered medium with allowance or mass transfer is constructed. The model is based on the methods of thermodynamics and on the equations of an electromagnetic field and is formulated without explicit separation of the charge carriers and the charge of an electric double layer.


---

## Body

## 1. Introduction

Let us consider the interfaceS between two media having different electrophysical properties. On each of its sides the magnetic-field and magnetic-inductance vectors as well as the electric-field and electric-displacement vectors are finite and continuous; however, at the surface S they can experience a discontinuity of the first kind. Furthermore, at the interface there arise induced surface charges σ and surface currents i (whose vectors lie in the plane tangential to the surface S) under the action of an external electric field.The existence of a surface charge at the interfaceS between the two media having different electrophysical properties is clearly demonstrated by the following example. We will consider the traverse of a direct current through a flat capacitor filled with two dielectric materials having relative permittivities ε1 and ε2 and electrical conductivities λ1 and λ2. A direct-current voltage U is applied to the capacitor plates; the total resistance of the capacitor is R (Figure 1). It is necessary to calculate the surface electric charge induced by the electric current.Figure 1
Dielectric media inside a flat capacitor.From the electric-charge conservation law follows the constancy of flow in a circuit; therefore, the following equation is fulfilled:(1)λ1En1=λ2En2=U(RS),
where En1 and En2 are the normal components of the electric-field vector.At the interface between the dielectrics, the normal components of the electric-inductance vector change spasmodically under the action of the electric field by a value equal to the value of the induced surface chargeσ:(2)ε0ε1En1-ε0ε2En2=σ.
Solving the system of (1) and (2), we obtain the expression for σ:(3)σ=(URS)ε0[(ε1λ1)-(ε2λ2)].
It follows from (3) that the charge σ is determined by the current and the multiplier accounting for the properties of the medium. If(4)(ε1λ1)-(ε2λ2)=0,
a surface charge σ is not formed. What is more, recent trends are toward increased use of micromachines and engines made from plastic materials, where the appearance of a surface charge is undesirable. For oiling of elements of such machines, it is best to use an oil with a permittivity εoil satisfying the relation(5)ε1<εoil<ε2.
This oil makes it possible to decrease the electrization of the moving machine parts made from dielectric materials. In addition to the charge σ, a contact potential difference arises always independently of the current.An electric field interacting with a material is investigated with the use of the Maxwell equation (1857). For convenience and correctness of further presentation we denote the axial vector of magnetic field strength, induction asH⇀ and B⇀(6)jtotal=∇×H⇀,∇⋅D=ρ,(7)-∂B⇀∂t=∇×E,∇⋅B⇀=0,
where  jtotal=λE+∂D⇀/∂t,B⇀=μμ0H⇀. In this case, at the interface S the above system of equations is supplemented with the boundary conditions [1, 2]:
(8)Dn1-Dn2=σ,(9)Eτ1-Eτ2=0,(10)Bn1-Bn2=0,(11)H⇀τ1-H⇀τ2=[ι⃗τ⋅n⃗].
The indices (subscripts) n and τ denote the normal and tangential components of the vectors to the surface S, and the indices 1 and 2 denote the adjacent media with different electrophysical properties. The index τ denotes any direction tangential to the discontinuity surface. At the same time, a closing relation is absent for the induced surface charge σ, which generates a need for the introduction of an impedance matrix [3–7] that is determined experimentally or, in some cases, theoretically from quantum representations [8–12]. We note that the Maxwell equation of electromagnetic field does not coincide completely with the modern form of them (6), (7) which was first given by Hertz (1884) and independently of him by Heaviside (1885). The equations of electromagnetic field with precise division of polar and axial vectors are of importance in considering the boundary conditions. One should bear in mind that electric field strength E is a polar vector. Gradients of concentration and temperature are also polar vectors. On the contrary, magnetic field strength H⇀, as well as H⇀, are axial vectors, or pseudovectors. This is well seen from the formula for the Lorentz force F=q[VH⇀]. Indeed, the behavior of H⇀ on reflection at the origin of coordinates and substitution r→(-r) is determined by the behavior of the polar vectors F and V and the properties of their vector product. On substitution r→(-r) the directions of the vectors F and V change to the opposite; the sign of the vector product also reverses. Consequently, on substitution r→(-r) the vector H⇀ must remain constant. This property is assumed to be an attribute of the axial vector.The induced surface chargeσ not only characterizes the properties of a surface, but also represents a function of the process, that is, σ(E(∂E/∂t,H⇀(∂H⇀/∂t))); therefore, the surface impedances [3–7] are true for the conditions under which they are determined. These impedances cannot be used in experiments conducted under other experimental conditions.The problem of determination of surface charge and surface current on metal-electrolyte boundaries becomes even more complicated in investigating and modeling nonstationary electrochemical processes, for example, pulse electrolysis, when lumped parametersL,C, and R cannot be used in principle.We will show thatσ can be calculated using the Maxwell phenomenological macroscopic electromagnetic equations and the electric-charge conservation law accounting for the special features of the interface between the adjacent media.Separate consideration will be given to ion conductors. In constructing a physicomathematical model, we take into account thatE and H⇀ are not independent functions; therefore, the wave equation for E or H⇀ is more preferable than the system of (6) and (7). It will be shown that in consideration of propagation, interference of waves in the heterogeneous medium it is preferable to use the wave equation for a polar vector of electric field strength.


## 2. Electron Conductors: New Closing Relations on the Boundaries of Adjacent Media

### 2.1. Generalized Wave Equation for E and Conditions on the Boundaries in the Presence of Strong Discontinuities of the Electromagnetic Field

#### 2.1.1. Physicomathematical Model

We will formulate a physicomathematical model of propagation of an electromagnetic field in a heterogeneous medium. Let us multiply the left and right sides of the equation for the total current (see (6)) by μμ0 and differentiate it with respect to time. Acting by the operator rot on the left and right sides of the first equation of (7) on condition that μ=const we obtain(12)∂jtotal∂t=1μμ0∇2E-1μμ0grad(div⁡E).
In Cartesian coordinates, (12) will take the form(13)∂jtotalx∂t=1μμ0(∂2Ex∂x2+∂2Ex∂y2+∂2Ex∂z2)-1μμ0∂∂x(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotaly∂t=1μμ0(∂2Ey∂x2+∂2Ey∂y2+∂2Ey∂z2)-1μμ0∂∂y(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotalz∂t=1μμ0(∂2Ez∂x2+∂2Ez∂y2+∂2Ez∂z2)-1μμ0∂∂z(∂Ex∂x+∂Ey∂y+∂Ez∂z).
At the interface, the following relation [2] is also true: (14)div⁡iτ+Iqx1-Iqx2=-∂σ∂t.
Let us write conditions (see (8)–(11)) in the Cartesian coordinate system: (15)Dx1-Dx2=σ,(16)Ey1-Ey2=0,(17)Ez1-Ez2=0,(18)Bx1-Bx2=0,(19)Hy1-Hy2=iz,(20)Hz1-Hz2=iy,
where iτ=iyj+izk is the surface-current density, and the coordinate x is directed along the normal to the interface. The densities iy and iz of the surface currents represent the electric charge carried in unit time by a segment of unit length positioned on the surface drawing the current perpendicularly to its direction.The order of the system of differential equations (13) is equal to (16). Therefore, at the interface S, it is necessary to set, by and large, nine boundary conditions. Moreover, the three additional conditions (15), (19), and (20) containing (prior to the solution) unknown quantities should be fulfilled at this interface. Consequently, the total number of conjugation conditions at the boundary S should be equal to (12) for a correct solution of the problem.As an example we consider a one-dimensional case of heat propagation in the layered material. As is known, to solve the problem we must specify two conditions at the interface: the equality of heat fluxes and the equality of temperatures, since in the heat conduction equation we have the second derivative with respect to the coordinate.Differentiating expression (15) with respect to time and using relation (14), we obtain the following condition for the normal components of the total current at the medium-medium interface:(21)div⁡iτ+jtotalx1=jtotalx2,
that allows one to disregard the surface charge σ. Let us introduce the arbitrary function f:[f]|x=ξ=f1|x=ξ+0-f2|x=ξ-0. In this case, expression (21) will take the form(22)[div⁡iτ+jtotalx1]|x=ξ=0.
It is assumed that, at the medium-medium interface, Ex is a continuous function of y and z. Then, differentiating Equation (21) with respect to y and z, we obtain(23)[∂∂yjtotalx1]|x=ξ=-∂(div⁡iτ)∂y,[∂∂zjtotalx1]|x=ξ=-∂(div⁡iτ)∂z.
Let us differentiate conditions (18)–(20) for the magnetic induction and the magnetic-field strength with respect to time. On condition that B⇀=μμ0H⇀:(24)[∂Bx∂t]|x=ξ=0,[1μμ0∂By∂t]|x=ξ=∂iz∂t,[1μμ0∂Bz∂t]|x=ξ=∂iy∂t.
Using (7) and expressing (24) in terms of projections of the electric-field rotor, we obtain(25)[rot⁡xE]|x=ξ=0,[∂Ez∂y-∂Ey∂z]|x=ξ=0,(26)[1μμ0rot⁡yE]|x=ξ=∂iz∂t,or[1μμ0(∂Ex∂z-∂Ez∂x)]|x=ξ=∂iz∂t,(27)[1μμ0rot⁡zE]|x=ξ=∂iz∂t,or[1μμ0(∂Ey∂x-∂Ex∂y)]|x=ξ=∂iz∂t.
Here, (25) is the normal projection of the electric-field rotor, (26) is the tangential projection of the rotor on y, and (27) is the rotor projection on z.Assuming thatEy and Ez are continuous differentiable functions of the coordinates y and z, from conditions (16) and (17) we find(28)[∂Ey∂y]|x=ξ=0,[∂Ey∂z]|x=ξ=0,[∂Ez∂y]|x=ξ=0,[∂Ez∂z]|x=ξ=0.
In accordance with the condition that the tangential projections of the electric field on z and y are equal and in accordance with conditions (16) and (17), the expressions for the densities of the surface currents iz and iy take the form(29)iz=λ̅Ez|x=ξ,iy=λ̅Ey|x=ξ,
where
(30)λ̅=12(λ1+λ2)|x=ξ
is the average value of the electrical conductivity at the interface between the adjacent media in accordance with the Dirichlet theorem for a piecewise-smooth, piecewise-differentiable function.Consequently, formulas (28)–(30) yield(31)[div⁡iτ⃗]|x=ξ=0.
Relations (31) and (32) and hence the equality of the normal components of the total current were obtained (in a different manner) by Grinberg and Fok [13]. In this work, it has been shown that condition (31) leads to the equality of the derivatives of the electric field strength along the normal to the surface(32)[∂Ex∂x]|x=ξ=0.
With allowance for the foregoing we have twelve conditions at the interface between the adjacent media that are necessary for solving the complete system of (13):(a)
the functionsEy and Ez are determined from (16) and (17);(b)
Ex is determined from condition (22);(c)
the values of∂Ex/∂y, ∂Ex/∂z, and ∂Ex/∂x are determined from relations (23) with the use of the condition of continuity of the total-current normal component at the interface (22) and the continuity of the derivative of the total current with respect to the coordinate x;(d)
the values of∂Ey/∂y, ∂Ey/∂z, and ∂Ez/∂z are determined from conditions (28) and (29) in consequence of the continuity of the tangential components of the electric field along y and z;(e)
the derivatives∂Ey/∂x and ∂Ez/∂x are determined from conditions (26) and (27) as a consequence of the equality of the tangential components of the electric-field rotor along y and z.Note that condition (21) was used by us in [14] in the numerical simulation of the pulsed electrochemical processes in the one-dimensional case. Condition (25) for the normal component of the electric-field rotor represents a linear combination of conditions (28) and (29); therefore, rot⁡xE=0 and there is no need to use it in the subsequent discussion. The specificity of the expression for the general law of electric-charge conservation at the interface is that the components ∂Ey/∂y and ∂Ez/∂z are determined from conditions (28) and (29) that follow from the equality and continuity of the tangential components Ey and Ez at the boundary S.Thus, at the interface between the adjacent media the following conditions are fulfilled: the equality of the total-current normal components; the equality of the tangential projections of the electric-field rotor; the electric-charge conservation law; the equality of the electric-field tangential components and their derivatives in the tangential direction; the equality of the derivatives of the total-current normal components in the direction tangential to the interface between the adjacent media, determined with account for the surface currents and without explicit introduction of a surface charge. They are true at each cross section of the sample being investigated.


#### 2.1.2. Features of Calculation of the Propagation of Electromagnetic Waves in Layered Media

The electromagnetic effects arising at the interface between different media under the action of plane electromagnetic waves have a profound impact on the equipment because all real devices are bounded by the surfaces and are inhomogeneous in space. At the same time, the study of the propagation of waves in layered conducting media and, according to [15], in thin films is reduced to the calculation of the reflection and transmission coefficients; the function E(x) is not determined in the thickness of a film, that is, the geometrical-optics approximation is used.The physicomathematical model proposed allows one to investigate the propagation of an electromagnetic wave in a layered medium without recourse to the assumptions used in [3–5, 7].Since conditions (21)–(29) are true at each cross section of a layered medium, we will use schemes of through counting without an explicit definition of the interface between the media. In this case, it is proposed to calculate Ex at the interface in the following way.In accordance with (15), Ex1≠Ex2, that is, Ex(x) experiences a discontinuity of the first kind. Let us determine the strength of the electric field at the discontinuity point x=ξ on condition that Ex(x) is a piecewise-smooth, piecewise-differentiable function having finite one-sided derivatives Ex+′(x) and Ex-′(x). At the discontinuity points xi,(33)Ex+′(xi)=lim⁡Δxi→+0E(xi+Δxi)-E(xi+0)Δxi,Ex-′(xi)=lim⁡Δxi→-0E(xi+Δxi)-E(xi-0)Δxi.
In this case, in accordance with the Dirichlet theorem [16], the Fourier series of the function E(x) at each point x, including the discontinuity point ξ, converges and its sum is equal to(34)Ex=ξ=12[E(ξ-0)+E(ξ+0)].
The Dirichlet condition (34) also has a physical meaning. In the case of contact of two solid conductors, for example, dielectrics or electrolytes in different combinations (metal-electrolyte, dielectric-electrolyte, metal-vacuum, and so on), at the interface between the adjacent media there always arises an electric double-layer (EDL) with an unknown (as a rule) structure that, however, substantially influences the electrokinetic effects, the rate of the electrochemical processes, and so on. It is significant that, in reality, the electrophysical characteristics λ, ε, and E(x) change uninterruptedly in the electric double-layer; therefore, (34) is true for the case where the thickness of the electric double-layer, that is, the thickness of the interphase boundary, is much smaller than the characteristic size of a homogeneous medium. In a composite, for example, in a metal with embedments of dielectric balls, where the concentration of both components is fairly large and their characteristic sizes are small, the interphase boundaries can overlap and condition (34) can break down.If the thickness of the electric double-layer is much smaller than the characteristic sizeL of an object, (34) also follows from the condition that E(x) changes linearly in the EDL region. In reality, the thickness of the electric double-layer depends on the kind of contacting materials and can comprise several tens of angstroms [17]. In accordance with the modern views, the outer coat of the electric double-layer consists of two parts, the first of which is formed by the ions immediately attracted to the surface of the metal (a “dense” or a “Helmholtz” layer of thickness h), and the second is formed by the ions separated by distances larger than the ion radius from the surface of the layer, and the number of these ions decreases as the distance between them and the interface (the “diffusion layer”) increases. The distribution of the potential in the dense and diffusion parts of the electric double-layer is exponential in actual practice [17], that is, the condition that E(x) changes linearly breaks down; in this case, the sum of the charges of the dense and diffusion parts of the outer coat of the electric double-layer is equal to the charge of its inner coat (the metal surface). However, if the thickness of the electric double-layer h is much smaller than the characteristic size of an object, the expansion of E(x) into a power series is valid and one can restrict oneself to the consideration of a linear approximation. In accordance with the more general Dirichlet theorem (1829), a knowledge of this function in the EDL region is not necessary to substantiate (34). Nonetheless, the above-indicated physical features of the electric double-layer lend support to the validity of condition (34).The condition at interfaces, analogous to (34), has been obtained earlier [18] for the potential field (where rot⁡E=0) on the basis of introduction of the surface potential, the use of the Green formula, and the consideration of the discontinuity of the potential of the double-layer. In [18], it is also noted that the consideration of the thickness of the double-layer and the change in its potential at h/L≪1 makes no sense in general; therefore, it is advantageous to consider, instead of the volume potential, the surface potential of any density. Condition (34) can be obtained, as was shown in [16], from the more general Dirichlet theorem for a nonpotential vorticity field [18].Thus, the foregoing and the validity of conditions (15)–(17) and (23)–(29) at each cross-section of a layered medium show that, for numerical solution of the problem being considered it is advantageous to use schemes of through counting and make the discretization of the medium in such a way that the boundaries of the layers have common points.The medium was divided into finite elements so that the nodes of a finite-element grid, lying on the separation surface between the media with different electrophysical properties, were shared by these media at a time. In this case, the total currents or the current flows at the interface should be equal if the Dirichlet condition (34) is fulfilled.


#### 2.1.3. Results of Numerical Simulation of the Propagation of Electromagnetic Waves in Layered Media

Let us analyze the propagation of an electromagnetic wave through a layered medium that consists of several layers with different electrophysical properties in the case where an electromagnetic-radiation source is positioned on the upper plane of the medium. It is assumed that the normal component of the electric-field vectorEx=0 and its tangential component Ey=asin⁡(ωt), where a is the electromagnetic-wave amplitude (Figure 2).Figure 2
Scheme of a layered medium: layers 1, 3, and 5 are characterized by the electrophysical parametersε1, λ1,and μ1, and layers 2, 4, and 6 by ε2, λ2, and μ2.In this example, for the purpose of correct specification of the conditions at the lower boundary of the medium, an additional layer is introduced downstream of layer 6; this layer has a larger conductivity and, therefore, the electromagnetic wave is damped out rapidly in it. In this case, the conditionEy=Ez=0 can be set at the lower boundary of the medium. The above manipulations were made to limit the size of the medium being considered because, in the general case, the electromagnetic wave is attenuated completely at an infinite distance from the electromagnetic-radiation source.Numerical calculations of the propagation of an electromagnetic wave in the layered medium with electrophysical parametersε1=ε2=1, λ1=100, λ2=1000, and μ1=μ2=1 were carried out. Two values of the cyclic frequency ω=2π/T were used: in the first case, the electromagnetic-wave frequency was assumed to be equal to ω=1014 Hz (infrared radiation), and, in the second case, the cyclic frequency was taken to be ω=109 Hz (radiofrequency radiation).As a result of the numerical solution of the system of (13) with the use of conditions S (22)–((31) and (32) at the interfaces, we obtained the time dependences of the electric-field strength at different distances from the surface of the layered medium (Figure 3).Figure 3
Time change in the tangential component of the electric-field strength at a distance of 1μm (1), 5 μm (2), and 10 μm (3) from the surface of the medium at λ1=100, λ2=1000, ε1=ε2=1, μ1=μ2=1, and ω=1014 Hz. t, sec.The results of our simulation (Figure4) have shown that a high-frequency electromagnetic wave propagating in a layered medium is damped out rapidly, whereas a low-frequency electromagnetic wave penetrates into such a medium to a greater depth. The model developed was also used for calculating the propagation of a modulated signal of frequency 20 kHz in a layered medium. As a result of our simulation (Figure 5), we obtained changes in the electric-field strength at different depths of the layered medium, which points to the fact that the model proposed can be used to advantage for calculating the propagation of polyharmonic waves in layered media; such a calculation cannot be performed on the basis of the Helmholtz equation.Figure 4
Distribution of the amplitude of the electric-field-strength at the cross section of the layered medium:ω=1014 (1) and 109 Hz (2). y,μm.Figure 5
Time change in the electric-field strength at a distance of 1 (1), 5 (2), and 10μm (3) from the surface of the medium. t, sec.The physicomathematical model developed can also be used to advantage for simulation of the propagation of electromagnetic waves in media with complex geometric parameters and large discontinuities of the electromagnetic field (Figure6).Distribution of the amplitude of the electric-field strength in the two-dimensional medium and in depth atε1=15, ε2=20, λ1=10-6, λ2=10, μ1=μ2=1, and ω=109 Hz (the dark background denotes medium 1, and the light background denotes medium 2). x,y, mm; E, V/m.
(a)
Distribution of the amplitude of the electric-field strength in the two-dimensional medium(b)
Distribution of the amplitude of the electric-field strength in depthFigure6(a) shows the cross-sectional view of a cellular structure representing a set of parallelepipeds with different cross-sections in the form of squares. The parameters of the materials in the large parallelepiped are denoted by Index 1, and the parameters of the materials in the small parallelepipeds (the squares in the figure) are denoted by Index 2.An electromagnetic wave propagates in the parallelepipeds (channels) in the transverse direction. It is seen from Figure6(b) that, in the cellular structure there are “silence regions,” where the amplitude of the electromagnetic-wave strength is close to zero, as well as inner regions where the signal has a marked value downstream of the “silence” zone formed as a result of the interference.


#### 2.1.4. Results of Numerical Simulation of the Scattering of Electromagnetic Waves in Angular Structures

It is radiolocation and radio-communication problems that are among the main challenges in the set of problems solved using radio-engineering devices.Knowledge of the space-time characteristics of diffraction fields of electromagnetic waves scattered by an object of location into the environment is necessary for solving successfully any radiolocation problem. Irradiated objects have a very intricate architecture and geometric shape of the surface consisting of smooth portions and numerous wedge-shaped for formations of different type-angular joints of smooth portions, surface fractures, sharp edges, and so forth, with rounded radii much smaller than the probing-signal wavelength. Therefore, solution of radiolocation problems requires the methods of calculation of the diffraction fields of electromagnetic waves excited and scattered by different surface portions of the objects, in particular, by wedge-shaped formations, be known, since the latter are among the main sources of scattered waves.For another topical problem, that is, radio communication effected between objects, the most difficult are the issues of designing of antennas arranged on an object, since their operating efficiency is closely related to the geometric and radiophysical properties of its surface.The issues of diffraction of an electromagnetic wave in wedge-shaped regions are the focus of numerous of the problems for a perfectly conducting and impedance wedge for monochromatic waves is representation of the diffraction field in an angular region in the form of a Sommerfeld integral [19].Substitution of Sommerfeld integrals into the system of boundary conditions gives a system of recurrence functional equations for unknown analytical integrands. The system’s coefficients are Fresnel coefficients defining the reflection of plane media or their refraction into the opposite medium. From the system of functional equations, one determines, in a recurrence manner, sequences of integrand poles and residues in these poles.The edge diffraction field in both media is determined using a pair of Fredholm-type singular integral equations of the second kind which are obtained from the above-indicated systems of functional equations with subsequent computation of Sommerfeld integrals by the saddle-point approximation. The branching points of the integrands condition the presence of creeping waves excited by the edge of the dielectric wedge.The proposed method is only true of monochromatic waves and of the approximate Leontovich boundary conditions, when the field of the electromagnetic wave slowly varies from point to point on a wavelength scale [20].We note that the existing approximate Leontovich conditions have a number of other constants and should be used with caution [20].In actual fact, the proposed calculation method does not work in the presence of, for example, two wedges, when the sharp angles are pointed at each other, that is, an optical knife, or in diffraction of the electromagnetic wave on a system of parallel lobes, when the gap between the lobes is in the region of microns, and the electromagnetic field is strongly “cut” throughout the space with a step much than the wavelength.(A) Optical Knife
Figure7 shows the field of an electromagnetic wave in its diffraction on the optical knife. The parameters of the wave at entry and at exit are Ex=104sin⁡(1010t),Ey=104cos⁡(1010t).
The electrophysical characteristics are as follows. The wedge is manufactured from aluminum:ε=1; μ=1; σ=3.774·107 S/m; the ambient medium is air.
The dimensions of the computational domain are0.1×0.05 m. The calculation time 10-9 sec, and the time step is 10-11 sec.
Numerical solution of the system of equations (13) yields the dependences of the distribution Ex,Ey on the optical knife. The calculation results are in good agreement with the existing experimental data and experiments specially conducted at the Department of the Physics and Chemistry of Nonequilibrium Media of the A.V. Luikov Heat and Mass Transfer Institute of the National Academy of Sciences by A. I. Bereznyak.
The experiments were carried out with an optical-range laser and were tentative in character but the obtained experimental photographs of diffraction fields, and the calculated results turned out to be in good qualitative agreement. The authors express their thanks to A. I. Bereznyak for the conducting of the experiments.Mesh, amplitudeEx  and isolines, and amplitude Ey  and isolines of the electromagnetic field strength.
(a)
Mesh(b)
AmplitudeEx(c)
AmplitudeEx(d)
IsolinesEx(e)
AmplitudeEy(f)
AmplitudeEy(g)
IsolinesEy(B) Diffraction Grating
The parameters of the wave and the interfacial conditions are the same, as those for the case “optical knife.” The electrophysical characteristics are as follows: 2D lobes,ε=12; and σ=100 S/m; the ambient medium is air; the characteristics of the prism and the square are identical to those of the lobes.
Figure8 corresponds to a calculation time of 10−10 sec; the time step is 10−12 sec. Figure 9 corresponds to a calculation time of 10−9 sec; the time step is 10−11 sec.
It is seen from the modeling results that the proposed “comb” cab be used as a filter of a high-frequency signal. Furthermore, we carried out numerical calculations of a modulated signal at a frequency of 20 kHz. The results of the modulated-signal calculations are not given. To analyze the difference scheme for stability was analyzed by the initial data. When the time and space steps are large there appear oscillations of the grid solution and of its “derivatives” (“ripple”) which strongly decrease the accuracy of the scheme. Undoubtedly, this issue calls for separate consideration. The proposed algorithm of solution of Maxwell equations allows circuitry-engineering modeling of high-frequency radio-engineering devices and investigation of the propagation of electromagnetic waves in media of intricate geometry in the presence of strong discontinuities of electromagnetic field.
The result of Para 2.1 were published in part [21].AmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEyAmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEy


#### 2.1.5. Conclusions

We were the first to construct a consistent physicomathematical model of propagation of electromagnetic waves in layered media without recourse to the matrices of the induced-surface-charge impedances. This model is based on the Maxwell equations, the electric-charge conservation law, the total-current continuity, and the Dirichlet theorem. Our numerical investigations have shown that the physical and mathematical model proposed can be used to advantage for simulation of the propagation of a high-frequency electromagnetic wave in a medium consisting of layers having different electrophysical properties.


### 2.2. Wave Equation forH⇀and Conditions on the Boundaries in the Presence of Strong Discontinuities of the Electromagnetic Field. Numerical Modeling of Electrodynamic Processes in the Surface Layer

#### 2.2.1. Introduction

During the interaction of an external magnetic field and magnetic abrasive particles, the particles are magnetized, and magnetic dipoles with the moment oriented predominantly along the field are formed. “Chains" along the force lines of the field [22, 23] appear that periodically act on the processable surface with a frequency ω=l/υ. A fixed elemental area of the material periodically experiences the effect of the magnetic field of one direction. Actually, the frequency and duration of the pulse will be still higher because of the rotation of the magnetic abrasive particle due to the presence of the moment of forces on contact and of the friction of the particle against the processable part. In what follows, we will not take into account the effect of rotation.We assume that the particle velocity on the polisher isυ. If the particle radius is r, then the angular frequency is ω=2πυ/r, and precisely this frequency determines the frequency of the effect of the variable magnetic field component due to the fact that for a ferromagnetic μ>1. The magnetic permeability μ of ferromagnetics, which are usually used in magnetoabrasive polishing, is measured by thousands of units in weak fields. However, in polishing, the constant external magnetic field is strong and amounts to 105-106 A/m, and in this case the value of μ for compounds of iron and nickel and for Heusler alloy decreases substantially.Because of the presence of a strong external magnetic fieldH0 the “small" absolute value of μ of an abrasive particle leads to a periodic “increase" and “decrease" in the normal component of the magnetic induction near the processable surface. In the present work we used neodymium magnets (neodymium-iron-boron) with H0>485,000 A/m. The magnetic permeability of a magnetic abrasive particle based on carbonyl iron was assumed in this case to be equal to μ1=100.Due to the continuity of the normal magnetic induction componentBn1=Bn2, where Bn1=μ1μ0H1; Bn2=μ2μ0H2. For example, in glasses (μ2=1; therefore at the boundary of contact of the glass with the magnetic abrasive particle an additional variable magnetic field of strength H1>H0 appears.In [24–27], magnetic field-induced effects in silicon are considered: a nonmonotonic change in the crystal lattice parameters in the surface layer of silicon, the gettering of defects on the surface, the change in the sorption properties of the silicon surface, and the change in the mobility of the edge dislocations and in the microhardness of silicon.In [28–30], the influence of an electromagnetic field on the domain boundaries, plasticity, strengthening, and on the reduction of metals and alloys was established.In view of the foregoing, it is of interest to find the relationship between the discrete-impulse actions of a magnetic field of one direction on the surface layer of the processable material that contains domains. According to [22], the size of domains is as follows: 0.05 μm in iron, 1.5 μm in barium ferrite; 8 μm in the MnBi compound, and 0.5–1 μm in the acicular gamma ferric oxide. According to [31], the size of a domain may reach 10−6 cm3 (obtained by the method of magnetic metallography).As a rule, an abrasive exhibits a distinct shape anisotropy, whereas the frequency of the effect is determined by the concentration of abrasive particles in a hydrophobic solution and by the velocity of its motion. We assume that on the surface of a processable crystal the magnetic field strengthH(t)=H1sin⁡4(ωt)+H0.It is required to find the value of the magnetic field strength in the surface layer that has the characteristicsλ1, ε1, and μ1 and contains domains with electrophysical properties λ2, ε2, and μ2. The domains may have the form of a triangular prism, a bar, a cylinder, and so forth.


#### 2.2.2. Physicomathematical Model: Wave Equation forH⇀

We consider interaction of a nonstationary magnetic field with substance on an example of a specific technology of magneto abrasive polishing based on use of magneto-rheological suspension. We will formulate a physicomathematical model of propagation of electromagnetic waves in a heterogeneous medium. The media in contact are considered homogeneous. We operate with the operator rot on the left- and right-hand sides of the first equation for the total current (6) and multiply by μ0μ; then we differentiate the second equation in (7) with respect to time. Taking into consideration the solenoidality of the magnetic field (7) and the rule of repeated application of the operator ∇ to the vector  H⇀, we obtain(35)μ0εε0∂2H⇀∂t2+λμ0∂H⇀∂t=1μ∇2H⇀.
In the Cartesian coordinates,(35) will have the form(36)εε0∂2Hx∂t2+λμ0∂Hx∂t=1μ(∂2Hx∂x2+∂2Hx∂y2+∂2Hx∂z2),εε0∂2Hy∂t2+λμ0∂Hy∂t=1μ(∂2Hy∂x2+∂2Hy∂y2+∂2Hy∂z2),εε0∂2Hz∂t2+λμ0∂Hz∂t=1μ(∂2Hz∂x2+∂2Hz∂y2+∂2Hz∂z2).
One fundamental electromagnetic field equation is the equation div⁡B⇀=0. The use of the Dirichlet theorem for approximation of the value of the magnetic field strength on the boundaries between adjacent media analogously to that of the electric field strength does not necessarily guarantees the observance of the condition of solenoidality of the magnetic field; furthermore, the magnetic properties of heterogeneous media were assumed constant in deriving generalized wave equations. The experience of numerical calculations has shown that when it is necessary to model nonstationary magnetic phenomena it is better in many cases to use a generalized wave equation for E⃗, accordingly expressing H⃗(t,r⃗) by E⃗(t,r⃗) and, if need be, to perform backward recalculation to H⃗(t,r⃗). This approach is difficult to apply to modeling of heterogeneous media with different magnetic properties, when the magnetic permeability μ is dependent on coordinates.In media with a weak heterogeneity whereμ(x,y,z) is a piecewise continuous quantity, the application of the proposed method of through counting is quite justified. Indeed, the system of (13) and (36) yields that the function’s discontinuity on the boundaries between adjacent media is determined by the complexes which will be called the generalized permeability ε*≡εμε0μ0 and the generalized conductivity λ*≡λμμ0. Using the Direchlet theorem for ε* and λ*, we obtain their values on the boundaries between adjacent media and the values for the electric field strength at the discontinuity point (see (34)); here, we note that the value of the electric field strength is obtained without solving Maxwell equations. In fact, at the discontinuity point, we use linear interpolation of the function to obtain the values of ε*, λ*, and Ex=ξ=(1/2)[E(ξ-0)+E(ξ+0)]. Consequently, for piecewise continuous quantity μ(x,y,z), the application of the proposed method of through counting is justified. We note that the equality of the derivatives of the electric field strength along the normal to the surface at the discontinuity point according to Equation (32) holds. When the wave equation for H⃗ is used for media with different magnetic permeabilities the condition of equality of the derivatives fails, that is,(37)∂Hx∂x|x=ξ-0≠∂Hx∂x|x=ξ+0,
which is a consequence of (10); therefore, the use of through-counting schemes for the wave equation for H⃗ is difficult.Moreover, reformulation of the boundary conditions for equations (6) and (7), with electric field strength being eliminated from the boundary conditions, does not even allow prediction of the surface current direction. Indeed, magnetic field strength H⇀ is an axial vector which on substitution r⃗→(-r⃗)  must remain constant. The value of the surface current also cannot be determined without knowledge of the polar vector of electric field strength. The generalized wave equation for E⃗ contains the term graddiv⁡E⃗ which directly allows for the influence of induced surface charges on the propagation of waves,  the right-hand side of the generalized wave equation has the same form as the equations of the linear elasticity theory, hydrodynamics. We note that the proposed method of calculation can be used on condition that there are no built-in space charges and extraneous electromotive forces [13].By virtue of what has been stated above, for modeling of the propagation of electromagnetic waves in glasses having roughness and defects, we used system (13) with boundary conditions (22)–(31) and (32).


#### 2.2.3. Results of Numerical Simulation

The physicomathematical model developed can also efficiently be used in modeling the propagation of electromagnetic waves in media with complex geometries and strong electromagnetic field discontinuities.The transverse cut of a cellular structure represents a set of parallelepipeds and triangular prisms of various cross-sections, as depicted in Figures 10(a) and11(a). An electromagnetic wave propagates across the direction of parallelepipeds and triangular prisms (channels) along the coordinate x.The size of the investigated two-dimensional object is14×20·10-6 m, and the sizes of the domains are 2–4 μm. The frequency of the influence of the magnetic field is ω=2π·106, and the strength of the field is(38)Hx=21⋅105sin4(2π⋅106t)Am.
The electrophysical properties are as follows: of the large parallelepiped, μ=1, ε=8, σ=10-9 Ω·m; of domains, μ=1, ε=6, σ=10-8 Ω·m. They correspond to the electrophysical properties of glasses.It was assumed that in a layer of thickness 15–20μm an electromagnetic wave propagates without attenuation; therefore, on all the faces of the large parallelepiped the fulfillment of condition (38) was considered valid. On the faces of the parallelepiped that are parallel to the OX axis condition (38) corresponded to the “transverse" tangential component of the wave; on the faces parallel to OY condition (38) corresponded to the normal component of the field.The calculations were carried out with a time step of 10−13 sec up to a time instant of 10−10 sec.Figures10(a) and 11(a) present the amplitude values of the magnetic field strength along Hx and Hy with a comparison scale, whereas Figures 11(b) and 11(b) present the corresponding isolines. An analysis of these figures shows that at the places of discontinuity, on the wedges, force lines of the electromagnetic field concentrate. According to [32], precisely wedges are often the sources and sinks of the vacancies that determine, for example, the hardness and plasticity of a solid body.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy and isolines of the magnetic field strength.
(a)
AmplitudeHy(b)
IsolinesAlso, we modeled the propagation of waves in media, when domains possess magnetic properties. We assumed, in the calculations, thatμ=100; the remaining parameters correspond to the previous example of solution (Figure 12).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyOf interest is the interaction of the electromagnetic wave with the rough surface shown in Figures13 and 14. As in the previous examples, we observe the concentration of electromagnetic energy on angular structures.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy  and isolines  of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesFrom Figures13 and 14, it is seen that electromagnetic heating of tapered structures may occur in addition to mechanical heating in magnetic abrasive machining.As we have mentioned above, for investigation of the propagation of electromagnetic waves in nonmagnetic materials, it is more expedient to use the generalized equation forE⃗. For the purpose of illustration we give an example of numerical calculation of an optical knife with the wave equation for H⃗ (Figure 15).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyFrom Figure15, it is seen that the actual problem of diffraction on the optical knife remains to be solved, that is, there is no “glow” on the optical-knife section, which is inconsistent with experimental data.As is known [33] in thermodynamically equilibrium systems the temperature T and the electrical φ and chemical μc potentials are constant along the entire system:(39)gradT=0,gradφ=0,gradμc=0.
If these conditions are not fulfilled (gradT≠0,gradφ≠0,gradμc≠0), irreversible processes of the transfer of mass, energy, electrical charge, and so forth appear in the system.The chemical potential of thejth component is determined, for example, as a change of the free energy with a change in the number of moles:(40)μcj=(∂F∂nj)T,V,
where
(41)dF=-SdT-PdV+H⇀dB⇀+EdD.
The last term in (41) takes into account the change in the free energy of a dielectric due to the change in the magnetic induction. The free energy of a unit volume of the dielectric in the magnetic field in this case has the form(42)F(T,D)=F0+μμ0H2⇀2+εε0E22.
We assume that changes in the temperature and volume of the dielectric are small. Then the mass flux is determined by a quantity proportional to the gradient of the chemical potential or, according to (41), we obtain(43)qi=-Dμcgrad(H⇀dB⇀+EdD)=-DμcgradW,
where W=μμ0(H2⇀/2)+εε0(E2/2) is the density of the magnetic field in the unit volume of the dielectric.In magnetic abrasive polishing on the sharp protrusions of domains the gradients of magnetic energy are great, which can lead to the origination of vacancy flows.An analysis of the results shows that the nonstationary component of the full electromagnetic energy is also concentrated in the region of fractures and wedges, that is, at the sharp angles of domains, which may lead to the improvement of the structure of the sublayer of the treated surface due to the “micromagnetoplastic" effect. Maximum values of the nonstationary part of the total electromagnetic energyWmax⁡ in the sublayer corresponds to a maximum value of the function sin (2π·106t)and occur for the time instants t=(n/4)10-6 sec, where n is the integer, with the value of Wmax⁡ for a neodymium magnet and a magnetoabrasive particle on the basis of carboxyl iron amounting to a value of the order of ((5)-(6))·106 J/m3. Having multiplied Wmax⁡ by the volume of a domain, vacancy, or atom, we may approximately obtain the corresponding energy. The density of the electromagnetic energy in all of the cases is much smaller than the bonding energy of atoms, 10−18-10−19 J. However, a periodic change in the magnetic field in one direction leads to a ponderomotive force that may influence the motion of various defects and dislocations to create a stable and equilibrium structure of atoms and molecules in magnetic abrasive polishing and, on the long run, in obtaining a surface with improved characteristics due to the “micromagnetoplastic" effect. The result of Para 2.2 were published in part [34].


## 2.1. Generalized Wave Equation for E and Conditions on the Boundaries in the Presence of Strong Discontinuities of the Electromagnetic Field

### 2.1.1. Physicomathematical Model

We will formulate a physicomathematical model of propagation of an electromagnetic field in a heterogeneous medium. Let us multiply the left and right sides of the equation for the total current (see (6)) by μμ0 and differentiate it with respect to time. Acting by the operator rot on the left and right sides of the first equation of (7) on condition that μ=const we obtain(12)∂jtotal∂t=1μμ0∇2E-1μμ0grad(div⁡E).
In Cartesian coordinates, (12) will take the form(13)∂jtotalx∂t=1μμ0(∂2Ex∂x2+∂2Ex∂y2+∂2Ex∂z2)-1μμ0∂∂x(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotaly∂t=1μμ0(∂2Ey∂x2+∂2Ey∂y2+∂2Ey∂z2)-1μμ0∂∂y(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotalz∂t=1μμ0(∂2Ez∂x2+∂2Ez∂y2+∂2Ez∂z2)-1μμ0∂∂z(∂Ex∂x+∂Ey∂y+∂Ez∂z).
At the interface, the following relation [2] is also true: (14)div⁡iτ+Iqx1-Iqx2=-∂σ∂t.
Let us write conditions (see (8)–(11)) in the Cartesian coordinate system: (15)Dx1-Dx2=σ,(16)Ey1-Ey2=0,(17)Ez1-Ez2=0,(18)Bx1-Bx2=0,(19)Hy1-Hy2=iz,(20)Hz1-Hz2=iy,
where iτ=iyj+izk is the surface-current density, and the coordinate x is directed along the normal to the interface. The densities iy and iz of the surface currents represent the electric charge carried in unit time by a segment of unit length positioned on the surface drawing the current perpendicularly to its direction.The order of the system of differential equations (13) is equal to (16). Therefore, at the interface S, it is necessary to set, by and large, nine boundary conditions. Moreover, the three additional conditions (15), (19), and (20) containing (prior to the solution) unknown quantities should be fulfilled at this interface. Consequently, the total number of conjugation conditions at the boundary S should be equal to (12) for a correct solution of the problem.As an example we consider a one-dimensional case of heat propagation in the layered material. As is known, to solve the problem we must specify two conditions at the interface: the equality of heat fluxes and the equality of temperatures, since in the heat conduction equation we have the second derivative with respect to the coordinate.Differentiating expression (15) with respect to time and using relation (14), we obtain the following condition for the normal components of the total current at the medium-medium interface:(21)div⁡iτ+jtotalx1=jtotalx2,
that allows one to disregard the surface charge σ. Let us introduce the arbitrary function f:[f]|x=ξ=f1|x=ξ+0-f2|x=ξ-0. In this case, expression (21) will take the form(22)[div⁡iτ+jtotalx1]|x=ξ=0.
It is assumed that, at the medium-medium interface, Ex is a continuous function of y and z. Then, differentiating Equation (21) with respect to y and z, we obtain(23)[∂∂yjtotalx1]|x=ξ=-∂(div⁡iτ)∂y,[∂∂zjtotalx1]|x=ξ=-∂(div⁡iτ)∂z.
Let us differentiate conditions (18)–(20) for the magnetic induction and the magnetic-field strength with respect to time. On condition that B⇀=μμ0H⇀:(24)[∂Bx∂t]|x=ξ=0,[1μμ0∂By∂t]|x=ξ=∂iz∂t,[1μμ0∂Bz∂t]|x=ξ=∂iy∂t.
Using (7) and expressing (24) in terms of projections of the electric-field rotor, we obtain(25)[rot⁡xE]|x=ξ=0,[∂Ez∂y-∂Ey∂z]|x=ξ=0,(26)[1μμ0rot⁡yE]|x=ξ=∂iz∂t,or[1μμ0(∂Ex∂z-∂Ez∂x)]|x=ξ=∂iz∂t,(27)[1μμ0rot⁡zE]|x=ξ=∂iz∂t,or[1μμ0(∂Ey∂x-∂Ex∂y)]|x=ξ=∂iz∂t.
Here, (25) is the normal projection of the electric-field rotor, (26) is the tangential projection of the rotor on y, and (27) is the rotor projection on z.Assuming thatEy and Ez are continuous differentiable functions of the coordinates y and z, from conditions (16) and (17) we find(28)[∂Ey∂y]|x=ξ=0,[∂Ey∂z]|x=ξ=0,[∂Ez∂y]|x=ξ=0,[∂Ez∂z]|x=ξ=0.
In accordance with the condition that the tangential projections of the electric field on z and y are equal and in accordance with conditions (16) and (17), the expressions for the densities of the surface currents iz and iy take the form(29)iz=λ̅Ez|x=ξ,iy=λ̅Ey|x=ξ,
where
(30)λ̅=12(λ1+λ2)|x=ξ
is the average value of the electrical conductivity at the interface between the adjacent media in accordance with the Dirichlet theorem for a piecewise-smooth, piecewise-differentiable function.Consequently, formulas (28)–(30) yield(31)[div⁡iτ⃗]|x=ξ=0.
Relations (31) and (32) and hence the equality of the normal components of the total current were obtained (in a different manner) by Grinberg and Fok [13]. In this work, it has been shown that condition (31) leads to the equality of the derivatives of the electric field strength along the normal to the surface(32)[∂Ex∂x]|x=ξ=0.
With allowance for the foregoing we have twelve conditions at the interface between the adjacent media that are necessary for solving the complete system of (13):(a)
the functionsEy and Ez are determined from (16) and (17);(b)
Ex is determined from condition (22);(c)
the values of∂Ex/∂y, ∂Ex/∂z, and ∂Ex/∂x are determined from relations (23) with the use of the condition of continuity of the total-current normal component at the interface (22) and the continuity of the derivative of the total current with respect to the coordinate x;(d)
the values of∂Ey/∂y, ∂Ey/∂z, and ∂Ez/∂z are determined from conditions (28) and (29) in consequence of the continuity of the tangential components of the electric field along y and z;(e)
the derivatives∂Ey/∂x and ∂Ez/∂x are determined from conditions (26) and (27) as a consequence of the equality of the tangential components of the electric-field rotor along y and z.Note that condition (21) was used by us in [14] in the numerical simulation of the pulsed electrochemical processes in the one-dimensional case. Condition (25) for the normal component of the electric-field rotor represents a linear combination of conditions (28) and (29); therefore, rot⁡xE=0 and there is no need to use it in the subsequent discussion. The specificity of the expression for the general law of electric-charge conservation at the interface is that the components ∂Ey/∂y and ∂Ez/∂z are determined from conditions (28) and (29) that follow from the equality and continuity of the tangential components Ey and Ez at the boundary S.Thus, at the interface between the adjacent media the following conditions are fulfilled: the equality of the total-current normal components; the equality of the tangential projections of the electric-field rotor; the electric-charge conservation law; the equality of the electric-field tangential components and their derivatives in the tangential direction; the equality of the derivatives of the total-current normal components in the direction tangential to the interface between the adjacent media, determined with account for the surface currents and without explicit introduction of a surface charge. They are true at each cross section of the sample being investigated.


### 2.1.2. Features of Calculation of the Propagation of Electromagnetic Waves in Layered Media

The electromagnetic effects arising at the interface between different media under the action of plane electromagnetic waves have a profound impact on the equipment because all real devices are bounded by the surfaces and are inhomogeneous in space. At the same time, the study of the propagation of waves in layered conducting media and, according to [15], in thin films is reduced to the calculation of the reflection and transmission coefficients; the function E(x) is not determined in the thickness of a film, that is, the geometrical-optics approximation is used.The physicomathematical model proposed allows one to investigate the propagation of an electromagnetic wave in a layered medium without recourse to the assumptions used in [3–5, 7].Since conditions (21)–(29) are true at each cross section of a layered medium, we will use schemes of through counting without an explicit definition of the interface between the media. In this case, it is proposed to calculate Ex at the interface in the following way.In accordance with (15), Ex1≠Ex2, that is, Ex(x) experiences a discontinuity of the first kind. Let us determine the strength of the electric field at the discontinuity point x=ξ on condition that Ex(x) is a piecewise-smooth, piecewise-differentiable function having finite one-sided derivatives Ex+′(x) and Ex-′(x). At the discontinuity points xi,(33)Ex+′(xi)=lim⁡Δxi→+0E(xi+Δxi)-E(xi+0)Δxi,Ex-′(xi)=lim⁡Δxi→-0E(xi+Δxi)-E(xi-0)Δxi.
In this case, in accordance with the Dirichlet theorem [16], the Fourier series of the function E(x) at each point x, including the discontinuity point ξ, converges and its sum is equal to(34)Ex=ξ=12[E(ξ-0)+E(ξ+0)].
The Dirichlet condition (34) also has a physical meaning. In the case of contact of two solid conductors, for example, dielectrics or electrolytes in different combinations (metal-electrolyte, dielectric-electrolyte, metal-vacuum, and so on), at the interface between the adjacent media there always arises an electric double-layer (EDL) with an unknown (as a rule) structure that, however, substantially influences the electrokinetic effects, the rate of the electrochemical processes, and so on. It is significant that, in reality, the electrophysical characteristics λ, ε, and E(x) change uninterruptedly in the electric double-layer; therefore, (34) is true for the case where the thickness of the electric double-layer, that is, the thickness of the interphase boundary, is much smaller than the characteristic size of a homogeneous medium. In a composite, for example, in a metal with embedments of dielectric balls, where the concentration of both components is fairly large and their characteristic sizes are small, the interphase boundaries can overlap and condition (34) can break down.If the thickness of the electric double-layer is much smaller than the characteristic sizeL of an object, (34) also follows from the condition that E(x) changes linearly in the EDL region. In reality, the thickness of the electric double-layer depends on the kind of contacting materials and can comprise several tens of angstroms [17]. In accordance with the modern views, the outer coat of the electric double-layer consists of two parts, the first of which is formed by the ions immediately attracted to the surface of the metal (a “dense” or a “Helmholtz” layer of thickness h), and the second is formed by the ions separated by distances larger than the ion radius from the surface of the layer, and the number of these ions decreases as the distance between them and the interface (the “diffusion layer”) increases. The distribution of the potential in the dense and diffusion parts of the electric double-layer is exponential in actual practice [17], that is, the condition that E(x) changes linearly breaks down; in this case, the sum of the charges of the dense and diffusion parts of the outer coat of the electric double-layer is equal to the charge of its inner coat (the metal surface). However, if the thickness of the electric double-layer h is much smaller than the characteristic size of an object, the expansion of E(x) into a power series is valid and one can restrict oneself to the consideration of a linear approximation. In accordance with the more general Dirichlet theorem (1829), a knowledge of this function in the EDL region is not necessary to substantiate (34). Nonetheless, the above-indicated physical features of the electric double-layer lend support to the validity of condition (34).The condition at interfaces, analogous to (34), has been obtained earlier [18] for the potential field (where rot⁡E=0) on the basis of introduction of the surface potential, the use of the Green formula, and the consideration of the discontinuity of the potential of the double-layer. In [18], it is also noted that the consideration of the thickness of the double-layer and the change in its potential at h/L≪1 makes no sense in general; therefore, it is advantageous to consider, instead of the volume potential, the surface potential of any density. Condition (34) can be obtained, as was shown in [16], from the more general Dirichlet theorem for a nonpotential vorticity field [18].Thus, the foregoing and the validity of conditions (15)–(17) and (23)–(29) at each cross-section of a layered medium show that, for numerical solution of the problem being considered it is advantageous to use schemes of through counting and make the discretization of the medium in such a way that the boundaries of the layers have common points.The medium was divided into finite elements so that the nodes of a finite-element grid, lying on the separation surface between the media with different electrophysical properties, were shared by these media at a time. In this case, the total currents or the current flows at the interface should be equal if the Dirichlet condition (34) is fulfilled.


### 2.1.3. Results of Numerical Simulation of the Propagation of Electromagnetic Waves in Layered Media

Let us analyze the propagation of an electromagnetic wave through a layered medium that consists of several layers with different electrophysical properties in the case where an electromagnetic-radiation source is positioned on the upper plane of the medium. It is assumed that the normal component of the electric-field vectorEx=0 and its tangential component Ey=asin⁡(ωt), where a is the electromagnetic-wave amplitude (Figure 2).Figure 2
Scheme of a layered medium: layers 1, 3, and 5 are characterized by the electrophysical parametersε1, λ1,and μ1, and layers 2, 4, and 6 by ε2, λ2, and μ2.In this example, for the purpose of correct specification of the conditions at the lower boundary of the medium, an additional layer is introduced downstream of layer 6; this layer has a larger conductivity and, therefore, the electromagnetic wave is damped out rapidly in it. In this case, the conditionEy=Ez=0 can be set at the lower boundary of the medium. The above manipulations were made to limit the size of the medium being considered because, in the general case, the electromagnetic wave is attenuated completely at an infinite distance from the electromagnetic-radiation source.Numerical calculations of the propagation of an electromagnetic wave in the layered medium with electrophysical parametersε1=ε2=1, λ1=100, λ2=1000, and μ1=μ2=1 were carried out. Two values of the cyclic frequency ω=2π/T were used: in the first case, the electromagnetic-wave frequency was assumed to be equal to ω=1014 Hz (infrared radiation), and, in the second case, the cyclic frequency was taken to be ω=109 Hz (radiofrequency radiation).As a result of the numerical solution of the system of (13) with the use of conditions S (22)–((31) and (32) at the interfaces, we obtained the time dependences of the electric-field strength at different distances from the surface of the layered medium (Figure 3).Figure 3
Time change in the tangential component of the electric-field strength at a distance of 1μm (1), 5 μm (2), and 10 μm (3) from the surface of the medium at λ1=100, λ2=1000, ε1=ε2=1, μ1=μ2=1, and ω=1014 Hz. t, sec.The results of our simulation (Figure4) have shown that a high-frequency electromagnetic wave propagating in a layered medium is damped out rapidly, whereas a low-frequency electromagnetic wave penetrates into such a medium to a greater depth. The model developed was also used for calculating the propagation of a modulated signal of frequency 20 kHz in a layered medium. As a result of our simulation (Figure 5), we obtained changes in the electric-field strength at different depths of the layered medium, which points to the fact that the model proposed can be used to advantage for calculating the propagation of polyharmonic waves in layered media; such a calculation cannot be performed on the basis of the Helmholtz equation.Figure 4
Distribution of the amplitude of the electric-field-strength at the cross section of the layered medium:ω=1014 (1) and 109 Hz (2). y,μm.Figure 5
Time change in the electric-field strength at a distance of 1 (1), 5 (2), and 10μm (3) from the surface of the medium. t, sec.The physicomathematical model developed can also be used to advantage for simulation of the propagation of electromagnetic waves in media with complex geometric parameters and large discontinuities of the electromagnetic field (Figure6).Distribution of the amplitude of the electric-field strength in the two-dimensional medium and in depth atε1=15, ε2=20, λ1=10-6, λ2=10, μ1=μ2=1, and ω=109 Hz (the dark background denotes medium 1, and the light background denotes medium 2). x,y, mm; E, V/m.
(a)
Distribution of the amplitude of the electric-field strength in the two-dimensional medium(b)
Distribution of the amplitude of the electric-field strength in depthFigure6(a) shows the cross-sectional view of a cellular structure representing a set of parallelepipeds with different cross-sections in the form of squares. The parameters of the materials in the large parallelepiped are denoted by Index 1, and the parameters of the materials in the small parallelepipeds (the squares in the figure) are denoted by Index 2.An electromagnetic wave propagates in the parallelepipeds (channels) in the transverse direction. It is seen from Figure6(b) that, in the cellular structure there are “silence regions,” where the amplitude of the electromagnetic-wave strength is close to zero, as well as inner regions where the signal has a marked value downstream of the “silence” zone formed as a result of the interference.


### 2.1.4. Results of Numerical Simulation of the Scattering of Electromagnetic Waves in Angular Structures

It is radiolocation and radio-communication problems that are among the main challenges in the set of problems solved using radio-engineering devices.Knowledge of the space-time characteristics of diffraction fields of electromagnetic waves scattered by an object of location into the environment is necessary for solving successfully any radiolocation problem. Irradiated objects have a very intricate architecture and geometric shape of the surface consisting of smooth portions and numerous wedge-shaped for formations of different type-angular joints of smooth portions, surface fractures, sharp edges, and so forth, with rounded radii much smaller than the probing-signal wavelength. Therefore, solution of radiolocation problems requires the methods of calculation of the diffraction fields of electromagnetic waves excited and scattered by different surface portions of the objects, in particular, by wedge-shaped formations, be known, since the latter are among the main sources of scattered waves.For another topical problem, that is, radio communication effected between objects, the most difficult are the issues of designing of antennas arranged on an object, since their operating efficiency is closely related to the geometric and radiophysical properties of its surface.The issues of diffraction of an electromagnetic wave in wedge-shaped regions are the focus of numerous of the problems for a perfectly conducting and impedance wedge for monochromatic waves is representation of the diffraction field in an angular region in the form of a Sommerfeld integral [19].Substitution of Sommerfeld integrals into the system of boundary conditions gives a system of recurrence functional equations for unknown analytical integrands. The system’s coefficients are Fresnel coefficients defining the reflection of plane media or their refraction into the opposite medium. From the system of functional equations, one determines, in a recurrence manner, sequences of integrand poles and residues in these poles.The edge diffraction field in both media is determined using a pair of Fredholm-type singular integral equations of the second kind which are obtained from the above-indicated systems of functional equations with subsequent computation of Sommerfeld integrals by the saddle-point approximation. The branching points of the integrands condition the presence of creeping waves excited by the edge of the dielectric wedge.The proposed method is only true of monochromatic waves and of the approximate Leontovich boundary conditions, when the field of the electromagnetic wave slowly varies from point to point on a wavelength scale [20].We note that the existing approximate Leontovich conditions have a number of other constants and should be used with caution [20].In actual fact, the proposed calculation method does not work in the presence of, for example, two wedges, when the sharp angles are pointed at each other, that is, an optical knife, or in diffraction of the electromagnetic wave on a system of parallel lobes, when the gap between the lobes is in the region of microns, and the electromagnetic field is strongly “cut” throughout the space with a step much than the wavelength.(A) Optical Knife
Figure7 shows the field of an electromagnetic wave in its diffraction on the optical knife. The parameters of the wave at entry and at exit are Ex=104sin⁡(1010t),Ey=104cos⁡(1010t).
The electrophysical characteristics are as follows. The wedge is manufactured from aluminum:ε=1; μ=1; σ=3.774·107 S/m; the ambient medium is air.
The dimensions of the computational domain are0.1×0.05 m. The calculation time 10-9 sec, and the time step is 10-11 sec.
Numerical solution of the system of equations (13) yields the dependences of the distribution Ex,Ey on the optical knife. The calculation results are in good agreement with the existing experimental data and experiments specially conducted at the Department of the Physics and Chemistry of Nonequilibrium Media of the A.V. Luikov Heat and Mass Transfer Institute of the National Academy of Sciences by A. I. Bereznyak.
The experiments were carried out with an optical-range laser and were tentative in character but the obtained experimental photographs of diffraction fields, and the calculated results turned out to be in good qualitative agreement. The authors express their thanks to A. I. Bereznyak for the conducting of the experiments.Mesh, amplitudeEx  and isolines, and amplitude Ey  and isolines of the electromagnetic field strength.
(a)
Mesh(b)
AmplitudeEx(c)
AmplitudeEx(d)
IsolinesEx(e)
AmplitudeEy(f)
AmplitudeEy(g)
IsolinesEy(B) Diffraction Grating
The parameters of the wave and the interfacial conditions are the same, as those for the case “optical knife.” The electrophysical characteristics are as follows: 2D lobes,ε=12; and σ=100 S/m; the ambient medium is air; the characteristics of the prism and the square are identical to those of the lobes.
Figure8 corresponds to a calculation time of 10−10 sec; the time step is 10−12 sec. Figure 9 corresponds to a calculation time of 10−9 sec; the time step is 10−11 sec.
It is seen from the modeling results that the proposed “comb” cab be used as a filter of a high-frequency signal. Furthermore, we carried out numerical calculations of a modulated signal at a frequency of 20 kHz. The results of the modulated-signal calculations are not given. To analyze the difference scheme for stability was analyzed by the initial data. When the time and space steps are large there appear oscillations of the grid solution and of its “derivatives” (“ripple”) which strongly decrease the accuracy of the scheme. Undoubtedly, this issue calls for separate consideration. The proposed algorithm of solution of Maxwell equations allows circuitry-engineering modeling of high-frequency radio-engineering devices and investigation of the propagation of electromagnetic waves in media of intricate geometry in the presence of strong discontinuities of electromagnetic field.
The result of Para 2.1 were published in part [21].AmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEyAmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEy


### 2.1.5. Conclusions

We were the first to construct a consistent physicomathematical model of propagation of electromagnetic waves in layered media without recourse to the matrices of the induced-surface-charge impedances. This model is based on the Maxwell equations, the electric-charge conservation law, the total-current continuity, and the Dirichlet theorem. Our numerical investigations have shown that the physical and mathematical model proposed can be used to advantage for simulation of the propagation of a high-frequency electromagnetic wave in a medium consisting of layers having different electrophysical properties.


## 2.1.1. Physicomathematical Model

We will formulate a physicomathematical model of propagation of an electromagnetic field in a heterogeneous medium. Let us multiply the left and right sides of the equation for the total current (see (6)) by μμ0 and differentiate it with respect to time. Acting by the operator rot on the left and right sides of the first equation of (7) on condition that μ=const we obtain(12)∂jtotal∂t=1μμ0∇2E-1μμ0grad(div⁡E).
In Cartesian coordinates, (12) will take the form(13)∂jtotalx∂t=1μμ0(∂2Ex∂x2+∂2Ex∂y2+∂2Ex∂z2)-1μμ0∂∂x(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotaly∂t=1μμ0(∂2Ey∂x2+∂2Ey∂y2+∂2Ey∂z2)-1μμ0∂∂y(∂Ex∂x+∂Ey∂y+∂Ez∂z),∂jtotalz∂t=1μμ0(∂2Ez∂x2+∂2Ez∂y2+∂2Ez∂z2)-1μμ0∂∂z(∂Ex∂x+∂Ey∂y+∂Ez∂z).
At the interface, the following relation [2] is also true: (14)div⁡iτ+Iqx1-Iqx2=-∂σ∂t.
Let us write conditions (see (8)–(11)) in the Cartesian coordinate system: (15)Dx1-Dx2=σ,(16)Ey1-Ey2=0,(17)Ez1-Ez2=0,(18)Bx1-Bx2=0,(19)Hy1-Hy2=iz,(20)Hz1-Hz2=iy,
where iτ=iyj+izk is the surface-current density, and the coordinate x is directed along the normal to the interface. The densities iy and iz of the surface currents represent the electric charge carried in unit time by a segment of unit length positioned on the surface drawing the current perpendicularly to its direction.The order of the system of differential equations (13) is equal to (16). Therefore, at the interface S, it is necessary to set, by and large, nine boundary conditions. Moreover, the three additional conditions (15), (19), and (20) containing (prior to the solution) unknown quantities should be fulfilled at this interface. Consequently, the total number of conjugation conditions at the boundary S should be equal to (12) for a correct solution of the problem.As an example we consider a one-dimensional case of heat propagation in the layered material. As is known, to solve the problem we must specify two conditions at the interface: the equality of heat fluxes and the equality of temperatures, since in the heat conduction equation we have the second derivative with respect to the coordinate.Differentiating expression (15) with respect to time and using relation (14), we obtain the following condition for the normal components of the total current at the medium-medium interface:(21)div⁡iτ+jtotalx1=jtotalx2,
that allows one to disregard the surface charge σ. Let us introduce the arbitrary function f:[f]|x=ξ=f1|x=ξ+0-f2|x=ξ-0. In this case, expression (21) will take the form(22)[div⁡iτ+jtotalx1]|x=ξ=0.
It is assumed that, at the medium-medium interface, Ex is a continuous function of y and z. Then, differentiating Equation (21) with respect to y and z, we obtain(23)[∂∂yjtotalx1]|x=ξ=-∂(div⁡iτ)∂y,[∂∂zjtotalx1]|x=ξ=-∂(div⁡iτ)∂z.
Let us differentiate conditions (18)–(20) for the magnetic induction and the magnetic-field strength with respect to time. On condition that B⇀=μμ0H⇀:(24)[∂Bx∂t]|x=ξ=0,[1μμ0∂By∂t]|x=ξ=∂iz∂t,[1μμ0∂Bz∂t]|x=ξ=∂iy∂t.
Using (7) and expressing (24) in terms of projections of the electric-field rotor, we obtain(25)[rot⁡xE]|x=ξ=0,[∂Ez∂y-∂Ey∂z]|x=ξ=0,(26)[1μμ0rot⁡yE]|x=ξ=∂iz∂t,or[1μμ0(∂Ex∂z-∂Ez∂x)]|x=ξ=∂iz∂t,(27)[1μμ0rot⁡zE]|x=ξ=∂iz∂t,or[1μμ0(∂Ey∂x-∂Ex∂y)]|x=ξ=∂iz∂t.
Here, (25) is the normal projection of the electric-field rotor, (26) is the tangential projection of the rotor on y, and (27) is the rotor projection on z.Assuming thatEy and Ez are continuous differentiable functions of the coordinates y and z, from conditions (16) and (17) we find(28)[∂Ey∂y]|x=ξ=0,[∂Ey∂z]|x=ξ=0,[∂Ez∂y]|x=ξ=0,[∂Ez∂z]|x=ξ=0.
In accordance with the condition that the tangential projections of the electric field on z and y are equal and in accordance with conditions (16) and (17), the expressions for the densities of the surface currents iz and iy take the form(29)iz=λ̅Ez|x=ξ,iy=λ̅Ey|x=ξ,
where
(30)λ̅=12(λ1+λ2)|x=ξ
is the average value of the electrical conductivity at the interface between the adjacent media in accordance with the Dirichlet theorem for a piecewise-smooth, piecewise-differentiable function.Consequently, formulas (28)–(30) yield(31)[div⁡iτ⃗]|x=ξ=0.
Relations (31) and (32) and hence the equality of the normal components of the total current were obtained (in a different manner) by Grinberg and Fok [13]. In this work, it has been shown that condition (31) leads to the equality of the derivatives of the electric field strength along the normal to the surface(32)[∂Ex∂x]|x=ξ=0.
With allowance for the foregoing we have twelve conditions at the interface between the adjacent media that are necessary for solving the complete system of (13):(a)
the functionsEy and Ez are determined from (16) and (17);(b)
Ex is determined from condition (22);(c)
the values of∂Ex/∂y, ∂Ex/∂z, and ∂Ex/∂x are determined from relations (23) with the use of the condition of continuity of the total-current normal component at the interface (22) and the continuity of the derivative of the total current with respect to the coordinate x;(d)
the values of∂Ey/∂y, ∂Ey/∂z, and ∂Ez/∂z are determined from conditions (28) and (29) in consequence of the continuity of the tangential components of the electric field along y and z;(e)
the derivatives∂Ey/∂x and ∂Ez/∂x are determined from conditions (26) and (27) as a consequence of the equality of the tangential components of the electric-field rotor along y and z.Note that condition (21) was used by us in [14] in the numerical simulation of the pulsed electrochemical processes in the one-dimensional case. Condition (25) for the normal component of the electric-field rotor represents a linear combination of conditions (28) and (29); therefore, rot⁡xE=0 and there is no need to use it in the subsequent discussion. The specificity of the expression for the general law of electric-charge conservation at the interface is that the components ∂Ey/∂y and ∂Ez/∂z are determined from conditions (28) and (29) that follow from the equality and continuity of the tangential components Ey and Ez at the boundary S.Thus, at the interface between the adjacent media the following conditions are fulfilled: the equality of the total-current normal components; the equality of the tangential projections of the electric-field rotor; the electric-charge conservation law; the equality of the electric-field tangential components and their derivatives in the tangential direction; the equality of the derivatives of the total-current normal components in the direction tangential to the interface between the adjacent media, determined with account for the surface currents and without explicit introduction of a surface charge. They are true at each cross section of the sample being investigated.


## 2.1.2. Features of Calculation of the Propagation of Electromagnetic Waves in Layered Media

The electromagnetic effects arising at the interface between different media under the action of plane electromagnetic waves have a profound impact on the equipment because all real devices are bounded by the surfaces and are inhomogeneous in space. At the same time, the study of the propagation of waves in layered conducting media and, according to [15], in thin films is reduced to the calculation of the reflection and transmission coefficients; the function E(x) is not determined in the thickness of a film, that is, the geometrical-optics approximation is used.The physicomathematical model proposed allows one to investigate the propagation of an electromagnetic wave in a layered medium without recourse to the assumptions used in [3–5, 7].Since conditions (21)–(29) are true at each cross section of a layered medium, we will use schemes of through counting without an explicit definition of the interface between the media. In this case, it is proposed to calculate Ex at the interface in the following way.In accordance with (15), Ex1≠Ex2, that is, Ex(x) experiences a discontinuity of the first kind. Let us determine the strength of the electric field at the discontinuity point x=ξ on condition that Ex(x) is a piecewise-smooth, piecewise-differentiable function having finite one-sided derivatives Ex+′(x) and Ex-′(x). At the discontinuity points xi,(33)Ex+′(xi)=lim⁡Δxi→+0E(xi+Δxi)-E(xi+0)Δxi,Ex-′(xi)=lim⁡Δxi→-0E(xi+Δxi)-E(xi-0)Δxi.
In this case, in accordance with the Dirichlet theorem [16], the Fourier series of the function E(x) at each point x, including the discontinuity point ξ, converges and its sum is equal to(34)Ex=ξ=12[E(ξ-0)+E(ξ+0)].
The Dirichlet condition (34) also has a physical meaning. In the case of contact of two solid conductors, for example, dielectrics or electrolytes in different combinations (metal-electrolyte, dielectric-electrolyte, metal-vacuum, and so on), at the interface between the adjacent media there always arises an electric double-layer (EDL) with an unknown (as a rule) structure that, however, substantially influences the electrokinetic effects, the rate of the electrochemical processes, and so on. It is significant that, in reality, the electrophysical characteristics λ, ε, and E(x) change uninterruptedly in the electric double-layer; therefore, (34) is true for the case where the thickness of the electric double-layer, that is, the thickness of the interphase boundary, is much smaller than the characteristic size of a homogeneous medium. In a composite, for example, in a metal with embedments of dielectric balls, where the concentration of both components is fairly large and their characteristic sizes are small, the interphase boundaries can overlap and condition (34) can break down.If the thickness of the electric double-layer is much smaller than the characteristic sizeL of an object, (34) also follows from the condition that E(x) changes linearly in the EDL region. In reality, the thickness of the electric double-layer depends on the kind of contacting materials and can comprise several tens of angstroms [17]. In accordance with the modern views, the outer coat of the electric double-layer consists of two parts, the first of which is formed by the ions immediately attracted to the surface of the metal (a “dense” or a “Helmholtz” layer of thickness h), and the second is formed by the ions separated by distances larger than the ion radius from the surface of the layer, and the number of these ions decreases as the distance between them and the interface (the “diffusion layer”) increases. The distribution of the potential in the dense and diffusion parts of the electric double-layer is exponential in actual practice [17], that is, the condition that E(x) changes linearly breaks down; in this case, the sum of the charges of the dense and diffusion parts of the outer coat of the electric double-layer is equal to the charge of its inner coat (the metal surface). However, if the thickness of the electric double-layer h is much smaller than the characteristic size of an object, the expansion of E(x) into a power series is valid and one can restrict oneself to the consideration of a linear approximation. In accordance with the more general Dirichlet theorem (1829), a knowledge of this function in the EDL region is not necessary to substantiate (34). Nonetheless, the above-indicated physical features of the electric double-layer lend support to the validity of condition (34).The condition at interfaces, analogous to (34), has been obtained earlier [18] for the potential field (where rot⁡E=0) on the basis of introduction of the surface potential, the use of the Green formula, and the consideration of the discontinuity of the potential of the double-layer. In [18], it is also noted that the consideration of the thickness of the double-layer and the change in its potential at h/L≪1 makes no sense in general; therefore, it is advantageous to consider, instead of the volume potential, the surface potential of any density. Condition (34) can be obtained, as was shown in [16], from the more general Dirichlet theorem for a nonpotential vorticity field [18].Thus, the foregoing and the validity of conditions (15)–(17) and (23)–(29) at each cross-section of a layered medium show that, for numerical solution of the problem being considered it is advantageous to use schemes of through counting and make the discretization of the medium in such a way that the boundaries of the layers have common points.The medium was divided into finite elements so that the nodes of a finite-element grid, lying on the separation surface between the media with different electrophysical properties, were shared by these media at a time. In this case, the total currents or the current flows at the interface should be equal if the Dirichlet condition (34) is fulfilled.


## 2.1.3. Results of Numerical Simulation of the Propagation of Electromagnetic Waves in Layered Media

Let us analyze the propagation of an electromagnetic wave through a layered medium that consists of several layers with different electrophysical properties in the case where an electromagnetic-radiation source is positioned on the upper plane of the medium. It is assumed that the normal component of the electric-field vectorEx=0 and its tangential component Ey=asin⁡(ωt), where a is the electromagnetic-wave amplitude (Figure 2).Figure 2
Scheme of a layered medium: layers 1, 3, and 5 are characterized by the electrophysical parametersε1, λ1,and μ1, and layers 2, 4, and 6 by ε2, λ2, and μ2.In this example, for the purpose of correct specification of the conditions at the lower boundary of the medium, an additional layer is introduced downstream of layer 6; this layer has a larger conductivity and, therefore, the electromagnetic wave is damped out rapidly in it. In this case, the conditionEy=Ez=0 can be set at the lower boundary of the medium. The above manipulations were made to limit the size of the medium being considered because, in the general case, the electromagnetic wave is attenuated completely at an infinite distance from the electromagnetic-radiation source.Numerical calculations of the propagation of an electromagnetic wave in the layered medium with electrophysical parametersε1=ε2=1, λ1=100, λ2=1000, and μ1=μ2=1 were carried out. Two values of the cyclic frequency ω=2π/T were used: in the first case, the electromagnetic-wave frequency was assumed to be equal to ω=1014 Hz (infrared radiation), and, in the second case, the cyclic frequency was taken to be ω=109 Hz (radiofrequency radiation).As a result of the numerical solution of the system of (13) with the use of conditions S (22)–((31) and (32) at the interfaces, we obtained the time dependences of the electric-field strength at different distances from the surface of the layered medium (Figure 3).Figure 3
Time change in the tangential component of the electric-field strength at a distance of 1μm (1), 5 μm (2), and 10 μm (3) from the surface of the medium at λ1=100, λ2=1000, ε1=ε2=1, μ1=μ2=1, and ω=1014 Hz. t, sec.The results of our simulation (Figure4) have shown that a high-frequency electromagnetic wave propagating in a layered medium is damped out rapidly, whereas a low-frequency electromagnetic wave penetrates into such a medium to a greater depth. The model developed was also used for calculating the propagation of a modulated signal of frequency 20 kHz in a layered medium. As a result of our simulation (Figure 5), we obtained changes in the electric-field strength at different depths of the layered medium, which points to the fact that the model proposed can be used to advantage for calculating the propagation of polyharmonic waves in layered media; such a calculation cannot be performed on the basis of the Helmholtz equation.Figure 4
Distribution of the amplitude of the electric-field-strength at the cross section of the layered medium:ω=1014 (1) and 109 Hz (2). y,μm.Figure 5
Time change in the electric-field strength at a distance of 1 (1), 5 (2), and 10μm (3) from the surface of the medium. t, sec.The physicomathematical model developed can also be used to advantage for simulation of the propagation of electromagnetic waves in media with complex geometric parameters and large discontinuities of the electromagnetic field (Figure6).Distribution of the amplitude of the electric-field strength in the two-dimensional medium and in depth atε1=15, ε2=20, λ1=10-6, λ2=10, μ1=μ2=1, and ω=109 Hz (the dark background denotes medium 1, and the light background denotes medium 2). x,y, mm; E, V/m.
(a)
Distribution of the amplitude of the electric-field strength in the two-dimensional medium(b)
Distribution of the amplitude of the electric-field strength in depthFigure6(a) shows the cross-sectional view of a cellular structure representing a set of parallelepipeds with different cross-sections in the form of squares. The parameters of the materials in the large parallelepiped are denoted by Index 1, and the parameters of the materials in the small parallelepipeds (the squares in the figure) are denoted by Index 2.An electromagnetic wave propagates in the parallelepipeds (channels) in the transverse direction. It is seen from Figure6(b) that, in the cellular structure there are “silence regions,” where the amplitude of the electromagnetic-wave strength is close to zero, as well as inner regions where the signal has a marked value downstream of the “silence” zone formed as a result of the interference.


## 2.1.4. Results of Numerical Simulation of the Scattering of Electromagnetic Waves in Angular Structures

It is radiolocation and radio-communication problems that are among the main challenges in the set of problems solved using radio-engineering devices.Knowledge of the space-time characteristics of diffraction fields of electromagnetic waves scattered by an object of location into the environment is necessary for solving successfully any radiolocation problem. Irradiated objects have a very intricate architecture and geometric shape of the surface consisting of smooth portions and numerous wedge-shaped for formations of different type-angular joints of smooth portions, surface fractures, sharp edges, and so forth, with rounded radii much smaller than the probing-signal wavelength. Therefore, solution of radiolocation problems requires the methods of calculation of the diffraction fields of electromagnetic waves excited and scattered by different surface portions of the objects, in particular, by wedge-shaped formations, be known, since the latter are among the main sources of scattered waves.For another topical problem, that is, radio communication effected between objects, the most difficult are the issues of designing of antennas arranged on an object, since their operating efficiency is closely related to the geometric and radiophysical properties of its surface.The issues of diffraction of an electromagnetic wave in wedge-shaped regions are the focus of numerous of the problems for a perfectly conducting and impedance wedge for monochromatic waves is representation of the diffraction field in an angular region in the form of a Sommerfeld integral [19].Substitution of Sommerfeld integrals into the system of boundary conditions gives a system of recurrence functional equations for unknown analytical integrands. The system’s coefficients are Fresnel coefficients defining the reflection of plane media or their refraction into the opposite medium. From the system of functional equations, one determines, in a recurrence manner, sequences of integrand poles and residues in these poles.The edge diffraction field in both media is determined using a pair of Fredholm-type singular integral equations of the second kind which are obtained from the above-indicated systems of functional equations with subsequent computation of Sommerfeld integrals by the saddle-point approximation. The branching points of the integrands condition the presence of creeping waves excited by the edge of the dielectric wedge.The proposed method is only true of monochromatic waves and of the approximate Leontovich boundary conditions, when the field of the electromagnetic wave slowly varies from point to point on a wavelength scale [20].We note that the existing approximate Leontovich conditions have a number of other constants and should be used with caution [20].In actual fact, the proposed calculation method does not work in the presence of, for example, two wedges, when the sharp angles are pointed at each other, that is, an optical knife, or in diffraction of the electromagnetic wave on a system of parallel lobes, when the gap between the lobes is in the region of microns, and the electromagnetic field is strongly “cut” throughout the space with a step much than the wavelength.(A) Optical Knife
Figure7 shows the field of an electromagnetic wave in its diffraction on the optical knife. The parameters of the wave at entry and at exit are Ex=104sin⁡(1010t),Ey=104cos⁡(1010t).
The electrophysical characteristics are as follows. The wedge is manufactured from aluminum:ε=1; μ=1; σ=3.774·107 S/m; the ambient medium is air.
The dimensions of the computational domain are0.1×0.05 m. The calculation time 10-9 sec, and the time step is 10-11 sec.
Numerical solution of the system of equations (13) yields the dependences of the distribution Ex,Ey on the optical knife. The calculation results are in good agreement with the existing experimental data and experiments specially conducted at the Department of the Physics and Chemistry of Nonequilibrium Media of the A.V. Luikov Heat and Mass Transfer Institute of the National Academy of Sciences by A. I. Bereznyak.
The experiments were carried out with an optical-range laser and were tentative in character but the obtained experimental photographs of diffraction fields, and the calculated results turned out to be in good qualitative agreement. The authors express their thanks to A. I. Bereznyak for the conducting of the experiments.Mesh, amplitudeEx  and isolines, and amplitude Ey  and isolines of the electromagnetic field strength.
(a)
Mesh(b)
AmplitudeEx(c)
AmplitudeEx(d)
IsolinesEx(e)
AmplitudeEy(f)
AmplitudeEy(g)
IsolinesEy(B) Diffraction Grating
The parameters of the wave and the interfacial conditions are the same, as those for the case “optical knife.” The electrophysical characteristics are as follows: 2D lobes,ε=12; and σ=100 S/m; the ambient medium is air; the characteristics of the prism and the square are identical to those of the lobes.
Figure8 corresponds to a calculation time of 10−10 sec; the time step is 10−12 sec. Figure 9 corresponds to a calculation time of 10−9 sec; the time step is 10−11 sec.
It is seen from the modeling results that the proposed “comb” cab be used as a filter of a high-frequency signal. Furthermore, we carried out numerical calculations of a modulated signal at a frequency of 20 kHz. The results of the modulated-signal calculations are not given. To analyze the difference scheme for stability was analyzed by the initial data. When the time and space steps are large there appear oscillations of the grid solution and of its “derivatives” (“ripple”) which strongly decrease the accuracy of the scheme. Undoubtedly, this issue calls for separate consideration. The proposed algorithm of solution of Maxwell equations allows circuitry-engineering modeling of high-frequency radio-engineering devices and investigation of the propagation of electromagnetic waves in media of intricate geometry in the presence of strong discontinuities of electromagnetic field.
The result of Para 2.1 were published in part [21].AmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEyAmplitudeEx and isolines and amplitude Ey and isolines of the electromagnetic field strength.
(a)
AmplitudeEx(b)
IsolinesEx(c)
AmplitudeEy(d)
IsolinesEy


## 2.1.5. Conclusions

We were the first to construct a consistent physicomathematical model of propagation of electromagnetic waves in layered media without recourse to the matrices of the induced-surface-charge impedances. This model is based on the Maxwell equations, the electric-charge conservation law, the total-current continuity, and the Dirichlet theorem. Our numerical investigations have shown that the physical and mathematical model proposed can be used to advantage for simulation of the propagation of a high-frequency electromagnetic wave in a medium consisting of layers having different electrophysical properties.


## 2.2. Wave Equation forH⇀and Conditions on the Boundaries in the Presence of Strong Discontinuities of the Electromagnetic Field. Numerical Modeling of Electrodynamic Processes in the Surface Layer

### 2.2.1. Introduction

During the interaction of an external magnetic field and magnetic abrasive particles, the particles are magnetized, and magnetic dipoles with the moment oriented predominantly along the field are formed. “Chains" along the force lines of the field [22, 23] appear that periodically act on the processable surface with a frequency ω=l/υ. A fixed elemental area of the material periodically experiences the effect of the magnetic field of one direction. Actually, the frequency and duration of the pulse will be still higher because of the rotation of the magnetic abrasive particle due to the presence of the moment of forces on contact and of the friction of the particle against the processable part. In what follows, we will not take into account the effect of rotation.We assume that the particle velocity on the polisher isυ. If the particle radius is r, then the angular frequency is ω=2πυ/r, and precisely this frequency determines the frequency of the effect of the variable magnetic field component due to the fact that for a ferromagnetic μ>1. The magnetic permeability μ of ferromagnetics, which are usually used in magnetoabrasive polishing, is measured by thousands of units in weak fields. However, in polishing, the constant external magnetic field is strong and amounts to 105-106 A/m, and in this case the value of μ for compounds of iron and nickel and for Heusler alloy decreases substantially.Because of the presence of a strong external magnetic fieldH0 the “small" absolute value of μ of an abrasive particle leads to a periodic “increase" and “decrease" in the normal component of the magnetic induction near the processable surface. In the present work we used neodymium magnets (neodymium-iron-boron) with H0>485,000 A/m. The magnetic permeability of a magnetic abrasive particle based on carbonyl iron was assumed in this case to be equal to μ1=100.Due to the continuity of the normal magnetic induction componentBn1=Bn2, where Bn1=μ1μ0H1; Bn2=μ2μ0H2. For example, in glasses (μ2=1; therefore at the boundary of contact of the glass with the magnetic abrasive particle an additional variable magnetic field of strength H1>H0 appears.In [24–27], magnetic field-induced effects in silicon are considered: a nonmonotonic change in the crystal lattice parameters in the surface layer of silicon, the gettering of defects on the surface, the change in the sorption properties of the silicon surface, and the change in the mobility of the edge dislocations and in the microhardness of silicon.In [28–30], the influence of an electromagnetic field on the domain boundaries, plasticity, strengthening, and on the reduction of metals and alloys was established.In view of the foregoing, it is of interest to find the relationship between the discrete-impulse actions of a magnetic field of one direction on the surface layer of the processable material that contains domains. According to [22], the size of domains is as follows: 0.05 μm in iron, 1.5 μm in barium ferrite; 8 μm in the MnBi compound, and 0.5–1 μm in the acicular gamma ferric oxide. According to [31], the size of a domain may reach 10−6 cm3 (obtained by the method of magnetic metallography).As a rule, an abrasive exhibits a distinct shape anisotropy, whereas the frequency of the effect is determined by the concentration of abrasive particles in a hydrophobic solution and by the velocity of its motion. We assume that on the surface of a processable crystal the magnetic field strengthH(t)=H1sin⁡4(ωt)+H0.It is required to find the value of the magnetic field strength in the surface layer that has the characteristicsλ1, ε1, and μ1 and contains domains with electrophysical properties λ2, ε2, and μ2. The domains may have the form of a triangular prism, a bar, a cylinder, and so forth.


### 2.2.2. Physicomathematical Model: Wave Equation forH⇀

We consider interaction of a nonstationary magnetic field with substance on an example of a specific technology of magneto abrasive polishing based on use of magneto-rheological suspension. We will formulate a physicomathematical model of propagation of electromagnetic waves in a heterogeneous medium. The media in contact are considered homogeneous. We operate with the operator rot on the left- and right-hand sides of the first equation for the total current (6) and multiply by μ0μ; then we differentiate the second equation in (7) with respect to time. Taking into consideration the solenoidality of the magnetic field (7) and the rule of repeated application of the operator ∇ to the vector  H⇀, we obtain(35)μ0εε0∂2H⇀∂t2+λμ0∂H⇀∂t=1μ∇2H⇀.
In the Cartesian coordinates,(35) will have the form(36)εε0∂2Hx∂t2+λμ0∂Hx∂t=1μ(∂2Hx∂x2+∂2Hx∂y2+∂2Hx∂z2),εε0∂2Hy∂t2+λμ0∂Hy∂t=1μ(∂2Hy∂x2+∂2Hy∂y2+∂2Hy∂z2),εε0∂2Hz∂t2+λμ0∂Hz∂t=1μ(∂2Hz∂x2+∂2Hz∂y2+∂2Hz∂z2).
One fundamental electromagnetic field equation is the equation div⁡B⇀=0. The use of the Dirichlet theorem for approximation of the value of the magnetic field strength on the boundaries between adjacent media analogously to that of the electric field strength does not necessarily guarantees the observance of the condition of solenoidality of the magnetic field; furthermore, the magnetic properties of heterogeneous media were assumed constant in deriving generalized wave equations. The experience of numerical calculations has shown that when it is necessary to model nonstationary magnetic phenomena it is better in many cases to use a generalized wave equation for E⃗, accordingly expressing H⃗(t,r⃗) by E⃗(t,r⃗) and, if need be, to perform backward recalculation to H⃗(t,r⃗). This approach is difficult to apply to modeling of heterogeneous media with different magnetic properties, when the magnetic permeability μ is dependent on coordinates.In media with a weak heterogeneity whereμ(x,y,z) is a piecewise continuous quantity, the application of the proposed method of through counting is quite justified. Indeed, the system of (13) and (36) yields that the function’s discontinuity on the boundaries between adjacent media is determined by the complexes which will be called the generalized permeability ε*≡εμε0μ0 and the generalized conductivity λ*≡λμμ0. Using the Direchlet theorem for ε* and λ*, we obtain their values on the boundaries between adjacent media and the values for the electric field strength at the discontinuity point (see (34)); here, we note that the value of the electric field strength is obtained without solving Maxwell equations. In fact, at the discontinuity point, we use linear interpolation of the function to obtain the values of ε*, λ*, and Ex=ξ=(1/2)[E(ξ-0)+E(ξ+0)]. Consequently, for piecewise continuous quantity μ(x,y,z), the application of the proposed method of through counting is justified. We note that the equality of the derivatives of the electric field strength along the normal to the surface at the discontinuity point according to Equation (32) holds. When the wave equation for H⃗ is used for media with different magnetic permeabilities the condition of equality of the derivatives fails, that is,(37)∂Hx∂x|x=ξ-0≠∂Hx∂x|x=ξ+0,
which is a consequence of (10); therefore, the use of through-counting schemes for the wave equation for H⃗ is difficult.Moreover, reformulation of the boundary conditions for equations (6) and (7), with electric field strength being eliminated from the boundary conditions, does not even allow prediction of the surface current direction. Indeed, magnetic field strength H⇀ is an axial vector which on substitution r⃗→(-r⃗)  must remain constant. The value of the surface current also cannot be determined without knowledge of the polar vector of electric field strength. The generalized wave equation for E⃗ contains the term graddiv⁡E⃗ which directly allows for the influence of induced surface charges on the propagation of waves,  the right-hand side of the generalized wave equation has the same form as the equations of the linear elasticity theory, hydrodynamics. We note that the proposed method of calculation can be used on condition that there are no built-in space charges and extraneous electromotive forces [13].By virtue of what has been stated above, for modeling of the propagation of electromagnetic waves in glasses having roughness and defects, we used system (13) with boundary conditions (22)–(31) and (32).


### 2.2.3. Results of Numerical Simulation

The physicomathematical model developed can also efficiently be used in modeling the propagation of electromagnetic waves in media with complex geometries and strong electromagnetic field discontinuities.The transverse cut of a cellular structure represents a set of parallelepipeds and triangular prisms of various cross-sections, as depicted in Figures 10(a) and11(a). An electromagnetic wave propagates across the direction of parallelepipeds and triangular prisms (channels) along the coordinate x.The size of the investigated two-dimensional object is14×20·10-6 m, and the sizes of the domains are 2–4 μm. The frequency of the influence of the magnetic field is ω=2π·106, and the strength of the field is(38)Hx=21⋅105sin4(2π⋅106t)Am.
The electrophysical properties are as follows: of the large parallelepiped, μ=1, ε=8, σ=10-9 Ω·m; of domains, μ=1, ε=6, σ=10-8 Ω·m. They correspond to the electrophysical properties of glasses.It was assumed that in a layer of thickness 15–20μm an electromagnetic wave propagates without attenuation; therefore, on all the faces of the large parallelepiped the fulfillment of condition (38) was considered valid. On the faces of the parallelepiped that are parallel to the OX axis condition (38) corresponded to the “transverse" tangential component of the wave; on the faces parallel to OY condition (38) corresponded to the normal component of the field.The calculations were carried out with a time step of 10−13 sec up to a time instant of 10−10 sec.Figures10(a) and 11(a) present the amplitude values of the magnetic field strength along Hx and Hy with a comparison scale, whereas Figures 11(b) and 11(b) present the corresponding isolines. An analysis of these figures shows that at the places of discontinuity, on the wedges, force lines of the electromagnetic field concentrate. According to [32], precisely wedges are often the sources and sinks of the vacancies that determine, for example, the hardness and plasticity of a solid body.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy and isolines of the magnetic field strength.
(a)
AmplitudeHy(b)
IsolinesAlso, we modeled the propagation of waves in media, when domains possess magnetic properties. We assumed, in the calculations, thatμ=100; the remaining parameters correspond to the previous example of solution (Figure 12).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyOf interest is the interaction of the electromagnetic wave with the rough surface shown in Figures13 and 14. As in the previous examples, we observe the concentration of electromagnetic energy on angular structures.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy  and isolines  of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesFrom Figures13 and 14, it is seen that electromagnetic heating of tapered structures may occur in addition to mechanical heating in magnetic abrasive machining.As we have mentioned above, for investigation of the propagation of electromagnetic waves in nonmagnetic materials, it is more expedient to use the generalized equation forE⃗. For the purpose of illustration we give an example of numerical calculation of an optical knife with the wave equation for H⃗ (Figure 15).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyFrom Figure15, it is seen that the actual problem of diffraction on the optical knife remains to be solved, that is, there is no “glow” on the optical-knife section, which is inconsistent with experimental data.As is known [33] in thermodynamically equilibrium systems the temperature T and the electrical φ and chemical μc potentials are constant along the entire system:(39)gradT=0,gradφ=0,gradμc=0.
If these conditions are not fulfilled (gradT≠0,gradφ≠0,gradμc≠0), irreversible processes of the transfer of mass, energy, electrical charge, and so forth appear in the system.The chemical potential of thejth component is determined, for example, as a change of the free energy with a change in the number of moles:(40)μcj=(∂F∂nj)T,V,
where
(41)dF=-SdT-PdV+H⇀dB⇀+EdD.
The last term in (41) takes into account the change in the free energy of a dielectric due to the change in the magnetic induction. The free energy of a unit volume of the dielectric in the magnetic field in this case has the form(42)F(T,D)=F0+μμ0H2⇀2+εε0E22.
We assume that changes in the temperature and volume of the dielectric are small. Then the mass flux is determined by a quantity proportional to the gradient of the chemical potential or, according to (41), we obtain(43)qi=-Dμcgrad(H⇀dB⇀+EdD)=-DμcgradW,
where W=μμ0(H2⇀/2)+εε0(E2/2) is the density of the magnetic field in the unit volume of the dielectric.In magnetic abrasive polishing on the sharp protrusions of domains the gradients of magnetic energy are great, which can lead to the origination of vacancy flows.An analysis of the results shows that the nonstationary component of the full electromagnetic energy is also concentrated in the region of fractures and wedges, that is, at the sharp angles of domains, which may lead to the improvement of the structure of the sublayer of the treated surface due to the “micromagnetoplastic" effect. Maximum values of the nonstationary part of the total electromagnetic energyWmax⁡ in the sublayer corresponds to a maximum value of the function sin (2π·106t)and occur for the time instants t=(n/4)10-6 sec, where n is the integer, with the value of Wmax⁡ for a neodymium magnet and a magnetoabrasive particle on the basis of carboxyl iron amounting to a value of the order of ((5)-(6))·106 J/m3. Having multiplied Wmax⁡ by the volume of a domain, vacancy, or atom, we may approximately obtain the corresponding energy. The density of the electromagnetic energy in all of the cases is much smaller than the bonding energy of atoms, 10−18-10−19 J. However, a periodic change in the magnetic field in one direction leads to a ponderomotive force that may influence the motion of various defects and dislocations to create a stable and equilibrium structure of atoms and molecules in magnetic abrasive polishing and, on the long run, in obtaining a surface with improved characteristics due to the “micromagnetoplastic" effect. The result of Para 2.2 were published in part [34].


## 2.2.1. Introduction

During the interaction of an external magnetic field and magnetic abrasive particles, the particles are magnetized, and magnetic dipoles with the moment oriented predominantly along the field are formed. “Chains" along the force lines of the field [22, 23] appear that periodically act on the processable surface with a frequency ω=l/υ. A fixed elemental area of the material periodically experiences the effect of the magnetic field of one direction. Actually, the frequency and duration of the pulse will be still higher because of the rotation of the magnetic abrasive particle due to the presence of the moment of forces on contact and of the friction of the particle against the processable part. In what follows, we will not take into account the effect of rotation.We assume that the particle velocity on the polisher isυ. If the particle radius is r, then the angular frequency is ω=2πυ/r, and precisely this frequency determines the frequency of the effect of the variable magnetic field component due to the fact that for a ferromagnetic μ>1. The magnetic permeability μ of ferromagnetics, which are usually used in magnetoabrasive polishing, is measured by thousands of units in weak fields. However, in polishing, the constant external magnetic field is strong and amounts to 105-106 A/m, and in this case the value of μ for compounds of iron and nickel and for Heusler alloy decreases substantially.Because of the presence of a strong external magnetic fieldH0 the “small" absolute value of μ of an abrasive particle leads to a periodic “increase" and “decrease" in the normal component of the magnetic induction near the processable surface. In the present work we used neodymium magnets (neodymium-iron-boron) with H0>485,000 A/m. The magnetic permeability of a magnetic abrasive particle based on carbonyl iron was assumed in this case to be equal to μ1=100.Due to the continuity of the normal magnetic induction componentBn1=Bn2, where Bn1=μ1μ0H1; Bn2=μ2μ0H2. For example, in glasses (μ2=1; therefore at the boundary of contact of the glass with the magnetic abrasive particle an additional variable magnetic field of strength H1>H0 appears.In [24–27], magnetic field-induced effects in silicon are considered: a nonmonotonic change in the crystal lattice parameters in the surface layer of silicon, the gettering of defects on the surface, the change in the sorption properties of the silicon surface, and the change in the mobility of the edge dislocations and in the microhardness of silicon.In [28–30], the influence of an electromagnetic field on the domain boundaries, plasticity, strengthening, and on the reduction of metals and alloys was established.In view of the foregoing, it is of interest to find the relationship between the discrete-impulse actions of a magnetic field of one direction on the surface layer of the processable material that contains domains. According to [22], the size of domains is as follows: 0.05 μm in iron, 1.5 μm in barium ferrite; 8 μm in the MnBi compound, and 0.5–1 μm in the acicular gamma ferric oxide. According to [31], the size of a domain may reach 10−6 cm3 (obtained by the method of magnetic metallography).As a rule, an abrasive exhibits a distinct shape anisotropy, whereas the frequency of the effect is determined by the concentration of abrasive particles in a hydrophobic solution and by the velocity of its motion. We assume that on the surface of a processable crystal the magnetic field strengthH(t)=H1sin⁡4(ωt)+H0.It is required to find the value of the magnetic field strength in the surface layer that has the characteristicsλ1, ε1, and μ1 and contains domains with electrophysical properties λ2, ε2, and μ2. The domains may have the form of a triangular prism, a bar, a cylinder, and so forth.


## 2.2.2. Physicomathematical Model: Wave Equation forH⇀

We consider interaction of a nonstationary magnetic field with substance on an example of a specific technology of magneto abrasive polishing based on use of magneto-rheological suspension. We will formulate a physicomathematical model of propagation of electromagnetic waves in a heterogeneous medium. The media in contact are considered homogeneous. We operate with the operator rot on the left- and right-hand sides of the first equation for the total current (6) and multiply by μ0μ; then we differentiate the second equation in (7) with respect to time. Taking into consideration the solenoidality of the magnetic field (7) and the rule of repeated application of the operator ∇ to the vector  H⇀, we obtain(35)μ0εε0∂2H⇀∂t2+λμ0∂H⇀∂t=1μ∇2H⇀.
In the Cartesian coordinates,(35) will have the form(36)εε0∂2Hx∂t2+λμ0∂Hx∂t=1μ(∂2Hx∂x2+∂2Hx∂y2+∂2Hx∂z2),εε0∂2Hy∂t2+λμ0∂Hy∂t=1μ(∂2Hy∂x2+∂2Hy∂y2+∂2Hy∂z2),εε0∂2Hz∂t2+λμ0∂Hz∂t=1μ(∂2Hz∂x2+∂2Hz∂y2+∂2Hz∂z2).
One fundamental electromagnetic field equation is the equation div⁡B⇀=0. The use of the Dirichlet theorem for approximation of the value of the magnetic field strength on the boundaries between adjacent media analogously to that of the electric field strength does not necessarily guarantees the observance of the condition of solenoidality of the magnetic field; furthermore, the magnetic properties of heterogeneous media were assumed constant in deriving generalized wave equations. The experience of numerical calculations has shown that when it is necessary to model nonstationary magnetic phenomena it is better in many cases to use a generalized wave equation for E⃗, accordingly expressing H⃗(t,r⃗) by E⃗(t,r⃗) and, if need be, to perform backward recalculation to H⃗(t,r⃗). This approach is difficult to apply to modeling of heterogeneous media with different magnetic properties, when the magnetic permeability μ is dependent on coordinates.In media with a weak heterogeneity whereμ(x,y,z) is a piecewise continuous quantity, the application of the proposed method of through counting is quite justified. Indeed, the system of (13) and (36) yields that the function’s discontinuity on the boundaries between adjacent media is determined by the complexes which will be called the generalized permeability ε*≡εμε0μ0 and the generalized conductivity λ*≡λμμ0. Using the Direchlet theorem for ε* and λ*, we obtain their values on the boundaries between adjacent media and the values for the electric field strength at the discontinuity point (see (34)); here, we note that the value of the electric field strength is obtained without solving Maxwell equations. In fact, at the discontinuity point, we use linear interpolation of the function to obtain the values of ε*, λ*, and Ex=ξ=(1/2)[E(ξ-0)+E(ξ+0)]. Consequently, for piecewise continuous quantity μ(x,y,z), the application of the proposed method of through counting is justified. We note that the equality of the derivatives of the electric field strength along the normal to the surface at the discontinuity point according to Equation (32) holds. When the wave equation for H⃗ is used for media with different magnetic permeabilities the condition of equality of the derivatives fails, that is,(37)∂Hx∂x|x=ξ-0≠∂Hx∂x|x=ξ+0,
which is a consequence of (10); therefore, the use of through-counting schemes for the wave equation for H⃗ is difficult.Moreover, reformulation of the boundary conditions for equations (6) and (7), with electric field strength being eliminated from the boundary conditions, does not even allow prediction of the surface current direction. Indeed, magnetic field strength H⇀ is an axial vector which on substitution r⃗→(-r⃗)  must remain constant. The value of the surface current also cannot be determined without knowledge of the polar vector of electric field strength. The generalized wave equation for E⃗ contains the term graddiv⁡E⃗ which directly allows for the influence of induced surface charges on the propagation of waves,  the right-hand side of the generalized wave equation has the same form as the equations of the linear elasticity theory, hydrodynamics. We note that the proposed method of calculation can be used on condition that there are no built-in space charges and extraneous electromotive forces [13].By virtue of what has been stated above, for modeling of the propagation of electromagnetic waves in glasses having roughness and defects, we used system (13) with boundary conditions (22)–(31) and (32).


## 2.2.3. Results of Numerical Simulation

The physicomathematical model developed can also efficiently be used in modeling the propagation of electromagnetic waves in media with complex geometries and strong electromagnetic field discontinuities.The transverse cut of a cellular structure represents a set of parallelepipeds and triangular prisms of various cross-sections, as depicted in Figures 10(a) and11(a). An electromagnetic wave propagates across the direction of parallelepipeds and triangular prisms (channels) along the coordinate x.The size of the investigated two-dimensional object is14×20·10-6 m, and the sizes of the domains are 2–4 μm. The frequency of the influence of the magnetic field is ω=2π·106, and the strength of the field is(38)Hx=21⋅105sin4(2π⋅106t)Am.
The electrophysical properties are as follows: of the large parallelepiped, μ=1, ε=8, σ=10-9 Ω·m; of domains, μ=1, ε=6, σ=10-8 Ω·m. They correspond to the electrophysical properties of glasses.It was assumed that in a layer of thickness 15–20μm an electromagnetic wave propagates without attenuation; therefore, on all the faces of the large parallelepiped the fulfillment of condition (38) was considered valid. On the faces of the parallelepiped that are parallel to the OX axis condition (38) corresponded to the “transverse" tangential component of the wave; on the faces parallel to OY condition (38) corresponded to the normal component of the field.The calculations were carried out with a time step of 10−13 sec up to a time instant of 10−10 sec.Figures10(a) and 11(a) present the amplitude values of the magnetic field strength along Hx and Hy with a comparison scale, whereas Figures 11(b) and 11(b) present the corresponding isolines. An analysis of these figures shows that at the places of discontinuity, on the wedges, force lines of the electromagnetic field concentrate. According to [32], precisely wedges are often the sources and sinks of the vacancies that determine, for example, the hardness and plasticity of a solid body.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy and isolines of the magnetic field strength.
(a)
AmplitudeHy(b)
IsolinesAlso, we modeled the propagation of waves in media, when domains possess magnetic properties. We assumed, in the calculations, thatμ=100; the remaining parameters correspond to the previous example of solution (Figure 12).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyOf interest is the interaction of the electromagnetic wave with the rough surface shown in Figures13 and 14. As in the previous examples, we observe the concentration of electromagnetic energy on angular structures.AmplitudeHx and isolines of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesAmplitudeHy  and isolines  of the magnetic field strength.
(a)
AmplitudeHx(b)
IsolinesFrom Figures13 and 14, it is seen that electromagnetic heating of tapered structures may occur in addition to mechanical heating in magnetic abrasive machining.As we have mentioned above, for investigation of the propagation of electromagnetic waves in nonmagnetic materials, it is more expedient to use the generalized equation forE⃗. For the purpose of illustration we give an example of numerical calculation of an optical knife with the wave equation for H⃗ (Figure 15).AmplitudeHx  and Hy of the magnetic field strength.
(a)
AmplitudeHx(b)
AmplitudeHyFrom Figure15, it is seen that the actual problem of diffraction on the optical knife remains to be solved, that is, there is no “glow” on the optical-knife section, which is inconsistent with experimental data.As is known [33] in thermodynamically equilibrium systems the temperature T and the electrical φ and chemical μc potentials are constant along the entire system:(39)gradT=0,gradφ=0,gradμc=0.
If these conditions are not fulfilled (gradT≠0,gradφ≠0,gradμc≠0), irreversible processes of the transfer of mass, energy, electrical charge, and so forth appear in the system.The chemical potential of thejth component is determined, for example, as a change of the free energy with a change in the number of moles:(40)μcj=(∂F∂nj)T,V,
where
(41)dF=-SdT-PdV+H⇀dB⇀+EdD.
The last term in (41) takes into account the change in the free energy of a dielectric due to the change in the magnetic induction. The free energy of a unit volume of the dielectric in the magnetic field in this case has the form(42)F(T,D)=F0+μμ0H2⇀2+εε0E22.
We assume that changes in the temperature and volume of the dielectric are small. Then the mass flux is determined by a quantity proportional to the gradient of the chemical potential or, according to (41), we obtain(43)qi=-Dμcgrad(H⇀dB⇀+EdD)=-DμcgradW,
where W=μμ0(H2⇀/2)+εε0(E2/2) is the density of the magnetic field in the unit volume of the dielectric.In magnetic abrasive polishing on the sharp protrusions of domains the gradients of magnetic energy are great, which can lead to the origination of vacancy flows.An analysis of the results shows that the nonstationary component of the full electromagnetic energy is also concentrated in the region of fractures and wedges, that is, at the sharp angles of domains, which may lead to the improvement of the structure of the sublayer of the treated surface due to the “micromagnetoplastic" effect. Maximum values of the nonstationary part of the total electromagnetic energyWmax⁡ in the sublayer corresponds to a maximum value of the function sin (2π·106t)and occur for the time instants t=(n/4)10-6 sec, where n is the integer, with the value of Wmax⁡ for a neodymium magnet and a magnetoabrasive particle on the basis of carboxyl iron amounting to a value of the order of ((5)-(6))·106 J/m3. Having multiplied Wmax⁡ by the volume of a domain, vacancy, or atom, we may approximately obtain the corresponding energy. The density of the electromagnetic energy in all of the cases is much smaller than the bonding energy of atoms, 10−18-10−19 J. However, a periodic change in the magnetic field in one direction leads to a ponderomotive force that may influence the motion of various defects and dislocations to create a stable and equilibrium structure of atoms and molecules in magnetic abrasive polishing and, on the long run, in obtaining a surface with improved characteristics due to the “micromagnetoplastic" effect. The result of Para 2.2 were published in part [34].


## 3. Ion Conductors: New Closing Relations on the Boundaries of Adjacent Media

### 3.1. Diffusion Electric Phenomena in Electrolytes

The medium under discussion is assumed to consist of a non ionized solvent, an electrolyte in the form of ions and uncharged components. At first we consider infinitely diluted electrolytes when molecules are completely dissociated into cations and anions.In a solution, ions are transferred by convection, diffusion, and migration in the presence of an electric field. Derivation of the equation for ion migration is based on the following considerations. In a solution with molar concentrationn and diffusion coefficient Di let there be ions with charge zi.When an ion is exposed to an electric field with intensity E applied to the solution, the ion experiences the force zieE, which brings it into motion. The ion velocity is related to this force by the usual expression known for the motion of particles in a viscous medium:(44)u=γi*zieE,
where γi* is the ion mobility. The latter may be expressed in terms of the diffusion coefficient using the known Einstein relation:(45)γi*=DiFRT.
The total flux of ions of the i th kind in a moving medium in the presence of diffusion and migration is determined by the Nernst-Planck equation:(46)qi=niυ-Di∇ni+DiziFERTni.
Formulas (44) and (45) have, in fact, a limited sphere of applicability. Indeed, Einstein's work [35] is concerned only with the diffusion of a neutral impurity with its small concentration in a solution when the usual relations of hydrodynamics are valid for a flow around a sphere.In the physics of plasma, formulas (44) and (45) are based on other considerations and provided the plasma is weakly ionized, that is, particles move independently of each other. Here, only the collision of charged particles with neutral ones is taken into account [36].In [37, 38], for describing the diffusion and migration of ions in a partially dissociated electrolyte it is suggested to take into consideration their transfer by neutral molecules. In the theory developed by Y. Kharkats [37] the expression for the total flux of ions is(47)qi=niυ-Di∇ni+DiziFERTni-Dzi∇nA.
The last term in (47) takes into account ion transfer by a flux of neutral molecules; here the concentrations of anions, cations, and neutral dissociated molecules are determined by the chemical equilibrium conditions(48)β=(n-)z+(n+)z-nA,
where β is the dissociation equilibrium constant. If the degree of dissociation at the prescribed total concentration n is α, then nA=(1-α)n and, consequently,(49)β=α2n1-α.
The dissociation constant β, unlike the degree of dissociation, must not depend on concentration. However, in real partially dissociated electrolytes those arguments for α and β are in rather poor agreement with experiment. In [39] the dissociation degree for HCl is listed in Table 1 as a function of concentrations, calculated in terms of measured electric conductivity α1 and e.m.f. α2.Table 1
The degree of HCl electrolytical dissociation in terms of measured electric conductivityα1 and E.m.f. α2.
CHCl, mole/liter0.0030.080.33.06.016.0α10.9860.9570.903———α20.990.880.7731.4023.413.2Inspection of the table reveals that the dissociation degrees obtained by different experimental methods coincide best in the case of dilute solutions. In the high-concentration range of the electrolyte,α2 even exceeds unity, which, naturally, has no physical sense.According to the Arrhenius theory the dissociation constantβ for the given electrolyte at the prescribed temperature and pressure must remain constant independently of the solution concentration. In [40], the dissociation constants of some electrolytes are given at their different concentrations. Only for very weak electrolytes (solutions of ammonia and acetic acid) does the dissociation constant remain more or less constant on dilution. For strong electrolytes (potassium chloride and magnesium sulfate), it changes by severalfold and in no way may be considered a constant.Of course, one may formally take into account the dependence of the dissociation constant on the electrolyte concentration and use it in the modified Nernst-Planck equation (47) and (48) but the main drawback in describing separately the diffusion and migration of ions in an electrolyte, in D. I. Mendeleev's opinion, lies in the fact that the interaction of particles of a dissolved substance between each other as well as with solvent molecules is neglected. Also, he has pointed out that not just the processes of formation of new compounds with solvent molecules are of importance for solutions. D. I. Mendeleev's viewpoints have been extended by Sakhanov [40],who believed that in addition to the usual dissociation reaction in an electrolyte solution there also proceeds associating of simple molecules. Molecular associations dissociate, in their turn, into complex and simple ions. In this case (47) will not hold mainly because of the fact that the expressions for diffusion, migration, and nondissociated flows of molecules must be determined relative to some mean liquid velocity. In concentrated solutions, this velocity does not coincide with the solvent velocity and must be determined from the fluid dynamics equations of a multicomponent mixture in which characteristics of its components (physical density, charge, diffusion coefficient of a complex ion) are, in fact, unknown. This is a reason why the theory of diffusion and migration of ions for a partially dissociated electrolyte solution, with a current traversing through it, encounters crucial difficulties. Besides, it is rather difficult to take an account of the force interaction of complex cations and anions between each other and with an external electric field. Therefore Einstein's formula (44) will also change its form.The main drawbacks of the theory of electrolytic dissociation are fully defined in collected papers “Fundamental Principles of Chemistry" by D. I. Mendeleev as well as in [40, 41].In our opinion, the principal disadvantage of the approaches considered above lies in the fact that the Nernst-Planck equation as well as its modified form for a partially dissociated electrolyte (see (47)) are based on the hydrodynamic theory of diffusion. We shall substantiate below that it is more reasonable to use the equations of thermodynamics of irreversible processes.When applying the Nernst-Planck equation (see (46)) and its modified forms (see (47)), it is implicitly postulated that cations and anions in the electrolyte solution are different components. The equations of transfer of cations and anions are written separately. On the other hand, the commonly used equation of electroneutrality does not ensure, according to Gibbs [39], the independence of differentials of cations ∂ni and anions ∂ne in a volumetric electrolyte solution their changes are functionally related and therefore they are identical, indiscernible components.


#### 3.1.1. The Case of a Current Traversing a Binary Fully Dissociated Electrolyte

Let us write system of equations (see (46)) in the form:(50)∂ne∂τ+vgradne=DeΔne+zeFDeRTdiv⁡(neE),∂ni∂τ+vgradni=DiΔni+ziFDiRTdiv⁡(niE),(51)zene-zini=0.
Analogously to [39], instead of concentrations ne and ni we introduce the molar concentration related to ne and ni by the expressions(52)C=nezi=nize.
Expressing ne and ni in (50) in terms of the molar concentration, we arrive at the following expression for the function C:(53)∂C∂τ+υ∇C=DefΔC,
where(54)Def=DeDi(ze+zi)zeDe+ziDi
is the diffusion coefficient of a salt or the effective diffusion coefficient of a binary electrolyte. Expressing the concentrations ne and ni in terms of the molar concentration C we find the expression for the vector of current density J [42]:(55)J=(Di-De)Fzize∇C+F2zizeRT(ziDi+zeDe)CE.
In [43] it is shown that by its meaning Def is a coefficient of molecular diffusion. For instance, the diffusion coefficients of copper sulfate have an intermediate value between the diffusion coefficients of copper and sulfate ions; if De=0.173·10-9 m2/sec and Di=1.065·10-9 m2/sec, then Def=0.854·10-9 m2/sec. However, expression (54) does not offer an explanation of the abnormally high mobility of H and OH ions when their displacement proceeds by a croquet-like mechanism [44].In the system of (53)–(55) no symmetry of cross terms is observed; therefore, it is difficult to take into account “the superposition effects." Indeed, from (55) it is evident that the molar concentration gradients VC exert an influence on the current J. On the other hand, one can see from (53) that an electric field has no effect on the diffusion of the molecules, and therefore the Onsager reciprocal relation is disturbed.One may raise the objection that according to the Nernst-Planck equations it follows that the concentration gradients of cations and anions exert an influence on the potential distribution of an electric field [42]: (56)∇(ρR∇φ)=-F∑iz∇(Di∇Ci).
However, relation (56) itself has been derived by Neumann [43] on the assumption that the resultant current does not depend on the mass flow.On derivation of the equations for diffusion and migration of ions in a volumetric electrolyte solution we shall employ the relations of the thermodynamics of irreversible processes.Finally, it is worthy of note that in plasma physics there is also no symmetry of “the superposition effects" in describing the ambipolar diffusion of ions, and the theory is based on the hydrodynamic model [36].In investigating the dynamics of an electron-hole semiconductor plasma, the hydrodynamic model (50)-(51) is also used.In our opinion, electric field-induced transfer of charges entails transfer of their kinetic energy as well as of heat and mass, and, conversely, mass or heat transfer may simultaneously cause, if we are concerned with the system of charged particles, charge transfer and give rise to an electromotive force.


#### 3.1.2. Conclusion

To phenomenologically describe diffusional-electrical phenomena in electrolytes it is suggested to use the methods of the thermodynamics of irreversible processes, rather than the hydrodynamic theory of diffusion of ions, for the experimentally observable thermodynamic flows and forces without clear discernment and identification ofrealion flows and mobilities.


### 3.2. Landau Model

In what follows we consider ion conductors that contain the electrolytic plasma with a fairly high concentration of charged particles. The interaction of the charged plasma particles with the external electromagnetic field and their collective interaction can theoretically be investigated completely only on the basis of the Boltzmann kinetic equation with a self-consistent field. In this probabilistic approach to the phenomena under study, one introduces characteristics that are average over a large ensemble of particles and are always related to the introduction of additional hypotheses on the properties of particles and their interactions and to the simplification of these properties and interactions. We note that in many cases, for example, when the electrolytic plasma in liquid electrolytes is considered, there is even no basis for constructing such methods. The developed methods are usually not an efficient means of solving problems by virtue of the excessive complexity of corresponding equations.Another general method is construction of a phenomenological macroscopic theory based on the general regularities and hypotheses obtained experimentally; this precisely method will be developed in the present work. For many problems of applied character, consideration of a plasma as a conducting gas is sufficient. Such approximation is strictly justified only in the case of a dense plasma when the free path of charged particles is much smaller than the characteristic dimension of the system and particle collisions are of crucial importance. The velocity distribution of particles is Maxwellian; at each point, it is fully determined by the local values of density, temperature, and macroscopic velocity. With these conditions being observed, we can carry out, in a unified context, macroscopic phenomenological description of a gas medium, a metal and dielectric plasma, a plasma in a liquid electrolyte, or an electrolytic plasma.The necessary condition of existence of the state of local thermodynamic equilibrium (LTE) of a plasma is high frequency of collisions of plasma particles (Maxwellization) so that the plasma state undergoes no marked changes over the period between collisions on the mean free path. The observance of this condition means the following: (1) the electrons have time to transfer the energy received from the electrostatic field to heavy particles; (2) ionization processes are nearly totally counterbalanced by recombination (the Saha equation holds true); (3) the overwhelming part of excited atoms gives up its energy in collisions; (4) energy exchange between particles dominates the processes in which the plasma energy markedly grows or diminishes.In the LTE plasma state, it is only the radiation field (photons) that is not in equilibrium with the plasma, but the energy loss by radiation is minor in a number of cases compared to other mechanisms of energy transfer. With allowance for the assumptions enumerated above, the plasma electrodynamics can be described by Maxwell continuity, motion, and energy equations [45]. For macroscopic description of transfer processes in ion conductors, we will use the equations of thermodynamics of irreversible processes without separating explicitly ion carriers.The possibility of such description was first suggested by L. D. Landau and E. M. Lifshits. This approach which has been substantiated above does not provide for the introduction and determination of the coefficients of molecular diffusion of cations and anions and of the degree of dissociation in an ion conductor. It is conditioned by the interaction of the mass and charge fluxes. The ion fluxes are not determined. The density of the conduction currentJq, the mass flux Jm, and the heat flux JT for the medium’s nonequilibrium state has the form [45]:(57)Jq=λ(E-βgradn)-λα(T)gradT,Jm=-Dmgradn-DA*λE-DTgradT,JT=-kgradT-Iq(Π+φ)-λDT*gradn,
where DM is the coefficient of molecular diffusion, DA* is the coefficient of ambipolar diffusion, DT* is the coefficient allowing for heat transfer by the impurity motion, β is the coefficient of specific electrical ambipolar conductivity due to ∇n, and Π is the Peltier coefficient.To investigate the interaction of electromagnetic impotential fields (rot⁡E⃗≠0) and heat and mass transfer rates it is advisable to use the system of (57) for polar vectors; if needed, the pseudovector of the magnetic field strength can be calculated from the Maxwell equation (7).


### 3.3. Interaction of Nonstationary Electric, Thermal, and Diffusion Fields with Allowance for Mass Transfer in a Layered Medium with the Example of an Electrochemical Cell

#### 3.3.1. Introduction

Investigation of the interaction of electric and thermal fields with allowance for mass transfer and contact phenomena is a complex and topical problem of the theory and practice of various fields of natural science and technology.The work seeks to construct a physicomathematical model of the interaction of nonstationary electric fields in a layered medium with allowance for nonstationary thermal phenomena and mass transfer without explicit separation of charge carriers. The media in contact are considered to be homogeneous. For the sake of clear representation, we will consider two-layer one-dimensional models.


#### 3.3.2. Interaction of Electric and Thermal Fields

In different substances, the processes of transfer of charge and energy are interrelated. The quantity of the released heat is determined by Joule heating and by the effects of Thomson and Peltier. The problem of interaction of nonstationary thermal phenomena has been considered in [46] without taking into account the Thomson effect.According to [45], the expressions for the conduction-current density and the energy-flux density in the absence of an external magnetic field or in the case of its slight influence have, respectively, the form
(58)Iq=-λ(α(T)gradT+gradφ),IT=-k(T)gradT+Iq(Π+φ),
where a(T) is the specific thermoelectromotive force, Π=αT is the Peltier coefficient, φ is the potential, and k(T) is the thermal-conductivity coefficient. We note that the problem is always nonlinear.Having eliminated the magnetic-field strength from system (13) according to one method of [47–49], we obtain the equation for the electric-field strength(59)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E.
The heat-balance equation has the form [45] (60)cpρ∂T∂t=div⁡(k(T)gradT)+Iq⋅(E-(α(T)+Tβ(T))gradT).
At the interface, the following relation [50] holds:(61)div⁡i+Iqn1-Iqn2=-∂σ∂t.
Here ε is the permittivity, c2=1/(ε0·μ0), μ is the magnetic permeability, ε0 is the electric constant, μ0 is the magnetic constant, cpρ is the product of the specific heat of the medium by its density, and β(T)=∂α(T)/∂T.In deriving (60), use has been made of the condition of local electroneutrality of the substance. We give a differential formulation of the problem and the method of its solution in greater detail.Let us investigate, on the segment[0≤x<l]=[0≤x<ξ]∪[ξ≤x≤l],the contact of homogeneous media 1 and 2 with different electrophysical properties. The quantities ε,λ,μ,E, and φ have discontinuities of the first kind at the point of the interface x=ξ. We will consider the case of plane contact where the influence of surface currents can be disregarded and the thickness of the electric double-layer is much smaller than the characteristic dimension of the object. We set E=-gradφ. Equations (59) and (60) for the one-dimensional problem will take the form(62)εc2∂2E∂t2+μ0∂∂t(-λ(α(T)∂T∂x-E))=1μ∂2E∂x2,cpρ∂T∂t=∂∂x(k(T)∂T∂x)-λ(α(T)∂T∂x-E)(E-(α(T)+Tβ(T))∂T∂x).
Conditions (39) and (61) will be written, respectively, in the form(63)ε1ε0E1|x=ξ-0-ε2ε0E2|x=ξ+0=σ|x=ξ,(64)λ1(E1-α1(T)∂T∂x)|x=ξ-0-λ2(E2-α2(T)∂T∂x)|x=ξ+0=-∂σ∂t|x=ξ,
By differentiating (63) with respect to time and taking into account (64), at the phase boundary we obtain the condition of equality of the total currents:(65)[λ(E-α(T)∂T∂x)+εε0∂E∂t]|x=ξ=0.
Here and in what follows, for the arbitrary function f we adopt the notation
(66)[f]|x=ξ=f|x=ξ+0-f|x=ξ-0.
Equality (63) is a corollary of the relation lim⁡x→ε±0div⁡D=ρ [50]. Taking into consideration the finiteness of the value of the space electric charge and the continuity of it, in the one-dimensional case we obtain(67)[εε0∂E∂x]|x=ξ=0.
Relations (65) and (67) reflect the laws of conservation and continuity of electric charge [45, 50]. At the contact point, we also have the equality of the temperatures and the energy fluxes [45]:(68)[T]|x=ξ=0,(69)[k(T)∂T∂x+λ(α(T)∂T∂x-E)(α(T)T+φ)]|x=ξ=0.
Thus, in the presence of the interaction of the electric and thermal fields in a layered medium, the equality of the charge fluxes (65), the equality of the charges (67) (when the conditions of quasineutrality of the contacting media beyond the electric double-layer are satisfied), the equality of the temperatures (68), and the equality of the energy fluxes must be fulfilled at the interface of the media. In the relations for the charge and energy fluxes, we take into account cross-thermoelectrical phenomena.The process of charge of the electric field is considered on the finite time interval [0≤t≤t0]. The initial conditions have the form(70)T(x,0)=T0(x),E(x,0)=f1(x),∂E(x,0)∂t=f2(x).
We give the boundary conditions, for example, for the value x=0. We set the value of the total-current density to be known(71)(λ(α(T)∂T∂x-E)-εε0∂E∂t)|x=0=j1(t),
and take into account the heat exchange at the boundary using the Newton relation(72)(k(T)∂T∂x+λ(α(T)∂T∂x-E)α(T)T+φ)|x=0=γ1(T-T*)|x=0,
where γ1 is the coefficient of heat exchange and T* is the ambient temperature.We consider the analogous relations at the right-hand boundary.


#### 3.3.3. Interaction of Electric and Thermal Fields with Allowance for Mass Transfer

Let us consider the nonstationary model of heat and mass transfer for electrochemical systems with the example of electrolysis. For this problem, an electric double-layer occurs at the site of contact of the electrolyte with the metallic anode and cathode. Let us set the density of the space charge beyond the electric double-layer of the metal-electrolyte contact to be zero at the initial instant and hence, according to [51], remain constant in the future while the voltage drop in the electrodes and leads is small. The influence of the electrodes on the temperature field of the electrolyte will be taken into account in terms of the coefficient of heat exchange.Two approaches to modeling of diffusion-electrical phenomena have currently been developed, each of which has certain disadvantages and advantages. The first approach [37, 38, 42] is characterized by consideration of the flows of ions of the corresponding sort in a completely or partially dissociated electrolyte. The equations derived contain many parameters, the reliable procedure of determination of which is absent in the majority of cases. Furthermore, the proposed theory fails to provide for taking into account the interaction of the cations and the anions of the dissolved substance with each other and with the molecules of the solvent. It must also be borne in mind that, apart from the dissociation reaction, we have the association of simple molecules which in turn are dissociated into complex and simple ions. It is difficult to determine the diffusion coefficients of complex ions; it is also difficult to take into account the force interaction of complex cations and anions with each other and with the external electric field. The difficulties in question become more serious when a multicomponent electrolyte is described. We note that such difficulties also arise in modeling, for example, certain problems of plasma physics.The second approach [52–56] fails to provide for the introduction and determination of the coefficients of molecular diffusion of cations and anions and the degree of dissociation of the electrolyte. It is based on the interaction of the mass and charge fluxes. The ionic flows are not determined. The model is formulated without explicit separation of the charge carriers and the charge of an electric double layer.The system of equations describing the electrodynamic processes in the electrolyte is as follows:(73)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E,∂n∂t=-div⁡(Im),c0ρ0∂T∂t=-div⁡(IT).
Here Iq, Im, and IT are the densities of the conduction current and of the mass and heat fluxes which (in the case of the nonequilibrium state of the medium) have, respectively, the form (57) [45], where DM is the molecular diffusion coefficient, DA* is the ambipolar diffusion coefficient, DT is the thermal diffusion coefficient, DT* is the coefficient taking into account the transfer of heat due to the motion of the impurity, and n is the concentration of the impurity.We give the initial(74)E(x,0)=E0(x),∂E(x,0)∂t=Ẽ(x),T(x,0)=T0(x),n(x,0)=n0(x),
and boundary conditions (anode-electrolyte, x=0):(75)Iq=λ(λA*∂n∂x-E)+λα(T)∂T∂x-εε0∂E∂t,γ1(T-T*)=k∂T∂x-Iq(Π+φ)+DT*∂n∂x,kaIq(t)=DM∂n∂x-DA*λE+DT∂T∂x,
in modeling the one-dimensional problem in the region G=[0≤x≤l]×[0≤t≤T].Here ka is the electrochemical equivalent of a substance which deposits at the anode. The analogous conditions hold for x=l as well.The constructed system of differential equations makes it possible to model the processes of transfer in the electrochemical system with allowance for the influence of an electric double-layer; this influence directly determines the thermoelectrical and ambipolar electrodiffusion phenomena at the boundary of the metal-electrolyte contact.


#### 3.3.4. Numerical Modeling of Electrochemical Systems. Heating of an Electrochemical Cell in the Case of Constant-Current and Pulse Electrolysis

In electrochemical systems, the electrodes are metals as a rule. The electrical conductivity of metals is hundreds of times higher than the specific electrical conductivity of electrolytes; therefore, one can disregard the voltage drop in electrodes and leads. In highly conductive metals, one can also disregard the displacement currents and the coefficients of ambipolar diffusion and thermal diffusion and restrict oneself to consideration of the system of equations (73)–(75) just in the electrolyte itself.Let us consider the process of heating of an electrochemical medium in the case of passage of an electric current with densityIq(t) through a solution of copper sulfate CuSO4·5H2O. We use copper (99.78%) as the anode. The copper from the solution is deposited at the cathode. We set the current efficiencies for the copper equal to 100% and the electrochemical equivalent ke equal to 0.6588·10−6 kg/C. The dependence of the specific electrical conductivity λ(n) of the copper-plating electrolyte on the concentration of the copper sulfate in water is given in [57]; the dependence of λ and k on the temperature will be considered to be insignificant. In the calculations for CuSO4, we set DM=5·10-10 m2/sec, DA*=10-11λ,λA*=10-4λ,ε=70, and μ=1 and disregard the effects of thermal diffusion. The distance between the cathode and the anode is L=0.05 m. The heat capacity, the density, and the thermal-conductivity coefficient of the electrolyte are taken to be 4.2·103 J/(kg·K), 103 kg/m3, and 0.6 V/(m·K) respectively. We note that the regimes of pulse action have been considered in [58] as applied to the processes of electrodeposition of alloys. This investigation has been carried out on the basis of a separate description of the transport of ions and the employment of the Kirchhoff law for quasistationary currents in a cell. With such an approach, the displacement currents are taken into account indirectly through introduction of the capacitive current of the electrode and determination of the experimental dependence of the polarization capacitance of the electrode on the character of the pulse action. In [58], it has been noted that neglect of the capacitive current in pulse electrolysis involves significant errors. The employment of system (73)–(75) to model the process of copper plating makes it possible to consider the problems without introducing the concepts of inductance and capacitance. The obtained nonlinear system of equations is solved by the finite-difference method analogously to [54].We have modeled numerically the process of copper plating for a constant current and a pulse current (Figure16). The density of the constant current was equal to 300 A/m2. The maximum density of the pulse current was also 300 A/m2, while the period of traversal of the current and the break were equal to 0.01 sec.Distribution of the normalized concentration of CuSO4  and the temperature in an electrochemical cell within 60 sec after switching on the current. The solid curve, at constant current; the dashed curve, at pulse current. T,°C; x,mm.
(a)
Distribution of the normalized concentration of CuSO4(b)
Distribution of the temperatureFigure16 gives the results of modeling the distributions of the concentration of CuSo4 and the temperature for different regimes of electrolysis. The measurements were carried out within 60 sec after switching on the current. It is seen from the figure that the concentration gradients near the cathode surface (x=1) are different and depend on the regime of copper plating. In the case of nonstationary action they are substantially smaller than in the case of electrolysis in the regime of constant current. As has been shown in [58], this tendency also holds in the case where one and the same total charge passes through the electrochemical cell.The analysis of the temperature field (Figure16(b)) shows that the temperatures of the electrolyte in the vicinity of the anode and the cathode differ significantly. Nonsymmetric heating of the electrochemical cell attributed to the am-bipolar diffusion and electrical ambipolar conductivity occurs. The modeling results are in agreement with the data of [59].With the aim of measuring the temperature in the surface layer of the electrolyte solution near the cathode and comparing the results to the calculated data, we developed and manufactured an experimental setup. Copper anode and cathode were placed in an electroplating bath with a sulfuric-acid copper-plating electrolyte. In the 1 mm thick cathode, a hole was made into which we placed a temperature-sensitive element. As the latter we employed a thermal resistor. The current strength in the bath was prescribed by a variable resistor.Figure17 gives the results of numerical modeling and the experimental data on heating of the cathode region in relation to the regime of current. It is seen that the increase in the density of the current (both constant and pulse) causes the temperature near the cathode surface to increase. In the anode, no increase in the temperature was observed either in the calculations or in the experiment, which is in agreement with the data of [59]. The surface temperature obtained in solving (73)–(75) in the regime of constant current is 10–15% higher than the temperature recorded in the experiments.Figure 17
Heating of the electrolyte in the cathode region as a function of the value and kind of polarizing current: (1) constant current; (2) pulse current; (3, 4) corresponding experimental data.ΔT,°C; Iq, A/m2.A more complicated situation occurs when we compare the experimental data and the results of numerical modeling in the regime of pulse electrolysis. The results of the experiments demonstrate that in the regime of pulse current the temperature of the electrolyte in the cathode region increases more than in the regime of constant current. At the same time, the numerical calculations of the pulse electrolysis show that the heating of the cathode region here is smaller than for the regime of constant current. The difference is, apparently, attributed to concentration convection [60], whose influence is substantial for high densities of the current for pulse electrolysis.In a longer term, one must take into account natural and concentration convection and heat exchange with the ambient medium. Nonetheless, despite the errors in determination of the coefficients and the assumptions made in the model, the calculation results are in qualitative agreement with the experimental data, which makes it possible, in certain cases, to use the proposed approach for modeling of nonstationary processes.Results were published in part [52, 53, 55, 56, 61].


## 3.1. Diffusion Electric Phenomena in Electrolytes

The medium under discussion is assumed to consist of a non ionized solvent, an electrolyte in the form of ions and uncharged components. At first we consider infinitely diluted electrolytes when molecules are completely dissociated into cations and anions.In a solution, ions are transferred by convection, diffusion, and migration in the presence of an electric field. Derivation of the equation for ion migration is based on the following considerations. In a solution with molar concentrationn and diffusion coefficient Di let there be ions with charge zi.When an ion is exposed to an electric field with intensity E applied to the solution, the ion experiences the force zieE, which brings it into motion. The ion velocity is related to this force by the usual expression known for the motion of particles in a viscous medium:(44)u=γi*zieE,
where γi* is the ion mobility. The latter may be expressed in terms of the diffusion coefficient using the known Einstein relation:(45)γi*=DiFRT.
The total flux of ions of the i th kind in a moving medium in the presence of diffusion and migration is determined by the Nernst-Planck equation:(46)qi=niυ-Di∇ni+DiziFERTni.
Formulas (44) and (45) have, in fact, a limited sphere of applicability. Indeed, Einstein's work [35] is concerned only with the diffusion of a neutral impurity with its small concentration in a solution when the usual relations of hydrodynamics are valid for a flow around a sphere.In the physics of plasma, formulas (44) and (45) are based on other considerations and provided the plasma is weakly ionized, that is, particles move independently of each other. Here, only the collision of charged particles with neutral ones is taken into account [36].In [37, 38], for describing the diffusion and migration of ions in a partially dissociated electrolyte it is suggested to take into consideration their transfer by neutral molecules. In the theory developed by Y. Kharkats [37] the expression for the total flux of ions is(47)qi=niυ-Di∇ni+DiziFERTni-Dzi∇nA.
The last term in (47) takes into account ion transfer by a flux of neutral molecules; here the concentrations of anions, cations, and neutral dissociated molecules are determined by the chemical equilibrium conditions(48)β=(n-)z+(n+)z-nA,
where β is the dissociation equilibrium constant. If the degree of dissociation at the prescribed total concentration n is α, then nA=(1-α)n and, consequently,(49)β=α2n1-α.
The dissociation constant β, unlike the degree of dissociation, must not depend on concentration. However, in real partially dissociated electrolytes those arguments for α and β are in rather poor agreement with experiment. In [39] the dissociation degree for HCl is listed in Table 1 as a function of concentrations, calculated in terms of measured electric conductivity α1 and e.m.f. α2.Table 1
The degree of HCl electrolytical dissociation in terms of measured electric conductivityα1 and E.m.f. α2.
CHCl, mole/liter0.0030.080.33.06.016.0α10.9860.9570.903———α20.990.880.7731.4023.413.2Inspection of the table reveals that the dissociation degrees obtained by different experimental methods coincide best in the case of dilute solutions. In the high-concentration range of the electrolyte,α2 even exceeds unity, which, naturally, has no physical sense.According to the Arrhenius theory the dissociation constantβ for the given electrolyte at the prescribed temperature and pressure must remain constant independently of the solution concentration. In [40], the dissociation constants of some electrolytes are given at their different concentrations. Only for very weak electrolytes (solutions of ammonia and acetic acid) does the dissociation constant remain more or less constant on dilution. For strong electrolytes (potassium chloride and magnesium sulfate), it changes by severalfold and in no way may be considered a constant.Of course, one may formally take into account the dependence of the dissociation constant on the electrolyte concentration and use it in the modified Nernst-Planck equation (47) and (48) but the main drawback in describing separately the diffusion and migration of ions in an electrolyte, in D. I. Mendeleev's opinion, lies in the fact that the interaction of particles of a dissolved substance between each other as well as with solvent molecules is neglected. Also, he has pointed out that not just the processes of formation of new compounds with solvent molecules are of importance for solutions. D. I. Mendeleev's viewpoints have been extended by Sakhanov [40],who believed that in addition to the usual dissociation reaction in an electrolyte solution there also proceeds associating of simple molecules. Molecular associations dissociate, in their turn, into complex and simple ions. In this case (47) will not hold mainly because of the fact that the expressions for diffusion, migration, and nondissociated flows of molecules must be determined relative to some mean liquid velocity. In concentrated solutions, this velocity does not coincide with the solvent velocity and must be determined from the fluid dynamics equations of a multicomponent mixture in which characteristics of its components (physical density, charge, diffusion coefficient of a complex ion) are, in fact, unknown. This is a reason why the theory of diffusion and migration of ions for a partially dissociated electrolyte solution, with a current traversing through it, encounters crucial difficulties. Besides, it is rather difficult to take an account of the force interaction of complex cations and anions between each other and with an external electric field. Therefore Einstein's formula (44) will also change its form.The main drawbacks of the theory of electrolytic dissociation are fully defined in collected papers “Fundamental Principles of Chemistry" by D. I. Mendeleev as well as in [40, 41].In our opinion, the principal disadvantage of the approaches considered above lies in the fact that the Nernst-Planck equation as well as its modified form for a partially dissociated electrolyte (see (47)) are based on the hydrodynamic theory of diffusion. We shall substantiate below that it is more reasonable to use the equations of thermodynamics of irreversible processes.When applying the Nernst-Planck equation (see (46)) and its modified forms (see (47)), it is implicitly postulated that cations and anions in the electrolyte solution are different components. The equations of transfer of cations and anions are written separately. On the other hand, the commonly used equation of electroneutrality does not ensure, according to Gibbs [39], the independence of differentials of cations ∂ni and anions ∂ne in a volumetric electrolyte solution their changes are functionally related and therefore they are identical, indiscernible components.


### 3.1.1. The Case of a Current Traversing a Binary Fully Dissociated Electrolyte

Let us write system of equations (see (46)) in the form:(50)∂ne∂τ+vgradne=DeΔne+zeFDeRTdiv⁡(neE),∂ni∂τ+vgradni=DiΔni+ziFDiRTdiv⁡(niE),(51)zene-zini=0.
Analogously to [39], instead of concentrations ne and ni we introduce the molar concentration related to ne and ni by the expressions(52)C=nezi=nize.
Expressing ne and ni in (50) in terms of the molar concentration, we arrive at the following expression for the function C:(53)∂C∂τ+υ∇C=DefΔC,
where(54)Def=DeDi(ze+zi)zeDe+ziDi
is the diffusion coefficient of a salt or the effective diffusion coefficient of a binary electrolyte. Expressing the concentrations ne and ni in terms of the molar concentration C we find the expression for the vector of current density J [42]:(55)J=(Di-De)Fzize∇C+F2zizeRT(ziDi+zeDe)CE.
In [43] it is shown that by its meaning Def is a coefficient of molecular diffusion. For instance, the diffusion coefficients of copper sulfate have an intermediate value between the diffusion coefficients of copper and sulfate ions; if De=0.173·10-9 m2/sec and Di=1.065·10-9 m2/sec, then Def=0.854·10-9 m2/sec. However, expression (54) does not offer an explanation of the abnormally high mobility of H and OH ions when their displacement proceeds by a croquet-like mechanism [44].In the system of (53)–(55) no symmetry of cross terms is observed; therefore, it is difficult to take into account “the superposition effects." Indeed, from (55) it is evident that the molar concentration gradients VC exert an influence on the current J. On the other hand, one can see from (53) that an electric field has no effect on the diffusion of the molecules, and therefore the Onsager reciprocal relation is disturbed.One may raise the objection that according to the Nernst-Planck equations it follows that the concentration gradients of cations and anions exert an influence on the potential distribution of an electric field [42]: (56)∇(ρR∇φ)=-F∑iz∇(Di∇Ci).
However, relation (56) itself has been derived by Neumann [43] on the assumption that the resultant current does not depend on the mass flow.On derivation of the equations for diffusion and migration of ions in a volumetric electrolyte solution we shall employ the relations of the thermodynamics of irreversible processes.Finally, it is worthy of note that in plasma physics there is also no symmetry of “the superposition effects" in describing the ambipolar diffusion of ions, and the theory is based on the hydrodynamic model [36].In investigating the dynamics of an electron-hole semiconductor plasma, the hydrodynamic model (50)-(51) is also used.In our opinion, electric field-induced transfer of charges entails transfer of their kinetic energy as well as of heat and mass, and, conversely, mass or heat transfer may simultaneously cause, if we are concerned with the system of charged particles, charge transfer and give rise to an electromotive force.


### 3.1.2. Conclusion

To phenomenologically describe diffusional-electrical phenomena in electrolytes it is suggested to use the methods of the thermodynamics of irreversible processes, rather than the hydrodynamic theory of diffusion of ions, for the experimentally observable thermodynamic flows and forces without clear discernment and identification ofrealion flows and mobilities.


## 3.1.1. The Case of a Current Traversing a Binary Fully Dissociated Electrolyte

Let us write system of equations (see (46)) in the form:(50)∂ne∂τ+vgradne=DeΔne+zeFDeRTdiv⁡(neE),∂ni∂τ+vgradni=DiΔni+ziFDiRTdiv⁡(niE),(51)zene-zini=0.
Analogously to [39], instead of concentrations ne and ni we introduce the molar concentration related to ne and ni by the expressions(52)C=nezi=nize.
Expressing ne and ni in (50) in terms of the molar concentration, we arrive at the following expression for the function C:(53)∂C∂τ+υ∇C=DefΔC,
where(54)Def=DeDi(ze+zi)zeDe+ziDi
is the diffusion coefficient of a salt or the effective diffusion coefficient of a binary electrolyte. Expressing the concentrations ne and ni in terms of the molar concentration C we find the expression for the vector of current density J [42]:(55)J=(Di-De)Fzize∇C+F2zizeRT(ziDi+zeDe)CE.
In [43] it is shown that by its meaning Def is a coefficient of molecular diffusion. For instance, the diffusion coefficients of copper sulfate have an intermediate value between the diffusion coefficients of copper and sulfate ions; if De=0.173·10-9 m2/sec and Di=1.065·10-9 m2/sec, then Def=0.854·10-9 m2/sec. However, expression (54) does not offer an explanation of the abnormally high mobility of H and OH ions when their displacement proceeds by a croquet-like mechanism [44].In the system of (53)–(55) no symmetry of cross terms is observed; therefore, it is difficult to take into account “the superposition effects." Indeed, from (55) it is evident that the molar concentration gradients VC exert an influence on the current J. On the other hand, one can see from (53) that an electric field has no effect on the diffusion of the molecules, and therefore the Onsager reciprocal relation is disturbed.One may raise the objection that according to the Nernst-Planck equations it follows that the concentration gradients of cations and anions exert an influence on the potential distribution of an electric field [42]: (56)∇(ρR∇φ)=-F∑iz∇(Di∇Ci).
However, relation (56) itself has been derived by Neumann [43] on the assumption that the resultant current does not depend on the mass flow.On derivation of the equations for diffusion and migration of ions in a volumetric electrolyte solution we shall employ the relations of the thermodynamics of irreversible processes.Finally, it is worthy of note that in plasma physics there is also no symmetry of “the superposition effects" in describing the ambipolar diffusion of ions, and the theory is based on the hydrodynamic model [36].In investigating the dynamics of an electron-hole semiconductor plasma, the hydrodynamic model (50)-(51) is also used.In our opinion, electric field-induced transfer of charges entails transfer of their kinetic energy as well as of heat and mass, and, conversely, mass or heat transfer may simultaneously cause, if we are concerned with the system of charged particles, charge transfer and give rise to an electromotive force.


## 3.1.2. Conclusion

To phenomenologically describe diffusional-electrical phenomena in electrolytes it is suggested to use the methods of the thermodynamics of irreversible processes, rather than the hydrodynamic theory of diffusion of ions, for the experimentally observable thermodynamic flows and forces without clear discernment and identification ofrealion flows and mobilities.


## 3.2. Landau Model

In what follows we consider ion conductors that contain the electrolytic plasma with a fairly high concentration of charged particles. The interaction of the charged plasma particles with the external electromagnetic field and their collective interaction can theoretically be investigated completely only on the basis of the Boltzmann kinetic equation with a self-consistent field. In this probabilistic approach to the phenomena under study, one introduces characteristics that are average over a large ensemble of particles and are always related to the introduction of additional hypotheses on the properties of particles and their interactions and to the simplification of these properties and interactions. We note that in many cases, for example, when the electrolytic plasma in liquid electrolytes is considered, there is even no basis for constructing such methods. The developed methods are usually not an efficient means of solving problems by virtue of the excessive complexity of corresponding equations.Another general method is construction of a phenomenological macroscopic theory based on the general regularities and hypotheses obtained experimentally; this precisely method will be developed in the present work. For many problems of applied character, consideration of a plasma as a conducting gas is sufficient. Such approximation is strictly justified only in the case of a dense plasma when the free path of charged particles is much smaller than the characteristic dimension of the system and particle collisions are of crucial importance. The velocity distribution of particles is Maxwellian; at each point, it is fully determined by the local values of density, temperature, and macroscopic velocity. With these conditions being observed, we can carry out, in a unified context, macroscopic phenomenological description of a gas medium, a metal and dielectric plasma, a plasma in a liquid electrolyte, or an electrolytic plasma.The necessary condition of existence of the state of local thermodynamic equilibrium (LTE) of a plasma is high frequency of collisions of plasma particles (Maxwellization) so that the plasma state undergoes no marked changes over the period between collisions on the mean free path. The observance of this condition means the following: (1) the electrons have time to transfer the energy received from the electrostatic field to heavy particles; (2) ionization processes are nearly totally counterbalanced by recombination (the Saha equation holds true); (3) the overwhelming part of excited atoms gives up its energy in collisions; (4) energy exchange between particles dominates the processes in which the plasma energy markedly grows or diminishes.In the LTE plasma state, it is only the radiation field (photons) that is not in equilibrium with the plasma, but the energy loss by radiation is minor in a number of cases compared to other mechanisms of energy transfer. With allowance for the assumptions enumerated above, the plasma electrodynamics can be described by Maxwell continuity, motion, and energy equations [45]. For macroscopic description of transfer processes in ion conductors, we will use the equations of thermodynamics of irreversible processes without separating explicitly ion carriers.The possibility of such description was first suggested by L. D. Landau and E. M. Lifshits. This approach which has been substantiated above does not provide for the introduction and determination of the coefficients of molecular diffusion of cations and anions and of the degree of dissociation in an ion conductor. It is conditioned by the interaction of the mass and charge fluxes. The ion fluxes are not determined. The density of the conduction currentJq, the mass flux Jm, and the heat flux JT for the medium’s nonequilibrium state has the form [45]:(57)Jq=λ(E-βgradn)-λα(T)gradT,Jm=-Dmgradn-DA*λE-DTgradT,JT=-kgradT-Iq(Π+φ)-λDT*gradn,
where DM is the coefficient of molecular diffusion, DA* is the coefficient of ambipolar diffusion, DT* is the coefficient allowing for heat transfer by the impurity motion, β is the coefficient of specific electrical ambipolar conductivity due to ∇n, and Π is the Peltier coefficient.To investigate the interaction of electromagnetic impotential fields (rot⁡E⃗≠0) and heat and mass transfer rates it is advisable to use the system of (57) for polar vectors; if needed, the pseudovector of the magnetic field strength can be calculated from the Maxwell equation (7).


## 3.3. Interaction of Nonstationary Electric, Thermal, and Diffusion Fields with Allowance for Mass Transfer in a Layered Medium with the Example of an Electrochemical Cell

### 3.3.1. Introduction

Investigation of the interaction of electric and thermal fields with allowance for mass transfer and contact phenomena is a complex and topical problem of the theory and practice of various fields of natural science and technology.The work seeks to construct a physicomathematical model of the interaction of nonstationary electric fields in a layered medium with allowance for nonstationary thermal phenomena and mass transfer without explicit separation of charge carriers. The media in contact are considered to be homogeneous. For the sake of clear representation, we will consider two-layer one-dimensional models.


### 3.3.2. Interaction of Electric and Thermal Fields

In different substances, the processes of transfer of charge and energy are interrelated. The quantity of the released heat is determined by Joule heating and by the effects of Thomson and Peltier. The problem of interaction of nonstationary thermal phenomena has been considered in [46] without taking into account the Thomson effect.According to [45], the expressions for the conduction-current density and the energy-flux density in the absence of an external magnetic field or in the case of its slight influence have, respectively, the form
(58)Iq=-λ(α(T)gradT+gradφ),IT=-k(T)gradT+Iq(Π+φ),
where a(T) is the specific thermoelectromotive force, Π=αT is the Peltier coefficient, φ is the potential, and k(T) is the thermal-conductivity coefficient. We note that the problem is always nonlinear.Having eliminated the magnetic-field strength from system (13) according to one method of [47–49], we obtain the equation for the electric-field strength(59)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E.
The heat-balance equation has the form [45] (60)cpρ∂T∂t=div⁡(k(T)gradT)+Iq⋅(E-(α(T)+Tβ(T))gradT).
At the interface, the following relation [50] holds:(61)div⁡i+Iqn1-Iqn2=-∂σ∂t.
Here ε is the permittivity, c2=1/(ε0·μ0), μ is the magnetic permeability, ε0 is the electric constant, μ0 is the magnetic constant, cpρ is the product of the specific heat of the medium by its density, and β(T)=∂α(T)/∂T.In deriving (60), use has been made of the condition of local electroneutrality of the substance. We give a differential formulation of the problem and the method of its solution in greater detail.Let us investigate, on the segment[0≤x<l]=[0≤x<ξ]∪[ξ≤x≤l],the contact of homogeneous media 1 and 2 with different electrophysical properties. The quantities ε,λ,μ,E, and φ have discontinuities of the first kind at the point of the interface x=ξ. We will consider the case of plane contact where the influence of surface currents can be disregarded and the thickness of the electric double-layer is much smaller than the characteristic dimension of the object. We set E=-gradφ. Equations (59) and (60) for the one-dimensional problem will take the form(62)εc2∂2E∂t2+μ0∂∂t(-λ(α(T)∂T∂x-E))=1μ∂2E∂x2,cpρ∂T∂t=∂∂x(k(T)∂T∂x)-λ(α(T)∂T∂x-E)(E-(α(T)+Tβ(T))∂T∂x).
Conditions (39) and (61) will be written, respectively, in the form(63)ε1ε0E1|x=ξ-0-ε2ε0E2|x=ξ+0=σ|x=ξ,(64)λ1(E1-α1(T)∂T∂x)|x=ξ-0-λ2(E2-α2(T)∂T∂x)|x=ξ+0=-∂σ∂t|x=ξ,
By differentiating (63) with respect to time and taking into account (64), at the phase boundary we obtain the condition of equality of the total currents:(65)[λ(E-α(T)∂T∂x)+εε0∂E∂t]|x=ξ=0.
Here and in what follows, for the arbitrary function f we adopt the notation
(66)[f]|x=ξ=f|x=ξ+0-f|x=ξ-0.
Equality (63) is a corollary of the relation lim⁡x→ε±0div⁡D=ρ [50]. Taking into consideration the finiteness of the value of the space electric charge and the continuity of it, in the one-dimensional case we obtain(67)[εε0∂E∂x]|x=ξ=0.
Relations (65) and (67) reflect the laws of conservation and continuity of electric charge [45, 50]. At the contact point, we also have the equality of the temperatures and the energy fluxes [45]:(68)[T]|x=ξ=0,(69)[k(T)∂T∂x+λ(α(T)∂T∂x-E)(α(T)T+φ)]|x=ξ=0.
Thus, in the presence of the interaction of the electric and thermal fields in a layered medium, the equality of the charge fluxes (65), the equality of the charges (67) (when the conditions of quasineutrality of the contacting media beyond the electric double-layer are satisfied), the equality of the temperatures (68), and the equality of the energy fluxes must be fulfilled at the interface of the media. In the relations for the charge and energy fluxes, we take into account cross-thermoelectrical phenomena.The process of charge of the electric field is considered on the finite time interval [0≤t≤t0]. The initial conditions have the form(70)T(x,0)=T0(x),E(x,0)=f1(x),∂E(x,0)∂t=f2(x).
We give the boundary conditions, for example, for the value x=0. We set the value of the total-current density to be known(71)(λ(α(T)∂T∂x-E)-εε0∂E∂t)|x=0=j1(t),
and take into account the heat exchange at the boundary using the Newton relation(72)(k(T)∂T∂x+λ(α(T)∂T∂x-E)α(T)T+φ)|x=0=γ1(T-T*)|x=0,
where γ1 is the coefficient of heat exchange and T* is the ambient temperature.We consider the analogous relations at the right-hand boundary.


### 3.3.3. Interaction of Electric and Thermal Fields with Allowance for Mass Transfer

Let us consider the nonstationary model of heat and mass transfer for electrochemical systems with the example of electrolysis. For this problem, an electric double-layer occurs at the site of contact of the electrolyte with the metallic anode and cathode. Let us set the density of the space charge beyond the electric double-layer of the metal-electrolyte contact to be zero at the initial instant and hence, according to [51], remain constant in the future while the voltage drop in the electrodes and leads is small. The influence of the electrodes on the temperature field of the electrolyte will be taken into account in terms of the coefficient of heat exchange.Two approaches to modeling of diffusion-electrical phenomena have currently been developed, each of which has certain disadvantages and advantages. The first approach [37, 38, 42] is characterized by consideration of the flows of ions of the corresponding sort in a completely or partially dissociated electrolyte. The equations derived contain many parameters, the reliable procedure of determination of which is absent in the majority of cases. Furthermore, the proposed theory fails to provide for taking into account the interaction of the cations and the anions of the dissolved substance with each other and with the molecules of the solvent. It must also be borne in mind that, apart from the dissociation reaction, we have the association of simple molecules which in turn are dissociated into complex and simple ions. It is difficult to determine the diffusion coefficients of complex ions; it is also difficult to take into account the force interaction of complex cations and anions with each other and with the external electric field. The difficulties in question become more serious when a multicomponent electrolyte is described. We note that such difficulties also arise in modeling, for example, certain problems of plasma physics.The second approach [52–56] fails to provide for the introduction and determination of the coefficients of molecular diffusion of cations and anions and the degree of dissociation of the electrolyte. It is based on the interaction of the mass and charge fluxes. The ionic flows are not determined. The model is formulated without explicit separation of the charge carriers and the charge of an electric double layer.The system of equations describing the electrodynamic processes in the electrolyte is as follows:(73)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E,∂n∂t=-div⁡(Im),c0ρ0∂T∂t=-div⁡(IT).
Here Iq, Im, and IT are the densities of the conduction current and of the mass and heat fluxes which (in the case of the nonequilibrium state of the medium) have, respectively, the form (57) [45], where DM is the molecular diffusion coefficient, DA* is the ambipolar diffusion coefficient, DT is the thermal diffusion coefficient, DT* is the coefficient taking into account the transfer of heat due to the motion of the impurity, and n is the concentration of the impurity.We give the initial(74)E(x,0)=E0(x),∂E(x,0)∂t=Ẽ(x),T(x,0)=T0(x),n(x,0)=n0(x),
and boundary conditions (anode-electrolyte, x=0):(75)Iq=λ(λA*∂n∂x-E)+λα(T)∂T∂x-εε0∂E∂t,γ1(T-T*)=k∂T∂x-Iq(Π+φ)+DT*∂n∂x,kaIq(t)=DM∂n∂x-DA*λE+DT∂T∂x,
in modeling the one-dimensional problem in the region G=[0≤x≤l]×[0≤t≤T].Here ka is the electrochemical equivalent of a substance which deposits at the anode. The analogous conditions hold for x=l as well.The constructed system of differential equations makes it possible to model the processes of transfer in the electrochemical system with allowance for the influence of an electric double-layer; this influence directly determines the thermoelectrical and ambipolar electrodiffusion phenomena at the boundary of the metal-electrolyte contact.


### 3.3.4. Numerical Modeling of Electrochemical Systems. Heating of an Electrochemical Cell in the Case of Constant-Current and Pulse Electrolysis

In electrochemical systems, the electrodes are metals as a rule. The electrical conductivity of metals is hundreds of times higher than the specific electrical conductivity of electrolytes; therefore, one can disregard the voltage drop in electrodes and leads. In highly conductive metals, one can also disregard the displacement currents and the coefficients of ambipolar diffusion and thermal diffusion and restrict oneself to consideration of the system of equations (73)–(75) just in the electrolyte itself.Let us consider the process of heating of an electrochemical medium in the case of passage of an electric current with densityIq(t) through a solution of copper sulfate CuSO4·5H2O. We use copper (99.78%) as the anode. The copper from the solution is deposited at the cathode. We set the current efficiencies for the copper equal to 100% and the electrochemical equivalent ke equal to 0.6588·10−6 kg/C. The dependence of the specific electrical conductivity λ(n) of the copper-plating electrolyte on the concentration of the copper sulfate in water is given in [57]; the dependence of λ and k on the temperature will be considered to be insignificant. In the calculations for CuSO4, we set DM=5·10-10 m2/sec, DA*=10-11λ,λA*=10-4λ,ε=70, and μ=1 and disregard the effects of thermal diffusion. The distance between the cathode and the anode is L=0.05 m. The heat capacity, the density, and the thermal-conductivity coefficient of the electrolyte are taken to be 4.2·103 J/(kg·K), 103 kg/m3, and 0.6 V/(m·K) respectively. We note that the regimes of pulse action have been considered in [58] as applied to the processes of electrodeposition of alloys. This investigation has been carried out on the basis of a separate description of the transport of ions and the employment of the Kirchhoff law for quasistationary currents in a cell. With such an approach, the displacement currents are taken into account indirectly through introduction of the capacitive current of the electrode and determination of the experimental dependence of the polarization capacitance of the electrode on the character of the pulse action. In [58], it has been noted that neglect of the capacitive current in pulse electrolysis involves significant errors. The employment of system (73)–(75) to model the process of copper plating makes it possible to consider the problems without introducing the concepts of inductance and capacitance. The obtained nonlinear system of equations is solved by the finite-difference method analogously to [54].We have modeled numerically the process of copper plating for a constant current and a pulse current (Figure16). The density of the constant current was equal to 300 A/m2. The maximum density of the pulse current was also 300 A/m2, while the period of traversal of the current and the break were equal to 0.01 sec.Distribution of the normalized concentration of CuSO4  and the temperature in an electrochemical cell within 60 sec after switching on the current. The solid curve, at constant current; the dashed curve, at pulse current. T,°C; x,mm.
(a)
Distribution of the normalized concentration of CuSO4(b)
Distribution of the temperatureFigure16 gives the results of modeling the distributions of the concentration of CuSo4 and the temperature for different regimes of electrolysis. The measurements were carried out within 60 sec after switching on the current. It is seen from the figure that the concentration gradients near the cathode surface (x=1) are different and depend on the regime of copper plating. In the case of nonstationary action they are substantially smaller than in the case of electrolysis in the regime of constant current. As has been shown in [58], this tendency also holds in the case where one and the same total charge passes through the electrochemical cell.The analysis of the temperature field (Figure16(b)) shows that the temperatures of the electrolyte in the vicinity of the anode and the cathode differ significantly. Nonsymmetric heating of the electrochemical cell attributed to the am-bipolar diffusion and electrical ambipolar conductivity occurs. The modeling results are in agreement with the data of [59].With the aim of measuring the temperature in the surface layer of the electrolyte solution near the cathode and comparing the results to the calculated data, we developed and manufactured an experimental setup. Copper anode and cathode were placed in an electroplating bath with a sulfuric-acid copper-plating electrolyte. In the 1 mm thick cathode, a hole was made into which we placed a temperature-sensitive element. As the latter we employed a thermal resistor. The current strength in the bath was prescribed by a variable resistor.Figure17 gives the results of numerical modeling and the experimental data on heating of the cathode region in relation to the regime of current. It is seen that the increase in the density of the current (both constant and pulse) causes the temperature near the cathode surface to increase. In the anode, no increase in the temperature was observed either in the calculations or in the experiment, which is in agreement with the data of [59]. The surface temperature obtained in solving (73)–(75) in the regime of constant current is 10–15% higher than the temperature recorded in the experiments.Figure 17
Heating of the electrolyte in the cathode region as a function of the value and kind of polarizing current: (1) constant current; (2) pulse current; (3, 4) corresponding experimental data.ΔT,°C; Iq, A/m2.A more complicated situation occurs when we compare the experimental data and the results of numerical modeling in the regime of pulse electrolysis. The results of the experiments demonstrate that in the regime of pulse current the temperature of the electrolyte in the cathode region increases more than in the regime of constant current. At the same time, the numerical calculations of the pulse electrolysis show that the heating of the cathode region here is smaller than for the regime of constant current. The difference is, apparently, attributed to concentration convection [60], whose influence is substantial for high densities of the current for pulse electrolysis.In a longer term, one must take into account natural and concentration convection and heat exchange with the ambient medium. Nonetheless, despite the errors in determination of the coefficients and the assumptions made in the model, the calculation results are in qualitative agreement with the experimental data, which makes it possible, in certain cases, to use the proposed approach for modeling of nonstationary processes.Results were published in part [52, 53, 55, 56, 61].


## 3.3.1. Introduction

Investigation of the interaction of electric and thermal fields with allowance for mass transfer and contact phenomena is a complex and topical problem of the theory and practice of various fields of natural science and technology.The work seeks to construct a physicomathematical model of the interaction of nonstationary electric fields in a layered medium with allowance for nonstationary thermal phenomena and mass transfer without explicit separation of charge carriers. The media in contact are considered to be homogeneous. For the sake of clear representation, we will consider two-layer one-dimensional models.


## 3.3.2. Interaction of Electric and Thermal Fields

In different substances, the processes of transfer of charge and energy are interrelated. The quantity of the released heat is determined by Joule heating and by the effects of Thomson and Peltier. The problem of interaction of nonstationary thermal phenomena has been considered in [46] without taking into account the Thomson effect.According to [45], the expressions for the conduction-current density and the energy-flux density in the absence of an external magnetic field or in the case of its slight influence have, respectively, the form
(58)Iq=-λ(α(T)gradT+gradφ),IT=-k(T)gradT+Iq(Π+φ),
where a(T) is the specific thermoelectromotive force, Π=αT is the Peltier coefficient, φ is the potential, and k(T) is the thermal-conductivity coefficient. We note that the problem is always nonlinear.Having eliminated the magnetic-field strength from system (13) according to one method of [47–49], we obtain the equation for the electric-field strength(59)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E.
The heat-balance equation has the form [45] (60)cpρ∂T∂t=div⁡(k(T)gradT)+Iq⋅(E-(α(T)+Tβ(T))gradT).
At the interface, the following relation [50] holds:(61)div⁡i+Iqn1-Iqn2=-∂σ∂t.
Here ε is the permittivity, c2=1/(ε0·μ0), μ is the magnetic permeability, ε0 is the electric constant, μ0 is the magnetic constant, cpρ is the product of the specific heat of the medium by its density, and β(T)=∂α(T)/∂T.In deriving (60), use has been made of the condition of local electroneutrality of the substance. We give a differential formulation of the problem and the method of its solution in greater detail.Let us investigate, on the segment[0≤x<l]=[0≤x<ξ]∪[ξ≤x≤l],the contact of homogeneous media 1 and 2 with different electrophysical properties. The quantities ε,λ,μ,E, and φ have discontinuities of the first kind at the point of the interface x=ξ. We will consider the case of plane contact where the influence of surface currents can be disregarded and the thickness of the electric double-layer is much smaller than the characteristic dimension of the object. We set E=-gradφ. Equations (59) and (60) for the one-dimensional problem will take the form(62)εc2∂2E∂t2+μ0∂∂t(-λ(α(T)∂T∂x-E))=1μ∂2E∂x2,cpρ∂T∂t=∂∂x(k(T)∂T∂x)-λ(α(T)∂T∂x-E)(E-(α(T)+Tβ(T))∂T∂x).
Conditions (39) and (61) will be written, respectively, in the form(63)ε1ε0E1|x=ξ-0-ε2ε0E2|x=ξ+0=σ|x=ξ,(64)λ1(E1-α1(T)∂T∂x)|x=ξ-0-λ2(E2-α2(T)∂T∂x)|x=ξ+0=-∂σ∂t|x=ξ,
By differentiating (63) with respect to time and taking into account (64), at the phase boundary we obtain the condition of equality of the total currents:(65)[λ(E-α(T)∂T∂x)+εε0∂E∂t]|x=ξ=0.
Here and in what follows, for the arbitrary function f we adopt the notation
(66)[f]|x=ξ=f|x=ξ+0-f|x=ξ-0.
Equality (63) is a corollary of the relation lim⁡x→ε±0div⁡D=ρ [50]. Taking into consideration the finiteness of the value of the space electric charge and the continuity of it, in the one-dimensional case we obtain(67)[εε0∂E∂x]|x=ξ=0.
Relations (65) and (67) reflect the laws of conservation and continuity of electric charge [45, 50]. At the contact point, we also have the equality of the temperatures and the energy fluxes [45]:(68)[T]|x=ξ=0,(69)[k(T)∂T∂x+λ(α(T)∂T∂x-E)(α(T)T+φ)]|x=ξ=0.
Thus, in the presence of the interaction of the electric and thermal fields in a layered medium, the equality of the charge fluxes (65), the equality of the charges (67) (when the conditions of quasineutrality of the contacting media beyond the electric double-layer are satisfied), the equality of the temperatures (68), and the equality of the energy fluxes must be fulfilled at the interface of the media. In the relations for the charge and energy fluxes, we take into account cross-thermoelectrical phenomena.The process of charge of the electric field is considered on the finite time interval [0≤t≤t0]. The initial conditions have the form(70)T(x,0)=T0(x),E(x,0)=f1(x),∂E(x,0)∂t=f2(x).
We give the boundary conditions, for example, for the value x=0. We set the value of the total-current density to be known(71)(λ(α(T)∂T∂x-E)-εε0∂E∂t)|x=0=j1(t),
and take into account the heat exchange at the boundary using the Newton relation(72)(k(T)∂T∂x+λ(α(T)∂T∂x-E)α(T)T+φ)|x=0=γ1(T-T*)|x=0,
where γ1 is the coefficient of heat exchange and T* is the ambient temperature.We consider the analogous relations at the right-hand boundary.


## 3.3.3. Interaction of Electric and Thermal Fields with Allowance for Mass Transfer

Let us consider the nonstationary model of heat and mass transfer for electrochemical systems with the example of electrolysis. For this problem, an electric double-layer occurs at the site of contact of the electrolyte with the metallic anode and cathode. Let us set the density of the space charge beyond the electric double-layer of the metal-electrolyte contact to be zero at the initial instant and hence, according to [51], remain constant in the future while the voltage drop in the electrodes and leads is small. The influence of the electrodes on the temperature field of the electrolyte will be taken into account in terms of the coefficient of heat exchange.Two approaches to modeling of diffusion-electrical phenomena have currently been developed, each of which has certain disadvantages and advantages. The first approach [37, 38, 42] is characterized by consideration of the flows of ions of the corresponding sort in a completely or partially dissociated electrolyte. The equations derived contain many parameters, the reliable procedure of determination of which is absent in the majority of cases. Furthermore, the proposed theory fails to provide for taking into account the interaction of the cations and the anions of the dissolved substance with each other and with the molecules of the solvent. It must also be borne in mind that, apart from the dissociation reaction, we have the association of simple molecules which in turn are dissociated into complex and simple ions. It is difficult to determine the diffusion coefficients of complex ions; it is also difficult to take into account the force interaction of complex cations and anions with each other and with the external electric field. The difficulties in question become more serious when a multicomponent electrolyte is described. We note that such difficulties also arise in modeling, for example, certain problems of plasma physics.The second approach [52–56] fails to provide for the introduction and determination of the coefficients of molecular diffusion of cations and anions and the degree of dissociation of the electrolyte. It is based on the interaction of the mass and charge fluxes. The ionic flows are not determined. The model is formulated without explicit separation of the charge carriers and the charge of an electric double layer.The system of equations describing the electrodynamic processes in the electrolyte is as follows:(73)εc2∂2E∂t2+μ0∂Iq∂t=1μ∇2E,∂n∂t=-div⁡(Im),c0ρ0∂T∂t=-div⁡(IT).
Here Iq, Im, and IT are the densities of the conduction current and of the mass and heat fluxes which (in the case of the nonequilibrium state of the medium) have, respectively, the form (57) [45], where DM is the molecular diffusion coefficient, DA* is the ambipolar diffusion coefficient, DT is the thermal diffusion coefficient, DT* is the coefficient taking into account the transfer of heat due to the motion of the impurity, and n is the concentration of the impurity.We give the initial(74)E(x,0)=E0(x),∂E(x,0)∂t=Ẽ(x),T(x,0)=T0(x),n(x,0)=n0(x),
and boundary conditions (anode-electrolyte, x=0):(75)Iq=λ(λA*∂n∂x-E)+λα(T)∂T∂x-εε0∂E∂t,γ1(T-T*)=k∂T∂x-Iq(Π+φ)+DT*∂n∂x,kaIq(t)=DM∂n∂x-DA*λE+DT∂T∂x,
in modeling the one-dimensional problem in the region G=[0≤x≤l]×[0≤t≤T].Here ka is the electrochemical equivalent of a substance which deposits at the anode. The analogous conditions hold for x=l as well.The constructed system of differential equations makes it possible to model the processes of transfer in the electrochemical system with allowance for the influence of an electric double-layer; this influence directly determines the thermoelectrical and ambipolar electrodiffusion phenomena at the boundary of the metal-electrolyte contact.


## 3.3.4. Numerical Modeling of Electrochemical Systems. Heating of an Electrochemical Cell in the Case of Constant-Current and Pulse Electrolysis

In electrochemical systems, the electrodes are metals as a rule. The electrical conductivity of metals is hundreds of times higher than the specific electrical conductivity of electrolytes; therefore, one can disregard the voltage drop in electrodes and leads. In highly conductive metals, one can also disregard the displacement currents and the coefficients of ambipolar diffusion and thermal diffusion and restrict oneself to consideration of the system of equations (73)–(75) just in the electrolyte itself.Let us consider the process of heating of an electrochemical medium in the case of passage of an electric current with densityIq(t) through a solution of copper sulfate CuSO4·5H2O. We use copper (99.78%) as the anode. The copper from the solution is deposited at the cathode. We set the current efficiencies for the copper equal to 100% and the electrochemical equivalent ke equal to 0.6588·10−6 kg/C. The dependence of the specific electrical conductivity λ(n) of the copper-plating electrolyte on the concentration of the copper sulfate in water is given in [57]; the dependence of λ and k on the temperature will be considered to be insignificant. In the calculations for CuSO4, we set DM=5·10-10 m2/sec, DA*=10-11λ,λA*=10-4λ,ε=70, and μ=1 and disregard the effects of thermal diffusion. The distance between the cathode and the anode is L=0.05 m. The heat capacity, the density, and the thermal-conductivity coefficient of the electrolyte are taken to be 4.2·103 J/(kg·K), 103 kg/m3, and 0.6 V/(m·K) respectively. We note that the regimes of pulse action have been considered in [58] as applied to the processes of electrodeposition of alloys. This investigation has been carried out on the basis of a separate description of the transport of ions and the employment of the Kirchhoff law for quasistationary currents in a cell. With such an approach, the displacement currents are taken into account indirectly through introduction of the capacitive current of the electrode and determination of the experimental dependence of the polarization capacitance of the electrode on the character of the pulse action. In [58], it has been noted that neglect of the capacitive current in pulse electrolysis involves significant errors. The employment of system (73)–(75) to model the process of copper plating makes it possible to consider the problems without introducing the concepts of inductance and capacitance. The obtained nonlinear system of equations is solved by the finite-difference method analogously to [54].We have modeled numerically the process of copper plating for a constant current and a pulse current (Figure16). The density of the constant current was equal to 300 A/m2. The maximum density of the pulse current was also 300 A/m2, while the period of traversal of the current and the break were equal to 0.01 sec.Distribution of the normalized concentration of CuSO4  and the temperature in an electrochemical cell within 60 sec after switching on the current. The solid curve, at constant current; the dashed curve, at pulse current. T,°C; x,mm.
(a)
Distribution of the normalized concentration of CuSO4(b)
Distribution of the temperatureFigure16 gives the results of modeling the distributions of the concentration of CuSo4 and the temperature for different regimes of electrolysis. The measurements were carried out within 60 sec after switching on the current. It is seen from the figure that the concentration gradients near the cathode surface (x=1) are different and depend on the regime of copper plating. In the case of nonstationary action they are substantially smaller than in the case of electrolysis in the regime of constant current. As has been shown in [58], this tendency also holds in the case where one and the same total charge passes through the electrochemical cell.The analysis of the temperature field (Figure16(b)) shows that the temperatures of the electrolyte in the vicinity of the anode and the cathode differ significantly. Nonsymmetric heating of the electrochemical cell attributed to the am-bipolar diffusion and electrical ambipolar conductivity occurs. The modeling results are in agreement with the data of [59].With the aim of measuring the temperature in the surface layer of the electrolyte solution near the cathode and comparing the results to the calculated data, we developed and manufactured an experimental setup. Copper anode and cathode were placed in an electroplating bath with a sulfuric-acid copper-plating electrolyte. In the 1 mm thick cathode, a hole was made into which we placed a temperature-sensitive element. As the latter we employed a thermal resistor. The current strength in the bath was prescribed by a variable resistor.Figure17 gives the results of numerical modeling and the experimental data on heating of the cathode region in relation to the regime of current. It is seen that the increase in the density of the current (both constant and pulse) causes the temperature near the cathode surface to increase. In the anode, no increase in the temperature was observed either in the calculations or in the experiment, which is in agreement with the data of [59]. The surface temperature obtained in solving (73)–(75) in the regime of constant current is 10–15% higher than the temperature recorded in the experiments.Figure 17
Heating of the electrolyte in the cathode region as a function of the value and kind of polarizing current: (1) constant current; (2) pulse current; (3, 4) corresponding experimental data.ΔT,°C; Iq, A/m2.A more complicated situation occurs when we compare the experimental data and the results of numerical modeling in the regime of pulse electrolysis. The results of the experiments demonstrate that in the regime of pulse current the temperature of the electrolyte in the cathode region increases more than in the regime of constant current. At the same time, the numerical calculations of the pulse electrolysis show that the heating of the cathode region here is smaller than for the regime of constant current. The difference is, apparently, attributed to concentration convection [60], whose influence is substantial for high densities of the current for pulse electrolysis.In a longer term, one must take into account natural and concentration convection and heat exchange with the ambient medium. Nonetheless, despite the errors in determination of the coefficients and the assumptions made in the model, the calculation results are in qualitative agreement with the experimental data, which makes it possible, in certain cases, to use the proposed approach for modeling of nonstationary processes.Results were published in part [52, 53, 55, 56, 61].


## 4. Interaction of Nonstationary Electric and Thermal Fields with Allowance for Relaxation Processes

We investigate electric and thermal fields created by macroscopic charges and currents in continuous media. Of practical interest is modeling of local heat releases in media on exposure to a high-frequency electromagnetic field. We should take into account the influence of the energy absorption on the propagation of an electromagnetic wave, since the transfer processes are interrelated.In an oscillatory circuit with continuously distributed parameters, the energy dissipation is linked [62] to the dielectric loss due to the dependence of the relative permittivity ε(ω) on frequency. In the general case ε is also complex, and the relationship between the electric displacement and electric field vectors has the form D=ε(ω)E, where ε(ω)=ε′(ω)-iε′′(ω); here, ε′ and ε′′ are determined experimentally. As of now, the problems of dielectric heating of a continuous medium are reduced in many cases to consideration of an equivalent circuit based on lumped parameters, such as capacitance, inductance, loss angle, and relative-dielectric-loss factor [63, 64], that are established experimentally.With this approach, there arise substantial difficulties in determining the temperature field of equivalent circuits. Also, we have polarization and the occurrence of an electric double-layer of a prescribed electric moment in contact of media with different properties. Equivalent circuits in lamellar media additionally involve empirical lumped parameters: surface capacitance and surface resistance [65]. The total current can always be separated into a dissipative or conduction current which is in phase with the applied voltage and a displacement current shifted in time relative to the voltage. The exact physical meaning of these components of the current is largely dependent on selection of an equivalent electric circuit. A unique equivalent circuit—series or parallel connection of the capacitor, the resistor, and the inductor—does not exist; it is determined by a more or less adequate agreement with experimental data.In the case of electrolytic capacitors, the role of one plate is played by the electric double-layer with a specific resistance much higher than the resistance of metallic plates. Therefore, decrease in the capacitance with frequency is observed, for such capacitors, even in the acoustic-frequency range [65]. Circuits equivalent to an electrolytic capacitor are very bulky: up to 12 R,L, and C elements can be counted in them; therefore, it is difficult to obtain a true value of, for example, the electrolyte capacitance. In [65] experimental methods of measurement of the dielectric properties of electrolyte solutions at different frequencies are given and ε′ and ε′′are determined. The frequency dependence of dispersion and absorption are essentially different consequences of one phenomenon: “dielectric-polarization inertia” [65]. In actual fact, the dependence ε(ω) is attributable to the presence of the resistance of the electric double-layer and to the electrochemical cell in the electrolytic capacitor being a system with continuously distributed parameters, in which the signal velocity is a finite quantity.Actually,ε′ and ε′′ are certain integral characteristics of a material at a prescribed constant temperature, which are determined by the geometry of the sample and the properties of the electric double-layer. It is common knowledge that in the case of a field arbitrarily dependent on time any reliable calculation of the absorbed energy in terms of ε(ω) turns out to be impossible [45]. This can only be done for a specific dependence of the field E on time. For a quasimonochromatic field, we have [45](76)E(t)=12[E0(t)e-iωt+E0*(t)eiωt],H(t)=12[H0(t)e-iωt+H0*(t)eiωt].
The values of E0(t) and H0(t), according to [66], must very slowly vary over the period T=2π/ω. Then, for absorbed energy, on averaging over the frequency ω, we obtain the expression [66]:(77)∂D(t)∂tE(t)=14d(ωε′(ω))d(ω)∂∂t[E0(t)E0*(t)]+ωε′′2E0(t)E0*(t)+id(ωε′′(ω))4dω(∂E0(t)∂tE0*(t)-E0*(t)∂tE0(t)),
where the derivatives with respect to frequency are taken at the carrier frequency ω. We note that for an arbitrary function E(t), it is difficult to represent it in the form(78)E(t)=a(t)cos⁡φ(t),
since we cannot unambiguously indicate the amplitude a(t) and the phase φ(t). The manner in which E(t) is decomposed into factors a and cos⁡φ is not clear. Even greater difficulties appear in the case of going to the complex representation W(t)=U(t)+iV(t) when the real oscillation E(t) is supplemented with the imaginary part V(t). The arising problems have been considered in [67] in detail. In the indicated work, it has been emphasized that certain methods using a complex representation and claiming higher-than-average accuracy become trivial without an unambiguous determination of the amplitude, phase, and frequency.Summing up the aforesaid, we can state that calculation of the dielectric loss is mainly empirical in character. Construction of the equivalent circuit and allowance for the influence of the electric double-layer and for the dependence of electrophysical properties on the field’s frequency are only true of the conditions under which they have been modeled; therefore, these are fundamental difficulties in modeling the propagation and absorption of electromagnetic energy.As we believe, the release of heat in media on exposure to nonstationary electric fields can be calculated on the basis of allowance for the interaction of electromagnetic and thermal fields as a system with continuously distributed parameters from the field equation and the energy equation which take account of the distinctive features of the boundary between adjacent media. When the electric field interacting with a material medium is considered we use Maxwell equations (6)-(7). We assume that space charges are absent from the continuous medium at the initial instant of time and they do not appear throughout the process. The energy equation will be represented in the form (79)ρCpdTdt=div⁡[k(T)grad(T)]+Q,
where Q is the dissipation of electromagnetic energy.According to [68], the electromagnetic energy converted to heat is determined by the expression(80)Q=ρ[Eddt(Dρ)+Hddt(Bρ)]+JqE.
In deriving this formula, we used the nonrelativistic approximation of Minkowski’s theory. If ε,μ, and ρ=const, there is no heat release; therefore, the intrinsic dielectric loss is linked to the introduction of ε′(ω) and ε′′(ω). The quantity Q is affected by the change in the density of the substance ρ(T).A characteristic feature of high frequencies is the lag of the polarization field behind the charge in the electric field in time. Therefore, the electric-polarization vector is expediently determined by solution of the equationP(t+τe)=(ε-1)ε0E(t) with allowance for the time of electric relaxation of dipoles τe. Restricting ourselves to the first term of the expansion P(t+τe) in a Taylor series, from this equation, we obtain(81)P(t)+τedP(t)dt=(ε-1)ε0E(t).
The solution (81), on condition that P=0 at the initial instant of time, will take the form(82)P=(ε-1)ε0τe∫t0tE(τ)e-(t-τ)/τedτ.
It is noteworthy that (81) is based on the classical Debay model. According to this model, particles of a substance possess a constant electric dipole moment. The indicated polarization mechanism involves partial arrangement of dipoles along the electric field, which is opposed by the process of disorientation of dipoles because of thermal collisions. The restoring “force”, in accordance with (81), does not lead to oscillations of electric polarization. It acts as if constant electric dipoles possessed strong damping.Molecules of many liquids and solids possess the Debay relaxation polarizability. Initially polarization aggregates of Debay oscillators turn back to the equilibrium stateP(t)=P(0)exp⁡(-t/τe).A dielectric is characterized, as a rule, by a large set of relaxation times with a characteristic distribution function, since the potential barrier limiting the motion of weakly coupled ions may have different values [63]; therefore, the mean relaxation time of the ensemble of interacting dipoles should be meant by τe in (81).To eliminate the influence of initial conditions and transient processes we sett0=-∞, E(∞)=0, H(∞)=0, as it is usually done. If the boundary regime acts for a fairly long time, the influence of initial data becomes weaker with time owing to the friction inherent in every real physical system. Thus, we naturally arrive at the problem without the initial conditions:(83)P=(ε-1)ε0τe∫-∞tE(τ)e-(t-τ)/τedτ.
Let us consider the case of the harmonic field E=E0sin⁡ωt; then, using (83) we have, for the electric induction vector(84)D=ε0E+P=(ε-1)ε0τe∫-∞tE(τ)e-(t-τ)/τedτ+ε0E0sin⁡ωt=E0(ε-1)ε01+ω2τe2(sin⁡ωt-ωτecos⁡ωt)+ε0E0sin⁡ωt.
The electric induction vector is essentially the sum of two absolutely different physical quantities: the field strength and the polarization of a unit volume of the medium.If the change in the density of the substance is small, we obtain, from formula (80), for the local instantaneous heat release(85)Q=EdDdt=E02(ε-1)ε01+ω2τe2(ωsin⁡ωtcos⁡ωt+ω2τesin2ωt),
when we write the mean value of Q over the total period T:(86)Q=12E02(ε-1)ε01+ω2τe2ω2τe+λE22.
For high frequencies (ω→∞), heat release ceases to be dependent on frequency, which is consistent with formula (86) and experiment [63].When the relaxation equation for the electric field is used we must also take account of the delay of the magnetic field, when the magnetic polarization lags behind the change in the strength of the external magnetic field:(87)I(t)+τidI(t)dt=μμ0H(t).
Formula (86) is well known in the literature; it has been obtained by us without introducing complex parameters. In the case of “strong” heating of a material where the electrophysical properties of the material are dependent on temperature expression (81) will have a more complicated form and the expression for Q can only be computed by numerical methods. Furthermore, in the presence of strong field discontinuities, we cannot in principle obtain the expression for Q because of the absence of closing relations for the induced surface charge and the surface current on the boundaries of adjacent media; therefore, the issue of energy relations in macroscopic electrodynamics is difficult, particularly, with allowance for absorption.Energy relations in a dispersive medium have repeatedly been considered; nonetheless, in the presence of absorption, the issue seems not clearly understood (or at least not sufficiently known), particularly in the determination of the expression of released heat on the boundaries of adjacent media.Indeed, it is known from the thermodynamics of dielectrics that the differential of the free energyF has the form(88)dF=-SdT-PdV+EdD.
If the relative permittivity and the temperature and volume of the dielectric are constant quantities, from (88) we have(89)F(T,D)=F0+εε0E22,
where F0 is the free energy of the dielectrics in the absence of the field.The change of the internal energy of the dielectrics during its polarization at constant temperature and volume can be found from the Gibbs-Helmholtz equation. DisregardingF0 which is independent of the field strength, we can obtain(90)U(T,D)=F(T,D)-T(dFdT)D.
If the relative dielectric constant is dependent on temperature (ε(T)), we obtain(91)U(T,D)=ε0E22+(ε+T(dFdT)D).
Expression (91) determines the change in the internal energy of the dielectric in its isothermal polarization but with allowance for the energy transfer to a thermostat, if the polarization causes the dielectric temperature to change. In the works on microwave heating, that we know, (91) is not used.A characteristic feature of high frequencies is that the polarization field lags behind the change in the external field in time; therefore, the polarization vector is expediently determined by solution of the equation:(92)P(t+τe)=(ε-1+T(dεdT)D)ε0E(t).
With allowance for the relaxation time, that is, restricting ourselves to the first term of the expansion P(t+τe) in a Taylor series, we obtain(93)P(t)+τedP(t)dT=(ε-1+T(dεdT)D)ε0E(t).
In the existing works on microwave heating with the use of complex parameters, they disregard the dependence ε′′(T). In [69], consideration has been given to the incidence of a one-dimensional wave from a medium with arbitrary complex parameters on one or two boundaries of media whose parameters are also arbitrary. The amplitudes of waves reflected from and transmitted by each boundary have been found. The refection, transmission, and absorption coefficients have been obtained from the wave amplitudes. The well-known proposition that a traditional selection of determinations of the reflection, transmission, and absorption coefficients from energies (reflectivity, transmissivity, and absorptivity) in the case of complex parameters of media comes into conflict with the law of conservation of energy has been confirmed and exemplified. The necessity of allowing for ε′′(T) still further complicates the problem of computation of the dissipation of electromagnetic energy in propagation of waves through the boundaries of media with complex parameters.The proposed method of computation of local heat release is free of the indicated drawbacks and makes it possible, for the first time, to construct a consistent model of propagation of nonmonochromatic waves in a heterogeneous medium with allowance for frequency dispersion without introducing complex parameters.In closing, we note that a monochromatic wave is infinite in space and time, has infinitesimal energy absorption in a material medium, and transfers infinitesimal energy, which is the idealization of real processes. However with these stringent constraints, too, the problem of propagation of waves through the boundary is open and far from being resolved even when the complex parameters of the medium are introduced and used. In reality, the boundary between adjacent media is not infinitely thin and has finite dimensions of the electric double-layers; therefore, approaches based on through-counting schemes for a hyperbolic equation without explicit separation of the boundary between adjacent media are promising.


## 5. Ponderomotive Forces in Heterogeneous Lamellar Media with Absorption

It is common knowledge that the electric field strengthEis equal to the ponderomotive mechanical force per unit charge acting at a given point of the field on a test charge (charged body). The issue of the force acting on induced surface charges on the boundaries remains to be solved, since there are no closing equations for σ and the charge is not determined on the interface itself. According to [45], the tension tensor of ponderomotive forces for media with dispersion and absorption due to the dipole relaxation cannot be found in general form in a macroscopic manner at all. As we believe, solution of this problem lies in computing the induced surface charge on the boundaries and the resultants of all tensions applied externally to the surface without the necessity of reducing bulk ponderomotive forces to tensions. The reduction of bulk forces to tensions is not necessarily possible; therefore, this issue will be considered in greater detail.


### 5.1. Reduction of Bulk Ponderomotive Forces to Tensions

Iff is the bulk density of ponderomotive forces, the resultant of all forces applied to the bodies within the volume V is(94)F=∫VfdV.
On the other hand, if bulk forces can be reduced to tensions at all, the set of tensions acting externally of a closed surface S must be equal to the same quantity. We denote by Tn a force acting externally per unit surface the external normal to which is directed along n; the components of this force along the coordinate axes are Txn,Tyn and Tzn. Then the resultant of all tensions applied externally to the surface will be(95)F=⨖sTndS.
Using the Ostrogradsky-Gauss theorem and equating (94) and (95), we obtain(96)fx=∂Txx∂x+∂Txy∂y+∂Txz∂z,fy=∂Tyx∂x+∂Tyy∂y+∂Tyz∂z,fz=∂Tzx∂x+∂Tzy∂y+∂Tzz∂z.
For the bulk forces and tensions to be equivalent, it is necessary that not only the resultant of forces applied to an arbitrary volume remains constant but the moment of these forces be invariant as well [70]. This circumstance imposes an additional constraint on the tension-tensor components. The moment of bulk forces applied to an arbitrary volume V is(97)N=∫⌊Rf⌋dV,
where R is the distance from point O for which the moment of forces is determined to the element dV. Consequently, the component Nx must satisfy relation (97):(98)Nx=∫V(yfz-zfy)dV=⨖s(yTzn-zTyn)dS+∫V(Tyz-Tzy)dV.
The surface integral on the right is equal to the moment of tension forces Tn applied to the surface S of volume V. The moment of these tension forces will be equal to the moment of bulk forces only in the case where the last integral on the right is equal to zero. In view of the arbitrariness of the volume V, this will occur only when the integrand
(99)Tyz=Tzy
is equal to zero.By repeating the same considerations forNy and Nz, we obtain Tzx=Txy, Txy=Tyx. Consequently, for the bulk forces and surface ponderomotive forces to be equivalent, the tension tensor must be symmetric. In anisotropic media, the tension tensor is asymmetric; therefore, direct computation of the resultant of forces in terms of tensions from formula (95) is expedient. We note that Maxwell was the first to reduce the field’s ponderomotive forces to tensions. Replacement of ponderomotive forces by an equivalent system of tensions makes it much easier to determine forces applied to the volume of the electric field.The force applied to a charged layer, given the uncompensated surface charge, has the form(100)FijS=σ(Eij1-Eij2).
If the layer is not charged, we have σ=0 and FijS=0. The resultant force applied to the boundary of a cell is(101)F=⨖FijdS.
The projections of fx, fy, and fz are determined according to (96), and in a variable electromagnetic field, we must take into account that fx,fy, and fz are expended in changing the density of the momentum flux of the substance and the electromagnetic field alike. If the forces acting only on the substance are meant by f̃x,f̃y,f̃z, according to [70],  we have(102)f̃x=fx-(εμ-1)∂∂t[EH⇀]x,f̃y=fy-(εμ-1)∂∂t[EH⇀]y,f̃z=fz-(εμ-1)∂∂t[EH⇀]z.
The expression (εμ-1)(∂/∂t)[EH⇀] is known to be called the Abraham force [70].If the medium is heterogeneous and the transmission of the electromagnetic wave involves its absorption and heating, we propose the following scheme of calculating ponderomotive forces:(1)
the domain under study is subdivided into N cells with constant properties;(2)
by the method presented in Para 1 and Para 2, we compute the surface chargeσ and the surface ponderomotive force  for each cell.In this approach, we disregard the influence of striction forces and tensions acting only on the force distribution over the body’s volume but exerting no influence on the value of the resultant of all forces and on their moment [70]. Allowance for ponderomotive forces is important, for example, for biomechanics problems.The history of the problem of ponderomotive forces developed into a paradoxical situation where the issues of certain basic propositions of classical electrodynamics have been discussed to date [66].


## 5.1. Reduction of Bulk Ponderomotive Forces to Tensions

Iff is the bulk density of ponderomotive forces, the resultant of all forces applied to the bodies within the volume V is(94)F=∫VfdV.
On the other hand, if bulk forces can be reduced to tensions at all, the set of tensions acting externally of a closed surface S must be equal to the same quantity. We denote by Tn a force acting externally per unit surface the external normal to which is directed along n; the components of this force along the coordinate axes are Txn,Tyn and Tzn. Then the resultant of all tensions applied externally to the surface will be(95)F=⨖sTndS.
Using the Ostrogradsky-Gauss theorem and equating (94) and (95), we obtain(96)fx=∂Txx∂x+∂Txy∂y+∂Txz∂z,fy=∂Tyx∂x+∂Tyy∂y+∂Tyz∂z,fz=∂Tzx∂x+∂Tzy∂y+∂Tzz∂z.
For the bulk forces and tensions to be equivalent, it is necessary that not only the resultant of forces applied to an arbitrary volume remains constant but the moment of these forces be invariant as well [70]. This circumstance imposes an additional constraint on the tension-tensor components. The moment of bulk forces applied to an arbitrary volume V is(97)N=∫⌊Rf⌋dV,
where R is the distance from point O for which the moment of forces is determined to the element dV. Consequently, the component Nx must satisfy relation (97):(98)Nx=∫V(yfz-zfy)dV=⨖s(yTzn-zTyn)dS+∫V(Tyz-Tzy)dV.
The surface integral on the right is equal to the moment of tension forces Tn applied to the surface S of volume V. The moment of these tension forces will be equal to the moment of bulk forces only in the case where the last integral on the right is equal to zero. In view of the arbitrariness of the volume V, this will occur only when the integrand
(99)Tyz=Tzy
is equal to zero.By repeating the same considerations forNy and Nz, we obtain Tzx=Txy, Txy=Tyx. Consequently, for the bulk forces and surface ponderomotive forces to be equivalent, the tension tensor must be symmetric. In anisotropic media, the tension tensor is asymmetric; therefore, direct computation of the resultant of forces in terms of tensions from formula (95) is expedient. We note that Maxwell was the first to reduce the field’s ponderomotive forces to tensions. Replacement of ponderomotive forces by an equivalent system of tensions makes it much easier to determine forces applied to the volume of the electric field.The force applied to a charged layer, given the uncompensated surface charge, has the form(100)FijS=σ(Eij1-Eij2).
If the layer is not charged, we have σ=0 and FijS=0. The resultant force applied to the boundary of a cell is(101)F=⨖FijdS.
The projections of fx, fy, and fz are determined according to (96), and in a variable electromagnetic field, we must take into account that fx,fy, and fz are expended in changing the density of the momentum flux of the substance and the electromagnetic field alike. If the forces acting only on the substance are meant by f̃x,f̃y,f̃z, according to [70],  we have(102)f̃x=fx-(εμ-1)∂∂t[EH⇀]x,f̃y=fy-(εμ-1)∂∂t[EH⇀]y,f̃z=fz-(εμ-1)∂∂t[EH⇀]z.
The expression (εμ-1)(∂/∂t)[EH⇀] is known to be called the Abraham force [70].If the medium is heterogeneous and the transmission of the electromagnetic wave involves its absorption and heating, we propose the following scheme of calculating ponderomotive forces:(1)
the domain under study is subdivided into N cells with constant properties;(2)
by the method presented in Para 1 and Para 2, we compute the surface chargeσ and the surface ponderomotive force  for each cell.In this approach, we disregard the influence of striction forces and tensions acting only on the force distribution over the body’s volume but exerting no influence on the value of the resultant of all forces and on their moment [70]. Allowance for ponderomotive forces is important, for example, for biomechanics problems.The history of the problem of ponderomotive forces developed into a paradoxical situation where the issues of certain basic propositions of classical electrodynamics have been discussed to date [66].


## 6. Conclusion

The consistent physicomathematical model of propagation of an electromagnetic wave in a heterogeneous medium has been constructed using the generalized wave equation and the Dirichlet theorem. Twelve conditions at the interfaces of adjacent media were obtained and justified without using a surface charge and surface current in explicit form. The conditions are fulfilled automatically in each section of the heterogeneous medium and are conjugate, which made it possible to use through-counting schemes for calculations. For the first time the effect of concentration of “medium-frequency" waves with a length of the order of hundreds of meters at the fractures and wedges of domains of size 1–3μm has been established. Numerical calculations of the total electromagnetic energy on the wedges of domains were obtained. It is shown that the energy density in the region of wedges is maximum and in some cases may exert an influence on the motion, sinks, and the source of dislocations and vacancies and, in the final run, improve the near-surface layer of glass due to the “micromagnetoplastic" effect.The results of these calculations are of special importance for medicine, in particular, when microwaves are used in the therapy of various diseases. For a small, on the average, permissible level of electromagnetic irradiation, the concentration of electromagnetic energy in internal angular structures of a human body (cells, membranes, neurons, interlacements of vessels, etc.) is possible.For the first time, we have constructed a consistent physicomathematical model of interaction of nonstationary electric and thermal fields in a layered medium with allowance for mass transfer. The model is based on the methods of thermodynamics and on the equations of an electromagnetic field and is formulated without explicit separation of the charge carriers and the charge of an electric double-layer. We have obtained the relations for the electric-field strength and the temperature, which take into account the equality of the total currents and the energy fluxes, to describe the electric and thermal phenomena in layered media where the thickness of the electric double-layer is small compared to the dimensions of the object under study.We have modeled numerically the heating of an electrochemical cell with allowance for the influence of the electric double-layer at the metal-electrolyte interface. The calculation results are in satisfactory agreement with experimental data.



---

*Source: 185647-2012-04-08.xml*