能指点一下我的二维burgers方程weno程序吗?

morxio

[这个贴子最后由morxio在 2007/04/24 07:33am 第 1 次编辑]

大家好，我用了两个月的时间学习二维weno格式的程序，
方程是 u_t + (u*u/2.0)_x + (u*u/2.0)_y =0
初值为：在四个象限里，从一象限到四象限分别为 -1.0, 0.5, -0.2, -1.0.
边值为：周期性边值条件。
用下面程序算出来，结果在 d:\x.txt , d:\y.txt , d:\pu.txt .
用matlab画图：plot3(x,y,pu)得到的图形没有之间的相互作用，只是沿着x方向描线,我改了很多遍都是这样。
前辈能不能给点建议或指点啊？乞求帮助！我算了两个月了，都快崩溃了！
谢谢了！！！

morxio

[这个贴子最后由morxio在 2007/04/24 09:15am 第 1 次编辑]

我用 tecplot10 画了一下图，出现了网格图，
但是问题是各个象限之间没有明显的相互作用,没有像
Efficient Implementation of Essentially Non-oscillatory Shock-Capturing Schemes II 里面Fig6-Fig7那样的图。

sphinx01

学习一下啊！！呵呵…………