OpenFOAM当中监测力和阻力系数

首先准备好我们自己的平常算例文件，本次我们以圆柱绕流的算例来说明用法

我们找到constant文件夹

打开其中的transportProperties文件

我们将其中的：

nu              [0 2 -1 0 0 0 0] 0.01;

修改为：

nu              nu [0 2 -1 0 0 0 0] 0.01;

的形式，然后保存

接下来我们找到system目录并打开

打开其中的controlDict文件

在文件后面追加一下语句（每条语句的含义见语句后面的注释）

functions

{

//力的输出

Forces

{

type forces; //名称

functionObjectLibs ("libforces.so"); //包含共享库

patches (CYLINDER); //此处填写需要监测的边界名称，如果有多个边界采用(patch1 patch2 ...)的形式

log true; //选择是否输出log文件，这里我们输出，其实影响不大

rho rhoInf; //参考密度，对于不可压问题，直接填写流体的密度即可

rhoInf 0.001; //这里填写流体密度的数值

CofR (0 0 0);

liftDir (0 1 0); //升力的方向

dragDir (1 0 0); //阻力方向

writeControl    timeStep; //按照timeStep的方式输出，其他形式，请参考OpenFOAM用户手册

writeInterval   5; //间隔我们设置每5步输出一次

}

//阻力系数的输出

forceCoeffs

{

type forceCoeffs; //名称

functionObjectLibs ("libforces.so"); //包含共享库

patches (CYLINDER); //此处填写需要监测的边界名称，如果有多个边界采用(patch1 patch2 ...)的形式

log true; //选择是否输出log文件，这里我们输出，其实影响不大

rho rhoInf; //参考密度，对于不可压问题，直接填写流体的密度即可

rhoInf 0.001; //这里填写流体密度的数值

CofR (0 0 0);

liftDir (0 1 0); //升力的方向

dragDir (1 0 0); //阻力方向

pitchAxis (0 0 1);

magUInf 1; //参考速度

lRef 2; //参考长度

Aref 3.14; //参考面积

writeControl    timeStep; //按照timeStep的方式输出，其他形式，请参考

OpenFOAM用户手册

writeInterval   5; //间隔我们设置每5步输出一次

}

};

到这里修改文件的配置就算是结束了，接下来我们按照正常流程开始求解，本次我们用icoFoam求解器，输入命令：

icoFoam | tee log

在算例的目录下我们可以看到多出了一个目录postProcessing目录：

然后我们打开该目录，会发现有多出了两个文件夹，分别名为Forces和forceCoeffs（是我们刚才在controlDict下面命名的）

其中Forces文件夹就包含着监测升力和阻力的数据，而forceCoeffs包含相关力的系数：

原视频下载地址为：https://yunpan.cn/cv2FtdIwUvsjT （提取码：9391）