@李东岳 谢谢李老师！A fifth-order shock capturing scheme.pdf 李老师，我是想把这个文献里的第二部分模块数值算法植入进去，我昨天上传的是第二部分前面一块，后面那一块主要是进行重构，去计算q_L和q_R，然后通过前面的公式去进行计算，前面那块我看和您发的rhocentralfoam解析中很像，但后面的公式不知道应该怎么植入，没有思路，想请教一下老师:chouchou:

@李东岳 因为换个地方继续算，重新分块的时候就根据最后时刻的网格来划分，所以用scotch 就是每个块基本同样网格数（上下几十浮动），而用simple划就会是那种XYZ方向的拉升

lz解决了吗？我也遇到同样的问题了，高压气流冲击真空环境，使用rhoCentralFoam，算了一段时间之后就出现negative initial temperature t0，改了好久，可以推迟出现的时间，但总会出现，不知道怎么弄了

@李东岳 非常感谢东岳老师的回复，这个问题出现在均匀网格，有的时候甚至同样的位置，但是气域的速度不同，他最后采样的范围也不一样。。。

@turbulence 02525037-654f-4a52-bc51-dcef3cc03688-image.png

@李东岳
感谢李老师的回复，前段时间做了其他的事情，就搁置了这个。哈哈
这个问题解决了，是因为边界条件的设置。
另外请教李老师另外的问题：

OpenFOAM中有求解瞬态不可压多孔介质的求解器吗？ 或者关于二相流中造单纯的流，关于相界面处的边界条件设置问题，李老师有更好的建议吗？
谢谢。:146:

@shiyu 在 BoundaryFoam求解近壁面流动（倾斜壁面） 中说：

已知local坐标系下的值，想转换成global坐标系下，方便给定openfoam的边界值，openfoam里有可以直接利用的tools吗？

没有现成的。不过你写代码的话，也不是特别费事。

http://www.dyfluid.cn/theory.pdf

合适的内容我会更新在《无痛苦NS方程笔记》里面。

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

捕获.PNG

libs ("libmyfvOptions.so"); 这句不能像一些官方提供的案例那样放在fvoptions文件里，得放在controldict文件下
Screenshot from 2021-12-03 17-31-56.png

@zhouxu 我就是设置的滑移边界，存在槽道流特征

测力时出现很强烈的震荡 想请教下是什么原因导致的 黑色的是想要验证的图（没有震荡），红色的是我计算的。用的Openfoam自带的测力工具。
17c6ea16-98db-4fc2-8b0f-25c32d22a333-image.png
ea932f4f-4999-469d-bdb4-d254f739078e-image.png a83847ec-c340-453c-a6c1-37a0913808f8-image.png

@李东岳 好的，感谢李老师回复，我们服务器系统是centos6.5，好像不太支持新版本:136:

@caijinjin 是计算停在这里还是终端会停止，方便的话可以交流一下吗，我微信15098932699

按说，那个 pdf 的大部分勘误记录都会跳转到这里
https://cfd-china.com/topic/2338/无痛苦n-s方程笔记-勘误/38

“错”和“误”不一样，（所以某些理论上无法消除、工具条件限制带来的与事实不符的差距叫“误差”不叫“错差”）（我对你措辞没意见，只是看你挺有互助热情，也“互助”你一下，说给你听听，参考参考。）

可以直接说出笔误位置，细微的笔误难以发现（容易发现也不会误了）。比如：

"接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
这句第二个坐标（或者第二个（x1, y1, z1））标号应该是(x2, y2, z2)

我看你之前也发过类似的帖子，可以看完了一遍书，把发现的笔误统一发一下，也方便修改。

欢迎继续:146:

@chaoscfd 我用了 但还是mesh.findcell最快

我猜是这样的，
process 里必须有动网格，这样计算的过程只涉及网格如何动，

如果 process 里没有动网格，那么当动网格运动到该区域，还要造动网格。

至于 process 里没有动网格也能算，那是动网格没有运动到该 process 覆盖的区域。

可以测试，在计算区域内设置动网格，但是不让它动，比如平静的水面。应该可以一直算。

paraview 里可以看每个 process 的区域，算不动的那步估计就是跨区域运动了

顶一下

这个是安装 docker 的问题，不是 OF 的问题 :wocao:
去 docker 相关社区问更容易得到答案:chigua:

@tens 再次感谢您！重新尝试了一下，可以输出了！