@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++知识严重不足，查了问题都看不懂，需要自己补补基础了。谢谢东岳老师的回复！

确实有很大改善，可能上方空气部分需要设置一个比较大的高度

@东岳 用interDyMFoam算潮流能水轮机，没有仔细看，计算速度大概有不到10%的提高

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？

解决了，把1换成1.0就好了，感谢 @麦迪文 师兄的帮助。

@东岳 是的VOF多相流，计算域三层，最上面是空气，红色是油，下面是水2020-10-22 21-10-44屏幕截图.png

@bestucan 你是行家呀，厉害厉害，终于知道原因了

如楼上所说，不指定就是六个自由度，重叠网格也不会出现网格变形的问题

晚上又测试了一下，在linux系统下是可以顺利的进行计算，谢谢东岳老师的帮助:xinxin: