@李东岳 谢谢李老师的回复，期待:xiexie:

多谢老师！

谢谢啦

多谢老铁，$\frac{1}{\Delta t}$已更正为$\frac{\rho_\rP^{*}}{\Delta t}$

想复现溃坝一物理模型的数据，需要保持上游水头一直不变，请问各位大神这怎么操作

@麦迪文 你好，请问能参考下您的代码吗: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涡环模拟 中说：

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

写出来看看

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

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