有限差分法：Dirichlet与Neumann边界处的求解问题咨询

我来梳理下你用有限差分法处理非齐次热方程的思路，尤其是边界条件的实现细节：

网格与边界条件定义 #

你的求解网格采用了混合边界规则：

• 标记为N的边界采用齐次Neumann边界条件（法向导数为0）
• 标记为D的边界采用非齐次Dirichlet边界条件（直接给定边界上的温度值）

齐次Neumann边界条件的虚点法处理 #

为了处理Neumann边界，你选择的虚点法是数值计算里非常经典的方案——把网格向外扩展一层"虚点"，通过中心差分近似一阶导数来推导虚点和实际网格点的关系，从而将边界条件融入离散方程。

以左侧边界为例，齐次Neumann条件要求边界处的法向导数为0：

$\frac{\partial T}{\partial x}\bigg|_{x=0} = 0$

假设网格步长为$\Delta x$，左侧边界的实际网格点是$T_{0,j}$，向外扩展的虚点为$T_{-1,j}$，相邻的内部网格点是$T_{1,j}$，用中心差分近似上述导数：

$\frac{T_{1,j} - T_{-1,j}}{2\Delta x} = 0$

直接整理就能得到虚点与内部点的等价关系：
$T_{-1,j} = T_{1,j}$

把这个关系代入边界点$T_{0,j}$的有限差分方程中，就能消去虚点变量，得到只包含实际网格点的求解方程，这样就顺利把Neumann边界条件整合到整个离散系统里了。

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