【电磁场计算】磁偶极子在外磁场中的受力

背景

最近在总结电磁场计算公式的时候发现，关于磁偶极子受力的计算公式，在不同文献中的说明有点出入，在此特地写一篇博客澄清一下。

受力推导

通过磁偶极子在外磁场中的磁势能求梯度得出其受力：
$$\begin{array}{rl}F& =-\nabla (U)\\ & =-\nabla (-m\cdot B)\\ & =m\times (\nabla \times B)+(m\cdot \nabla )B\end{array}$$
在计算的区域内（磁偶极子所在的微小空间）内可认为无外电流，即$\nabla \times B=0$，故有：
$$F=(m\cdot \nabla )B\text{\hspace{1em}(1)}$$

以上是标准的计算公式，很多教材中也是这样写的。而有的文献中直接描述为$F$等于磁矩$m$与磁场梯度$\nabla B$相乘——$F=m\cdot \nabla B$，这是不对的。

矢量展开

式(1)右侧是磁矩$m$右点乘梯度算子$\nabla$，再对磁感应强度$B$求梯度，展开为:
$$\begin{array}{rl}F& =(m_x\cdot \frac{\partial }{\partial x}+m_y\cdot \frac{\partial }{\partial y}+m_z\cdot \frac{\partial }{\partial z})B\\ & =\left[\begin{array}{c}{m}_x\frac{\partial B_x}{\partial x}+m_y\frac{\partial B_x}{\partial y}+m_z\frac{\partial B_x}{\partial z}\\ m_x\frac{\partial B_y}{\partial x}+m_y\frac{\partial B_y}{\partial y}+m_z\frac{\partial B_y}{\partial z}\\ m_x\frac{\partial B_z}{\partial x}+m_y\frac{\partial B_z}{\partial y}+m_z\frac{\partial B_z}{\partial z}\end{array}\right]\\ & =\left[\begin{array}{ccc}\frac{\partial B_x}{\partial x}& \frac{\partial B_x}{\partial y}& \frac{\partial B_x}{\partial z}\\ \frac{\partial B_y}{\partial x}& \frac{\partial B_y}{\partial y}& \frac{\partial B_y}{\partial z}\\ \frac{\partial B_z}{\partial x}& \frac{\partial B_z}{\partial y}& \frac{\partial B_z}{\partial z}\end{array}\right]\left[\begin{array}{c}{m}_x\\ m_y\\ m_z\end{array}\right]\\ & =J(B)m\end{array}$$

上式中，$J(B)$为磁场的雅可比矩阵，仔细观察，它与磁场的梯度区别如下，根据向量的梯度的定义(数值计算之 梯度向量和梯度矩阵，雅可比矩阵，海森矩阵)，有：
$$\nabla B=\left[\begin{array}{ccc}\frac{\partial B_x}{\partial x}& \frac{\partial B_y}{\partial x}& \frac{\partial B_z}{\partial x}\\ \frac{\partial B_x}{\partial y}& \frac{\partial B_y}{\partial y}& \frac{\partial B_z}{\partial y}\\ \frac{\partial B_x}{\partial z}& \frac{\partial B_y}{\partial z}& \frac{\partial B_z}{\partial z}\end{array}\right]$$
即有：${\nabla B}^{\top }=J(B)$

结论

根据以上推导，拟清楚了B的雅可比矩阵和梯度矩阵之间的区别，就可以把磁偶极子的受力写成如下简洁的表达式：
$$F=(m\cdot \nabla )B=J(B)m={\nabla B}^{\top }m=\nabla B\cdot m$$