OpenFOAM常用linux命令

1 删除所有计算结果文件

• OpenFOAM自带命令

  foamListTimes -rm
  

  删除除了0文件之外的结果文件。

  如果0文件也要删除，则添加 -withZero 选项：

  foamListTimes -withZero -rm
  

  另一些方法则利用一些linux自带的命令

• 方法二：通配符

  rm -rf [:digit:].[:digit:]
  rm -rf [0-9].[0-9]
  

• 方法三：花括号扩展

  rm -rf {0..9}.{0..9}
  

• 方法四：删除所有求解产生的文件。假定求解前改动的最后一个文件是 controlDict，那么先找到所有比controlDict时间更新的文件，然后删除之。

  find . -newer controlDict -type d -exec rm -r '{}' +
  

2 进程相关

jobs -l # 查看作业列表
bg %<n> # 把序号n的作业放到后台
fg %<n> # 把序号n的作业放到前台
Ctrl + Z # 把前台正在执行的作业放到后台并暂停它
paraFoam & # 在后台运行paraFoam
icoFoam & # 在后台运行icoFoam

当用户退出远程系统的时候，nohup 会让进程继续运行。-n 是 nice 的选项，可以设置为19，最低优先级。> log 表示把结果输出到 log。& 表示后台运行 icoFoam。

nohup nice -n 19 icoFoam > log & 

3 查找

3.1 查找可用字段

比如说我想知道压力场可以用什么单位，此时可以用 foamSearch 查找教程目录下所有 p 文件的 dimension 字段：

foamSearch $FOAM_TUTORIALS/ p dimension

dimensions [ 0 2 -2 0 0 0 0 ];
dimensions [ 1 -1 -2 0 0 ];
dimensions [ 1 -1 -2 0 0 0 0 ];
dimensions [ m^2 s^-2];

如果指定 -c ，可以查看这些字段出现的次数。例如，我想查看 fvSchemes 中的 gradSchemes 有哪些写法。

此时可以用：

foamSearch -c $FOAM_TUTORIALS fvSchemes gradSchemes

如果我们想进一步查看 grad(P) 的，可以在 gradSchemes 后面加个斜杠，再加上 grad(P)，即 gradSchmes/grad(P)。不过此时需要注意，因为字段里面有括号，属于 shell 的元字符，所以必须用引号包裹这个字段：

foamSearch -c $FOAM_TUTORIALS fvSchemes "gradSchemes/grad(p)"

 13 grad(p)         Gauss linear;1 grad(p)         leastSquares;

我们看到，大部分都是 Gauss linear，只有一个是 leastSquares。如果我们想找出这个leastSqures 的特例，就需要查找这个文件。

3.2 查找文件

接续前文，我们想找到教程所有fvSchems 文件中的 gradSchems/grad(p) 为 leastSquares 的是哪个文件，它完整路径是什么。

find $FOAM_TUTORIALS -name fvSchemes | xargs grep -El "grad\(p\)\s+leastSquares"

/home/jasminefeng/OpenFOAM/OpenFOAM-9/tutorials/incompressible/nonNewtonianIcoFoam/offsetCylinder/system/fvSchemes

xargs 不能少。比如 | xargs 前的结果是 file1 file2，本来后面的语句应该解析成 grep -El "grad\(p\)\s+leastSquares" file1 file2，file1 file2 是作为参数的。但没有 xargs 的时候，后面 file1 file2 就不能作为参数了。注意：args 就是参数的英文。

再解释一下 grep 命令。这里用 -E 表示使用正则表达式，引号中的\s+ 表示一个或多个空格。grep 的 -l 选项表示列出所有含关键字的文件名（还有一个孪生选项 -L ，表示列出所有不含关键字的文件名）。

3.3 查询命令

我们想找到一个命令（例如 mapFields）的用法，除了使用帮助命令：

mapFields -help

还可以使用 foamInfo。这两个是有区别的。-help 提供的主要是 Usage

而 foamInfo 提供的是命令的源文件位置、简述和使用案例。罗列出的案例默认情况下不超过十个，我们可以用 -a 来列出所有案例（不过这里 mapFields 只有4处使用了）。

不过我很好奇 foamInfo 是怎么判断这些案例是否确实使用了 mapFields 的，难道是在哪里记录了吗？我扒了一遍 foamInfo 的源代码后，发现它原来是根据 Allrun 中有没有 mapFields 这个命令来判断的。

Allrun 就是把运行一个案例所需要的所有命令都集成在一起的一个文件。如果这个里面有 mapFields，那就说明它要使用，就可以参考这个案例。原来如此。

P.S. VSCode 有一个插件 Bash Debug 可以用来对 shell 脚本进行断点调试。在 launch.json 中添加 args 可以设置运行命令的参数。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-HkDMNI4h-1680234744104)(https://fengimages-1310812903.cos.ap-shanghai.myqcloud.com/20230331110626.png)]

3.4 查找字典文件

我们知道 OpenFOAM 的命令大都依赖于一些字典文件，这些字典文件是用户可编辑的，有点像 OpenFOAM 的前端。那我们如何找到这些字典文件呢？为此 OpenFOAM 提供了 foamGet 命令来获取这些字典文件。在使用 foamGet 时我们必须确保当前目录下有一个 system 的子目录。例如我们想获取 mapFieldsDict 的文件，则

foamGet mapFieldsDict

Multiple files with “mapFieldsDict” prefix found:

1. /home/jasminefeng/OpenFOAM/OpenFOAM-9/etc/caseDicts/annotated/mapFieldsDict
2. /home/jasminefeng/OpenFOAM/OpenFOAM-9/etc/caseDicts/preProcessing/mapFieldsDict
   ** Note: it is easier to use files NOT in the “annotated” directory

Enter file number (1-2) to obtain description (suggest 2) :

此时命令告诉我们找到了两个 mapFieldsDict 文件，并且还贴心地提示我们 annotated 目录下地文件不太好用，然后问我们选择哪个，那我们就选择 postProcessing 目录下的，然后键入 2，此时该命令就会把这个文件复制到 system 目录下。

Copying /home/jasminefeng/OpenFOAM/OpenFOAM-9/etc/caseDicts/preProcessing/mapFieldsDict to system

查看它的内容，可以看到这其实是一个字典的模板，被尖括号包围的内容就是我们需要替换的。

4 复制算例

foamCloneCase -latestTime cavity cavityHighRe

指定 lastestTime 就会拷贝时间最新的文件夹（如 0.5，以作为初始条件的文件夹）。如果不指定 lastestTime 就会拷贝0文件夹。

5 把日志同时输出到终端和文件

用 tee 命令

icoFoam | tee log.txt &

在 $WM_PROJECT_DIR/bin 下还有不少好用的小工具，以后找机会再介绍~