我从事电磁建模，我使用有限体积法作为工具，所以我想通过改进的有限体积法（三角网格）开发代码，我研究的问题是在可渗透环境（铁磁材料）中。

首先，我尝试 进行离散化，（见图） $-\mathrm{div}(\mathrm{grad}(A)=\mu_0(J + \mathrm{rot}(M))$

我现在的问题

我搜索以确定三角形重心处的磁矢量势$A$

我的问题是如何像矩形网格一样通过 FVM 三角形离散化上面的这个方程

$$-\left(\frac{\partial^2 A}{\partial x^2} + \frac{\partial^2 A}{\partial y^2}\right) = \mu_0 \left(J_S + \frac{\partial M_y}{\partial x} - \frac{\partial M_x}{\partial y}\right)\, .$$