cfd理論過渡到編程的傻瓜入門教程-資料下載頁
2025-06-28 08:13本頁面
【正文】 。j3。j++) cell[i].Wp[j]=cell[i].W[j]。 } for(j=0。j3。j++) cell[i].W[j]=cell[i].Wp[j]+ar[ir]*dt*cell[i].R[j]。 SolveWtoU(cell[i].W,cell[i].U)。 } }}//RungeKutta方法void SolveRungeKutta(){ double ar[4]={1/,1/,}。 int it,ir。 timesum=。 for(it=0。it++)//迭代步數(shù) { for(ir=0。ir4。ir++)//四步RougerKutta法 { boundarycondition()。 SolveResidual()。 SolveNextstep(ar,ir)。 } timesum+=dt。 printf(it=%d,timesum=%f\n,it,timesum)。 if(timesum=) { output()。 getchar()。 printf(Please enter any key to continue...)。 break。 } }}//結果輸出void output(){ int i。 FILE *fpd,*fpu,*fpp。 if((fpd=fopen(,w))==NULL) { printf(connot open infile\n)。 return。 } fprintf(fpd,TITLE = \TestCase\\n)。 fprintf(fpd,VARIABLES = \x\, \density\\n)。 fprintf(fpd,ZONE T=\Only Zone\, I=%d, F=POINT\n,Nc)。 if((fpp=fopen(,w))==NULL) { printf(connot open infile\n)。 return。 } fprintf(fpp,TITLE = \TestCase\\n)。 fprintf(fpp,VARIABLES = \x\, \pressure\\n)。 fprintf(fpp,ZONE T=\Only Zone\, I=%d, F=POINT\n,Nc)。 if((fpu=fopen(,w))==NULL) { printf(connot open infile\n)。 return。 } fprintf(fpu,TITLE = \TestCase\\n)。 fprintf(fpu,VARIABLES = \x\, \velocity\\n)。 fprintf(fpu,ZONE T=\Only Zone\, I=%d, F=POINT\n,Nc)。 for(i=0。iNc。i++) { fprintf(fpd,%f\t%f\n,cell[i].xc,cell[i].U[0])。 fprintf(fpu,%f\t%f\n,cell[i].xc,cell[i].U[1])。 fprintf(fpp,%f\t%f\n,cell[i].xc,cell[i].U[2])。// printf(x=%f,d=%f,u=%f,p=%f\n,cell[i].xc,cell[i].U[0],cell[i].U[1],cell[i].U[2])。getchar()。 } fclose(fpd)。 fclose(fpp)。 fclose(fpu)。}//網格類型判斷void SetMeshFlag(){ int i。 for(i=0。iNc。i++) { if(cell[i].xc=XS1amp。amp。cell[i].xc=XS2) cell[i].flag=0。 else if(cell[i].xcXS2) cell[i].flag=2。 else cell[i].flag=1。// printf(cell[%d].flag=%d\n,i,cell[i].flag)。getchar()。 }}//熵的計算void SovleEntropy(){ int i。 for(i=0。iNc。i++) { cell[i].S=cell[i].U[2]/pow(cell[i].U[0],gama)。 }}//主函數(shù)void main(void){ initialsolve()。 SetMeshFlag()。 SolveRungeKutta()。 output()。
試題試卷相關推薦
鄂ICP備17016276號-1