@李东岳 谢谢师兄，因为模型有辐射模型，流固耦合，还有多物种传输，所以一直觉得是自己设置有问题，或者部分模型需要重新编译，忘了最基本的，刚才把所有收敛系数都调小了.可以跑了，应该就是前期不稳定。但是，也只是多跑了3步，我回头再看看

@李东岳 哦哦，懂了，那我都试试吧。谢谢李老师及时回复。:xiexie: 感恩的心:xinxin:

@李东岳 李老师，这是我自己编译的一个求解器，我在constant文件夹下将g从（0 0 -9.8）改为（0 0 0），后消失了，但是不符合正常要求呀

fvSolution里面的pFinal，relTol应该调到0吧。
矩阵最终求解的这些*Final一般都强制要求满足容差。
应该是这个原因导致你某一步没收敛，或者误差逐渐积累，导致发散。

找到数值发散的原因了，是热流密度设置过大

谢谢啦，主要想找一个带相变的可压缩气液两相流的solver，我再找找试试~

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

@小鲨鱼 在 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)

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