0_1515589858123_捕获.JPG

采用相同的求解器，上面是并行计算，下面是单核计算，计算结果有偏差。

感谢关注，过年的时候应该可以 :cheeky:

在自己的笔记本上跑算例的时候，如果打开Chrome，经常由于内存不足导致网页自动崩溃。回国赶紧换工作站。

我觉得可能是win10系统占用的内存更高了，或者是新版本Chrome内存占用过高。

@李东岳 我最后还是把这个打开了。我模拟的是潮流能水轮机，开之前功率系数不过30%出头，开了之后到了45%左右（TSR=6），功率系数-叶尖速比曲线基本上与实验对上了。而且开了之后对时间步长就不敏感了，时间步从0.0001s到0.002s，CoN从40多一直到700多，功率系数变化在5%以内。

而且，OpenFOAM/tutorials/multiphase/interPhaseChangeDyMFoam/propeller这个算例也是开启了动量预测，我考虑是不是固定转速的旋转机械都需要开这个选项。

@流卜 的确是接口差别很大，目前准备把extend上的部分内容移到v1706，感觉崩溃了快

@李东岳 0_1515462127425_2018-01-09 09-40-58屏幕截图.png
我按照您说的，我把整个算例搬过来，只是改了一下0文件里面各个变量，在system文件加了几个字典文件，结果又出现了上图的问题。我怕是我哪个地方，出错了，又重新操作了一遍，还是这个问题。我又去运行了一下原来的算例，没有出现任何报错。我看了一下，除了data这个我没找到，其他的都能找到。

对于单气泡的曳力，升力更倾向于是一种面力。但是由于这个面不是网格单元的面。处理起来存在一定困境。我们把这些力当做体积力？有一些模拟结果表明存在震荡。插值做面力（类似基金会这种处理）？但是结果是错的。

需要深入研究力的处理方式。

可以使用inletOutlet或者outletInlet，满足你的要求。

非常感谢 @程迪 的讨论，我重新回顾一下看看能不能分享一些看法。

有关ddtPhi()这个函数的形式大家已经明白了，更重要的是明白为什么要这么做。MULES限制器和ddtPhi都是Henry Weller发明的，MULES限制器的植入以及原因目前都清楚了。但是ddtPhi的原因尚不明了。

@dyj19901127
宇老师大作被翻出来了，哈哈

我大略看了一下：

#include "createFluidMeshes.H" #include "createSolidMeshes.H"

划分两个网格的时候就存在了边界场，在求解fluid之后，边界场信息会更新，你不需要显性的传递边界场。

fvMesh& mesh = solidRegions[i]; solidThermo& thermo = thermos[i]; tmp<volScalarField> trho = thermo.rho(); const volScalarField& rho = trho(); tmp<volScalarField> tcp = thermo.Cp(); const volScalarField& cp = tcp();

这几行代码自动包含边界场，其中给固体的是网格共享的边界场信息。

应该是这个思路， 你可以确认下。

@程迪
rhoE应该是$\frac{\partial \rho e}{\partial t} + \frac{\partial \rho K}{\partial t}$
phiEp应该是$\nabla \cdot (\rho \mathbf{U} e)+\nabla \cdot (\rho \mathbf{U} K)+\nabla\cdot(p\mathbf{U})$
sigmaDotU应该是$\nabla \cdot(\tau \cdot \mathbf{U})$
感觉拆不了。

另外，请问对定常的求解器也有时间步限制吗？如果没有的话，是不是定常的求解器中，比如buoyanySimpleFoam，就不用管隐式还是显式了？

@lxjcfd
同求第三篇文献，在网上确实找不到。

@liuyf
截面jump是 p(x=0) == p(x=1) + dp
截面两侧单位jump 是 p(x=0+dx/2) == p(x=1-dx/2) + dp

边界网格越粗，这个差异应该会越大吧。

楼主现在有没有做过delta的测试？

@李东岳 谢谢东岳老师，明天再试一下

@搬运工不好当 您好，请问您对cavitatingFoam研究的怎么样？可以指教一下吗？

你可以先用topoSet的nearestToPoint生成pointSet，再由pointSet生成faceSet，在生成的若干面里面剔除无关的面即可。

@李东岳
工程上有意义的是达到指定误差的计算代价，一般理论上认为误差可以从某种定义的残差推算出的不确定来衡量，时间以taubench用25000网格迭代10步作为1个时间单位来衡量。

这样的假设下做横向比较是可能的。但是，具体情况嘛，呵呵。湍流模型的不确定度就说不清楚。

@yhdthu 好的，多谢您的回答

我猜测你这个是源代码没写对，或者拷贝复制的时候丢了代码，undefined reference这种错误只是在编译类后捆绑到求解器的时候函数没定义导致。可能你代码里面有些函数代码块丢失了。

以上只是猜测，因为不清楚你安装的详细过程不能确诊。