高等计算流体力学-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,i