非结构网格的粘性物面的处理方法

liujw512

我在编三维粘性非结构网格的程序时，遇到了一点问题，就是粘性物面的处理不对，压力等值线是锯齿性的，其他的参数也感觉不是很对，我怀疑是导数求的不对。我在处理边界时引入了虚拟网格，内网格的各物理量的导数可以由高斯定理得到，如何得到物面对应的虚拟网格的导数呢？不知道大家有什么看法？非常谢谢。

brittops

我做非结构时面的导数是由左右两侧单元的导数得到的，对边界面，没有使用虚网格，根据单侧单元的导数得到。