The dynamics of incompressible fluids are established within the frame work of Newtonian mechanics and must satisfy the conservation laws of Newtonian mechanics: Conservation of mass, conservation of charge, and conservation of energy. These conservation laws can be derived from Reynolds’ transport theorem, which states that the rate of a physical function in the local volume $V$ of space with time is equal to the rate of change of the physical function in the volume with time plus the flux of the function through the surface $S$ of the $V$ at a certain time interval $\Delta t$ [5] . The physical function mentioned here can be any physical properties of the system, for the electromagnetic field, a composite system composed of $E$ and $B$ fields, the physical functions include electric field strength and magnetic induction intensity. The Reynolds theorem can be formulated as

$\frac{\text{d}}{\text{d}t}{\displaystyle {\oint }_{V}f\text{d}v}={\displaystyle {\oint }_V\left[\frac{\partial f}{\partial t}+\nabla \cdot \left(fv\right)\right]\text{d}v}$(2)

where $f$ is the physical function. The theorem implies that only the physical properties corresponding to the two items in the same equation can participate in the transport process at the same time and in the same place, otherwise it is impossible to participate in them.

Suppose there be only one type of energy transport in a local volume $V$ of space, $\varphi$ be the energy density (nothing except it), $\phi v$ be the energy flow density, $v$ be the transporting velocity. The total rate of change of $\varphi$ in the $V$ can be expressed by using Equation (2)

$\frac{\text{d}}{\text{d}t}{\displaystyle {\oint }_V\varphi \text{d}v}={\displaystyle {\oint }_V\left[\frac{\partial \varphi }{\partial t}+\nabla \cdot \left(\varphi v\right)\right]\text{d}v}$(3)

If there is no energy loss, the integrand on the right-hand of Equation (3) is equal to zero, and the energy is conserved:

$\frac{\partial \varphi }{\partial t}+\nabla \cdot \left(\varphi v\right)=0$(4)

This equation is in a fully Eulerian form. The energy of electromagnetic fields is composed of the E-field energy and the B-field energy. There is no electromagnetic wave without F law and M law, and these two laws run through the entire wave propagation process. From this point of view, they describe the energy transport characteristics of electromagnetic fields. The forms of F law and M law are [4]

$\nabla \times E+\frac{\text{1}}{c}\frac{\partial B}{\partial t}=0$(5)

$\nabla \times B-\frac{1}{c}\frac{\partial E}{\partial t}=0$(6)

Equations (5) and (6) are basic laws for unstable electromagnetic fields, which indicate that $E$ and $B$ are not independent, they are influence and restrict each other. Consider the electric field and magnetic field as fluid fields by Maxwell’s analogy, and notice that the materiality of the field indicates that the two formulas are certainly energy transport equations.

By definition, the direction of the curl field $\nabla \times E$ is the direction in which the circulation density is maximized, which can ensure the conservation of energy in the energy transition process [6] :

$\nabla \times E=\underset{\Delta S\to 0}{\lim }\frac{1}{\Delta S}n\cdot {\left\{{\displaystyle {\oint }_{L}E}\cdot \text{d}l\right\}}_{\max }$(7)

where $L$ is the closed path surrounding $S$, the $S$ is the surface of the volume $V$, $n$ is the unit normal vector of the $\Delta S$. The choice of $n$ ensures that the circulating current has a maximum value. Because of the conservation of electromagnetic energy, the change of energy in the volume $V$ should equal to the change of energy passing through the surface S of the $V$, namely that

${\displaystyle {\int }_S\left[\left(\nabla \times E\right)+\frac{1}{c}\frac{\partial B}{\partial t}\right]\cdot n\cdot \text{d}s=0}$(8)

There is an analogy between equations (8) and (3), the correspondences are

$\frac{\partial \varphi }{\partial t}\leftrightarrow \frac{\text{1}}{c}\frac{\partial B}{\partial t}\cdot n$, $\nabla \cdot \left(\varphi v\right)\leftrightarrow \left(\nabla \times E\right)\cdot n$(9)

Similarly, another correspondences are

$\frac{\partial \varphi }{\partial t}\leftrightarrow -\frac{\text{1}}{c}\frac{\partial E}{\partial t}\cdot n$, $\nabla \cdot \left(\varphi v\right)\leftrightarrow \left(\nabla \times B\right)\cdot n$(10)

These show that F law and M law are respectively the representations of Reynolds’ theorem for free electromagnetic fields under the condition of energy conservation.

Biot-Savart law is an experimental law, the current intensity $I_0$ is actually contributed by the negative charge of electrons, which current intensity is $I_e$, and $-I_0=I_e$, the formula is rewritten as

$\text{d}B=-\frac{{\mu }_0}{4\pi }\frac{I_e\text{d}l\times r_0}{r^2}$(11)

where $r_0=r/r$, $r$ is radius vector, ${\mu }_0$ is the permeability of vacuum, $I_e\text{d}l$ is the current-carrying element for the electrons. The mass flow intensity of electrons formed by the masses of electrons are $I_m$, which is in the same direction as $I_e$, and the resulting $W$ can be expressed as

$\text{d}W=-\frac{{\xi }_0}{4\pi }\frac{I_m\text{d}l\times r_0}{r^2}$(12)

where ${\xi }_0/4\pi$ is proportional constant. The negative sign in the right side of this formula indicates that the wagnetic field obeys the left-handed rule. On one hand, Newton’s equation of gravitational field and Coulomb’s equation of electrostatic field have similar algebraic forms (parity); on the other hand, $\text{d}W$ and $\text{d}B$ share a similar algebraic form (parity). This suggests that we can describe the properties of the gravitowagnetic field using force lines, which can be likened to an incompressible fluid. Helmholtz theorem provides a theoretical reason for the existence of wagnetic field. By vector properties, the gravitational field is a longitudinal field like electrostatic field, and wagnetic field is a transverse field like magnetic field.

Helmholtz’s theorem is only applicable vector fields in an infinite space [4] . It also indicates that: 1) a vector field in an infinite space is uniquely determined by its divergence and rotation. If both the divergence and rotation of a vector are zero, the vector field also vanish; 2) any vector field can be expressed by the sum of an irrotational field and a non-divergence field. The existence of electromagnetic waves indicates that the electromagnetic field conforms to this theorem. In recent years, gravitational waves have been detected, reminding us that the vector field excited by mass may also conform to this theorem. Introduce the scalar potential and vector potential, and represent the field excited by mas as

$G=\nabla {\phi }_m$, $W=-\nabla \times A_m$(13)

where $G$ is static gravitational field, $W$ is steady-state wagnetic field, ${\phi }_m$ is scalar potential, and $A_m$ is vector potential. For fields excited by electric charge, there are interactions between charges of the same sign as well between charges of opposite signs. For fields excited by mass, there is no interaction between negative mass and mass; there is only interaction between the masses of two objects. Therefore, the interaction between the masses of two objects is similar to the interaction two electric charges of the same sign. The difference lies in the fact that the force related to mass is gravitational, while the force related to electric charge is repulse. This means that the gravitowagnetic field follows left-handed rules.

The experimental law of magnetic induction intensity states

$\nabla \times B=\frac{4\pi }{c}{j}_0=-\frac{4\pi }{c}{j}_e$(14)

where $j_e=-j_0$, $j_0$ is current density, $j_e$ is current density contributed by electronic charges, the corresponding electronic mass flux density is $j_m$, which has the same direction as $j_e$, So

$\nabla \times W=-\frac{4\pi }{c}{j}_m$(15)

During the transition from a high energy orbit to a low energy orbit, an electron emits electromagnetic wave. Since the speed of the wave is limited, the electromagnetic wave continues to propagate after the transition is complete, and the received radiation wave can be considered as a source-free wave. After the process of gravitowagnetic waves generated in the distant depths of the universe has ended, the Earth receives this radiation waves, and therefore it is also considered a source-free waves. Since the properties of $G$ and $W$ are similar to those of $E$ and $B$, we can write similar wave equations,

${\nabla }^{2}G-\frac{1}{c^2}\frac{{\partial }^{2}G}{\partial t^2}=0$(16)

${\nabla }^{2}W-\frac{1}{c^2}\frac{{\partial }^{2}W}{\partial t^2}=0$(17)

The corresponding Poynting vector (obeys the left-handed rule) is

$S_m=-\frac{c}{4\pi }G\times W=\frac{c}{4\pi }W\times G$(18)

$S_m$ is parallel to z-axis. The plane wave solutions fro equations (16) and (17) are

$G=\left(G_{\max }\cdot j\right)\sin \left(k_z\cdot z-\omega t\right)=G_y\cdot j$(19)

$W=\left(W_{\max }\cdot i\right)\sin \left(k_z\cdot z-\omega t\right)=W_x\cdot i$(20)

where $G_{\max }=W_{\max }$, $i$ and $j$ are basis vectors of the x-axis and y-axis, $k_z$ is wave vector in the z-axis, $\omega$ is frequency, $k_z=\omega c$,

$\nabla \times G=\left(-i\right)G_{\max }\cos \left(k_z\cdot z-\omega t\right)\cdot k_z$(21)

$\frac{\partial W}{\partial t}=\left(W_{\max }\cdot i\right)\cos \left(k_z\cdot z-\omega t\right)\cdot \left(-\omega \right)$(22)

From Equations (21) and (22), we get the F law for the gravitowagnetic field

$\nabla \times G=\frac{1}{c}\frac{\partial W}{\partial t}$(23)

With the same reason, the relevant M law is

$\nabla \times W=-\frac{1}{c}\frac{\partial G}{\partial t}$(24)

In contrast to the F law and M law of electromagnetic field, the F law and M law embody compliance and resistance, respectively.

General Maxwell’s equations include electromagnetic field equations and gravitowagnetic field equations,

$\nabla \cdot E=4\pi {\rho }_{\text{0}}$, $-\nabla \cdot G=4\pi {\rho }_m$(25)

$\nabla \times E=-\frac{1}{c}\frac{\partial B}{\partial t}$, $-\nabla \times G=-\frac{1}{c}\frac{\partial W}{\partial t}$(26)

$\nabla \cdot B=0$, $\nabla \cdot W=0$(27)

$\nabla \times B=\frac{1}{c}\frac{\partial E}{\partial t}+\frac{4\pi }{c}{j}_{\text{0}}$, $-\nabla \times W=\frac{1}{c}\frac{\partial G}{\partial t}+\frac{4\pi }{c}{j}_m$(28)

The equations indicate that electromagnetic field and gravitowagnetic field maintain the parity of the excited vacuum, with no overall time reversal, meaning time maintains a forward direction. Locally, there is time reversal and parity violation, where a weak force field exists. Such states exist only momentarily and do not interfere with the overall parity and unidirectionality of time. As T. D. Lee pointed out: symmetry breaking is observable and there should exist field and force that disrupt symmetry [7] . To date, no force or field has been observed that violates the vacuum excitation-state’s parity on a large scale. Here it needs to be emphasized that parity is a property of topological space, and it’s not related to numerical value because measure is not a topological quantity, a topological space can exist without a metric.

The F law and M law should exhibit transport characteristics, from which the wave equations can be derived.

Take the time derivative of the second equation in (28), then substitute the second equation of (26) into the differential expression to eliminate $W$,

$\frac{1}{c^2}\frac{{\partial }^{2}G}{\partial t^2}=-\nabla \times \left(\nabla \times G\right)-\frac{4\pi }{c^2}\frac{\partial j_m}{\partial t}$(29)

Using the formula of vectors $\nabla \times \left(\nabla \times G\right)=\nabla \left(\nabla \cdot G\right)-{\nabla }^{2}G$, get active wave equation,

${\nabla }^{2}G-\frac{1}{c^2}\frac{{\partial }^{2}G}{\partial t^2}=4\pi \left(-\nabla {\rho }_m+\frac{1}{c^2}\frac{\partial j_m}{\partial t}\right)$(30)

The radiation wave equation is

${\nabla }^{2}G-\frac{1}{c^2}\frac{{\partial }^{2}G}{\partial t^2}=0$(31)

Similarly, for the $W$

${\nabla }^{2}W-\frac{1}{c^2}\frac{{\partial }^{2}W}{\partial t^2}=\frac{4\pi }{c}\nabla \cdot j_m$(32)

${\nabla }^{2}W-\frac{1}{c^2}\frac{{\partial }^{2}W}{\partial t^2}=0$(33)

Apart from the plane wave solutions, Equations (31) and (33) also have other forms of solutions. This indicates that the F law and M law are universally applicable, even though they were derived from the solutions for plane waves. Using the second formula of (13) and Equation (23), we get $\nabla \times \left(G+\frac{1}{c}\frac{\partial A_m}{\partial t}\right)=0$, the brackets enclose a longitudinal field, referring to the first formula of Equation (13) and noticing the parity of $G$ and $E$:

$E=-\nabla {\phi }_e-\frac{1}{c}\frac{\partial A_e}{\partial t}$, it can be represented as $\nabla {\phi }_m$, and

$G=\nabla {\phi }_m-\frac{1}{c}\frac{\partial A_m}{\partial t}$(34)

Substituting Equation (34) into the second formula of Equation (25), we obtain the wave equation of the scalar vector,

${\nabla }^2{\phi }_m-\frac{1}{c^2}\frac{{\partial }^2{\phi }_m}{\partial t^2}=-4\pi {\rho }_m$(35)

The radiation wave equation is

${\nabla }^2{\phi }_m-\frac{1}{c^2}\frac{{\partial }^2{\phi }_m}{\partial t^2}=0$(36)

On the other hand, Equation (35) can be rewritten as

$\frac{1}{c}\frac{\partial }{\partial t}\left[-\nabla \cdot A_m+\frac{1}{c}\frac{\partial {\phi }_m}{\partial t}\right]=0$(37)

The sum of the two terms in brackets of the formula is constant, let it be zero,

$-\nabla \cdot A_m+\frac{1}{c}\frac{\partial A_m{}_m}{\partial t}=0$(38)

This is the Lorenz gauge for the gravitowagnetic field. Below, we prove the rationality of the Lorenz gauge. Let the left side of Equation (38) be

${H|}_{t=0}={\left(-\nabla \cdot A_m+\frac{1}{c}\frac{\partial A_m{}_m}{\partial t}\right)}_{t=0}=0$(39)

If $H$ meets this initial condition, then at any given moment $H=0$. From the derivation of the Lorenz gauge, it is known that the initial values must satisfy the following equation

${\left[-\frac{1}{c}\frac{\partial }{\partial t}\left(\nabla \cdot A_m\right)+{\nabla }^2{\phi }_m+4\pi {\rho }_m\right]}_{t=0}=0$(40)

Combine (39), get

${\left(\nabla \cdot G+4\pi {\rho }_m\right)}_{t=0}=0$(41)

Substitute Equation (35) into (40), get

${\frac{\partial H}{\partial t}|}_{t=0}=0$(42)

Take the derivation of the wave equation of ${\phi }_m$ with respect to $t$, and then take the divergence of the wave equation of $A_m$, get

${\nabla }^{2}H-\frac{1}{c^2}\frac{{\partial }^{2}H}{\partial t^2}=\frac{4\pi }{c}\left(\nabla \cdot j_m+\frac{\partial {\rho }_m}{\partial t}\right)$(43)

The two terms within the parentheses on the right side of Equation (43) precisely form the mass continuity equation, which equals zero when mass is conservation. $H$ satisfies the homogeneous wave equation

${\nabla }^{2}H-\frac{1}{c^2}\frac{{\partial }^{2}H}{\partial t^2}=0$(44)

Combining Equations (38) and (39), it can be derived that $H=0$ at any time. Therefore, under the condition of mass conservation, the Lorenz gauge holds true.

The generalized Maxwell’s equations indicate that the wagnetic field and the magnetic field have similar properties. Therefore, we believe that there is an action of the wagnetic field on moving mass, analogous to the action of the magnetic field on moving charges. We collectively refer to this force as the Lorentz force $f_g$,

$f_g=-mv\times W=-p\times W$(45)

where $v$ is velocity, $p$ is momentum, and the negative sign indicates that the force conforms the let-handed rule.

The diameter of the Sun is about 1.39 million kilometers, which is 109 times that of the Earth’s diameter. The volume of the Sun is 1.3 million times that of the Earth, and its mass is 330,000 times that of the Earth. The Sun is a star located at the center of the solar system and is almost an ideal sphere. Due to the presence of flowing hot plasma, the wagnetic field generated by its rotation fluctuates unpredictably in space, resembling an alternating field. The photons emitted from the Sun may exhibit a precession similar to that electrons in an alternating magnetic field, under the action of the Lorentz force of the alternating wagnetic field. In addition, the influence of the Sun’s wagnetic field on Mercury’s motion cannot be ignored when calculating the position of Mercury’s perihelion.

When photons approach a massive star in the universe, they are subjected to the Lorentz force of the wagnetic field of the rotating star, causing the path of light to bend. The more higher the frequency, the more the momentum, and the more curved the path.

Consider the phenomena observed in superfluid experiments. From equation (23), F law, it is known that the wagnetic field induced by the superfluid is aligned with the Earth’s wagnetic field, and the Lorentz force $f_g$ acting on the upper layer of the superfluid molecules is as illustrated in Figure 1 , where $W_a$ is the wagnetic field of the Earth, $v$ is the velocity of the molecules in the upper fluid layer. In the experiment, on the cross-section perpendicular to $W_a$, the liquid circulate in a counterclockwise direction: It is this force that drives the superfluid molecules in the upper layer to climb up the walls of the container. As shown in Figure 1 , the climbing phenomenon exhibits a distinct directional nature, with no such effect observed on the opposite side of the container. If the force were due to molecular interactions with the container material, the wall climbing phenomenon would be no directionality. Similarly, at the bottom of the container, the superfluid molecules experience a downward Lorentz force from the Earth’s wagnetic field because their velocity direction is opposite to that of the upper layer molecules, causing them to penetrate through the bottom and exit the container. The osmosis also occurs at the side walls parallel to $W_a$, with no osmosis at the side walls perpendicular to $W_a$. In summary, the phenomena of wall climbing and penetration demonstrate directionality [8] .