@麦迪文 你好，请问能参考下您的代码吗:xiexie:

@chengan-wang 怎么修正啊

@dzw05 @溯光 @李东岳 您好，我想请问一下这个interface里面的C表征的是流体对粒子一个什么样的力呀？我最近在验证波浪作用下一个粒子的输移运动，发现改变不同的C会使得粒子在水平和垂向位移上与stokes质量输移有较大的偏差

image.png

各位老师好，想请教一个问题，我现在在controlDict文件里写了一个function，里面需要读取一个字典中的值。字典中的值为一个标量，读取之后想改变这个值并重新写入字典中，这样下次读取之后就是新的值。
其中字典文件定义如下
1.png
function里面读取的代码如下
2.png
我尝试用write函数来改写，但是报错是这样
3.png
请问各位老师这个该怎么解决呢？谢谢！

@Wayne 你好，请问一下这个问题怎么解决？

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