DG法后处理问题

zqb138

各位做DG的是怎么把间断解连续光滑化的，否则这个等值线没法画出来啊？还有一个画流线的问题，我们只有流场的流速数据，如何用matlab画流线图呢？

glandetian

用软件Tecplot， 按它的格式输出计算结果 不需要把间断解连续光滑，按单元依次输出结果就可以了。流线它自己会画出来。

zqb138

我是这样理解的：间断有限元的是间断解，同一个物理量在同一个顶点会有好几个值，画等值线那就有矛盾了，是否是tecplot已经将其连续化了然后才画出的图？我自己一直是用matlab画图的，对tecplot不是了解。看到一篇文献上是将分片间断有限元空间向连续有限元空间做L2投影得到连续化结果。

zqb138

tecplot中一个点的一个物理量只能给一个数值啊，怎么按单元输出数据啊，望能给个简单例子参考一下啊！

youy624

用matlab画流线我没有试过，但是画过等温线，希望通用。命令contour。当然前提是知道流函数这个数组吧。

glandetian

Title=  " Euler-DG "
VARIABLES= " X "," Y "," rho "," u "," v "," E "," p "," March "
Zone n=4,e=1,f=fepoint, et=quadrilateral
5.000000e+000  0.000000e+000  1.000000e+000  3.800000e-001  -8.820019e-015  1.857914e+000  7.142857e-001  1.000000e+000
4.962000e+000  4.341200e-001  1.000000e+000  3.800000e-001  -1.380393e-014  1.857914e+000  7.142857e-001  1.000000e+000
4.364975e+000  3.818850e-001  1.000000e+000  3.800000e-001  7.254945e-015  1.857914e+000  7.142857e-001  1.000000e+000
4.398400e+000  0.000000e+000  1.000000e+000  3.800000e-001  5.081015e-015  1.857914e+000  7.142857e-001  1.000000e+000
1 2 3 4
Title=  " Euler-DG "
VARIABLES= " X "," Y "," rho "," u "," v "," E "," p "," March "
Zone n=4,e=1,f=fepoint, et=quadrilateral
4.398400e+000  0.000000e+000  1.000000e+000  3.800000e-001  5.081015e-015  1.857914e+000  7.142857e-001  1.000000e+000
4.364975e+000  3.818850e-001  1.000000e+000  3.800000e-001  7.254945e-015  1.857914e+000  7.142857e-001  1.000000e+000
3.767950e+000  3.296500e-001  1.000000e+000  3.800000e-001  -5.726223e-014  1.857914e+000  7.142857e-001  1.000000e+000
3.796800e+000  0.000000e+000  1.000000e+000  3.800000e-001  -4.114518e-014  1.857914e+000  7.142857e-001  1.000000e+000
1 2 3 4
Title=  " Euler-DG "
VARIABLES= " X "," Y "," rho "," u "," v "," E "," p "," March "
Zone n=4,e=1,f=fepoint, et=quadrilateral
4.962000e+000  4.341200e-001  1.000000e+000  3.800000e-001  -1.380393e-014  1.857914e+000  7.142857e-001  1.000000e+000
4.924000e+000  8.682400e-001  1.000000e+000  3.800000e-001  -2.586421e-014  1.857914e+000  7.142857e-001  1.000000e+000
4.331550e+000  7.637700e-001  1.000000e+000  3.800000e-001  1.356063e-014  1.857914e+000  7.142857e-001  1.000000e+000
4.364975e+000  3.818850e-001  1.000000e+000  3.800000e-001  7.254945e-015  1.857914e+000  7.142857e-001  1.000000e+000
1 2 3 4
Title=  " Euler-DG "
VARIABLES= " X "," Y "," rho "," u "," v "," E "," p "," March "
Zone n=4,e=1,f=fepoint, et=quadrilateral
4.364975e+000  3.818850e-001  1.000000e+000  3.800000e-001  7.254945e-015  1.857914e+000  7.142857e-001  1.000000e+000
4.331550e+000  7.637700e-001  1.000000e+000  3.800000e-001  1.356063e-014  1.857914e+000  7.142857e-001  1.000000e+000
3.739100e+000  6.593000e-001  1.000000e+000  3.800000e-001  -1.026514e-013  1.857914e+000  7.142857e-001  1.000000e+000
3.767950e+000  3.296500e-001  1.000000e+000  3.800000e-001  -5.726223e-014  1.857914e+000  7.142857e-001  1.000000e+000
1 2 3 4




上面给出了一个例子，（为了提高显示的精度）把一个单元分成了4个四边形单元输出

zqb138

又有一个问题了，画这些数据的二维平面图倒是没问题，我用这些数据画三维立体图时总是只能画出rho的，后面的都出不来，也不知道在哪里设置改变物理量，麻烦了。

eryun05

用tecplot画图方便，建议使用它

zqb138

试了用tecplot，感觉没有matlab方便，fortran可以与matlab混合编程，在fortran中的计算结果图可以随算随显随保存，命令可以在fortran中编好，感觉更方便。不爽的就是matlab无法根据速度场画流线图，还有就是等值线无法自动标注数值。可能还是各有优缺点吧。

[ 本帖最后由 zqb138 于 2012-6-7 15:12 编辑 ]