流过平板的超声速流动.pdfVIP

D:\A\CFD\Flat_Plate\Flat_Plate\MacCormack.F901

!平板上的有粘超声速流动：完整的N-S方程的MacCormack数值解法

!变量说明：

!IMAX:网格的X方向个数

!JMAX:网格的Y方向个数

!MAXIT:循环推进限制的最大次数

!ITER：循环推进控制变量

!MACH_FAR:远场的马赫数

!V_FAR:远场的来流速度

!TEMP_FAR:远场的温度

!A_FAR:远场的声

!P_FAR:远场的压力

!ROU_FAR:远场的密度

!MIU_FAR:远场的动力粘性系数

!K_FAR:远场的热传导系数

!E_FAR:远场的单位微团内能

!LHORI:平板长度

!TW:物面温度（假设物面为常壁面温度条件）

!GAMMA:比热比

!PR:普朗特数

!MIU_SEA:标准海平面的动力粘性系数

!TEMP_SEA:标准海平面的温度

!R:气体常数

!CV：定容比热容

!CP:定压比热容

!REL:初始雷诺数

!DELTA:层流边界层厚度

!LVERT：计算域Y方向的高度

!DX,DY:X,Y方向步长

!IE:标示调用E/F中粘性应力项及热传导项的时刻

!U、E、F：N-S程原始的通量项的列向量形式

!U1:U1项初始时刻值及DT时间后校正值(U/E/F中的其他各分量同)，作为循环推进不断更新

!U1_ES:U1项T+DT时刻的预测值(U/E/F中的其他各分量同)

!VX,VY,ROU,P,E,TEMP,MIU,K:原始变量，作为循环不断推进更新变化

!ROU_PRE:用来存储前一时刻的密度值，与当前时刻的密度值比较后来作为判断是否达到收敛的条件

!定义常量模块

MODULEGLOBAL

IMPLICITNONE

REAL::TW0,GAMMA,PR,MIU_SEA,R,TEMP_SEA,CV,CP

PARAMETER(TW0=1.0,GAMMA=1.4,PR=0.71,MIU_SEA=1.7894E-5,R=287.,TEMP_SEA=288.16,CV=717.5,CP=1004.5)

ENDMODULE

!主程序

PROGRAMMAIN

USEGLOBAL

IMPLICITNONE

INTEGER::IMAX,JMAX,MAXIT,ITER,I,J

REAL::MACH_FAR,TEMP_FAR,A_FAR,P_FAR,ROU_FAR,V_FAR,LHORI,TW,REL,E_FAR,MIU_FAR,K_FAR,DELTA,LVERT,DX,DY,

DT

REAL,ALLOCATABLE::VX(:,:),VY(:,:),A(:,:),ROU(:,:),P(:,:),E(:,:),TEMP(:,:),MIU(:,:),K(:,:),U1(:,:),U2(

:,:),U3(:,:),U5(:,:),

E1(:,:),E2(:,:),E3(:,:),E5(:,:),F1(:,:),F2(:,:),F3(:,:),F5(:,:),ROU_PRE(:,:)

REAL,EXTERNAL::DYNVIS,THERMC

INTEGER,EXTERNAL::CONVER

OPEN(6,FILE=OUTPUT.DAT)

OPEN(5,FILE=INPUT.DAT)

!WRITE(*,(请输入划分网格X方向个数：))

READ(5,*)IMAX

!WRITE(*,(请输入划分网格Y方向个数：))

READ(5,*)JMAX

WRITE(6,101)IMAX,JMAX

101FORMAT(5X,*****输入的网格信息*****/,8x,IMAX=,I4/,8X,JMAX=,I4/)

!WRITE(*,(请输入循环推进最大次数限制：))

READ(5,*)MAXIT

WRITE(6,102)MAXIT

102FORMAT(5X,*****循环推进最大次数限制*****/,8x,MAXIT=,I4/)

D:\A\CFD\Flat_Plate\Flat_Plate\MacCormack.F902

!WRITE(*,(请输入来流马赫数：))

READ(5,*)MACH_FAR

!WRITE(*,(请输入