Mechanics of electromagnetic interactions
Valery P. Dmitriyev Lomonosov University P.O.Box 160, Moscow 117574, Russia∗ (Dated: November 13, 2018)
We consider an elastic-plastic medium whose motion equations are isomorphic to Maxwell’s equations. Electrical charges are modeled by pressure centers of the medium. The electric interaction is shown to be concerned with the conservation law in the torsion ﬁeld of the medium. The Lorentz force may correspond to the Coriolis driving force due to the entrainment of the pressure center in the medium’s ﬂow.

arXiv:physics/0612210v4 [physics.gen-ph] 28 Oct 2009

I. SUBSTRATUM FOR PHYSICS

The hypothetical substratum underlying classical electromagnetism was shown [1] to be an elastic-plastic jelly-like medium. The peculiar feature of the jelly is that it is stiﬀ to compression yet liable to shear deformations. Its motion is governed by the following linear equations. The continuity equation for the incompressible medium is

∇·u=0

(1)

where u(x, t) is the velocity ﬁeld of the medium. Dynamics of the jelly is described in the Lagrange representation by the type of the Lame equation

ς∂tu + ςc2∇ × (∇ × s) + ∇p − ςf = 0

(2)

where s(x, t) is the small displacement of the element of the medium,

u = ∂ts,

(3)

ς the density of the medium, c the rate of the transverse wave spreading in the medium and f is an external force. Expression (2) diﬀers from the standard Lame equation by that for the incompressible medium we have the pressure p term instead of ∇ · s. In the medium under consideration the pressure gradient, i.e. the longitudinal stress, is created by a stationary inclusion that is described by the static external force. While the transverse stress given by the solenoidal term c2∇ × (∇ × s) is of a transitory character. In order to describe this feature we must reformulate (2) into the kind of the Prandtl-Reuss model of the plastic-elastic medium which speciﬁes the joint evolution of the velocity u(x, t) and external force −R(x, t) ﬁelds:

ς∂tu + ∇p + ςR = 0,

(4)

∂tR − c2∇ × (∇ × u) + J = 0

(5)

where J(x, t) is a current for the external force. By (5) u is a small quantity though according to (3) the displacement s may be large. This expresses the material ﬂow in the elastic-plastic medium. The medium in question is realized in the ideal turbulent ﬂuid with the elasticity c2∇ × (∇ × u) and volume force R arising due to Reynolds stresses [2]. Now we have three vector equations, (1), (4) and (5), for three vector and one scalar unknown ﬁelds, u, R, J and p. The set of equations (1), (4) and (5) will be closed if we know the dependence of J on u and p. In the case of the electromagnetic substratum such a closure can be performed in the following way.
The inhomogeneous ﬁeld R(x, t) is generated in the medium by the point defect that is described by the singularity

∇ · R = 4πbδ(x − ξ)

(6)

where b is the strength of the defect and ξ its location [1]. We suppose that the discontinuity (6) can move freely in the medium its motion velocity being deﬁned as

v

=

dξ dt

.

(7)

∗Electronic address: aether@yandex.ru

2

The motion of defects underlies the phenomenology of the plasticity. A moving with the velocity (7) singularity (6) forms in the medium the current of the defect’s strength [1]:

J = 4πbvδ(x − ξ).

(8)

We assign to the defect (6) the mass m and assume that it moves in the medium as a classical particle:

m

dv dt

=

F.

(9)

Now the problem of the closure for the set of equations (1), (4)-(9) will be to ﬁnd the dependence of F on b and the ﬁelds u and p, and also on v. In physical terms we are searching for the law that defects interact with each other by means of the ﬁelds that they generate in the medium.

II. INTERACTION DUE TO CONSERVATION LAW IN THE TORSION FIELD

Take notice that multiplying the current (8) by any vector ﬁeld R(x, t) and integrating all over the volume we will obtain the expression

J · Rd3x = 4πbv · R(ξ, t)

(10)

that depends only on coordinates of the defect (6) referred here as the particle. The form (10) can be interpreted as the rate of the work done on the particle

κ2 4π

J · Rd3x = v · F(ξ, t)

(11)

if we will assume that the force acted on the particle is

F = κ2bR

(12)

where κ2 is a positive constant. Basing himself on this observation we will multiply (5) by the term R of the external ﬁeld and integrate it over the
volume:

1 2

∂t

R2d3x − c2

R · [∇ × (∇ × u)]d3x +

J · Rd3x = 0.

(13)

Using (4) in the second term of (13) and integrating it by parts gives

−

R · [∇ × (∇ × u)]d3x =

(∂tu

+

∇p)

·

[∇

×

(∇

×

u)]d3x

=

1 2

∂t

(∇ × u)2d3x

(14)

where vector relation B · (∇ × A) = A · (∇ × B) + ∇ · (A × B) and ∂tu × (∇ × u) → 0 at inﬁnity were used. Substituting (11) and (14) into (13):

κ2 8π

∂t

[R2 + c2(∇ × u)2]d3x + F · v = 0.

(15)

Using (9) in (15) gives

d dt

κ2 8π

[R2

+

c2(∇

×

u)2]d3x

+

mv2 2

= 0.

(16)

So, with a reasonable assumption (12) for the force on the particle we have balance (16) of the particle’s kinetic energy within the integral in the torsion ﬁeld of the medium. The terms with R2 and c2(∇ × u)2 in (16) can be interpreted
as the energy of respective ﬁelds concerned with the energy of the particle. However, as we see, they have no relation to the energy of the substratum which is proportional to u2 and an integral of R over a space coordinate. Using (4)
in (12) gives for the ﬁeld force

F = −κ2b(∂tu + ∇p/ς).

(17)

3

III. DRIVING FORCE DUE TO ENTRAINMENT OF THE PARTICLE IN THE FLOW OF THE MEDIUM

We will assume that the singularity (6) is entrained in the ﬂow of the medium. So that the full velocity of the defect can be represented as the sum of the defect’s own velocity and the medium’s velocity. In a more general case the defect may be a composite structure. We will assume that it moves as a classical particle of the mass m with µ being its entrained part, m ≥ µ. Then the full velocity of the particle v can be represented as the sum of the particle’s own velocity υ and the part of the medium’s velocity

v

=

υ

+

µ u
m

(18)

where µ is the mass coeﬃcient of the pressure center (6). Suppose that the kinematics of the ﬁeld u(x, t) is known. Then we may thus or otherwise compute the driving force on the particle. In the absence of a speciﬁc force it is formally given by

F

=

m

dv dt

.

(19)

As can be shown (see Appendix) the driving force on the particle consists of the inertial force plus the force due to rotation of the particle’s velocity by the medium’s ﬂow:

F

=

µ

du dt

−

µυ

×

(∇

×

u).

(20)

In order to interpret terms in the right-hand part of (20) we will apply (20) to the well-known particular case, the rotational ﬁeld

u

=

dx dt

=

ω

×

x,

(21)

We have for (20) with (21) in two dimensions

du dt

=

−ω2x,

(22)

∇ × u = 2ω.

(23)

Substituting (23) into the last term of (20) we obtain the expression for the Coriolis driving force on the particle. Relation (22) used in (20) gives the centripetal driving force on the particle.
In the linear-elastic medium the full time derivative in (20) can be changed to partial derivative and the full particle’s velocity v substituted for its relative velocity υ:

F = µ∂tu − µv × (∇ × u).

(24)

With some suggestions concerning the mass µ of the pressure center (6) relation (24) provides an addition to the speciﬁc force in (17) found from ﬁeld equations.

IV. ENTRAINMENT IN VARIOUS SYSTEMS

We will consider three variants for the entrainment of a particle in a medium’s ﬂow. Let an impurity be advected by the laminar ﬂow of an ideal ﬂuid according to (18). Using in (20) the Euler equation

ς

du dt

+

∇p

=

0

(25)

we obtain

F

=

−

µ ς

∇p

−

µυ

×

(∇

×

u).

(26)

Relations (22)-(23) provide a particular case to (26).

4

Next we consider transport in the elastic-plastic medium. Using (4) in (24) we obtain

F

=

−µ(

1 ς

∇p

+

R)

−

µv

×

(∇

×

u).

(27)

We see that there may be a regime when the external force fully equilibrates the pressure in (27). The result will be a purely rotation of the particle’s velocity.
We have a special case in the electromagnetic substratum. Combining (24) with (17) gives with the account of (4) for the full force on the singularity

F = κ2bR + κ2bv × (∇ × u)

(28)

where the mass coeﬃcient µ in (24) is related with the strength b of the pressure center (6) by

µ = −κ2b.

(29)

As in contrast with (26) and (25) we may have in (28) and (4) separately the Coriolis force and the central force. This situation becomes possible due to the excluding in (25) of non-linear terms when passing to Lagrange coordinates. The addition of the Coriolis driving force to (17) does not change the result of the last term in (15).

V. ELECTROMAGNETIC ANALOGY

We deﬁne vector and scalar electromagnetic potentials

A = κcu,

(30)

ςϕ = κ(p − p0),

(31)

E = κR

(32)

where p0 is the background pressure. Substituting (30)-(32) into (4) we will obtain the Maxwell’s equation

1 c

∂tA

+

E

+

∇ϕ

=

0.

(33)

With the substitution (30) continuity equation (1) takes the form of the Coulomb gauge

∇ · A = 0.

(34)

We will deﬁne the electric charge by

q = κb.

(35)

Substituting (32) and (35) into (6) we obtain the source equation

∇ · E = 4πqδ(x − ξ).

(36)

Substituting (30), (32) and (35) into (5) with (8) gives another Maxwell’s equation

∂tE − c∇ × (∇ × A) + 4πqvδ(x − ξ) = 0.

(37)

Substituting (32), (30) and (35) into (28) we come to the expression for the force acted on the electric charge

v

F = q[E + c × (∇ × A)].

(38)

The proton q > 0 is modeled mesoscopically by a cavity stabilized via the perturbation of the ﬁeld R [2]. Through (6) and (4) with (1), the decreased pressure in the center of the defect corresponds to b < 0 . Hence by (35) κ < 0. Comparing (35) with (29) we may conclude that

κ

=

−

µ e

(39)

where µ is the mass coeﬃcient and e the charge of the electron. By (39) with κ < 0 we have µ < 0 for e < 0. This means that the electron will move counter to the aether ﬂow. The
loop on a vortex ﬁlament microscopic mechanism that underlies the upstream motion of the particle was considered in my earlier publication [3].

5

VI. REALIZATION

The medium that provides substratum for physics can be realized in the ideal turbulent ﬂuid [2]. Following Osborne Reynolds we expand the ﬂuid velocity into the sum of the averaged u over a short time interval and ﬂuctuation u′ components:

u = u + u′

(40)

where

u′ = 0.

(41)

Substituting (40) into (25), averaging with the account of (41) and linearizing gives [4]

ς∂t u + ς∂j u′u′j + ∇ p = 0

(42)

where ∂iu′i = 0 was used. The summation over recurrent index is implied here and throughout. Comparing (42) with (4) we see that all the ﬁelds mentioned can be related to averaged perturbations of the ideal turbulence, and the external force R can be rendered in terms of Reynolds stresses:

R = ∂j u′u′j .

(43)

Using (43) in (32) we ﬁnd that the electric ﬁeld is concerned with the turbulence force acted on the element of the medium:

E = κ∂j u′u′j

(44)

where with the account of (39) κ < 0. The meaning of (44) is that the electromotive force is proportional to the diﬀerence of the turbulence energy density u′iu′i /2 across two open-circuited terminals. Alike in (44), deﬁnitions of vector (30) and scalar (31) potentials should be changed to respective averaged values. Using (43) in the ﬁrst term of (28) we obtain for the speciﬁc force acted on the particle

Fi = −|κ|q∂j u′iu′j

(45)

where (35) and κ < 0 were used. As was above said, q > 0 corresponds to a cavity stabilized in the disturbed turbulence. According to (45) we may expect that the cavity will move in the region of the decreased turbulence intensity. This tendency just corresponds to the phenomenon in macroscopic transport known as turbophoresis. When q > 0 the dependence (45) coincides with the expression for the force on a large particle derived theoretically in [5]. The largeness of the pressure center (6) consists in that it generates in the medium the voluminous perturbation ﬁeld −b/r. The inversion of the turbophoretic force with the change of the sign of the pressure center suggests that the turbophoresis may have the same vortex loop microscopic mechanism as the Lorentz force [3]
In order to reproduce properly the Lorenz force in (38) we must assume in (18) that the pressure center is entrained by the averaged ﬂow, i.e.

v

=

υ

+

µ m

u

,

(46)

and then average the force (24) acted on the particle.

APPENDIX A: THE ENTRAINMENT FORCE

Let the particle at ξ is driven by the ﬂow of a medium. The kinematics of the medium is speciﬁed by the velocity ﬁeld u(x, t). We will assume that at any instant of time the full velocity v of the particle can be represented as the sum of the particle’s own velocity (relative velocity) υ and the medium’s velocity:

v = υ + u(ξ, t).

(A1)

The kinematics of the particle can be calculated if we know the dependence on time of the particle’s relative velocity in the reference frame xˇ accompanying the medium’s velocity ﬁeld at ξ:

υˇ = ddξtˇ.

(A2)

6

We will assume that the relative motion of the particle is governed by a speciﬁc force ℜ:

dυˇ dt

=

ℜˇ .

(A3)

Considering Lagrange dynamics in the non-inertial reference frame xˇ we will compute the driving force F − ℜ on

the particle due to the entrainment (A1). The kinetic energy of the particle in the reference frame xˇ is

where

Kˇ

=

1 2

m[υˇ(t)

+

uˇ(ξˇ, t)]2

uˇ(ξˇ(t), t) = u(ξˇ(t), t) − u(xˇ).

(A4) (A5)

Here we do not indicate explicitly the time dependence in the velocity ﬁeld u(x) of the medium since this term will

not be diﬀerentiated with respect to time. Lagrange equations for the evolution of the particle’s relative velocity (A2)

are

d ∂Kˇ dt ∂υˇ

=

∂Kˇ ∂ξˇ

+ Fˇ.

(A6)

Substituting (A4) in (A6) we have with the account of (A2)

d ∂Kˇ dt ∂υˇ

=

m

d dt

(υˇ

+

uˇ)

=

m

dυˇ dt

+

υˇ · ∇ˇ

uˇ

+

duˇ dt

,

(A7)

where and

duˇ dt

=

∂uˇ ∂t

+

uˇ · ∇ˇ

uˇ,

∂Kˇ ∂ξˇ

=

1 2

m∇ˇ

(υˇ

+

uˇ)2

=

m

υˇ · ∇ˇ

uˇ

+

mυˇ

×

(∇ˇ

×

uˇ)

+

1 2

m∇ˇ uˇ2,

(A8) (A9)

where the vector relation ∇ (a · b) = (a · ∇) b + (b · ∇) a + a × (∇ × b) + b × (∇ × a) was used. Substituting (A7) and (A9) into (A6) we obtain

m

dυˇ dt

+

m

duˇ dt

=

mυˇ

×

(∇ˇ

×

uˇ)

+

1 2

m∇ˇ uˇ2

+

Fˇ .

Using (A3) in (A10) gives for the driving force of the entrainment

(A10)

Fˇ

−

ℜˇ

=

m

duˇ dt

−

1 2

m∇ˇ uˇ2

−

mυˇ

×

(∇ˇ

×

uˇ).

(A11)

Next we may convert (A11) into the inertial reference frame ξˇ → ξ and xˇ → x taking advantage of that ∇ˇ → ∇, Fˇ → F and υˇ → υ. Substituting (A5) into (A11) gives

F

=

ℜ

+

m

du(ξ(t), dt

t)

−

1 2

m∇[u(ξ(t),

t)

−

u(x)]2

−

mυ

×

[∇

×

u(ξ(t),

t)].

(A12)

Taking x = ξ the gradient of the quadratic term in (A12) will be nulliﬁed, and we shall have ﬁnally for the force on the particle entrained in the medium

F

=

ℜ

+

m

du dt

−

mυ

×

(∇

×

u).

(A13)

Specifying dynamical law for ℜ and u in the right-hand part of (A13) we shall ﬁnd the force exerted on a particle by

a particular medium.

[1] V.P.Dmitriyev, Elasticity and electromagnetism, Meccanica 39, No 6, 511-520 (2004) [2] V.P.Dmitriyev, Mechanical realization of the medium imitating electromagnetic ﬁelds and particles, Meccanica 42, No 3,
283-291 (2007) [3] V.P.Dmitriyev, Mechanical model of the Lorentz force and Coulomb interaction, CEJP 6, No 3, 711-716 (2008) [4] O.V.Troshkin, On wave properties of an incompressible turbulent ﬂuid, Physics A 168, No 2, 881-899 (1990) [5] M.W.Reeks, The transport of descrete particles in inhomogeneous turbulence, J. Aerosol Science 14, No 6, 729-739 (1983)

