前台阶数值计算结果，求分析问题。

逐冰格

我格式精度不高，计算得到前台阶如图。我想问，只是根据这幅图可以知道我计算过程有哪些问题吗？并不是有一隐瞒格式，只是初学，想听听大家的分析。学习一下。用的是laxwendroff格式。谢谢大家了。

逐冰格

再另发一幅计算结果的图。不同算法啊。结果又有什么问题呢？

zhujiejie

两幅图所示结果计算精度不同，跟计算网格已经空间、时间离散格式精度都有关系，没有明显错误

glandetian

前台阶问题的计算结果是随时间变化的。
第一幅图的结果是时间为t=2.5左右的结果吗？
如果是，我的结果跟你稍微有点不同，我的结果是已经在上壁面出现了马赫杆。

[ 本帖最后由 glandetian 于 2014-5-7 16:20 编辑 ]

逐冰格

我的结果是时间t=3的结果了。但和2.5几乎没有差别。应该是收敛了吧。

逐冰格

你说我继续算下去也没有马赫杆出现，这是为什么？或者出现了都是很短的。没有继续发展变长。这是怎么回事呀？

逐冰格

你也是用的laxwendroff格式吗？

glandetian

请参看文献《高精度算法与小波多分辨分析》，第131页。
结果是随时间变化的，如下图所示，这是我计算的结果，是密度等值线。
当然，最后会得到定常的结果，上述文献中计算到了t=10的结果，大概就是收敛的结果了。

t=2.5和t=3.0时，结果是不同的。
我用的龙格库塔-间断Galerkin方法，数值通量是Lax - Friedrichs型的。

[ 本帖最后由 glandetian 于 2014-5-9 10:45 编辑 ]

逐冰格

谢谢了。辛苦了。十分感动啊。