CFD中文网

@Tong 好的，真知灼见。大概理解了下，您强调的是充分利用物理性质来构造预条件。。确实偏数学。。“针对体系结构和并行的优化可能不同”，的确如此。感觉如果是做优化的话，还是要针对特定的硬件架构去做实现层面的优化，想要在数学（算法）层面优化太难了。。。

@田畔的风 在 有人在mac上装of么？ 中说：

@一颗鸭蛋 默认情况下，你编译求解器生成的可执行文件会保存到$FOAM_APPBIN，比如我的默认在<上级文件夹>/OpenFOAM-v2206/platforms/darwinARM64ClangDPInt32Opt/bin，但是这个App的加载磁盘是只读的，所以会报错。

解决方案有两个：

修改在Make/files中你的求解器保存路径 EXE = $(FOAM_APPBIN)/xxxFoam -> EXE = <新路径>/xxxFoam

然后在环境变量中为$FOAM_APPBIN附加这个路径，以便让OpenFOAM在执行时能搜索到这个文件

export $FOAM_APPBIN=<新路径>:$FOAM_APPBIN 创建一个可读写的分区，直接编译安装。参考 https://github.com/BrushXue/OpenFOAM-AppleM1。

您的意思是把openfoam安装到创建的可读写分区吗？我用的也是M1的Mac，现在想在Mac上编译自己修改的求解器，大佬能不能给点详细的指导？谢谢🙏

@东岳 过两天我测试一下6142

@veen 好的。等你的教程，嘿嘿

@bestucan Thanks for your attention.

@东岳 在 讨论 Fluent计算效率优于OpenFOAM 中说：

你需要深挖一下

给你举个例子，下面是运算的一个log文件，时间步长非常小，2e-5，算2天之后，结果完全可以是正确的，但我觉得有问题

捕获.JPG

在自适应调节时间步长的情况下，时间步长会自动跳转满足CFL标准。其实对于我这个，如果在代码里加上几行，输出Ur的话，会发现下图中全场计算域中，就这么一个地方Ur的值特别大，导致局部库郎数过大。

捕获2.JPG

我建议你看看你的算例是不是也是这样，我怀疑跟我差不多

想请教一下遇到这种情况应该怎么办呢？特别算vof的时候，界面某些点出经常出现很大的速度，导致时间步长小的无法接受。

@yhdthu 现在很多实验研究都证实了这种压力波动，并且影响到空泡的周期行为

@李东岳 最近在跑作业，后面停歇时试试

@胡同小巷 你那个算例我下载了之后没法运行。好像不是openfoam8的。不过你这个压力的图，看起来比较正常。

@小火人 3年过去了，不知道你是否愿意分享你的这个研究算例？可以的话我们聊一下细节

最近把理论这部分弄了一下，

捕获.PNG

@东岳 等我今天晚上也跑一下，我也是9750H。
好像我的6142鸽了好久了。。。

@李东岳 我的模型两侧还有静止区域，只是这个简单模型还没加上，所以SRF是不适用的。

@Haoooo 已关注您，下面回复有我联系方式

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

@李东岳 好的，谢谢老师。

不好讲故事是比较困境了，哎..

@李东岳 非常感谢李老师，终于搞定，能输出数据了，但还需要进一步验证是否正确。差点放弃想升级新版本of。

volScalarField& summass = const_cast<volScalarField&>(p.mesh().objectRegistry::lookupObject<volScalarField>("summass")); const label cellIp = p.cell(); scalar& PPC = summass.internalField() [cellIp]; PPC += p.nParticle()*p.mass();

希望对以后需要的朋友有点借鉴作用。特别感谢 @星星星星晴 不厌其烦的解答:146:

@李东岳 好的好的，谢谢老师！

这段用来统计网格单元内粒子个数的代码，编译的路径是什么，求东岳老师指导