好贴，谢谢李老师分享！！！

@小鲨鱼 在 OpenFOAM涡环模拟 中说：

数学上是有一个积分公式实现涡量场到速度场的转化的

写出来看看

如果可以的话，可以有偿指挥。

老哥，可以分享一下算例吗？我正在学习

这个就是发散了，目前看不出来什么原因

各位大神好，我在做一个关于可压缩流体的动网格算例时遇到了问题：我的模型是一个长方体，当我将x，y方向上的边界都设置为slip时，z方向的边界沿着z轴移动，从而实现压缩这个模型，这样是可以计算的。
但是当我想要实现沿x和z两个方向都压缩时则出现了问题，此时x和z的边界都沿着各自的轴移动，这样在x和z边界的接触的位置网格产生了畸变。想请教一下应该如何设置？

@李东岳 感谢东岳老师的回复。已解决。需要在system里有一个单独的forceCoeffs文件，然后在controlDict里的添加functions{#include "forceCoeffs"}，这样会在输出文件里输出每一步的Cd等相关信息。感谢老师！

@Dingcy 在 如何输出出口escape的颗粒数量随时间变化的曲线？ 中说：

如何输出出口escape的颗粒数量随时间变化的曲线？之前算单向流的时候我输出出口浓度随时间变化曲线是通过paraview那个plot selection over time, 但是现在这个颗粒好像不能那么整了....:136:
1c05a757-a68a-499f-bfc9-1582a3c49a26-image.png
前辈您好，请问这个操作怎么实现，得到下图:134:

我在无痛苦ns方程笔记里面增加了对这一行代码的解释：

fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY??? - fvm::laplacian(alpha1alpha2f()*pPrimeByA_(), alpha1) ); alpha1Eqn.relax(); alpha1Eqn.solve();

尤其是其中的

fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY???

@Tens 恭喜恭喜！！！！隔壁师兄又要嫉妒死了！:mihu: :mihu: :mihu:

@李东岳 多谢李老师:146:

@李东岳 感谢东岳老师

@pengdi LES这么添加是不对的。你那个公式是RANS的k，不是k_sgs。LES那面如果你用smagorinsky，不需要给k

Distribute the map ... [cu01:53769] *** An error occurred in MPI_Bsend [cu01:53769] *** reported by process [140195840131073,17] [cu01:53769] *** on communicator MPI COMMUNICATOR 3 SPLIT FROM 0 [cu01:53769] *** MPI_ERR_BUFFER: invalid buffer pointer [cu01:53769] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort, [cu01:53769] *** and potentially your MPI job)

@李东岳 谢谢李老师，我再回去看看是不是我公式的问题

我看了你的一些帖子，大部分都是基本的使用问题。这一类问题在学习用户手册后都可以被解决。 http://www.dyfluid.com/OpenFOAM用户指南-9.pdf

先处理出来涡量 ： http://dyfluid.com/code.html#id29 使用streamLine函数可以处理某个平面的速度，压力，涡量等值： http://dyfluid.com/code.html#id40

建立两个网格就行了。P一套网格。Q一套网格。

思路具体是需要对OpenFOAM中某一求解器中的头文件进行修改和变化，并修改其中其中一些项。是自己建立一个求解器然后将原来的自带求解器中代码粘贴到建立的求解器中修改，然后编译的过程吗？有什么注意事项吗？

@李东岳 谢谢东哥

同问一下，我怎样能够提取某一截面中某一部分的速度，感谢