求助关于二维前台阶超声速无粘流编程

zhanglonglin

最近在编个二维前台阶超声速无粘流的问题，试过ＮＮＤ格式和MUSCL格式，但都出现了第一道激波迁移，反射激波逐渐消失的现象，好像后面通道堵住了。估计原因应该是边界处理的问题。入口和出口都是外推边界，问题可能在台阶拐角处几个点的处理（我用的是张德良的那本书上的处理方法）。不知哪位大侠能够指出可能的问题和解决办法。先谢过了