【paraview进行MFiX后处理】如何计算某一气体组分的质量流量/气体产量

MFiX版本：20.1.2
paraview版本：5.8
目的：在一个反应容器中，发生反应并产生H2，我们需要算出H2的质量流量，并且测定出气体产量。

当前MFiX只能输出H2气体质量分数，但是没有办法输出总质量或者质量流量，所以问题归结于求质量流量，并且进行时间积分。

求总产量，相当于求对某一切面流过的质量流量积分。

1. 导入流场。

2. 做切面。

点Y normal可以直接把切面与Y轴垂直
如果没显示切面，勾选show plane
在 origin中输入切面所在位置

3. 用calculator 计算密流

2020年7月30日更正

$$\dot{m}=\rho VA$$
需要先计算密度rho，然后计算rho*V，最后在面上积分，就相当于乘以了A

点击左上角的计算器

先用如下公式计算出密度（如上图所示）

P_G/((8314/Gas_Mixture_MW)*Gas_temperature)

注意：只有高版本的MFiX才可以输出相应的气体常数，并且要在模拟前勾选上

然后再加一个calculator，计算出rho*V
如图所示

4. 对切面积分

在filter > data analysis中有一个名为integrate variables的filter
加上以后别忘了在跳出的表格视图中选择显示attribute > cell data

此时显示的一行数据就是切面积分后的数据

【2020-12-9补充】
intgrate variables这个过滤器其实是对面做积分之后再除以总面积，相当于是积分后再平均。
（见附注：验证integrate variables过滤器的作用[^1]）

5. 作流量随时间变化图

先选中表格中那行数据
然后使用plot variable over time
显示以后，取消其他无关的变量，只勾选Result
就会出现流量随时间变化图

补充：求4种气体的产量

在面积分和求rhoV之间添加四个计算器即可。
计算的值是rhoV*相应的气体的质量分数

如图（可以在过滤器上右键重命名）

四个计算器分别是

得到的intgrateVariable1显示的值为

最后再plot over time

得到结果

注意，得到结果后需要除以相应的分子量才得到摩尔分数

验证MFiX中的monitor与paraview

设定paraview中的monitor为

得到的结果为
paraview得到的结果为（摩尔分数）

两者在绝对值上是不同的，但是变化趋势是差不多的。
用origin积分后发现
这是monitor得到的

这是paraview得到的

两者差距还可以，但是还是有些误差的。

附注：验证integrate variables过滤器的作用[^1]

核算下面积分到底是不是积分再平均
上面integrate得到的rhoV_CH4是0.00075

CH4显示的值为
导出后发现rhoV*CH4平均值是0.0125
而切面的面积是0.06

0.0125*0.06恰好是0.0125.

故integrate variables这个过滤器其实是面积分然后再平均。