555c7ea9-b5fe-484b-ae06-64ec9cbd376d-image.png 这个是我想求解的方程形式，这段时间看了点相关东西，请问各位前辈，foam-extend里的blockMatrix能实现我想要的求解方式吗，然后就是曳力项应该OF里应该怎么写呢？是fvm : : Sp ( K, Ug ) - fvm : : Sp ( K, Ul )吗？希望各位大佬不吝赐教啊:mianmo:

https://www.cfd-online.com/Forums/openfoam-solving/181258-vof-model-mixture-model-ras-les-model.html

https://www.openfoam.com/documentation/guides/v2112/doc/

右上方有个小窗口，从那里搜相关关键词，会快一些定位错误

也可以编译本地的这种文档，比网页更快一些

有可能。非结构网格有的时候算出来是这样。这算的查的还不是特别多。纯粹是数值问题。

由于相变的存在，导致interPhaseChangeFoam的压力方程中额外添加了由相变产生的速度散度项，这一项是直接由空化模型得到的。如果我并没有使用空化模型而是使用别的方法的得到了每个时刻的alpha场，如何通过新的alpha与旧的alpha构造出压力方程中的速度散度呢

@李老师和各位老师：
LES 的 inlet 的初始条件，好像就是先运行 （同条件的） RANS，然后 把下游某处 map 给 LES 的 inlet

刚接触 LES，还有个问题要和老师们确认：
波浪的 LES，还有必要人为添加扰动吗？

谢谢

@Samuel-Tu 我在用拉格朗日方法计算颗粒的时候也出现了相似的问题，单核能算，但是并行算几步就不行了，想问一下你当初怎么解决了，你楼上的帖子说“往scalarField里面加数据的时候超过了定义的大小”我没看懂。
920e3c52-5310-4117-85d0-015244998d03-image.png

@李东岳
明白了，谢谢李老师。。我另外那问题，还要再麻烦您或其他老师帮指点下

@李东岳 好吧，谢谢大佬！

@Foamer 可以设置不同的边界，然后设置不同的润湿角？

@lys 你可以看下你的solver中有没有添加storeGlobalPositions这个函数，我对你使用的这个solver不是很了解，拉格朗日是你自己添加进去的吗

请问大家，一般对风力发电机方向的仿真一般基于哪个标准求解器？然后怎么做出过那种涡流线的结果图？

壁面的速度确实不能为0，壁面的surface是0. 从这个图来看，【猜测】是粒子直径太大了导致的。这个速度场不太正确。

@李东岳 嗷嗷，好的，谢谢李老师，我再调研一下。

看不出具体原因,只能一步一步排查,最后减小时间步长,可以运行.

@李东岳 好的好的，谢谢大家

在OpenFOAM5.X之后的版本中都改了这个position文件的记录方式，那如果想得到粒子的位置信息从而得到轨迹，是不是需要使用p.position函数得到坐标矢量，存到新的矢量变量中，这样就可以得到粒子坐标轨迹了

@cccrrryyy
谢谢老师，关键2点总结如下：

1 . 您的 “。。。同一个case里面网格如果有变化，那filter体积就是在变化的。。。” 意思是：
网格密（即小）的地方，vortices 模拟的量，要少于网格稀（即大）的地方

您的 “。。。他只是在举个例子。。。”， 意思是：
（尽管如此），但那仍是个 "laminar case"， (as no turbulence variables, e.g. vorticity etc. are computed / genertated during case run)

再次感谢老师的指点

我也遇到这个问题,版本v1912,调用了表面采样,截面穿过了AMI,计算的时候把界面关了,reconstructPar成功了

尊敬的李老师:
现有几个东岳流体网interFoam解析部分的几个小的建议，请您过目:
1.公式（24）是否应该为 \begin{equation*}\nabla\cdot\tau=\nabla\cdot\left(\mu\left(\nabla\bfU+\nabla\bfU^{\mathbf{T}}\right)\right)
\end{equation*}
2.公式（28）之上的那一句话，“方程(28）”是否应该写为"公式（3)"
3.公式(29）是否可以写作
\begin{gather*}
\frac{\mathrm{D}\rho}{\mathrm{D}t} = \frac{\mathrm{D} \left(\alpha_1\left(\rho_1 - \rho_2 \right)+\rho_2 \right)}{\mathrm{D}t}= (\rho_1-\rho_2)\frac{\mathrm{D}\alpha_1}{\mathrm{D}t}
\end{gather*}
$$\frac{\mathrm{D}\alpha}{\mathrm{D}t}=\frac{\partial \alpha}{\partial t}+\bfU\cdot\nabla\alpha=0$$

Many thanks.