DSMC时间平均

[kuangyonglin]

dsmc求解中，由于是使用统计学相关的算法进行求解，计算结果具有较大的统计噪音，使用时间平均法是降噪的一种有效方法.

作者在学习相关软件时，发现tecplot具有pytceplot拓展模块，能够实现此功能.

首先需要在系统中装python环境，并且安装pytecplot模块

直接使用pip install pytecplot命令即可安装

在官方github仓库中找到tputils.py, tpmath.py, TimeAverage.py3个文件，将其放到OpenFOAM的结果文件夹下

使用tecplot打开contralDict文件，并且打开软件的PyTecplot Connections功能

在终端使用python环境运行TimeAverage.py脚本，待程序运行完成，在zone style中将会出现时间平均后的计算结果

以某微纳通道为例，进行时间平均，并做对比：

时间平均前的结果：

时间平均后的结果：

显然，时间平均起到了降噪的作用

[Nonocore]

您好，请问这种方法是针对DSMC瞬态仿真结果处理的吗？
我在使用DSMC进行稳态算例仿真时，收敛后利用dsmcfoam自带的时间平均功能，采样平均了很多时间步后，最终流场图结果仍然有很多明显的噪声，请问这种情况的原因是什么呢？是因为我的DSMC粒子数太少了吗？我只保证了每个网格大概10个DSMC粒子。

[kuangyonglin]

@Nonocore 从我的经验来看，若是统计宏观量，10个DSMC粒子数是够用了，但是要观察流场细节，尤其是在Kn数较大的情况，10个DSMC粒子是不够的，需要更多的粒子 和 更多的时间平均才能较好的降噪.
有问题也可以随时邮件联系

[Nonocore]

@kuangyonglin 好的好的，谢谢您的建议。

[jeffx]

@kuangyonglin 您好，请问在使用dsmcfoam+进行计算时，初始运行一段时间步数，计算就会停止是什么原因（并没有到达设定的endTime）

[kuangyonglin]

@jeffx 在 DSMC时间平均 中说：

@kuangyonglin 您好，请问在使用dsmcfoam+进行计算时，初始运行一段时间步数，计算就会停止是什么原因（并没有到达设定的endTime）

看不出来啊，看着 controlDict 看着挺正常的

[jeffx]

@kuangyonglin 请问您是使用的hyStrath版本的dsmcfoam+求解器(https://hystrath.github.io/solvers/fleming/dsmcfoam/)还是另外一个MNF版本的求解器(https://github.com/MicroNanoFlows/OpenFOAM-2.4.0-MNF)。

[kuangyonglin]

@jeffx 在 DSMC时间平均 中说：

hyStrath版本

[jeffx]

@kuangyonglin 好的 多谢大佬

[jeffx]

@kuangyonglin 大佬如果不考虑化学反应的话，分子性质需要定义红框内参数吗？因为我看官方tutorial里supersonicFlatPlate案例里没有定义那些参数，因为我需要模拟水分子，很多参数查不到。

[kuangyonglin]

@jeffx 在 DSMC时间平均 中说：

@kuangyonglin 大佬如果不考虑化学反应的话，分子性质需要定义红框内参数吗？因为我看官方tutorial里supersonicFlatPlate案例里没有定义那些参数，因为我需要模拟水分子，很多参数查不到。

不需要，你可以看一下bird的书，也可以扒一下代码，不考虑化学反应的话，看你选择的粒子模型，根据模型的需要定义相关参数就可以了