二维差值实验报告.docVIP

中山大学 本科生实验报告书 二维差值的fortran编程 院 系 工学院应用力学与工程系 专业班级 理论与应用力学10级 实验课程 计算流体力学 姓 名 程 彬 学 号 实验地点 中山大学东校区实验中心 实验时间 2012年 4 月 5日 指导教师 詹杰民 一、实验问题： 1、编程用迎风C格式、A格式在=0.5, 1, 2三种情况下编程计算。 2、计算范围为在，取计算至90个的时间步长。 3、讨论两种格式，在=0.5, 1, 2下的差分格式的稳定性。选择另外一种稳定地二阶精度格式，进行编程计算和对比分析。 4、分别绘制出三种格式在不同时间t下的图，并选择若干将其制作成为动画。 二、实验程序 1.主程序 !------------------主程序（变量声明赋值）------------------- !A=u(x,y,t)定义一个三维变量 !i=x坐标变量，j=y坐标变量，k=时间变量，m=x,y初值边界 !dt/dx=P,dt/dy=Q program main implicit none integer i,j,k,m character(len=2)::time real A(-106:106,-106:106,0:90) real P,Q P=.5 Q=.5 m=16 !-------------------------------------调用子程序----------------------------------- call subl(a,p,q,m) !-------------------------------------输出txt格式文件------------------------------ open(10,file=test0425L05.txt) do k=0,90 write(10,100)((A(i,j,k),j=-16,16),i=-16,16), !按照33行，33列输出每个时层的u(x,y,t),并空出一行 enddo 100 format(1x,33G10.2) !将结果存储到test0425A05.txt !----------------------输出dat格式文件----------------------------- do k=0,90 write(time,(i2))k open(unit=20,file=testL05//trim(time)//.dat) write(20,*) variables=,x,,,y,,,z do i=-16,16 do j=-16,16 write(20,*) i,j,a(i,j,k) enddo enddo enddo end program main 2.A格式子程序 subroutine suba(A,P,Q,M) real A(-106:106,-106:106,0:90) A(-106:106,-106:106,0)=0 A(-4:4,-4:4,0)=1 !给三维变量A(x,y,t)赋初值:先全部赋初值为0，再对-4=x=4,-4=y=4赋初值1 !按照A格式计算1到90层的值 do k=1,90 !按时间层循环 do i=-106+k,106-k !循环x变量，从-105到105 do j=-106+k,106-k !循环y变量，从-105到105 A(i,j,k)=A(i,j,k-1)-P/2.0*(A(i+1,j,k-1)-A(i-1,j,k-1))-Q/2.0*(A(i,j+1,k-1)-A(i,j-1,k-1)) enddo enddo enddo end subroutine 3.拉克斯格式子程序subroutine subl(A,P,Q,M) real A(-106:106,-106:106,0:90) A(-106:106,-106:106,0)=0 A(-4:4,-4:4,0)=1 !给三维变量A(x,y,t)赋初值:先全部赋初值为0，再对-4=x=4,-4=y=4赋初值1 !按照A格式计算1到90层的值 do k=1,90 !按时间层循环 do i=-106+k,106-k !循环x变量，从-105到105 do j=-106+k,106-k !循环y变量，从-105到105 A(i,j,k)=A(i,j,k-1)-0.5*P*(A(i+1,j,k-1)-A(i-1,j,k-1))+0.5*P**2*(A(i