拉普拉斯方程的数值解.pdfVIP

金文璨 2007201230 作业6.1 一． 题目 有一个中空的无限长棱柱，其横截面如图所示，棱柱壁是金属的，棱柱中是金属 棒，在金属棒和外壁上分别加电压，计算金属棒和棱柱壁之间的电场，画出等势 线，电势和电场的分布。 二． 算法 金属棒为等势体，金属棒和外壁之间的区域内，电势满足Laplace 方程。解电势 φ(x,y)的Laplace 边值问题： 2 2 ∂ φ ∂ φ ∆φ= 2 + 2 = 0 ∂x ∂y φ(x=±1,y)=0；φ(x,y= ±1) = 0 φ(x, y) = 1;x∈ [−0.3.0.3],y∈ [−0.3.0.3] 用差分法表示： u + u + u + u i+1,j i−1, j i, j−1 i, j+1 u = i,j 4 （*） u = 0,u = 0,u