求助双马赫数反射van leer方法出现问题

一壶酒

1.出现了可以迭代求解ρ，ρu,ρv,ρE，再利用E=e+0.5*(u^2+v^2)和 p=（γ-1）ρe,求出p,但是会出现负压，即无法求出声速c(c=sqrt(γp/ρ))，无法判断Ma的大小。现在我用p<1则p=1来算 不知道合不合理（我这里在想我一直用当地的ρ和p来求c会不会其实用无穷远的参数来算或者声速c不是这么求得）
2.算出来的结果密度很大且最大点不是在x=1/6,y=0的地方 感觉出错了

以下是问题和我算的结果
希望有人指导下我

一壶酒

具体问题如下

LDK8

计算的初值是不是没有设置好啊，看你的视频里面初始时刻正激波前后的区域密度似乎是一样的，不过我也是初学