用5阶weno格式处理前台阶问题

binghen

这几日lz本科毕设，用5阶weno格式处理前台阶问题，作为一个cfd菜鸟，真是大问题小问题层出不穷。
我做出的结果如图1：
而标准图如图3：
我的初始变量以及边界条件的设置和论文中一模一样，另外为了验证我的边界条件，我把自己的边界条件移到2阶MUSCL算法上，可以取得非常好的效果，但是我自己写出的WENO-JS（5阶）的结果却如图一那样，和标准论文里的结果相去甚远。
我的左边界为入口，右边界为出口条件，上下边界为反射条件。
算法首先是用steger-warming矢通量分裂进行通量分解，然后用5阶WENO算法进行计算，对于带有虚拟点的模板我没有另外处理，我见有一份公共课件上说，对于带有虚拟店的模板，把其权重置为0，我试过把带有虚拟点的模板权重置为0，出现了更加匪夷所思的结果如图2：
看图片的感觉就好象该反射的地方没有反射，该出流的地方没有出流，好像边界特性都没有纳入考虑，所以我就删除了这一程序段，保留了原有的带有虚拟点的模板权重，并且没有做任何处理。就得到了图1，貌似还不错的图，但是对比标准图还是有些差距。
我试着重新调试程序，但是着实找不出原因，但是隐约感觉是反射边界处的通量处理的问题。
期待各位大神们的回复，帮帮菜鸟～

[ 本帖最后由 binghen 于 2014-5-10 15:47 编辑 ]