A New Formulation of Maxwell’s Equations

: In this paper, new forms of Maxwell’s equations in vector and scalar variants are presented. The new forms are based on the use of Gauss’s theorem for magnetic induction and electrical induction. The equations are formulated in both differential and integral forms. In particular, the new forms of the equations relate to the non-stationary expressions and their integral identities. The indicated methodology enables a thorough analysis of non-stationary boundary conditions on the behavior of electromagnetic ﬁelds in multiple continuous regions. It can be used both for qualitative analysis and in numerical methods (control volume method) and optimization. The last Section introduces an application to equations of magnetic ﬂuid in both differential and integral forms.


Introduction
Magnetic and magnetorheological fluids have been increasingly used in recent years [1][2][3].Their properties, especially viscosity, can be significantly changed by the effects of the magnetic field [4][5][6][7][8][9][10].The mathematical model is, in such cases, composed of Maxwell's and Navier-Stokes equations [11], which are mostly solved by numerical methods [12][13][14][15].One of these models is presented at the end of this paper.The solution behavior of the mentioned models can be assessed on the basis of the symmetry or asymmetry of the operators forming the mathematical model [16,17].Let us write, for example, Maxwell's equations in the following operator form: From here, the symmetry of the problem for a nonconductive environment, where j = 0, ρ = 0, is quite obvious.Furthermore, it follows that the non-stationary field is rotational, because it is not possible to set E = gradφ.It is therefore necessary to use new methods for non-stationary problems solution.The operator equations, written in the summation symbolics, also show that The expression .
While the Levi-Civit tensor ε ijk is antisymmetric, only the antisymmetric part E kj is used in Expression (3).Thus: (curlE These modifications could be deepened by a more detailed study of the symmetry of the mathematical model operators depending on the boundary conditions.A new formulation of Maxwell's equations, which is essentially based on Gauss's divergence theorem [17], can also contribute to this analysis. According to the authors, the mathematical formulation of Gauss's divergence theorem is underestimated in terms of the physical substantiality of the problem.Gauss's divergence theorem explains the following important finding: every spatial change in the physical variable f(x, t), regardless of its tensor character, has its response at the boundary of a closed region: where f can be understood as a symbolic variable.Of course, the term (5) has the opposite meaning.If we act on the boundary of the region S with the effect of the variable f, this variable changes within the volume V.In brief, every change within the volume V can be measured at the boundary S, where it is reflected in its functional value; see Figure 1.This principle can be used in qualitative analysis, numerical methods, and optimization of the electromagnetic field.The main reason is that this principle gives a very good idea, both qualitative and quantitative, of the influence of boundary conditions on the problem, which assignment is often in hands of the researcher.
After all, in electromagnetic fields solution, the engineering practice successfully uses integral identities both on the principle of Gaussian divergence theorem for the closed area S and Stokes theorem for the open area Ŝ bounded by a curve k [15][16][17][18][19].
For example: and others [20].Since this paper is focused on solving non-stationary problems, it is appropriate to modify Maxwell's equations into a more suitable form for optimization, qualitative analysis, and numerical methods.This can be achieved by using Gauss's divergence theorem and the symmetry of the Kronecker delta operator δ ij = δ ji .The adjustment applies to all non-stationary terms of Maxwell's equations for non-conductive media and ∂B ∂t for conductive media, where j = 0, ρ = 0.The modification, the proof of which is given in Section 2, can be expressed by a non-stationary term in a more suitable form for the application of Gauss's divergence theorem: The scalar variant of Maxwell's equations (Section 4) is also presented in the work, where new functions are introduced: Modified intensity of imprinted forces: Modified stress of printed forces: Based on Gauss's divergence theorem, they are again modified to a form suitable for analysis: Based on (7), by integrating (6) over the domain V surrounded by the surface S, new integral identities can be found, this time over the closed surface, where all boundary conditions appear.
The proof is given in Section 2. It follows from the above-mentioned equations that the article focuses on the appropriate modification of Maxwell's equations so that the influence of non-stationary terms in the field V is expressed by their values on the boundary S of the closed region.The solution is based on the use of symmetry conditions and Gauss's divergence theorem.The use of this procedure for qualitative analysis, numerical methods, and optimization are presented in the individual Sections for both the vector variant and the scalar variant of Maxwell's equations.Both conductive and non-conductive environments are considered in the solution.The last Section presents a mathematical model of the interaction of a magnetorheological fluid with a magnetic field.Even for this interdisciplinary problem, Gauss's divergence theorem can be used to redefine the mathematical model of Navier-Stokes equations.An example is given in Section 6.
A special part is devoted to the finite volume method for non-stationary problems [13][14][15][16][17][18][19][20][21][22].In the classical method, a non-stationary term is identified through the control volume in terms of the mean values of the integral calculus.This method does not allow the use of the finite volume method, while the new variant will allow it; see Section 2.

Symmetry in Principles of the Solution
In the technical sciences, symmetry conditions play a special role, especially in the stability conditions of the system [13,14,23].They have the same importance in the solution of electromagnetism tasks, where we can find in the term ∂B i ∂x j and ∂D i ∂x j .These terms conclude the symmetric S ij and antisymmetric A ij part.It can be written: Based on the symmetry principles, it is possible to find new shapes of Maxwell's equations using the Gauss divergence theorem.The principle can be explained, for example, on magnetic induction.Let us consider the following equations: The same formulated in the index symbolic (Einstein summation symbolics): Einstein's summation symbolic is used in the mentioned relation and the following text.We note that these relations depend on the antisymmetric operator ε ijk and the expression ∂E k ∂x j , which can be decomposed into symmetric and antisymmetric parts (see Section 1).Therefore, in the left part of Equation ( 14), only the antisymmetric part of the expression ∂E k ∂x j manifests.
After these remarks, let us proceed to a derivation of a new variant of the equation ∂t by modifying its right part containing a non-stationary term.We start from the validity of the equation ∂B i ∂x i = 0.For this purpose let us put: The relationship uses per partes integration in 3D space.After using a Kronecker delta symmetry adjustment: Thus, B j It is possible to write Equation (15) in the shape: Because it also holds that ∂ ∂x j ∂B j ∂t = 0, it can be derived by analogy: Expressions (19) and (20) are very important, because they point out the fact that unsteady states of the magnetic flux density are generated on the borders of the area and vice versa.From the Expression (17) follows the next important result based on the divergence theorem: Using Expression (21), it is possible to correct the form of Equation ( 14) on the principle of symmetry.After the substitution of the above-mentioned criteria, we obtain a new shape of Maxwell equations.
In Equation ( 22), the first term in bracelets represents the antisymmetric tensor of second grade, and in the second term, it is possible to decompose into the symmetric and antisymmetric parts.
Equation ( 19) can be used to modify the control volume method.
In the current control volumes method, the integration of the non-stationary term given in Maxwell's equations is expressed on the basis of the mean value of the integral calculus [5,[11][12][13][14]20,22,24].Thus, in the form (see Figure 1): where In the newly proposed method, the integration is performed directly by using Relation (21) as follows: The same in the vector form: After the integration over the volume V, (22) can be easily written in the vector variant: Note that from the obtained results (20), it is observable that the non-stationary change in magnetic induction within the field V can be determined by integration only at the system boundary.Conversely, time changes in the magnetic induction in the field V are generated at the system boundary.This fact can be advantageously used for both qualitative analysis of non-stationary boundary conditions and optimization of non-stationary tasks by selecting a suitable target function depending on the boundary conditions.For the non-conductive space, the derivations are described in Section 3.

A Non-Conductive Environment
We assume σ = 0, ρ e = 0.In this case, Maxwell's equations will have the index symbolic form [1,3]: Due to the validity of Equations ( 17) and (18) using transform (21), new forms of Maxwell's equations can be written without evidence, since: From the above, it is visible that in the case of a non-conductive environment, it is possible to derive a new variant for all Maxwell's equations.All the conclusions given in Section 3 remain valid, including the control volumes method.
All these results can be used to solve the interdisciplinary problem of the motion of an incompressible fluid with the effects of a nonconductive magnetic field using the Maxwell stress tensor for this case; see Section 6.

The Scalar Variants of Maxwell's Equations
By the scalar variant of Maxwell's equations [1,3], we mean the product of the multiplication of Equation ( 13) and the position vector x.
The results of the solution of the scalar variant can be used again for the boundary conditions analysis and in the optimization area.For the solution, the divergence theorem ( 13) is beneficially used: By analogy to (15), we apply the multiplication: From (37) follows the important knowledge: In the vector form written as: Considering the divergence theorem in the shape div ∂B ∂t = 0, Equations ( 38) and (39) can be written: V ∂B ∂t where From Equation (40) and using the divergence theorem, it is possible to derive the following important dependence that allows one to reformulate Maxwell's Equation (35): When we implement (42) into (35), we obtain: The term in the bracelet can again be divided into the symmetric and antisymmetric parts as well as in (22).In the integral form: For the vector form, it holds: One of the scalar variants of Maxwell's equations was again derived under the assumption of Gauss's divergence theorem validity.Derived relationships can be used to evaluate the results obtained by numerical methods.Even in this case, non-stationary changes in magnetic induction are reflected at the boundary of the region, and, here, it is possible to determine their values as a function of time.The resulting equations can also be easily used for optimization because the target function is scalar in this case.

A Non-Conductive Environment-Scalar Variant
If we return to the problem of a non-conductive environment, we can rewrite Equations ( 28) and (29) in a differential form: Original equations: New variant: Integral form: Original equations: New variant: The new form of Maxwell's equations is useful, whether for analysis or numerical solution.Comparing the left sides of Equations ( 49) and (50), it is obvious that nonstationary variables D(x, t), B(x, t), are generated only at the boundary of the system, and, therefore, it is possible to influence their process within the volume V.This can be essential in optimizing the non-stationary problems of electromagnetism.
Here, scalar variants of Maxwell's equations were also derived under the assumption of Gauss's divergence theorem validity.Derived relations can be used to evaluate the results obtained by numerical methods.Even in this case, non-stationary changes in magnetic induction and electrical induction are reflected at the boundary of the region, and here, their values can be determined as a function of time.The resulting equations can also be easily used for optimization because the target function is scalar in this case.

An Interaction of a Non-Conductive Magnetic Liquid with a Magnetic Field
In Section 3, it is shown that in the case of a non-conductive environment, it is possible to derive a new variant for all Maxwell's equations.All the conclusions given in the Section remain valid, including the control volumes method.All these results can be used in solving the interdisciplinary problem of the motion of an incompressible fluid with the effects of a non-conductive magnetic field using the Maxwell stress tensor for this case [3].In the presented case, we assume: The density of the volumetric magnetic force that acts on the elementary volume can be written in the form [8]: Navier-Stokes equations of the magnetic liquid in the presented case are in the form [2,11,[25][26][27]: Considering then, Equation (53) can be written in a more transparent form: Because the liquid is considered to be incompressible, the continuity equation is in the form: Now, if we consider Equation ( 21), the Navier-Stokes equation for the incompressible magnetic liquid can be written in a new form: This equation can be, using the Divergence theorem [16,17], rewritten in the new integral form: By comparing the original equation, Equation (53), and the new equation, Equation (57), the advantage of the new variant is evident both for the numerical solution by the finite volume method and for the analysis of the influence of boundary conditions.Since in the above-mentioned case, assuming diva = 0, with respect to (20), it holds: and concurrently for the result of the continuity equation:

Discussion
A new formulation of Maxwell's equations was derived, both in differential and integral variants.The basis for the derivation was Gauss's divergence theorem, used for magnetic flux density B and electric flux density D. By the use of Gauss's divergence theorem, Maxwell's equations were transformed.This resulted in a tool that can be used in the numerical finite volume method and optimization.The obtained equations will also allow the qualitative analysis of the influence of boundary conditions.The mentioned changes concern the non-stationary terms of the type ∂B ∂t , resp.∂D ∂t .This resulted in a new form of Maxwell's equations that can be used in solving the interdisciplinary problem of the motion of an incompressible fluid with the effects of a non-conductive magnetic field using the Maxwell stress tensor for this case.For example:

Differential forms
Original Maxwell s equation:

Integral forms
Original variant: In these relations, it is interesting that the effect of non-stationary members ∂B ∂t within the region V is reflected at the system boundary only by its normal component ∂B ∂t •n .

Conclusions
The work was focused on the analysis of non-stationary Maxwell equations.A new shape of non-stationary magnetic flux density was derived.This made the analyses of Maxwell's equations possible by using the Gaussian divergence theorem.Maxwell's equations were defined in both vector and scalar variants.The new shape of the Maxwell equations simplifies the analyses of the solution quality depending on the boundary conditions, considering the non-stationary magnetic induction.It also allows the numerical solution of Maxwell's equations to be extended to the large control volume method.Using the Gaussian divergence theorem, the new method allows the region to be optimized depending on the non-stationary field of magnetic induction.
A special part was devoted to the finite volume method for non-stationary problems.In the classical method, a non-stationary term is identified through the control volume in terms of the mean values of the integral calculus.This method does not allow the use of large control volumes, while the new variant allows it.
Both conductive and non-conductive environments were considered in the solution.The last Section presents a mathematical model of the interaction of a magnetorheological fluid with a magnetic field.Even for this interdisciplinary problem, Gauss's divergence theorem can be used to redefine the mathematical model of Navier-Stokes equations.

Figure 1 .
Figure 1.General volume of the liquid V surrounded by continuous closed area S that consists of open areas Ŝ (outlined by closed curve k).Volume V can be split into control volumes ∆V, closed by control areas ∆S.
x 1 , x 2 , x 3 ) spatial vector n = (n 1 , n 2 , n 3 ) unit normal vector E fluid velocity v = (v 1 , v 2 , v 3 ) g gravity acceleration y•z = y i z iscalar product of two vectors y, b