请教一个二阶迎风格式的问题

runheart

我自己做了个二维同位网格simpler算法的程序，包含了，一阶迎风，指数，quick，和二阶迎格式。但计算发先，前三种格式都能算好层流blasius解，只有二阶迎有问题：
   二阶迎风是这样给出的：
   Qi+1/2=Qi＋DQi.DRi         Ui+1/2 > 0
   Qi+1/2=Qi+1＋DQi+1.DRi+1   Ui+1/2 < 0
  Q 为（u ,v) ；i为节点 ；DQ为梯度； DR为节点到界面的距离矢量。
附图：二阶迎风的结果：
  请教各位了，谢谢！！