OpenFOAM显示残差

本文主要讲解两种方法用来显示OpenFOAM的计算残差，一种是采用OpenFOAM自带的foamMonitor来输出残差，另一种就是大家经常看见的采用pyFoam来输出残差。不管采用哪一种方法都必须安装gnuplot。

安装gnuplot的方法为：在Linux（以乌班图为例）的终端中输入

sudo apt-get install gnuplot

方 法一：采用自带的foamMonitor来输出残差我们首先需要将OpenFOAM中输出残差的文件拷贝到算例的system目录下，我们可以采用以下命 令在OpenFOAM的安装目录下查找该文件

find /opt/openfoam4/ -name residuals

然后我们使用cp命令将查找到的residuals文件拷贝到算例system文件夹下

cp “/opt/openfoam4/etc/caseDicts/postProcessing/numerical/residuals”  “/home/dell/new/system”

接着我们在算例的system/controlDict文件中添加后处理函数

functions

{

#includeFunc            residuals

}

文件为：

FoamFile

{

version     2.0;

format      ascii;

class       dictionary;

location    "system";

object      controlDict;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     simpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         20000;

deltaT          1;

writeControl    timeStep;

writeInterval   50;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

functions

{

#includeFunc residuals

}

我们运行算例

simpleFoam &

我们在该窗口输入以下命令：

cd postProcessing

cd residuals

cd 0

foamMonitor -l residuals.dat

上面那么繁琐的命令可以用用一句命令表示：

foamMonitor -l postProcessing/residuals/0/residuals.dat

接着见证奇迹的时刻到来了！残差出现了

​未尽事宜，请参考视频操作

原版视频下载地址：https://yunpan.cn/ckffDej3CXzQJ（提取码：1e27）

方法二：应用Python的pyFoam来输出残差

可参看博文：

http://blog.sina.com.cn/s/blog_5fdfa7e60100g4zw.html

http://blog.sina.com.cn/s/blog_4a0a8b5d0100s832.html

两篇博文均采用svn的方式安装pyFoam，但是不只是时间太长还是其他什么原因，本人无法通过svn的方式下载PyFoam，但我们可以下载已经编译为二进制的PyFoam包，采用以下命令安装pyFoam：

sudo pip install ./PyFoam-0.6.6-py2.py3-none-any.whl

pyFoam二进制包的下载地址为：

https://pypi.python.org/pypi/PyFoam

需要特别指出的一点是两篇博文在利用log日志文件进行残差显示时都采用的

pyFoamPlotWatch.py cavity.log，而正确的命令应该为：

pyFoamPlotWatcher.py cavity.log

未尽事宜，请参考视频操作

原版视频下载地址：https://yunpan.cn/ckff4DQRC5CiK（提取码：8e55）