最近在用OpenFOAM和SU2提交计算时，反复出现"Authorization required, but no authorization protocol specified"报错信息。尽管不理会它最终还是能完成计算，但总是出现的报错信息让人心里不安。
怀疑是OpenMPI的问题，通过以下命令验证：
mpirun -np 2 hostname
果然出现报错信息：
翻阅网上信息，比较靠谱的解决方案是这个1：在Slurm脚本中增加以下环境变量：
export HWLOC_COMPONENTS=-gl
测试一下，报错问题解决了：
前期编译安装的是com版的OpenFOAM，版本号是v2412，没有附带cfmesh的源码。根据官方文档的提示，需要手动下载cfmesh的源码文件：
cd $WM_PROJECT_DIR
git submodule update --init --recursive plugins/cfmesh
出现以下错误：
fatal: 不是 git 仓库（或者任何父目录）：.git
无语了，直接通过git下载源码到指定文件夹，文件夹和URL路径可以查看.gitmodules文件：
cd $WM_PROJECT_DIR
git clone https://develop.openfoam.com/Community/integration-cfmesh.git plugins/cfmesh
进入cfmesh文件夹，开始编译：
cd plugins/cfmesh
./Allwmake -jN # N替换成CPU核心数
运行pMesh，验证是否安装成功：
pMesh -help
OpenFOAM的ccm工具包含ccmToFoam和foamToCcm，前者用于将STAR-CCM+输出的ccm格式网格转换成OpenFOAM的网格，后者相反。需要先编译安装第三方库libccmio：
cd $WM_THIRD_PARTY_DIR
# 下载libccmio
wget http://portal.nersc.gov/project/visit/third_party/libccmio-2.6.1.tar.gz
# 备选链接
# wget https://sourceforge.net/projects/foam-extend/files/ThirdParty/libccmio-2.6.1.tar.gz
tar xvzf libccmio-2.6.1.tar.gz # 解压下载的压缩包
./makeCCMIO # 运行libccmio的编译程序
接下来编译ccmToFoam：
# 先编译libccm
cd $WM_PROJECT_DIR/src/conversion/ccm
./Allwmake
# 然后编译ccmToFoam和foamToCcm
cd $WM_PROJECT_DIR/applications/utilities/mesh/conversion/ccm
./Allwmake
运行ccmToFoam，验证是否安装成功：
ccmToFoam -help
本文采用不同CFD求解器，对层流流动问题进行验证。通过列举不同求解器的操作过程和输出结果的差异，验证各个求解器的精度。
如下图，建立两个同心圆柱间的定常层流模型。流动由内筒以恒定角速度旋转而引起，而外筒保持静止。使用周期性边界，只需要对流域的一部分进行建模。物理模型和输入数据如下表所示。
| 来源 | Ansys验证算例 |
|---|---|
| 参考文献 | F. M. White. Viscous Fluid Flow. Section 3-2.3. McGraw-Hill Book Co., Inc.. New York, NY. 1991. |
| 物理模型 | 层流，旋转壁面 |
流体物性参数、几何尺寸和边界条件如下表所示。
| 特征 | 单位 | 参数 |
|---|---|---|
| 流体密度 | kg/m^3 | 1.0 |
| 流体粘度 | kg/m-s | 0.0002 |
| 内圆柱半径 | mm | 17.8 |
| 外圆柱半径 | mm | 46.28 |
| 内圆柱转速 | rad/s | 1.0 |
本文所用求解器和输入模型如下。
| 求解器 | 输入模型 |
|---|---|
| Ansys Fluent 2020R2 | VMFL001_rot_conc_cyl_2D.msh |
| Ansys CFX 2020R2 | VMFL001_rot_conc_cyl_3D.msh |
| Siemens STAR-CCM+ 2020.2.1(15.04.010-R8) | VMFL001_rot_conc_cyl_2D.msh |
| OpenFOAM v2006 | VMFL001_rot_conc_cyl_2D.msh |
| SU2 v7.1.1 | VMFL001_rot_conc_cyl_2D.cgns |
打开Fluent，选择2D求解器。网格数量较少，只用1个处理器核心运行。