船舶力学

海洋运动物体尾流感应电磁场计算模型

PDF下载

洪方文 ¹^,² , 石杰栋 ¹^,² , 周根水 ¹^,² , 张军 ¹^,²

船舶力学 | 流体力学 2026,30(6): 843-852

收起

船舶力学 |流体力学 2026 , 30 (6) : 843 -852

海洋运动物体尾流感应电磁场计算模型

全屏

洪方文¹^,², 石杰栋¹^,², 周根水¹^,², 张军¹^,²

作者信息

¹中国船舶科学研究中心，江苏无锡 214082

²船舶振动噪声重点实验室，江苏无锡 214082

A Wake-Induced Electromagnetic Fields Calculation Model for Marine Moving Objects

Fang-wen HONG¹^,², Jie-dong SHI¹^,², Gen-shui ZHOU¹^,², Jun ZHANG¹^,²

Affiliations

¹China Ship Scientific Research Center, Wuxi 214082, China

²National Key Laboratory on Ship Vibration & Noise, Wuxi 214082, China

HONG Fang-wen (1971−), male, Ph.D., professor, corresponding author, E-mail: hongfangwen@sina.com

出版时间: 2026-06-15 doi: 10.3969/j.issn.1007-7294.2026.06.001

文章导航

摘要

收起

为计算海洋运动物体绕流场感应的电磁场，本研究基于相对性原理推导了运动导电介质中的麦克斯韦方程组，并解决了欧姆定律的洛伦兹协变性问题。通过洛伦兹变换将问题转换至随体坐标系，尾流感应电磁场的计算简化为稳态问题。同时，本研究给出了点源与偶极子流场感应电磁场的解析表达式，证明涡量场是产生感应电荷的必要条件。建立的航行体尾流感应电磁场近似计算模型表明，在物体后方1公里处感应磁感应强度约为0.0006纳特斯拉。

关键词

尾流 / 感应电磁场 / 洛伦兹变化

Abstract

收起

To calculate the electromagnetic fields induced by the flow around moving objects in the ocean, Maxwell’s equations for a moving conductive medium were derived based on the principle of relativity, and the Lorentz covariant problem of Ohm's law was addressed. By applying the Lorentz transformation to shift the problem into a reference frame moving with the object, the calculation of electromagnetic fields induced by the wake becomes a steady-state problem, significantly simplifying the process.Formulas for electromagnetic fields induced by point sources and dipole flow fields are provided, revealing that a vortex field is a necessary condition for generating induced charge. An approximate calculation model for the induced electromagnetic fields in the wake of a marine moving object is established, with which it is estimated that the induced magnetic intensity at a distance of 1 kilometer behind the object is approximately 0.0006 nT.

Key words

wake / induced electromagnetic fields / Lorentz transformation

引用本文

洪方文, 石杰栋, 周根水, 张军. 海洋运动物体尾流感应电磁场计算模型. 船舶力学, 2026 , 30 (6) : 843 -852 . DOI: 10.3969/j.issn.1007-7294.2026.06.001

Fang-wen HONG, Jie-dong SHI, Gen-shui ZHOU, Jun ZHANG. A Wake-Induced Electromagnetic Fields Calculation Model for Marine Moving Objects[J]. Journal of Ship Mechanics, 2026 , 30 (6) : 843 -852 . DOI: 10.3969/j.issn.1007-7294.2026.06.001

正文

收起

0　Introduction

收起

Moving objects in the ocean generate large-scale and long-lasting wakes. Due to the conductivity of seawater, the motion of seawater in the wake across the Earth's magnetic field induces electromagnetic fields. These induced electromagnetic fields have become important subjects of study in marine physics^[1].

In the 19th century, while formulating the law of electromagnetic induction, Michael Faraday also investigated the relationship between seawater motion and the geomagnetic field. From the 1960s and 1970s onward, research on electromagnetic fields induced by seawater flow has been conducted from an ocean physics perspective^[2], focusing on ocean waves^[3−5], tides^[6], currents^[7−8], and internal waves^[9]. After Fraser^[3] noted that the flow induced by ship motion can also disturb the geomagnetic field, studies on calculation models^[10−12], simulation methods^[13−16], experimental technologies^[17−19], and the characteristics of wake-induced electromagnetic fields^[20−23] have become key areas in the marine science and technology.

The calculation method for the magnetic fields induced by seawater motion is based on simplified Maxwell equations. For quasi-steady flows, a steady-state approximation is used, such as for ocean circulation^[7]. For flows where time variation is not dominant, the displacement current term in Maxwell’s equations is ignored (e.g., waves and swell^[4−5]); For the wake of ships, Maxwell’s equations in the ground frame^[10−11,13] are solved, where the wake's time-dependent nature complicates the problem. Ohm's law is essential for these calculations^[23−24]. For the current form of Maxwell's equations, Ohm's law does not satisfy the covariance of coordinate transformation, and it is necessary to solve the unsteady Maxwell's equations in the ground coordinate system or to use some approximations^[12]. By adopting a frame moving with the ship, the wake field becomes time-invariant, simplifying the problem significantly.

This paper analyzes the covariance of Ohm's law and proposes an improved form of Maxwell's equations to ensure Lorentz covariance. Maxwell’s equations are then presented in the ship's moving frame, along with solution methods for steady-state wake-induced electromagnetic fields.

1　The covariance of Ohm's law

收起

According to the electromagnetic theory, the electric field strength

$ \boldsymbol{E} $

and magnetic induction strength

$ \boldsymbol{B} $

induced by seawater flow in the Earth's magnetic field satisfy Maxwell's equations. In the ground frame

$ (x,{y},{z}) $

, there are

$ \left\{\begin{aligned}& \nabla \cdot \boldsymbol{E}=\frac{{\rho }_{{\mathrm{e}}}}{{\varepsilon }_{0}}\\&\nabla \times \boldsymbol{E}=-\frac{\partial \boldsymbol{B}}{\partial t}\\&\nabla \cdot \boldsymbol{B}=0\\&\nabla \times \boldsymbol{B}={\mu }_{0}\boldsymbol{J}\mathbf+{\mu }_{0}{\varepsilon }_{0}\frac{\partial \boldsymbol{E}}{\partial t}\end{aligned}\right. $

(1)

where

$ \boldsymbol{J} $

is the electric current density,

$ {\rho }_{{\mathrm{e}}} $

is the charge density，

$ {\varepsilon }_{0} $

and

$ {\mu }_{0} $

are the permittivity and permeability of free space, respectively,

$ t $

is time. For conductive seawater, the current density relates to the electric and magnetic fields as

$ \boldsymbol{J}\mathbf=\sigma \left(\boldsymbol{E}\mathbf+\boldsymbol{u}\times \boldsymbol{B}\right) $

(2)

where,

$ \boldsymbol{u} $

is the velocity vector of seawater, and

$ \sigma $

is the conductivity of seawater. Eq.(2) represents Ohm's law.

In special relativity, Maxwell's equations exhibit Lorentz covariance. The relationship between physical quantities in a moving frame with velocity

$ ({v},0,0) $

and those in the ground frame is given by

$ \begin{aligned}{\boldsymbol{E}}^{\prime}&=\left({E}_{x},\gamma \left({E}_{y}-{v}{B}_{z}\right),\gamma \left({E}_{z}+{v}{B}_{y}\right)\right)\\{\boldsymbol{B}}^{\prime}&=\left({B}_{x},\gamma \left({B}_{y}+\frac{{v}}{{c}^{2}}{E}_{z}\right),\gamma \left({B}_{z}-\frac{{v}}{{c}^{2}}{E}_{y}\right)\right)\end{aligned} $

(3)

$ \begin{aligned}\rho _{{\mathrm{e}}}^{\prime}&=\gamma \left({\rho }_{{\mathrm{e}}}-\frac{{v}{J}_{x}}{{c}^{2}}\right)\\{\boldsymbol{J}}^{\prime}&=\left(\gamma \left({J}_{x}-{v}{\rho }_{{\mathrm{e}}}\right),{J}_{y},{J}_{z}\right)\end{aligned} $

(4)

$ {\boldsymbol{u}}^{\prime}=\frac{1}{1-\dfrac{{v}{u}_{x}}{{c}^{2}}}\left({u}_{x}-{v},{u}_{y}/\gamma ,{u}_{z}/\gamma \right) $

(5)

where the superscript apostrophe represents the physical quantity in the moving frame, and

$ \gamma =1/\sqrt{1-{{v}}^{2}/{c}^{2}} $

is the Lorentz factor, c is the speed of light.

In the moving frame, according to Lorentz covariance, Ohm's law should also hold and have the same form as Ohm's law in the ground frame, that is

$ {\boldsymbol{J}}^{\prime}\mathbf=\sigma \left({\boldsymbol{E}}^{\prime}\mathbf+{\boldsymbol{u}}^{\prime}\times {\boldsymbol{B}}^{\prime}\right) $

(6)

In addition, by substituting the inverse transformation of Eqs. (3) and (4) into Eq.(2) , the calculation formula for current density in the moving frame can also be obtained

$ {\boldsymbol{J}}^{\prime}=\frac{\sigma }{\gamma \left(1+\dfrac{{v}u_{x}^{\prime}}{{c}^{2}}\right)}\left({\boldsymbol{E}}^{\prime}+{\boldsymbol{u}}^{{{}_{'}}}\times {\boldsymbol{B}}^{\prime}+{\boldsymbol{i}}^{\prime}\frac{{v}}{{c}^{2}}{\boldsymbol{u}}^{{{}_{'}}}\cdot {\boldsymbol{E}}^{\prime}\right)-{\boldsymbol{i}}^{\prime}{v}\rho _{{\mathrm{e}}}^{\prime} $

(7)

where

$ {\boldsymbol{i}}^{\prime} $

is the unit vector of the coordinate

$ x $

components in the moving frame, which is consistent with the direction of the frame's motion velocity. Eq.(7) is significantly different from Eq.(6) , which means that Eq.(6) has no covariance and requires a new form to ensure the covariance of Ohm's law.

2　Maxwell's equations for moving conductor

收起

The issue of Ohm's law covariance arises from the incomplete electromagnetic constitutive relations for moving conductive media. By introducing two auxiliary quantities, electric displacement vector

$ \boldsymbol{D} $

and magnetic field strength

$ \boldsymbol{H} $

, the polarization and magnetization effects of electromagnetic medium can be solved^[25]. For a moving medium, the influence of motion is reflected by reconstructing the relationships among the electric displacement vector

$ \boldsymbol{D} $

, magnetic field strength

$ \boldsymbol{H} $

, polarization strength

$ \boldsymbol{P} $

, and magnetization strength

$ \boldsymbol{M} $

^[26−28]. For the case of moving conductors, there have also been studies^[29−30], but there is no complete and concise form proposed, mainly low-speed approximation results^[31−33]. In this section, Maxwell’s equations for moving conductors will be presented through the Lorentz transformation between the frame that moves with the conductive medium and the ground frame.

For a conductive medium, assuming its polarization and magnetization are both zero, at a point in space, Maxwell’s equations in the frame moving with the medium are

$ \left\{\begin{aligned}& {\nabla }^{\prime}{}^{\prime}\cdot {\boldsymbol{E}}^{\prime}{}^{\prime}=\frac{{\rho }_{0}}{{\varepsilon }_{0}}\\&{\nabla }^{''}\times {\boldsymbol{E}}^{''}=-\frac{\partial {\boldsymbol{B}}^{''}}{\partial {t}^{''}}\\&{\nabla }^{''}\cdot {\boldsymbol{B}}^{''}=0\\&{\nabla }^{''}\times {\boldsymbol{B}}^{''}={\mu }_{0}\sigma {\boldsymbol{E}}^{''}\mathbf+{\mu }_{0}{\varepsilon }_{0}\frac{\partial {\boldsymbol{E}}^{''}}{\partial {t}^{''}}\end{aligned}\right. $

(8)

where the superscript with two strokes represents the physical quantity in the frame that moves with the medium, and

$ {\rho }_{0} $

is the charge density measured in a frame that is stationary relative to the medium. Convert the first equation of Eqs. (8) to the ground frame without converting the charge density, resulting in

$ \nabla \cdot \boldsymbol{E}\mathbf+u\left(\frac{1}{{c}^{2}}\frac{\partial {E}_{x}}{\partial t}-{\left(\nabla \times \boldsymbol{B}\right)}_{x}\right)=\frac{{\rho }_{0}}{{\gamma }_{{\mathrm{u}}}{\varepsilon }_{0}} $

(9)

where

$ {\gamma }_{{\mathrm{u}}}=1/\sqrt{1-{u}^{2}/{c}^{2}} $

. Convert the fourth equation of Eq.(8) to the ground frame and write it as a component in the form of

$\left\{\begin{aligned}& \frac{u}{{c}^{2}}\nabla \cdot \boldsymbol{E}\mathbf+\left({\mu }_{0}{\varepsilon }_{0}\frac{\partial {E}_{x}}{\partial t}-{\left(\nabla \times \boldsymbol{B}\right)}_{x}\right)=-\frac{{\mu }_{0}\sigma }{{\gamma }_{{\mathrm{u}}}}{E}_{x}\\&\left({\mu }_{0}{\varepsilon }_{0}\frac{\partial {E}_{y}}{\partial t}-{\left(\nabla \times B\right)}_{y}\right)=-{\gamma }_{{\mathrm{u}}}{\mu }_{0}\sigma \left({E}_{y}-u{B}_{z}\right)\\&\left({\mu }_{0}{\varepsilon }_{0}\frac{\partial {E}_{z}}{\partial t}-{\left(\nabla \times \boldsymbol{B}\right)}_{z}\right)=-{\gamma }_{{\mathrm{u}}}{\mu }_{0}\sigma \left({E}_{z}+u{\boldsymbol{B}}_{y}\right)\end{aligned}\right. $

(10)

Multiply the first equation in Eqs. (10) by

$ u $

and subtract it from Eq.(9) to obtain

$ \nabla \cdot \boldsymbol{E}={\gamma }_{{\mathrm{u}}}\left(\frac{{\rho }_{0}}{{\varepsilon }_{0}}+{\mu }_{0}\sigma u{E}_{x}\right) $

(11)

Multiply Eq.(8) by

$ \dfrac{u}{{c}^{2}} $

and subtract the first equation in Eqs. (10) to obtain

$ {\left(\nabla \times \boldsymbol{B}\right)}_{x}={\gamma }_{{\mathrm{u}}}{\mu }_{0}\left(\sigma u{E}_{x}+u{\rho }_{0}\right)+{\mu }_{0}{\varepsilon }_{0}\frac{\partial {E}_{x}}{\partial t} $

(12)

The formula

$ {c}^{2}=1/{\mu }_{0}{\varepsilon }_{0} $

for calculating the speed of light is utilized in the above calculations. Considering

$ \boldsymbol{u}\mathbf{=}\left(u,0,0\right) $

, the following related statement extensions can be made as

$\begin{split}& u{E}_{x}\rightarrow \boldsymbol{u}\cdot \boldsymbol{E}\\& {\rho }_{0}u\rightarrow {\rho }_{0}\boldsymbol{u}\\& -u{B}_{z}\rightarrow {\left(\boldsymbol{u}\times \boldsymbol{B}\right)}_{y}\\& u{B}_{y}\rightarrow {\left(\boldsymbol{u}\times \boldsymbol{B}\right)}_{z}\end{split} $

(13)

At the same time, similar transformations and calculation methods are adopted for the middle two equations of Eq.(8) , and Maxwell’s equations for a moving conductive medium in the ground frame are obtained as follows

$\left\{\begin{aligned}& \nabla \cdot \boldsymbol{E}={\gamma }_{{\mathrm{u}}}\left(\frac{{\rho }_{0}}{{\varepsilon }_{0}}+{\mu }_{0}\sigma \boldsymbol{u}\cdot \boldsymbol{E}\right)\\&\nabla \times \boldsymbol{E}=-\frac{\partial \boldsymbol{B}}{\partial t}\\&\nabla \cdot \boldsymbol{B}=0\\&\nabla \times \boldsymbol{B}={\gamma }_{{\mathrm{u}}}{\mu }_{0}\left[\sigma \left(\boldsymbol{E}\mathbf+\boldsymbol{u}\times \boldsymbol{B}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right]+{\mu }_{0}{\varepsilon }_{0}\frac{\partial \boldsymbol{E}}{\partial t}\end{aligned}\right. $

(14)

Compared with Eq.(1) , we can obtain

$ {\rho }_{{\mathrm{e}}}={\gamma }_{{\mathrm{u}}}{\rho }_{0}+{\gamma }_{{\mathrm{u}}}{\mu }_{0}{\varepsilon }_{0}\sigma \boldsymbol{u}\cdot \boldsymbol{E} $

(15)

$ \boldsymbol{J}\mathbf={\gamma }_{{\mathrm{u}}}\left[\sigma \left(\boldsymbol{E}+\boldsymbol{u}\times \boldsymbol{B}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right] $

(16)

Eq.(15) indicates that the charge of a moving conductive medium includes the net charge in a stationary state and the charge induced by the medium's motion. Eq.(16) is the form of Ohm's law, which includes the current generated by the movement of net charges and the expansion effect of medium motion. We substitute the expressions for charge density and current density into Eqs. (15) and (16) , and also substitute the inverse transformations of Eqs. (3) and (5) into Eq.(4), to recalculate the expressions for current density and charge density in the moving frame

$ {{{\rho }^{\prime}_{{\mathrm{e}}}}}={\gamma }_{{{{\mathrm{u}}}^{\prime}}}{\rho }_{0}+{\gamma }_{{{{\mathrm{u}}}^{\prime}}}{\mu }_{0}{\varepsilon }_{0}\sigma {\boldsymbol{u}}^{\prime}\cdot {\boldsymbol{E}}^{\prime} $

(17)

$ {\boldsymbol{J}}^{\prime}\mathbf={\gamma }_{{{{\mathrm{u}}}^{\prime}}}\left[\sigma \left({\boldsymbol{E}}^{\prime}+{\boldsymbol{u}}^{\prime}\times {\boldsymbol{B}}^{\prime}\right)\mathbf+{\rho }_{0}{\boldsymbol{u}}^{\prime}\right] $

(18)

The expressions of Eqs. (17) and (18) are exactly the same as those of Eqs. (15) and (16) , indicating that the newly defined charge density and Ohm's law are covariant under Lorentz transformation, thus ensuring that the expression of the equation system has an invariant form under various inertial frames, and thereby providing a consistent mathematical and physical description of the variation of the electromagnetic fields in the ship frame proposed earlier.

3　Induced electromagnetic fields of steady-state wave and flow

收起

The induced electromagnetic fields by seawater motion, as the name suggests, originate from the movement of seawater.When the motion remains constant in an appropriate reference frame, the induced electromagnetic fields also remain constant, thereby transforming the problem into a steady-state problem. The induced electromagnetic fields from various steady marine motions—such as ship waves generated by vessels moving at constant speed, internal waves caused by submerged bodies (the so-called body effect), and slowly varying ocean currents—all exhibit this characteristic. These are collectively categorized as steady-state wave- and flow-induced electromagnetic problems. Based on Eq.(14) , the time derivative term can be eliminated by considering the realistic condition where the flow velocity is significantly lower than the speed of light (

$ {\gamma }_{{\mathrm{u}}}\cong 1 $

) . The governing equations for this class of electromagnetic fields therefore become

$\left\{\begin{aligned}& \nabla \cdot \boldsymbol{E}=\frac{{\rho }_{0}}{{\varepsilon }_{0}}+{\mu }_{0}\sigma \boldsymbol{u}\cdot \boldsymbol{E}\\&\nabla \times \boldsymbol{E}=0\\&\nabla \cdot \boldsymbol{B}=0\\&\nabla \times \boldsymbol{B}={\mu }_{0}\left[\sigma \left(\boldsymbol{E}+\boldsymbol{u}\times \boldsymbol{B}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right]\end{aligned}\right. $

(19)

Performing a divergence operation on the fourth equation in Eq.(19) leads to

$ \sigma \nabla \cdot \boldsymbol{E}\mathbf+\sigma \nabla \cdot \left(\boldsymbol{u}\times \boldsymbol{B}\right)+\nabla \cdot \left({\rho }_{0}\boldsymbol{u}\right)=0 $

(20)

Applying the vector operation formula

$ \nabla \cdot \left(\boldsymbol{u}\times \boldsymbol{B}\right)=\boldsymbol{B}\cdot \left(\nabla \times \boldsymbol{u}\right)-\boldsymbol{u}\cdot \left(\nabla \times \boldsymbol{B}\right) $

(21)

By substituting the fourth equation of Eq.(19) into Eq.(21), we can obtain

$ \nabla \cdot \left(\boldsymbol{u}\times \boldsymbol{B}\right)=\boldsymbol{B}\cdot \left(\nabla \times \boldsymbol{u}\right)-{\mu }_{0}\left[\sigma \boldsymbol{u}\cdot \boldsymbol{E}+{\rho }_{0}\boldsymbol{u}\cdot \boldsymbol{u}\right] $

(22)

Substitute the first equation of Eqs. (19) and (22) into Eq.(20), and use

$ {\gamma }_\text{u}\cong 1 $

to obtain

$ \frac{{\rho }_{0}}{{\varepsilon }_{0}}+\frac{1}{\sigma }\nabla \cdot \left({\rho }_{0}\boldsymbol{u}\right)=-\boldsymbol{B}\cdot \left(\nabla \times \boldsymbol{u}\right) $

(23)

Due to

$ {\varepsilon }_{0}=8.85\times {10}^{-12}\;{\mathrm{F/m}} $

，

$ \sigma =5\;{\mathrm{s/m}} $

，and the

$ \boldsymbol{u} $

is generally on the order of m/s, the above equation can be approximately expressed as

$ {\rho }_{0}\cong -{\varepsilon }_{0}\boldsymbol{B}\cdot \left(\nabla \times \boldsymbol{u}\right) $

(24)

As shown in Eq.(24) , the presence of vorticity in the oceanic velocity field generates induced charges, while irrotational motion yields zero charge density when measured in the medium's reference frame.

Furthermore, from Eq.(19) , it can be seen that the curl of the electric field

$ \boldsymbol{E} $

and the divergence of the magnetic field

$ \boldsymbol{B} $

are equal to 0, which can be expressed as

$ \boldsymbol{E}\mathbf=\nabla \phi $

(25)

$ \boldsymbol{B}\mathbf=\nabla \times \boldsymbol{A } $

(26)

where,

$ \phi $

is a scalar function,

$ \boldsymbol{A } $

is a vector function, and

$ \nabla \cdot \boldsymbol{A}=0 $

, in this way

$\left\{\begin{aligned}& \Delta \phi =\frac{{\rho }_{0}}{{\varepsilon }_{0}}+{\mu }_{0}\sigma \boldsymbol{u}\cdot \boldsymbol{E}\\&\Delta \boldsymbol{A}=-{\mu }_{0}\left[\sigma \left(\boldsymbol{E}+\boldsymbol{u}\times \boldsymbol{B}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right]\end{aligned}\right. $

(27)

The electric and magnetic fields in the right-hand expression of Eq.(27) can be decomposed into their background and induced components as follows

$ \boldsymbol{E}\mathbf={\boldsymbol{E}}_{{\mathrm{e}}}+{\boldsymbol{E}}_{{\mathrm{u}}} $

(28)

$ \boldsymbol{B}\mathbf={\boldsymbol{B}}_{{\mathrm{e}}}+{\boldsymbol{B}}_{{\mathrm{u}}} $

(29)

where,

$ {\boldsymbol{E}}_{{\mathrm{e}}} $

and

$ {\boldsymbol{B}}_{{\mathrm{e}}} $

are the electric field and magnetic induction intensity of the Earth in the moving frame, both being spatially uniform and temporally constant within the studied domain.

$ {\boldsymbol{E}}_{{\mathrm{u}}} $

and

$ {\boldsymbol{B}}_{{\mathrm{u}}} $

are the electric field and magnetic induction intensity induced by seawater motion, generally much smaller than the electric field and magnetic field of the Earth. Therefore, Eqs. (24) and (27) can be approximated as

$ \left\{\begin{aligned}& {\rho }_{0}=-{\varepsilon }_{0}{\boldsymbol{B}}_{{\mathrm{e}}}\cdot \left(\nabla \times \boldsymbol{u}\right)\\&\Delta \phi =\frac{{\rho }_{0}}{{\varepsilon }_{0}}+{\mu }_{0}\sigma \boldsymbol{u}\cdot {\boldsymbol{E}}_{{\mathrm{e}}}\\& \Delta \boldsymbol{A}=-{\mu }_{0}\left[\sigma \left({\boldsymbol{E}}_{{\mathrm{e}}}+\boldsymbol{u}\times {\boldsymbol{B}}_{{\mathrm{e}}}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right]\end{aligned}\right. $

(30)

Eq.(30) constitute an approximate governing system for steady-state wave- and flow-induced electromagnetic phenomena, forming a well-posed boundary value problem when supplemented with appropriate electromagnetic boundary conditions. For infinite space, the solutions of the steady-state wave- and flow-induced electric field and magnetic induction intensity obtained using scalar and vector Poisson equations are

$ \left\{\begin{aligned}& {\rho }_{0}=-{\varepsilon }_{0}{\boldsymbol{B}}_\text{e}\cdot \left(\nabla \times \boldsymbol{u}\right)\\&{\boldsymbol{E}}_\text{u}=-\nabla {\iiint }_{v}\frac{\left({\rho }_{0}/{\varepsilon }_{0}+{\mu }_{0}\sigma \boldsymbol{u}\cdot {\boldsymbol{E}}_\text{e}\right)}{4\text{π} r}{\mathrm{d}}V\\&{\boldsymbol{B}}_\text{u}=\nabla \times {\iiint }_{v}\frac{{\mu }_{0}\left[\sigma \left({\boldsymbol{E}}_\text{e}+\boldsymbol{u}\times {\boldsymbol{B}}_\text{e}\right)\mathbf+{\rho }_{0}\boldsymbol{u}\right]}{4\text{π} r}{\mathrm{d}}V\end{aligned}\right. $

(31)

The third equation in Eq.(31) is the Biot Savart formula, where

$ \boldsymbol{r} $

is the vector from the field point to the spatial integration point, and

$ r $

is its magnitude.

4　Induced electromagnetic source and sink flow fields

收起

4.1　Induced electromagnetic a source flow field

The velocity field of the unit point source located at

$ {\boldsymbol{r}}_{0} $

can be expressed as

$ \boldsymbol{u}\left({\boldsymbol{r}}^{\prime}\right)=-\nabla \frac{1}{4\text{π} \left| {\boldsymbol{r}}^{\prime}-{\boldsymbol{r}}_{0}\right| } $

(32)

According to Eq.(31) , the induced electromagnetic source field is

$\left\{\begin{aligned}& {\rho }_{0}=0\\&{\boldsymbol{E}}_\text{u}=\frac{{\mu }_{0}\sigma }{16{\text{π} }^{2}}\nabla \left({\boldsymbol{E}}_\text{e}\cdot \nabla {\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime}-{\boldsymbol{r}}_{0}\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }{\mathrm{d}}{V}^{\prime}\right)\\&{\boldsymbol{B}}_\text{u}=-\frac{{\mu }_{0}\sigma }{16{\text{π} }^{2}}\nabla \times \left({\boldsymbol{B}}_\text{e}\times \nabla {\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime}-{\boldsymbol{r}}_{0}\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }{\mathrm{d}}{V}^{\prime}\right)\end{aligned}\right.$

(33)

The kernel integral of Eq.(33) is

$ {\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime}-{\boldsymbol{r}}_{0}\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }{\mathrm{d}}{V}^{\prime}={\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime\prime}+\boldsymbol{D}\right| }\frac{1}{\left| {\boldsymbol{r}}^{\prime\prime}-\boldsymbol{D}\right| }{\mathrm{d}}{V}^{\prime\prime} $

(34)

where

$ \boldsymbol{D}\mathbf{=}\dfrac{\boldsymbol{r}\mathbf{-}{\boldsymbol{r}}_{0}}{2} $

$ {\boldsymbol{r}}^{\prime\prime}\mathbf{=}{\boldsymbol{r}}^{\prime}\mathbf{-}\dfrac{\boldsymbol{r}\mathbf{+}{\boldsymbol{r}}_{0}}{2} $

, equivalent to translating the coordinate origin to

$ \dfrac{\boldsymbol{r}\mathbf{+}{\boldsymbol{r}}_{0}}{2} $

. Define the angle between

$ {\boldsymbol{r}}^{\prime\prime} $

and

$ \boldsymbol{D} $

$ \phi $

, then

${\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime\prime}+\boldsymbol{D}\right| }\frac{1}{\left| {\boldsymbol{r}}^{\prime\prime}-\boldsymbol{D}\right| }{\mathrm{d}}{V}^{\prime\prime}=\int \limits_{0}^{\mathrm{\infty }}\int \limits_{0}^{\text{π} }\frac{-2\text{π} {{{r}^{\prime\prime}}}^{2}{\mathrm{d}}{r}^{\prime\prime}{\mathrm{d}}{\mathrm{cos}}\phi }{\sqrt{{\left({{{r}^{\prime\prime}}}^{2}+{D}^{2}\right)}^{2}-{\left(2{r}^{\prime\prime}D\right)}^{2}{{\mathrm{cos}}}^{2}\phi }}=D\int \limits_{0}^{\mathrm{\infty }}2\text{π} \eta {{\mathrm{sin}}}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta $

(35)

where

$ \eta ={r}^{\prime\prime}/D $

. When

$ \eta \rightarrow \mathrm{\infty } $

, the integral of Eq.(35) tends towards infinity, which does not converge, can be written as a two-stage integral as follows

$ D\int \limits_{0}^{\mathrm{\infty }}2\text{π} \eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta=D\int \limits_{0}^{\Lambda }2\text{π} \eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta+D\int \limits_{\Lambda }^{\mathrm{\infty }}2\text{π} \eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta $

(36)

where

$ \Lambda $

is a large but finite value, and the second term on the right side of Eq.(36) is approximately

$ D\int \limits_{\Lambda }^{\mathrm{\infty }}2\text{π} \eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta =4\text{π} D\int \limits_{\Lambda }^{\mathrm{\infty }}{\mathrm{d}}\eta =4\text{π} \int \limits_{\Lambda D}^{\mathrm{\infty }}{\mathrm{d}}{r}^{\prime\prime} $

(37)

let

$ \mathit{\Pi }=\int \limits_{0}^{\Lambda }2\text{π} \eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}{\mathrm{d}}\eta $

(38)

The value of Eq.(37) is infinite, but

$ \mathit{\Pi } $

is a finite value. The integral of Eq.(34) is written as follows

$ {\iiint }_{v}\frac{1}{\left| {\boldsymbol{r}}^{\prime}-{\boldsymbol{r}}_{0}\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }{\mathrm{d}}{V}^{\prime}=D\mathit{\Pi }+4\text{π} \int \limits_{\Lambda D}^{\mathrm{\infty }}{\mathrm{d}}{r}^{\prime\prime} $

(39)

Substituting Eq.(39) into Eq.(33) yields the induced electromagnetic fields of the point source induced flow field, but its value is infinite and has no practical significance.

4.2　Induced electromagnetic dipole flow field

Arrange the source and sink of dipole at positions

$ \pm d $

on the x-axis, respectively. According to the calculation formula for the induced electromagnetic fields of a point source induced flow field in the previous section, a kernel integral similar to Eq.(33) can be obtained as follows

$ {\iiint }_{v}\left(\frac{1}{\left| {\boldsymbol{r}}^{\prime}-\boldsymbol{r}_{0}^+\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }-\frac{1}{\left| {r}^{\prime}-\boldsymbol{r}_{0}^-\right| }\frac{1}{\left| \boldsymbol{r}-{\boldsymbol{r}}^{\prime}\right| }\right){\mathrm{d}}{V}^{\prime}=\left({D}^+-{D}^-\right)\mathit{\Pi }-4\text{π} \int \limits_{\Lambda {D}^-}^{\Lambda {D}^+}{\mathrm{d}}{\boldsymbol{r}}^{\prime\prime}=\left({D}^+-{D}^-\right)\left(\mathit{\Pi }-4\text{π} \Lambda \right) $

(40)

where

$ \boldsymbol{r}_{0}^{+} $

is the location

$ \left(-d,0,0\right) $

of point source,

$ \boldsymbol{r}_{0}^{-} $

is the location

$ \left(d,0,0\right) $

of point sink, and

$ {D}^{\pm }\mathbf{=}\left| \dfrac{\boldsymbol{r}\mathbf{-}\boldsymbol{r}_{0}^{\pm }}{2}\right| $

. Using Eq.(38) , Eq.(40) can be expressed as

$ \left({D}^+-{D}^-\right)\left(\mathit{\Pi }-4\text{π} \Lambda \right)=\left({D}^+-{D}^-\right)\int \limits_{0}^{\Lambda }2\text{π} \left(\eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}-2\right){\mathrm{d}}\eta $

(41)

By using numerical integration, we can obtain

$ \int \limits_{0}^{\Lambda }2\text{π} \left(\eta {\sin}^{-1}\frac{2\eta }{{\eta }^{2}+1}-2\right){\mathrm{d}}\eta =-12.56 $

(42)

This leads to

(43)

The induced electromagnetic fields of the flow field induced by the point source and point sink arranged at position

$ \pm d $

on the x-axis are as follows

$\left\{\begin{aligned}& {\boldsymbol{E}}_\text{u}=\frac{-12.56{\mu }_{0}\sigma }{16{\text{π} }^{2}}\nabla \left[{\boldsymbol{E}}_\text{e}\cdot \nabla \left({D}^+-{D}^-\right)\right]\\&{\boldsymbol{B}}_\text{u}=\frac{12.56{\mu }_{0}\sigma }{16{\text{π} }^{2}}\nabla \times \left[{\boldsymbol{B}}_\text{e}\times \nabla \left({D}^+-{D}^-\right)\right]\end{aligned}\right.$

(44)

4.3　The electromagnetic field induced by wake of a marine moving object

This section demonstrates a simple application of the derived theory: calculating the magnetic fields induced by the flow field of a moving vehicle with constant speed in an unbounded fluid domain. The flow field is expressed using potential of source and sink. A ground-fixed reference frame (x, y, z) is established with the positive x-axis pointing towards the east, the positive y-axis pointing towards the north, the positive z-axis pointing upwards. The length, diameter, and volume of the vehicle are

$ L=100\; {\mathrm{m}} $

$ D=10\; {\mathrm{m}} $

, and

$ V=5000\;{{\mathrm{m}}}^{3} $

, respectively. The speed is

$ U=10\; {\mathrm{m/s}} $

, and its direction is (−1, 0, 0), moving from east to west. The geomagnetic intensity is

$ {\boldsymbol{B}}_\text{e} $

$ 5\times {10}^{-5}\;{\mathrm{T}} $

, which is parallel to the horizontal plane and points to the north (0, 1, 0) . These are shown in Fig. 1.

In the geodetic frame, the flow field induced by the moving vehicle varies with time, and the problem studied has unsteady characteristics, requiring the solution of the unsteady Maxwell’s equations. If the reference frame is established on the moving vehicle, the flow field around it will not change over time, and the problem will become a steady-state problem, requiring only the solution of the stationary Maxwell’s equations. The frame moving with the vehicle is

$ \left({x}^{\prime},\;{y}^{\prime},\;{z}^{\prime}\right) $

, whose coordinate axes are aligned with the geodetic frame, but the origin is located at the midpoint of the vehicle. In frame

$ \left({x}^{\prime},\;{y}^{\prime},\;{z}^{\prime}\right) $

, the vehicle is fixed, but there is an incoming flow of

$ \boldsymbol{U}\mathbf{=} $

(10, 0, 0) from far ahead. According to the electromagnetic fields transformation relationship between reference frames, under the approximation of

$ {\gamma }_\text{u}\cong 1 $

, the Earth's magnetic field remains unchanged

$ \boldsymbol{B}_\text{e}^{\prime}={\boldsymbol{B}}_\text{e} $

, and there is a uniform electric field with a direction

$ {z}^{\prime} $

in space, and its magnitude is

$ U{B}_\text{e} $

, namely

$ \left(0,0,U{B}_\text{e}\right) $

For the simplification of the problem, without considering the viscous effect, the induced flow field of the vehicle is expressed using a source and a sink with a distance of L and an intensity of

$ UV/L $

. The calculation formula for the velocity field is^[23]

$ {\boldsymbol{u}}^{\prime}\mathbf{=-}\frac{UV}{4\text{π} L}{\nabla }^{\prime}\left(\frac{1}{\sqrt{{\left({x}^{\prime}+L/2\right)}^{2}+{{{y}^{\prime}}}^{2}+{{{z}^{\prime}}}^{2}}}-\frac{1}{\sqrt{{\left({x}^{\prime}-L/2\right)}^{2}+{{{y}^{\prime}}}^{2}+{{{z}^{\prime}}}^{2}}}\right) $

(45)

The superscripts in the above equation represent the physical quantities in the reference frame that move with the vehicle,

$ {\boldsymbol{i}}^{\prime} $

$ {\boldsymbol{j}}^{\prime} $

$ {\boldsymbol{k}}^{\prime} $

are the unit directions vectors of

$ {x}^{\prime} $

$ {y}^{\prime} $

$ {z}^{\prime} $

respectively. To simplify the calculation of the electromagnetic field strength at the far rear of the vehicle, only the points on the axis are considered, namely

$ {y}^{\prime}={z}^{\prime}=0 $

. Using Eq.(44) , we obtain

$ {\rho }_{0}=0 $

(46)

$ {{{\boldsymbol{E}}^{\prime}_\text{u}}}=\frac{12.56{\mu }_{0}\sigma {U}^{2}V{B}_\text{e}}{32{\text{π} }^{2}\left({{{x}^{\prime}}}^{2}-{L}^{2}/4\right)}{\boldsymbol{k}}^{\prime} $

(47)

$ {{{\boldsymbol{B}}^{\prime}_\text{u}}}=\frac{-12.56{\mu }_{0}\sigma UV{B}_\text{e}}{32{\text{π} }^{2}\left({{{x}^{\prime}}}^{2}-{L}^{2}/4\right)}{\boldsymbol{j}}^{\prime} $

(48)

Using the inverse transformation of Eq.(3) , considering that

$ U \lt < c $

, transform Eqs. (47) - (48) into the ground frame, we have

$ {\boldsymbol{E}}_\text{u}=0 $

(49)

$ {\boldsymbol{B}}_\text{u}=\frac{-12.56{\mu }_{0}\sigma UV{B}_\text{e}}{32{\text{π} }^{2}\left({\left(x-Ut\right)}^{2}-{L}^{2}/4\right)}\boldsymbol{j} $

(50)

From Eq.(50) , the variation of magnetic induction intensity with distance in the wake of a sailing vehicle is shown in Fig. 2. At a distance of 1 kilometer behind the vehicle, the magnetic induction intensity is about 0.0006 nT, which can be detected by the most sensitive magnetic detector currently available.

5　Concluding remarks

收起

The electromagnetic induction calculation and characteristics of marine moving object’s wake are complex interdisciplinary problems of electromagnetics and fluid mechanics. In this article, the Lorentz covariance of Ohm's law is analyzed and presents Maxwell's equations in moving conductors through the principle of relativity, thereby solving the Lorentz covariance problem of Ohm's law. By using Lorentz transformation, the problem of the induced magnetic fields in the wake of a vehicle is transformed into a reference frame that follows the vehicle's motion and the unsteady problem is transformed into a steady-state problem, greatly simplifying the solving process. At the same time, an approximate calculation method for the induced magnetic fields of vehicle’s wake is provided.

However, the rationality of several issues in the theoretical analysis of electromagnetic fields in moving conductors still needs further research: 1) the relationship between Ohm's law and polarization and magnetization is unclear, and whether Ohm's law includes polarization and magnetization needs further analysis; 2) the relationship between the evolution equation of charge and the velocity of conductor motion is unclear. The use of moving coordinates seems to indicate that the evolution of charge is independent of the velocity of the medium. However, further research is needed to investigate what happens when there is a spatial gradient in the velocity of the medium.

基金

收起

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Wei B J, Chen T, Liu J, et al. Characteristics of ship’s wake induction magnetic field under background of waves[J]. Journal of China University of Petroleum (Edition of Natural Science), 2024, 48(1): 104-114.

Huang B, Liu Z Y, Xu Y J, et al. Characteristics of magnetic fields induced by the wake of an underwater vehicle[J]. Applied Sciences, 2022, 12(16): 7964.

Fraser D C. Magnetic field of ocean waves[J]. Nature, 1965, 206: 605–606.

Crews A, Futterman J. Geomagnetic micro pulsations due to the motion of ocean waves[J]. Journal of Geophysical Research, 1962, 67(5): 299-306.

Weaver J T. Magnetic variations associated with ocean waves and swell[J]. Journal of Geophysical Research, 1965, 70(8): 1921-1929.

Zhang Z L, Wei W B, Liu B H, et al. Theoretical calculation of electromagnetic field generated by ocean waves[J]. Acta Oceanologica Sinica, 2008, 30(1): 42-46.

Larsen J C. Electric and magnetic fields induced by deep sea tide[J]. Geophysical Journal International, 1968, 16(1): 47-70.

Sanford T B. Motional induced electric and magnetic fields in the sea[J]. Journal of Geophysical Research, 1971, 76(15): 3476-3492.

Liu J, Wang J T, Wang Z M, et al. Comparasion on induced electromagnetic fields generated by ocea current and tidal wave[J]. Journal of China University of Petroleum(Edition of Natural Science), 2024, 46(2): 52-62.

Chave A D. On the electromagnetic fields induced by oceanic internal waves[J]. Journal of Geophysical Research Oceans, 1984, 89(C6): 10519-10528.

Madurasinghe D, Tuck E O. The induced electromagnetic fields associated with submerged moving bodies in an unstratified conducting fluid[J]. IEEE Journal of Oceanic Engineering, 1994, 19(6): 193-199.

Madurasinghe D. Induced electromagnetic fields associated with large ship wakes[J]. Wave Motion, 1994, 20(3): 283-292.

Liu J B, Zheng X L, Du X Z. Computing method of magnetic an0maly induced by underwater vehicle’ motion in geomagnetic field[J]. Journal of Naval Academy of Engineering, 1994, 69(4): 40-45.

Yaakobi O, Zilman G, Miloh T. Detection of the electromagnetic field induced by the wake of a ship moving in a random sea of finite depth[J]. Journal of Engineering Mathematics, 2011, 70(3): 17-27.

Zhu X, Xia M. Magnetic field induced by wake of moving body in wind waves[J]. Progress in Electromagnetics Research, 2014, 149: 109-118.

Zhang J X, Xiong L, Jing R X. Modeling and simulation of electromagnetic field induced by wake of a submerged vehicle moving in shallow sea[J]. Systems Engineering and Electronics, 2016, 38(5): 1004-1009.

Wang X K, Li Y G, Li J K, et al. Numerical simulation of electromagnetic fields induced by ship wakes in ocean waves of finite depth water[J]. Feriodical of Ocean University Of China (Edition of Natural Science), 2020, 50(12): 107-114.

Zou N, Nehorai A. Detection of ship wakes using an airborne magnetic transducer[J]. IEEE Transactions on Geoscience & Remote Sensing, 2000, 38(1): 532-539.

Hirota M, Furuse T, Ebana K. Magnetic detection of a surface ship by an airborne LTS SQUID MAD[J]. IEEE Transactions on Applied Superconductivity, 2001, 11(1): 884-887.

Huang C G, Luo S R, Wu W F. Design scheme of magnetometer based on detecting electromagnetic field characteristics induced by submarine wake[C]//Proceedings of the 21^st Comprehensive Academic Symposium on Marine Surveying and Mapping, Tianjin, 2009.

Fallah M A, Abiri H. Electromagnetic fields induced by the motion of di-hull bodies in a conducting fluid[J]. IEEE Transactions on Magnetics, 2013, 49(10): 5257-5263.

Zhang J W, Jiang R X, Gong S G. Study of the electric field induced by the wake of a moving ship[J]. Journal of Harbin Engineering University, 2014, 35(8): 931-935.

Xu Z H, Du C P, Xia M Y. Modeling of electromagnetic fields induced by moving seawater due to a undersea vehicle[C]//IEEE International Conference on Computational Electromagnetic (ICCEM), 2017.

Huang B, Liu Z Y, Xu Y J, et al. Far-field characteristics and evolutions of electromagnetic field induced by the wake of underwater vehicles[J]. Applied Ocean Research, 2024, 145: 103933.

Guo S H. Electrodynamics[M]. Peking: People’s Education Press, 1979.

Penfield P, Hermann M. Electrodynamics of moving media[M]. Massachusetts: The M. I. T. Press, 1967.

Wang W Y, Xu Y. Lorentz covariant form of electromagnetism in the moving media[J]. Journal of Physics and Engineering, 2018, 28(2): 13-25.

Wang Z L, Shao J J. Maxwell’s equations for a mechano-driven varying-speed motion media system under slow motion and nonrelativistic approximations[J]. Sci Sin Tech, 2022, 52(8): 1198-1211 .

Boffi L V. Electrodynamics of moving media[D]. Massachusetts Institute of Technology, 1957.

Landau L D, Lifshitz E M. Electrodynamics of continuous media[J]. New York: Pergamon Press, 1984.

Wilhelm H E. Galilei covariant electrodynamics of moving media with applications to the experiments of Fizeau and Hoek[J]. Apeiron: Studies in Infinite Nature, 1993, 15: 1-6.

Rozov A. Maxwell equations for slow-moving media[J]. Z. Naturforsch, 2015, 70(12)a: 1019–1024.

Chen J W. A new model of covariant electrodynamics of moving media based on an alternative coordinate transformation relation instead of the usual Lorentz transformation[J]. The European Physical Journal Puls, 2022, 137: 920-930.

2026年第30卷第6期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2026.06.001

接收时间：2025-12-11
首发时间：2026-11-03
出版时间：2026-06-15

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-12-11

基金

Supported by the National Natural Science Foundation of China (51779232)

作者信息

¹中国船舶科学研究中心，江苏无锡 214082

²船舶振动噪声重点实验室，江苏无锡 214082

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2026.06.001

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT