电势电场的分布.docVIP

数值模拟作业 姓名：杜海龙 学号 专业：等离子体 (一)求解在轴对称圆柱体中的电势跟电场的分布 一：数值方法 所研究的问题是在柱坐标系下，根据求解问题的轴对称性，由原来的三维问题转化为一个一维问题。 (一)数学模型： 求解在柱坐标系下的泊松方程： (1) 根据电场跟电势之间的关系： (2)边界条件是： (3) （二）模型离散化 ：离散的方法是有限差分法 （均匀网格差分） (4) (5) (6) 其中离子的密度假设是： （7） 将（5），(6),(7)以及边界条件带入（4）得三对角线性方程组，因此用追赶法（线性方程组的一种特殊的行形式）来求解； 求解电场的分布：（一阶三点法计算） () （非均匀网格差分）基本过程如下： （8） （9） 其中h为公比 将(8),(9),带入(1)，同样是得到对三角线性方程组，用追赶法求解； 二 具体程序如下：（参考计算物理中求解二阶线性微分方程程序） subroutine cw(N,XA,XB,A0,A1,A2,B0,B1,B2,A,B,EPS) DIMENSION A(N,N),B(N),x(n),e(n) H=(XB-XA)/(N-1) H2=H*H DO 1 I=1,N IF(I.EQ.1)B(I)=2.0*H*A2 IF(I.EQ.N)B(I)=2.0*H*B2 IF(I.NE.1.AND.I.NE.N)B(I)=2.0*H2*F(XA+(I-1)*H) 1 CONTINUE DO 3 I=1,N DO 3 J=1,N 3 A(I,J)=0.0 A(1,1)=2.0*H*A1-3.0*A0 A(1,2)=4.0*A0 A(1,3)=-A0 A(N,N-2)=B0 A(N,N-1)=-4*B0 A(N,N)=3*B0+2*H*B1 DO 2 I=2,N-1 i1=i-1 do 2 j=i-1,i+1 IF(I-1.EQ.J)A(I,J)=2.0-P(XA+(I-1)*H)*H IF(I.EQ.J)A(I,J)=2*H*H*Q(XA+(I-1)*H)-4 IF(I+1.EQ.J)A(I,J)=2+P(XA+(I-1)*H)*H 2 CONTINUE CALL GSL(N,A,B,EPS) do 10 i=1,n x(i)=XA+(I-1)*H WRITE(10,*) x(i),b(i) open(10,file=fu.dat) write(*,*)x(i),b(i) 10 continue e(1)=(-3*b(1)+4*b(2)-b(3))/2*h e(100)=(b(98)-4*b(99)+3*b(100))/2*h do i=2,99 e(i)=(b(i-1)-b(i+1))/2*h end do open(11,file=ni.dat) do i=1,100 write(*,*)x(i),e(i),b(i) write(11,*)x(i),e(i) enddo RETURN END PROGRAM TCWB DIMENSION A(100,100),B(100) open(10,file=fu.dat) A0=0.0 A1=1.0 A2=-1 B0=0.0 B1=1.0 B2=100 XA=3.0 XB=5.0 N=100 EPS=1.0e-10 CALL cwb(N,XA,XB,A0,A1,A2,B0,B1,B2,A,B,EPS) STOP END function q(x) q=0.0 return end FUNCTION P(X)P=1.0/x RETURN END FUNCTION F(X) pr=8.8