cfd理論過渡到編程的傻瓜入門教程(參考版)
2025-07-01 08:13本頁(yè)面
【正文】 }。 SolveRungeKutta()。 }}//主函數(shù)void main(void){ initialsolve()。iNc。 }}//熵的計(jì)算void SovleEntropy(){ int i。// printf(cell[%d].flag=%d\n,i,cell[i].flag)。 else if(cell[i].xcXS2) cell[i].flag=2。amp。iNc。}//網(wǎng)格類型判斷void SetMeshFlag(){ int i。 fclose(fpp)。getchar()。 fprintf(fpp,%f\t%f\n,cell[i].xc,cell[i].U[2])。i++) { fprintf(fpd,%f\t%f\n,cell[i].xc,cell[i].U[0])。 for(i=0。 fprintf(fpu,VARIABLES = \x\, \velocity\\n)。 return。 fprintf(fpp,ZONE T=\Only Zone\, I=%d, F=POINT\n,Nc)。 } fprintf(fpp,TITLE = \TestCase\\n)。 if((fpp=fopen(,w))==NULL) { printf(connot open infile\n)。 fprintf(fpd,VARIABLES = \x\, \density\\n)。 return。 FILE *fpd,*fpu,*fpp。 break。 getchar()。 printf(it=%d,timesum=%f\n,it,timesum)。 SolveNextstep(ar,ir)。ir++)//四步RougerKutta法 { boundarycondition()。it++)//迭代步數(shù) { for(ir=0。 timesum=。 } }}//RungeKutta方法void SolveRungeKutta(){ double ar[4]={1/,1/,}。j++) cell[i].W[j]=cell[i].Wp[j]+ar[ir]*dt*cell[i].R[j]。 } for(j=0。j3。i=Nc。 } } }}//流場(chǎng)值更新void SolveNextstep(double ar[],int ir){ int i,j。j3。 //二階格式 // SolveAUSMFlux(cell[i1].U,cell[i].U,Fm)。//一階格式 //i1/2 SolveReconstruction(cell[i2].U, cell[i1].U, cell[i].U, cell[i+1].U, UL, UR)。 //MSUCL重構(gòu) SolveAUSMFlux(UL,UR,Fp)。iNc。 double UL[3],UR[3],Fp[3],Fm[3]。 Fc[2]=*Mn*((WL[2]+PL)*cl+(WR[2]+PR)*cr)*fMn*((WR[2]+PR)*cr(WL[2]+PL)*cl)。 Fc[0]=*Mn*(WL[0]*cl+WR[0]*cr)*fMn*(WR[0]*crWL[0]*cl)。 if(fabs(Mn)deta) fMn=fabs(Mn)。} Mn=Ml+Mr。} else {Mr=MR。} else if(fabs(MR)) {Mr=*(MR1)*(MR1)。} if(MR=) {Mr=。} else {Ml=。} else if(fabs(ML)) {Ml=*(ML+1)*(ML+1)。 if(ML=) {Ml=ML。 ML=UL[1]/cl。 PR=UR[2]。 PL=UL[2]。 SolveUtoW(WL,UL)。 double pl,pr,pn,PL,PR。 double deta,fMn。//u W[2]=U[2]/()+*U[0]*U[1]*U[1]。}//基本變量轉(zhuǎn)化為守恒變量void SolveUtoW(double W[3],double U[3]){ W[0]=U[0]。//d U[1]=W[1]/W[0]。 UL[i]=UI[i]+*sL[i]*((1+KAKA*sL[i])*aL[i]+(1KAKA*sL[i])*bL[i])。i3。 sL[i]=(*aL[i]*bL[i]+)/(aL[i]*aL[i]+bL[i]*b
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1