backward和C-N没有二阶精度？

[yhdthu]

偶然看到一个韩国人写的，说明了通量一阶处理无法使高阶时间格式达到应有的效果，大家怎么看。文章在这

[程迪]

@yhdthu
你管他呢，就说时间格式有二阶精度完了呗。解有没有二阶精度谁知道。而且时间精度还和时间步长有关呢。

[程迪]

@yhdthu 看了文章，似乎是对的。本来CN格式是2阶，但是解压力方程的时候阶数下降了。

感觉这个和压力方程的本质有关系，压力本来没有非定常项，都是通过速度算出来的。

所以里面的没修正的结果都是时间：速度2阶精度，压力1阶精度，最后到达一个平台。

我的理解是他是玩了个插值修正。把压力方程源项弄得更准了。

[yhdthu]

@程迪 是的，改起来也不是很麻烦，挺好

[程迪]

@yhdthu 呵呵，OF的精度问题不止如此，湍流模型如果是后加的，算子分裂可能会导致精度下降到1阶和稳态无法达到的问题。所以Jasak在FE中玩耦合求解是另外的套路。