python matplot包绘制有限体积解

python matplotlib包绘制非结构(三角)网格有限体积解

在偏微分方程数值解的有限体积方法中，方程的解是控制体上变量的均值。可以使用Python matplotlib中的函数完成有限体积解的绘制。

matplotlib.tripcolor

该函数是用于绘制三角网格上伪色彩图，可以用于表示有限体积解，官方说明文档点击这里.

plt.tripcolor(triangle,rho,edgecolors='none',cmap=plt.cm.rainbow,alpha=0.8)

输入参数如下：

1. triangle 三角网格信息。包含三角形顶点的坐标以及三角形和顶点的对应关系。详见matplotlib.tri模块。

也可以分开输入。注意三角形的顶点编号是从0开始的，而MATLAB函数中对应是从1开始的。

1. rho是要绘制的变量。或者写成facecolors=rho.
2. edgecolors 是三角网格边的颜色。
3. cmap 是matplotlib库中绘图colormap.
4. alpha是透明度。

下面的数值解结果是Euler方程的sod激波管问题。假定流体是理想气体，比热比$\gamma =1.4$,方程的初值如下，
$$({\rho }_L,u_L,v_L，p_L)=(1.0,0,0,1.0)$$

$$({\rho }_R,u_R,v_L,p_R)=(0.125,0,0,0.1)$$

数值结果如下：

在$y=0.25$处的密度分布如下：

在$x=0.7$处有接触间断。