南京大学计算物理课程课件5.pptVIP

计算物理 李文飞 南京大学物理系 E-mail: wfli@nju.edu.cn 流程图 开始 输入N,M,V,H 初始边界 条件赋值 K=1,2,…,M I=1,2,…,N 重新赋边值 输出 * * 2011-3-29 E=mc2 热传导方程 ---偏微分方程的数值解法(抛物型) 对内部无热源物体： 热传导方程 Q u(x,y,z) ds dV 一维热传导方程 l 思路：用差分代替微分 计算步骤： x t k+1 k i-1 i i+1 稳定条件 边界条件差分格式： 第一类边界条件 第二类边界条件 第三类边界条件 举例： 计算过程： 1 u=0 u=0 x t k i-1 i i+1 网格法 DIMENSION U(0:200,0:200) OPEN(10,FILE=CONDUCT1.DAT) A=1/6. H=0.005 DO K=0,200 U(0,K)=0. U(200,K)=0. ENDDO DO I=1,199 U(I,0)=4*I*H*(1-I*H) ENDDO DO K=0,200 DO I=1,199 U(I,K+1)=A*U(I+1,K)+(1-2*A)*U(I,K)+A*U(I-1,K) ENDDO ENDDO DO I=0,200 WRITE(10,*)I,U(I,200) ENDDO CLOSE(10) END 计算结果： 二维扩散方程的差分解法 砷 J 线积分变面积分： 若粒子流仅仅由扩散导致，则 二维扩散方程 建立差分格式： M1 M2 (i, j) M N O 对点(i, j),在k时刻有： 代入扩散方程： 整理得递推公式： i-1 i+1 i k k+1 j-1 j j+1 t x 边界条件： 流程图 开始 输入N,M,M1,M2,D,H,MT 初始边界 条件赋值 S= TD / H2 K=1,2,…,MT 重新赋边值 输出 I=2,3,…,M-1 J=2,3,…,M-1 DIMENSION C1(51,71),C2(51,71) OPEN(10,FILE=DIFFUSION.DAT) READ *, N,M,M1,M2,D,H,T,MT N1=N+1 MX=M+1 S=T*D/(H*H) IF(S.GT.0.25)STOP IF(N1.GT.51.OR.MX.GT.71)STOP DO I=1,N1 DO J=1,MX C1(I,J)=0.0 ENDDO ENDDO DO I=M1,M2 C1(1,I)=1.0 ENDDO DO K=1,MT DO I=2,N DO J=2,M C2(I,J)=(1-4.0*S)*C1(I,J)+S*(C1(I+1),J)+C1(I-1,J)+C1(I,J+1)+C1(I,J-1)) ENDDO ENDDO DO J=2,M DO I=2,N C1(I,J)=C2(I,J) ENDDO C1(N1,J)=C2(N,J) ENDDO DO I=2,M1 C1(1,I)=C2(2,I) ENDDO DO I=M2,M C1(1,I)=C2(2,I) ENDDO WRITE(10,*)((I,J,C1(I,J),I=1,N1),J=1,MX) ENDDO END 计算结果： 波动方程 ---偏微分方程的数值解法(双曲型) 弦线的横振动方程 线密度 张力 外力 对均匀弦线，无外力的自