CFD中文网

@姜磊 大佬 你有动态接触角的UDF了嘛 求大佬分享 2529693979@qq.com 万分感谢

@linhan-ge 我觉得也应该是这样，但我还没理解李老师说的，您理解了吗？

@李东岳 delta函数

@Samuel-Tu
嗯，好的。毕竟这个也没有实验数据可参考，也只能这样了。

@李东岳 感谢李老师，问题已解决

@李东岳 图方便，我只带了一个商务本，，，周一去实验室试试
笔记本买不啊，这头的实在是贵:xinlei:

如果单纯的mules求解每一相是无法保证alpha1+alpha2=1的

现在我就想确定，alpha2通过对应方程推导，用MULES求解不可行的吗？

MULES的出发点是保证变量的有界。所以如果你用MULES求解alpha1，那么理论上alpha1是有界的。alpha2的求解没有调用MULES，他是通过1-alpha1算出来的，如果alpha1有界，alpha2也有界。

如果你用算法同时求解alpha1和alpha2，如何处理耦合？这俩个变量是耦合在一起的，就像速度和压力。你不能单独的去分离求解，如果分离求解就就需要迭代。迭代就导致计算速度变慢。因此现存大厂据我所知都是只求解alpha1，然后alpha2=1-alpha1。

我在去你年底验证了MULES，理论上是可以保证有界，但是真实计算的时候，还是会越界。尤其是物理模型比较复杂的时候。个人觉得这方面内容搞出来，绝对是个好文。

@程迪
您好，有个问题想请教您。看到您在《OpenFOAM的残差定义》一文中提到了残差定义方式：https://chengdi123000.github.io/2018/01/04/OpenFOAM的残差定义/
804487c1-8422-49c6-a9b1-116b1949a77f-image.png
关于分母normFactor，有些疑问，不知道为什么要这样定义，有什么数学上的原理吗？
谢谢！

好的老师，我算完给您回复。

@李东岳 对的，大概有十几公里，那谢谢东岳老师，我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了，虽然能看懂离散和求解过程，但是想在openfoam里面实现觉得难度挺大的，当时老师您带我们植入的单相可压缩的，目前这是两相可压缩，我努力努力，再次感谢东岳老师的持续关心:xinxin:

If you would prefer to write Python, and can afford to write Python, we recommend using the Python interface to PyTorch. However, if you would prefer to write C++, or need to write C++ (because of multithreading, latency or deployment requirements), the C++ frontend to PyTorch provides an API that is approximately as convenient, flexible, friendly and intuitive as its Python counterpart. The two frontends serve different use cases, work hand in hand, and neither is meant to unconditionally replace the other.

我理解的应该是对的。pytorch，libtorch看个人爱好

https://pytorch.org/cppdocs/frontend.html

偶尔误入这个帖子，还挺有意思。

鼠标垫基本都发货了。Ricardo，你的地址信息不全

你好，你现在解决这个问题了吗？我现在也遇到了这个问题，不过是因为pressureDirectedInletVelocity 这个进口条件的问题。
我也在做类似的研究，可以交流一下吗

如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

@李东岳 模型缺陷吗.....那就有点麻烦了。我不单纯写文章，是老师的项目。我再考虑考虑

@尚善若水 应该可以吧。RR[speciei][celli]是celli网格里speciei的反应速率，那该网格里该组分的Qdot应该可以这么写：

scalar Qdot-speciei-celli = hi * RR_[speciei][celli];

@Dingcy 哈哈我是在用低雷诺数湍流模型算超临界CO2流动传热的例子，物性变化很剧烈，不好收敛。

@exthan 您好，我刚刚也下载了OF7版本的makeAxialmesh,但好像没找到编译文件，没有Allmake文件或者wmake文件，请问您是怎么编译的呢d27bf862-45db-4120-8ba1-050e18d46128-image.png

@东岳 对，我那个5118是学校的计算集群，我记得我在哪儿看到过说这堆集群都是把超线程关了的，所以不可能使用的线程数量超过物理核心数量；然后睿频的话，应该是说如果你不关他，他默认是开的，但学校集群可能是关了的，你的如果是自己用应该是开着的，就是说如果睿频开着的，CPU温度和功率不超过他的限值的时候他可以运行在更高的频率上。比如你的5120 是14cores 28 threads 2.2G, 但是如果散热比较好功耗不超过某个限制的时候，他是可以整体跑在更高的频率上，这就能更快。
话说我觉得最惊奇的是我的9750H跑6线程居然比服务器CPU慢这么多，理论上说不过去啊，因为我看到跑的时候主频都跑到4G了，服务器才2.3G。得研究一下构架了