高等计算流体力学-08.pptxVIP

第八讲求解二维Euler方程的计算程序基本特征结构网格二阶有限体积方法重构:准一维重构+TVD-MUSCL插值Riemann Solver: Lax, Roe时间推进: 二阶TVD Runge-Kutta边界处理:基于Roe格式的简化特征边界条件 数据结构定义数组维数时imax=im+s,jmax=jm+s, s为虚拟网格层数im,jmib=3,jb=3i+1/2,j+1/2 i+1,j+1i+1/2,j i+1,ji-1/2,j+1/2 i,j+1i,ji-/2,ji,ji,j-1/2i,ji-1/2,j-1/2i,j变量定义module main parameter(iq=405,jq=105,ngmax=max(iq,jq)+5) !数组维数 parameter(small=1.0d-20,small1=1.0d-10,alarge=1.0d20) ! 小量和大量 parameter(one=1.0,zero=0.0) ! 1和0 logical :: steady1 ！steady1: logical variable, ==.true. compute steady flow; ==.false. the flow is unsteady integer :: ib,jb,im,jm,irsolver ！ib,jb,im,jm网格指标；irsolver， irsolver=0 Roe; =1 Lax -riemann solver real :: basel,basev,baser,baset,basem ！无量刚基准，长度-速度-密度-温度-马赫数Real:: ro(iq,jq),p(iq,jq),roe(iq,jq),tp(iq,jq),rovx(iq,jq),rovy(iq,jq),vx(iq,jq),vy(iq,jq),ho(iq,jq), ！ρ，p， ρE（单位体积总能），T（温度），ρu， ρv，u， v，总焓 , hro(iq,jq),hre(iq,jq),hrx(iq,jq),hry(iq,jq) ！ 通量代数和（方程右端项）连续-能量-动量,rom1(iq,jq),roem1(iq,jq),rovxm1(iq,jq),rovym1(iq,jq) ！N时刻守恒变量real :: x(iq,jq),y(iq,jq),xc(iq,jq),yc(iq,jq),step(iq,jq),vol(iq,jq)！控制体角点和中心点坐标，时间步长，控制体面积 real :: aix(iq,jq),aiy(iq,jq),ajx(iq,jq),ajy(iq,jq),！控制体面积矢量 Snx，Snyaim(iq,jq),ajm(iq,jq) ! 控制体面积大小 ,alagmx(iq,jq),alagmy(iq,jq)控制体长度和宽度控制体i-界面左右两侧的速度、压力和密度控制体j-界面左右两侧的速度、压力和密度real :: vib(iq,jq),vjb(iq,jq) ！控制体界面中点运动的法向速度（动网格计算中有用） real :: uil (iq,jq),uir(iq,jq), vil (iq,jq),vir(iq,jq), pil (iq,jq),pir(iq,jq), ril (iq,jq),rir(iq,jq), ujl (iq,jq),ujr(iq,jq), vjl (iq,jq),vjr(iq,jq), pjl (iq,jq),pjr(iq,jq), rjl (iq,jq),rjr(iq,jq) integer ::nmax,n,irest,ioup,iinp, icgl, iimplicit,idgstart,nts1p, nunsloop,ninnerloop,ip,iprint,istart,nb！ nmax: maximum advancing steps；!irest: start from very beginning(0) or start from a previous computation(1);!iprint: output every iprint steps！ icgl=1: gloal time step; =0 local time step！ip=1输出计算结果； n：current time-step（累积）；istart： current time-step（当前计算） real :: cfl,ttime,tstep,epsall ! Cfl数，计算时间，时间步长，收敛标准 real ::cp,ga,pi,fga,ga1,rcp,rcpcv,rfga,cv, prl,prt,prte！Cp，γ，π