FiPy 中一维扩散方程方程的默认边界条件实现

Question

我正在尝试使用 FiPy 求解一维传热方程。 文档指出，如果没有在外表面上指定边界条件，则默认边界条件相当于零梯度

如果我想对左脸应用零梯度，我应该消除左脸的边界条件代码吗？下面的代码代码会自动对左面应用零约束边界条件吗？

# User defined Boundary Condition for the right face
valueRight = 100
phi.constrain(valueRight, mesh.facesRight)

这是我第一次在 python 和 FiPy 工作。任何帮助将不胜感激。

Answer 1

默认边界条件为零通量。这通常会转化为零梯度，但并非必须如此，而且假设它们是同一件事可能会产生误导。

如果您指定的是

valueRight = 100
phi.constrain(valueRight, mesh.facesRight)

那么你将在左边有一个 Neumann（零通量）条件，在右边有一个 Dirichlet（固定值）条件。