资源描述：
《二阶迎风型差分格式求解一维可压缩黏性流动问题fortran语言》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、Fortran77语言源程序UpwindTVD_lD.f二阶迎风型TVD差分格式求解一维可压缩黏性流动问题（Fbrtran77语言版本）programUPWIND_TVD_1Dimplicitreal氺8(a-h，o-z)parameter(mx二201，Tt=5.0)dimensionQ(3,mx),Qold(3,mx)!Q:[rou，rou*u，E]dimensionrou(mx)，u(mx),p(mx),T(mx),E(mx)real氺8Macommon/paradef/Sf,Re,Ma,pr,gama,cp,cv,
2、dx,dt，imim=mxcallInitialize⑼timc=0.0n=01n=n+ltime=time+dtdoi=l,imdo1=1,3Qold(l,i)=Q(l,i)enddoenddocallUpwindTVD_lD_Solver(Q)if(n/10000*10000.eq.n)thenwrite(*,20)n,time,error(Q,Qold)callOutPut(Q)endif20format(lx,i10,’time=’，el6.9,’error=’，el6.9,lx)if((time.It.Tt).a
3、nd.(error(Q,Qold).gt.1.0e~8))thengoto1endifcallOutPut(Q)write(氺，氺)’Programend!’endisubroutineTnitialize(Q)implicitreal氺8(a-h,o-z)real氺8Macommon/paradef/Sf,Re,Ma,pr,gama,cp,cv,dx,dt，imdimensionrou(im)，u(ini)，p(ini),T(im),E(im),Q(3,im)Re=50.0!雷诺数Ma=2.0!马赫数pr=0.72!普朗
4、特数gama=l-4»气体常数cv=l.0/(gama*(gama-1.0)*Ma**2)!比定容热容cp=gama氺cv!比定压热容dt=l.0e~6dx二1.0/(1.0*(im—1))!空间步长xl=O.0xr=l.0rouO=l.0uO=l.0T0=1.0temp=(gama+1.)/(gama~l.)ul=(2./(gama-1.)+Ma**2)/(temp*Ma**2)Tl=(2.*gama/(gama+1)/temp)*(l./temp+2./(gama+1.)/Ma**2)doi=l,imrou(i)=ro
5、uOxi=(i-l)本dxu(i)=(xi-xl)/(xr-xl)*(ul-uO)+u0T(i)=(xi-xl)/(xr-xl)*(Tl-TO)+TOP(i)=rou(i)*T(i)/(gama*Ma**2)E(i)=rou(i)*(cv*T(i)+0.5*u(i)**2)enddodoi=l,imQ(l,i):rou(i)Q(2,i)=rou(i)*u(i)Q(3,i)=E(i)enddoendsubroutineQ2U(Q,rou,u,p,T,E,a)implicitreal氺8(a~h,o_z)real*8Macom
6、mon/paradef/Sf,Re,Ma,pr,gama,cp,cv,dx,dt,imdimensionQ(3,im),rou(im),u(im)dimensionp(im),T(im),E(im),a(im)doi=l,imrou(i)=Q(l，i)u(i)=Q(2,i)/Q(l,i)E(i)=Q(3，i)T(i)=(E(i)/rou(i)-0.5*u(i)**2)/cvp(i)=rou(i)*T(i)/(gama*Ma**2)a(i)=dsqrt(T(i))/MaenddoendsubroutineUpwindTVDI
7、DSolver(Q)implicitreal氺8(a~h,o~z)real*8Ma,lamda,minmodcommon/para_def/Sf,Re,Ma,pr,gama,cp,cv,dx,dt,imdimensionQ(3,im),Qold(3,im),dq(3,im),rou(im),u(im),p(im),T(im)dimensionE(im),a(im),f(3,im),fwave(3,im),fr(3),fl(3)dimensionGv(3,im),g(3,im),lamda(3,im),sigma(3,iin
8、)dimensionalfa(3,im),Rn(3,3,im),R(3,3,im),beta(3,im)callQ2U(Q,rou,u,p,T,E,a)doi=l,imf(l,i)=rou(i)*u(i)f(2,i)=rou(i)*u(i)**2+p(i)f(3，i)=u(i)*(E⑴+p⑴)enddo
