采用lax格式用c语言模拟一维激波管问题

《采用lax格式用c语言模拟一维激波管问题》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、采用lax格式用C语言模拟一维激波管问题#include#include#include#defineN101#definej0intmain(){FILE*fp;doubletimeOut=0.006,lendom=10.0,dx=lendom/(N-1),DiniL,DiniR,Tini=300.0,PiniL=10.0,PiniR=1.0,dt0,Eini,aini,Hini,uini=0.0,gamm=1.4,R=287.0,CFL=0.98;//初

2、始条件DiniL=100000.0*PiniL/(R*Tini);DiniR=100000.0*PiniR/(R*Tini);//求初始密度Eini=R*Tini/(gamm-1.0);//求初始eaini=sqrt(gamm*R*Tini);//求初始声速Hini=gamm*R*Tini/(gamm-1.0);//求初始hdt0=CFL*dx/aini;inti,k;doubleD[2][N],E[2][N],a[2][N],H[2][N],U1[2][N],U2[2][N],dt[2][N],U3[2][N

3、],F1[2][N],F2[2][N],F3[2][N],u[2][N],C[2][N],P[2][N];doubletime,dt1;time=0;for(i=0;i

4、][i]*u[j][i];U3[j][i]=D[j][i]*E[j][i];F1[j][i]=D[j][i]*u[j][i];F2[j][i]=D[j][i]*u[j][i]*u[j][i]+P[j][i];F3[j][i]=D[j][i]*u[j][i]*H[j][i];}do{for(i=1;i

5、[i-1])/2-dt1/dx*(F2[j][i+1]-F2[j][i-1])/2;U3[j+1][i]=(U3[j][i+1]+U3[j][i-1])/2-dt1/dx*(F3[j][i+1]-F3[j][i-1])/2;u[j+1][i]=U2[j+1][i]/U1[j+1][i];E[j+1][i]=U3[j+1][i]/U1[j+1][i];D[j+1][i]=U1[j+1][i];P[j+1][i]=R*D[j+1][i]*(E[j+1][i]-u[j+1][i]*u[j+1][i]/2)/717.5

6、;H[j+1][i]=(E[j+1][i]-u[j+1][i]*u[j+1][i]/2)*1.4+u[j+1][i]*u[j+1][i]/2;//h=e/cv*cp,t=e/cvC[j+1][i]=sqrt(gamm*R*(E[j+1][i]-u[j+1][i]*u[j+1][i]/2)/717.5);F1[j+1][i]=D[j+1][i]*u[j+1][i];F2[j+1][i]=D[j+1][i]*u[j+1][i]*u[j+1][i]+P[j+1][i];F3[j+1][i]=D[j+1][i]*u[j

7、+1][i]*H[j+1][i];a[j+1][i]=(u[j+1][i]>0?C[j+1][i]+u[j+1][i]:C[j+1][i]-u[j+1][i]);//a=/c/+/u/dt[j+1][i]=CFL*dx/a[j+1][i];dt1=(dt[j+1][i-1]>dt[j+1][i])?dt[j+1][i]:dt[j+1][i-1];}time=time+dt1;U1[j+1][0]=U1[j+1][1];U1[j+1][N-1]=U1[j+1][N-2];U2[j+1][N-1]=U2[j+1][

8、N-2];U2[j+1][0]=U2[j+1][1];U3[j+1][N-1]=U3[j+1][N-2];U3[j+1][0]=U3[j+1][1];F1[j+1][0]=F1[j+1][1];F1[j+1][N-1]=F1[j+1][N-2];F2[j+1][N-1]=F2[j+1][N-2];F2[j+1][0]=F2[j+1][1];F3[j+1][N-1]=F3[j+1][N-