各向异性介质电流模拟问询：基于有限体积法的石墨类金属恒定电流建模

嘿，刚好之前做过各向异性导体的恒定电流模拟，结合有限体积法（FVM）给你理清楚思路，应该能解决你的问题：

咱们先从基础定律推导，核心逻辑和各向同性是通的，只是把标量电导率换成张量形式：

1. 恒定电流场的电流连续性方程：$\nabla \cdot \mathbf{J} = 0$，这个不管介质是否各向异性都成立，意思是控制体内没有电荷积累。
2. 各向异性的欧姆定律：不再是标量的 $J = \sigma E$，而是张量形式 $\mathbf{J} = \boldsymbol{\sigma} \cdot \mathbf{E}$。对于石墨这类层状材料，只要把坐标系对齐它的各向异性轴（比如层内x/y方向，层间z方向），电导率张量可以简化为对角矩阵：

   \boldsymbol{\sigma} = \begin{pmatrix}
   \sigma_{xx} & 0 & 0 \\
   0 & \sigma_{yy} & 0 \\
   0 & 0 & \sigma_{zz}
   \end{pmatrix}
   

   石墨的层内电导率（$\sigma_{xx}/\sigma_{yy}$）远大于层间（$\sigma_{zz}$），数值上可能差3-5个数量级，这个后续求解要注意。
3. 电场和电势的关系还是 $\mathbf{E} = -\nabla \phi$，把上面两个式子代入连续性方程，就能得到各向异性下的广义拉普拉斯方程：
   $$\nabla \cdot (\boldsymbol{\sigma} \cdot \nabla \phi) = 0$$
   对比各向同性的 $\sigma \nabla^2 \phi = 0$，这个方程里电导率是和梯度做张量点乘后再取散度，是模拟的核心。

FVM的核心是对控制方程做积分，避免微分的数值误差，针对这个广义拉普拉斯方程，步骤如下：

1. 控制体积分与高斯定理转换：对任意控制体$V$积分方程：
   $$\int_V \nabla \cdot (\boldsymbol{\sigma} \cdot \nabla \phi) dV = 0$$
   用高斯定理把体积分转成控制体表面$S$的面积分：
   $$\oint_S (\boldsymbol{\sigma} \cdot \nabla \phi) \cdot \mathbf{n} dS = 0$$
   这里$\mathbf{n}$是控制体面的外法向量，这个式子的物理意义是：流入控制体的总电流等于流出的总电流，符合恒定电流的要求。
2. 面通量的离散计算：
   • 如果你的网格坐标系和材料各向异性轴对齐，那每个面的通量计算很直接。比如二维中，垂直x方向的面（法向量$\mathbf{n}=(1,0)$），电流通量就是 $J_n = \sigma_{xx} \cdot \frac{\partial \phi}{\partial x}$；垂直y方向的面则是 $J_n = \sigma_{yy} \cdot \frac{\partial \phi}{\partial y}$。
   • 如果坐标系没对齐，需要先把电导率张量做坐标变换，转换到网格局部方向上，再计算通量，这个稍微麻烦一点，但原理一样。
3. 组装线性方程组：用相邻控制体的电势近似梯度（比如 $\frac{\partial \phi}{\partial x} \approx \frac{\phi_{邻体} - \phi_{本控制体}}{\Delta x}$），代入面积分后，每个控制体都会得到一个关于电势的线性方程，最后把所有方程组装成大型线性方程组，用求解器（比如CG、LU分解）求解即可。

你提到的特定边界设置固定电势（模拟电接触），这个处理逻辑和各向同性基本一致：

• 对于固定电势的边界控制体，直接把它的电势设为给定值$\phi_0$，在组装方程组时，把对应行的系数矩阵设为1（对应自身电势的项），右端项设为$\phi_0$，其他项设为0。
• 非固定电势的边界如果是绝缘的，那电流通量为0，这个在FVM里会自然处理——因为边界外没有控制体，对应的面通量项直接设为0即可。

• 刚性方程组的求解：因为石墨层内和层间电导率差异极大，组装出来的线性方程组会很“刚性”，普通的求解器可能收敛很慢甚至不收敛。建议用带预条件子的共轭梯度法（PCG），或者如果网格规模不大，直接用LU分解求解。
• 网格分辨率：层间方向电导率小，电流在层间的变化可能更剧烈，建议在层间方向适当加密网格，保证数值精度。

内容的提问来源于stack exchange，提问作者lmr