12.jpg
用重叠网格求解器，overPimpleDyMFoam计算了一个平板前移，在物面网格和背景网格之间有非常明显的界线。想请教各位大佬，生成网格的时候，内外层网格之间需要注意什么才可以让两层网格间的界限不那么明显。

@bestucan :146: 感谢感谢

@东岳 好的东岳老师！ :ok3:
最近也是刚开始试着改写求解器，顿感困难重重，确实需要多锻炼多积累。

只要你没有明确的格式化你的硬盘，你的数据就不会丢失。只是从一个带系统的硬盘变成了大U盘。
Ctrl + alt + F1~F7 都算可以调出来终端，一般前两个或最后一个是图形的。
即使上面七个都挂了，grub 里还有一个，Linux 没那么容易崩的:mihu:

这就是 linux 难推广的地方，很多经验并不相通，ubuntu-drivers 是 ubuntu 自己搞的，debian 系的并不通用，这是开始夹私货了。

@calrel是这样的，我的仿真时间是从0到50ms，计算时间间隔是e-6s，写出的时间间隔是e-4s，我想把这些写出的所有的时间步，从0.0001到0.0499，都重构处理，这样在paraView里就可以看到完整的模型场嘛，如果直接输入reconstructParMesh的话，报错如下图， dae0adad-a66c-4b4b-ac26-83d9091ef2f1-image.png 即使后面加了-case，-fullMatch，-noZero都一样的报错信息，跟上面一个大佬的回答不太一样，不知道你是不是也一样？

经过尝试发现这个reconstructParMesh指令需要指定单独的时间步（reconstructParMesh -time 0.001）进行操作，后面紧跟一个reconstructPar -time 0.001，然后再去重构下一个时间步，如此往复循环直至所有的时间步都被重构，不知道是不是都要这样来，比较麻烦，如下图所示b021e43b-7f85-49f0-b2e7-0cc3f7fc7e6e-image.png 0b1b8e86-1887-4b71-9483-ff6d8e918477-image.png9c139540-068d-45ac-8513-e1941ad04ede-image.png 随手截的图，时间没有对的上号

谢谢老师

在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 学到了，感谢

phi，phiHbyA是通量，通量的边界是固定的，不需要插值

@东岳 请问李老师您能帮忙看一下这个问题吗？之前也了解到您做了许多这方面的很好的工作。另外，我在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 初始化一定高度的颗粒堆积在流化床底部试试？现在这个设置，侧面全是周期边界，入口固定速度，假设是单相流情形，因为没有壁面剪切力，速度场会是完全均匀的。针对多相流情形，因为你入口设定的固相体积分率也是固定值，所以这个流动里面一切都是那么均匀，没有任何会引起不均匀的因素，所以算出来是均匀的也可以理解。

@dzw05 耦合计算和单域计算的流体域是使用上面提到的模型求解，求解器则是不同的，一个是壁面导电耦合求解器，一个壁面是绝缘求解器。我的耦合也只是壁面和流体的电势耦合。

谢谢老师

各位老师好，我想在OF里植入一个与流动直径有关的边界条件。通过sHM后得到的网格，其各流动通道直径不一且比较复杂，能否获取各网格所在流动通道的当量直径呢？在网上查阅wallDist(mesh).y()可以计算网格到壁面的最近的距离，但不知道能否用来计算出直径。

楼主能不能把算例分享一下，我目前三维的rhoCentralFoam基本都是能量方程发散，rhoSimpleFoam都是一步就发散了。。。

谢谢东岳老师指点，我刚刚测试过在空化流里面 R设置为7255或者3000之类的数值 对最后结果确实影响很小:chitang:

@RebelYoung http://dyfluid.com/rhoSimpleFoam.html 这面有一些rhoSimpleFoam的算例，可以先研究研究

@xpqiu 感谢。我试下你的方法

自己顶一下，有人临幸一下么？

@东岳 就是自己琢磨的，:xiezuoye: 写出来一是分享，二来有错的，大家能指出来，老师不用太在意，兴许是逻辑相当自恰的错误连篇:143: :chigua2:

四壁导电矩形管的对流传热case，现在想实现在其中一个壁面添加热源（实际上就是在其中一个流固交界面施加固定的热流密度q）。该怎么样设置边界条件呢？

@东岳 稠密的.sediment transportation.

老师们晚上好，我想问一下OpenFOAM里liquidEvaporation的模型是怎样控制水分蒸发的速率的？教程里面只指定了物质种类H2O，可不可以指定蒸发速率？非常感谢！

liquidEvaporationCoeffs { enthalpyTransfer enthalpyDifference; activeLiquids ( H2O ); }

谢谢东岳大佬:146:

一个网格一个粒子的话，$\alpha=\frac{v_i}{V_{cell}}$，一个网格多个粒子的话，$\alpha=\frac{\sum v_i}{V_{cell}}$，$v_i$是颗粒体积

在FLUENT里面用欧拉-拉格朗日求解时，欧拉方程里是认为α=1的，所以不适合用于高密度颗粒流的求解，而OpenFOAM中是考虑和计算体积分数α的。

理论上式对的，不过你要确认下在FLUENT里面用欧拉-拉格朗日求解时，欧拉方程里是认为α=1的？

各位大佬，请问无法并行计算，出现这个提示是怎么回事呀？

ERROR: unable to allocate launching block

IB Method的优势就是背景网格不用变动，可以节约很多dynamic mesh的资源。
如果想快速测试一下IB Method，可以考虑使用LS Dyna.

@马乔 好的，谢谢。 顺便问一下，这个auto是指什么？数据类型自动根据后面的指定嘛？

多谢，已更新:duang:

@bestucan 对的，我正在监测，我看了自带算例的残差，0-2秒残差下降，2秒这一刻突然集体升高，然后2-4秒一直上升，4-6秒上下波动趋于稳定。然后我现在主要就是两秒这个时候在那疯狂迭代，我两秒之前和自带算例一样，残差是下降的

internalField就是算例的初始值

@东岳 在 LES中 周期性边界条件的问题 中说：

segregationVector

李老师您好，这个关键词应该按照什么编译，麻烦了

@xpqiu 在 请教如何减小yPlus值，在简单的dam break算例里 中说：

@Sloan 多相流的不太好估计，只能尝试不同的值看看了。

好的好的，我再试试看，谢谢！