Poisson方程五点差分格式例题及解答

写在前面

本文针对偏微分方程数值解中出现的一道例题进行分析，详细介绍了五点差分格式的公式推导及应用。

例题

在单位正方形$\overline{\Omega }:0⩽x⩽1,0⩽y⩽1$上，
$$\{\begin{array}{rl}& \Delta u=-16,(x,y)\in \Omega \\ & u{|}_{x=1}=0,\frac{\partial u}{\partial y}{|}_{y=1}=-u\\ & \frac{\partial u}{\partial x}{|}_{x=0}=\frac{\partial u}{\partial y}{|}_{y=0}=0\end{array}$$

在$\overline{\Omega }$上布置以$h=\frac{1}{4}$为步长的正方形网格，写出五点差分格式解上述问题的方程组。

解题过程

由于其步长为$h=\frac{1}{4}$，再根据边值条件，可知
$x=1$这条轴上的五个点其值均为$0$，故可得图如下，一共有20个未知量，其中有9个正则内点和11个边界点。

首先计算内点：
对于9个正则内点$u_{i,j},i=1,2,3;j=1,2,3$，由五点差分格式公式可以得到：
$$u_{i+1,j}+u_{i,j+1}+u_{i,j-1}+u_{i-1,j}-4u_{i,j}=-16\times \frac{1}{4^2}=-1$$
即：

$$4u_{i,j}-(u_{i+1,j}+u_{i,j+1}+u_{i,j-1}+u_{i-1,j})=1,i,j=1,2,3\text{\hspace{1em}(1)}$$

再计算边界点：

①对于下边界${\frac{\partial u}{\partial y}|}_{y=0}=0$，用一阶中心差商代替一阶偏导数，得到
$$\frac{u_{i,1}-u_{i,-1}}{2h}=0,$$
即$u_{i,1}=u_{i,-1}$，又因为$j=0$时，
$$4u_{i,0}-(u_{i+1,0}+u_{i,1}+u_{i,-1}+u_{i-1,0})=1,$$
代入$u_{i,1}=u_{i,-1}$得到
$$2u_{i,0}-\left(\frac{1}{2}{u}_{i+1,0}+u_{i,1}+\frac{1}{2}{u}_{i-1,0}\right)=\frac{1}{2},i=1,2,3\text{\hspace{1em}(2)}$$

②对左边界${\frac{\partial u}{\partial x}|}_{x=0}=0$，同理得到：
$$2u_{0,j}-\left(\frac{1}{2}{u}_{0,j+1}+u_{1,j}+\frac{1}{2}{u}_{0,j-1}\right)=\frac{1}{2},j=1,2,3\text{\hspace{1em}(3)}$$

③对上边界${\frac{\partial u}{\partial y}|}_{y=1}=-u$，用一阶中心差商代替一阶偏导数，得到$u_{i,5}=u_{i,3}-\frac{1}{2}{u}_{i,4}$，所以得到：
$$\frac{9}{4}{u}_{i,4}-\left(\frac{1}{2}{u}_{i+1,4}+u_{i,3}+\frac{1}{2}{u}_{i-1,4}\right)=\frac{1}{2},i=1,2,3\text{\hspace{1em}(4)}$$

④对点$u_{0,0}=(0,0)$，根据下边界和左边界的式子，可以分别得到：
$$2u_{0,0}-\left(\frac{1}{2}{u}_{1,0}+u_{0,1}+\frac{1}{2}{u}_{-1,0}\right)=\frac{1}{2},$$

$$2u_{0,0}-\left(\frac{1}{2}{u}_{0,1}+u_{1,0}+\frac{1}{2}{u}_{0,-1}\right)=\frac{1}{2},$$
以及内点差分方程$(1)$式
$$4u_{0,0}-\left(u_{1,0}+u_{0,1}+u_{0,-1}+u_{-1,0}\right)=1$$
可得：
$$u_{0,0}-\frac{1}{2}{u}_{1,0}-\frac{1}{2}{u}_{0,1}=\frac{1}{4}.$$

⑤对点$u_{0,4}$，同理可由左边界和上边界的式子分别得到：
$$2u_{0,4}-\left(\frac{1}{2}{u}_{0,5}+u_{1,4}+\frac{1}{2}{u}_{0,3}\right)=\frac{1}{2},$$

$$\frac{9}{4}{u}_{0,4}-\left(\frac{1}{2}{u}_{1,4}+u_{0,3}+\frac{1}{2}{u}_{-1,4}\right)=\frac{1}{2},$$
以及内点差分方程$(1)$式
$$4u_{0,4}-\left(u_{1,4}+u_{0,5}+u_{0,3}+u_{-1,4}\right)=1.$$
可得：
$$\frac{9}{8}{u}_{0,4}-\frac{1}{2}{u}_{1,4}-\frac{1}{2}{u}_{0,3}=\frac{1}{4}.$$

总结

在运用五点差分格式时，要把握好对初值及边值条件的求解，其步骤一般为：

1. 作草图，找正则内点和边界点；
2. 代入五点差分格式公式，得到正则内点表达式；
3. 分别根据边值及初值条件计算边界上的点和角位置的点；
4. 如果需要，可将方程组转化为矩阵形式。

在进行计算时，还可以整理成矩阵的形式，便于计算机求解。

参考

[1] 常微分数值方法部分课后习题答案