计算流体力学作业二

《计算流体力学作业二》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、计算流体力学椭圆型方程作业学院：海洋科学学院专业及班级：12级1班姓名:何敏学号：20121382014一、解析解程序programmain!主程序开始implicitnoneinteger,parameter::ii=31,jj=31!声明变量及赋值real,parameter::pi=3.1415926535898real,parameter::dx=0.1,dy=0.1real::u(0:ii,0:jj),x(0:ii,0:jj),y(0:ii,0:jj),f(0:ii,0:jj)integer::i,j,count,on

2、off,kopen(10,file='E:jsltlxzy2outputcs.txt')f=1.0forall(i=0:ii,j=0:jj)x(i,j)=dx*real(i)y(i,j)=dy*real(j)endforallf(0,0)=1.0doi=0,iif(i,0)=exp(x(i,0))enddodoj=0,jjf(0,j)=exp(y(0,j))enddodoi=1,iidoj=1,jjf(i,j)=exp(x(i,j)+y(i,j))enddoenddowrite(10,*)xwrite(10,*)ywri

3、te(10,*)fclose(10)end一、高斯赛德尔-点迭代程序programmain!主程序开始implicitnoneinteger,parameter::ii=31,jj=31!声明变量及赋值real,parameter::pi=3.1415926535898real,parameter::dx=0.1,dy=0.1real,parameter::TV_limit=0.001real::time,beta_x,beta_y,omegareal::u(0:ii,0:jj),x(0:ii,0:jj),y(0:ii,0:jj

4、),analyze(0:ii,0:jj),temp(0:ii,0:jj)integer::i,j,count,onoff,kcallinitial()!调用子程序，赋值计算初始值(以及边界条件)onoff=0count=0dowhile(count<50000)temp=ucallgauss_seidel(ii,jj,beta_x,beta_y,beta_xy,omega,u,dx,dy)!调用高斯-赛德尔迭代子程序count=count+1callTV_cal(onoff)!总变差计算,如果满足条件，输出计算结果，并停止程序e

5、nddocontainssubroutineinitial()!初值赋值子程序implicitnonebeta_x=dx**2/(dx**2+dy**2)/2.0beta_y=dy**2/(dx**2+dy**2)/2.0forall(i=0:ii,j=0:jj)x(i,j)=real(i)*dxy(i,j)=real(j)*dyendforalldoi=0,iidoj=0,jju(i,j)=exp(x(i,j)+y(i,j))enddoenddoendsubroutineinitialsubroutineTV_cal(onof

6、f)!总变差计算子程序implicitnonereal::TVinteger,intent(inout)::onoffonoff=0TV=0.0doj=1,jj-1doi=1,ii-1TV=TV+abs(u(i,j)-temp(i,j))!全场变差累加enddoenddoif(TV

7、racter*6::char_countcharacter*4::char_iicharacter*4::char_jjcallint_to_char(6,count,char_count)!调用整形数字符串转换子程序callint_to_char(4,ii+1,char_ii)callint_to_char(4,jj+1,char_jj)open(10,file="output/ADI3.plt")!数据文件名write(10,*)'title=T'!titlewrite(10,*)'variables="x","y","T"

8、'!variableswrite(10,*)'zonei='//char_ii//',j='//char_jj//',k=1,f=block'write(10,*)xwrite(10,*)ywrite(10,*)uwrite(10,*)countclose(10