也是实测。OpenFOAM自带算例。256万网格。虚拟机16核跑了1772s，真系统跑了548s。

虚拟机大约要慢3倍

@队长别开枪 专业 学习了学习了:140: 我看朋友用blender做过一些美术动画，看起来还是挺舒适的:papa: 目测不太适合工程:xinlei:

@同学博 是的，基于 freeCAD，专门给 OF 开发的前处理工具

@李东岳 谢谢李老师！我发现问题了，这个报错其实就是最后的correctADDStress的问题，所以看这个报错提示就是看最后一个信息，错误的原因是我用wclean去移除编译的库的时候其实没有移除，所以后面不管我怎么改，都是调用的原来的库，解决办法就是手动删除openfoam-7的platform中生成的自己编译的库，重新编译就好了:136:

可能是超算上有不只一个节点，所以节点间通过IP地址交流信息。
而linux系统里有一套“权限”和“环境变量”的继承规则。你加载 openfoam 环境变量的的那句应该是放在家目录的.bashrc了（就是source ~/OpenFOAM/etc/bashrc，大概是这句）

而 miprun 再通过IP调用其他节点的时候，这个网络通信的用户权限和环境变量，应该没有继承你的用户名的环境变量；而是系统的环境变量。所以不包含你的用户目录下.bashrc里的东西。

把你家目录那句 source 放到 /etc/profile 试试 。这是系统范围的环境变量。得是绝对路径：/home/<用户名>/OpenFOAM/……

当然也可能得在各个节点上都配上环境变量、装上 openfoam，这细致活得慢慢调……

部分求解器出现上述状况（编译成功，但是不能算），还有一部分求解器完全正常。

没用过多节点的，但问题是出在环境变量上没跑了。
到这个路径下
~/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin
这都是可运行的二进制文件。

直接./<文件名>，比如
./icoFoam
就知道是没编译好还是环境变量没配置好了。

一年多了还是没有解决这个算例，不知道有没有人做过相类似的，重新描述下这个算例，如图所示一个后台阶流动，但是是在图中左下（0,0,0）位置为原点绕Z轴旋转，主要流动参数，
雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$，其中$U_b$是体积平均速度；
旋转数$U_b=\frac{\Omega(H-h)}{U_b}$，$\Omega$是角速度；
台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$，这个与槽道流DNS结果是完全一致的，也就是说台阶前是完全发展的槽道流。目前存在以下问题：
1、求解器选择，SRF类还是需要自己修改求解器？
2、合理的进口条件，不旋转的相对好给，但是旋转的case，速度进口边界条件和初始值不知道如何给合理。
3、无量纲化的流向速度为什么会小于1，如下图所示。以我的经验，槽道流里面最大速度比体积平均速度肯定是大的，也就是$\frac{U_{max}}{U_b}$是大于1的，但是为什么在旋转流动里面就会于1呢？
6c444575-46fa-460b-bb57-8d655a396d45-image.png
参考文献
1、DNS原始文献：Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation
2、图片来自：Development and Application of Rotation and Curvature Correction，是篇博士论文

http://c.biancheng.net/view/2226.html

this就是当前对象的别称，有哪些成员，就看它所在的函数是哪个对象的成员函数了。
去定义这个对象的头文件查，与之同名的.c文件有各成员函数的定义。

感谢回复，我研究研究。我把体积分转化成了面积分，然后用python calculator进行了求和。感觉结果不太理想

@dzw05 你遇到的极值局部点距离给定压力条件的边界接近吗？我今天下午算的时候发现时间步过大时，距离fixedValue p_rgh边界最近的cell中心也会出现温度极值....调小了时间步后就没有了。我算的这段时间一直感觉温度和压强的联系太密切了...

不好判断 但感觉就是设置的问题 如果可以上传算例并且网格数量小的话 我可以给你看看

电势是相对的，上下机板电势差不变，就像电池正负两极保持电势差一样。不影响是否绝缘。
电荷密度零梯度，仅在极板与空气交界线上。看这种情况：紧挨着交接线两边各有10个电荷，所以零梯度。但是一边诗歌电荷后面站着100个电荷，另一边十个电荷后面站着20个电荷。

不同于流场的变化，仅分析相邻区域即可。电荷之间的相互作用是“超距”的，不相邻也可以用场去作用。
总的来看，电极板之间已有电势，只要有载流子就有电流。就像在电容间加水，就存不住电势差了。
交界处电荷梯度为零可以两边同时增加或者减少，或者增加减少动态平衡。

如果要绝缘，按照物理意义的绝缘最靠谱：没有载流子。
如果强行绝缘，应该是电势梯度为0，没有电势差，就没有电荷受力。但是电势场是场内所有电荷分布造成的结果，应该不能直接设定。

具体怎么实现得看源代码是怎么算电场的了

@veen 多谢大佬！^O^/

@史浩 你好，抱歉，之前居然没注意到这个回复，刚刚才发现。现在“导出该球表面积随时间的变化”这个问题应该解决了吧？
两种方法，一个是save data输出成csv，多个csv合并成一个文件，在excel里处理数据。另一个是plot data over time 也可以

另外，请问，表征形状是不是好的球体这个问题，解决了么？
如果是少量的液滴，可以求一下质心坐标的偏差，但是你这个图的液滴数量比较多，取坐标什么的应该很麻烦。所以我很好奇这个问题有什么解决办法么？

@史浩 您好，请问对于自适应网格的设置讲解有什么资料可以帮助我学习参数设置的吗

![Y`IW6LQFCDZXZ5RHH%K69V.png**图为剖面图，两部分的网格没有connect，这里有两个内部面，如何将其缝合使网格节点对应
**

@赵一铭
你好，请问对于并行计算的案例该如何用gnuplot实时监测残差呢？

@红豆沙 fluent meshing画的多面体of是支持的，多面体和六面体核心结合的新型网格形式of不支持，不知道新版的有没有更新

回来填下坑
为了在dragModels中引入一个新的区间判断函数between()，除了需要在

\src\OpenFOAM\primitives\Scalar\Scalar.H

中定义函数主题，还需要将函数进行重载，这是因为写入Scalar.H的文件只是单纯意义上的scalar类型，并不是计算所需要的

Foam::tmp<Foam::GeometricField<double, Foam::FvPatchField, Foam::volMesh

类型，为了使函数返回的类型符合所需数据类型需要对函数进行重载操作，具体操作的方法是在下列地方写入重载的定义，分别是

\src\OpenFOAM\dimensionedTypes\dimensionedType.H \src\OpenFOAM\dimensionedTypes\dimensionedScalar.H \src\OpenFOAM\dimensionedSet\dimensionedSet.H \src\OpenFOAM\fields\DimensionedFields\DimensionedFieldFunctions.H \src\OpenFOAM\fields\FieldFields\scalarFieldField.H \src\OpenFOAM\fields\GeometricFields\GeometricField\GeometricFieldsFunctions.H \src\OpenFOAM\fields\GeometricFields\GeometricScalarField\GeometricScalarField.H

重载之后可以使between()函数返回的数值类型为

Foam::tmp<Foam::GeometricField<double, Foam::FvPatchField, Foam::volMesh

之后重新编辑一下就可以运行了

@Foamer24 你好可以分享下您关于setSet的具体操作吗

@李东岳 没用过。因为要做耦合 所以用的foamextend4.0

楼主搞定了吗，我最近分析一个圆柱加速运动，也想计算一下其附加质量力

好像有个wsl和wsl2，我也没太研究，我直接关了..

@李东岳 老师好 我用的这个求解器虽然基于OF 不过还是有些差别 我用的是该求解器自带的网格处理工具，做的是结构网格，由“点→线→面”生成不同的block，然后block组合成整个几何，您说的交界面我不太理解是什么意思。网格代码比较长，如果老师对该求解器划分网格的方式感兴趣的话我可以把网格代码发上来。

这是我的一些参数选择 双破折号后面的为注释
代码.png

另外像您说的，壁面去掉粘性后，计算结果正常多了
无粘性密度-tuya-tuya.png

忽略左下角哈哈哈 那个block的壁面我忘记去掉粘性了:chouchou:

@bestucan 谢谢你~ 我找到问题了，就是ADD范围不对，我用了最值函数限制了就好了，这一点论文里面没有提，很不好找，谢谢你~

@李东岳 谢谢李老师，我再详细看一下推导过程

@史浩史老师，help，帮忙看一下呗:140: https://www.cfd-china.com/topic/4247/流电耦合模型中边界条件的设置?_=1604991413128

或者如果不能导出网格连接关系的话怎么导出网格面积以及法向量呢

@Calf-Z-DNS 好久没联系 你的DNS算的怎么样了

还是看CFD文献有意思 回归正轨 回归正轨:mianmo:

有偿请教，问一下有关拉格朗日颗粒喷射的问题，我往喷射模型（manualInjection模型）动态的传递喷射点坐标，Positions文件坐标每一个时间步都会发生变化，怎么实现时刻读取新的喷射点坐标并往出喷粒子呢？@东岳老师

@benqing 谢谢你的解答！我再想想其他解决办法

@hotairballon1997 看着我真诚的眼神，告诉我，你看了我的帖子

or per-case override in the $FOAM_CASE/system/controlDict file

你懂我的意思吧.jpg

具体是非固体边界第一层的网格的yplus与uplus是怎么计算的。

就按照壁面的计算的。我那个图出的非常麻烦，没有代码一键搞出来，还得结合paraview导出数据，所以就没放。反正算法就是跟壁面的一样，你看看你能搞出来个一键运行不 Good luck :baobao:

请教一下各位前辈，用OF算完旋转机械之后，利用paraview进行后处理，在画流线的时候发现在动网格区域显示的是相对速度的流线，怎么才能让它显示绝对速度曲线呢？

@Marcllory 对，氢气才两个质子的大小。空气中氧气、氮气、二氧化碳二三十个质子中字。当空气是个窟窿眼很大的网，氢气在里面穿梭。

如果想确定这种简化误差有多大，可以关注这两个数
Screenshot_20201106_221545.jpg
Screenshot_20201107_095249.jpg

忽略的就是氢气的密度（动量运输）和氢气浓度差造成的扩散（扩散速率）。如果这两个占的比重跟对流比起来很小那就可以忽略。

以上都是我猜的，我没做过这个:chigua2:

有偿请教，问一下有关拉格朗日颗粒喷射的问题，我往喷射模型（manualInjection模型）动态的传递喷射点坐标，Positions文件坐标每一个时间步都会发生变化，怎么实现时刻读取新的喷射点坐标并往出喷粒子呢？

@刘雄国 赞，太感谢了。我赶紧去补一下vector<>:146:

这个问题感觉和我前几天提的问题一样，我已经解决了，原因就是低粘度液体在半径大的管道中就会出现这种速度分布的平台状。
解决办法有两个：提高流体的粘度，或者减小管道的半径尺寸。
你看一下两个模型的尺寸设置是不是一样的或者流体粘度设置的是不是一样。

tahts difficult to answer, I would suggest to get to know the algorithm first, if you want to know the settings, the knowledge should come to the first place. So, get some books to have a look at the algorithm, afterwards the settings shouldn't be a problem.

@东岳 是的，粘度调高也可以出现抛物线，感谢 感谢。我把网格尺度变小了，也能观察到。

:143:

我的上述计算中没有用到估算的方法，那论文里的是不是因为边界层网格采用的估算的行为，所以才会导致friction Reynolds numbe因人而异，出现微小的差别呢！

@东岳 李博，我准备试试，之前写过fluent多面体的网格转换程序，不过这个发不了论文:mihu:

@fangyuanaza :146: 非常感谢！

@fubianhanshu 好的我自己试一下，谢谢各位前辈 @星星星星晴 @fubianhanshu

例如下图所示的波浪-圆柱体耦合计算

波浪圆柱耦合.png

根据理论波浪力方程，圆柱所受到的波浪力（水平方向、垂直方向）为

理论波浪力.png

Fy表达式中的最后一项ρgV0为圆柱体的静水浮力（initial buoyancy in still water）；

请问：OF输出的Fy（垂直方向上圆柱受到的Force）里面包含ρgV0这部分吗？想要获得圆柱受到的垂直方向上的波浪力，是不是需要从OF获得的计算结果中手动减去ρgV0这部分？

@yfclark 谢谢！

十分感谢！@东岳 :146:

http://dyfluid.com/icoFoam.html 另外 这个例子你看过么 最后那个扰流

谢谢李老师一针见血的点评，新版本的komegaSST模型是有一个共同的Base的，然后衍生出不同的komegaSST湍流模型，komega湍流模型也是这样的结构，我修改的时候这部分没有处理好。之前觉得就是给动量方程加一项，不懂C++应该也不难，现在发现C++知识严重不足，查了问题都看不懂，需要自己补补基础了。谢谢东岳老师的回复！