python - 您如何在 Fipy 中指定 Neumann(垂直于面的固定通量)边界条件？

Question

如何在不限制面内通量分量的情况下，将 fipy 网格中边界面的法向通量明确设置为特定值？

Neumann 边界条件可以指定为:(1) 固定的垂直于边界面的通量分量，或 (2) 作为面处通量的完整规范。默认的fipy条件是前者(value = 0)，但是显式方法(faceGrad.constrain)是后者。可以通过在 fipy diffusion.mesh20x20 末尾添加以下代码来理解该问题示例并注意不同的面部梯度结果。

facesNeumann = mesh.exteriorFaces & ~facesTopLeft & ~facesBottomRight
print 'grad(phi) with default Neumann BC...'
print phi.faceGrad.value.T[facesNeumann.value]
phi.faceGrad.constrain(0.,facesNeumann)
DiffusionTerm().solve(var=phi)
print 'and with explicit Neumann BC...'
print phi.faceGrad.value.T[facesNeumann.value]

Answer 1

请参阅关于 fixed flux boundary conditions 的讨论.基本上，您将包含所需边界通量散度的源添加到 FiPy 的默认无通量条件。