我是静磁学领域的新手。我想编写一个二维有限元代码来获取螺线管线圈内的磁场。我从二维码开始，并遵循了J.Bastos 和 N Sadowski书中描述的方法。我设法从FEMM（磁性和静电的开源软件包）中获得了解决方案并对其进行了验证。我还使用Nathan Ida书中描述的方法构建了二维轴对称问题的代码。但是，当我尝试使用 FEMM 验证我的代码时（因为书中没有验证问题），我无法匹配结果。

我正在尝试的轴对称问题是边 = 60mm 的正方形域。对称轴是正方形的底边。中心有一个边长为 20mm 的小方块，其电流密度为$10^6$ $A/m^2$. 外部正方形边界的所有边缘都在 A=0 处，其中 A 是磁矢量势。具有电流密度的内部正方形的材料是一种$\mu_r$= 1 并且外部正方形有$\mu_r$= 10.我正在求解更大正方形内的矢量势。我正在求解的方程是

$\vec{\nabla} \times \vec{H}=\vec{J}$

$\vec{\nabla}\cdot\vec{B}=0$

$\vec{B}=\mu \vec{H}$

我正在求解磁矢量势$A$在哪里

$\vec{\nabla} \times \vec{A}=\vec{B}$

我的问题是

1. 我已经坚持了一段时间。是否有任何论文或书籍存在验证问题，并且清楚地提出了该方法？

2. 是否有我可以参考的轴对称案例的开源代码？或者任何可以帮助我研究磁矢量势轴对称有限元算法的专家？

PS：我不确定我是否已经清楚地解释了这个问题。让我知道是否需要任何澄清