在creatFields.H文件中看到了如下代码，没有查到相关解读，创建类fieldTable的对象fields，然后调用了add(形参)函数，我在multivariateSurfaceInterpolationScheme.H找到了该函数，感觉是往hashTable中插入该变量的相关信息？麻烦有了解的大佬指导一番:xinxin:

Info<< "Creating multi-variate interpolation scheme\n" << endl; multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields; fields.add(h); fields.add(rho);

另外，在multivariateSurfaceInterpolationScheme.H文件中发现如下的类派生关系，两组类的派生关系这么声明是什么含义？第一次见:136: ，还望大佬解读一下:

namespace Foam { /*---------------------------------------------------------------------------*\ Class multivariateSurfaceInterpolationScheme Declaration \*---------------------------------------------------------------------------*/ template<class Type> class multivariateSurfaceInterpolationScheme : public refCount { public: //- fieldTable class fieldTable : public HashTable<const GeometricField<Type, fvPatchField, volMesh>*> { public: fieldTable() {} void add(const GeometricField<Type, fvPatchField, volMesh>& f) { this->insert(f.name(), &f); } }; ...... }；

OpenFOAM 当然更配 paraview。paraview 有个数据筛叫 celltopoint ，相反的也有，pointtocell。转了之后再输出。
tecplot用的不多，不清楚。

@bestucan 学到了，感谢

@东岳 请问李老师您能帮忙看一下这个问题吗？之前也了解到您做了许多这方面的很好的工作。另外，我在openfoam里面试着对比了一下无滑移和自由滑移的速度边界条件，发现二者的差别非常小。
另外不知道openfoam里面是不是也按照式（2）的方式，通过接触角算出界面的法向量，以此来植入接触角边界的？鄙人目前还看不懂openfoam源码，只是简单的应用，惭愧。所以向您请教。

速度场
U.png 压力场
p.png 粘度场
nut.png

粘度边界：进出口calculated，机翼表面nutUSpaldingWallFunction，上下边界symmetry简化问题

计算收敛，domainminandmax监测稳定无震荡

求解器：default interPhaseChangeFoam，ofv1912

湍流模型：kEqn LES

粘度场异常原因？

好像是资源不足，
虽然是超算，但是你用了slurmd，可能分配到的资源不够，把算的内容改的特别简单试试

当然以上建立在你描述的操作都无误的情况下

@xpqiu 嗯嗯，非常感谢~:xiexie:

@bestucan 在 fvc::DDt是什么意思？ 中说：

@xpqiu 是啊，我的签名就声明了我就是来搞笑的:chitang: 滚来滚去，滚去滚来，薛定谔的滚蛋状态:mihu: 布朗滚蛋运动:mihu:

像您这样的能答上来的人没有出现，只有先搞笑一下，不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究，一般不愿意回答的问题要是有个错误回答，就会引来很多纠错。因为纠错的成就感更强一些。

:146:

:146: :146:

Thanks a lot. The problem was solved.

@东岳 初步判定可能是of2.3.1版本对应的mpi环境配置问题，当输入mpirun命令后会提示：It looks like opal_init failed for some reason; your parallel process is likely to abort. There are many reasons that a parallel process can fail during opal_init; some of which are due to configuration or environment problems. This failure appears to be an internal failure; here's some additional information (which may only be relevant to an Open MPI developer):
然后具体解决办法在摸索中

搞定，输出叶片上的压力数据（不要选择internal mesh），用slice切一下，最后用plot data即可

@bestucan 太感谢了:xiexie:

@千里之行 求问大神 如何改变状态方程 编写新的psi,比如说我的状态方程是dac29137-1ecb-4648-88ec-e27adc2aac6e-image.png 那我的psi应该是什么呢

@random_ran OK 感谢! 很好的分享！

求问大神们，我计划使用rhoCentralFoam来模拟液体的弱可压缩性和热膨胀性，状态方程准备使用 0f46d7a5-4daf-440b-821b-24f146dbbe63-image.png，包含压缩系数和膨胀系数，为常数，这个表达式跟Boussinesq假定有点像 。但是，rhoCentralFoam只允许psiThermo类型，更具体地说是理想气体(称为perfectGas)和Peng Robinson EOS。1.我这个换状态方程想法可以实现吗，都需要改哪些文件呢？2.我知道Boussinesq假定是在rhoThermo里面，psiThermo和rhoThermo本质区别是什么呢，为什么rhoCentralFoam不能调用rhoThermo？3.psiThermo里面我知道FOAM是假设了rho = psi * p，psi是可压缩性为1/RT，相当于是理想气体，如果是我的这种形式的状态方程那这个psi应该是什么呢？

根据横向力做傅里叶变化即可求出频率

@BlookCFD 解决了，受到下面代码的启发，即重复利用代码，更新边界名称即可

p { solver GAMG; tolerance 1e-6; relTol 0; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 100; mergeLevels 1; } pFinal { $p; relTol 0; } 最后的要实现结果如下图： forceCoeffs_object { // rhoInf - reference density // CofR - Centre of rotation // dragDir - Direction of drag coefficient // liftDir - Direction of lift coefficient // pitchAxis - Pitching moment axis // magUinf - free stream velocity magnitude // lRef - reference length // Aref - reference area type forceCoeffs; functionObjectLibs ("libforces.so"); //patches ("body1" "body2" "body3"); patches (C0); pName p; Uname U; rho rhoInf; rhoInf 1.0; //// Dump to file log true; CofR (0.0 0 0); liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1); magUInf 1.0; lRef 1.0; // reference lenght for moments!!! Aref 2.0; // reference area 1 for 2d writeControl timeStep; writeInterval 1; } //监测新的边界条件，只需要重复利用代码，更新边界名称即可 C1 { $forceCoeffs_object patches (C1); }

感谢李老师点评！

@tidedrinker 感谢分享！:xiexie:

@Tian 初始化一定高度的颗粒堆积在流化床底部试试？现在这个设置，侧面全是周期边界，入口固定速度，假设是单相流情形，因为没有壁面剪切力，速度场会是完全均匀的。针对多相流情形，因为你入口设定的固相体积分率也是固定值，所以这个流动里面一切都是那么均匀，没有任何会引起不均匀的因素，所以算出来是均匀的也可以理解。