Shock Fronts in Non-Polar Electro-Conducting Media

We analyze the propagation of electromagnetic fronts in unbounded electric conductors. Our model is based on the Maxwell model of electromagnetism, which includes the displacement current and Ohm’s law in its simplest forms. The shock-like electromagnetic front is a propagating surface, across which the electric and magnetic fields, as well as their higher temporal and spatial derivatives, experience finite jumps. The shock-like fronts are essentially different as compared with the weak fronts; in particular, the bulk Maxwell equations are essentially insufficient for the analysis of the shock-like fronts, and they should be amended we the physical jump conditions. We choose these additional conditions by using conditions similar to those suggested by Heaviside. We derive the basic shock intensity relationships implied by this model.


Introduction
Establishing the self-consistent mathematical model of electromagnetic waves (and of light, in particular) is the culmination of the Maxwell model of electromagnetism [1]. Maxwell himself has not seen recognition of his model by the greatest electricians of his time, including Kelvin and Helmholtz. However, forthcoming developments in science eventually proved to be instrumental in thousands of applications.
In the absence of electric and magnetic polarization, the bulk system of Maxwell equations includes (in its contemporary form) 4 partial differential equa-tions for the distributed electric current and the electric and magnetic fields.
This system is basically linear (the only source of nonlinearity can appear in the generalized Ohms law for the high amplitude currents.) Many approaches to solving Maxwell equations and the explicit exact solutions can be found in thousands of textbooks and monographs.
The most close to our approach is the approaches of the monographs of Luneburg [2] and Born and Wolfe [3]. These two monographs contain the sections dealing with the electromagnetic wave-fronts treated as surfaces of discontinuities of electromagnetic fields. For other media, similar methods were used by Hadamard [4], Levi-Civita [5], Thomas [6], and Keller [7]. In their publications, the interested readers can find further references and historic discussions.
The key instrument in analysis of wavefronts is the so-called compatibility conditions. In 20th century, the biggest contributions in their usage belong to the outstanding geometers and mathematical physicists first usage of this approach in electromagnetism belongs to Luneburg [2], as well as Born and Wolf [3].
In the recent paper Grinfeld and Grinfeld [8], we explored an evolution of the so-call wek fronts in electric conductors. What still remains unexplored, is the dynamics of the shock-like wave fronts in dissipative electro-conducting media in the framework of the full Maxwell model. In this paper, we fill this gap of the classical analysis.
Compatibility conditions have been discussed in detail in the monographs of Michael Grinfeld [9], and Pavel Grinfeld [10]. From these monographs, we borrow required results and notation. We use standard tensorial notation and operations in the form, presented in the same monographs [9] [10]. The Latin (spatial) indices , , , i j k  run the values 1, 2, 3); the Greek (surface) indices , , , α β γ  run the values 1, 2. The space is referred to the immobile coordinates i z ; j i δ is the spatial Kronecker delta; ij z and ij z are the co-and contravariant components of the metrics, which is used for raising and lowering ("juggling") spatial indices, and for definition of covariant differentiation i ∇ ; ijk z is the Levi-Civita skew-symmetric tensor.

The Bulk Master System of Equations
The bulk Maxwell equations read where c is a speed of light, i E is an electric field, i H is a magnetic field; Q and i I are the bulk free charge and electric current.
If there were non-zero charges Q at 0 t = , they charges will decay exponen-

The Heaviside Shock Conditions and Their Implications
The bulk system of Maxwell Equations (1.1) -(1.4) is sufficient when analyzing smooth solutions. However, when the electromagnetic field and its derivative experience discontinuities, Equations (1.1) -(1.4) become meaningless in some points of space-time and the description, based on these partial differential equations, appears to be essentially incomplete. This incompleteness should be compensated by some additional relationships.
The additional relationships are different for different types of singularities.
For the so-called weak discontinuities, the fields i E and i H are continuous everywhere but not their first temporal and spatial derivatives. On the contrary, the first derivatives experience finite jumps on some moving singular surfaces, called wavefronts. No additional physical ideas should be used for analysis of weak discontinuities. The physical Equations (1.1) -(1.6), in principle, fully describe the weak discontinuities. But still some additional relationships should be used at the fronts. They are called the compatibility relationships. We explored this sort of weak fronts in our earlier publication [8].
When dealing with electromagnetic shock waves, not only first derivatives but also the fields themselves experience finite jumps. Therefore, geometric and geometric boundary conditions significantly change as compared with the weak fronts, as it follows from comparison of the compatibility relations of [8], and the relations below. Moreover, for shock fronts, it is not sufficient to use only the geometric and kinematic compatibility conditions. For this sort of singular surfaces, we have to amend the bulk Maxwell equations with their analogies at the front.
There are different approaches to establishing additional conditions across the shock waves. We will dwell on the additional conditions known to Heaviside (Born and Wolfe [3]). We will call them the Heaviside conditions; they read: The Equations (2.3) -(2.6) are very important for the future and self-explanatory.

The Compatibility Conditions for Shock Fronts
In addition to the physical shock-front conditions (2.1), (2.2), the following condition is implied by the geometry and kinematics only: . .
where the zero and first order jump vectors are defined as We notice that the jumps vectors are defined on the moving shock fronts only.
Let us present the zero order jump vectors , i i a b in the following form:

Shock-Front Velocity and the Zeroth Order Jump Vector
Let us rewrite Equations (2.1), (2.2) as follows

Evolution of Shock-Wave Intensity along Rays
Calculating the jumps of the terms in the bulk master system (1.1) -(1.4) and using the compatibility conditions (3.1) -(3.3), we arrive at the following relationships: . .